Skip to main content

Comparative genomics of Toll-like receptor signalling in five species

Abstract

Background

Over the last decade, several studies have identified quantitative trait loci (QTL) affecting variation of immune related traits in mammals. Recent studies in humans and mice suggest that part of this variation may be caused by polymorphisms in genes involved in Toll-like receptor (TLR) signalling. In this project, we used a comparative approach to investigate the importance of TLR-related genes in comparison with other immunologically relevant genes for resistance traits in five species by associating their genomic location with previously published immune-related QTL regions.

Results

We report the genomic localisation of TLR1-10 and ten associated signalling molecules in sheep and pig using in-silico and/or radiation hybrid (RH) mapping techniques and compare their positions with their annotated homologues in the human, cattle and mouse whole genome sequences. We also report medium-density RH maps for porcine chromosomes 8 and 13. A comparative analysis of the positions of previously published relevant QTLs allowed the identification of homologous regions that are associated with similar health traits in several species and which contain TLR related and other immunologically relevant genes. Additional evidence was gathered by examining relevant gene expression and association studies.

Conclusion

This comparative genomic approach identified eight genes as potentially causative genes for variations of health related traits. These include susceptibility to clinical mastitis in dairy cattle, general disease resistance in sheep, cattle, humans and mice, and tolerance to protozoan infection in cattle and mice. Four TLR-related genes (TLR1, 6, MyD88, IRF3) appear to be the most likely candidate genes underlying QTL regions which control the resistance to the same or similar pathogens in several species. Further studies are required to investigate the potential role of polymorphisms within these genes.

Background

The innate immune system is the first line of defence against invading pathogens and is activated by conserved pathogen associated molecular patterns (PAMPs). Toll-like receptors (TLRs), a family of signalling molecules that bind to PAMPs and consequently trigger an immune response [1], play a major role within the innate immune system. TLRs are found in all animals and even plant homologues have been described [2], illustrating the ancient origin of this gene family. Most mammalian species share ten TLR genes (TLR1-10), each detecting PAMPs with different molecular structures.

TLRs bind their ligands in a horseshoe-shaped leucine rich repeat (LRR) domain, which enables a Toll/interleukin-1 receptor (TIR) domain to associate with adapter proteins like the Toll/interleukin-1 receptor domain-containing adapter protein (TIRAP), lymphocyte antigen 96 (LY96 or MD2), or myeloid differentiation primary response protein (MyD88) which binds with the interleukin-1 receptor-associated kinase 1 (IRAK-1). This binding activates the tumour necrosis factor receptor-associated factor 6 (TRAF6), triggering a cascade which finally results in nuclear factor-kappa B (NF-κB) liberation, activating the expression of pro-inflammatory genes (reviewed by Werling & Jungi [3]). An additional molecule, the Toll-interacting protein (TOLLIP), is involved in the regulation of this process [4]. MyD88, TIRAP, IRAK-1 and TRAF6 are also involved in TLR-induced apoptosis mediated by caspase-8 (CASP8) (reviewed by Bannerman & Goldblum [5]). The toll-like receptor adaptor molecules (TICAM-1 or TRIF and TICAM-2 or TRAM) have been shown to activate TRAF6 and also to trigger interferon α or β (IFN-α/β) responses [6, 7]. The transcriptional regulation of type I interferons is coordinated, at least in part, by interferon regulatory factors 3 and 7 (IRF3/7). IRF3 and IRF7 can also be activated by kinases which are regulated by MyD88/TRAF6 [8]. The pre-eminence of TLRs and these associated signalling molecules in the initial recognition of pathogens suggests that they could be strong candidates for animal health traits.

In humans, polymorphisms within genes coding for TLR and associated signalling molecules are associated with a predisposition to several diseases [911]. There is increasing evidence pointing to the strong possibility that polymorphisms in livestock TLR genes might affect immune related traits [1214] and might explain at least part of the observed variation in disease resistance. A number of immune-related quantitative trait loci (QTL) studies have been conducted in the major livestock species and the data are made publicly available [15, 16]. However, the causal genes underlying these QTLs have not been identified. Consequently, the TLR genes and their related signalling molecules which are located within these QTLs should be considered as potential candidates for explaining phenotypic variation in disease related traits and could therefore be exploited through genetic selection for desirable alleles.

Another approach to identify genes underlying variation in immune responses is the analysis of gene expression patterns in populations with divergent resistance status pre or post infection. In mice, several differential gene expression studies involving a multitude of traits have been conducted, resulting in large datasets which are publicly available [17, 18]. However, this type of information is more limited for livestock species. Differential expression of TLR s and related genes has been analysed in the gastrointestinal tract of sheep infected with Haemonchus contortus and Trichostronglyus colubriformis [19]. In cattle, expression differences have been investigated in breeds of different susceptibility to Theileria annulata [20] and Trypanosoma congolense (Kemp, personal communication). Studies in pig have mostly addressed the role of specific TLRs during host pathogen interaction and have been reviewed recently [21]. However, to date, no studies in pigs have been undertaken to investigate TLR gene expression differences in phenotypically divergent lines.

Although genes involved in TLR signalling have been annotated in the mouse and human genomes and successfully mapped in cattle [22, 23], only a subset of TLR s and no TLR-associated signalling molecules have been localised in other livestock species. Only TLR2, TLR4, TLR6 and TLR9 have been mapped in the porcine genome [24, 25], while the locations of the sheep TLR and associated genes are currently unknown.

Here we report the genomic locations of ten TLR genes (TLR1-10) and a further ten associated signalling molecules in sheep, pig, cattle, human and mouse and compare their positions with previously published health related QTLs. We identify TLR-related genes which are located in homologous regions that are associated with similar health related traits in several species and investigate their importance by functional comparison with other linked immune related genes.

Results

Localisation of TLR and signalling genes in the pig genome

The Roslin-Cambridge porcine RH panel [26] was screened with six TLR-related genes, 20 other genes and 72 microsatellite markers, all predicted to be on porcine chromosomes 8 and 13, based on comparative analysis of the pig fingerprinted contig (FPC) map [27].

RH map of pig chromosome 8

Fifty-nine markers (23 genes and 36 microsatellites) were assigned to five linkage groups at LOD4 on porcine chromosome 8 (SSC8), which themselves were ordered into two groups corresponding to both arms (SSC8a with 24 markers and SSC8b with 35 markers) of the chromosome using markers in common with the MARC v2 porcine linkage map [28] as the scaffold (Additional file 1: Pig chromosome 8). The length of SSC8a and b was 915.2 centiray (cR) and 1312.9 cR, respectively. The resulting RH maps showed a very consistent marker order when compared to the MARC v2 map (Additional file 1: Pig chromosome 8). The here created maps are publicly available in the Arkdb database [29]. Four TLR genes were assigned to this chromosome; TLR1, 6 and 10 are closely linked between 446.9 and 490.5 cR on SSC8a, whereas TLR2 maps at 310.5 cR on SSC8b (Additional file 1: Pig chromosome 8).

RH map of pig chromosome 13

Thirty-nine markers (three genes and 36 microsatellites) were assigned to 14 linkage groups at LOD4 on porcine chromosome 13 (SSC13), which themselves were ordered along the chromosome using markers in common with the MARC v2 map [28] as the scaffold. The total map length was 2669.0 cR.

Comparison with the MARC v2 map displays a very consistent marker order (Additional file 2: Pig chromosome 13). The resulting RH map is now publicly available at the ArkDB database [29]. MyD88 was located on this chromosome at 354.4 cR and TLR9 at 588.8 cR.

Comparison with the porcine FPC map

Eighteen of the 20 TLR and associated signalling genes could be localised using comparative information between the porcine FPC map and the human whole genome sequence [30]. Two genes (TICAM1 and TOLLIP) could not be assigned to a position in the porcine FPC map because a 3 Mb human sequence fragment surrounding the localisation of the genes produced no significant alignment with any porcine clone mapped on the FPC map. Of the 18 genes with predicted locations, six (TLR1, 2, 6, 10 on SSC8 and TLR9 and MyD88 on SSC13) were mapped on the porcine RH map using the Cambridge-Roslin RH panel (Table 1). The positions for TLR4 on SSC1 and TLR9 on SSC13 on the FPC map agree with another study [25]. Thus, in-silico positions were confirmed by lab based mapping techniques for seven of the 18 TLR-related genes (Table 1). In addition, an alignment of publicly available porcine mRNA sequences of the genes against the current pre-assembled HTGS (high throughput genomic sequence) pig sequence database [31] resulted in 14 alignments which all confirmed the positions predicted by the FPC map (Table 1).

