- Research article
- Open Access
Comparative genomics of Toll-like receptor signalling in five species
BMC Genomics volume 10, Article number: 216 (2009)
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.
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.
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.
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 , play a major role within the innate immune system. TLRs are found in all animals and even plant homologues have been described , 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 ). An additional molecule, the Toll-interacting protein (TOLLIP), is involved in the regulation of this process . MyD88, TIRAP, IRAK-1 and TRAF6 are also involved in TLR-induced apoptosis mediated by caspase-8 (CASP8) (reviewed by Bannerman & Goldblum ). 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 . 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 [9–11]. There is increasing evidence pointing to the strong possibility that polymorphisms in livestock TLR genes might affect immune related traits [12–14] 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 . In cattle, expression differences have been investigated in breeds of different susceptibility to Theileria annulata  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 . 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.
Localisation of TLR and signalling genes in the pig genome
The Roslin-Cambridge porcine RH panel  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 .
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  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 . 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  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 . 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 . 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 . 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  resulted in 14 alignments which all confirmed the positions predicted by the FPC map (Table 1).
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  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).
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, [35–71]). 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 .
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).
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  (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).
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  (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  (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).
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 , 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  and the BLAST analysis of the pig HTGS sequence database  we also determined the location of several of the genes using the Roslin-Cambridge porcine radiation hybrid panel . 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 . 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  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  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  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 . For example in mice 12 TLRs are known , 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 . 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 , 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  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 . FCGRT haplotypes have been associated with the capacity to transfer IgG from cow to calf in beef cattle , 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 . 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 . In addition multiple other TLR independent pathways which are activated by pathogen recognition can result in activation of IRF3 . 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  and Plasmodium chabaudi infections in mice  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 . 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  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) .
There are multiple chemokine receptors and ligands which are involved in the trafficking of leukocytes . 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 . The chemokine receptor genes have also been investigated as possible candidates for health traits in livestock . 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. ). MyD88 has been associated with a protective effect during infection with Trypanosoma  and Toxoplasma  strains. Interestingly, during malaria infections MyD88 signalling is involved in an excessive cytokine production which is responsible for most of the clinical symptoms . 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)  and mice (Listeria moncytogenes)  is consistent with the function of TLR1 and 6 and polymorphisms within these genes have been associated with tuberculosis  and malaria  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.
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  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.
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 . 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  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  against the emerging pig genome sequence .
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  using the virtual sheep genome browser . 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.  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  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  and the PiGMaP consortium linkage map .
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  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 .
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 . 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 . Groups were then ordered and orientated along the chromosomes using the order of common markers with the porcine MARC v2 map .
Ovine radiation hybrid panel
The 88 cell lines of the USUoRH 5000rad ovine radiation hybrid panel  were amplified as described above. The ovine RH maps were constructed using the rh_tsp_map 3.0 software package  and CONCORDE  linked with the QSopt package  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  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 .
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  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.
Takeda K, Kaisho T, Akira S: Toll-like receptors. Annual Review of Immunology. 2003, 21: 335-376. 10.1146/annurev.immunol.21.120601.141126.
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.
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
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.
Bannerman DD, Goldblum SE: Mechanisms of bacterial lipopolysaccharide-induced endothelial apoptosis. Am J Physiol Lung Cell Mol Physiol. 2003, 284 (6): L899-914. Review
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
Takeda K, Akira S: Toll-like receptors in innate immunity. Int Immunol. 2005, 17 (1): 1-14. 10.1093/intimm/dxh186. Review
Barton GM, Medzhitov R: Linking Toll-like receptors to IFN-alpha/beta expression. Nat Immunol. 2003, 4 (5): 432-433. 10.1038/ni0503-432.
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
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
Misch EA, Hawn TR: Toll-like receptor polymorphisms and susceptibility to human disease. Clin Sci (Lond). 2008, 114 (5): 347-360. Review
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.
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.
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.
Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/projects/geo/]
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Porcine BES Search. [http://www.sanger.ac.uk/cgi-bin/Projects/S_scrofa/BESsearch.cgi]
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.
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.
Virtual sheep genome browser. [http://www.livestockgenomics.csiro.au/perl/gbrowse.cgi/vsheep1.2/]
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.
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.
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.
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.
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.
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.
Holmberg M, Andersson-Eklund L: Quantitative trait loci affecting health traits in Swedish dairy cattle. J Dairy Sci. 2004, 87 (8): 2653-2659.
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.
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.
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.
Gerald PS, Bruns GA: Genetic determinants of viral susceptibility. Birth Defects. 1978, 14 (6A): 1-7.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
InnateDB non-redundant list (Curated & Gene expression). [http://www.innatedb.com/resources.jsp]
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.
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.
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.
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.
Leulier F, Lemaitre B: Toll-like receptors – taking an evolutionary approach. Nature Reviews Genetics. 2008, 9: 165-178. 10.1038/nrg2303.
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.
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.
Kacskovics I: Fc receptors in livestock species. Vet Imm Immunopath. 2004, 102 (4): 351-362. 10.1016/j.vetimm.2004.06.008.
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.
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.
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
Gene expression Viewer. [http://www.genomics.liv.ac.uk/tryps/GeneExpressionViewer/PublicExpressionForm.V1.html]
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.
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.
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.
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.
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.
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.
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.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.
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.
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.
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.
QSopt package. [http://www.isye.gatech.edu/~wcook/qsopt]
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.
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.
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.
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
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. , 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)
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. , 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)
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 . 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)
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)
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
- Quantitative Trait Locus
- Bacterial Artificial Chromosome
- Quantitative Trait Locus Region
- Quantitative Trait Locus Effect
- Radiation Hybrid