- Research article
- Open Access
Insights into the European rabbit (Oryctolagus cuniculus) innate immune system: genetic diversity of the toll-like receptor 3 (TLR3) in wild populations and domestic breeds
BMC Genetics volume 14, Article number: 73 (2013)
Toll-like receptors (TLRs) belong to the innate immune system and are a major class of pattern recognition receptors representing the first line of the innate immune response. The TLR molecule is structurally composed by an ectodomain that contains leucine rich repeats (LRRs) that interact with pathogen associated molecular patterns (PAMPs), a transmembrane domain and a conserved cytoplasmic domain designated TIR (Toll-IL1 receptor) that is responsible for the intracellular signaling. TLR3 has been associated with the direct recognition of double-stranded viral RNA resulting from viral replication, while TLR7 and TLR8 target single-stranded viral RNA. In the European rabbit (Oryctolagus cuniculus), TLR7 and TLR8 were reported to be absent and pseudogenised, respectively, making TLR3 the only available TLR for the recognition of viral RNA. Thus, the levels of diversity of TLR3 were evaluated in the European rabbit by analysing different genetic backgrounds and exposure to pathogen pressures.
We detected 41 single nucleotide polymorphisms (SNPs) in the coding sequence of TLR3. The highest diversity was observed in the wild populations of Iberian Peninsula, between 22–33 polymorphic positions. In the French population, 18 SNPs were observed and only 4 polymorphic positions were detected in the domestic breeds. 14 non-synonymous substitutions were observed, most of them in the LRR molecules. The remaining were scattered across the transmembrane and TIR domains.
The study of TLR3 in European rabbit populations might be relevant to understand the interplay between RNA viruses and innate immunity. Wild rabbit populations presented more diversity than domestic breeds and other mammals previously studied. This might be linked to the absence of population bottlenecks during their evolution and to the almost inexistence of man-mediated selection. The observed variability might have also been potentiated by the contact of the wild populations with various pathogens. The study of these patterns of variability might reveal scenarios of host-pathogen interaction and identify TLR3 polymorphisms’ that arose due to viral pathogens affecting wild populations.
The Toll like receptor (TLR) molecules are pattern recognition receptors (PRR) that mediate the recognition of pathogens by the innate immune system and can also lead to the activation of the adaptive immune response. To date, thirteen mammalian TLRs have been described that recognise pathogen associated molecular patterns (PAMPs), including bacterial lipoproteins, lipopolysaccharide (LPS) and viral nucleic acids . TLR3 is triggered directly by double-stranded RNA (dsRNA) viruses or dsRNA resulting from their replication process . The human TLR3 molecule is structurally characterised by the presence of an ectodomain that encompasses 23 leucine-rich repeats (LRRs) capped at each end by a LRR-NT and a LRR-CT molecule [3, 4]. The LRRs main function is the interaction with PAMPs which in TLR3 occurs in the concave face of the TLR3-TLR3 homodimer . In addition, human TLR3 has a transmembrane domain and a conserved cytoplasmic domain homologous to interleukin-1 receptor (IL1R) and IL-18 receptor (IL18R) designated TIR domain that is responsible for the intracellular signalling [5–7].
The lack of this receptor has been associated with higher susceptibility and mortality in infections by cytomegalovirus (CMV) and herpes simplex virus (HSV) in mice and humans [8, 9]. On the other hand, some viruses, like West Nile Virus (WNV), benefit from the interaction with TLR3 that facilitates infection and dissemination in the organism . In line with this disease association, several studies have pointed out a link between polymorphisms in TLR3 and resistance/susceptibility to diseases [11, 12]. For instance, in humans, the TLR3 single nucleotide polymorphism (SNPs) involved in the amino acid replacement N284I is responsible for the partial or total loss of activity of the receptor possibly by destabilising LRR architecture and the P554S SNP is associated with encephalitis caused by HSV .
In the present study a general approach was made to assess TLR3 diversity in wild and domestic European rabbits (Oryctolagus cuniculus) that have different genetic backgrounds and are exposed differently to pathogens. Viral TLRs (TLR3, 7, 8 and 9) act as the first line of defence against viral pathogens. However, in rabbit, TLR7 and TLR8 which are responsible for the recognition of ssRNA viruses are absent or pseudogenised, respectively , making TLR3 the only available TLR for the recognition of viral RNA.
The European rabbit originated in the Iberian Peninsula (IP) where two subspecies that diverged about 2 million years ago co-exist: Oryctolagus cuniculus algirus (OCA) in the Southwest part and Oryctolagus cuniculus cuniculus (OCC) in the Northeast part of IP  (Figure 1). More recently, the Northeast populations expanded to France and from these the species has been domesticated in the last 1,500 years with more than 200 breeds currently being recognised . Thus, rabbit populations found outside the IP represent a subset of the original gene pool. Here, we have characterised genetically the TLR3 coding sequence in wild rabbits from the Iberian Peninsula and from France, and from domestic animals from 4 different breeds (Table 1) that represent different genetic backgrounds and which are exposed to different pathogen pressures aiming at finding polymorphisms that might interfere with the TLR3 function.
Results and discussion
According to Ensembl (http://www.ensembl.org; ENSOCUG00000017763), TLR3 is assigned to chromosome 2 of the European rabbit genome, with a total length of 11762 bp. The gene is organised into 4 or 5 exons with the coding sequence (CDS) comprising 2718 bp and 2655 bp, respectively. These transcripts translate into 905 and 884 amino acids, respectively, with the shorter transcript having a 21 amino acid deletion at exon 3 corresponding to amino acids 783–803. The analysis of the TLR3 structure revealed the presence of an ectodomain with a LRR-NT, twenty four LRRs and a LRR-CT, a transmembrane and TIR domains (Figure 2). This structure is in agreement to that of other TLRs.
Several studies addressing the TLR3 genetic variation showed low levels of diversity in many mammals. Indeed, in nine breeds of bovines, seven SNPs were found of which three were non-synonymous . In addition, in 84 male pigs from eleven breeds nine synonymous SNPs and six non-synonymous were detected  while in five horse breeds four SNPs were described of which two are non-synonymous . Here we detected 41 SNPs in the coding sequence of TLR3 of which 14 are non-synonymous (Additional file 1: Table S1). This discrepancy in values could be due to the inclusion of wild animals in the present study which were not considered in the species previously described, and which in turn might have contributed for the higher levels of diversity found here. In particular, the results showed higher diversity in the OCA populations where 33 polymorphic positions were observed, whereas only 22 were observed in OCC from IP, 18 in OCC from France and only 4 in the domestic breeds (Table 2). These results are further confirmed by the population genetic parameters summarised in Table 2. Even though no diagnostic positions were observed, some were exclusive to each population (Table 2). In addition, a gradient of genetic diversity loss in the same direction as the domestication process (from the Iberian populations to the French populations; Table 2) seems to occur as expected due to the bottleneck effect associated with this process . However, the Tajima D’s test did not reject the null hypothesis that the populations evolve under neutrality (Table 2).
Seventy-five haplotypes were recovered of which the majority (64/75 ≈ 85%) are rare occurring only once in heterozigoty and only two (haplotypes 39 and 40) are shared between different populations (see Additional file 2: Table S2). Indeed, haplotype 39 is common in OCC from IP and France and all domestic breeds, which is expected as these are the closest populations. Haplotype 40 is more restricted and shared only by Iberian OCC from Tarragona and France and the English Spot breed. It is also interesting to note that OCA populations do not share haplotypes with the other subspecies neither among the two populations which reveal a high degree of genetic variability across individuals. With respect to the domestic individuals, all specimen from French Lop, Argent Champagne and New Zealand are characterised by presenting a single haplotype while in the English Spot breed 3 different haplotypes are present. The estimation of haplotype diversity showed a pattern of loss of genetic diversity from the Iberian to the French populations (Table 3). As expected, the lowest values were observed for the domestic breeds. The higher variability in the wild populations may be due to the absence of man-mediated selection and a greater exposure to pathogens. Additionally, the low genetic diversity observed for the French populations when comparing with the IP populations might be explained because they are a subsample of the Iberian OCC which has recently colonised France [20, 25, 26].
Of the total number of polymorphisms, ≈ 34% (14/41) are non-synonymous substitutions that could be a source of protein variability that might interfere in the ability of TLR3 to recognise ligands and lead to resistance/susceptibility to viral diseases (Table 4). In particular, the variation observed at residues 480, 560 and 864 is unique to the European rabbit as these amino acids are generally conserved among mammals (data not shown). Although no data associate the observed substitutions with disease, it is likely that populations that are in contact with different pathogenic agents co-evolve with them and display more genetic diversity in order to confer resistance. Examples of two viral diseases responsible for high mortalities in domestic and wild European rabbits are RHDV and myxomatosis [27, 28]. This variation could however be just the result of stochastic events that occurred during evolution of the protein.
Regarding the location of the amino acid changes, most (11) are in the LRR molecules and the remaining (3) are distributed across the transmembrane and TIR domains (Figure 2). The higher variation in the LRRs is expected as these molecules are responsible for the direct interaction with PAMPs and should evolve to maintain recognition and trigger the innate immune response. The interaction between TLR3 and dsRNA occurs in two binding sites localised in the region comprising LRR-NT and LRR3 and LRR19-20 . Of the non-synonymous substitutions observed in the LRRs, three are located within these interaction points. Nevertheless, none of these alterations (p.I46V in the LRR-NT, p.L101S in LRR2 and p.K560E in LRR20; Table 4) is conserved among vertebrates (data not shown) or essential for recognition [29, 30]. LRR20, along with LRR21 and LRR22, are known for their importance in binding and signalling , but the p.L575F alteration observed in the LRR20 is conservative.
Other regions in the ectodomain have also been described as being important for ligand binding such as those between R65 to K163 and K330 to K493 that were described by Choe, Kelker et al. (2005) as a patch of positively charged residues . In this region, of the four amino acid changes observed in this study, only one (p.Q353H; Figure 2 and Table 4) results in changes in the polarity. The importance of the ectodomain in the PAMPs recognition as well as the role of the possible coevolution with the pathogens they recognise in the evolution of this domain has been a recurrent matter particularly shown in the high level of positive selection in this domain across mammal TLRs . This supports that the changes observed in the ectodomain of TLR3 might be the result of adaptation of TLR3 to maintain recognition of the evolving viral PAMPs. Nonetheless, none of the non-synonymous substitutions here detected is located within the sites detected previously as being under positive selection in mammals .
In the transmembrane domain, two amino acid changes were found, but both are conservative in terms of polarity and charge. In contrast, the alteration in the TIR domain, A864S (Figure 2 and Table 4), results in change of a hydrophobic residue for one hydrophilic corresponding to a possible site of N-glycosylation that may have consequences in the functionality of this domain. The TLR3 molecule is highly glycosylated in its ectodomain where there are 15 predicted N-glycosylation sites that affect ligand binding or control the access to the interaction surfaces . Ligand binding in TLR3 only occurs at the glycan-free surface . Despite this, the functional consequences of a possible N-glycosylation site in the TIR domain are unknown.
The biological consequences of each of the described alterations in the TLR3 molecule of the European rabbit also remains unknown, but further studies should be conducted to assess their impact since in this species only this TLR is available for detection of RNA viruses.
The differences in the TLR3 genetic variability detected between the wild rabbit populations and the domestic breeds reflects the domestication process of this species but might also result from the different exposure to pathogens. Indeed, while wild populations are usually in contact with a higher number of pathogens, domestic breeds are usually bred in highly controlled facilities. The higher number of non-synonymous substitutions found in the ectodomain, which corresponds to the binding surface to the PAMPs, further suggests a role of the host-pathogen interaction in promoting variability. In addition, the limited degree of variability found in the domestic breeds might compromise survival in face of a viral infection putting at risk the studied breeds.
In this study we analysed 20 individuals of the subspecies Oryctolagus cuniculus cuniculus from Spain (Tarragona (OCCIP_T) and Zaragoza (OCCIP_Z)) and 20 of the subspecies Oryctolagus cuniculus algirus from Portugal (Pancas (OCA_P) and Mertola (OCA_M)). These natural populations are located within the species’ original range. Twenty individuals from a French population (Cerisay (OCCF)) and 20 domestic animals (5 individuals per breed (DOM)) belonging to the French Lop, English Spot, Argent Champagne and New Zealand breeds were also studied (Table 1). DNA was extracted with the phenol-chloroform method . Amplification was carried out with Taq DNA Polymerase (Taq PCR Master Mix Kit, Qiagen) with several sets of primers (see Additional file 3: Table S3). The thermal profile consisted of a denaturation of 95°C for 15 min and then, 30 cycles of 95°C for 30s; optimal annealing temperature for 30s; 72°C for 45 s and a final amplification at 60°C for 20 min. PCR amplicons of the expected size were sequenced with BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City) in the following temperatures: 94°C for 3 min, 23 cycles at 96°C for 10s, 55°C for 5 s and 60°C for 4 min. The electropherogram analysis was performed in 3130xl Genetic Analyzer (Applied Biosystems, Foster City).
Sequences were edited, aligned and analysed with CLUSTAL W  as implemented in the BioEdit software . Haplotype reconstruction was obtained for each individual using PHASE (Version 2.1)  as implemented in DnaSP (Version 5.10.01) . Only haplotypes with P > 0.90 were considered as inferred reliably. Population genetic parameters (S, π), Tajima’s D, number of haplotypes, haplotype diversity, variance of haplotype diversity and standard deviation of haplotype diversity were estimated with DnaSP. The sequences obtained were also aligned and compared to the reference sequence available in Ensembl (ENSOCUT00000017763).
Availability of supporting data
Data obtained in this study was submitted to GenBank under the following accession numbers: KC963181 to KC963340.
Ethics for field studies
All sampling collection was conducted in accordance with local legislation and with the permissions and licences of the National institutions that supervise the hunting activities. The Convention on Biological Diversity and the Convention on the Trade in Endangered Species of Wild Fauna were respected.
Pathogen associated molecular pattern
Single nucleotide polymorphism
Pattern recognition receptor
Herpes simplex virus
West Nile virus
Oryctolagus cuniculus algirus
Oryctolagus cuniculus cuniculus
Takeuchi O, Akira S: Pattern recognition receptors and inflammation. Cell. 2010, 140 (6): 805-820. 10.1016/j.cell.2010.01.022.
Alexopoulou L, Holt AC, Medzhitov R, Flavell RA: Recognition of double-stranded RNA and activation of NF-kappaB by toll-like receptor 3. Nature. 2001, 413 (6857): 732-738. 10.1038/35099560.
Botos I, Liu L, Wang Y, Segal DM, Davies DR: The toll-like receptor 3:dsRNA signaling complex. Biochim Biophys Acta. 2009, 1789 (9–10): 667-674.
Bell JK, Botos I, Hall PR, Askins J, Shiloach J, Segal DM, Davies DR: The molecular structure of the toll-like receptor 3 ligand-binding domain. Proc Natl Acad Sci USA. 2005, 102 (31): 10976-10980. 10.1073/pnas.0505077102.
Bell JK, Mullen GE, Leifer CA, Mazzoni A, Davies DR, Segal DM: Leucine-rich repeats and pathogen recognition in toll-like receptors. Trends Immunol. 2003, 24 (10): 528-533. 10.1016/S1471-4906(03)00242-4.
Beutler B, Jiang Z, Georgel P, Crozat K, Croker B, Rutschmann S, Du X, Hoebe K: Genetic analysis of host resistance: toll-like receptor signaling and immunity at large. Annu Rev Immunol. 2006, 24: 353-389. 10.1146/annurev.immunol.24.021605.090552.
Medzhitov R, Janeway CA: Innate immunity: the virtues of a nonclonal system of recognition. Cell. 1997, 91 (3): 295-298. 10.1016/S0092-8674(00)80412-2.
Tabeta K, Georgel P, Janssen E, Du X, Hoebe K, Crozat K, Mudd S, Shamel L, Sovath S, Goode J, et al: Toll-like receptors 9 and 3 as essential components of innate immune defense against mouse cytomegalovirus infection. Proc Natl Acad Sci USA. 2004, 101 (10): 3516-3521. 10.1073/pnas.0400525101.
Zhang J, Ugolini S, Smahi A, Elain G, Romero P, Segal D, Sancho-Shimizu V, Lorenzo L, Puel A, et al: TLR3 Deficiency in patients with herpes simplex encephalitis. Science. 2007, 317 (5844): 1522-1527. 10.1126/science.1139522.
Wang T, Town T, Alexopoulou L, Anderson JF, Fikrig E, Flavell RA: Toll-like receptor 3 mediates West Nile virus entry into the brain causing lethal encephalitis. Nat Med. 2004, 10 (12): 1366-1373. 10.1038/nm1140.
Cargill EJ, Womack JE: Detection of polymorphisms in bovine toll-like receptors 3, 7, 8, and 9. Genomics. 2007, 89 (6): 745-755. 10.1016/j.ygeno.2007.02.008.
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.
Astakhova NM, Perelygin AA, Zharkikh AA, Lear TL, Coleman SJ, MacLeod JN, Brinton MA: Characterization of equine and other vertebrate TLR3, TLR7, and TLR8 genes. Immunogenetics. 2009, 61 (7): 529-539. 10.1007/s00251-009-0381-z.
Carneiro M, Ferrand N, Nachman MW: Recombination and speciation: loci near centromeres are more differentiated than loci near telomeres between subspecies of the European rabbit (Oryctolagus cuniculus). Genetics. 2009, 181 (2): 593-606.
Whitman B: Domestic rabbits and their histories. 2004, Leawodd (Kansas): Leathers Publishing
Carmo CR, Esteves PJ, Ferrand N, van der Loo W: Genetic variation at chemokine receptor CCR5 in leporids: alteration at the 2nd extracellular domain by gene conversion with CCR2 in Oryctolagus, but not in Sylvilagus and Lepus species. Immunogenetics. 2006, 58 (5-6): 494-501. 10.1007/s00251-006-0095-4.
Carneiro M, Blanco-Aguiar JA, Villafuerte R, Ferrand N, Nachman MW: Speciation in the European rabbit (Oryctolagus cuniculus): islands of differentiation on the X chromosome and autosomes. Evolution. 2010, 64 (12): 3443-3460. 10.1111/j.1558-5646.2010.01092.x.
Guillon P, Ruvoen-Clouet N, Le Moullac-Vaidye B, Marchandeau S, Le Pendu J: Association between expression of the H histo-blood group antigen, alpha1,2fucosyltransferases polymorphism of wild rabbits, and sensitivity to rabbit hemorrhagic disease virus. Glycobiology. 2009, 19 (1): 21-28.
Marchandeau S, Le Gall-Recule G, Bertagnoli S, Aubineau J, Botti G, Lavazza A: Serological evidence for a non-protective RHDV-like virus. Vet Res. 2005, 36 (1): 53-62. 10.1051/vetres:2004049.
Carneiro M, Afonso S, Geraldes A, Garreau H, Bolet G, Boucher S, Tircazes A, Queney G, Nachman MW, Ferrand N: The genetic structure of domestic rabbits. Mol Biol Evol. 2011, 28 (6): 1801-1816. 10.1093/molbev/msr003.
Letunic I, Doerks T, Bork P: SMART 6: recent updates and new developments. Nucleic acids research. 2009, 37 (Database issue): D229-232.
Morozumi T, Uenishi H: Polymorphism distribution and structural conservation in RNA-sensing toll-like receptors 3, 7, and 8 in pigs. Biochim Biophys Acta. 2009, 4: 267-274.
Nei M: Molecular evolutionary genetics. 1987, COLUMBIA: University Press
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.
Esteves Lanning D, Ferrand N, Knight KL, Zhai SK, Van der Loo W: Allelic variation at the VHa locus in natural populations of rabbit (Oryctolagus cuniculus, L.). J Immunol. 2004, 172 (2): 1044-1053.
Surridge AK, van der Loo W, Abrantes J, Carneiro M, Hewitt GM, Esteves PJ: Diversity and evolutionary history of the MHC DQA gene in leporids. Immunogenetics. 2008, 60 (9): 515-525. 10.1007/s00251-008-0309-z.
Kerr PJ: Myxomatosis in Australia and Europe: a model for emerging infectious diseases. Antiviral Res. 2012, 93 (3): 387-415. 10.1016/j.antiviral.2012.01.009.
Abrantes J, van der Loo W, Le Pendu J, Esteves PJ: Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review. Vet Res. 2012, 43 (1): 12-10.1186/1297-9716-43-12.
Liu L, Botos I, Wang Y, Leonard JN, Shiloach J, Segal DM, Davies DR: Structural basis of toll-like receptor 3 signaling with double-stranded RNA. Science. 2008, 320 (5874): 379-381. 10.1126/science.1155406.
Bell JK, Askins J, Hall PR, Davies DR, Segal DM: The dsRNA binding site of human toll-like receptor 3. Proc Natl Acad Sci USA. 2006, 103 (23): 8792-8797. 10.1073/pnas.0603245103.
Choe J, Kelker MS, Wilson IA: Crystal structure of human toll-like receptor 3 (TLR3) ectodomain. Science. 2005, 309 (5734): 581-585. 10.1126/science.1115253.
Areal H, Abrantes J, Esteves PJ: Signatures of positive selection in toll-like receptor (TLR) genes in mammals. BMC Evol Biol. 2011, 11 (1): 368-10.1186/1471-2148-11-368.
Offord V, Coffey TJ, Werling D: LRRfinder: a web application for the identification of leucine-rich repeats and an integrative toll-like receptor database. Dev Comp Immunol. 2010, 34 (10): 1035-1041. 10.1016/j.dci.2010.05.004.
Chomczynski P, Sacchi N: The single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction: twenty-something years on. Nat Protoc. 2006, 1 (2): 581-585. 10.1038/nprot.2006.83.
Thompson JD, Higgins DG, Gibson TJ, Clustal W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucl Acids Symp. 1999, 41: 95-98.
Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73 (5): 1162-1169. 10.1086/379378.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics (Oxford, England). 2009, 25 (11): 1451-1452. 10.1093/bioinformatics/btp187.
The Portuguese Foundation for Science and Technology supported the doctoral fellowship of Helena Areal (SFRH/BD/74948/2010) and the post-doctoral fellowship of Joana Abrantes (SFRH/BPD/73512/2010) as well as the project PTDC/CVT/108490/2008. We are grateful to the Fédération Française de Cuniculture, Union des Sélectionneurs Avicoles du Sud Ouest, Société central d’Aviculture de France, Union des Aviculteurs du Haut-Rhin, Groupement Avicole Pyrénées Atlantiques et Landes, Lapin Club du Sud-Ouest, J. Ruesche, D. Allain, and G. Saleil for providing the domestic rabbit samples.
The authors declare that they have no competing interests.
JA and PJE conceived of the study, participated in its design and coordination and helped to interpret data and draft the manuscript. HA carried out the experimental work and analysed the data. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Abrantes, J., Areal, H. & Esteves, P.J. Insights into the European rabbit (Oryctolagus cuniculus) innate immune system: genetic diversity of the toll-like receptor 3 (TLR3) in wild populations and domestic breeds. BMC Genet 14, 73 (2013). https://doi.org/10.1186/1471-2156-14-73
- Toll-like receptor 3 (TLR3)
- European rabbit (Oryctolagus cuniculus)
- RNA virus
- Genetic diversity