Table 1 Location of TLR and related signalling genes on the porcine FPC map compared to the porcine RH map

Localisation of TLR and signalling genes in the sheep genome

In-silico and radiation hybrid mapping in sheep

All 20 TLR-related genes could be localised in the virtual sheep genome [32, 33]. In order to confirm these in-silico positions (Table 2), primers for the 20 genes (Additional file 3: Primers used for RH mapping of TLR and signalling molecules) were used to screen the USU oRH5000 ovine radiation hybrid panel [34] to analyse linkage with previously assigned markers on the ovine RH map. Significant linkage was demonstrated for 16 of the 20 loci by LOD scores greater than 5, allowing an assignment of the loci to the ovine RH map and a comparison to the in-silico position predicted by the virtual sheep genome (Table 2).

Table 2 Location of TLR and related signalling genes on the virtual sheep genome compared to the ovine RH map

The remaining four genes were linked to markers on the ovine RH map but with LOD scores of ≤ 5.0. One of these four genes (TOLLIP) was tentatively linked to a marker on chromosome 21 (LOD = 2.49), the same location predicted by the virtual sheep genome. Therefore, while the LOD score for the RH mapping was not significant, the RH analysis supported the in-silico position (Table 2). Three genes were tentatively assigned on the RH map to locations other than predicted by the virtual sheep genome (TICAM2 and CASP8) or unlinked to any other marker (IRAK-1) but all three had non-significant LOD scores (2.97, 3.7, and 2.97, respectively), suggesting that the location predicted by the virtual sheep genome was more plausible than the RH location (Table 2).

In summary, the predicted positions of 17 genes on the virtual sheep genome were confirmed (LOD ≥ 5.0) or supported (LOD = 2.49) by RH mapping. The remaining three genes were not positioned with confidence on the RH map so the positions predicted by the virtual sheep genome could not be confirmed.

Homologous regions affecting related traits in several species

Genomic coordinates of the TLR-related genes were compared with the locations of health-related QTLs in pig, sheep, cattle, human and mouse (Table 3, [3571]). Six of the analysed genes are located in homologous QTL regions which control the susceptibility to the same or a closely related pathogen in several species. Five of them (TLR1, 6, 9, MyD88 and IRF3) are functionally involved in immune responses against the QTL associated pathogens (Table 3). In addition association studies have linked polymorphic variants of human and murine TLR1, 6 and IRF3 with susceptibility to relevant diseases. Further evidence arises also for MyD88 by differential expression in mouse strains of divergent resistance post infection with Trypanosoma congolense which is of particular interest because of the ambiguous involvement of MyD88 into the control of protozoan infections [72].

Table 3 Comparative localisation of TLR and related signalling molecules

These four genes are located in QTL regions which harbour further immunologically relevant genes. Assuming that homologous QTLs are controlled by the same genes in several species, the QTL span was narrowed down to the common block of conserved gene synteny among the species (Figures 1, 2 and 3).

Figure 1
figure 1

Position of IRF3 and overlap of QTLs in mouse, human, cattle and sheep. QTL positions are indicated by bold blue lines. Green boxes indicate the localisation of syntenic blocks conserved between species. Inversions of the gene order are indicated by red arrows. Markers located on the boundaries of the QTLs in mouse (susceptibility to Mycobaterium tuberculosis), human (Coxsackie virus resistance), cattle (susceptibility to clinical mastitis) and sheep (Nematodirus egg count) or the blocks of conserved synteny are indicated in blue. Under the assumption that the indicated QTLs are caused by the same loci, the significant region can be narrowed to two segments with a combined length of less than 7 Mb (brown line in syntenic blocks). Immunologically relevant genes located in these regions are listed in additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3.

Figure 2
figure 2

Position of MyD88 and overlap of QTLs in cattle and mouse. QTL positions are indicated by bold blue lines. Green boxes indicate the localisation of syntenic blocks conserved among species. Inversions of the gene order are indicated by red arrows. Loci located on the boundaries of the overlap between the QTL in cattle (Tryanosoma resistance) and in mice (Plasmodium chabaudi malaria) are indicated in blue. Under the assumption that the indicated QTLs are caused by the same loci, the significant region can be narrowed to segments with a combined length of approximately 10 Mb (brown line in syntenic blocks). Immunologically relevant genes located in these regions are listed in additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3.

Figure 3
figure 3

Position of the TLR1 family cluster and overlap of QTLs in cattle, mouse and human. QTL positions are indicated by bold blue lines. Green boxes indicate the localisation of the syntenic block conserved between species. Under the assumption that the indicated QTLs are caused by the same loci, the significant region can be narrowed to a segments with a length of approximately 20 Mb (brown line in syntenic blocks). Two loci (GPR125 and YIP7, in blue) limit the overlap of the QTL for susceptibility to clinical mastitis in cattle with the QTL for Listeria monocytogenes susceptibility in mice. Polymorphisms in human TLR6 (red, within the QTL overlap) have been associated with susceptibility to tuberculosis. Immunologically relevant genes located in this region are listed in additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3.

IRF3 is located in a region affecting health traits in all five species, but the QTL controls host responses for a wide range of pathogens (Table 3). The homologous QTL overlap among mouse, human, cattle and sheep comprises two blocks of conserved gene synteny between these species and has a combined extent of approximately 7 Mb (Figure 1). The region contains in human, mouse and cattle 241, 263 and 210 genes, respectively. Seventy-seven of them are listed in the innatedb non-redundant gene list of immune-related murine or humane genes [73] (Additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3). Eight genes were considered as functionally relevant according to their gene ontology (GO) annotation (Table 4).

Table 4 Potential QTL related candidate genes with functional relevance, differential expression in divergent phenotypes, and localization within QTL regions

MyD88 is located in a QTL related to protozoan infections in cattle and mice (Table 3). There are four regions of conserved synteny between bovine chromosome 22 and murine chromosome 9 which are differentially ordered and orientated between both species. Together they comprise approximately 10 Mb, within which the parasite-related QTLs in mouse and cattle overlap (Figure 2). Both QTLs share an area which in cattle and mouse comprises 97 and 100 genes, respectively. Thirty-eight genes are listed in the innatedb non-redundant gene list of immunologically relevant murine or human genes [73] (Additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3) and six were considered as functionally relevant according to their GO annotation (Table 4).

The homologous QTL regions overlaying the TLR1 family cluster controls bacterial infections in three species (Table 3). Their 20 Mb overlap region (Figure 3) contains in human, cattle and mouse 68, 63 and 71 genes, respectively. Out of those 16 genes are listed in the innatedb non-redundant gene list of immunologically relevant murine or bovine genes [73] (Additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3).

The GO annotation indicates that four of them are involved in immune responses and therefore might be functionally relevant (Table 4).

Discussion

Reliability of the pig FPC map

It was the aim of this study to use information from different sources to infer the location of 20 porcine TLR-related genes. The gene content of the pig bacterial artificial chromosome (BAC) clones predicted on the basis of BES (BAC end sequence) alignments with the human genome has been validated by subsequent sequencing of BACs in the pig genome project [31], suggesting that the FPC map is a solid tool to identify gene locations. In addition to the in-silico information deduced from the pig FPC map [27] and the BLAST analysis of the pig HTGS sequence database [31] we also determined the location of several of the genes using the Roslin-Cambridge porcine radiation hybrid panel [26]. We performed the RH analysis on porcine chromosomes 8 and 13 because these two chromosomes were expected to harbour six genes of which the position of five is of particular interest. The common location of TLR2 with the TLR1 family cluster on one chromosome is unique to pig and human. The molecules of the TLR1 family (TLR1, 6, 10) broaden their ligand spectrum by heterodimerisation with TLR2 which is then signalled via a MyD88 dependent pathway [74]. It is striking that in human and pig these closely interacting molecules are linked together, whereas in other species they are on different chromosomes. This merits a more detailed analysis of the involved genomic region.

The high consistency of the marker order between the MARC v2 [28] and RH maps (additional files 1 and 2) and the confirmation of the predicted positions with all six genes mapped using RH techniques confirms that the FPC map is a reliable source of mapping information.

Reliability of the virtual sheep genome

