- Open Access
Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds
© Vinkler et al.; licensee BioMed Central Ltd. 2014
- Received: 11 June 2014
- Accepted: 10 October 2014
- Published: 12 November 2014
Toll-like receptors (TLR) are essential activators of the innate part of the vertebrate immune system. In this study, we analysed the interspecific variability of three TLR (bacterial-sensing TLR4 and TLR5 and viral-sensing TLR7) within the Galloanserae bird clade, investigated their phylogeny, assessed their structural conservation and estimated site-specific selection pressures.
Physiochemical properties varied according to the TLR analysed, mainly with regards to the surface electrostatic potential distribution. The predicted ligand-binding features (mainly in TLR4 and TLR5) differed between the avian proteins and their fish and mammalian counterparts, but also varied within the Galloanserae birds. We identified 20 positively selected sites in the three TLR, among which several are topologically close to ligand-binding sites reported for mammalian and fish TLR. We described 26, 28 and 25 evolutionarily non-conservative sites in TLR4, TLR5 and TLR7, respectively. Thirteen of these sites in TLR4, and ten in TLR5 were located in functionally relevant regions. The variability appears to be functionally more conserved for viral-sensing TLR7 than for the bacterial-sensing TLR. Amino-acid positions 268, 270, 343, 383, 444 and 471 in TLR4 and 180, 183, 209, 216, 264, 342 and 379 in TLR5 are key candidates for further functional research.
Host-pathogen co-evolution has a major effect on the features of host immune receptors. Our results suggest that avian and mammalian TLR may be differentially adapted to pathogen-derived ligand recognition. We have detected signatures of positive selection even within the Galloanserae lineage. To our knowledge, this is the first study to depict evolutionary pressures on Galloanserae TLR and to estimate the validity of current knowledge on TLR function (based on mammalian and chicken models) for non-model species of this clade.
- Root Mean Square Deviation
- Human TLR4
- Phenetic Analysis
- Electrostatic Potential Distribution
- Datamonkey Server
Toll-like receptors (TLR) are part of the pattern recognition receptor machinery and play a key role in initial pathogen recognition in vertebrates . Since TLR are responsible for the recognition of microbe-associated molecular patterns that are present on pathogens , the efficiency of innate immunity in vertebrates is conditioned by their optimal functioning. Substantial variability has been described for TLR, both at the interspecific and intraspecific level . Variability in the structure and binding features of TLR could significantly influence host resistance to diseases and vulnerability to autoimmune damage. TLR evolution has been intensively studied in vertebrates in general - and within the mammalian clade in particular -. In birds, although the number of studies on the evolution of TLR has steadily increased -, there is still a very limited understanding of the functional significance of the putatively adaptive variability observed.
There are ten TLR genes in birds ,,, including those encoding predominantly bacterial-sensing TLR, such as TLR4 and TLR5, and viral-sensing TLR, such as TLR7. Although the TLR and their corresponding genes have been well characterised in the domestic chicken (Gallus gallus domesticus-), our knowledge on TLR of other avian species remains rather fragmentary ,-. Based mainly on data for the domestic chicken, we assume that the ligands recognized by avian TLR are similar to their mammalian counterparts , i.e. TLR4 binds lipopolysaccharide (LPS) , TLR5 binds flagellin  and TLR7 binds viral single-stranded RNA and synthetic antivirals ,. The structure of bacterial ligands, such as LPS or flagellin, varies between species -, which may exert selective pressure on TLR and lead to host-pathogen co-evolution of these molecules . Previously, it has been shown that TLR-ligand binding is, in many aspects, species-specific and that TLR protein variation may have a functional significance ,,- that could affect resistance to pathogen infections ,.
The function of the receptor is determined by its tertiary structure and surface features that confer specificity for ligands. TLR are type I integral membrane glycoproteins that are characterized by an extracellular ligand-binding domain (ECD) and a cytoplasmic signalling Toll/IL-1R homology (TIR) domain. Three-dimensional (3D) molecular structures of the ECD - and TIR domains , have been described for several human and mouse TLR. In several other cases, the ECD structures are predicted -. Most recently, the TLR5 ECD structure has been resolved in zebrafish (Danio rerio). These studies have identified TLR ligand-binding sites. All these protein structures, however, represent mammalian or fish TLR only, no avian TLR have been characterised or even predicted in 3D as yet.
In this study, we provide a comprehensive view on the evolution of TLR in Galloanserae birds. Galloanserae is an ancient clade of avian species that includes the orders Galliformes and Anseriformes which are largely separated from all other modern birds of the Neoaves lineage . Currently, this is the most extensively studied avian taxon regarding TLR immunogenetics. Based on published data, we investigated the sequence and structural variability of TLR4, TLR5 and TLR7 within this bird clade. These three TLR were chosen as representatives of the bacterial-sensing and viral-sensing TLR based on published sequence data. Since comparison of 3D structures and protein surface features may reveal biologically interesting similarities not detectable by sequence analysis , the 3D tertiary structure of these three proteins was modelled. We then carried out a structural comparison of functionally important regions, a comparison of surface electrostatic potentials and four independent analyses of positive selection. Superposition of the 3D structures allowed us, for the first time, to conduct a phenetic (not only phylogenetic) analysis of avian TLR evolution.
Coding DNA sequences (CDS) for the selected TLR were downloaded from the National Centre for Biotechnology Information (NCBI) GenBank for all Galloanserae species currently available in full length (n = 13) and for humans and mouse. GenBank ID are given in Additional file 1: Table S1. CDS translation was performed using BioEdit Sequence Alignment Editor (Tom Hall, Ibis Biosciences, Carlsbad, California, USA) and protein sequences for each gene were aligned using ClustalW multiple-sequence alignment software. Alignment of nucleotide codons and their corresponding protein sequences was undertaken using the PAL2NAL webtool .
To predict the distribution of structural domains in the proteins, we applied SMART . Since the analysis revealed differences between species in the number and position of leucine-rich repeats (LRR) we tested LRR distribution by an independent approach using the LRRfinder, with upper and lower boundaries fixed at 95% and 80%, respectively . Because the results of these two approaches differed slightly, SMART predictions were only used to identify N-terminal LRR (LRRNT) and C-terminal LRR (LRRCT) motifs which in most cases were not identified by LRRfinder, while other LRR were identified based on LRRfinder predictions. Molecular weight and charge at pH = 7 were also calculated for each predicted protein. Presence of a transmembrane (TM) region in each protein was checked for on the DAS-TMfilter server . When the TIR domain was not detected by SMART, we used PFAM comparison  to check for its presence. Signal peptides were identified using SignalP 4.0 . Finally, secondary and 3D tertiary structures of the three TLR were predicted by applying a comparative modelling approach (see e.g. , or ) using I-TASSER , which is currently the leading protein structure prediction server (see ) and http://predictioncenter.org/casp10/)). The I-TASSER server uses a hierarchical protein-structure modelling approach based on secondary-structure enhanced profile-profile threading alignment and iterative implementation of the threading assembly refinement program . For the prediction of secondary structures, we used the whole CDS, whereas ECD and TIR domains were modelled separately in order to compare tertiary structures. Since it has been shown that the sequence for TLR7 up to LRR14 is cleaved in the endoplasmic reticulum , we modelled the sequence starting at amino-acid 417 only. Models with the highest C-scores and conformation similarities to other modelled structures were used for further analysis. In these models, we excluded regions with limited structural stability (i.e. regions with high modelling errors: signal peptides and regions following LRRCT). Modelling errors in the regions of interest were estimated using ModFOLD ; all models had high levels of confidence with P-values less than 0.002 and Global model quality scores greater than 0.37. Although the accuracy of our models may still have been limited, the aim of this study was not to describe the proteins’ tertiary structures precisely but to assess average structural similarity between the receptors. The error estimates obtained indicate that the models constructed represent reasonably reliable inputs for further phenetic analysis. Images of the predicted protein 3D structures were visualised using PyMOL software v. 1.5 (http://pymol.org/). Protein electrostatic potentials were calculated using PDB2PQR v. 1.9.0  (PDB2PQR Server, http://nbcr-222.ucsd.edu/pdb2pqr_1.9.0/) based on the PARSE force-field and electrostatic calculation in the APBS web solver  (http://www.poissonboltzmann.org/). Surface charge distribution was visualised using Jmol v. 12.2 (http://jmol.sourceforge.net/).
Phylogenetic and phenetic analysis
Alignments of TLR4, TLR5 and TLR7 CDS were used for phylogenetic analysis using a maximum likelihood (ML) method. As outgroups, we used orthologous human and mouse sequences. FindModel (http://www.hiv.lanl.gov/content/sequence/findmodel/findmodel.html) was used to evaluate the fit of 28 nested nucleotide substitution models to the data, the best model for each alignment being selected on the basis of the Akaike information criterion. ML analyses were performed using PHYML , with the NNI algorithm and BIONJ distance-based tree as the starting tree. Bootstrap analysis (with 1000 replicates) was performed to estimate the robustness of internal nodes. The results were visualised in FigTree v. 1.3.1 (http://tree.bio.ed.ac.uk/software/figtree/). A consensus phylogenetic tree including all investigated species (see Additional file 2: Figure S1) was constructed using the avian phylogenetic tool available at http://birdtree.org/. We used the Hackett backbone  as the source tree with 1000 randomly generated trees. The maximum clade credibility tree was produced using the TreeAnnotator v. 1.8.0 tool in BEAST v. 1.8.0 software . Phenetic similarity analysis of the predicted protein secondary and tertiary structures was performed to detect conserved structures in avian TLR. Secondary structures were compared based on alignments obtained using the EMBOSS Needle pairwise alignment tool (http://www.ebi.ac.uk/Tools/psa/emboss_needle/), using chicken (Gallus gallus) GaGaTLR sequences as references. To predict 3D structures, we used the adjusted I-TASSER pdb models for structural superposition in the DALI pairwise comparison tool . The pair-wise root mean square deviations (RMSD) metric was used to compare protein structures , and to construct distance matrices that subsequently served as matrices of Euclidean distances in cluster analysis using an unweighted pair group method with arithmetic mean (UPGMA) method in STATISTICA v. 6.0 (StatSoft, Inc., Tulsa, OK, USA; ; for a similar approach see ).
Before testing for selection, all codons containing gaps in any species in the alignment were removed (this applied to only six codons in TLR5; throughout the text, codon positions are numbered according to the chicken GaGaTLR5 sequence (see Additional file 1: Table S1). We used two methods to test for positive selection on individual residues at the interspecific level within the Galloanserae clade, i.e. (1) the hierarchical Bayes (Bayes Empirical Bayes, BEB) approach with implemented Markov chain Monte Carlo routine - PAML (Phylogenetic Analysis by Maximum Likelihood ); and (2) FUBAR (A Fast, Unconstrained Bayesian AppRoximation for inferring selection ). For PAML (v. 4.6), we used codon-based substitution models (codeml) to identify amino acid sites under positive selection in the CDS comparing the neutral M8a (beta&ω = 1) model with the alternative M8 (beta&ω) model. The likelihood ratio test (LRT) for the comparison of two nested models was calculated using chi-square approximation: Chi2 = 2 × (lnLM8 - lnLM8a), where LM8 and LM8a are likelihood values. The degrees of freedom (df) were defined as the difference in the number of parameters in the models used (see Additional file 1: Table S6). If the LRT was less than 0.05, positive selection was considered significant. The BEB approach  was used to determine site-specific posterior probabilities of positive selection (≥0.9). FUBAR analysis was performed on the Datamonkey server (http://www.datamonkey.org/, ) using a default significance level of posterior probability set at 0.9. In this study, we applied the FUBAR algorithm because it is more robust and much faster than selection analysis based on random effect likelihood (REL methods ).
We tested the degree of dissimilarity of amino acid substitutions according to their physiochemical properties using the new PRIME (PRoperty Informed Model of Evolution) tool available on the Datamonkey server . We used the set of five composite physiochemical properties proposed by Atchley et al. , i.e. polarity index, secondary structure factor, volume, refractivity/heat capacity and charge/iso-electric point. A change in these properties was considered significant if the posterior probability was greater than 0.9. The evolutionary conservation of amino acid positions was predicted using the ConSurf tool , with the assumption that positively selected residues (functionally important for pathogen binding) were the least conserved. For all ConSurf analyses, we used GaGaTLR protein 3D models (obtained as described above) and the LG substitution matrix . A phylogenetic tree of the three TLR genes studied including all investigated species in the PAML and ConSurf analyses was constructed as described above (see Additional file 2: Figure S1).
We were able to verify the homology of all sequences examined (see Additional file 2: Figure S2) and subsequently to assess amino acid identity and similarity (see Additional file 1: Table S2). The phylogeny of the TLR genes was consistent with known phylogeny for the Galloanserae clade  (see Additional file 2: Figures S1 and S2). While the size of both TLR4 and TLR7 was uniform (843 aa and 1047 aa, respectively), TLR5 varied in length between 859 aa and 862 aa (see Additional file 1: Table S3 and Additional file 3: Section S1). The orthologues showed little variation in molecular weights, although some differed markedly in their charge at pH = 7 (see Additional file 1: Table S3). Since localised charge variability can influence protein conformation and domain composition and produce variation in ligand-binding features, we examined charge and structural variation in more detail.
Protein structure evolution
TLR proteins are likely to be involved in host-parasite co-evolution, and thus shaped by parasite-mediated natural selection . Interspecific differences in TLR protein structure, therefore, may exhibit imprints of structural evolutionary convergence due to selection. In this study, for the first time in birds, we were able to model secondary and tertiary structures of TLR4, TLR5 and TLR7 for all Galloanserae species with currently known CDS. Secondary structures for the proteins predicted by I-TASSER revealed low levels of interspecific structural variability within individual TLR (see Additional file 2: Figure S3). Although two regions of potential functional interest in TLR4 were polymorphic (see Additional file 3: Section S2), the most important structural motifs in all three receptors were conserved in all species analysed, resulting in more than 90% interspecific identity in secondary structure distribution (see Additional file 1: Table S4). As a result, 3D extracellular domain models invariably had a horseshoe-like shape, in which the concave surface comprised β-sheets and the convex surface parallel loops and short helices. The TIR domain (which was modelled separately) had a globular shape. A phenetic analysis of RMSD distances obtained by superposition of the modelled structures revealed that structural variability in the TIR domain of all three receptors was lower than ECD variability (see Additional file 2: Figures S4d, S4e and S4f). This may be a result of higher conservation in the TIR domain when compared to ECD , although domain size may also have played a role as RMSD tends to increase with protein size. Avian TLR4 ECD showed a stronger structural resemblance to human TLR4 than mouse TLR4 (see Additional file 2: Figure S4a), which suggests that avian TLR4 may exhibit similar binding features to human TLR4. At present, however, this must remain a hypothesis since current experimental data provide no support . As the RMSD of our superposed models were below 2.5 Å (i.e. deviation between models was lower than the accuracy of individual models), we conclude that TLR4 and TLR5 ECD structural variability within the Galloanserae is generally low and probably unimportant (see Additional file 2: Figures S4a and S4b). Of more interest is the ECD of TLR7, with the phenogram indicating that Anseriform TLR7 exhibits a close structural relationship to human TLR7, while Galliform TLR7 clusters with murine TLR7 (see Additional file 2: Figure S4c). Despite this, the RMSD were too low to indicate any meaningful structural variation. The low TLR structural variability observed in this study is consistent with recent findings for rodents .
Surface electrostatic potential
Changes in physiochemical features of predicted ligand-binding residues
To further explore potential interspecific variability in the ligand-binding properties of avian TLR proteins, we compared basic chemical features of residues known or predicted to be functionally important in other vertebrate TLR ,,,,,. These included residues directly involved in ligand binding, TLR homodimerization and MD-2 heterodimerization. At most of the predicted functional sites in avian TLR, similar chemical properties were preserved as in mammalian TLR (see Additional file 1: Table S5). Such conservation was particularly visible in TLR7, which suggests that most receptor-ligand interactions in TLR7 are highly conserved (see Additional file 3: Section S3). In TLR5, we identified substantial residue differences at the binding region previously described in zebrafish  as separating fish from amniotes (see Additional file 1: Table S5). Primary binding interfaces A and B, as well as both dimerization interfaces, appear to be only modestly conserved, indicating that flagellin binding probably differs in fish and amniotes. This is further supported by the fact that avian TLR5 sites generally retain their physiochemical properties at those binding residues predicted for mammals. Intriguingly, we identified potentially important changes in amino acid properties at several TLR4 functional sites, which suggests variability in LPS binding and MD-2 dimerization (see Additional file 1: Table S5). In a number of cases, residue changes precluded the existence of charge (positions 268 and 397) or hydrophobic interactions (positions 449 and 472) known for human TLR4 ,. Interestingly, all avian species possess arginine at position 393; in this respect, avian TLR4 is identical to equine TLR4  but distinct from human TLR4. This suggests that lipid IVa in birds serves as an agonist triggering a TLR4-MD-2-mediated immune response, just as it does in horses . Galliformes also display systematic differences from Anseriformes at several positions (see Additional file 3: Section S3). Thus, our results suggest that the ligand-binding features of avian TLR differ not only from mammalian TLR but also between the various avian taxa. This is in concordance with previously reported experimental results ,.
Analysis of positive selection
While the vast majority of residues in TLR were evolutionarily highly conservative, consistent with purifying selection on TLR ,,, ConSurf analysis revealed 26 reliable evolutionarily non-conservative sites in TLR4, 28 in TLR5 and 25 in TLR7 (Figures 2, 3 and 4 and (see Additional file 1: Table S9 for detailed information). These may also represent positively selected sites since functionally important positively-selected residues are also predicted to be the least conserved. This is supported by the location of six of the non-conservative sites in the functionally important TLR4 signal peptide, four in the predicted TLR4-MD-2-binding regions ,, three in the LPS-binding regions ,, two in the TLR4-TLR4-dimerization regions  and one in the trans-membrane region (see Additional file 3: Section S4). One of the non-conservative positions identified in our Galloanserae dataset (383) has previously been identified as a possible single nucleotide polymorphism (SNP) responsible for Salmonella resistance in the domestic chicken . In TLR5, one non-conservative site was situated in the signal peptide, five in the flagellin-binding residues identified in fish , two in flagellin-binding residues predicted in mammals  and two in the trans-membrane region. Most of the non-conservative sites in TLR7 were situated in the excised region  (and hence likely to be functionally unimportant), with just one in the potentially relevant signal peptide and one in the TIR domain.
Intriguingly, both the positively-selected and non-conservative sites identified in this study partly match the positively selected sites detected in previous studies in other vertebrate taxa ,,,,,,, (Figures 2, 3 and 4, Additional file 1: Table S10 and Additional file 3: Section S4). In TLR4 and TLR5, 15 and 20 candidate sites, respectively, were either precisely the same or lay in close proximity (up to 2 aa) to a site identified in previous studies. In contrast, only seven sites were close to, or the same as, previously detected sites in TLR7. Most importantly, there was consensus for positive selection on site 270 in the predicted MD-2-binding region in TLR4 in birds, and on sites 183, 209, 264 and 342 in TLR5 both in birds and mammals, all of which are predicted flagellin binding sites or sites in close proximity to these binding sites (see Additional file 3: Section S4). No consensus site was located in the predicted ligand-binding region in TLR7. Independent concordance between the results of selection analysis in several studies that cover different taxa, and topological agreement between the identified sites and predicted functionally important regions, strongly support the proposed importance of diversifying selection or positive selection in the evolution of these residues. In human TLR4, for example, it has been demonstrated that even SNP located far from MD-2-dimerization or LPS-binding sites can modify responsiveness to LPS  or influence the binding of other ligands. Thus, other consensus selected sites identified in TLR4 (205, 246, 273–275, 301 and 406, Figure 2) and TLR5 (281, 293, 314, 422, 508, 525 and 532, Figure 3) may play an important role in the evolution of TLR-ligand binding in birds (see Additional file 1: Table S10). However, the potentially important consensus sites identified in TLR7, are located mainly in the excised region  (residues at positions 39, 99 and 383) and do not appear to influence receptor function. The only exception to this is site 665 which lies, however, outside the predicted ligand-binding region. Given the generally low concordance in predicted positively-selected sites between studies, we suggest that many of the sites reported for TLR7 represent false-positive predictions resulting from the method applied (most sites were identified by REL or FEL; see  for discussion).
In concordance with previous studies, we found stronger signatures of positive selection acting on the ligand-binding regions in TLR4  and TLR5  than in TLR7 . There are several possible explanations for this. First, bacterial-sensing TLR such as TLR4 and TLR5 recognise structurally variable ligands -, while viral-sensing TLR7 detects structurally invariant RNA molecules regardless of their precise sequence . Hence, TLR7 is likely to have evolved mainly under purifying selection, while TLR4 and TLR5 evolved mainly under diversifying selection. Furthermore, since TLR4 is capable of recognising several unrelated ligands , it is also possible that this receptor evolved more rapidly. Second, there appears to be greater redundancy in bacterial recognition than viral recognition. In TLR5, for example, several studies have proposed relaxed purifying selection ,, while nonsense stop-codons have been described in both birds  and mammals . This may be a result of the presence of other flagellin receptors which may compensate for malfunction of TLR5 (discussed in ). Finally, the limited number of positively selected sites detected in the TLR7 ligand-binding region may result from the limited knowledge about its precise location. While the precise mechanism of ligand-binding has been described for both TLR4 - and TLR5 , only crude predictions are available for TLR7 . Although all these explanations may be relevant to some extent, the relatively low concordance between results of multiple studies involving TLR7 (see Additional file 1: Table S10) tends to support the view that there is, indeed, only weak positive selection acting on TLR7.
Interspecific comparisons within the Galloanserae clade revealed relatively high sequence variability in all three TLR investigated. Such variation has been shown to influence the physiochemical properties of proteins. Despite high tertiary-structure conservation, evolutionary changes were manifested by alterations to protein surface characteristics, such as changes in electrostatic potential distribution. Importantly, not only does surface charge in Galloanserae birds differ distinctly from that of mammals, to some extent there are also distinct differences observed within the clade. These variations most likely affect receptor binding features, a theory that is consistent with the idea of a host-pathogen evolutionary arms race  in which any adaptation enabling a pathogen to escape host immunity leads inevitably to a counter-adaptation in host receptors that again enables pathogen detection. Co-evolution of this kind has previously been described in both human TLR4 and TLR4 of the bacterium Pseudomonas aeruginosa. In fact, there are now numerous examples of individual host TLR adaptations to pathogens known (summarised, for example, in ). What we presently lack, however, is a functional understanding of evolutionarily-tested beneficial innovations present at the interspecific level in vertebrates. In this study, selection analysis identified a number of positively selected sites, mostly in the ligand-binding ECD. As in previous studies, we observed stronger positive selection acting on the ligand-binding regions of TLR4  and TLR5  compared to TLR7 . Sites subjected to selection in TLR4 and TLR5 were frequently located either precisely in, or in close topological proximity to, ligand-binding sites known or predicted in mammalian or fish TLR. This suggests that, although avian TLR may be differentially adapted to pathogen-derived ligand recognition compared to that of other vertebrate species, identical regions are responsible for ligand binding. We suggest that future investigations in this field should focus on functional testing of evolutionarily relevant substitutions detected by selection analyses. Based on the evidence summarised above, we propose several sites for more detailed investigation. In particular, sites 270, 343, 444 and 471 in TLR4 and sites 183, 209, 216, 264, 342 and 379 in TLR5 represent key candidates for further research on the functional significance of selection acting on TLR in birds. We suggest that site 268 (and possibly also site 383) in TLR4 and site 180 in TLR5 may be of particular evolutionary importance in Galloanserae birds since no selection on these sites has previously been observed in either mammals or other birds. Furthermore, while we do not know the functional importance of the sites, concordance between selection analyses in both birds and mammals suggests that special attention should be paid to positions 244–246, 273–275, 301 in TLR4 and 293, 294, 314, 342, 422, 525, 532, 533 in TLR5.
Taken together, our results have identified variability in Galloanserae birds that very likely results from pathogen-mediated evolution of species-specific TLR binding features. To our knowledge, this is the first study to depict evolutionary pressures on Galloanserae TLR and to estimate the validity of current knowledge on TLR function (based on mammalian and chicken models) for non-model species of this clade. Functional testing of the importance of individual sites (such as that performed by Walsh et al.  in mammals) should provide novel understanding of evolutionary mechanisms increasing resistance to pathogens in avian species. Any knowledge gained would be of great practical relevance, with applications in animal breeding for increased resistance to diseases.
MV conceived and designed the study, participated in the protein structure analysis and drafted the manuscript. HB carried out the selection analysis and participated in the protein structure analysis. JB performed the phylogenetic and phenetic analyses. All authors read and approved the final manuscript.
We are grateful to Marian Novotný, Zuzana Bainová and Anežka Pojezdná for their help and advice, and Dagmar Vinklerová and Kevin Roche for language correction. This study was supported through Czech Science Foundation Projects P505/10/1871 and P502/12/P179 and through Institutional Research Support (Grant SVV 260 087/2014). The work was also partially supported by the Ministry of Education, Youth and Sport of the Czech Republic within the NextGenProject (CZ.1.07/2.3.00/20.0303).
- Palsson-McDermott EM, O'Neill LAJ: Building an immune system from nine domains. Biochem Soc Trans. 2007, 35: 1437-1444. 10.1042/BST0351437.View ArticlePubMedGoogle Scholar
- Akira S, Uematsu S, Takeuchi O: Pathogen recognition and innate immunity. Cell. 2006, 124: 783-801. 10.1016/j.cell.2006.02.015.View ArticlePubMedGoogle Scholar
- Vinkler M, Albrecht T: The question waiting to be asked: Innate immunity receptors in the perspective of zoological research. Folia Zool. 2009, 58: 15-28.Google Scholar
- Roach JC, Glusman G, Rowen L, Kaur A, Purcell MK, Smith KD, Hood LE, Aderem A: The evolution of vertebrate Toll-like receptors. Proc Natl Acad Sci U S A. 2005, 102: 9577-9582. 10.1073/pnas.0502272102.PubMed CentralView ArticlePubMedGoogle Scholar
- Temperley ND, Berlin S, Paton IR, Griffin DK, Burt DW: Evolution of the chicken Toll-like receptor gene family: A story of gene gain and gene loss. BMC Genomics. 2008, 9: 62-10.1186/1471-2164-9-62.PubMed CentralView ArticlePubMedGoogle Scholar
- Mikami T, Miyashita H, Takatsuka S, Kuroki Y, Matsushima N: Molecular evolution of vertebrate Toll-like receptors: Evolutionary rate difference between their leucine-rich repeats and their TIR domains. Gene. 2012, 503: 235-243. 10.1016/j.gene.2012.04.007.View ArticlePubMedGoogle Scholar
- Wlasiuk G, Khan S, Switzer WM, Nachman MW: A history of recurrent positive selection at the Toll-like receptor 5 in primates. Mol Biol Evol. 2009, 26: 937-949. 10.1093/molbev/msp018.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang YH, Temperley ND, Ren LM, Smith J, Li N, Burt DW: Molecular evolution of the vertebrate TLR1 gene family - a complex history of gene duplication, gene conversion, positive selection and co-evolution. BMC Evol Biol. 2011, 11: 149-10.1186/1471-2148-11-149.PubMed CentralView ArticlePubMedGoogle Scholar
- Areal H, Abrantes J, Esteves PJ: Signatures of positive selection in Toll-like receptor (TLR) genes in mammals. BMC Evol Biol. 2011, 11: 368-10.1186/1471-2148-11-368.PubMed CentralView ArticlePubMedGoogle Scholar
- Fornuskova A, Vinkler M, Pages M, Galan M, Jousselin E, Cerqueira F, Morand S, Charbonnel N, Bryja J, Cosson J-F:Contrasted evolutionary histories of two Toll-like receptors (Tlr4andTlr7) in wild rodents (MURINAE).BMC Evol Biol. 2013, 13: 194-10.1186/1471-2148-13-194.PubMed CentralView ArticlePubMedGoogle Scholar
- Cormican P, Lloyd AT, Downing T, Connell SJ, Bradley D, O'Farrelly C: The avian Toll-Like receptor pathway-Subtle differences amidst general conformity. Dev Comp Immunol. 2009, 33: 967-973. 10.1016/j.dci.2009.04.001.View ArticlePubMedGoogle Scholar
- Alcaide M, Edwards SV: Molecular evolution of the Toll-like receptor multigene family in birds. Mol Biol Evol. 2011, 28: 1703-1715. 10.1093/molbev/msq351.View ArticlePubMedGoogle Scholar
- Grueber CE, Wallis GP, Jamieson IG: Episodic positive selection in the evolution of avian toll-like receptor innate immunity genes. PLoS One. 2014, 9: e89632-10.1371/journal.pone.0089632.PubMed CentralView ArticlePubMedGoogle Scholar
- Brownlie R, Allan B: Avian toll-like receptors. Cell Tissue Res. 2011, 343: 121-130. 10.1007/s00441-010-1026-0.View ArticlePubMedGoogle Scholar
- Iqbal M, Philbin VJ, Smith AL: Expression patterns of chicken Toll-like receptor mRNA in tissues, immune cell subsets and cell lines. Vet Immunol Immunopathol. 2005, 104: 117-127. 10.1016/j.vetimm.2004.11.003.View ArticlePubMedGoogle Scholar
- Leveque G, Forgetta V, Morroll S, Smith AL, Bumstead N, Barrow P, Loredo-Osti JC, Morgan K, Malo D: Allelic variation in TLR4 is linked to susceptibility to Salmonella enterica serovar typhimurium infection in chickens. Infect Immun. 2003, 71: 1116-1124. 10.1128/IAI.71.3.1116-1124.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Fukui A, Inoue N, Matsumoto M, Nomura M, Yamada K, Matsuda Y, Toyoshima K, Seya T: Molecular cloning and functional characterization of chicken Toll-like receptors - A single chicken toll covers multiple molecular patterns. J Biol Chem. 2001, 276: 47143-47149. 10.1074/jbc.M103902200.View ArticlePubMedGoogle Scholar
- Lynn DJ, Lloyd AT, O'Farrelly C: In silico identification of components of the Toll-like receptor (TLR) signaling pathway in clustered chicken expressed sequence tags (ESTs). Vet Immunol Immunopathol. 2003, 93: 177-184. 10.1016/S0165-2427(03)00058-8.View ArticlePubMedGoogle Scholar
- Philbin VJ, Iqbal M, Boyd Y, Goodchild MJ, Beal RK, Bumstead N, Young J, Smith AL: Identification and characterization of a functional, alternatively spliced Toll-like receptor 7 (TLR7) and genomic disruption of TLR8 in chickens. Immunology. 2005, 114: 507-521. 10.1111/j.1365-2567.2005.02125.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Yilmaz A, Shen SX, Adelson DL, Xavier S, Zhu JJ: Identification and sequence analysis of chicken Toll-like receptors. Immunogenetics. 2005, 56: 743-753. 10.1007/s00251-004-0740-8.View ArticlePubMedGoogle Scholar
- Smith J, Speed D, Law AS, Glass EJ, Burt DW: In-silico identification of chicken immune-related genes. Immunogenetics. 2004, 56: 122-133. 10.1007/s00251-004-0669-y.View ArticlePubMedGoogle Scholar
- Vinkler M, Bryjova A, Albrecht T, Bryja J: Identification of the first Toll-like receptor gene in passerine birds: TLR4 orthologue in zebra finch (Taeniopygia guttata). Tissue Antigens. 2009, 74: 32-41. 10.1111/j.1399-0039.2009.01273.x.View ArticlePubMedGoogle Scholar
- Gopinath VP, Biswas M, Raj GD, Raja A, Kumanan K, Elankumaran S: Molecular cloning and tissue-specific expression of Toll-like receptor 5 gene from turkeys. Avian Dis. 2011, 55: 480-485. 10.1637/9590-102710-ResNote.1.View ArticlePubMedGoogle Scholar
- Elfeil WK, Abouelmaatti RR, Sun CJ, Han WY, Li XK, Ma JS, Lei LC, Liu SS, Yang YJ, Wang Y, Mandour M, Fawzy M, Shalaby MN: Identification, cloning, expression of a novel functional Anas platyrhynchos mRNA TLR4. J Anim Vet Adv. 2012, 11: 1727-1733. 10.3923/javaa.2012.1727.1733.View ArticleGoogle Scholar
- MacDonald MRW, Xia JG, Smith AL, Magor KE: The duck Toll like receptor 7: Genomic organization, expression and function. Mol Immunol. 2008, 45: 2055-2061. 10.1016/j.molimm.2007.10.018.View ArticlePubMedGoogle Scholar
- Jia H, Li G, Li J, Tian Y, Wang D, Shen J, Tao Z, Xu J, Lu L: Cloning, expression and bioinformatics analysis of the duck TLR 4 gene. Br Poult Sci. 2012, 53: 190-197. 10.1080/00071668.2012.674208.View ArticlePubMedGoogle Scholar
- Wang F, Lu LZ, Yuan H, Tian Y, Li JJ, Shen JD, Tao ZR, Fu Y: cDNA cloning, characterization and expression analysis of Toll-like receptor 4 gene in goose. Can J Anim Sci. 2011, 91: 371-377. 10.4141/cjas2011-002.View ArticleGoogle Scholar
- Ruan WK, Wu YH, Zheng SJJ: Different genetic patterns in avian Toll-like receptor (TLR)5 genes. Mol Biol Rep. 2012, 39: 3419-3426. 10.1007/s11033-011-1113-7.View ArticlePubMedGoogle Scholar
- Ramasamy KT, Verma P, Reddy MR, Murugesan S: Molecular characterization of coding sequence and mRNA expression pattern of Toll-like receptor 15 in Japanese quail (Coturnix japonica) and indigenous chicken breeds (Aseel and Kadaknath). J Poult Sci. 2011, 48: 168-175. 10.2141/jpsa.011008.View ArticleGoogle Scholar
- Uno Y, Usui T, Fujimoto Y, Ito T, Yamaguchi T: Quantification of interferon, interleukin, and Toll-like receptor 7 mRNA in quail splenocytes using real-time PCR. Poult Sci. 2012, 91: 2496-2501. 10.3382/ps.2012-02283.View ArticlePubMedGoogle Scholar
- Downing T, Lloyd AT, O'Farrelly C, Bradley DG: The differential evolutionary dynamics of avian cytokine and TLR gene classes. J Immunol. 2010, 184: 6993-7000. 10.4049/jimmunol.0903092.View ArticlePubMedGoogle Scholar
- Keestra AM, van Putten JPM: Unique properties of the chicken TLR4/MD-2 complex: Selective lipopolysaccharide activation of the MyD88-dependent pathway. J Immunol. 2008, 181: 4354-4362. 10.4049/jimmunol.181.6.4354.View ArticlePubMedGoogle Scholar
- Keestra AM, de Zoete MR, van Aubel RAMH, van Putten JPM: Functional characterization of chicken TLR5 reveals species-specific recognition of flagellin. Mol Immunol. 2008, 45: 1298-1307. 10.1016/j.molimm.2007.09.013.View ArticlePubMedGoogle Scholar
- Brownlie R, Zhu JZ, Allan B, Mutwiri GK, Babiuk LA, Potter A, Griebel P: Chicken TLR21 acts as a functional homologue to mammalian TLR9 in the recognition of CpG oligodeoxynucleotides. Mol Immunol. 2009, 46: 3163-3170. 10.1016/j.molimm.2009.06.002.View ArticlePubMedGoogle Scholar
- DeMarco ML, Woods RJ: From agonist to antagonist: Structure and dynamics of innate immune glycoprotein MD-2 upon recognition of variably acylated bacterial endotoxins. Mol Immunol. 2011, 49: 124-133. 10.1016/j.molimm.2011.08.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Marr N, Novikov A, Hajjar AM, Caroff M, Fernandez RC: Variability in the lipooligosaccharide structure and endotoxicity among Bordetella pertussis strains. J Infect Dis. 2010, 202: 1897-1906. 10.1086/657409.View ArticlePubMedGoogle Scholar
- Andersen-Nissen E, Smith KD, Strobe KL, Barrett SLR, Cookson BT, Logan SM, Aderem A: Evasion of Toll-like receptor 5 by flagellated bacteria. Proc Natl Acad Sci U S A. 2005, 102: 9247-9252. 10.1073/pnas.0502040102.PubMed CentralView ArticlePubMedGoogle Scholar
- Walsh C, Gangloff M, Monie T, Smyth T, Wei B, McKinley TJ, Maskell D, Gay N, Bryant C: Elucidation of the MD-2/TLR4 interface required for signaling by lipid IVa. J Immunol. 2008, 181: 1245-1254. 10.4049/jimmunol.181.2.1245.View ArticlePubMedGoogle Scholar
- Hajjar AM, Ernst RK, Tsai JH, Wilson CB, Miller SI: Human Toll-like receptor 4 recognizes host-specific LPS modifications. Nat Immunol. 2002, 3: 354-359. 10.1038/ni777.View ArticlePubMedGoogle Scholar
- Lien E, Means TK, Heine H, Yoshimura A, Kusumoto S, Fukase K, Fenton MJ, Oikawa M, Qureshi N, Monks B, Finberg RW, Ingalls RR, Golenbock DT: Toll-like receptor 4 imparts ligand-specific recognition of bacterial lipopolysaccharide. J Clin Invest. 2000, 105: 497-504. 10.1172/JCI8541.PubMed CentralView ArticlePubMedGoogle Scholar
- Poltorak A, Ricciardi-Castagnoli P, Citterio S, Beutler B: Physical contact between lipopolysaccharide and Toll-like receptor 4 revealed by genetic complementation. Proc Natl Acad Sci U S A. 2000, 97: 2163-2167. 10.1073/pnas.040565397.PubMed CentralView ArticlePubMedGoogle Scholar
- Calenge F, Kaiser P, Vignal A, Beaumont C: Genetic control of resistance to salmonellosis and to Salmonella carrier-state in fowl: a review. Genet Sel Evol. 2010, 42: 11-10.1186/1297-9686-42-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik SG, Lee HY, Lee JO: Crystal structure of the TLR1-TLR2 heterodimer induced by binding of a tri-acylated lipopeptide. Cell. 2007, 130: 1071-1082. 10.1016/j.cell.2007.09.008.View ArticlePubMedGoogle Scholar
- Kang JY, Nan X, Jin MS, Youn SJ, Ryu YH, Mah S, Han SH, Lee H, Paik SG, Lee JO: Recognition of lipopeptide patterns by Toll-like receptor 2-Toll-like receptor 6 heterodimer. Immunity. 2009, 31: 873-884. 10.1016/j.immuni.2009.09.018.View ArticlePubMedGoogle Scholar
- Bell JK, Botos I, Hall PR, Askins J, Shiloach J, Segal DM, Davies DR: The molecular structure of the Toll-like receptor 3 ligand-binding domain. Proc Natl Acad Sci U S A. 2005, 102: 10976-10980. 10.1073/pnas.0505077102.PubMed CentralView ArticlePubMedGoogle Scholar
- Choe J, Kelker MS, Wilson IA: Crystal structure of human Toll-like receptor 3 (TLR3) ectodomain. Science. 2005, 309: 581-585. 10.1126/science.1115253.View ArticlePubMedGoogle Scholar
- Liu L, Botos I, Wang Y, Leonard JN, Shiloach J, Segal DM, Davies DR: Structural basis of Toll-like receptor 3 signaling with double-stranded RNA. Science. 2008, 320: 379-381. 10.1126/science.1155406.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim HM, Park BS, Kim JI, Kim SE, Lee J, Oh SC, Enkhbayar P, Matsushima N, Lee H, Yoo OJ, Lee JO: Crystal structure of the TLR4-MD-2 complex with bound endotoxin antagonist eritoran. Cell. 2007, 130: 906-917. 10.1016/j.cell.2007.08.002.View ArticlePubMedGoogle Scholar
- Park BS, Song DH, Kim HM, Choi BS, Lee H, Lee JO: The structural basis of lipopolysaccharide recognition by the TLR4-MD-2 complex. Nature. 2009, 458: 1191-1195. 10.1038/nature07830.View ArticlePubMedGoogle Scholar
- Ohto U, Fukase K, Miyake K, Shimizu T: Structural basis of species-specific endotoxin sensing by innate immune receptor TLR4/MD-2. Proc Natl Acad Sci U S A. 2012, 109: 7421-7426. 10.1073/pnas.1201193109.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu YW, Tao X, Shen BH, Horng T, Medzhitov R, Manley JL, Tong L: Structural basis for signal transduction by the Toll/interleukin-1 receptor domains. Nature. 2000, 408: 111-115. 10.1038/35047056.View ArticlePubMedGoogle Scholar
- Nyman T, Stenmark P, Flodin S, Johansson I, Hammarstrom M, Nordlund P: The crystal structure of the human Toll-like receptor 10 cytoplasmic domain reveals a putative signaling dimer. J Biol Chem. 2008, 283: 11861-11865. 10.1074/jbc.C800001200.View ArticlePubMedGoogle Scholar
- Mizel SB, West AP, Hantgan RR: Identification of a sequence in human Toll-like receptor 5 required for the binding of Gram-negative flagellin. J Biol Chem. 2003, 278: 23624-23629. 10.1074/jbc.M303481200.View ArticlePubMedGoogle Scholar
- Andersen-Nissen E, Smith KD, Bonneau R, Strong RK, Aderem A: A conserved surface on Toll-like receptor 5 recognizes bacterial flagellin. J Exp Med. 2007, 204: 393-403. 10.1084/jem.20061400.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei TD, Gong J, Rossle SC, Jamitzky F, Heckl WM, Stark RW: A leucine-rich repeat assembly approach for homology modeling of the human TLR5-10 and mouse TLR11-13 ectodomains. J Mol Model. 2011, 17: 27-36. 10.1007/s00894-010-0697-5.View ArticlePubMedGoogle Scholar
- Wei TD, Gong J, Jamitzky F, Heckl WM, Stark RW, Rossle SC: Homology modeling of human Toll-like receptors TLR7, 8, and 9 ligand-binding domains. Protein Sci. 2009, 18: 1684-1691. 10.1002/pro.186.PubMed CentralView ArticlePubMedGoogle Scholar
- Kubarenko AV, Ranjan S, Colak E, George J, Frank M, Weber ANR: Comprehensive modeling and functional analysis of Toll-like receptor ligand-recognition domains. Protein Sci. 2010, 19: 558-569.PubMed CentralPubMedGoogle Scholar
- Zhou KF, Kanai R, Lee P, Wang HW, Modis Y: Toll-like receptor 5 forms asymmetric dimers in the absence of flagellin. J Struct Biol. 2012, 177: 402-409. 10.1016/j.jsb.2011.12.002.View ArticlePubMedGoogle Scholar
- Yoon SI, Kurnasov O, Natarajan V, Hong MS, Gudkov AV, Osterman AL, Wilson IA: Structural basis of TLR5-flagellin recognition and signaling. Science. 2012, 335: 859-864. 10.1126/science.1215584.PubMed CentralView ArticlePubMedGoogle Scholar
- Hackett SJ, Kimball RT, Reddy S, Bowie RCK, Braun EL, Braun MJ, Chojnowski JL, Cox WA, Han KL, Harshman J, Huddleston CJ, Marks BD, Miglia KJ, Moore WS, Sheldon FH, Steadman DW, Witt CC, Yuri T: A phylogenomic study of birds reveals their evolutionary history. Science. 2008, 320: 1763-1768. 10.1126/science.1157704.View ArticlePubMedGoogle Scholar
- Holm L, Kaariainen S, Rosenstrom P, Schenkel A: Searching protein structure databases with DaliLite v. 3. Bioinformatics. 2008, 24: 2780-2781. 10.1093/bioinformatics/btn507.PubMed CentralView ArticlePubMedGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006, 34: W609-W612. 10.1093/nar/gkl315.PubMed CentralView ArticlePubMedGoogle Scholar
- Letunic I, Doerks T, Bork P: SMART 7: recent updates to the protein domain annotation resource. Nucleic Acids Res. 2012, 40: D302-D305. 10.1093/nar/gkr931.PubMed CentralView ArticlePubMedGoogle Scholar
- Offord V, Coffey TJ, Werling D: LRRfinder: A web application for the identification of leucine-rich repeats and an integrative Toll-like receptor database. Dev Comp Immunol. 2010, 34: 1035-1041. 10.1016/j.dci.2010.05.004.View ArticlePubMedGoogle Scholar
- Cserzo M, Eisenhaber F, Eisenhaber B, Simon I: TM or not TM: transmembrane protein prediction with low false positive rate using DAS-TMfilter. Bioinformatics. 2004, 20: 136-137. 10.1093/bioinformatics/btg394.View ArticlePubMedGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8: 785-786. 10.1038/nmeth.1701.View ArticlePubMedGoogle Scholar
- Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5: 725-738. 10.1038/nprot.2010.5.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008, 9: 40-10.1186/1471-2105-9-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Ewald SE, Lee BL, Lau L, Wickliffe KE, Shi G-P, Chapman HA, Barton GM: The ectodomain of Toll-like receptor 9 is cleaved to generate a functional receptor. Nature. 2008, 456: 658-662. 10.1038/nature07405.PubMed CentralView ArticlePubMedGoogle Scholar
- McGuffin LJ: The ModFOLD server for the quality assessment of protein structural models. Bioinformatics. 2008, 24: 586-587. 10.1093/bioinformatics/btn014.View ArticlePubMedGoogle Scholar
- Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA: PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 2004, 32: W665-W667. 10.1093/nar/gkh381.PubMed CentralView ArticlePubMedGoogle Scholar
- Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA: Electrostatics of nanosystems: Application to microtubules and the ribosome. Proc Natl Acad Sci U S A. 2001, 98: 10037-10041. 10.1073/pnas.181342398.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online - a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33: W557-W559. 10.1093/nar/gki352.PubMed CentralView ArticlePubMedGoogle Scholar
- Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO: The global diversity of birds in space and time. Nature. 2012, 491: 444-448. 10.1038/nature11631.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMedGoogle Scholar
- Cozzetto D, Kryshtafovych A, Fidelis K, Moult J, Rost B, Tramontano A: Evaluation of template-based models in CASP8 with standard measures. Proteins. 2009, 77: 18-28. 10.1002/prot.22561.PubMed CentralView ArticlePubMedGoogle Scholar
- Hasegawa H, Holm L: Advances and pitfalls of protein structural alignment. Curr Opin Struct Biol. 2009, 19: 341-348. 10.1016/j.sbi.2009.04.003.View ArticlePubMedGoogle Scholar
- Kalinowski ST: How well do evolutionary trees describe genetic relationships among populations?. Heredity. 2009, 102: 506-513. 10.1038/hdy.2008.136.View ArticlePubMedGoogle Scholar
- Yang Z: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Murrell B, Moola S, Mabona A, Weighill T, Sheward D, Kosakovsky Pond SL, Scheffler K: FUBAR: a fast, unconstrained bayesian approximation for inferring selection. Mol Biol Evol. 2013, 30: 1196-1205. 10.1093/molbev/mst030.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang ZH, Wong WSW, Nielsen R: Bayes empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.View ArticlePubMedGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21: 2531-2533. 10.1093/bioinformatics/bti320.View ArticleGoogle Scholar
- Atchley WR, Zhao JP, Fernandes AD, Druke T: Solving the protein sequence metric problem. Proc Natl Acad Sci U S A. 2005, 102: 6395-6400. 10.1073/pnas.0408677102.PubMed CentralView ArticlePubMedGoogle Scholar
- Ashkenazy H, Erez E, Martz E, Pupko T, Ben-Tal N: ConSurf 2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res. 2010, 38: W529-W533. 10.1093/nar/gkq399.PubMed CentralView ArticlePubMedGoogle Scholar
- Le SQ, Gascuel O: An improved general amino acid replacement matrix. Mol Biol Evol. 2008, 25: 1307-1320. 10.1093/molbev/msn067.View ArticlePubMedGoogle Scholar
- Eo SH, Bininda-Emonds ORP, Carroll JP: A phylogenetic supertree of the fowls (Galloanserae, Aves). Zoologica Scripta. 2009, 38: 465-481. 10.1111/j.1463-6409.2008.00382.x.View ArticleGoogle Scholar
- Resman N, Vasl J, Oblak A, Pristovsek P, Gioannini TL, Weiss JP, Jerala R: Essential roles of hydrophobic residues in both MD-2 and Toll-like receptor 4 in activation by endotoxin. J Biol Chem. 2009, 284: 15052-15060. 10.1074/jbc.M901429200.PubMed CentralView ArticlePubMedGoogle Scholar
- Smirnova I, Poltorak A, Chan EKL, McBride C, Beutler B: Phylogenetic variation and polymorphism at the Toll-like receptor 4 locus (TLR4). Genome Biol. 2000, 1: research002.1-research002.10. 10.1186/gb-2000-1-1-research002.View ArticleGoogle Scholar
- 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. Proc Natl Acad Sci U S A. 2003, 100: 10364-10369. 10.1073/pnas.1333957100.PubMed CentralView ArticlePubMedGoogle Scholar
- Wlasiuk G, Nachman MW: Adaptation and constraint at Toll-like receptors in primates. Mol Biol Evol. 2010, 27: 2172-2186. 10.1093/molbev/msq104.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith SA, Jann OC, Haig D, Russell GC, Werling D, Glass EJ, Emes RD: Adaptive evolution of Toll-like receptor 5 in domesticated mammals. BMC Evol Biol. 2012, 12: 122-10.1186/1471-2148-12-122.PubMed CentralView ArticlePubMedGoogle Scholar
- Arbour NC, Lorenz E, Schutte BC, Zabner J, Kline JN, Jones M, Frees K, Watt JL, Schwartz DA: TLR4 mutations are associated with endotoxin hyporesponsiveness in humans. Nat Genet. 2000, 25: 187-191. 10.1038/76048.View ArticlePubMedGoogle Scholar
- Diebold SS, Massacrier C, Akira S, Paturel C, Morel Y, Reis e Sousa C: Nucleic acid agonists for Toll-like receptor 7 are defined by the presence of uridine ribonucleotides. Eur J Immunol. 2006, 36: 3256-3267. 10.1002/eji.200636617.View ArticlePubMedGoogle Scholar
- Bainova H, Kralova T, Bryjova A, Albrecht T, Bryja J, Vinkler M: First evidence of independent pseudogenization of Toll-like receptor 5 in passerine birds. Dev Comp Immunol. 2014, 45: 151-155. 10.1016/j.dci.2014.02.010.View ArticlePubMedGoogle Scholar
- Woolhouse MEJ, Webster JP, Domingo E, Charlesworth B, Levin BR: Biological and biomedical implications of the co-evolution of pathogens and their hosts. Nat Genet. 2002, 32: 569-577. 10.1038/ng1202-569.View ArticlePubMedGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.