Molecular evolution of the vertebrate TLR1 gene family - a complex history of gene duplication, gene conversion, positive selection and co-evolution
© Huang et al; licensee BioMed Central Ltd. 2011
Received: 10 December 2010
Accepted: 28 May 2011
Published: 28 May 2011
The Toll-like receptors represent a large superfamily of type I transmembrane glycoproteins, some common to a wide range of species and others are more restricted in their distribution. Most members of the Toll-like receptor superfamily have few paralogues; the exception is the TLR1 gene family with four closely related genes in mammals TLR1, TLR2, TLR6 and TLR10, and four in birds TLR1A, TLR1B, TLR2A and TLR2B. These genes were previously thought to have arisen by a series of independent gene duplications. To understand the evolutionary pattern of the TLR1 gene family in vertebrates further, we cloned the sequences of TLR1A, TLR1B, TLR2A and TLR2B in duck and turkey, constructed phylogenetic trees, predicted codons under positive selection and identified co-evolutionary amino acid pairs within the TLR1 gene family using sequences from 4 birds, 28 mammals, an amphibian and a fish.
This detailed phylogenetic analysis not only clarifies the gene gains and losses within the TLR1 gene family of birds and mammals, but also defines orthologues between these vertebrates. In mammals, we predict amino acid sites under positive selection in TLR1, TLR2 and TLR6 but not TLR10. We detect co-evolution between amino acid residues in TLR2 and the other members of this gene family predicted to maintain their ability to form functional heterodimers. In birds, we predict positive selection in the TLR2A and TLR2B genes at functionally significant amino acid residues. We demonstrate that the TLR1 gene family has mostly been subject to purifying selection but has also responded to directional selection at a few sites, possibly in response to pathogen challenge.
Our phylogenetic and structural analyses of the vertebrate TLR1 family have clarified their evolutionary origins and predict amino acid residues likely to be important in the host's defense against invading pathogens.
KeywordsGene Duplication Gene Conversion Positive Selection Co-evolution TLR1 gene family
The innate immune response is the first line of defense against invading pathogens and is required for an efficient adaptive immune response. Toll-like receptors (TLRs) play a crucial role in innate immunity by recognizing specific pathogen-associated molecular patterns including lipoproteins, lipopeptides, lipopolysaccharide, flagellin, dsRNA, ssRNA and CpG DNA motifs (Akira and Takeda 2004; West et al. 2006). Phylogenetic analysis has classified vertebrate TLRs into six major gene families [1–3]: TLR1, TLR3, TLR4, TLR5, TLR7 and TLR11. Most vertebrates have a single gene for each TLR gene family [1, 2]; the major exception is the TLR1 gene family, which appears to have multiple, paralogous genes. In mammals, the TLR1 gene family consists of four members: TLR1, TLR2, TLR6 and TLR10. In birds, there are also four paralogues with multiple names in the literature, two TLR1-like genes: TLR1La (TLR1 type 1, TLR1.1, TLR1/6/10 and TLR16) and TLR1Lb (TLR1 type 2 and TLR1.2), and two TLR2-like genes: TLR2a (TLR2 type 1 and TLR2.1) and TLR2b (TLR2 type 2 and TLR2.2). More distantly related members of the TLR1 gene family are TLR14 and TLR15 found in mammals and birds, respectively. Phylogenetic and genome analyses suggest that TLR1- and TLR2-like genes diverged from a common ancestor early in vertebrate evolution [1, 2]. The genes of the TLR1 gene family are usually in tandem in both avian and mammalian genomes and are likely to be the product of successive rounds of tandem gene duplications from an ancestral gene.
Gene duplication followed by sequence divergence is one of the most important mechanisms for generating new genes with novel functions [4, 5]. Genetic drift and positive selection have both played a role in the fixation and early evolution of duplicated genes [6–9]. Many studies have classified genes showing evidence of positive selection , where genes with roles in host defense/immunity, chemosensory perception and reproduction are particularly over represented. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving region in the TLR-related genes in primates[11–13]. A recent comparative sequence analysis from a panel of domestic chickens and wild jungle fowl suggest that TLRs show positive selection . More specifically, an earlier study suggested that the leucine-rich repeats of chicken TLR1La and TLR1Lb may have been subject to positive selection. In contrast, the sequences of all other avian TLRs were under purifying selection , which appears to be the trend for all TLRs in mammals .
In mammals, TLR2 dimerizes with either TLR1 to recognize microbial triacyl lipoproteins or TLR6 to recognize diacyl lipopeptides found in Mycoplasma, lipoteichoic acid of Gram-positive bacteria or Zymosan of yeasts [16, 17]. Recently, human TLR10 has also been shown to dimerize with TLR2 and to recognize triacyl lipoproteins , although signaling by this complex appears to differ from the other TLR1/6/2 complexes. In chicken [19, 20] expression of TLR1La or b, or TLR2a or b alone fail to be activated by exogenous agonists; however heterodimers of all four TLR1/TLR2 combinations were activated by both diacyl (Malp-2) and triacyl (Pam3) lipopeptides, with the exception of TLR2a/TLR1Lb, which was activated by Pam3 but not Malp-2. In addition, TLR2a/TLRL1b was also activated by peptidoglycan. The interaction of these TLR1- and TLR2-like proteins in the absence of any agonist was also shown by immunoprecipitation. Thus, chicken TLR1-like proteins interact with TLR2-like proteins and recognize agonists identical to those used in mammals by heterodimers between TLR2 and TLR1, 6 or 10.
A recent study  showed sequence divergence in the C-terminal regions of TLR1 and TLR6 in mammals to be severely limited by gene conversion. This 300 amino acid region encodes the LRR16-19 motifs, the C-terminal cap motif, the transmembrane domain and most of the intracellular TIR domain. The need to co-evolve with binding partners, such as MyD88 and TIRAP may explain the selective advantage for gene conversion in the TIR domains of TLR1 and TLR6. A recent report also found evidence for gene conversion in the N- and C-terminal regions of the duplicated TLR2A and TLR2B genes found in the genomes of chicken and zebra finch . Thus, gene conversion may be a general mechanism that limits the divergence of critical regions in these proteins, while allowing other domains, such as the ligand binding domain, to evolve diversified functions.
The aim of this work was to seek possible explanations for the apparent multiple, independent gene duplications proposed to account for the diversity and size of the TLR1 gene family . From a phylogenetic analysis of many avian and mammalian TLR1-like and TLR2-like genes we have defined the likely role of positive selection, gene duplication, gene conversion and gene co-evolution in the evolutionary dynamics of this family. Our phylogenetic and structural analysis of vertebrate TLR1 gene family members also predicts amino acid residues likely to be important in the primary function of TLRs in the host's defense against invading pathogens.
Nucleic acid extraction
Total RNA from turkey spleen was extracted by TRIZOL (Invitrogen, Paisley, UK) according to the manufacture's protocol. The sample was re-suspended in a final volume of 50 μl of RNAse-free water. A routine phenol/chloroform extraction method was used to purify turkey genomic DNA. Both the quantity and quality of total RNA and genomic DNA were assessed at OD260 and OD280 using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, USA).
DNA cloning and sequencing
BAC clones containing TLR1A/TLR1B or TLR2A/TLR2B sequence were screened by 4D, two-step PCR from a HindIII, female Beijing duck genomic BAC library with inserts from vector pIndig-5 (Epi-Center, USA). The primers used to screen the duck BAC library, were designed from the orthologous sequences in other species such as chicken, pig, human or the partial sequence of duck (AY838880) (Additional file 1, Table S1) . Positive clones were confirmed by sequence analysis of PCR products, then subcloned into pUC118 and sequenced by the shotgun method, (with an average of 6.38Х coverage) using Big Dye Terminator V3.1 (ABI, Foster city, USA) on an ABI3130X Sequencer (ABI, Columbia, MD). Base calling and quality assessment were performed using Phred [25, 26], Assembly was performed using Phred and Consed [27, 28]. Sequence identity between duck and chicken genomic sequences was analyzed using zPicture .
To clone the TLR1 gene family sequences from turkey, PCR reactions were carried out in 100ųl total volume, containing 100 ng turkey genomic DNA, 50 mM Tris/HCl (pH8.3), 10 mM KCl, 50 mM (NH4)2SO4, 2.0 mM MgCl2, 200 mM dNTP, 10 U FastStart Taq DNA polymerase (Roche, Mannheim, Germany), and 100 pmol each primer designed from the orthologous sequence in chicken and duck (Additional file 1, Table S1). PCR reaction conditions were: denaturing for 5 min at 94°C, followed by 94°C for 30s, annealing for 30s, and 72°C for 1 min per 1 kb, 35 cycles, with a final 7 min elongation step at 72°C. PCR products were purified from agarose gels with the QIAEXII gel extraction kit (Qiagen, Hilden, Germany), ligated into pCR4-TOPO (Invitrogen, Carlsbad, USA), transformed into E. coli TOP10 competent cells and plasmid DNA was sequenced using T3 or T7 universal primer. In addition, the sequences at the 3`end of TLR2A and TLR2B were cloned from turkey spleen RNA by 3`-RACE using the FirstChoice RLM-RACE kit (Ambion, Austin, USA) and confirmed against the genomic sequence of turkey.
Sequences from the TLR1 gene family were collected from the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov) using the Blast search program , retrieved through searches in ENSEMBL (http://www.ensembl.org/index.html) and identified by Blat searches on the UCSC Genome Bioinformatics database (http://genome.ucsc.edu/). The alignments of amino acid and nucleotide sequence used for the analysis of gene conversion and the construction of phylogenetic trees, and for the calculation of percentage identities were made using ClustalX (version 2.0) . Sequence homologies were displayed using JalView 2.5.1.
Gene conversion analysis
Recombination analysis was performed using GENECONV (version 1.8) with default settings . The more conservative and accurate p-values from permutations (10,000 pseudo-replicates) were used and p-values (p < 0.05) from global fragments were corrected for multiple comparisons. The aligned sequences were examined for possible gene conversion events using the sliding-window genetic diversity plot (Simplot) (version 3.5.1) by bootscan with a neighbor-joining (NJ) tree, maximum-likelihood distance model, with 1,000 pseudo-replicates in a 200 bp window and with a step size of 20 bp .
Molecular phylogenies based on sequence alignments are only as accurate as the alignment data from which they are produced; consequently it is important that the alignment data quality is determined . In order to assess the alignment's tree like structure, likelihood mapping was carried out using Tree-Puzzle (version 5.2) , which generally showed alignments to have strong phylogenetic signals with 89-94% of quartets (1,000 pseudo-replicates) support a unique phylogeny (data not shown).
Phylogenetic relationships among orthologues and paralogues of TLR1 and TLR2 subfamilies were carried out on the protein alignment using PHYML (version 2.4.4) . After selection of the best substitution model for the aligned protein or DNA sequence by comparing the likelihood values of different models, two rounds of analysis were used to search the maximum-likelihood trees under the JTT model of molecular evolution with 4 substitution rates classes for protein sequences or HKY85 for DNA sequences. The phylogenetic tree, the proportion of invariable sites and Gamma distribution for protein sequence or transition/transversion ratio for DNA sequence parameters were estimated under the above model in the first round. Then, the robustness of the inferred tree was assessed using bootstrapping (1,000 pseudo-replicates) by fixing the proportion of invariable sites and Gamma distribution or transition/transversion parameters in the second round. Gene Trees were displayed using FigTree (version 1.3.1) (http://tree.bio.ed.ac.uk/).
Estimation of the time of gene duplication events
The fossil calibration time of the orthologous TLRs was inferred from species divergence times obtained from the TIMETREE web site (http://www.timetree.org/index.php). Divergence times of gene conversion and duplication were estimated under global and local clock models with two runs using CODEML or BASEML within the PAML package (version 4) . Firstly, the "kappa" (the transition/transversion rate ratio) and "omega" (the non-synonymous/synonymous rate ratio) for CODEML or "alpha" (gamma-distributed rate) for BASEML were estimated, and then the divergence times were estimated with fixed values for "kappa" and "omega" for CODEML or "alpha" for BASEML.
Maximum likelihood tests of positive selection
Multiple sequence alignment of coding sequences were created using DAMBE , first on the translated DNA sequences and then the nucleotide sequence of codon triplets was recovered by back translation. DAMBE was also used to check for saturation of nucleotide substitutions using a plot of the number of transitions and transversions for each pairwise comparison against the genetic distance calculated with the F84 model of nucleotide substitution , which allows different equilibrium nucleotide frequencies and a transition rate-transversion rate bias. Multiple sequence alignments not showing saturation of nucleotide substitutions were then analysed further using CODEML. A series of nested likelihood ratio tests (LRTs) were performed using CODEML to investigate whether some sites were under positive selection in each TLR group. The first test compared M1A (nearly neutral: p0, p1, ω0 < 1, ω1 = 1, NS sites = 1) against M2A (positive selection: p0, p1, p2, ω0 < 1, ω1 = 1, ω2 < 1, NS sites = 2). The second test compared M7 (beta: p, q, NS sites = 7) with M8 (beta & ω: p0, p1, p, q, ωs > 1, NS sites = 8). The third test was between M8A (beta & ωs = 1: fix omega = 1, omega = 1, NS sites = 8) and M8.
Intra- and intermolecular co-evolution analysis
To identify co-evolutionary patterns, we used the parametric method based on correlated evolutionary patterns among amino-acid sites  implemented in CAPS (Version 1.0) . The probabilities and significance of the correlated evolutionary patterns among amino-acid sites were estimated using a large number (10,000) of random samplings and a small α value (0.001) to minimize a false-positive rate (type I error). CAPS implements the step-down permutation procedure  to correct for multiple testing. The scores for the amino-acid substitutions were obtained using the appropriate blocks substitution matrix (BLOSUM80) , depending on the similarity of the protein sequences. Possible interactions were also filtered based on cellular compartments: extracellular, transmembrane and cytoplasmic. All amino-acid sites reported in the co-evolutionary analyses refer to the positions in the protein sequences of human TLR1, TLR2, TLR6 and TLR10.
Visualization of co-evolutionary networks
Cytoscape (Version 2.7.0)  was used to visualize the co-evolutionary networks identified by CAPS, and to generate the networks of correlation between co-evolving amino acids. The correlation coefficients generated in CAPS were used to determine the color of the lines between nodes (amino-acid residues).
Protein structure analysis
The Uniprot protein database (http://www.uniprot.org) was used as an initial resource to define the protein domain locations of specific amino acid residues in TLR proteins. If there were no annotated protein domains, these were determined by aligning genes with their human orthologues, and then inferring domains from Uniprot. To examine potential functional or structural significance of specific amino acid residues, their location was mapped onto the known or predicted three-dimensional structures of TLR proteins. The known structures of human TLR1/2/6 (Jin et al. 2007) were used directly and as templates to predict the structures of human TLR10 using CPHmodels (version 3.0), a protein homology modeling server . Template recognition was based on profile-profile alignment guided by secondary structure and exposure predictions. The atomic coordinates for the PDB models of human TLR1 and TLR2 proteins suggested the secondary structure of human TLR10 was similar to human TLR1 (E value 7.0 × 10-89). Amino acid residues were mapped onto these structures using the protein BLAST alignment of each group, and a PyMOL script file was generated for visualization using PyMOL (version 1.1) (http://www.pymol.org).
Isolation of Toll-like receptor-like 1 and 2 genes from the turkey and duck genomes
The coding sequences for chicken and zebra finch members of the TLR1 gene family were extracted from sequence databases. The coding sequences for turkey TLR1-like and TLR2-like genes were cloned from genomic DNA by PCR using primers based on the chicken orthologues. The 3`-end mRNA sequences of turkey TLR2A and TLR2B were cloned by 3`-RACE with the primers shown in Additional file 1, Table S1. Duck BAC clones containing genes for TLR1A, TLR1B, TLR2A and TLR2B were isolated by four-dimensional, two-step PCR and sequenced (see Materials and Methods). A 67,954 bp BAC clone (DBS1203G04, GenBank accession FJ477859) contained both TLR1A and TLR1B genes, and a 75,080 bp BAC clone (DBS1405N01, GenBank accession FJ477862) contained TLR2A and TLR2B (Additional file 2). The coding sequences of the TLR1 gene family in duck and turkey are encoded by single exons, as found in the chicken. Sequence alignment revealed the average percentage identity for pair-wise comparisons between the orthologues in chicken/duck to be 85.3 ± 1.5% for nucleotide and 79.8 ± 2.5% for amino acid sequences, and in chicken/turkey to be 93.0 ± 1.2% for nucleotide and 90.3 ± 2.2% for protein sequences. Sequences for the duck and turkey TLR1 genes were submitted to GenBank (Accessions FJ477857-FJ477862).
Multiple gene conversion events during the evolution of the TLR1 gene family
Gene conversion analysis for the TLR1 gene family using coding sequences from birds and mammals
BC KA P-value
Number of Polymorphisms
Number of Differences
We divided the sequences of groups 1, 2 and 3 into regions predicted to have undergone gene conversion or not (Table 1). Amino acid sequences encoded by these regions were used to create separate ML phylogenies using PHYML  for group 1, 2 or 3 genes. A TLR1 or TLR2 orthologue from zebrafish (Danio rerio) was used as an outgroup gene to root these trees. The ML trees of group 1, 2 and 3 genes calculated using the gene conversion free regions (Additional file 8, Figure S6A, Additional file 9, Figure S7B, Additional file 10, Figure S8A) were as expected from the species tree. In contrast, using group 1, 2 and 3 sequences in the regions predicted to be susceptible to gene conversion, showed paralogues to be more similar to each other than the corresponding orthologues (Additional file 8, Figure S6B, Additional file 9, Figure S7A, S7C). For most species, except primates, group 3 paralogues were more similar than orthologues and clustered together in TLR1/6 pairs within each species (Additional file 10, Figure S8B). However, for most primates, except the common marmoset (Callithrix Jacchus), we find orthologues to be more similar, indicating the last gene conversion event occurred prior to the Callithric Jacchus/Macaca mulatta split, 45 million years ago (Mya).
The sequence identities in specific regions between orthologues and paralogues are summarized in Additional file 1, Table S2. For groups 1 and 3 sequence identities in the C-terminal regions of paralogues (within the same species) were significantly greater than the corresponding identities between orthologues (p < 0.05). In contrast, sequence identities between the N-terminal regions of paralogues were lower than orthologues (p < 0.05). For group 2 sequence identities of N- and C-terminal regions of paralogues were higher than orthologues (p < 0.05). However, the central region was more conserved between orthologues (p < 0.05).
Amino acid sequences that remained invariant throughout evolution or became identical through convergent evolution are less likely to be encoded by the same codons than sequences that became identical after gene conversion. To discriminate between these possibilities as the cause for the greater sequence identity between paralogues, we compared the codon usage of conserved (and identical) amino acids in group 1, 2 and 3 genes. Overall, paralogues share 71.5%, 82.5% and 62.0% amino acid identity in groups 1, 2 and 3, respectively (Additional file 1, Table S2). Analysis of codon usage for conserved amino acids (Additional file 1, Table S3) in groups 1 and 3 indicated that identical codons were more frequently used in the C-terminal (96-99%) than in the N-terminal (57-66%) region. In group 2 identical codons in both N- and C-terminal regions were used more frequently (99%) than those in the central region (56%).
Taken together these results support a process of gene conversion rather than gene convergence, as the major cause for increased sequence identity between paralogues.
Identification of TLR2 pseudogenes in mammalian genomes
Two functional copies of TLR2 (TLR2A and TLR2B) were found in the genomes of chicken, turkey, duck and zebra finch, extending the earlier observations . To examine the evolutionary history of the TLR2 gene duplication, we searched all available vertebrate genome sequences for possible functional or non-functional sequences. Two functional copies of TLR2 were also detected in Xenopus tropicalis but only one functional gene was identified in 5 species of fish and 14 species of mammals. Previous work  detected a TLR2 pseudogene (TLR2P) in tandem with the functional copy in the genomes of opossum, dog and human. We have extended this observation to include an additional five mammalian genomes: chimpanzee, orangutan, rhesus monkey, marmoset and horse (Additional file 4). The identification of multiple TLR2 genes and pseudogenes on many mammals and Xenopus tropicalis suggests an early gene duplication event in the history of the TLR2 gene family. In most cases, mammalian TLR2 functional genes and pseudogenes formed distinct clusters (Additional file 11) as expected from the species tree. The exception was the horse, where its functional gene was more similar to its pseudogene, suggesting a recent gene conversion event between these paralogues. This was supported by finding high sequence homology between these paralogues using SIMPLOT (Additional file 6, Figure S4D).
A complex history of gene duplication and gene conversion during the evolution of the vertebrate TLR1 gene family
The dates for gene duplication were estimated from an analysis of the sequences free of gene conversion using CODEML  and summarized in Additional file 1, Table S4. Either global or local clock models of nucleotide substitution were used and calibrated using dates from the fossil record. For the TLR1 subfamily, we estimate an initial gene duplication event 359 Mya that gave rise to the TLR1A-TLR10 and TLR1B-TLR1/6 lineages, which was prior to the common ancestor for birds and mammals 310-325 Mya. Subsequently a further gene duplication 270-282 Mya gave rise to the TLR1 and TLR6 genes in mammals. The gene duplication event that gave rise to the TLR2 subfamily was estimated to be 341-356 Mya, which was prior to the common ancestor for birds and mammals 310-325 Mya. Similarly, the dates for gene conversion were inferred from an analysis of sequences predicted to have undergone gene conversion. Predicted dates were very variable between all species, varying from 0-47 Mya and generally older in mammals (Additional file 1, Table S4).
Positive selection in the coding regions of the TLR1 gene family in birds and mammals
The estimation of substitution saturation from either the TLR1 or TLR2 subfamilies suggested the number of transitions reached a plateau when d > 0.7 (higher than the number of 160 Mya). Beyond this time, positions are saturated for pairs of nucleotide sequences (Additional file 13). Based on these comparisons we then used site-specific CODEML analysis to detect evidence of positive selection on the full length sequences of TLR1 gene family orthologues[38, 46]: four groups of mammals (A: 12 TLR1, B: 12 TLR6, C: 15 TLR10 and D: 27 TLR2 species) and four birds (E: 4 TLR1A, F: 4 TLR1B, G: 4 TLR2A, H: 4 TLR2B species).
Positive selection in the TLR1 gene family in birds and mammals
No. of sequences
Proportion of sites ω> 1
Positively selected sites2
M1A vs. M2A
ω= 3.46 (p = 0.02)
M7 vs. M8
ω= 2.60 (p = 0.04)
174T, 238A, 286Q*, 599S**
M8A vs. M8
ω= 2.84 (p = 0.02)
M1A vs. M2A
ω = 2.84 (p = 0.02)
M7 vs. M8
ω = 2.02 (p = 0.04)
M8A vs. M8
M1A vs. M2A
M7 vs. M8
ω = 2.47 (p = 0.02)
M8A vs. M8
M1A vs. M2A
M7 vs. M8
ω = 2.07 (p = 0.02)
M8A vs. M8
M1A vs. M2A
M7 vs. M8
M8A vs. M8
M1A vs. M2A
M7 vs. M8
M8A vs. M8
M1A vs. M2A
ω = 2.85 (p = 0.04)
M7 vs. M8
ω = 2.19 (p = 0.09)
184T, 304N, 506S
M8A vs. M8
M1A vs. M2A
ω = 2.00 (p = 0.11)
M7 vs. M8
ω = 1.83 (p = 0.14)
175T, 260V, 496S
M8A vs. M8
Avian TLR1A and TLR1B parameter estimates under M2A or M8 models did not suggest any sites to be under positive selection. The lack of evidence for positive selection in the avian TLR1A/TLR1B genes may have been due to lack of power with only four sequences to compare, however this was not the case for the four avian TLR2A/TLR2B sequences. TLR2A estimates under M2A suggest 4% of sites to be under positive selection with ŵ2 = 2.9. Estimates under M8 suggest 9% of sites to be under positive selection with ŵs = 2.2. Three sites were inferred to have ω > 1 (p > 0.9) under M8: residues 184T, 304N and 506S. TLR2B estimates under M2A suggest 11% of sites to be under positive selection with ŵ2 = 2.0. Estimates under M8 suggest 14% of sites to be under positive selection with ŵs = 1.8. Three sites were inferred to have ω > 1 (p > 0.9) under M8: residues 175T, 260V and 496S.
Co-evolution of heterodimers between TLR1/6/10 and TLR2 in mammals
The CAPS analysis between 11 pairs of mammalian TLR1/2 sequences identified 28 groups of co-evolving amino acid residues (Additional file 14). Within TLR1, 35 residues were identified as co-evolving with 141 in TLR2 with correlation coefficients between 0.50 and 0.82 (P < 0.001). Most co-evolving residues were located in the extracellular region (33 in TLR1 and 130 in TLR2) and far fewer (2 in TLR1 and 8 in TLR2) in cytoplasmic regions. Similarly, we identified 63 groups of co-evolving amino acid residues using 12 pairs of TLR6/2 sequences (Additional file 15). Within TLR6, 76 amino acid residues were identified as co-evolving with 300 amino acid residues within TLR2 with correlation coefficients between 0.50 and 0.88 (P < 0.001). All of the inter-co-evolving sites of TLR1 and TLR6 were distributed in extracellular regions. Finally, we identified 76 groups of co-evolving amino acid residues using 12 pairs of TLR10/2 sequences (Additional file 16). Within TLR10, 91 amino acid residues were identified as co-evolving with 146 amino acid residues within TLR2 with correlation coefficients between 0.50 and 0.95 (P < 0.001). Most co-evolving residues within TLR10 were located in the extracellular region (56 in TLR10 and 102 in TLR2). However, inter-co-evolving sites were also detected in transmembrane (1 in TLR10 and 1 in TLR2) and cytoplasmic (2 in TLR10 and 4 in TLR2) regions.
We also looked for evidence of intra-molecular co-evolution using intra-protein CAPS analysis with the sequences of TLR1 family genes. Within TLR1, 111 amino acids from the extracellular region and 12 from the cytoplasmic region were identified as co-evolving with correlation coefficients between 0.72 and 1.00 (p < 0.001). However, there was no evidence of co-evolving amino acids from an analysis of 27 mammalian TLR2 sequences (p < 0.001). Within TLR6, 130 amino acids from the extracellular region and 14 from the cytoplasmic region which were identified as co-evolving with correlation coefficients between 0.61 and 0.99 (p < 0.001). Within TLR10, 29 amino acids from the extracellular region, 1 from the transmembrane domain and 6 from the cytoplasmic region were identified as co-evolving with correlation coefficients between 0.75 and 1.00 (p < 0.001).
TLRs play a crucial role in the host's early defense against invading pathogens [47, 48]. In the past few years our understanding of the mechanisms by which TLRs recognize pathogens and induce distinct signaling pathways has increased significantly, particularly in mammals [49, 50]. The tertiary structures of various TLRs and ligand complexes have also been determined [16, 51]. In contrast, studies on avian TLRs have been very limited[2, 15, 19, 20, 52–54]. In this paper, we examined the evolution and selective constraints on members of the TLR1 gene family in both birds and mammals. Here we discuss the details of these evolutionary origins and the significance of selection constraints on the structure and function of this family of cell surface receptors.
Evolutionary origins of the TLR1 gene family
Previous phylogenetic analyses of the TLR1 gene family suggested there had been many independent gene duplications in birds and mammals [1, 2, 23]. This contrasts with the TLR3, TLR4, TLR5 and TLR7 families where gene trees generally recapitulate species trees . Detailed analysis presented in this paper, together with the previous studies by Kruithof et al. and Cormican et al, [21, 22] suggests gene conversion events between paralogues had occurred during the evolution of the TLR1 gene family in both birds and mammals, giving rise to ambiguous orthologous relationships. We show that gene conversion events have occurred independently in most avian and mammalian lineages during the last 45 million years. Generally gene conversion events have been more recent in birds (0-20 Mya) than mammals (5-45 Mya). Gene conversion occurred in the mammalian TLR1/TLR6 genes 42-44 Mya in the common ancestor of four primate species (human, chimpanzee, orangutan and rhesus monkey) and independently in other mammals. Interestingly, gene conversion in the N- and C-terminal regions of avian TLR2A/TLR2B may have occurred independently. For example, in the duck we predict these events to be separated by 14 million years. These results suggest that N- and C-terminal regions of avian TLR2 genes may be under different selective constraints.
To understand more about the evolutionary origins of the TLR1 gene family, we constructed gene trees based on sequences from the TLR1 or TLR2 subfamilies free from gene conversion. For the TLR1 subfamily, two groups of orthologues were defined in birds and mammals as avian TLR1A/mammalian TLR10 and avian TLR1B/mammalian TLR1/TLR6, predicted to be the product of a gene duplication event 359 Mya. Subsequently in mammals, a further gene duplication event 270-282 Mya gave rise to the paralogues TLR1/6. This orthology between TLR1A/TLR10 and TLR1B/TLR1/TLR6 was strengthened by finding additional sequence identity in the flanking, non-coding regions of avian and mammalian genes (Additional file 12). Similarly, phylogenetic analysis strongly indicated that TLR2A to be orthologous to the functional TLR2 gene in mammals, sharing a common ancestor with TLR2B 341-356 Mya, prior to the split of reptiles and mammals 310-325 Mya. This analysis prompted us to search the genomes of mammals for TLR2 paralogues, possibly missed by current annotations. TLR2P was found in tandem with the functional TLR2 gene in human, chimpanzee, orangutan, rhesus monkey, marmoset, horse and dog. Together with the predicted time of divergence of the TLR2A/B genes (341-356 Mya), these results supposed that the TLR2B/2P genes were orthologues. In mammals, TLR2P may have recently been functional and independently pseudogenised in each species, possibly due to selection pressure against duplicate TLR2 genes. If true, this would suggest that different selective constraints are acting on the duplicate TLR2 genes in birds. This would be consistent with the more broad responses found in avian TLR1/2 heterodimers, possibly in response to a more diverse set of pathogens.
Functional and structural significance of selection constraints within the TLR1 family
Many genes of the immune system and those involved in host-pathogen interactions, such as the MHC complex, TRIM5alpha and HAVCR1, and show significant evidence of positive selection [10, 55, 56]. However, work so far suggests that most amino acid sites in TLRs are subject to purifying selection . Domain-switching experiments between TLR1/6 showed the LRR9-12 modules to play crucial roles in lipopeptide recognition . These modules constitute a critical portion of the central domain and form the ligand binding domain in TLR1/6/10. They also represent the most divergent regions and are likely to be subject to positive selection in response to pathogen infection. In mammals, our analysis detected positively selected (PS) sites in exposed regions in both TLR1 (238A, 286Q Figure 3A) and TLR6 (296D Figure 3B) at similar regions (suggesting both regions are functionally important, and possibly share a common function) in or near the central domain but not in TLR2. There were other residues inferred to be subject to positive selection, in the extracellular domain of TLR1 (174T in LRR6 Figure 3A), TLR6 (90L in LRR2 Figure 3B) and TLR2 (453G), and in the transmembrane domain of TLR1 (599S p > 0.99) and the TIR domain of TLR2 (766M). Changes in these sites may modify the ligand-binding activity of the TLR1/2 and TLR6/2 complexes. These results suggest that in mammals, changes in the ligand-binding activity of TLR1 and TLR6 may be modified through positive selection. In contrast to mammals, PS sites were detected in both avian TLR2A (304N in LRR10) and TLR2B (260V in LRR9), which map onto an exposed region (Figure 3C) in or near the ligand-binding domain and shows species-specific response to lipopeptides. The other sites are located at homologous positions in TLR2A (184T, 506S) and TLR2B (175T, 496S) in the extracellular domain at LRR6 and LRR19, respectively (Figure 3C), also suggesting those regions shared important functions. These results suggest that in birds, changes in the ligand-binding activity of TLR2A/B may be modified through positive selection, which contrasts with mammals where this appears not to be case. Also, unlike mammals, we did not find any sites showing high posterior probabilities in either the TLR1A/B molecules. These results could be due to lack of power in the small samples or they may indicate a difference in the selective constraints between TLR1 family proteins in birds and mammals.
TLR2 initiates a potent immune response by recognizing diacylated and triacylated lipopeptides. Its ligand specificity is controlled by whether it heterodimerizes with TLR1/6/10 [16–18]. Therefore extensive protein-protein interactions between TLR1/6/10 and TLR2 are to be expected, which should be reflected in a network of co-evolving amino acids. Intra-molecular CAPS analysis detected extensive co-evolution between residues within TLR1, 6 and 10 consistent with their ability to form homodimers. TLR2 does not form homodimers, which was reflected in the lack of evidence for any co-evolution between amino acid residues even with the alignment of sequences from 27 species. Inter-molecular CAPS analysis found extensive co-evolution between amino acids in TLR2 complexes with TLR1/2/10. Interestingly, amino acid residues from different domains within TLR2 were observed as co-evolving with the same residues in TLR1, TLR6 or TLR10, suggesting that there is a high codependence of evolution between the functional domains of TLR2.
TLR6 recognizes diacylated lipopeptides by the formation of a heterodimer with TLR2. Previous analysis suggested that TLR6 did not have a lipid-binding channel as found in TLR1 . This correlates with the lack of co-evolving sites in the homologous region in TLR6 but extensive co-evolution within the lipid binding pocket in TLR2 (29/90 residues). Protein-protein interactions are stronger in TLR6/TLR2 than TLR1/TLR2 where the dimerization interface in TLR6/TLR2 is 80% larger than TLR1/TLR2 . This is consistent with there being more co-evolving sites in TLR2 (35/81) than TLR6 (6/79), which has a similar number to that found in TLR1 (7/76). There are also two other small clusters of co-evolving amino acid residues, the largest also involving residues in the lipid pocket of TLR2 and residues (459, 468, 470 and 490) elsewhere in the TLR6 protein (Figure 4B), which may represent other sites of functional and/or structural significance, but this awaits further investigation.
The TLR1 gene family was previously thought to have arisen by a series of independent gene duplications in multiple vertebrates. However, we show here that the true history of gene duplication events has been masked by gene conversion between paralogous genes. Our detailed phylogenetic analysis not only clarifies the gene gains and losses within the TLR1 families of birds and mammals, but also defines orthologues between these vertebrates. TLR1 and 6 in mammals are the orthologues of avian TLR1B, and TLR10 in mammals is the orthologue of TLR1A. The functional TLR2 gene in mammals is the orthologue of the avian TLR2A genes. The mammalian orthologue of avian TLR2B has been lost in most mammals, with pseudogene fragments (TLR2P) present in a few species. In mammals, we predict amino acid sites under positive selection in TLR1, TLR2 and TLR6 but not TLR10 and detect co-evolution of amino acid residues between TLR1/TLR2 proteins that may be required to maintain their ability to form functional TLR1/2 complexes. In birds, we were able to predict positive selection in the TLR2A/B genes at functionally significant amino acid residues but not in TLR1A/B. In both birds and mammals these amino acid residues include known structural and functional sites, involved in ligand binding and dimerization. This analysis also predicts many other sites of functional and/or structural significance, but further investigations are required to define these roles. This analysis also suggests divergence in the function of avian and mammalian TLR1 gene family proteins, each under different selective constraints. In mammals, it appears that TLR1 and 6 are subject to positive selection more significant than TLR2, where in birds the reverse seems to be the case for TLR2-like genes. In contrast, the TLR10/TLR2 complex does not appear to be subject to positive selection. Our phylogenetic and structural analyses of vertebrate TLR1 family members has clarified their evolutionary history and predicts amino acid residues likely to be important in their primary function in the host's defense against invading pathogens.
Sequence clone and identities: Additional file 1, Table S1-2, Additional file 2, 5, 12 Gene conversion: Additional file 1, Table S3-4, Additional file 6, 7 Phylogenetic topology: Additional file 3, 8, 9, 10, 11 Positive selection: Additional file 13 Co-evolution: The groups, correlation coefficients and distribution on human or mouse tertiary structures of the inter-co-evolving sites in mammalian TLR1/2, TLR6/2 and TLR10/2 pairs are listed in Additional file 14, 15, 16. The co-evolutionary network identified by CAPS in mammalian TLR1/2, TLR6/2 and TLR10/2 pairs can be visualized with Cytoscape version 2.7.0 software using Additional file 14, 15, 16.
Antidocas marsupialis (springbok)
Anas platyrhynchos (duck)
Bison bison (american bison)
Bubalus bubalis (water buffalo)
Bos indicus (gudali zebu)
Bos taurus (cow)
Boselaphus tragocamelus (nilgai)
Canis familiaris (dog)
Cricetulus griseus (chinese hamster)
Capra hircus (goat)
Capra ibex (ibex)
Callithrix Jacchus (marmoset)
Dasypus novemcinctus (nine-banded armadillo)
Pongo pygmaeus (phillipsi blesbok)
Danio rerio (zebra fish)
Equus caballus (horse)
Erinaceus europaeus (hedgedog)
Giraffa camelopardalis (giraffe)
Gyps fulvus (griffon vulture)
Gallus gallus (chicken)
Homo sapiens (human)
Monodelphis domestica (gray short-tailed opossum)
Macaca fascicularis (crab-eating macaque)
Meleagris gallopavo (turkey)
Macaca mulatta (rhesus monkey)
Microcebus murinus (mouse lemur)
Mus musculus (mouse)
Ornithorhynchus anatinus (platypus)
Ovis aries (sheep)
Oryctolagus cuniculus (rabbit)
Ochotona princeps (pika)
Pongo pygmaeus (orangutan)
Pan troglodytes (chimpanzee)
Rattus norvegicus (rat)
Sus scrofa (pig)
Taeniopygia guttata (zebra finch)
Xenopus tropicalis (western clawed frog).
This work was supported by grants from the Biotechnology and Biological Sciences Research Council (BBSRC), the European Community, the National High Technology Research and Development Program of China and Natural Scientific Foundation of China.
- 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 USA. 2005, 102 (27): 9577-9582. 10.1073/pnas.0502272102.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Matsushima N, Tanaka T, Enkhbayar P, Mikami T, Taga M, Yamada K, Kuroki Y: Comparative sequence analysis of leucine-rich repeats (LRRs) within vertebrate toll-like receptors. BMC Genomics. 2007, 8: 124-10.1186/1471-2164-8-124.View ArticlePubMedPubMed CentralGoogle Scholar
- Long M, Betran E, Thornton K, Wang W: The origin of new genes: glimpses from the young and old. Nat Rev Genet. 2003, 4 (11): 865-875.View ArticlePubMedGoogle Scholar
- Roth C, Rastogi S, Arvestad L, Dittmar K, Light S, Ekman D, Liberles DA: Evolution after gene duplication: models, mechanisms, sequences, systems, and organisms. J Exp Zool B Mol Dev Evol. 2007, 308 (1): 58-73.View ArticlePubMedGoogle Scholar
- Clark AG: Invasion and maintenance of a gene duplication. Proc Natl Acad Sci USA. 1994, 91 (8): 2950-2954. 10.1073/pnas.91.8.2950.View ArticlePubMedPubMed CentralGoogle Scholar
- Lynch M, Force A: The probability of duplicate gene preservation by subfunctionalization. Genetics. 2000, 154 (1): 459-473.PubMedPubMed CentralGoogle Scholar
- Walsh B: Population-genetic models of the fates of duplicate genes. Genetica. 2003, 118 (2-3): 279-294.View ArticlePubMedGoogle Scholar
- Ortiz M, Kaessmann H, Zhang K, Bashirova A, Carrington M, Quintana-Murci L, Telenti A: The evolutionary history of the CD209 (DC-SIGN) family in humans and non-human primates. Genes Immun. 2008, 9 (6): 483-492. 10.1038/gene.2008.40.View ArticlePubMedPubMed CentralGoogle Scholar
- Kosiol C, Vinar T, da Fonseca RR, Hubisz MJ, Bustamante CD, Nielsen R, Siepel A: Patterns of positive selection in six Mammalian genomes. PLoS Genet. 2008, 4 (8): e1000144-10.1371/journal.pgen.1000144.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakajima T, Ohtani H, Satta Y, Uno Y, Akari H, Ishida T, Kimura A: Natural selection in the TLR-related genes in the course of primate evolution. Immunogenetics. 2008, 60 (12): 727-735. 10.1007/s00251-008-0332-0.View ArticlePubMedGoogle Scholar
- Barreiro LB, Ben-Ali M, Quach H, Laval G, Patin E, Pickrell JK, Bouchier C, Tichit M, Neyrolles O, Gicquel B, Kidd JR, Kidd KK, Alcais A, Ragimbeau J, Pellegrini S, Abel L, Casanova JL, Quintana-Murci L: Evolutionary dynamics of human Toll-like receptors and their different contributions to host defense. PLoS Genet. 2009, 5 (7): e1000562-10.1371/journal.pgen.1000562.View ArticlePubMedPubMed CentralGoogle Scholar
- Wlasiuk G, Nachman MW: Adaptation and constraint at Toll-like receptors in primates. Mol Biol Evol. 2010, 27 (9): 2172-2186. 10.1093/molbev/msq104.View ArticlePubMedPubMed CentralGoogle 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 (12): 6993-7000. 10.4049/jimmunol.0903092.View ArticlePubMedGoogle Scholar
- Yilmaz A, Shen S, Adelson DL, Xavier S, Zhu JJ: Identification and sequence analysis of chicken Toll-like receptors. Immunogenetics. 2005, 56 (10): 743-753. 10.1007/s00251-004-0740-8.View ArticlePubMedGoogle Scholar
- Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik SG, Lee H, Lee JO: Crystal structure of the TLR1-TLR2 heterodimer induced by binding of a tri-acylated lipopeptide. Cell. 2007, 130 (6): 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 (6): 873-884. 10.1016/j.immuni.2009.09.018.View ArticlePubMedGoogle Scholar
- Guan Y, Ranoa DR, Jiang S, Mutha SK, Li X, Baudry J, Tapping RI: Human TLRs 10 and 1 share common mechanisms of innate immune sensing but not signaling. J Immunol. 2010, 184 (9): 5094-5103. 10.4049/jimmunol.0901888.View ArticlePubMedGoogle Scholar
- Keestra AM, de Zoete MR, van Aubel RA, van Putten JP: The central leucine-rich repeat region of chicken TLR16 dictates unique ligand specificity and species-specific interaction with TLR2. J Immunol. 2007, 178 (11): 7110-7119.View ArticlePubMedGoogle Scholar
- Higuchi M, Matsuo A, Shingai M, Shida K, Ishii A, Funami K, Suzuki Y, Oshiumi H, Matsumoto M, Seya T: Combinational recognition of bacterial lipoproteins and peptidoglycan by chicken Toll-like receptor 2 subfamily. Dev Comp Immunol. 2008, 32 (2): 147-155. 10.1016/j.dci.2007.05.003.View ArticlePubMedGoogle Scholar
- Kruithof EK, Satta N, Liu JW, Dunoyer-Geindre S, Fish RJ: Gene conversion limits divergence of mammalian TLR1 and TLR6. BMC Evol Biol. 2007, 7: 148-10.1186/1471-2148-7-148.View ArticlePubMedPubMed CentralGoogle 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 (9): 967-973. 10.1016/j.dci.2009.04.001.View ArticlePubMedGoogle Scholar
- Brownlie R, Allan B: Avian toll-like receptors. Cell Tissue Res. 2011, 343 (1): 121-30. 10.1007/s00441-010-1026-0.View ArticlePubMedGoogle Scholar
- Yuan X, Zhang M, Ruan W, Song C, Ren L, Guo Y, Hu X, Li N: Construction and characterization of a duck bacterial artificial chromosome library. Anim Genet. 2006, 37 (6): 599-600. 10.1111/j.1365-2052.2006.01526.x.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.View ArticlePubMedGoogle Scholar
- Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8 (3): 195-202.View ArticlePubMedGoogle Scholar
- Gordon D, Desmarais C, Green P: Automated finishing with autofinish. Genome Res. 2001, 11 (4): 614-625. 10.1101/gr.171401.View ArticlePubMedPubMed CentralGoogle Scholar
- Ovcharenko I, Loots GG, Hardison RC, Miller W, Stubbs L: zPicture: dynamic alignment and visualization tool for analyzing conservation profiles. Genome Res. 2004, 14 (3): 472-477. 10.1101/gr.2129504.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.View ArticlePubMedPubMed CentralGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948.Google Scholar
- Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ: Jalview Version 2--a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25 (9): 1189-1191. 10.1093/bioinformatics/btp033.View ArticlePubMedPubMed CentralGoogle Scholar
- Sawyer S: Statistical tests for detecting gene conversion. Mol Biol Evol. 1989, 6 (5): 526-538.PubMedGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999, 73 (1): 152-160.PubMedPubMed CentralGoogle Scholar
- Smith VS, Page RDM, Johnson KP: Data incongruence and the problem of avian louse phylogeny. Zool Scr. 2004, 33: 239-259. 10.1111/j.0300-3256.2004.00149.x.View ArticleGoogle Scholar
- Strimmer K, von Haeseler A: Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc Natl Acad Sci USA. 1997, 94 (13): 6815-6819. 10.1073/pnas.94.13.6815.View ArticlePubMedPubMed CentralGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92 (4): 371-373. 10.1093/jhered/92.4.371.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Rannala B: Phylogenetic methods come of age: testing hypotheses in an evolutionary context. Science. 1997, 276 (5310): 227-232. 10.1126/science.276.5310.227.View ArticlePubMedGoogle Scholar
- Fares MA, Travers SA: A novel method for detecting intramolecular coevolution: adding a further dimension to selective constraints analyses. Genetics. 2006, 173 (1): 9-23. 10.1534/genetics.105.053249.View ArticlePubMedPubMed CentralGoogle Scholar
- Fares MA, McNally D: CAPS: coevolution analysis using protein sequences. Bioinformatics. 2006, 22 (22): 2821-2822. 10.1093/bioinformatics/btl493.View ArticlePubMedGoogle Scholar
- Henikoff S, Henikoff JG: Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci USA. 1992, 89 (22): 10915-10919. 10.1073/pnas.89.22.10915.View ArticlePubMedPubMed CentralGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.View ArticlePubMedPubMed CentralGoogle Scholar
- Lund O, Nielsen M, Lundegaard C, Worning P: CPHmodels 2.0: X3M a computor Program to Extract 3D Models. CASP5 conference. 2002, A102-Google Scholar
- Yang Z, Bielawski JP: Statistical methods for detecting molecular adaptation. Trends Ecol Evol. 2000, 15 (12): 496-503. 10.1016/S0169-5347(00)01994-7.View ArticlePubMedGoogle Scholar
- Janeway CA, Medzhitov R: Innate immune recognition. Annu Rev Immunol. 2002, 20: 197-216. 10.1146/annurev.immunol.20.083001.084359.View ArticlePubMedGoogle Scholar
- Beutler B, Rietschel ET: Innate immune sensing and its roots: the story of endotoxin. Nat Rev Immunol. 2003, 3 (2): 169-176. 10.1038/nri1004.View ArticlePubMedGoogle Scholar
- West AP, Koblansky AA, Ghosh S: Recognition and signaling by toll-like receptors. Annu Rev Cell Dev Biol. 2006, 22: 409-437. 10.1146/annurev.cellbio.21.122303.115827.View ArticlePubMedGoogle Scholar
- Akira S, Takeda K: Toll-like receptor signalling. Nat Rev Immunol. 2004, 4 (7): 499-511. 10.1038/nri1391.View 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 (5): 906-917. 10.1016/j.cell.2007.08.002.View ArticlePubMedGoogle Scholar
- Higgs R, Cormican P, Cahalane S, Allan B, Lloyd AT, Meade K, James T, Lynn DJ, Babiuk LA, O'Farrelly C: Induction of a novel chicken Toll-like receptor following Salmonella enterica serovar Typhimurium infection. Infect Immun. 2006, 74 (3): 1692-1698. 10.1128/IAI.74.3.1692-1698.2006.View ArticlePubMedPubMed CentralGoogle 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 (1-2): 117-127. 10.1016/j.vetimm.2004.11.003.View 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 (50): 47143-47149. 10.1074/jbc.M103902200.View ArticlePubMedGoogle Scholar
- Vallender EJ, Lahn BT: Positive selection on the human genome. Hum Mol Genet. 2004, 13 (Spec No 2): R245-254.View ArticlePubMedGoogle Scholar
- Sawyer SL, Wu LI, Emerman M, Malik HS: Positive selection of primate TRIM5alpha identifies a critical species-specific retroviral restriction domain. Proc Natl Acad Sci USA. 2005, 102 (8): 2832-2837. 10.1073/pnas.0409853102.View ArticlePubMedPubMed CentralGoogle Scholar
- Omueti KO, Beyer JM, Johnson CM, Lyle EA, Tapping RI: Domain exchange between human toll-like receptors 1 and 6 reveals a region required for lipopeptide discrimination. J Biol Chem. 2005, 280 (44): 36616-36625. 10.1074/jbc.M504320200.View ArticlePubMedGoogle Scholar
- Sanghavi SK, Shankarappa R, Reinhart TA: Genetic analysis of Toll/Interleukin-1 Receptor (TIR) domain sequences from rhesus macaque Toll-like receptors (TLRs) 1-10 reveals high homology to human TLR/TIR sequences. Immunogenetics. 2004, 56 (9): 667-674. 10.1007/s00251-004-0734-6.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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.