The virtual sheep genome [33] has been established by aligning BAC end sequence data from the CHORI-243 ovine BAC library against the sequences from the human (build hg17), bovine (build 2.0) and canine (build canFam2) genomes, and anchoring those with the ovine linkage map (version 4.6). The three inconsistencies between the virtual sheep genome and RH positions in this study (CASP8, IRAK-1, TICAM2) are likely due to limited loci on the RH map resulting in non-significant linkage, which could be resolved in the near future by adding additional markers to the RH map. Thus, the confirmation of 17 of the 20 genes by RH mapping using the USUoRH5000 panel [34] suggests that the in-silico approach for predicting gene positions using the virtual sheep genome holds great promise.

Combined approach to identify health related candidate genes

Particularly in mice a multitude of health related QTLs cover relatively large proportions of the chromosomes so that the localisation of the genes within such a region provides rather limited evidence for an involvement into the mechanisms shaping the variation of the trait. However, homologous QTL regions can be narrowed by including comparative information from other species [75, 76]. The more the size of the inter-specific QTL region can be reduced, the greater becomes the support for the candidate genes within this region.

The relationships of a gene to a phenotype can also be indicated by its functional relevance or by expression patterns which differ among phenotypes. Hence we pursue a combined approach which can provide much stronger evidence.

Genes within homologous QTL overlaps can be selected based on their ontology. Functional relevance has often been deduced from other species. This approach is not always reliable, as gene functions might change during evolution leading to limited differences among species [77]. For example in mice 12 TLRs are known [78], but only 9 of the TLR1-10 which are common in most mammals are functional in mice. However, the functional differences between mammals are relatively small and it can therefore be assumed that the gene function established in one mammalian species can in most cases be extrapolated to the others.

In addition differential expression patterns indicate that the genes might be involved in the mechanism(s) resulting in phenotypic differences, either as a consequence of a polymorphism in an upstream gene or in the gene itself. However, genes can also be associated with a divergent phenotype without a related difference in expression. Also, differential expression patterns do not necessarily indicate a direct involvement of a gene in the phenotype.

Therefore a combination of approaches is necessary and can provide much stronger evidence for or against the involvement of candidate genes in variations of disease resistance traits (Table 4).

QTL regions and candidate genes

IRF3 and linked genes

The wide range of pathogens controlled by the QTL seems to suggest that different genes might be responsible for the QTL effect in the different species. In the pig, it is likely that the QTL affecting Escherichia (E.) coli resistance is caused by a polymorphism in the FUT1 (fucosyltransferase 1) gene which is closely linked to IRF3 [36]. The FUT1 enzyme modifies a structure that enables specific binding of E. coli (ECF18) to the intestinal mucosa and therefore could not explain the QTL effects on the non-bacterial pathogens in the other species (Table 3).

For TYROBP (TYRO protein tyrosine kinase binding protein) some evidence for adaptive selection within cattle populations has been found [79], indicating that polymorphism might influence resistance traits in cattle. TYROBP activates natural killer (NK) cells and therefore plays an important role in anti-viral defence [80] and could explain the human Coxsackie virus resistance locus, but not the overlapping QTLs in the other species.

To date only polymorphisms in IRF3 and FCGRT (IgG Fc fragment receptor transporter alpha chain) have been associated with relevant health traits. The FCGRT binds IgG, and serves to transfer IgG to mucosal surfaces. In ruminants and pigs it is likely to be particularly important in colostral immunoglobulin transfer to newborns as it is expressed in the newly lactating mammary gland. However, it is also expressed in adult mammalian tissues [81]. FCGRT haplotypes have been associated with the capacity to transfer IgG from cow to calf in beef cattle [82], which is of upmost importance in newborn, but not in adult animals on which the underlying phenotypic data of the discussed QTL studies are based on. Hence there is not sufficient evidence for the involvement of FCGRT polymorphism into the variation caused by the QTLs. An IRF3 polymorphism in mice alters induction of IFN-β response and affects resistance to Listeria infections [71]. Pathogens which use the same underlying mechanism of a pathogen-driven induction of IFN-β transcription to reduce the host's defence would probably also be affected by a similar polymorphism in IRF3.

All pathogens related to the QTLs in all analysed species can potentially be recognized by TLR3 or TLR4 which can activate an immune response via a MyD88 independent pathway resulting in activation of IRF3 [83]. In addition multiple other TLR independent pathways which are activated by pathogen recognition can result in activation of IRF3 [84]. Polymorphisms altering IRF3 transcript levels could therefore affect the resistance to a range of pathogens. This indeed was observed in resistant C57BL/6 mice compared to susceptible A/J mice nine days post infection with Trypanosoma congolense (Table 4).

For the other potentially relevant genes located in the homologous QTL regions (Table 4) to our knowledge no results suggesting an involvement of these genes in the QTL effects have been reported. Assuming, that one common gene is underlying the same QTL in sheep, cattle, mouse and human, IRF3 would be a compelling candidate.

MyD88 and linked genes

The QTL effect on the susceptibility to Trypanosoma congolense infections in cattle [43] and Plasmodium chabaudi infections in mice [54] might be related, as both diseases are the result of protozoan infections which presumably carry similar PAMPs and activate the same pathways. The chromosomal overlap of these QTLs suggests that they could be caused by the same genes in both species, while a connection with an overlapping QTL for Borrelia burgdorferi in mice is less obvious [51]. To date no evidence of differential expression has been reported for any of the six potentially QTL related genes (Table 4) in response to infections in pig, cattle or sheep. However, expression studies in mice show that three genes are differentially expressed in divergent mouse phenotypes post infection with Trypanosoma congolense [85] which includes MyD88, chemokine (C-X3-C motif) receptor 1 (CX3CR1), and chemokine binding protein 2 (CCBP2) and one post infection with Yersinia enterocolitica (CX3CR1) [86].

There are multiple chemokine receptors and ligands which are involved in the trafficking of leukocytes [87]. Although several of them are coded in the region around MyD88 in the so-called chemokine receptor cluster, the comparative genomic approach (Figure 2 and Table 4) excluded several of them due to their localisation. Polymorphisms in a number of chemokine receptors are associated with susceptibility and resistance to human immunodeficiency virus (HIV) infection [87]. The chemokine receptor genes have also been investigated as possible candidates for health traits in livestock [88]. However, to date no significant associations of the chemokine receptors located in the homologous QTL with protozoan infections have been detected.

In contrast, MyD88 is due to its central position as an adaptor molecule involved in the immune responses to many different pathogens, including protozoa (reviewed by Ropert et al. [72]). MyD88 has been associated with a protective effect during infection with Trypanosoma [89] and Toxoplasma [90] strains. Interestingly, during malaria infections MyD88 signalling is involved in an excessive cytokine production which is responsible for most of the clinical symptoms [72]. Thus, a hypothetical MyD88 polymorphism affecting the gene function could balance protection against different protozoan parasites. However, to date no evidence for such polymorphism is available. It can therefore concluded, that further investigations are required to elucidate the role of MyD88, CX3CR1 or CCBP2 in the variation caused by the QTL.

TLR1 gene family cluster and linked genes

The association of this chromosomal location with the susceptibility to bacterial infections in cattle (clinical mastitis) [38] and mice (Listeria moncytogenes) [48] is consistent with the function of TLR1 and 6 and polymorphisms within these genes have been associated with tuberculosis [70] and malaria [66] in humans. The association with malaria suggests together with the differential expression in divergent mouse, sheep and cattle phenotypes post infection with protozoan or other parasites (Table 4), that the TLR1 family cluster might also be involved in the recognition of further yet unknown ligands. The ligand for TLR10 is still unknown. However, TLR10 is not functional in mice and must therefore be excluded as a common candidate for both species, although it remains a possible candidate gene for the mastitis related QTL in cattle.

Another relevant gene, RFC1 (replication factor 1), had higher transcript levels in bone marrow-derived macrophages (BMDM) isolated from disease susceptible BALB/c mice than from resistant C57BL/6 mice post infection with Yersinia enterocolitica. This indicates that different variants might play a divergent role in the disease response. The RFC1 GO annotation points among others to its involvement in the defence [GO:0006952], which includes recovery functions such as DNA repair. However, to our knowledge no gene functions linking RFC1 directly with a mechanism which could be responsible for the trait variations are known and its low expression in resistant mice might simply reflect reduced requirement for DNA repair in more resistant animals. It can therefore be concluded that the TLR1 family gene cluster is the most likely candidate for the overlaying QTLs.

