Diversification of the expanded teleost-specific toll-like receptor family in Atlantic cod, Gadus morhua
© Sundaram et al.; licensee BioMed Central Ltd. 2012
Received: 1 October 2012
Accepted: 26 December 2012
Published: 29 December 2012
Toll-like receptors (Tlrs) are major molecular pattern recognition receptors of the innate immune system. Atlantic cod (Gadus morhua) is the first vertebrate known to have lost most of the mammalian Tlr orthologues, particularly all bacterial recognising and other cell surface Tlrs. On the other hand, its genome encodes a unique repertoire of teleost-specific Tlrs. The aim of this study was to investigate if these duplicate Tlrs have been retained through adaptive evolution to compensate for the lack of other cell surface Tlrs in the cod genome.
In this study, one tlr21, 12 tlr22 and two tlr23 genes representing the teleost-specific Tlr family have been cloned and characterised in cod. Phylogenetic analysis grouped all tlr22 genes under a single clade, indicating that the multiple cod paralogues have arisen through lineage-specific duplications. All tlrs examined were transcribed in immune-related tissues as well as in stomach, gut and gonads of adult cod and were differentially expressed during early development. These tlrs were also differentially regulated following immune challenge by immersion with Vibrio anguillarum, indicating their role in the immune response. An increase in water temperature from 4 to 12°C was associated with a 5.5-fold down-regulation of tlr22d transcript levels in spleen. Maximum likelihood analysis with different evolution models revealed that tlr22 genes are under positive selection. A total of 24 codons were found to be positively selected, of which 19 are in the ligand binding region of ectodomain.
Positive selection pressure coupled with experimental evidence of differential expression strongly support the hypothesis that teleost-specific tlr paralogues in cod are undergoing neofunctionalisation and can recognise bacterial pathogen-associated molecular patterns to compensate for the lack of other cell surface Tlrs.
Toll-like receptors (TLRs) are an integral part of the innate immune system in all organisms and form one of the first lines of defence against invading pathogens. They are a class of pathogen recognition receptors (PRRs) that elicit specific responses against pathogens upon recognising pathogen-associated molecular patterns (PAMPs) . Most TLRs are type-I transmembrane proteins that are composed of three domains: an intracellular Toll/interleukin-1 receptor (TIR) domain, a transmembrane region and an extracellular domain. The TIR domain is highly conserved across all transmembrane TLRs and initiates signal transduction, while the variable extracellular domain is composed of leucine-rich repeats (LRR) motifs that are involved in recognising specific PAMPs . To date, 21 different TLRs have been identified across numerous vertebrates . Based on phylogenetic analyses, they are organised in six major families: TLR1 (TLRs 1, 2, 6, 10 and 14), TLR3, TLR4, TLR5, TLR7 (TLRs 7, 8, 9) and TLR11 (TLRs 11 to 13 and TLRs 21 to 23) . Avian, amphibian and teleost genomes encode for most of the mammalian orthologues, as well as additional TLRs [4–6]. Tlr15 has been identified only in birds, whereas Tlr18, Tlr19 and Tlr20 are found in teleosts. Tlr21, Tlr22 and Tlr23 are generally termed as ‘teleost-specific Tlrs’, since they are present in several teleost taxa . Nevertheless, a putative Tlr21 has been identified in chicken (Gallus gallus) , while Tlr21 and Tlr22 have been found in Xenopus tropicalis.
Even though the role of TLRs in detecting pathogens is well documented, these molecules are also known to be activated by endogenous agonists and to be involved in other biological functions. Early studies in Drosophila melanogaster have demonstrated that they control the formation of the dorso-ventral axis during embryogenesis . Heat shock proteins, inflammatory mediators and fragments of molecules from extracellular matrix, which are mainly generated in response to stress or as a consequence of tissue injury, have the potential to activate TLRs .
In spite of a large degree of conservation between teleost TLRs and their mammalian orthologues, there are some differences in signalling and their ability to recognise PAMPs . Unlike in mammals, there is not always a one to one relationship between teleost Tlr families and the PAMPs that they recognise. Immunostimulation experiments have revealed that several teleost TLRs respond to PAMPs from bacterial and viral origin . In particular, the teleost-specific Tlr22 is known to recognise dsRNA in tiger pufferfish, Takifugu rubripes, but it also responds to other PAMPs from Gram-positive and Gram-negative bacteria in other teleosts [4, 11–13].
The recently published Atlantic cod (Gadus morhua) genome draft has uncovered a unique feature of its immune system: the absence of the genes encoding for major histocompatibility complex (MHC) II, CD4 and invariant chain, which are key components of the adaptive immune system in jawed vertebrates . However, this fish has a large number of MHC I genes and a unique repertoire of TLR families in its genome. The cod genome encodes four of the mammalian homologues (tlr3, tlr7, tlr8 and tlr9), and all three teleost-specific tlrs (tlr21, tlr22 and tlr23), representing three of the six TLR families. It has lost all cell surface receptors as well as bacterial recognising mammalian homologues from the TLR1, TLR4 and TLR5 families. A single copy of tlr14 has been identified in the cod genome, but the ligand specificity of this Tlr family member is still unknown.
Gene duplication is a major force of adaptive genome evolution, since it allows duplicate genes to explore different aspects of the multidimensional functional space . Even if most duplicates degenerate into pseudogenes (nonfunctionalisation or pseudogenisation) within 50 million years following the duplication event, a remarkable number of gene duplicates are found in vertebrate genomes . One of the main mechanisms that account for the increased probability of retaining duplicate genes is the acquisition of a novel function (neofunctionalisation) by one of the copies, which is no longer required to maintain the original functions . An alternative model, which is not incompatible with subsequent neofunctionalisation, is the sharing of ancestral functions between gene duplicates (subfunctionalisation), namely partitioning of spatio-temporal expression domains . The relative contribution of neofunctionalisation and subfunctionalisation in early vertebrate evolution is still a matter of controversial debate and little is known about the role of adaptive and/ or non-adaptive pressures in the maintenance of duplicate genes (reviewed in ). One of the factors that make it difficult to distinguish these processes is the long divergence time, which clouds direct tests of selection on ancient evolutionary events. Analyses of more recent duplications, such as the ones found in some teleost lineages, may prove useful to overcome this issue .
The expanded teleost-specific Tlr family in cod is so far unique amongst teleosts and provides a good model to better understand how and why so many duplicate genes have been retained during vertebrate evolution. It is plausible that these multiple teleost-specific paralogues are retained through adaptive evolution to compensate for the lack of other cell surface Tlrs in the cod genome. To address this hypothesis, we have examined the molecular evolution and differential expression of all teleost-specific Tlrs present in the current cod genome assembly.
Sources of biological samples
Tissue and embryo samples from naïve fish
Two-year old Atlantic cod (Codfarmers ASA, Norway), reared in land based tanks at Mørkvedbukta research station (University of Nordland, Norway) were used for this study. The flow-through rearing system was supplied with sea water at 7–8°C and the fish were fed daily with a commercial diet (Amber Neptun, Skretting AS, Stavanger, Norway). Adult fish were humanely killed by immersion in an anaesthetic bath containing 0.5 g·L-1 tricaine methanesulfonate (Sigma) in accordance with the national guidelines detailed in the “Norwegian Regulation on Animal Experimentation” (Forsøksdyrutvalget, Norway). Head-kidney, kidney, spleen, liver, stomach, gut, heart, gills, muscle, skin, brain, blood and gonads were collected, snap-frozen in liquid nitrogen and stored at −80°C for subsequent RNA extraction.
Cod eggs for this study were kindly provided by Codfarmers ASA (Norway). Unfertilised eggs were immediately frozen in liquid nitrogen and stored at −80°C until RNA extraction. Eggs from individual cod spawning pairs were artificially fertilised in drum-filtered (30 μm) UV treated seawater (7°C) and maintained without aeration at a density of 10 mL·L-1. Up to one third of the seawater was replaced on a daily basis, so as to keep the oxygen concentration above 6.5 mg·L-1. Embryos at different developmental stages (4-cells, 16-cells, oblong, 25% epiboly, 75% epiboly, 10-somite, 30-somite and golden eye) and larvae (hatched, bladder stage, hindgut stage and first feeding) were observed under an optical microscope and approximately 50 specimens from each stage were collected, snap-frozen in liquid nitrogen and stored at −80°C for further analysis.
Tissue samples from fish exposed to a bacterial pathogen
This experiment was conducted at the Institute of Marine Research, Norway. One hundred and twenty adult fish with an average weight of 60 g were equally distributed in three 250 L tanks, which were part of a flow-through system that was supplied with sea water at 7–8°C. The fish were maintained in this system for a period of five weeks prior to the challenge experiment. They were fed daily with a fishmeal based feed  at 1.5% (w/w) of their body weight every day. Prior to bacterial challenge, initial control samples were collected from six fish, two per tank. Thereafter, the water flow was stopped and fish in all three tanks were subjected to bath challenge with V. anguillarum strain H610 at a concentration of 2.6·107 cfu·ml-1 for 1 h . Post-challenge samples were collected at 4 (4 hpc) and 48 (48 hpc) h after exposure. The samples collected included head-kidney, gills and spleen, which were immediately snap-frozen in liquid nitrogen and maintained at −80°C for further analysis.
The temperature stress experiment was conducted at the indoor facilities of Mørkvedbukta research station. Fifty adult cod with a mean weight of 263 ± 50 g were evenly distributed in two 500 L tanks and fed daily (Amber Neptun, Skretting AS, Norway) to 1.5% (w/w) of their body weight. Seawater at 4°C was supplied to the rearing tanks and the fish were allowed to acclimatise for a period of one week prior to the temperature stress experiment. Initial control samples were collected at the start of the experiment. Water temperature was then increased from 4°C to 12°C at a rate of 2°C·h-1 and the first post-stress samples were collected at 4 h (4 hps) when the water temperature reached 12°C. Fish were further maintained at 12°C and the final sample was collected after 72 h (72 hps). Three fish were taken from each tank at each sampling point (n=6) and humanely killed as above. Head-kidney and spleen were immediately dissected, snap-frozen in liquid nitrogen and stored at −80°C prior to RNA extraction.
RNA extraction and cDNA synthesis
The above samples were lysed in Lysing Matrix D (MP Biomedicals, USA) and total RNA extracted using QIAzol (Qiagen, Netherlands) according to the manufacturer’s instructions. Quality and quantity of total RNA were assessed by agarose electrophoresis and spectrophotometry (NanoDrop, Thermo Scientific, USA), respectively. Complementary DNA was synthesised using the Quantitect reverse transcriptase kit (Qiagen, Netherlands). Total RNA was treated with gDNA wipeout buffer provided in the reverse transcriptase kit to remove any traces of genomic DNA. Luciferase mRNA (Promega, USA) was used as an external control, as previously reported .
Cloning of Atlantic cod tlr21, tlr22 and tlr23 genes
Tlr21, Tlr22 and Tlr23 protein sequences from zebrafish (Danio rerio), stickleback (Gasterosteus aculeatus), green-spotted pufferfish (Tetraodon nigroviridis) and tiger pufferfish were used as queries in Ensembl BLAST searches (http://www.ensembl.org) against the cod genome (gadMor1 v67.1). In order to predict gene sequences, contigs and scaffolds, the above BLAST hits were further analysed using the AUGUSTUS gene prediction server at University of Greifswald . Based on predicted coding sequences, primers were designed to amplify partial coding regions of the respective paralogues (Additional file 1). Total RNA from head-kidney, kidney, spleen and gills were pooled, reverse transcribed as above and used as PCR template. Following amplification by PCR, the products of interest were analysed using gel electrophoresis, purified, cloned and sequenced as described elsewhere . The GeneRacer kit with SuperScript III RT (Invitrogen, USA) was used to perform RACE PCR in order to obtain full length cDNA sequences. Outer and inner gene specific primers for both 5′ and 3′ RACE were designed based on the partial sequences obtained above. RACE cDNA was synthesised as per the manufacturer’s protocol using total RNA pooled from head-kidney, kidney, spleen and gills. PCR products were cloned and sequenced using the primers listed in Additional file 1.
Teleost-specific tlr s of Atlantic cod
Sequence length (bp)
Protein length (aa)
A total of 41 sequences from 14 teleosts (Additional file 2) were used to perform the phylogenetic analysis to elucidate the evolution of teleost tlrs. MUSCLE (http://www.ebi.ac.uk/Tools/msa/muscle) was used to align cDNA sequences and the best nucleotide substitution model was identified using MrModelTest v2.3  and PAUP* v4.0b10 , as reported . The best model to describe the data was identified based on the Akaike information criterion (AIC). Maximum likelihood phylogenetic analysis was carried out with PhyML  and Bayesian inference was performed as detailed elsewhere . The multiple sequence alignment used for phylogenetic reconstruction and corresponding tree have been submitted to TreeBASE (http://www.treebase.org/) under the accession ID 13554.
Quantification of gene expression
Primers used for semi quantitative (RT-PCR) and real-time PCR (qPCR) of teleost-specific tlr s in Atlantic cod
qPCR primer (Forward and Reverse) (5′-3′)
RT-PCR/qPCR annealing (°C)
Qualitative RT-PCR (RT-PCR)
Gene expression across tissues and developmental stages for Atlantic cod tlr21, tlr22 and tlr23 was determined using RT-PCR. Recombinant Taq DNA polymerase (VWR, USA) was used for RT-PCR with the following thermocycling parameters: 95°C: 2 min, 35 cycles of (95°C: 15 sec, annealing temperature (Table 1): 30 sec and 72°C: 2 min) and 72°C: 7 min. Amplification was carried out in Bio-Rad C1000 thermocycler (Bio-Rad, USA). Samples were analysed by electrophoresis on 1.5% (w/v) gels and then visualised and photographed using the Kodak Gel Logic 200 Imaging System (Carestream, USA).
Real-time PCR (qPCR)
Quantification of gene expression was performed by real-time PCR with SYBR green chemistry on a LightCycler 480 (Roche, USA), as detailed elsewhere . A dissociation step with a gradient from 65°C to 97°C was performed to check the specificity of the qPCR reaction and the absence of primer dimers. Specificity was further confirmed by Sanger sequencing of qPCR products. CT values were calculated with a fluorescence threshold of 0.5 and the average of two technical replicates was used to calculate relative gene expression. Data were normalised against eef1a and ubi expression using geometric normalisation factors obtained from GeNorm (http://medgen.ugent.be/genorm/), as previously described . Relative gene expression against the initial control sample was determined and statistical analysis was performed by one-way ANOVA with Tukey’s HSD post-hoc tests using the SigmaPlot 12.0 (Systat Software Inc., USA). When the data did not meet normality or equal variance requirements, a Kruskal-Wallis one-way ANOVA by ranks and median tests was performed. Significance levels were set at P < 0.05. The sample size was too small to exclude a tank effect but there was no obvious pattern of differential gene expression in one particular tank.
Tests of selection pressure and divergence
All complete and partial tlr22 paralogues were used for selection pressure analysis, except tlr22a, tlr22e and tlr22k, since these genes had only partial sequences of 1654, 1360 and 293 bp, respectively. Coding sequences of the other nine tlr22 paralogues were aligned with MUSCLE and a codon alignment was performed using the Codon Align software (http://www.hiv.lanl.gov). The N-terminal portion of the codon aligned sequences was too variable and hence 210 bp of this region were removed prior to positive selection tests. Similarly, the C-terminal region coding for TIR domain was not included in the analysis, as it is highly conserved across all known transmembrane TLRs. Instead, a codon alignment comprising 75% (2169 bp) of the total CDS and without stop codons was used. The best nucleotide substitution model was selected using MrModelTest v2.3  and PAUP* v4.0b10  based on AIC. Differences in sequence diversity between the regions that code for different domain structures were examined by calculating the average number of synonymous (dS) and non-synonymous (dN) substitutions, insertions and deletions in the codon alignments using SNAP .
Codon based Z-tests of selection were performed to test the hypothesis of positive selection in MEGA4  using the modified Nei-Gojobori method (Jukes-Cantor) and calculating the variance with 1000 bootstrap replicates . Evolutionary distances between the nine Tlr22 paralogues were estimated by Tajima’s relative rate test . Each pair of paralogues was compared with Tlr22b as outgroup, since it was the most distant Tlr22 paralogue for which the complete sequence was available. In addition, tests for positive selection were performed using the maximum likelihood methods implemented in the CODEML program of PAML, as detailed elsewhere . The dN/dS ratio (ω) was calculated using models M0 (neutral), M1 (nearly neutral), M2 (positive selection), M7 (beta) and M8 (beta & ω). Models were compared against each other using likelihood ratio tests (LRTs). Bayesian posterior probabilities (p) were calculated for positively selected sites using naive empirical Bayes (NEB) and Bayes empirical Bayes (BEB). REL, FEL and SLAC analyses were carried out in Datamonkey (http://www.datamonkey.org) to calculate ω values for each codon, along with the corresponding probability values .
Expanded teleost-specific tlrs in Atlantic cod
Synteny and phylogenetic analysis of teleost-specific tlrs in cod
Expression profiles of teleost-specific tlrs in adult cod tissues and during early ontogeny
Tlr22c, tlr22h, tlr22j and tlr22k transcripts were found in unfertilised eggs (Figure 5B). Tlr22k was the only tlr22 paralogue to be expressed throughout early development and its transcripts were detected at epiboly, somite stage, golden eye, hatching, bladder and hindgut stages. Low expression of tlr21 and tlr22a was detected at later stages from hatching until first feeding, while tlr23a and tlr23b were not present in any of the developmental stages examined.
Differential expression following pathogen challenge
Response to temperature stress
Molecular evolution of the cod tlr22 family
Tests of selection and relative rate tests
Codon based Z-test of positive selection analysis between Atlantic cod tlr22 paralogues
Atlantic cod paralogues
Identification of positively selected sites in Atlantic cod tlr22 paralogues by maximum likelihood analysis
Positively Selected sites
ω = 1.12
M1: nearly neutral
ω0 = 0.081, ω1 = 1
p0 = 0.39, p1 = 0.61
M2: positive selection
ω0 = 0.05, p0 = 0.28
M2 vs M1
4, 6, 30, 41, 73, 126, 170, 224, 250, 274, 279, 318, 326, 333, 369, 371, 427, 443, 452, 455, 458, 478, 484, 501, 503, 505, 507, 509, 528, 531, 553, 577, 674
ω1 = 1, p1 = 0.54
2ΔlnL = 192.35,
ω 2 = 4.08, p2 = 0.18
df = 2, p = 0
ω0 = 0.18, p0 = 0.35
M3 vs M0
1, 3, 4, 6, 9, 11, 12, 13, 16, 18, 23, 25, 27, 28, 30, 33, 37, 40, 41, 43, 44, 49, 51, 52, 54, 56, 57, 59, 68, 70, 71, 73, 75, 76, 78, 80, 81, 82, 83, 94, 95, 97, 108, 110, 112, 115, 116, 122, 124, 126, 127, 129, 132, 134, 147, 152, 154, 156, 157, 159, 163, 169, 170, 172, 174, 175, 177, 178, 179, 180, 181, 189, 194, 196, 197, 210, 211, 212, 215, 224, 227, 230, 231, 233, 234, 236, 237, 240, 241, 245, 246, 247, 250, 251, 252, 253, 257, 260, 262, 267, 268, 271, 273, 274, 276, 277, 279, 281, 284, 287, 288, 294, 295, 297, 298, 301, 302, 303, 305, 307, 311, 313, 315, 316, 318, 320, 326, 330, 331, 333, 334, 335, 336, 338, 339, 340, 341, 342, 344, 345, 347, 349, 350, 352, 354, 355, 357, 358, 368, 369, 371, 375, 376, 378, 379, 382, 383, 387, 388, 392, 393, 395, 397, 400, 402, 403, 406, 410, 413, 414, 416, 417, 419, 421, 426, 427, 428, 429, 430, 432, 434, 438, 440, 441, 443, 445, 448, 452, 453, 455, 457, 458, 460, 471, 472, 474, 475, 476, 478, 480, 481, 484, 493, 495, 498, 499, 501, 503, 504, 505, 506, 507, 508, 509, 510, 517, 524, 526, 528, 529, 530, 531, 532, 533, 535, 538, 539, 548, 549, 550, 551, 552, 553, 555, 556, 557, 572, 574, 577, 578, 579, 595, 596, 598, 600, 601, 603, 613, 616, 619, 620, 622, 624, 639, 642, 643, 645, 650, 662, 666, 670, 673, 674, 680, 681, 691, 702, 712, 719
ω1 = 1.19 , p1 = 0.49
2ΔlnL = 413.43,
ω 2 = 4.36, p2 = 0.16
df = 4, p = 0
p = 0.02, q = 0.01
M8: β + ωS>1
p = 0.1, q = 0.05
M8 vs M7
1, 4, 6, 30, 41, 56, 73, 126, 170, 174, 224, 245, 250, 274, 279, 287, 295, 318, 326, 333, 342, 344, 355, 369, 371, 378, 393, 400, 427, 443, 452, 455, 457, 458, 460, 471, 474, 478, 484, 501, 503, 505, 506, 507, 509, 528, 529, 530, 531, 553, 577, 619, 674
ω = 4.06
2ΔlnL = 201.48,
p0 = 0.81, p1 = 0.19
df = 2, p = 0
4, 6, 30, 41, 49, 73, 76, 126, 147, 157, 170, 224, 246, 274, 279, 295, 318, 320, 326, 333, 342, 369, 371, 397, 400, 427, 452, 455, 478, 484, 501, 503, 507, 509, 553, 578, 613
We have characterised the full-repertoire of the highly expanded teleost-specific tlr family in Atlantic cod, which includes one tlr21, twelve tlr22 and two tlr23 genes encoded by its genome. Phylogenetic analysis of tlr paralogues from 15 teleost species recovered monophyly of all tlr22 paralogues, suggesting their origin from a common teleost ancestor. All cod tlr22 paralogues were grouped under a single clade, which indicates that they have likely arisen through tandem duplications. Cod is the first sequenced vertebrate identified to have lost all the mammalian cell surface and bacterial recognising TLR orthologues . Based on the knowledge of the functional coverage of the vertebrate TLRs, 10 TLRs are predicted to be present in the common vertebrate ancestor, namely, TLR2, 3, 4, 5, 7, 8, 9, 11, 21 and 22 . Genes encoding Tlr2, Tlr4, Tlr5 and Tlr11 are absent, while Tlr3, Tlr7-9 are intracellular Tlrs. Hence, Tlr21 and Tlr22 are the only plausible cell surface Tlrs encoded by the cod genome.
Partial synteny analysis based on the current genome build revealed conservation between Tlr22 encoding genes in cod and those in stickleback, tiger pufferfish and green-spotted pufferfish, within the genomic region containing sh3kbp1 and map3k15 genes. Sh3kbp1 (SH3-domain kinase binding protein 1) is an adapter protein involved in regulating diverse signal transduction pathways, while Map3k15 (mitogen-activated protein kinase kinase kinase 15) plays a key role in signal transduction and is essential for stress-induced apoptosis . Several tlr22 paralogues are in close proximity within the cod genome and seem to have arisen through tandem duplications. There is no uniform exon-intron structure within these tlr22 paralogues. Full length CDS of cod tlr22b, tlr22d, tlr22g and tlr22i are encoded by 5, 3, 3 and 3 exons, respectively. In the case of goldfish (Carassius auratus), zebrafish and rainbow trout (22 and 22 l) tlr22 has a single exon, while the tiger pufferfish and large yellow croaker orthologues are encoded by four, three and three exons, respectively [10, 39, 40]. Tlr22 genes in basal teleosts such as Cyprinidae and Salmonidae are represented by a single exon, while their orthologues in higher teleosts (Sciaenidae, Tetraodontidae and Gadidae) contained multiple exons. This suggests that Tlr22 may have been encoded by an uninterrupted exon in the common vertebrate ancestor and has acquired additional introns during the evolution. According to homology, synteny and phylogenetic analyses, tlr22a encoded by a single exon (based on partial sequence) seems to be the ancestral Tlr22 encoding gene and the remaining eleven paralogues have arisen through tandem duplications. It was not possible to perform a synteny analysis for tlr21 and tlr23 paralogues, since the corresponding genomic scaffolds were short and did not contain more than two genes. Cod tlr21 is represented by an uninterrupted exon in the genome, sharing this gene structure with zebrafish, tiger pufferfish and stickleback (tlr21a) homologues, while stickleback (tlr21b), channel catfish (Ictalurus punctatus) and medaka homologues are encoded by multiple exons. Tlr23 has been identified in two more teleosts, tiger pufferfish and green-spotted pufferfish, both comprising three exons each, while cod tlr23a had five exons and tlr23b is composed of three exons in its genome. Completion of current genome build will provide a better understanding of the origin of the various paralogues as well as synteny with other teleosts.
TLRs have cysteine clusters flanking either side of the LRR region with two to five cysteine residues, which are denoted LRRCT and LRRNT domains . While LRRNT regions are variable among TLRs, LRRCT contains a highly conserved consensus sequence and is known to play a crucial role in TLR signalling. The LRRCT forms a compact structure stabilised by disulphide bridges positioning the extracellular domain of the TLR relative to the membrane, as seen in the structure of human TLR3 protein . Similar to other known teleost Tlr21s, the Atlantic cod Tlr21 protein has a CxCx24Cx15C motif at its LRRCT domain, while Tlr22 and Tlr23 had a CxCx24Cx18C motif in their LRRCT, characteristic of teleost Tlr22 proteins [11, 25]. The vertebrate TLR N-terminal ectodomain is made up of several LRRs and is involved in recognising PAMPs. The ectodomain of the teleost-specific Tlrs in cod is made of up to 27 LRR repeats . Full length CDS of two tlr23 paralogues encoded for proteins containing 27 (Tlr23a) and 14 (Tlr23b) LRRs within their N-terminal domain. It is noteworthy that cod Tlr23b contains such a low number of LRRs, since vertebrate Tlrs contain 16 to 28 LRRs. As cod is the first vertebrate known to encode for two tlr23 paralogues, it is likely that tlr23a is the ancestral gene and Tlr23b has lost LRRs during evolution. Homology modelling of Tlr22b based on human TLR3 ectodomain (PDB ID: 2A0Z)  revealed a characteristic horseshoe-shaped structure. The human TLR3 ectodomain is composed of 23 LRRs forming the classical horseshoe-shaped structure and the concave inner surface is composed of 21 parallel beta sheets with the hydrophobic residues pointing inwards forming a hydrophobic core. In the cod Tlr22b model, LRR22 formed an external protrusion similar to human LRR20. LRR11 formed a very large regular alpha helix and protruded outwards similar to human LRR12. These two LRRs may be involved in the recognition of PAMPs as observed for LRR12 and LRR20 in human TLR3.
A similar pattern of cod tlr21, tlr22 and tlr23 expression was observed in zebrafish tlr21 and tlr22, channel catfish tlr21, rainbow trout tlr22 and tlr22l, large yellow croaker (Larimichthys crocea) tlr22, grass carp (Ctenopharyngodon idella) tlr22, goldfish tlr22 and orange-spotted grouper tlr21. Nevertheless, the differential expression pattern observed across tissues for tlr22 paralogues indicates that this gene may have diversified to attain specific roles in different tissues of cod. To date, cod and grass carp  are the only two teleosts known to express a tlr22 paralogue in fast muscle. Cod tlr22k was expressed in skin, similarly to channel catfish tlr22. The skin is an important mucosal defence organ  and the presence of tlr22k transcripts may trigger the innate immune response by detecting PAMPs, once they cross the mucosal layer into the skin. Cod testis expressed most tlrs, similarly to zebrafish tlr22. Several mammalian TLRs in mouse are reported to be involved in the testicular innate immune response especially in Sertoli cells . Thus, tlr22 expression in testis suggests that it may be involved in protecting the male reproductive tract in cod and other teleosts. Teleost-specific tlrs showed varied developmental expression patterns and unfertilised eggs had tlr22c, tlr22h and tlr22k transcripts, possibly derived from maternal source. Historically, Drosophila toll was identified as a key player in specification of the dorso-ventral axis during embryonic development and several toll genes were found to be expressed throughout the developmental stages . The main focus of mammalian TLR research is on the immune function of the gene and less evidence of their role in embryogenesis is established in vertebrates. A recent study on mouse brain has identified specific expression patterns of TLR7 and TLR9 expression in developing brain, which has been linked to the development of the central nervous system of vertebrates . In grass carp, tlr22 transcripts were also found during late developmental stages . This study corroborates our data, suggesting that teleost-specific tlrs may also play a role in embryogenesis.
Tiger pufferfish Tlr22 was originally thought to be a functional substitute of human TLR3, as it responds to dsRNA and may therefore promote antiviral protection in teleosts . Several in vivo and in vitro studies have shown that teleost-specific tlrs do respond to a wide variety of PAMPs originating from bacteria and parasites . An increase in expression of tlr22 was found in LPS stimulated macrophages as well as in LPS, Aeromonas salmonicida or Mycobacterium cheloni stimulated leucocytes in goldfish . LPS, peptidoglycan and poly(I:C) as well as M. marinum up-regulated expression of tlr22 in larvae and adult zebrafish, respectively [11, 48]. Continuous exposure of rainbow trout PBL, spleen and kidney to inactivated A. salmonicida, induced up to 8-fold increase in expression of tlr22 and tlr22l after 24 h. Following stimulation with poly(I:C), high levels of tlr22 transcripts in spleen of large yellow croaker . A similar effect was seen on the expression of tlr22 in grass carp infected with reovirus . In the present study, bath challenge with V. anguillarum induced a 2.1-fold increase of tlr22f, tlr22h and tlr22k transcript levels in spleen at 48 hpc compared to the initial control. In general, most of the genes analysed in this study responded to bacterial bath challenge across the three tissues that were examined, but with some tissue-specific responses. In particular, tlr22f and tlr22k were down-regulated in gills but up-regulated in spleen at 48 hpc. Our data revealed that in addition to recognising dsRNA, teleost-specific tlrs respond to PAMPs from bacterial origin.
Tlr22d was down-regulated by 5.5-fold in spleen following exposure to high temperature, indicating that it may be involved directly in the immune response to heat shock. Heat stress is known to induce an innate immune response by activating the overexpression of various heat shock proteins. In mammals, TLR2 and TLR4 up-regulation is mediated by p38-kinase and might be involved in the enhanced response of PAMP in humans monocytes induced by head shock . In the thermal shock experiment, heat shock protein 70 (hsp70) was up-regulated in skin by 3-fold at 72 hps (data not shown), thus confirming the effect of the temperature stress. Human HSP70 released into the extracellular milieu binds to TLR2 and TLR4 and exerts immunoregulatory effects through its chaperokine activity . HSP60 is also found to be dependent on TLR4 for induction of specific cytokines . Most cod teleost-specific tlr genes were differentially regulated in this temperature stress experiment, indicating that they may be involved in regulating the heat shock as well as the immune response.
In our previous study, we demonstrated that tlr22 genes from several teleost taxa (Cyprinidae, Gasterosteidae, Salmonidae, Adrianichthyidae, Tetraodontidae) are under adaptive selection pressure . Pairwise comparison of N-terminal LRR domains amongst cod tlr22 paralogues showed that they are evolving at different rates. Significant dN-dS values greater than one were observed for most comparisons, the highest being between tlr22c and tlr22i. Similar results were obtained from Tajima’s relative rate test, which in fact revealed that tlr22d is evolving considerably faster than all other paralogues. This accelerated divergence may account for its involvement in the heat shock response, since tlr22d was significantly down-regulated following thermal stress.
Average non-synonymous nucleotide substitutions within Atlantic cod tlr22 paralogues were generally much higher than synonymous ones, especially within the LRR coding region. Five PSCs at positions 73, 170, 274, 369 and 371 were found within the β sheets forming the hydrophobic core of the TLR ectodomain . At positions 73, 170 and 369, all nine cod Tlr22 paralogues predominantly contained a hydrophilic residue, while 274 and 371 were mostly hydrophobic in nature. PSCs 427, 452, 455 and 509 are present after LRR19, LRR20, LRR20 and LRR22, respectively, and are generally represented by different hydrophilic amino acids, where hydrophobic residues are normally found . LRR11 and LRR22, which protrude outwards from the horseshoe-shaped domain that recognises PAMPs, have one and four PSCs, respectively. Ligand specificity may be based on variations in the amino acids in the solvent-exposed beta sheets or on variations in the convex surface of the horseshoe-shaped domain . The substitution of hydrophilic amino acids for hydrophobic ones in and around the beta sheets will affect the polarity of the core. Also, changes i in other PSCs may be altering the polarity and the structure of the ectodomain, thus producing striking variations in the PAMP recognising sites of Tlr22 paralogues in cod. Four PSCs were found to be unique to Tlr22d, which seems to be involved in the heat shock response. Unlike other Tlr22 isoforms, the first three sites of Tlr22d had a negative charge (E318, E427 and E452) while the fourth was positive (H455).
Positive selection within duplicate genes has been related to their functional diversification through neofunctionalisation [54, 55]. Our study revealed that several PSCs in cod tlr22 genes may produce striking changes in critical protein sites and may therefore be associated with adaptation to evolving pathogens or acquisition of additional functions, such as the heat shock response. Hence, it is likely that these duplicate tlr22 genes are undergoing neofunctionalisation. Taken together with the observed asymmetric evolution rates amongst cod tlr22 paralogues, our data favour the adaptation model, as opposed to the Dykhuizen-Hartl model (reviewed in ).
We have identified and annotated 15 tlr genes representing all the members of the highly expanded teleost-specific Tlr family in Atlantic cod, which includes 12 tlr22 paralogues. They seem to have evolved through lineage-specific tandem duplications, perhaps to compensate for the absence of bacterial recognising and other cell surface Tlrs. The various tlr22 paralogues are evolving at different molecular rates and several codons in the region coding for their ligand binding domain are under adaptive selection, which may contribute to their functional diversification through neofunctionalisation. This conclusion is corroborated by experimental evidence of differential expression upon thermal shock and bacterial challenge.
This work was supported by the Faculty of Biosciences and Aquaculture, University of Nordland, Bodø, Norway and the project ‘Mucosal Immune System of Atlantic Cod - Creating a Knowledge Base (184703)’ funded by the Research Council of Norway. Arvind Sundaram’s PhD is funded by a scholarship from the University of Nordland (Norway). We are grateful to Codfarmers ASA, Norway for providing some of the biological samples used in this study. We also thank Dr. Kazue Nagasawa, Ms. Catarina Campos, Mr. Jep Lokesh and Ms. Hilde Ribe for their assistance with the fish husbandry and laboratory work at University of Nordland. The support of Dr. Øivind Bergh, Mr. Hari Rudra, Ms. Ann Cathrine Einen and Mr. Rolf Hetlelid of the Institute of Marine Research (Norway) for performing the pathogen challenge study is acknowledged.
- Akira S, Uematsu S, Takeuchi O: Pathogen recognition and innate immunity. Cell. 2006, 124: 783-801. 10.1016/j.cell.2006.02.015.PubMedView ArticleGoogle Scholar
- Imler JL, Zheng L: Biology of toll receptors: lessons from insects and mammals. J Leukoc Biol. 2004, 75: 18-26.PubMedView ArticleGoogle 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 USA. 2005, 102: 9577-9582. 10.1073/pnas.0502272102.PubMedPubMed CentralView ArticleGoogle Scholar
- Rebl A, Goldammer T, Seyfert HM: Toll-like receptor signaling in bony fish. Vet Immunol Immunopathol. 2010, 134: 139-150. 10.1016/j.vetimm.2009.09.021.PubMedView ArticleGoogle Scholar
- Ishii A, Kawasaki M, Matsumoto M, Tochinai S, Seya T: Phylogenetic and expression analysis of amphibian Xenopus toll-like receptors. Immunogenetics. 2007, 59: 281-293. 10.1007/s00251-007-0193-y.PubMedView ArticleGoogle Scholar
- Brownlie R, Allan B: Avian toll-like receptors. Cell Tissue Res. 2011, 343: 121-130. 10.1007/s00441-010-1026-0.PubMedView ArticleGoogle Scholar
- Palti Y: Toll-like receptors in bony fish: from genomics to function. Dev Comp Immunol. 2011, 35: 1263-1272. 10.1016/j.dci.2011.03.006.PubMedView ArticleGoogle Scholar
- Brownlie R, Zhu J, 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.PubMedView ArticleGoogle Scholar
- Lemaitre B, Nicolas E, Michaut L, Reichhart JM, Hoffmann JA: The dorsoventral regulatory gene cassette spatzle/Toll/cactus controls the potent antifungal response in Drosophila adults. Cell. 1996, 86: 973-983. 10.1016/S0092-8674(00)80172-5.PubMedView ArticleGoogle Scholar
- Matsuo A, Oshiumi H, Tsujita T, Mitani H, Kasai H, Yoshimizu M, Matsumoto M, Seya T: Teleost TLR22 recognizes RNA duplex to induce IFN and protect cells from birnaviruses. J Immunol. 2008, 181: 3474-3485.PubMedView ArticleGoogle Scholar
- Sundaram AYM, Consuegra S, Kiron V, Fernandes JMO: Positive selection within teleost toll-like receptors tlr21 and tlr22 subfamilies and their response to temperature stress and microbial components in zebrafish. Mol Biol Rep. 2012, 39: 8965-8975. 10.1007/s11033-012-1765-y.PubMedView ArticleGoogle Scholar
- Rebl A, Siegl E, Köllner B, Fischer U, Seyfert HM: Characterization of twin toll-like receptors from rainbow trout (Oncorhynchus mykiss): evolutionary relationship and induced expression by Aeromonas salmonicida salmonicida. Dev Comp Immunol. 2007, 31: 499-510. 10.1016/j.dci.2006.08.007.PubMedView ArticleGoogle Scholar
- Stafford JL, Ellestad KK, Magor KE, Belosevic M, Magor BG: A toll-like receptor (TLR) gene that is up-regulated in activated goldfish macrophages. Dev Comp Immunol. 2003, 27: 685-698. 10.1016/S0145-305X(03)00041-7.PubMedView ArticleGoogle Scholar
- Star B, Nederbragt AJ, Jentoft S, Grimholt U, Malmstrom M, Gregers TF, Rounge TB, Paulsen J, Solbakken MH, Sharma A, Wetten OF, Lanzen A, Winer R, Knight J, Vogel JH, Aken B, Andersen O, Lagesen K, Tooming-Klunderud A, Edvardsen RB, Tina KG, Espelund M, Nepal C, Previti C, Karlsen BO, Moum T, Skage M, Berg PR, Gjoen T, Kuhl H, et al: The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011, 477: 207-210. 10.1038/nature10342.PubMedPubMed CentralView ArticleGoogle Scholar
- Hughes AL: Gene duplication and the origin of novel proteins. Proc Natl Acad Sci USA. 2005, 102: 8791-8792. 10.1073/pnas.0503922102.PubMedPubMed CentralView ArticleGoogle Scholar
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290: 1151-1155. 10.1126/science.290.5494.1151.PubMedView ArticleGoogle Scholar
- Ohno S: Evolution by gene duplication. 1970, Berlin: Springer-VerlagView ArticleGoogle Scholar
- Force A, Lynch M, Pickett FB, Amores A, Yan Y, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999, 151: 1531-1545.PubMedPubMed CentralGoogle Scholar
- Hahn MW: Distinguishing among evolutionary models for the maintenance of gene duplicates. J Hered. 2009, 100: 605-617. 10.1093/jhered/esp047.PubMedView ArticleGoogle Scholar
- Mazet FM, Shimeld S: Gene duplication and divergence in the early evolution of vertebrates. Curr Opin Genet Dev. 2002, 12: 393-396. 10.1016/S0959-437X(02)00315-5.PubMedView ArticleGoogle Scholar
- Lokesh J, Fernandes JMO, Korsnes K, Bergh Ø, Brinchmann MF, Kiron V: Transcriptional regulation of cytokines in the intestine of Atlantic cod fed yeast derived mannan oligosaccharide or β-glucan and challenged with Vibrio anguillarum. Fish Shellfish Immunol. 2012, 33: 626-631. 10.1016/j.fsi.2012.06.017.PubMedView ArticleGoogle Scholar
- Campos C, Valente LMP, Fernandes JMO: Molecular evolution of zebrafish dnmt3 genes and thermal plasticity of their expression during embryonic development. Gene. 2012, 500: 93-100. 10.1016/j.gene.2012.03.041.PubMedView ArticleGoogle Scholar
- Stanke M, Diekhans M, Baertsch R, Haussler D: Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008, 24: 637-644. 10.1093/bioinformatics/btn013.PubMedView ArticleGoogle Scholar
- Campos C, Valente LMP, Borges P, Bizuayehu T, Fernandes JMO: Dietary lipid levels have a remarkable impact on the expression of growth-related genes in Senegalese sole (Solea senegalensis Kaup). J Exp Biol. 2010, 213: 200-209. 10.1242/jeb.033126.PubMedView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView ArticleGoogle Scholar
- Swofford DL: PAUP* version 4.0b10. Phylogenetic analysis using parsimony (*and other methods). 2003, Massachusetts: Sinauer AssociatesGoogle Scholar
- Fernandes JMO, Ruangsri J, Kiron V: Atlantic cod piscidin and its diversification through positive selection. PLoS One. 2010, 5: e9501-10.1371/journal.pone.0009501.PubMedPubMed CentralView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.PubMedView ArticleGoogle Scholar
- Nagasawa K, Giannetto A, Fernandes JMO: Photoperiod influences growth and mll (mixed-lineage leukaemia) expression in Atlantic cod. PLoS One. 2012, 7: e36908-10.1371/journal.pone.0036908.PubMedPubMed CentralView ArticleGoogle Scholar
- Korber B: HIV signature and sequence variation analysis. Computational analysis of HIV molecular sequences. Edited by: Rodrigo A, Learn G. 2000, Netherlands: Kluwer Academic Publishers, 55-72.Google Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Nei M, Kumar S: Molecular evolution and phylogenetics. 2000, New York: Oxford University PressGoogle Scholar
- Tajima F: Simple methods for testing the molecular evolutionary clock hypothesis. Genetics. 1993, 135: 599-607.PubMedPubMed CentralGoogle Scholar
- Sarropoulou E, Fernandes JMO, Mitter K, Magoulas A, Mulero V, Sepulcre MP, Figueras A, Novoa B, Kotoulas G: Evolution of a multifunctional gene: the warm temperature acclimation protein Wap65 in the European seabass Dicentrarchus labrax. Mol Phylogenet Evol. 2010, 55: 640-649. 10.1016/j.ympev.2009.10.001.PubMedView ArticleGoogle Scholar
- Delport W, Poon AFY, Frost SDW, Kosakovsky Pond SL: Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics. 2010, 26: 2455-2457. 10.1093/bioinformatics/btq429.PubMedPubMed CentralView ArticleGoogle Scholar
- Oshiumi H, Matsuo A, Matsumoto M, Seya T: Pan-vertebrate toll-like receptors during evolution. Curr Genomics. 2008, 9: 488-493. 10.2174/138920208786241234.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaji T, Yoshida S, Kawai K, Fuchigami Y, Watanabe W, Kubodera H, Kishimoto T: ASK3, a novel member of the apoptosis signal-regulating kinase family, is essential for stress-induced cell death in HeLa cells. Biochem Biophys Res Commun. 2010, 395: 213-218. 10.1016/j.bbrc.2010.03.164.PubMedView ArticleGoogle Scholar
- Xiao X, Qin Q, Chen X: Molecular characterization of a toll-like receptor 22 homologue in large yellow croaker (Pseudosciaena crocea) and promoter activity analysis of its 5′-flanking sequence. Fish Shellfish Immunol. 2011, 30: 224-233. 10.1016/j.fsi.2010.10.014.PubMedView ArticleGoogle Scholar
- Su J, Heng J, Huang T, Peng L, Yang C, Li Q: Identification, mRNA expression and genomic structure of TLR22 and its association with GCRV susceptibility/resistance in grass carp (Ctenopharyngodon idella). Dev Comp Immunol. 2012, 36: 450-462. 10.1016/j.dci.2011.08.015.PubMedView ArticleGoogle 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 USA. 2005, 102: 10976-10980. 10.1073/pnas.0505077102.PubMedPubMed CentralView ArticleGoogle Scholar
- Baoprasertkul P, Xu P, Peatman E, Kucuktas H, Liu Z: Divergent toll-like receptors in catfish (Ictalurus punctatus): TLR5S, TLR20, TLR21. Fish Shellfish Immunol. 2007, 23: 1218-1230. 10.1016/j.fsi.2007.06.002.PubMedView ArticleGoogle Scholar
- Li YW, Luo XC, Dan XM, Qiao W, Huang XZ, Li AX: Molecular cloning of orange-spotted grouper (Epinephelus coioides) TLR21 and expression analysis post Cryptocaryon irritans infection. Fish Shellfish Immunol. 2012, 32: 476-481. 10.1016/j.fsi.2011.11.021.PubMedView ArticleGoogle Scholar
- Riccioli A, Starace D, Galli R, Fuso A, Scarpa S, Palombi F, De Cesaris P, Ziparo E, Filippini A: Sertoli cells initiate testicular innate immune responses through TLR activation. J Immunol. 2006, 177: 7122-7130.PubMedView ArticleGoogle Scholar
- Kambris Z, Hoffmann JA, Imler JL, Capovilla M: Tissue and stage-specific expression of the Tolls in Drosophila embryos. Gene Expr Patterns. 2002, 2: 311-317. 10.1016/S1567-133X(02)00020-0.PubMedView ArticleGoogle Scholar
- Kaul D, Habbel P, Derkow K, Krüger C, Franzoni E, Wulczyn FG, Bereswill S, Nitsch R, Schott E, Veh R, Naumann T, Lehnardt S: Expression of toll-like receptors in the developing brain. PLoS One. 2012, 7: e37767-10.1371/journal.pone.0037767.PubMedPubMed CentralView ArticleGoogle Scholar
- Lv J, Huang R, Li H, Luo D, Liao L, Zhu Z, Wang Y: Cloning and characterization of the grass carp (Ctenopharyngodon idella) toll-like receptor 22 gene, a fish-specific gene. Fish Shellfish Immunol. 2012, 32: 1022-1031. 10.1016/j.fsi.2012.02.024.PubMedView ArticleGoogle Scholar
- Meijer AH, Gabby Krens SF, Medina Rodriguez IA, He S, Bitter W, Ewa Snaar-Jagalska B, Spaink HP: Expression analysis of the toll-like receptor and TIR domain adaptor families of zebrafish. Mol Immunol. 2004, 40: 773-783. 10.1016/j.molimm.2003.10.003.PubMedView ArticleGoogle Scholar
- Zhou J, An H, Xu H, Liu S, Cao X: Heat shock up-regulates expression of toll-like receptor-2 and toll-like receptor-4 in human monocytes via p38 kinase signal pathway. Immunology. 2005, 114: 522-530. 10.1111/j.1365-2567.2004.02112.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Asea A: Heat shock proteins and toll-like receptors. Toll-like receptors (TLRs) and innate immunity. Edited by: Bauer S, Hartmann G. 2008, Berlin/Heidelberg: Springer, 111-127. [Bauer S, Hartmann G (Series Editor): Handbook of Experimental Pharmacology]View ArticleGoogle Scholar
- Ohashi K, Burkart V, Flohé S, Kolb H: Cutting edge: heat shock protein 60 is a putative endogenous ligand of the toll-like receptor-4 complex. J Immunol. 2000, 164: 558-561.PubMedView ArticleGoogle Scholar
- Jin MS, Lee JO: Structures of the toll-like receptor family and its ligand complexes. Immunity. 2008, 29: 182-191. 10.1016/j.immuni.2008.07.007.PubMedView ArticleGoogle Scholar
- Bell JK, Mullen GED, Leifer CA, Mazzoni A, Davies DR, Segal DM: Leucine-rich repeats and pathogen recognition in toll-like receptors. Trends Immunol. 2003, 24: 528-533. 10.1016/S1471-4906(03)00242-4.PubMedView ArticleGoogle Scholar
- Gojobori J, Innan H: Potential of fish opsin gene duplications to evolve new adaptive functions. Trends Genet. 2009, 25: 198-202. 10.1016/j.tig.2009.03.008.PubMedView ArticleGoogle Scholar
- Rodríguez-Trelles F, Tarrío R, Ayala FJ: Convergent neofunctionalization by positive Darwinian selection after ancient recurrent duplications of the xanthine dehydrogenase gene. Proc Natl Acad Sci USA. 2003, 100: 13413-13417. 10.1073/pnas.1835646100.PubMedPubMed CentralView ArticleGoogle 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.