Conclusion

A comparative approach enabled us to identify TLR-related genes in regions of conserved synteny among mammals that affect related traits in several species. We investigated their functional relevance for the trait in question, reviewed expression studies and analysed further immune related genes located in the regions. With the increasing availability of QTL and expression data, this approach could be extended to identify additional genes of economic interest in livestock and also to provide new insights into complex phenotypes in humans.

The genes involved in TLR signalling are suggested to be candidates for health traits in mammalian species. The most compelling evidence for involvement in pathogen susceptibility traits has been demonstrated for TLR1, TLR6, MyD88 and IRF3. Due to their close linkage and their functions or expression patterns some evidence suggests in addition FCGRT, CX3CR1, CCBP2 and TLR10 as further potential candidate genes. For FCGRT, TLR1, TLR6 and TLR10 SNPs have been established in pig [14] and cattle [82, 91, 92]. The other genes could be screened for SNPs which could then be tested for associations with health related traits in livestock.

The other TLR-related genes and further closely linked genes might be involved in mechanisms shaping immune related traits, although they were not considered here due to the limited availability of evidence. Additional investigations of polymorphisms in these genes should be pursued.

Methods

In-silico mapping using the pig FPC map

Positions of TLR-related genes in the pig genome were predicted using information of the porcine FPC map [27]. This integrated physical BAC map contains contigs constructed by fingerprinting and BAC end sequencing and is ordered using landmark maps and alignments with the human genome. The in-silico position for each locus was predicted by an alignment of the human genome sequence surrounding the localisation of the TLR-related gene [30] with the BAC end sequences in the FPC map. The reliability of this in-silico method was tested by RH mapping (see below) and BLAST analysis [93] against the emerging pig genome sequence [31].

In-silico mapping using the virtual sheep genome

Positions of the TLR-related genes were predicted in sheep by identifying the gene sequences within the virtual sheep genome [32] using the virtual sheep genome browser [33]. These in-silico positions were tested by RH mapping (see below).

Primer design for pig and sheep genes

Primers for the porcine signalling molecules and TLR s were designed from published sequences, including genomic and cDNA sequences (Additional file 3: Primers used for RH mapping of TLR and signalling molecules). Intron-exon boundaries were determined by aligning porcine cDNA sequences against either the partial pig genome sequence assembly (build SScrofa5) or against the bovine whole genome sequence assembly (build Btau 4.0), assuming conserved gene structures between both species.

Primer sequences derived from Connor et al. [23] for all TLR signalling genes but MyD88 were used for RH mapping in sheep. Oligonucleotides for the ovine MyD88 and the TLR genes were designed from ovine cDNA sequences. In order to identify intron-exon boundaries to facilitate primer design, bovine or ovine cDNA sequences were aligned with the bovine whole genome sequence assembly (build Btau 4.0).

All new primers were designed using Primer3 [94] with a targeted amplicon length of 300 bp. Other primers used for the development of RH maps for porcine chromosomes 8 and 13 were derived from the MARC v2 [28] and the PiGMaP consortium linkage map [95].

RH maps for pig and sheep

Porcine radiation hybrid panel

DNAs from 94 cell lines of the 3000 rad porcine Cambridge-Roslin Radiation Hybrid panel [26] were amplified in order to establish presence or absence of the gene in each cell line. PCR was performed with the same touchdown program for all markers: 13 cycles with an initial annealing temperatures of 67°C, dropping by 0.5°C each cycle, followed by 24 further cycles with an annealing temperature of 60°C. Genomic ovine and hamster DNA were used as positive and negative controls, respectively. The amplification of each cell line was assessed by electrophoresis in 2.8% agarose gels. All reactions were conducted twice and scored independently by eye and/or by using GelScore software [96].

Resulting vectors (Additional file 5: RH vectors of markers used for mapping in pig and sheep) were assigned to chromosomes and two- and multi-point analysis were performed using Carthagene software [97]. Fifty-nine and 39 markers were included in the RH maps for SSC8 and SSC13, containing five and 14 linkage groups (LOD4), respectively. The marker order within each group was determined using the Default algorithm of Carthagene [97]. Groups were then ordered and orientated along the chromosomes using the order of common markers with the porcine MARC v2 map [28].

Ovine radiation hybrid panel

The 88 cell lines of the USUoRH 5000rad ovine radiation hybrid panel [34] were amplified as described above. The ovine RH maps were constructed using the rh_tsp_map 3.0 software package [98] and CONCORDE [99] linked with the QSopt package [100] as described [101, 102]. Two-point RH linkage groups were constructed with a LOD of at least 5.0.

Positions of genes in human, mouse and cattle

Positions of the analysed genes in the human (NCBI 36), mouse (NCBI m37) and cattle (Btau 4.0) genomes were retrieved from the ENSEMBL website [103] by name string-search.

Definition of QTL overlaps

Markers limiting the significant QTL boundaries were identified in the relevant studies (Table 3) and their positions identified as described before. The genes limiting the QTL regions were then used to identify the homologous region in species with related QTLs. Genes located within the resulting homologous QTL overlaps were retrieved from the ENSEMBL database [103].

Selection of candidate genes based on gene ontology annotation

Genes located within the homologous QTL overlaps and listed within the InnateDB non-redundant gene list [73] were considered as functionally relevant if their GO annotations contained the keywords "immune response", "cellular defence", "response to...(any pathogen)" or "defence to ... (any pathogen)".

Analysis of expression data

Gene transcript data were retrieved from the corresponding databases and analysed for differential expression by calculating the ratio of transcript levels between populations. Differences between Means were tested by a two-tailed t-test using the corresponding Excel function. Only significantly different transcript levels (p < 0.05) were considered further.

References

  1. Takeda K, Kaisho T, Akira S: Toll-like receptors. Annual Review of Immunology. 2003, 21: 335-376. 10.1146/annurev.immunol.21.120601.141126.

    Article  CAS  PubMed  Google Scholar 

  2. Baker B, Zambryski P, Staskawicz B, Dinesh-Kumar SP: Signalling in plant-microbe interactions. Science. 1997, 276: 726-733. 10.1126/science.276.5313.726.

    Article  CAS  PubMed  Google Scholar 

  3. Werling D, Jungi TW: TOLL-like receptors linking innate and adaptive immune response. Vet Immunol Immunopathol. 2003, 91 (1): 1-12. 10.1016/S0165-2427(02)00228-3. Review

    Article  CAS  PubMed  Google Scholar 

  4. Zhang G, Ghosh S: Negative regulation of toll-like receptor-mediated signaling by Tollip. J Biol Chem. 2002, 277 (9): 7059-7065. 10.1074/jbc.M109537200.

    Article  CAS  PubMed  Google Scholar 

  5. Bannerman DD, Goldblum SE: Mechanisms of bacterial lipopolysaccharide-induced endothelial apoptosis. Am J Physiol Lung Cell Mol Physiol. 2003, 284 (6): L899-914. Review

    Article  CAS  PubMed  Google Scholar 

  6. Vogel SN, Fitzgerald KA, Fenton MJ: TLRs: differential adapter utilization by toll-like receptors mediates TLR-specific patterns of gene expression. Mol Interv. 2003, 3 (8): 466-477. 10.1124/mi.3.8.466. Review

    Article  CAS  PubMed  Google Scholar 

  7. Takeda K, Akira S: Toll-like receptors in innate immunity. Int Immunol. 2005, 17 (1): 1-14. 10.1093/intimm/dxh186. Review

    Article  CAS  PubMed  Google Scholar 

  8. Barton GM, Medzhitov R: Linking Toll-like receptors to IFN-alpha/beta expression. Nat Immunol. 2003, 4 (5): 432-433. 10.1038/ni0503-432.

    Article  CAS  PubMed  Google Scholar 

  9. Schröder NW, Schumann RR: Single nucleotide polymorphisms of Toll-like receptors and susceptibility to infectious disease. Lancet Infect Dis. 2005, 5 (3): 156-164. Review

    Article  PubMed  Google Scholar 

  10. Tapping RI, Omueti KO, Johnson CM: Genetic polymorphisms within the human Toll-like receptor 2 subfamily. Biochem Soc Trans. 2007, 35 (Pt 6): 1445-1448. 10.1042/BST0351445. Review

    Article  CAS  PubMed  Google Scholar 

  11. Misch EA, Hawn TR: Toll-like receptor polymorphisms and susceptibility to human disease. Clin Sci (Lond). 2008, 114 (5): 347-360. Review

    Article  Google Scholar 

  12. Leveque G, Forgetta V, Morroll S, Smith AL, Bumstead N, Barrow P, Loredo-Osti JC, Morgan K, Malo D: Allelic variation in TLR4 is linked to susceptibility to Salmonella enterica serovar Typhimurium infection in chickens. Infect Immun. 2003, 71 (3): 1116-1124. 10.1128/IAI.71.3.1116-1124.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Sharma BS, Leyva I, Schenkel F, Karrow NA: Association of toll-like receptor 4 polymorphisms with somatic cell score and lactation persistency in Holstein bulls. J Dairy Sci. 2006, 89 (9): 3626-3635.

    Article  CAS  PubMed  Google Scholar 

  14. Shinkai H, Tanaka M, Morozumi T, Eguchi-Ogawa T, Okumura N, Muneta Y, Awata T, Uenishi H: Biased distribution of single nucleotide polymorphisms (SNPs) in porcine Toll-like receptor 1 (TLR1), TLR2, TLR4, TLR5, and TLR6 genes. Immunogenetics. 2006, 58 (4): 324-330. 10.1007/s00251-005-0068-z.

    Article  CAS  PubMed  Google Scholar 

  15. SheepQTLdb. [http://sphinx.vet.unimelb.edu.au/cgi-bin/QTLdb/OA/browse]

  16. AnimalQTLdb. [http://www.animalgenome.org/QTLdb/]

  17. ArrayExpress. [http://www.ebi.ac.uk/microarray-as/ae/]

  18. Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/projects/geo/]

  19. Ingham A, Reverter A, Windon R, Hunt P, Menzies M: Gastrointestinal nematode challenge induces some conserved gene expression changes in the gut mucosa of genetically resistant sheep. Int J Parasitol. 2008, 38 (3–4): 431-442. 10.1016/j.ijpara.2007.07.012.

    Article  CAS  PubMed  Google Scholar 

  20. Jensen K, Paxton E, Waddington D, Talbot R, Darghouth MA, Glass EJ: Differences in the transcriptional responses induced by Theileria annulata infection in bovine monocytes derived from resistant and susceptible cattle breeds. Int J Parasitol. 2008, 38 (3–4): 313-325. 10.1016/j.ijpara.2007.08.007.

    Article  CAS  PubMed  Google Scholar 

  21. Uenishi H, Shinkai H: Porcine Toll-like receptors: The front line of pathogen monitoring and possible implications for disease resistance. Dev Comp Immunol. 2009, 33 (3): 353-361. 10.1016/j.dci.2008.06.001.

    Article  CAS  PubMed  Google Scholar 

  22. McGuire K, Jones M, Werling D, Williams JL, Glass EJ, Jann O: Radiation hybrid mapping of all 10 characterized bovine Toll-like receptors. Anim Genet. 2006, 37 (1): 47-50. 10.1111/j.1365-2052.2005.01364.x.

    Article  CAS  PubMed  Google Scholar 

  23. Connor EE, Cates EA, Williams JL, Bannerman DD: Cloning and radiation hybrid mapping of bovine toll-like receptor-4 (TLR-4) signalling molecules. Veterinary Immunology and Immunopathology. 2006, 112 (3–4): 302-308. 10.1016/j.vetimm.2006.03.003.

    Article  CAS  PubMed  Google Scholar 

  24. Muneta Y, Uenishi H, Kikuma R, Yoshihara K, Shimoji Y, Yamamoto R, Hamashima N, Yokomizo Y, Mori Y: Porcine TLR2 and TLR6: identification and their involvement in Mycoplasma hyopneumoniae infection. Journal of Interferon and Cytokine Research. 2003, 23: 583-590. 10.1089/107999003322485080.

    Article  CAS  PubMed  Google Scholar 

  25. Franceschi A, Cassini P, Scalabrini D, Botti S, Bandi CM, Giuffra E: Radiation hybrid mapping of two members of the Toll-like receptor gene family in pigs. Animal Genetics. 2004, 35: 245-264. 10.1111/j.1365-2052.2004.01134.x.

    Article  Google Scholar 

  26. Lopez-Corrales NL, Mungall C, McCarthy L, McDowall S, Goodfellow PN, Archibald AL: A radiation hybrid map of pig chromosome 4. Proceedings of the 13th European Colloquium on Cytogenetics of Domestic Animals. Állattenyésztés Takarmányozás [Hungarian Journal of Animal Production]. 1999, 48: 65-66.

    Google Scholar 

  27. Humphray SJ, Scott CE, Clark R, Marron B, Bender C, Camm N, Davis J, Jenks A, Noon A, Patel M, Sehra H, Yang F, Rogatcheva MB, Milan D, Chardon P, Rohrer G, Nonneman D, de Jong P, Meyers SN, Archibald A, Beever JE, Schook LB, Rogers J: A high utility integrated map of the pig genome. Genome Biol. 2007, 8 (7): R139-10.1186/gb-2007-8-7-r139.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Rohrer GA, Alexander LJ, Hu Z, Smith TP, Keele JW, Beattie CW: A comprehensive map of the porcine genome. Genome Res. 1996, 6 (5): 371-391. 10.1101/gr.6.5.371.

    Article  CAS  PubMed  Google Scholar 

  29. ArkDB. [http://www.thearkdb.org/]

  30. Porcine BES Search. [http://www.sanger.ac.uk/cgi-bin/Projects/S_scrofa/BESsearch.cgi]

  31. Schook LB, Beever JE, Rogers J, Humphray S, Archibald A, Chardon P, Milan D, Rohrer G, Eversole K: Swine Genome Sequencing Consortium (SGSC): a strategic roadmap for sequencing the pig genome. Comparative and Functional Genomics. 2005, 6: 251-255. 10.1002/cfg.479.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Dalrymple BP, Kirkness EF, Nefedov M, McWilliam S, Ratnakumar A, Barris W, Zhao S, Shetty J, Maddox JF, O'grady M, Nicholas F, Crawford A, Smith T, de Jong P, McEwan J, Oddy H, Cockett NE: Using comparative genomics to reorder the human genome sequence into a virtual sheep genome. Genome Bio. 2007, 8 (7): R152-10.1186/gb-2007-8-7-r152.

    Article  Google Scholar 

  33. Virtual sheep genome browser. [http://www.livestockgenomics.csiro.au/perl/gbrowse.cgi/vsheep1.2/]

  34. Wu CH, Nomura K, Goldammer T, Hadfield T, Womack JE, Cockett NE: An ovine whole-genome radiation hybrid panel used to construct an RH map of ovine chromosome 9. Anim Genet. 2007, 38 (5): 534-536. 10.1111/j.1365-2052.2007.01626.x.

    Article  CAS  PubMed  Google Scholar 

  35. Edfors-Lilja I, Wattrang E, Andersson L, Fossum C: Mapping quantitative trait loci for stress induced alterations in porcine leukocyte numbers and functions. Anim Genet. 2000, 31 (3): 186-193. 10.1046/j.1365-2052.2000.00628.x.

    Article  CAS  PubMed  Google Scholar 

  36. Meijerink E, Neuenschwander S, Fries R, Dinter A, Bertschinger HU, Stranzinger G, Vögeli P: A DNA polymorphism influencing alpha(1,2)fucosyltransferase activity of the pig FUT1 enzyme determines susceptibility of small intestinal epithelium to Escherichia coli F18 adhesion. Immunogenetics. 2000, 52 (1–2): 129-136. 10.1007/s002510000263.

    Article  CAS  PubMed  Google Scholar 

  37. Davies G, Stear MJ, Benothman M, Abuagob O, Kerr A, Mitchell S, Bishop SC: Quantitative trait loci associated with parasitic infection in Scottish blackface sheep. Heredity. 2006, 96 (3): 252-258. 10.1038/sj.hdy.6800788.

    Article  CAS  PubMed  Google Scholar 

  38. Klungland H, Sabry A, Heringstad B, Olsen HG, Gomez-Raya L, Våge DI, Olsaker I, Ødegård J, Klemetsdal G, Schulman N, Vilkki J, Ruane J, Aasland M, Rønningen K, Lien S: Quantitative trait loci affecting clinical mastitis and somatic cell count in dairy cattle. Mamm Genome. 2001, 12 (11): 837-842. 10.1007/s00335001-2081-3.

    Article  CAS  PubMed  Google Scholar 

  39. Schulman NF, Viitala SM, de Koning DJ, Virta J, Mäki-Tanila A, Vilkki JH: Quantitative trait Loci for health traits in Finnish Ayrshire cattle. J Dairy Sci. 2004, 87 (2): 443-449.

    Article  CAS  PubMed  Google Scholar 

  40. Holmberg M, Andersson-Eklund L: Quantitative trait loci affecting health traits in Swedish dairy cattle. J Dairy Sci. 2004, 87 (8): 2653-2659.

    Article  CAS  PubMed  Google Scholar 

  41. Heyen DW, Weller JI, Ron M, Band M, Beever JE, Feldmesser E, Da Y, Wiggans GR, VanRaden PM, Lewin HA: A genome scan for QTL influencing milk production and health traits in dairy cattle. Physiol Genomics. 1999, 1 (3): 165-175.

    CAS  PubMed  Google Scholar 

  42. Zhang Q, Boichard D, Hoeschele I, Ernst C, Eggen A, Murkve B, Pfister-Genskow M, Witte LA, Grignola FE, Uimari P, Thaller G, Bishop MD: Mapping quantitative trait loci for milk production and health of dairy cattle in a large outbred pedigree. Genetics. 1998, 149 (4): 1959-1973.

    PubMed Central  CAS  PubMed  Google Scholar 

  43. Hanotte O, Ronin Y, Agaba M, Nilsson P, Gelhaus A, Horstmann R, Sugimoto Y, Kemp S, Gibson J, Korol A, Soller M, Teale A: Mapping of quantitative trait loci controlling trypanotolerance in a cross of tolerant West African N'Dama and susceptible East African Boran cattle. Proc Natl Acad Sci USA. 2003, 100 (13): 7443-7448. 10.1073/pnas.1232392100.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Gerald PS, Bruns GA: Genetic determinants of viral susceptibility. Birth Defects. 1978, 14 (6A): 1-7.

    CAS  PubMed  Google Scholar 

  45. Badalová J, Svobodová M, Havelková H, Vladimirov V, Vojtísková J, Engová J, Pilcík T, Volf P, Demant P, Lipoldová M: Separation and mapping of multiple genes that control IgE level in Leishmania major infected mice. Genes Immun. 2002, 3 (4): 187-195. 10.1038/sj.gene.6363838.

    Article  PubMed  Google Scholar 

  46. Mitsos LM, Cardon LR, Fortin A, Ryan L, LaCourse R, North RJ, Gros P: Genetic control of susceptibility to infection with Mycobacterium tuberculosis in mice. Genes Immun. 2000, 1 (8): 467-477. 10.1038/sj.gene.6363712.

    Article  CAS  PubMed  Google Scholar 

  47. Sánchez F, Radaeva TV, Nikonenko BV, Persson AS, Sengul S, Schalling M, Schurr E, Apt AS, Lavebratt C: Multigenic control of disease severity after virulent Mycobacterium tuberculosis infection in mice. Infect Immun. 2003, 71 (1): 126-131. 10.1128/IAI.71.1.126-131.2003.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Boyartchuk VL, Broman KW, Mosher RE, D'Orazio SE, Starnbach MN, Dietrich WF: Multigenic control of Listeria monocytogenes susceptibility in mice. Nat Genet. 2001, 27 (3): 259-260. 10.1038/85812.

    Article  CAS  PubMed  Google Scholar 

  49. Graefe SE, Meyer BS, Müller-Myhsok B, Rüschendorf F, Drosten C, Laue T, Steeg C, Nürnberg P, Fleischer B: Murine susceptibility to Chagas' disease maps to chromosomes 5 and 17. Genes Immun. 2003, 4 (5): 321-325. 10.1038/sj.gene.6363972.

    Article  CAS  PubMed  Google Scholar 

  50. Butterfield RJ, Roper RJ, Rhein DM, Melvold RW, Haynes L, Ma RZ, Doerge RW, Teuscher C: Sex-specific quantitative trait loci govern susceptibility to Theiler's murine encephalomyelitis virus-induced demyelination. Genetics. 2003, 163 (3): 1041-1046.

    PubMed Central  CAS  PubMed  Google Scholar 

  51. Roper RJ, Weis JJ, McCracken BA, Green CB, Ma Y, Weber KS, Fairbairn D, Butterfield RJ, Potter MR, Zachary JF, Doerge RW, Teuscher C: Genetic control of susceptibility to experimental Lyme arthritis is polygenic and exhibits consistent linkage to multiple loci on chromosome 5 in four independent mouse crosses. Genes Immun. 2001, 2 (7): 388-397. 10.1038/sj.gene.6363801.

    Article  CAS  PubMed  Google Scholar 

  52. Bagot S, Campino S, Penha-Gonçalves C, Pied S, Cazenave PA, Holmberg D: Identification of two cerebral malaria resistance loci using an inbred wild-derived mouse strain. Proc Natl Acad Sci USA. 2002, 99 (15): 9919-9923. 10.1073/pnas.152215199.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Sebastiani G, Olien L, Gauthier S, Skamene E, Morgan K, Gros P, Malo D: Mapping of genetic modulators of natural resistance to infection with Salmonella typhimurium in wild-derived mice. Genomics. 1998, 47 (2): 180-186. 10.1006/geno.1997.5116.

    Article  CAS  PubMed  Google Scholar 

  54. Foote SJ, Burt RA, Baldwin TM, Presente A, Roberts AW, Laural YL, Lew AM, Marshall VM: Mouse loci for malaria-induced mortality and the control of parasitaemia. Nat Genet. 1997, 17 (4): 380-381. 10.1038/ng1297-380.

    Article  CAS  PubMed  Google Scholar 

  55. Ohno T, Ishih A, Tanaka S, Nishimura M, Terada M: Chromosomal mapping of host susceptibility loci to Angiostrongylus costaricensis nematode infection in mice. Immunogenetics. 2002, 53 (10–11): 925-929.

    CAS  PubMed  Google Scholar 

  56. Mulero-Marchese RD, Blank KJ, Sieck TG: Strain-dependent migration of lymphocytes to the vaginal mucosa after peripheral immunization. Immunogenetics. 1999, 49 (11–12): 973-980. 10.1007/s002510050581.

    Article  CAS  PubMed  Google Scholar 

  57. Bix M, Wang ZE, Thiel B, Schork NJ, Locksley RM: Genetic regulation of commitment to interleukin 4 production by a CD4(+) T cell-intrinsic mechanism. J Exp Med. 1998, 188 (12): 2289-2299. 10.1084/jem.188.12.2289.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Ogus AC, Yoldas B, Ozdemir T, Uguz A, Olcen S, Keser I, Coskun M, Cilli A, Yegin O: The Arg753GLn polymorphism of the human toll-like receptor 2 gene in tuberculosis disease. Eur Respir J. 2004, 23 (2): 219-223. 10.1183/09031936.03.00061703.

    Article  CAS  PubMed  Google Scholar 

  59. Bochud PY, Hawn TR, Aderem A: Cutting edge: a Toll-like receptor 2 polymorphism that is associated with lepromatous leprosy is unable to mediate mycobacterial signaling. J Immunol. 2003, 170 (7): 3451-3454.

    Article  CAS  PubMed  Google Scholar 

  60. Tabel Y, Berdeli A, Mir S: Association of TLR2 gene Arg753Gln polymorphism with urinary tract infection in children. Int J Immunogenet. 2007, 34 (6): 399-405. 10.1111/j.1744-313X.2007.00709.x.

    Article  CAS  PubMed  Google Scholar 

  61. Schröder NW, Diterich I, Zinke A, Eckert J, Draing C, von Baehr V, Hassler D, Priem S, Hahn K, Michelsen KS, Hartung T, Burmester GR, Göbel UB, Hermann C, Schumann RR: Heterozygous Arg753Gln polymorphism of human TLR-2 impairs immune activation by Borrelia burgdorferi and protects from late stage Lyme disease. J Immunol. 2005, 175 (4): 2534-2540.

    Article  PubMed  Google Scholar 

  62. Lorenz E, Mira JP, Cornish KL, Arbour NC, Schwartz DA: A novel polymorphism in the toll-like receptor 2 gene and its potential association with staphylococcal infection. Infect Immun. 2000, 68 (11): 6398-6401. 10.1128/IAI.68.11.6398-6401.2000.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Agnese DM, Calvano JE, Hahm SJ, Coyle SM, Corbett SA, Calvano SE, Lowry SF: Human toll-like receptor 4 mutations but not CD14 polymorphisms are associated with an increased risk of gram-negative infections. J Infect Dis. 2002, 186 (10): 1522-1525. 10.1086/344893.

    Article  CAS  PubMed  Google Scholar 

  64. Arcaroli J, Silva E, Maloney JP, He Q, Svetkauskaite D, Murphy JR, Abraham E: Variant IRAK-1 haplotype is associated with increased nuclear factor-kappaB activation and worse outcomes in sepsis. Am J Respir Crit Care Med. 2006, 173 (12): 1335-1341. 10.1164/rccm.200603-341OC.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Mockenhaupt FP, Cramer JP, Hamann L, Stegemann MS, Eckert J, Oh NR, Otchwemah RN, Dietz E, Ehrhardt S, Schröder NW, Bienzle U, Schumann RR: Toll-like receptor (TLR) polymorphisms in African children: Common TLR-4 variants predispose to severe malaria. Proc Natl Acad Sci USA. 2006, 103 (1): 177-182. 10.1073/pnas.0506803102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Leoratti FM, Farias L, Alves FP, Suarez-Mútis MC, Coura JR, Kalil J, Camargo EP, Moraes SL, Ramasawmy R: Variants in the toll-like receptor signaling pathway and clinical outcomes of malaria. J Infect Dis. 2008, 198 (5): 772-780. 10.1086/590440.

    Article  CAS  PubMed  Google Scholar 

  67. Genc MR, Vardhana S, Delaney ML, Onderdonk A, Tuomala R, Norwitz E, Witkin SS, MAP Study Group: Relationship between a toll-like receptor-4 gene polymorphism, bacterial vaginosis-related flora and vaginal cytokine responses in pregnant women. Eur J Obstet Gynecol Reprod Biol. 2004, 116 (2): 152-156. 10.1016/j.ejogrb.2004.02.010.

    Article  CAS  PubMed  Google Scholar 

  68. Hawn TR, Verbon A, Lettinga KD, Zhao LP, Li SS, Laws RJ, Skerrett SJ, Beutler B, Schroeder L, Nachman A, Ozinsky A, Smith KD, Aderem A: A common dominant TLR5 stop codon polymorphism abolishes flagellin signaling and is associated with susceptibility to legionnaires' disease. J Exp Med. 2003, 198 (10): 1563-1572. 10.1084/jem.20031220.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Khor CC, Chapman SJ, Vannberg FO, Dunne A, Murphy C, Ling EY, Frodsham AJ, Walley AJ, Kyrieleis O, Khan A, Aucan C, Segal S, Moore CE, Knox K, Campbell SJ, Lienhardt C, Scott A, Aaby P, Sow OY, Grignani RT, Sillah J, Sirugo G, Peshu N, Williams TN, Maitland K, Davies RJ, Kwiatkowski DP, Day NP, Yala D, Crook DW, Marsh K, Berkley JA, O'Neill LA, Hill AV: A Mal functional variant is associated with protection against invasive pneumococcal disease, bacteremia, malaria and tuberculosis. Nat Genet. 2007, 39 (4): 523-528. 10.1038/ng1976.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  70. Ma X, Liu Y, Gowen BB, Graviss EA, Clark AG, Musser JM: Full-exon resequencing reveals toll-like receptor variants contribute to human susceptibility to tuberculosis disease. PLoS ONE. 2007, 2 (12): e1318-10.1371/journal.pone.0001318.

    Article  PubMed Central  PubMed  Google Scholar 

  71. Garifulin O, Qi Z, Shen H, Patnala S, Green MR, Boyartchuk V: Irf3 polymorphism alters induction of interferon beta in response to Listeria monocytogenes infection. PLoS Genet. 2007, 3 (9): 1587-1597. 10.1371/journal.pgen.0030152.

    Article  CAS  PubMed  Google Scholar 

  72. Ropert C, Franklin BS, Gazzinelli RT: Role of TLRs/MyD88 in host resistance and pathogenesis during protozoan infection: lessons from malaria. Semin Immunopathol. 2008, 30 (1): 41-51. 10.1007/s00281-007-0103-2.

    Article  CAS  PubMed  Google Scholar 

  73. InnateDB non-redundant list (Curated & Gene expression). [http://www.innatedb.com/resources.jsp]

  74. Hajjar AM, O'Mahony DS, Ozinsky A, Underhill DM, Aderem A, Klebanoff SJ, Wilson CB: Cutting edge: functional interactions between toll-like receptor (TLR) 2 and TLR1 or TLR6 in response to phenol-soluble modulin. J Immunol. 2001, 166 (1): 15-19.

    Article  CAS  PubMed  Google Scholar 

  75. Barton A, Eyre S, Myerscough A, Brintnell B, Ward D, Ollier WE, Lorentzen JC, Klareskog L, Silman A, John S, Worthington J: High resolution linkage and association mapping identifies a novel rheumatoid arthritis susceptibility locus homologous to one linked to two rat models of inflammatory arthritis. Hum Mol Genet. 2001, 10 (18): 1901-1906. 10.1093/hmg/10.18.1901.

    Article  CAS  PubMed  Google Scholar 

  76. Xu C, Dai Y, Lorentzen JC, Dahlman I, Olsson T, Hillert J: Linkage analysis in multiple sclerosis of chromosomal regions syntenic to experimental autoimmune disease loci. Eur J Hum Genet. 2001, 9 (6): 458-463. 10.1038/sj.ejhg.5200653.

    Article  CAS  PubMed  Google Scholar 

  77. Werling D, Jann OC, Offord V, Glass EJ, Coffey TJ: Variation matters: TLR structure and species-specific pathogen recognition. Trends in Immunology. 2009, 30 (3): 124-130. 10.1016/j.it.2008.12.001.

    Article  CAS  PubMed  Google Scholar 

  78. Leulier F, Lemaitre B: Toll-like receptors – taking an evolutionary approach. Nature Reviews Genetics. 2008, 9: 165-178. 10.1038/nrg2303.

    Article  CAS  PubMed  Google Scholar 

  79. Freeman AR, Lynn DJ, Murray C, Bradley DG: Detecting the effects of selection at the population level in six bovine immune genes. BMC Genet. 2008, 9: 62-10.1186/1471-2156-9-62.

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  81. Kacskovics I: Fc receptors in livestock species. Vet Imm Immunopath. 2004, 102 (4): 351-362. 10.1016/j.vetimm.2004.06.008.

    Article  CAS  Google Scholar 

  82. Laegreid WW, Heaton MP, Keen JE, Grosse WM, Chitko-McKown CG, Smith TP, Keele JW, Bennett GL, Besser TE: Association of bovine neonatal Fc receptor alpha-chain gene (FCGRT) haplotypes with serum IgG concentration in newborn calves. Mamm Genome. 2002, 13 (12): 704-710. 10.1007/s00335-002-2219-y.

    Article  CAS  PubMed  Google Scholar 

  83. Kawai T, Takeuchi O, Fujita T, Inoue J, Mühlradt PF, Sato S, Hoshino K, Akira S: Lipopolysaccharide stimulates the MyD88-independent pathway and results in activation of IFN-regulatory factor 3 and the expression of a subset of lipopolysaccharide-inducible genes. J Immunol. 2001, 167 (10): 5887-5894.

    Article  CAS  PubMed  Google Scholar 

  84. Servant MJ, Grandvaux N, Hiscott J: Multiple signaling pathways leading to the activation of interferon regulatory factor 3. Biochem Pharmacol. 2002, 64 (5–6): 985-992. 10.1016/S0006-2952(02)01165-6. Review

    Article  CAS  PubMed  Google Scholar 

  85. Gene expression Viewer. [http://www.genomics.liv.ac.uk/tryps/GeneExpressionViewer/PublicExpressionForm.V1.html]

  86. Van Erp K, Dach K, Koch I, Heesemann J: Role of strain differences on host resistance and the transcriptional response of macrophages to infection with Yersinia enterocolitica. Physiol Genomics. 2006, 25 (1): 75-84. 10.1152/physiolgenomics.00188.2005.

    Article  CAS  PubMed  Google Scholar 

  87. Allen SJ, Crown SE, Handel TM: Chemokine: Receptor structure, interactions, and antagonism. Ann Rev Imm. 2007, 25: 787-820. 10.1146/annurev.immunol.24.021605.090529.

    Article  CAS  Google Scholar 

  88. Leyva-Baca I, Schenkel F, Sharma BS, Jansen GB, Karrow NA: Identification of single nucleotide polymorphisms in the bovine CCL2, IL8, CCR2 and IL8RA genes and their association with health and production in Canadian Holsteins. Anim Genet. 2007, 38 (3): 198-202. 10.1111/j.1365-2052.2007.01588.x.

    Article  CAS  PubMed  Google Scholar 

  89. Drennan MB, Stijlemans B, Abbeele Van den J, Quesniaux VJ, Barkhuizen M, Brombacher F, De Baetselier P, Ryffel B, Magez S: The induction of a type 1 immune response following a Trypanosoma brucei infection is MyD88 dependent. J Immunol. 2005, 175 (4): 2501-2509.

    Article  CAS  PubMed  Google Scholar 

  90. Sukhumavasi W, Egan CE, Warren AL, Taylor GA, Fox BA, Bzik DJ, Denkers EY: TLR adaptor MyD88 is essential for pathogen control during oral toxoplasma gondii infection but not adaptive immunity induced by a vaccine strain of the parasite. J Immunol. 2008, 181 (5): 3464-3473.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  91. Seabury CM, Womack JE: Analysis of sequence variability and protein domain architectures for bovine peptidoglycan recognition protein 1 and Toll-like receptors 2 and 6. Genomics. 2008, 92 (4): 235-245. 10.1016/j.ygeno.2008.06.005.

    Article  CAS  PubMed  Google Scholar 

  92. Seabury CM, Cargill EJ, Womack JE: Sequence variability and protein domain architectures for bovine Toll-like receptors 1, 5, and 10. Genomics. 2007, 90 (4): 502-515. 10.1016/j.ygeno.2007.07.001.

    Article  CAS  PubMed  Google Scholar 

  93. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.

    Article  CAS  PubMed  Google Scholar 

  94. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  95. Archibald AL, Haley CS, Brown JF, Couperwhite S, McQueen HA, Nicholson D, Coppieters W, Weghe Van de A, Stratil A, Winter AK, Fredholm M, Larsen NJ, Nielsen VH, Milan D, Woloszyn N, Robic A, Dalens M, Riquet J, Gellin J, Caritez JC, Burgaud G, Ollivier L, Bidanel JP, Vaiman M, Renard C, Geldermann H, Davoli R, Ruyter D, Verstege EJM, Groenen MAM, Davies W, Høyheim B, Keiserud A, Andersson L, Ellegren H, Johansson M, Marklund L, Miller , Anderson Dear DV, Signer E, Jeffreys AJ, Moran C, Le Tissier P, Muladno , Rothschild MF, Tuggle CK, Vaske D, Helm J, Liu HC, Rahman A, Yu TP, Larson RG, Schmitz CB: The PiGMaP consortium linkage map of the pig (Sus scrofa). Mammalian Genome. 1995, 6: 157-175. 10.1007/BF00293008.

    Article  CAS  PubMed  Google Scholar 

  96. GelScore. [http://www.wesbarris.com/GelScore/]

  97. Schiex T, Gaspin C: CARTHAGENE: constructing and joining maximum likelihood genetic maps. Proceedings of International Conference on Intellligent Systems for Molecular Biology. 1997, 5: 258-267.

    CAS  Google Scholar 

  98. Schäffer AA, Rice ES, Cook W, Agarwala R: rh_tsp_map 3.0: end-to-end radiation hybrid mapping with improved speed and quality control. Bioinformatics. 2007, 23 (9): 1156-1158. 10.1093/bioinformatics/btm077.

    Article  PubMed Central  PubMed  Google Scholar 

  99. CONCORDE. [http://www.isye.gatech.edu/~wcook/rh]

  100. QSopt package. [http://www.isye.gatech.edu/~wcook/qsopt]

  101. Agarwala R, Applegate DL, Maglott D, Schuler GD, Schäffer AA: A fast and scalable radiation hybrid map construction and integration strategy. Genome Res. 2000, 10 (3): 350-364. 10.1101/gr.10.3.350.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  102. Schäffer AA, Rice ES, Cook W, Agarwala R: rh_tsp_map 3.0: end-to-end radiation hybrid mapping with improved speed and quality control. Bioinformatics. 2007, 23 (9): 1156-1158. 10.1093/bioinformatics/btm077.

    Article  PubMed Central  PubMed  Google Scholar 

  103. Ensembl. [http://www.ensembl.org]

Download references

Acknowledgements

The study was financed by the Biotechnology and Biological Sciences Research Council and Pfizer Inc [grant number BBD5240401], the Biotechnology and Biological Sciences Research Council [BBSRC Institute Strategic Programme Grant, grant numbers EGA16307, PAG04437, BBE0105201] and the Wellcome Trust Host-Pathogen Project [grant number GR066764MA]. The authors wish to thank two referees for their helpful comments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Oliver C Jann.

Additional information

Authors' contributions

OJ designed the primers, screened the ovine and porcine RH panels, calculated the pig RH maps, built in-silico maps, performed the comparative QTL overlap study and prepared the draft manuscript. AK screened the porcine RH panel with microsatellite markers on SSC8. NLC screened the porcine RH panel with microsatellite markers on SSC13. SIA searched literature and murine microarray databases for evidence of TLR-related transcriptional response variation. KJ searched literature and bovine and ovine microarray databases for evidence of TLR-related transcriptional response variation. TAA searched literature and porcine microarray databases for evidence of TLR-related transcriptional response variation. HF built porcine in-silico maps. CW mapped RH vectors to the sheep RH map. NEC supervised ovine RH mapping and helped draft the manuscript. ALA supervised the pig RH project, conducted searches of the pig genome sequence and reviewed the manuscript, EJG designed and supervised the study and helped with drafting of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2008_2100_MOESM1_ESM.ppt

Additional file 1: Pig chromosome 8. The file contains an RH map of porcine chromosome 8a and 8b (left) compared to the MARC v2 linkage map (Rohrer et al. [28], right). Common markers are connected by red lines. RH linkage groups (LOD4) are indicated by blue lines and the outer most markers of each group are indicated. TLR-related genes are boxed. Distances on the RH maps are indicated in cR and on the linkage map in cM. (PPT 60 KB)

12864_2008_2100_MOESM2_ESM.ppt

Additional file 2: Pig chromosome 13. The file contains an RH map of porcine chromosome 13 (left) compared to the MARC v2 linkage map (Rohrer et al. [28], right). Markers common to both maps are connected by red lines. RH linkage groups (LOD4) are indicated by blue lines and the extreme markers of each group are indicated. TLR-related genes are boxed. Distances on the RH map are indicated in cR and on the linkage map in cM. (PPT 62 KB)

12864_2008_2100_MOESM3_ESM.xls

Additional file 3: Primers used for RH mapping of TLR and signalling molecules. The file contains the primer sequences for the mapped loci in sheep and pig. (XLS 35 KB)

12864_2008_2100_MOESM4_ESM.xls

Additional file 4: Immunologically relevant genes in regions of conserved synteny surrounding the TLR1 family cluster, MyD88 and IRF3. The file contains a list of genes located in the regions of conserved synteny which overlap with the discussed QTLs and which are listed in the innatedb gene list [73]. Ensembl IDs, gene names, murine orthologs, gene ontologies (GO term) and chromosomal localisation in human are given. Genes unique to mouse or murine genes for which the human orthologs are not listed in the innatedb gene list are itemized with their position in mouse. Genes considered as functionally relevant are highlighted by green background. (XLS 152 KB)

12864_2008_2100_MOESM5_ESM.xls

Additional file 5: RH vectors of markers used for mapping in pig and sheep. The file contains RH vectors for each locus mapped in sheep and pig. Each position in the vector represents a cell line with "0" indicating no retention, "1" indicating retention and "2" indicating ambiguous results. (XLS 96 KB)

Authors’ original submitted files for images

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

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Jann, O.C., King, A., Corrales, N.L. et al. Comparative genomics of Toll-like receptor signalling in five species. BMC Genomics 10, 216 (2009). https://doi.org/10.1186/1471-2164-10-216

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-216

Keywords