Polymorphism in pattern recognition receptor genes in pigs

Linnaeus University Dissertations No 8/2010

POLYMORPHISM IN PATTERN RECOGNITION RECEPTOR GENES IN PIGS

INGRID-MARIA BERGMAN VDM

LINNAEUS UNIVERSITY PRESS

Linnaeus University Dissertations No 8/2010

POLYMORPHISM IN PATTERN RECOGNITION RECEPTOR GENES IN PIGS

INGRID-MARIA BERGMAN VDM

LINNAEUS UNIVERSITY PRESS

Abstract

Bergman VDM, Ingrid-Maria (2010). Polymorphism in pattern recognition receptor genes in pigs Linnaeus University Dissertations No 8/2010. ISBN: 978-91-86491-10-9. Writ- ten in English. The mammalian immune defense consists of two systems, which are interconnected and co-operate to provide host defense. The innate immune system is always active and detects and responds to non-self without delay. The adaptive immune system has a lag phase, but is more specific and has got a memory. The innate immune system relies on pattern recognition receptors (PRRs) to detect molecular patterns signaling microbial presence. This thesis focuses on a centrally placed family of PRRs, namely the Toll-like receptors (TLRs), and on mannan-binding (MBL), a PRR which initiates the lectin activation pathway of complement. TLRs are expressed on the cell surface and in intracellular compartments, while MBL is a soluble protein present in most body fluids. Polymorphism – literally ’many forms’ – refers to variation between individuals, at DNA level as well as in traits. A single nucleotide polymorphism (SNP) implicates that alternative nucleotides are present at a particular position in the genome. Mutations, to- gether with phenomena like gene duplication and whole genome duplication, are the ultimate source of variation in nature and the fuel for evolution. Through natural selec- tion and breeding, i.e. artificial selection, species are shaped and change over time. Domestic animals are well suited for genetic studies, since they enable comparisons of populations exposed to different selection criteria and environmental challenges. Also, in the case of pigs, comparisons to the wild ancestor – i.e. the wild boar – can shed light on the evolutionary process. Moreover, pigs are large animal models for hu- mans. Paper I reports the refinement of previously identified quantitative trait loci for im- mune-related traits on pig chromosome 8. Papers II and III report differences in polymorphic patterns between wild boars and domestic pigs in the TLR1, TLR2, TLR6, and TLR10 genes. In TLR1 and TLR2, more SNPs were present in the domestic pigs than in the wild boars. In TLR6, SNP numbers were similar in both animal groups, but the level of heterozygosity was higher in the domestic pigs than in the wild boars. In TLR10, again, more SNPs were present in the domestic pigs, and a higher number of non-synonymous SNPs was detected in TLR10 compared to the other genes. This might suggest redundancy for TLR10 in pigs. Paper IV reports the presence of an SNP, previously detected in domestic pigs and assumed to affect MBL concentrations in serum, in European wild boars. Also, the POLYMORPHISM IN PATTERN RECOGNITION RECEPTOR GENES connection between the presumed low-producing allele and low MBL concentration in IN PIGS serum was confirmed. Moreover, a novel SNP, with potential to be functionally impor- Doctoral dissertation, School of Natural Sciences, Linnaeus University 2010. tant, was detected. Owing to the domestication process and differences in selection pressure, differ- ences in polymorphic patterns between wild boars and domestic pigs are not surprising. However, since breeding means choosing among genotypes, the opposite pattern – Series editor: Kerstin Brodén more SNPs in wild boars than in domestic pigs – would have been expected. However, ISBN: 978-91-86491-10-9 the result confirms other studies, which have shown that European wild boars went Printed by: Intellecta Infolog, Göteborg 2010 through a bottle neck before domestication started. The higher number of SNPs in domestic pigs may be due to relaxed purifying selection during the domestication proc- ess.

Abstract

Bergman VDM, Ingrid-Maria (2010). Polymorphism in pattern recognition receptor genes in pigs Linnaeus University Dissertations No 8/2010. ISBN: 978-91-86491-10-9. Writ- ten in English. The mammalian immune defense consists of two systems, which are interconnected and co-operate to provide host defense. The innate immune system is always active and detects and responds to non-self without delay. The adaptive immune system has a lag phase, but is more specific and has got a memory. The innate immune system relies on pattern recognition receptors (PRRs) to detect molecular patterns signaling microbial presence. This thesis focuses on a centrally placed family of PRRs, namely the Toll-like receptors (TLRs), and on mannan-binding lectin (MBL), a PRR which initiates the lectin activation pathway of complement. TLRs are expressed on the cell surface and in intracellular compartments, while MBL is a soluble protein present in most body fluids. Polymorphism – literally ’many forms’ – refers to variation between individuals, at DNA level as well as in traits. A single nucleotide polymorphism (SNP) implicates that alternative nucleotides are present at a particular position in the genome. Mutations, to- gether with phenomena like gene duplication and whole genome duplication, are the ultimate source of variation in nature and the fuel for evolution. Through natural selec- tion and breeding, i.e. artificial selection, species are shaped and change over time. Domestic animals are well suited for genetic studies, since they enable comparisons of populations exposed to different selection criteria and environmental challenges. Also, in the case of pigs, comparisons to the wild ancestor – i.e. the wild boar – can shed light on the evolutionary process. Moreover, pigs are large animal models for hu- mans. Paper I reports the refinement of previously identified quantitative trait loci for im- mune-related traits on pig chromosome 8. Papers II and III report differences in polymorphic patterns between wild boars and domestic pigs in the TLR1, TLR2, TLR6, and TLR10 genes. In TLR1 and TLR2, more SNPs were present in the domestic pigs than in the wild boars. In TLR6, SNP numbers were similar in both animal groups, but the level of heterozygosity was higher in the domestic pigs than in the wild boars. In TLR10, again, more SNPs were present in the domestic pigs, and a higher number of non-synonymous SNPs was detected in TLR10 compared to the other genes. This might suggest redundancy for TLR10 in pigs. Paper IV reports the presence of an SNP, previously detected in domestic pigs and assumed to affect MBL concentrations in serum, in European wild boars. Also, the POLYMORPHISM IN PATTERN RECOGNITION RECEPTOR GENES connection between the presumed low-producing allele and low MBL concentration in IN PIGS serum was confirmed. Moreover, a novel SNP, with potential to be functionally impor- Doctoral dissertation, School of Natural Sciences, Linnaeus University 2010. tant, was detected. Owing to the domestication process and differences in selection pressure, differ- ences in polymorphic patterns between wild boars and domestic pigs are not surprising. However, since breeding means choosing among genotypes, the opposite pattern – Series editor: Kerstin Brodén more SNPs in wild boars than in domestic pigs – would have been expected. However, ISBN: 978-91-86491-10-9 the result confirms other studies, which have shown that European wild boars went Printed by: Intellecta Infolog, Göteborg 2010 through a bottle neck before domestication started. The higher number of SNPs in domestic pigs may be due to relaxed purifying selection during the domestication proc- ess.

Till Martin

Till Martin

LIST OF PAPERS

This thesis is based on the following papers, referred to in the text by their numbers. Paper II was reprinted with permission from the publisher.

Paper I Bergman IM, Johansson A, Wattrang E, Fossum C, Andersson

L, Edfors I (2010) Refined analysis of quantitative trait loci

(QTLs) for immune capacity related traits on pig chromosome

8. Manuscript

Paper II Bergman IM, Rosengren JK, Edman K, Edfors I (2010) Euro-

pean wild boars and domestic pigs display different polymor-

phic patterns in the Toll-like receptor (TLR) 1, TLR2, and

TLR6 genes. Immunogenetics 62:49-58. doi:10.1007/s00251-

009-0409-4

Paper III Bergman IM, Edman K, Rosengren JK, Edfors I (2010) Single

nucleotide polymorphisms (SNPs) in the Toll-like receptor

(TLR) 10 gene in European wild boars and domestic pigs.

Submitted manuscript

Paper IV Bergman IM, Sandholm K, Juul-Madsen HR, Heegaard PM,

Nilsson Ekdahl K, Edfors I (2010) MBL-A concentrations and

MBL1 genotypes in European wild boars, Large White pigs,

and wild boar/Large White crossbreds. Manuscript

3

LIST OF PAPERS

This thesis is based on the following papers, referred to in the text by their numbers. Paper II was reprinted with permission from the publisher.

Paper I Bergman IM, Johansson A, Wattrang E, Fossum C, Andersson

L, Edfors I (2010) Refined analysis of quantitative trait loci

(QTLs) for immune capacity related traits on pig chromosome

8. Manuscript

Paper II Bergman IM, Rosengren JK, Edman K, Edfors I (2010) Euro-

pean wild boars and domestic pigs display different polymor-

phic patterns in the Toll-like receptor (TLR) 1, TLR2, and

TLR6 genes. Immunogenetics 62:49-58. doi:10.1007/s00251-

009-0409-4

Paper III Bergman IM, Edman K, Rosengren JK, Edfors I (2010) Single

nucleotide polymorphisms (SNPs) in the Toll-like receptor

(TLR) 10 gene in European wild boars and domestic pigs.

Submitted manuscript

Paper IV Bergman IM, Sandholm K, Juul-Madsen HR, Heegaard PM,

Nilsson Ekdahl K, Edfors I (2010) MBL-A concentrations and

MBL1 genotypes in European wild boars, Large White pigs,

and wild boar/Large White crossbreds. Manuscript

3

ABBREVIATIONS

APC antigen-presenting cell

BC2 second backcross generation

BCR B cell receptor

bp base pairs

C cystein-rich N-terminal domain

C1qRp C1q receptor for phagocytosis enhancement

C3aR C3a receptor

C5aR C5a receptor

CD cluster of differentiation

Caco-2 human cancer cell line

CLD collagen-like domain

cM centimorgan

ConA , a lymphocyte

CR3 complement receptor-3

CRD carbohydrate recognition domain

CRP C reactive protein

DC dendritic cell

DNA deoxyribonucleic acid

F1 first filial generation

F2 second filial generation

FLS2 flagellin-sensitive 2

Hb hemoglobin

HBB beta hemo globin gene

HbS the mutant HBB allele, causing sickle cell ane-

mia

HCA-7 human cancer cell line

Hem hematocrit

IBD identity-by-descent

Ig immunoglobulin IGF2 insulin-like growth factor 2 IL interleukin IFN interferon KIT mast/stem cell growth factor receptor LOD log of odds LPS lipopolysaccharide LRR leucin-rich repeat Mal MyD88-adapter like MAMP microorganism-associated molecular pattern MAP-1 MBL/ associated protein 1

4 5

ABBREVIATIONS

APC antigen-presenting cell

BC2 second backcross generation

BCR B cell receptor

bp base pairs

C cystein-rich N-terminal domain

C1qRp C1q receptor for phagocytosis enhancement

C3aR C3a receptor

C5aR C5a receptor

CD cluster of differentiation

Caco-2 human cancer cell line

CLD collagen-like domain

cM centimorgan

ConA concanavalin A, a lymphocyte mitogen

CR3 complement receptor-3

CRD carbohydrate recognition domain

CRP C reactive protein

DC dendritic cell

DNA deoxyribonucleic acid

F1 first filial generation

F2 second filial generation

FLS2 flagellin-sensitive 2

Hb hemoglobin

HBB beta hemo globin gene

HbS the mutant HBB allele, causing sickle cell ane-

mia

HCA-7 human cancer cell line

Hem hematocrit

IBD identity-by-descent

Ig immunoglobulin IGF2 insulin-like growth factor 2 IL interleukin IFN interferon KIT mast/stem cell growth factor receptor LOD log of odds LPS lipopolysaccharide LRR leucin-rich repeat Mal MyD88-adapter like MAMP microorganism-associated molecular pattern MAP-1 MBL/ficolin associated protein 1

4 5 MASP mannan-binding protein associated serine pro- TNF tumor necrosis factor tease TLR Toll-like receptor MHC I and II major histocompatibility complexes I and II TRAM TRIF-related adapter molecule MyD88 myeloid differentiation primary response pro- TRIF TIR domain-containing adapter inducing inter- tein 88 feron-β MBL mannan-binding lectin V-TLR vertebrate-like TLR MC1R melanocortin 1 receptor ms microsatellite N neck region Neu segmented neutrophils NK natural killer NFκB nuclear factor κB NOD 1 and 2 nucleotide-binding oligomerization domain- containing proteins 1 and 2 P parental generation p short chromosome arm PCR polymerase chain reaction PHA PKR protein kinase R PRR pattern recognition receptor PTX3 pentraxin 3 P-TLR protostome-like TLR PWM pokeweed mitogen q long chromosome arm OTG quantitative trait gene QTL quantitative trait loci QTN quantitative trait nucleotide R resistance RAG 1 and 2 recombination activating genes 1 and 2 REL protein domain found in a family of transcrip- tion factors, also known as RHD, rel homology domain RNA ribonucleic acid s a sine annu (without year) SARM sterile α and HEAT-Armadillo motifs SAP serum amyloid P component SSC8 Sus scrofa (pig) chromosome 8 SMART simple modular architecture research tool SNP single nucleotide polymorphism SP-A and -D surfactant proteins A and -D TCR T cell receptor TIR Toll/interleukin 1 receptor homology domain TICAM-1 and -2 TIR-containing adapter molecule-1 and -2 TIRAP TIR domain-containing adapter

6 7 MASP mannan-binding protein associated serine pro- TNF tumor necrosis factor tease TLR Toll-like receptor MHC I and II major histocompatibility complexes I and II TRAM TRIF-related adapter molecule MyD88 myeloid differentiation primary response pro- TRIF TIR domain-containing adapter inducing inter- tein 88 feron-β MBL mannan-binding lectin V-TLR vertebrate-like TLR MC1R melanocortin 1 receptor ms microsatellite N neck region Neu segmented neutrophils NK natural killer NFκB nuclear factor κB NOD 1 and 2 nucleotide-binding oligomerization domain- containing proteins 1 and 2 P parental generation p short chromosome arm PCR polymerase chain reaction PHA phytohaemagglutinin PKR protein kinase R PRR pattern recognition receptor PTX3 pentraxin 3 P-TLR protostome-like TLR PWM pokeweed mitogen q long chromosome arm OTG quantitative trait gene QTL quantitative trait loci QTN quantitative trait nucleotide R resistance RAG 1 and 2 recombination activating genes 1 and 2 REL protein domain found in a family of transcrip- tion factors, also known as RHD, rel homology domain RNA ribonucleic acid s a sine annu (without year) SARM sterile α and HEAT-Armadillo motifs SAP serum amyloid P component SSC8 Sus scrofa (pig) chromosome 8 SMART simple modular architecture research tool SNP single nucleotide polymorphism SP-A and -D surfactant proteins A and -D TCR T cell receptor TIR Toll/interleukin 1 receptor homology domain TICAM-1 and -2 TIR-containing adapter molecule-1 and -2 TIRAP TIR domain-containing adapter

6 7

TABLE OF CONTENTS

Sammanfattning

General introduction

Part I: Immune defense 1. IMMUNE DEFENSE IN PLANTS AND IN DROSOPHILA MELANOGASTER 2. THE MAMMALIAN IMMUNE DEFENSE 2.1 The innate immune defense 2.1.1 Epithelial barriers, effector cells, defense proteins, and cytokines 2.1.2 Recognition strategies 2.1.3 The complement system 2.1.4 Resistance against innate immunity and maintenance of homeostasis 2.2 Main features of the adaptive immune defense 2.3 Bridges between the innate and adaptive immune systems

Part II: Pattern recognition receptors 3. SELECTED PATTERN RECOGNITION RECEPTORS AND THEIR FUNCTIONS 4. THE TOLL-LIKE RECEPTORS 4.1 Evolutionary perspective 4.2 The TLR molecule 4.3 TLR signaling 4.4 Comparisons to Drosophila Toll and plant TIR 5. MANNAN-BINDING LECTIN

Part III: Quantitative traits and QTL analysis 6. QUANTITATIVE TRAITS 7. QTL ANALYSIS

8 9

TABLE OF CONTENTS

Sammanfattning

General introduction

Part I: Immune defense 1. IMMUNE DEFENSE IN PLANTS AND IN DROSOPHILA MELANOGASTER 2. THE MAMMALIAN IMMUNE DEFENSE 2.1 The innate immune defense 2.1.1 Epithelial barriers, effector cells, defense proteins, and cytokines 2.1.2 Recognition strategies 2.1.3 The complement system 2.1.4 Resistance against innate immunity and maintenance of homeostasis 2.2 Main features of the adaptive immune defense 2.3 Bridges between the innate and adaptive immune systems

Part II: Pattern recognition receptors 3. SELECTED PATTERN RECOGNITION RECEPTORS AND THEIR FUNCTIONS 4. THE TOLL-LIKE RECEPTORS 4.1 Evolutionary perspective 4.2 The TLR molecule 4.3 TLR signaling 4.4 Comparisons to Drosophila Toll and plant TIR 5. MANNAN-BINDING LECTIN

Part III: Quantitative traits and QTL analysis 6. QUANTITATIVE TRAITS 7. QTL ANALYSIS

8 9 Part IV: Mutations, selection, and evolution Acknowledgements 8. MUTATIONS 9. SELECTION AND EVOLUTION References 10. WILD BOARS AND DOMESTIC PIGS

Part V: Materials and methods 11. PIGS IN GENETIC ANALYSES 12. METHODOLOGICAL CONSIDERATIONS 12.1 Paper I 12.1.1 Duplex PCR reactions using a touch-down procedure 12.1.2 Microsatellites as markers 12.1.3 Interpretation of signals 12.1.4 CRI-MAP 12.1.5 QTL Express 12.1.6 Limitations in connection with QTL analysis 12.2 Papers II and III 12.2.1 Direct sequencing 12.2.2 Protein prediction 12.2.3 Mann-Whitney U-test 12.2.4 The dN/dS statistic 12.3 Paper IV 12.3.1 Sequencing subsequent to cloning 12.3.2 ELISA

Aims

Part VI: Papers, conclusions, future 13. PAPER I 13.1 Objectives 13.2 Results and discussion 14. PAPERS II AND III 14.1 Objectives 14.2 Results and discussion 15. PAPER IV 15.1 Objectives 15.2 Results and discussion 16. CONCLUSIONS 17. FUTURE

10 11 Part IV: Mutations, selection, and evolution Acknowledgements 8. MUTATIONS 9. SELECTION AND EVOLUTION References 10. WILD BOARS AND DOMESTIC PIGS

Part V: Materials and methods 11. PIGS IN GENETIC ANALYSES 12. METHODOLOGICAL CONSIDERATIONS 12.1 Paper I 12.1.1 Duplex PCR reactions using a touch-down procedure 12.1.2 Microsatellites as markers 12.1.3 Interpretation of signals 12.1.4 CRI-MAP 12.1.5 QTL Express 12.1.6 Limitations in connection with QTL analysis 12.2 Papers II and III 12.2.1 Direct sequencing 12.2.2 Protein prediction 12.2.3 Mann-Whitney U-test 12.2.4 The dN/dS statistic 12.3 Paper IV 12.3.1 Sequencing subsequent to cloning 12.3.2 ELISA

Aims

Part VI: Papers, conclusions, future 13. PAPER I 13.1 Objectives 13.2 Results and discussion 14. PAPERS II AND III 14.1 Objectives 14.2 Results and discussion 15. PAPER IV 15.1 Objectives 15.2 Results and discussion 16. CONCLUSIONS 17. FUTURE

10 11 Sammanfattning

Däggdjurs immunförsvar består av två system, ett äldre och ett yngre, som är sammanflätade och samarbetar för att försvara värden mot mikroorganismer. Det medfödda immunförsvaret är alltid aktivt och reagerar omedelbart på när- varon av icke-själv. Det evolutionärt yngre förvärvade immunförsvaret behöver tid för att komma igång, men är i gengäld mer specifikt och har också ett min- ne. Detta minne är den mekanism som utnyttjas vid vaccinering.

Det medfödda immunförsvarets funktion är beroende av receptorer som känner igen molekylära mönster som är typiska för mikroorganismer men inte förekommer hos värden. Den här avhandlingen fokuserar på en centralt place- rad familj av sådana receptorer, nämligen de Toll-lika receptorerna (TLR), och mannan-bindande lectin (MBL), som hör till en annan proteinfamilj men har samma slags funktion. TLR förekommer på cellytan och i vissa av cellens orga- neller, medan MBL är ett lösligt protein som finns i de flesta kroppsvätskor.

Ordet polymorfi betyder ’många former’ och betecknar variation mellan in- divider, dels på gennivå och dels i egenskaper. En SNP (’Single Nucleotide Po- lymorphism’) är variation på en enstaka position i genomet. Mutationer, till- sammans med fenomen som genduplikation och genomduplikation, är den yt- tersta orsaken till variation i naturen och grunden för evolution. Genom natur- ligt urval och avel, det vill säga artificiellt urval, formas och förändras arter över tid.

Husdjur är väl lämpade för genetiska studier, eftersom det är möjligt att jämföra populationer med olika bakgrund och göra försökskorsningar. När det gäller grisar kan jämförelser med vildsvinet dessutom kasta ljus över den evolu- tionära processen. Grisen är också en vanlig stordjursmodell för att studera sjukdomar hos människor.

Den första studien i avhandlingen är en QTL-analys. Med hjälp av en så- dan kan man avgöra var någonstans i genomet det finns gener som påverkar kvantitativa egenskaper, som t ex immunkapacitet.

Arbetet med TLR-generna visade att vildsvin och tamgrisar har olika poly- morfa mönster i generna som kodar för TLR1, TLR2, TLR6 och TLR10. An- talet SNPer i TLR1 och TLR2 var högre hos tamgrisar än hos vildsvin. I TLR6 var antalet SNPer ungefär det samma i båda djurgrupperna, men heterozygo- tigraden – dvs. förekomsten av båda versionerna av SNPn på den aktuella posi- tionen – var högre i tamgrisgruppen än i vildsvinsgruppen. Även i TLR10 fanns det fler SNPer hos tamgrisarna än hos vildsvinen, men antalet icke-synonyma (betydelsebärande) SNPer var dessutom fler i TLR10 jämfört med de andra ge- nerna. Detta skulle kunna tyda på att TLR10 inte är en nödvändig receptor,

12 13 Sammanfattning

Däggdjurs immunförsvar består av två system, ett äldre och ett yngre, som är sammanflätade och samarbetar för att försvara värden mot mikroorganismer. Det medfödda immunförsvaret är alltid aktivt och reagerar omedelbart på när- varon av icke-själv. Det evolutionärt yngre förvärvade immunförsvaret behöver tid för att komma igång, men är i gengäld mer specifikt och har också ett min- ne. Detta minne är den mekanism som utnyttjas vid vaccinering.

Det medfödda immunförsvarets funktion är beroende av receptorer som känner igen molekylära mönster som är typiska för mikroorganismer men inte förekommer hos värden. Den här avhandlingen fokuserar på en centralt place- rad familj av sådana receptorer, nämligen de Toll-lika receptorerna (TLR), och mannan-bindande lectin (MBL), som hör till en annan proteinfamilj men har samma slags funktion. TLR förekommer på cellytan och i vissa av cellens orga- neller, medan MBL är ett lösligt protein som finns i de flesta kroppsvätskor.

Ordet polymorfi betyder ’många former’ och betecknar variation mellan in- divider, dels på gennivå och dels i egenskaper. En SNP (’Single Nucleotide Po- lymorphism’) är variation på en enstaka position i genomet. Mutationer, till- sammans med fenomen som genduplikation och genomduplikation, är den yt- tersta orsaken till variation i naturen och grunden för evolution. Genom natur- ligt urval och avel, det vill säga artificiellt urval, formas och förändras arter över tid.

Husdjur är väl lämpade för genetiska studier, eftersom det är möjligt att jämföra populationer med olika bakgrund och göra försökskorsningar. När det gäller grisar kan jämförelser med vildsvinet dessutom kasta ljus över den evolu- tionära processen. Grisen är också en vanlig stordjursmodell för att studera sjukdomar hos människor.

Den första studien i avhandlingen är en QTL-analys. Med hjälp av en så- dan kan man avgöra var någonstans i genomet det finns gener som påverkar kvantitativa egenskaper, som t ex immunkapacitet.

Arbetet med TLR-generna visade att vildsvin och tamgrisar har olika poly- morfa mönster i generna som kodar för TLR1, TLR2, TLR6 och TLR10. An- talet SNPer i TLR1 och TLR2 var högre hos tamgrisar än hos vildsvin. I TLR6 var antalet SNPer ungefär det samma i båda djurgrupperna, men heterozygo- tigraden – dvs. förekomsten av båda versionerna av SNPn på den aktuella posi- tionen – var högre i tamgrisgruppen än i vildsvinsgruppen. Även i TLR10 fanns det fler SNPer hos tamgrisarna än hos vildsvinen, men antalet icke-synonyma (betydelsebärande) SNPer var dessutom fler i TLR10 jämfört med de andra ge- nerna. Detta skulle kunna tyda på att TLR10 inte är en nödvändig receptor,

12 13 utan kan förändras i högre grad än de övriga utan negativa konsekvenser för General introduction grisen. The topic of this thesis is polymorphic patterns in genes encoding pattern MBL-studien visar att en SNP, som tidigare detekterats i tamgrisar och recognition receptors (PRRs) in pigs. Polymorphisms in immune-related genes förmodligen har stor betydelse för koncentrationen av MBL i blodet, också fö- commonly affect the immune competence of the host and are therefore of grea- rekommer hos vildsvin. Sambandet mellan den version av SNPn som antas ge test interest. Domestic animals are well suited for genetic studies, since they låga MBL-koncentrationer och låga mätvärden hos individen kunde också be- enable comparisons of populations exposed to different selection criteria and kräftas. Dessutom detekterades en tidigare okänd SNP, som skulle kunna för- environmental challenges. Also, in the case of pigs, comparisons to the wild an- väntas ge upphov till ett protein med förändrad funktion. cestor – i.e. the wild boar – can shed light on the evolutionary process. More- over, pigs are large animal models for humans. Eftersom vildsvin och tamgrisar har olika bakgrund är det inte förvånande att de uppvisar olika polymorfa mönster i de analyserade generna. Däremot The immune defense is dependent on evolutionarily conserved PRRs to de- kunde mönstret ha förväntats vara det motsatta – fler SNPer hos vildsvin än tect the presence of invading microorganisms and altered self. These receptors hos tamgrisar – eftersom avel innebär att vissa genotyper väljs bort. Resultatet are present on a variety of cell types and exert a variety of functions. The focus bekräftar dock andra studier, som visar att de europeiska vildsvinen gick ige- of the work behind this thesis has been on a centrally placed family of PRRs, nom en evolutionär flaskhals – dvs. under en period fanns det mycket få indivi- namely the Toll-like receptors (TLRs), and mannan-binding lectin (MBL), der – innan domesticeringen startade. Evolutionära flaskhalsar resulterar i låg also a PRR, which induces the activity of a system of defense proteins. grad av variation. Det högre antalet SNPer i tamgrisar skulle kunna bero på att den renande selektionen, som skyddar ett proteins funktion, har arbetat mindre Traits influenced by several genes exhibit quantitative variation. In order to effektivt under den tid då grisen varit husdjur. determine which genes should be focused on when studying a particular quan- titative trait of interest, e.g. immune competence, a quantitative trait loci (QTL) analysis can be applied. Moreover, in order to interpret genetic snap- shots in time, like those presented in this thesis, insights concerning selective pressure and the evolutionary process are necessary.

Part I of this thesis focuses on immune defense. Paragraph 1 is concerned with the immune defense in plants and in Drosophila melanogaster, the fruit fly. These two ancient systems put the mammalian immune defense into perspec- tive and simplify the understanding of many of its features. Paragraph 2 de- scribes the most central features of the mammalian immune defense, focusing on the innate immune system and its main function: the discrimination be- tween self and non-self.

Part II describes some PRRs and their genes, focusing on TLRs and MBL.

Part III discusses quantitative traits and QTL analysis, part IV is focused on mutations, selection, and evolution, while part V is concerned with the us- age of pigs in this work and in genetic studies in general and with methodo- logical considerations.

Part VI introduces the papers included in the thesis, presents conclusions drawn from the work, and outlines possible future research themes.

14 15 utan kan förändras i högre grad än de övriga utan negativa konsekvenser för General introduction grisen. The topic of this thesis is polymorphic patterns in genes encoding pattern MBL-studien visar att en SNP, som tidigare detekterats i tamgrisar och recognition receptors (PRRs) in pigs. Polymorphisms in immune-related genes förmodligen har stor betydelse för koncentrationen av MBL i blodet, också fö- commonly affect the immune competence of the host and are therefore of grea- rekommer hos vildsvin. Sambandet mellan den version av SNPn som antas ge test interest. Domestic animals are well suited for genetic studies, since they låga MBL-koncentrationer och låga mätvärden hos individen kunde också be- enable comparisons of populations exposed to different selection criteria and kräftas. Dessutom detekterades en tidigare okänd SNP, som skulle kunna för- environmental challenges. Also, in the case of pigs, comparisons to the wild an- väntas ge upphov till ett protein med förändrad funktion. cestor – i.e. the wild boar – can shed light on the evolutionary process. More- over, pigs are large animal models for humans. Eftersom vildsvin och tamgrisar har olika bakgrund är det inte förvånande att de uppvisar olika polymorfa mönster i de analyserade generna. Däremot The immune defense is dependent on evolutionarily conserved PRRs to de- kunde mönstret ha förväntats vara det motsatta – fler SNPer hos vildsvin än tect the presence of invading microorganisms and altered self. These receptors hos tamgrisar – eftersom avel innebär att vissa genotyper väljs bort. Resultatet are present on a variety of cell types and exert a variety of functions. The focus bekräftar dock andra studier, som visar att de europeiska vildsvinen gick ige- of the work behind this thesis has been on a centrally placed family of PRRs, nom en evolutionär flaskhals – dvs. under en period fanns det mycket få indivi- namely the Toll-like receptors (TLRs), and mannan-binding lectin (MBL), der – innan domesticeringen startade. Evolutionära flaskhalsar resulterar i låg also a PRR, which induces the activity of a system of defense proteins. grad av variation. Det högre antalet SNPer i tamgrisar skulle kunna bero på att den renande selektionen, som skyddar ett proteins funktion, har arbetat mindre Traits influenced by several genes exhibit quantitative variation. In order to effektivt under den tid då grisen varit husdjur. determine which genes should be focused on when studying a particular quan- titative trait of interest, e.g. immune competence, a quantitative trait loci (QTL) analysis can be applied. Moreover, in order to interpret genetic snap- shots in time, like those presented in this thesis, insights concerning selective pressure and the evolutionary process are necessary.

Part I of this thesis focuses on immune defense. Paragraph 1 is concerned with the immune defense in plants and in Drosophila melanogaster, the fruit fly. These two ancient systems put the mammalian immune defense into perspec- tive and simplify the understanding of many of its features. Paragraph 2 de- scribes the most central features of the mammalian immune defense, focusing on the innate immune system and its main function: the discrimination be- tween self and non-self.

Part II describes some PRRs and their genes, focusing on TLRs and MBL.

Part III discusses quantitative traits and QTL analysis, part IV is focused on mutations, selection, and evolution, while part V is concerned with the us- age of pigs in this work and in genetic studies in general and with methodo- logical considerations.

Part VI introduces the papers included in the thesis, presents conclusions drawn from the work, and outlines possible future research themes.

14 15

PART I: IMMUNE DEFENSE

The integrity of a single cell is maintained by the cell membrane. Similarly,

multicellular organisms rely on physical and chemical barriers between them-

selves and the environment, to uphold their integrity and prevent the entry of

microorganisms. However, microorganisms frequently manage to bypass these

barriers. At these occasions, more elaborate defense mechanisms come into

play.

1. Immune defense in plants and in Drosophila melano- gaster Plants have a two-layered germ line encoded innate immune defense. The primary immune system relies on PRRs, which recognize microorganism- associated molecular patterns (MAMPs) and induce defense responses, such as cell wall alterations and accumulation of defense proteins. However, plant pathogens evade the primary immune system by means of effectors that enable them to cause disease. In response to this, the host further protects itself thro- ugh a secondary immune system relying on resistance (R) proteins. Further- more, plants have a mechanism known as systemic acquired resistance, which protects from subsequent pathogenic attacks (de Wit 2007).

Drosophila melanogaster, the fruit fly, possesses a systemic immune response which is based on antimicrobial peptides secreted from the fat body into the hemolymph subsequent to infection. The epithelial barrier is armed with effi- cient defense systems, including a lining of chitinous matrix, secretion of ly- sozymes, and local production of reactive oxygen species and antimicrobial pep- tides. Moreover, specialized hemocytes with phagocytic and encapsulating ca- pacity as well as clotting and melanization play important roles in host defense in Drosophila (Lemaitre and Hoffmann 2007).

The Drosophila Toll gene was originally discovered as a maternal-effect gene: females lacking Toll gene activity produce embryos in which all cells be- have like the dorsal cells of the wild-type embryo (Anderson et al. 1985). How- ever, Toll also plays a role in synaptogenesis and axon path finding and is es- sential in the immune defense against fungi and Gram-positive bacteria (Le- maitre 2004; Rose et al. 1997). Signaling through Toll subsequent to infection is initiated by the cleavage of the cytokine Spätzle, which is a ligand for Toll (Lemaitre 2004). Thus, Drosophila Toll does not interact directly with MAMPs and is not a PRR. Activation of Toll leads – via an intracellular path-

16 17

PART I: IMMUNE DEFENSE

The integrity of a single cell is maintained by the cell membrane. Similarly,

multicellular organisms rely on physical and chemical barriers between them-

selves and the environment, to uphold their integrity and prevent the entry of

microorganisms. However, microorganisms frequently manage to bypass these

barriers. At these occasions, more elaborate defense mechanisms come into

play.

1. Immune defense in plants and in Drosophila melano- gaster Plants have a two-layered germ line encoded innate immune defense. The primary immune system relies on PRRs, which recognize microorganism- associated molecular patterns (MAMPs) and induce defense responses, such as cell wall alterations and accumulation of defense proteins. However, plant pathogens evade the primary immune system by means of effectors that enable them to cause disease. In response to this, the host further protects itself thro- ugh a secondary immune system relying on resistance (R) proteins. Further- more, plants have a mechanism known as systemic acquired resistance, which protects from subsequent pathogenic attacks (de Wit 2007).

Drosophila melanogaster, the fruit fly, possesses a systemic immune response which is based on antimicrobial peptides secreted from the fat body into the hemolymph subsequent to infection. The epithelial barrier is armed with effi- cient defense systems, including a lining of chitinous matrix, secretion of ly- sozymes, and local production of reactive oxygen species and antimicrobial pep- tides. Moreover, specialized hemocytes with phagocytic and encapsulating ca- pacity as well as clotting and melanization play important roles in host defense in Drosophila (Lemaitre and Hoffmann 2007).

The Drosophila Toll gene was originally discovered as a maternal-effect gene: females lacking Toll gene activity produce embryos in which all cells be- have like the dorsal cells of the wild-type embryo (Anderson et al. 1985). How- ever, Toll also plays a role in synaptogenesis and axon path finding and is es- sential in the immune defense against fungi and Gram-positive bacteria (Le- maitre 2004; Rose et al. 1997). Signaling through Toll subsequent to infection is initiated by the cleavage of the cytokine Spätzle, which is a ligand for Toll (Lemaitre 2004). Thus, Drosophila Toll does not interact directly with MAMPs and is not a PRR. Activation of Toll leads – via an intracellular path-

16 17 way involving myeloid differentiation primary response protein 88 (MyD88) cells and have different areas of responsibility: TNF, IL-1, and chemokines and other proteins – to translocation of Dorsal-related immunity factor to the mediate inflammation, IFN-γ activates macrophages, while IL-15 and IL-12 nucleus, where it regulates genes coding for antimicrobial peptides. Similarly, stimulate proliferation and activity of NK cells. IL-12 is also a key inducer of the immune deficiency pathway, initiated in response to Gram-negative bacte- the adaptive immune response, while IL-6 has diverse functions in both innate rial infections, leads to the cleavage of Relish and translocation of its REL do- and adaptive immunity (Kabelitz and Medzhitov 2007; Palm and Medzhitov main to the nucleus (Lemaitre 2004). The Drosophila genome encodes eight 2009) Toll-related receptors. Most of these have roles in development, but Toll-9 is capable of activating the drosomycin promoter, possibly through the Toll signal- ing pathway (Bilak et al. 2003). 2.1.2 Recognition strategies

The innate immune defense uses three strategies of recognition: recognition of microbial non-self, recognition of missing self, and recognition of altered self 2. The mammalian immune defense (Fig. 1, Medzhitov and Janeway Jr 2002). During evolution, newly developed defense mechanisms have been added to the old ones, resulting in a layered structure of the immune system (Zänker Recognition of microbial non-self is based on PRR’s – such as TLRs and 2008). The mammalian immune defense consists of two sub-systems – the in- MBL – interaction with MAMPs. MAMPs are conserved structures, invariant nate and the adaptive immune systems – which are interconnected and co- in a particular class of microorganisms, which are not present within the host. operate to provide host defense. The innate immune system has previously Furthermore, these structures are essential for the viability/adaptive fitness of been considered to be unsophisticated and non-specific. However, the growing the microorganism and thus not easily discarded (Medzhitov 2001). Secreted awareness of the immense importance of TLRs has promoted innate immunity PRRs, like MBL, bind to microbial cells and flag them for phagocytosis or eli- research and led to new insights concerning its elaborate nature and specificity. mination by the complement system, while membrane-bound PRRs, like TLRs, activate signaling pathways that induce antimicrobial effector mecha- nisms and inflammation (Medzhitov and Janeway Jr 2002). Since infections are 2.1 The innate immune defense likely to activate TLRs and complement simultaneously, it is reasonable to as- sume that these two danger-sensing systems intersect (Hawlish and Köhl 2006). Indeed, binding of the complement-derived anaphylatoxins C5a and 2.1.1 Epithelial barriers, effector cells, defense proteins, C3a to their receptors – C5aR and C3aR, respectively – modulates TLR4, and cytokines TLR2/6, and TLR9 signaling (Zhang et al. 2007). Moreover, PTX3, a PRR belonging to the pentraxin family, is produced by a variety of cells and tissues in The mammalian innate immune defense consists of epithelial barriers, ef- response to TLR signaling and modulates complement activation through in- fector cells, defense proteins – notably those of the complement cascade – and teraction with C1q and factor H (Deban et al. 2009; Garlanda et al. 2009). In- cytokines. The skin and the mucosal surfaces of the gastrointestinal and respi- tracellular crosstalk between complement and TLRs is also known. The com- ratory tracts prevent microbial entry into the body. They also produce antim- plement receptor-3 (CR3), an integrin, can be transactivated by TLR2 via an icrobial substances and harbor intraephitelial T lymphocytes. In addition, B-1 inside-out signaling pathway which is distinct from the MyD88-dependent cells are present in the peritoneal cavity. Intraephitelial T lymphocytes and B-1 pro-inflammatory signaling pathway. Conversely, CR3 can initiate TLR2 and cells express a more limited repertoire of antigen receptors than other T and B TLR4 by promoting the recruitment of the adapter molecule Mal (see below) cell populations and serve as sentinels at common sites of microbial invasion. (Hajishengallis and Lambris 2010). Cross-communication between TLRs and The main effector cell populations of the innate immune defense are granulo- complement may serve to avoid misinterpretation of signals from non- cytes (in particular neutrophils), macrophages, and natural killer (NK) cells. dangerous non-self and help in the fine-tuning of the subsequent immune re- These cell types are all known to express TLRs, a family of PRRs which are sponse to a particular microbe (Friec and Kemper 2009). homologous to the Drosophila Toll receptor. Neutrophils specialize in early phagocytosis while macrophages are more effective, since they persist longer at The missing self recognition strategy is based on molecular markers ex- sites of inflammation. NK cells lyse infected cells and activate macrophages. pressed on healthy cells: if these markers are missing, the cell is targeted for de- Cytokines are produced mainly by the effector cells, but also by endothelial and struction. The major histocompatibility complex (MHC) I proteins, constitu- some epithelial cells. Cytokines serve to communicate information between

18 19 way involving myeloid differentiation primary response protein 88 (MyD88) cells and have different areas of responsibility: TNF, IL-1, and chemokines and other proteins – to translocation of Dorsal-related immunity factor to the mediate inflammation, IFN-γ activates macrophages, while IL-15 and IL-12 nucleus, where it regulates genes coding for antimicrobial peptides. Similarly, stimulate proliferation and activity of NK cells. IL-12 is also a key inducer of the immune deficiency pathway, initiated in response to Gram-negative bacte- the adaptive immune response, while IL-6 has diverse functions in both innate rial infections, leads to the cleavage of Relish and translocation of its REL do- and adaptive immunity (Kabelitz and Medzhitov 2007; Palm and Medzhitov main to the nucleus (Lemaitre 2004). The Drosophila genome encodes eight 2009) Toll-related receptors. Most of these have roles in development, but Toll-9 is capable of activating the drosomycin promoter, possibly through the Toll signal- ing pathway (Bilak et al. 2003). 2.1.2 Recognition strategies

The innate immune defense uses three strategies of recognition: recognition of microbial non-self, recognition of missing self, and recognition of altered self 2. The mammalian immune defense (Fig. 1, Medzhitov and Janeway Jr 2002). During evolution, newly developed defense mechanisms have been added to the old ones, resulting in a layered structure of the immune system (Zänker Recognition of microbial non-self is based on PRR’s – such as TLRs and 2008). The mammalian immune defense consists of two sub-systems – the in- MBL – interaction with MAMPs. MAMPs are conserved structures, invariant nate and the adaptive immune systems – which are interconnected and co- in a particular class of microorganisms, which are not present within the host. operate to provide host defense. The innate immune system has previously Furthermore, these structures are essential for the viability/adaptive fitness of been considered to be unsophisticated and non-specific. However, the growing the microorganism and thus not easily discarded (Medzhitov 2001). Secreted awareness of the immense importance of TLRs has promoted innate immunity PRRs, like MBL, bind to microbial cells and flag them for phagocytosis or eli- research and led to new insights concerning its elaborate nature and specificity. mination by the complement system, while membrane-bound PRRs, like TLRs, activate signaling pathways that induce antimicrobial effector mecha- nisms and inflammation (Medzhitov and Janeway Jr 2002). Since infections are 2.1 The innate immune defense likely to activate TLRs and complement simultaneously, it is reasonable to as- sume that these two danger-sensing systems intersect (Hawlish and Köhl 2006). Indeed, binding of the complement-derived anaphylatoxins C5a and 2.1.1 Epithelial barriers, effector cells, defense proteins, C3a to their receptors – C5aR and C3aR, respectively – modulates TLR4, and cytokines TLR2/6, and TLR9 signaling (Zhang et al. 2007). Moreover, PTX3, a PRR belonging to the pentraxin family, is produced by a variety of cells and tissues in The mammalian innate immune defense consists of epithelial barriers, ef- response to TLR signaling and modulates complement activation through in- fector cells, defense proteins – notably those of the complement cascade – and teraction with C1q and factor H (Deban et al. 2009; Garlanda et al. 2009). In- cytokines. The skin and the mucosal surfaces of the gastrointestinal and respi- tracellular crosstalk between complement and TLRs is also known. The com- ratory tracts prevent microbial entry into the body. They also produce antim- plement receptor-3 (CR3), an integrin, can be transactivated by TLR2 via an icrobial substances and harbor intraephitelial T lymphocytes. In addition, B-1 inside-out signaling pathway which is distinct from the MyD88-dependent cells are present in the peritoneal cavity. Intraephitelial T lymphocytes and B-1 pro-inflammatory signaling pathway. Conversely, CR3 can initiate TLR2 and cells express a more limited repertoire of antigen receptors than other T and B TLR4 by promoting the recruitment of the adapter molecule Mal (see below) cell populations and serve as sentinels at common sites of microbial invasion. (Hajishengallis and Lambris 2010). Cross-communication between TLRs and The main effector cell populations of the innate immune defense are granulo- complement may serve to avoid misinterpretation of signals from non- cytes (in particular neutrophils), macrophages, and natural killer (NK) cells. dangerous non-self and help in the fine-tuning of the subsequent immune re- These cell types are all known to express TLRs, a family of PRRs which are sponse to a particular microbe (Friec and Kemper 2009). homologous to the Drosophila Toll receptor. Neutrophils specialize in early phagocytosis while macrophages are more effective, since they persist longer at The missing self recognition strategy is based on molecular markers ex- sites of inflammation. NK cells lyse infected cells and activate macrophages. pressed on healthy cells: if these markers are missing, the cell is targeted for de- Cytokines are produced mainly by the effector cells, but also by endothelial and struction. The major histocompatibility complex (MHC) I proteins, constitu- some epithelial cells. Cytokines serve to communicate information between

18 19 tively expressed on all nucleated cells but often down-regulated as a result of vi- Kemper 2009). Moreover, it has recently been shown that properdin can act as ral infection, serves as a self marker and ligand for inhibitory receptors which a recognition molecule and initiate the alternative pathway (Kemper et al. block the lytic activity of NK cells. Conversely, recognition of altered self is 2009). The alternative pathway has previously been considered to be the origi- based on markers expressed only by abnormal or damaged cells, which thus are nal route of activation, since it does not require antibodies, a feature of the flagged for elimination (Medzhitov and Janeway Jr 2002). adaptive immune system, for its function. However, there is accumulating evi- dence that complement originated as a lectin-based opsonisation system (Fujita et al. 2004). The activating pathways of the complement system lead to the generation of the two C3 convertases, which, in turn, generate C3b, the major mammalian opsonin, and the anaphylatoxin C3a. Further activation initiates the formation of the C5 convertases. These generate the potent anaphyla- toxin/chemotaxin C5a as well as C5b, which initiates the formation of the C5b-9 complex, which, in turn, promotes cell lysis (Friec and Kemper 2009).

The lectin pathway is initiated by the binding of MBLs or to man- nose and N-acetyl glucose amine on the surface of pathogens (Fig. 2). Ficolin-3 has been found to be the most potent of the lectin pathway initiators in hu- mans, followed by ficolin-2 and MBL (Hummelshoj et al. 2008). MBL is similar in structure to C1q, a subunit of the classical pathway C1 complex, and is thought to have the same ability as C1q to stimulate phagocytosis through the C1qRp receptor (Friec and Kemper 2009; Homskov et al. 2003; Phatsara et al. 2007). MBL and ficolins interact with MASPs, the function of which is equivalent to that of the C1 complex: C4 and C2 are cleaved and the classi- cal/lectin pathway C3 and C5 convertases – C4b2a and C4b2aC3b, respec- tively – are formed (Friec and Kemper 2009; Lillie et al. 2005). Recently, a novel MBL/ficolin associated protein, denoted MAP-1, has been detected in humans. This protein is expressed in myocardial and skeletal muscle and Fig. 1 Three recognition strategies used by the innate immune defense. A. Recognition of microbial thought to inhibit the complement system by preventing the cleavage of C4 non-self induces immune response. B. NK cells interact with target cells through activating and in- (Skjoedt et al. 2010). Furthermore, there are indications that MASP-1 might hibitory receptors. When both types of receptors are engaged, the inhibitory receptors are dominant be the initiator of the alternative complement activation pathway (Takahashi et and the NK cell is not activated. However, if self marker molecules are missing, the NK cell is re- al. 2010). leased from its state of inhibition. C. Expression of markers of altered self flags the cell for destruc- tion.

2.1.3 The complement system

The mammalian complement system consists of more than 30 membrane- bound and soluble proteins and traces its origin back at least to Cnidaria: two C3, two factor B, and one mannan-binding protein associated serine protease (MASP) is present in Nematostella vectensis, the sea anemone (Kimura et al. 2009). The mammalian complement can be activated through three major pathways – the classical, lectin, and alternative pathways – and direct cleavage of C3 or C5 by proteases from the clotting cascade is also known (Friec and

20 21 tively expressed on all nucleated cells but often down-regulated as a result of vi- Kemper 2009). Moreover, it has recently been shown that properdin can act as ral infection, serves as a self marker and ligand for inhibitory receptors which a recognition molecule and initiate the alternative pathway (Kemper et al. block the lytic activity of NK cells. Conversely, recognition of altered self is 2009). The alternative pathway has previously been considered to be the origi- based on markers expressed only by abnormal or damaged cells, which thus are nal route of activation, since it does not require antibodies, a feature of the flagged for elimination (Medzhitov and Janeway Jr 2002). adaptive immune system, for its function. However, there is accumulating evi- dence that complement originated as a lectin-based opsonisation system (Fujita et al. 2004). The activating pathways of the complement system lead to the generation of the two C3 convertases, which, in turn, generate C3b, the major mammalian opsonin, and the anaphylatoxin C3a. Further activation initiates the formation of the C5 convertases. These generate the potent anaphyla- toxin/chemotaxin C5a as well as C5b, which initiates the formation of the C5b-9 complex, which, in turn, promotes cell lysis (Friec and Kemper 2009).

The lectin pathway is initiated by the binding of MBLs or ficolins to man- nose and N-acetyl glucose amine on the surface of pathogens (Fig. 2). Ficolin-3 has been found to be the most potent of the lectin pathway initiators in hu- mans, followed by ficolin-2 and MBL (Hummelshoj et al. 2008). MBL is similar in structure to C1q, a subunit of the classical pathway C1 complex, and is thought to have the same ability as C1q to stimulate phagocytosis through the C1qRp receptor (Friec and Kemper 2009; Homskov et al. 2003; Phatsara et al. 2007). MBL and ficolins interact with MASPs, the function of which is equivalent to that of the C1 complex: C4 and C2 are cleaved and the classi- cal/lectin pathway C3 and C5 convertases – C4b2a and C4b2aC3b, respec- tively – are formed (Friec and Kemper 2009; Lillie et al. 2005). Recently, a novel MBL/ficolin associated protein, denoted MAP-1, has been detected in humans. This protein is expressed in myocardial and skeletal muscle and Fig. 1 Three recognition strategies used by the innate immune defense. A. Recognition of microbial thought to inhibit the complement system by preventing the cleavage of C4 non-self induces immune response. B. NK cells interact with target cells through activating and in- (Skjoedt et al. 2010). Furthermore, there are indications that MASP-1 might hibitory receptors. When both types of receptors are engaged, the inhibitory receptors are dominant be the initiator of the alternative complement activation pathway (Takahashi et and the NK cell is not activated. However, if self marker molecules are missing, the NK cell is re- al. 2010). leased from its state of inhibition. C. Expression of markers of altered self flags the cell for destruc- tion.

2.1.3 The complement system

The mammalian complement system consists of more than 30 membrane- bound and soluble proteins and traces its origin back at least to Cnidaria: two C3, two factor B, and one mannan-binding protein associated serine protease (MASP) is present in Nematostella vectensis, the sea anemone (Kimura et al. 2009). The mammalian complement can be activated through three major pathways – the classical, lectin, and alternative pathways – and direct cleavage of C3 or C5 by proteases from the clotting cascade is also known (Friec and

20 21 as a sensor to maintain homeostasis (Portnoy 2005). In such a scenario, avoid- ing and initiating inflammation are of equal importance. Commensal bacteria play an important role in maintaining tolerance and active stability of the intes- tinal epithelial barrier by suppressing inflammation (Cario et al. 2007). Intesti- nal epithelial cells constitutively express TLRs, but it seems that this expression is restricted to specific cell lineages and, in some cases, to the basolateral surface of the cell. Thus, pathogenic bacteria, which can penetrate the epithelial bar- rier, are recognized as a threat and elicit a pro-inflammatory response, while commensal bacteria, which remain on the apical side, elicit a homeostatic anti- inflammatory response (Abreau 2010). In polarized human HCA-7 and Caco- 2 cells (cancer cell lines), TLR9 is expressed on both the apical and basolateral surfaces, but only basolaterally localized TLR9 induce a pro-inflammatory re- sponse (Lee et al. 2008). Moreover, it has been shown that TLR2 signaling protects tight junctions against stress-induced damage in vivo (Cario et al. 2007). Fig. 2 Activation of complement by MBL. MBL (and ficolins) recognizes carbohydrate residues on pathogen cell surfaces. Through interaction with MASPs, C4 and C2 are cleaved and the C3 con- vertase is formed. sMAP is also part of the MBL/MASP complex. C3a and C4a are anaphylatox- ins, while the function of C2b is unknown. Subsequent to the primary cleavage of C3 by the C3 con- 2.2 Main features of the adaptive immune defense vertase, factor I (together with co-factors) cleaves C3b in two steps, thus producing iC3b and C3d. The mammalian adaptive immune system traces its origin to the jawed ver- The downstream effector functions of complement are opsonisation, cell lysis, and inflammation. tebrates about 500 million years ago. The emergence of this more specific, di- verse, and specialized defense system is conventionally linked to the emergence of the recombination activating genes (RAG), which encode recombinase pro- teins essential for somatic recombination (Friec and Kemper 2009). However, 2.1.4 Resistance against innate immunity and mainte- albeit with unknown function, there is a gene cluster with similarity to the ja- wed vertebrate RAG 1 and 2 genes present in the sea urchin genome (Hibino nance of homeostasis et al. 2006), and the sea lamprey possesses variable lymphocyte receptors, with Many pathogens have developed strategies to avoid or exploit the innate somatically rearranged extracellular domains (Pancer et al. 2004). Thus, the ad- immune response to promote their pathogenesis. The pathogenicity of Yersinia aptive immune system seen in today’s mammals may not be the only existing pestis relies on several mechanisms. It is known that this bacterium is capable of form of adaptive immunity (Huang et al. 2008). producing an altered lipid A structure, which does not fully stimulate TLR4, when grown under conditions resembling those in the human host. Moreover, The mammalian adaptive immune defense has a cell-mediated and a hu- it can also stimulate production of IL-10, which is the main inhibitor of acti- moral branch, based on the properties of B and T cells and on antibodies, re- vated macrophages and dendritic cells (DCs), in a TLR2-dependent manner spectively. Somatic recombination accounts, to a large extent, for the diversity and cause depletion of NK cells. Similarly, Mycobacterium tuberculosis is capable of the adaptive immune defense. In this process, the sequences in the immu- of down-regulating IL-12 expression as well as inhibiting macrophage re- noglobulin (Ig) and T cell receptor (TCR) loci in immature B and T cells are γ sponses to IFN- in a TLR2-dependent manner (Portnoy 2005). Moreover, recombined through enzymatic activity, to form functional antigen receptors. Candida albicans, a pathogenic yeast, has got at least two surface proteins which Unlike the innate immune system, the adaptive immune system has got a bind the complement regulators factor H and factor H like protein-1; thus, memory, which enables a quick response to a specific pathogen on the second Candida disguises itself and avoids attention from the complement system (Luo encounter. This property of the adaptive immune defense is the basis for vacci- et al. 2009). nation. Moreover, self-limitation and non-reactivity to self ensures that host injury is prevented (Abbas and Janeway 2000). Since many microorganisms are able to avoid and/or exploit the mammalian innate immune system, it has been suggested that one of its main functions is

22 23 as a sensor to maintain homeostasis (Portnoy 2005). In such a scenario, avoid- ing and initiating inflammation are of equal importance. Commensal bacteria play an important role in maintaining tolerance and active stability of the intes- tinal epithelial barrier by suppressing inflammation (Cario et al. 2007). Intesti- nal epithelial cells constitutively express TLRs, but it seems that this expression is restricted to specific cell lineages and, in some cases, to the basolateral surface of the cell. Thus, pathogenic bacteria, which can penetrate the epithelial bar- rier, are recognized as a threat and elicit a pro-inflammatory response, while commensal bacteria, which remain on the apical side, elicit a homeostatic anti- inflammatory response (Abreau 2010). In polarized human HCA-7 and Caco- 2 cells (cancer cell lines), TLR9 is expressed on both the apical and basolateral surfaces, but only basolaterally localized TLR9 induce a pro-inflammatory re- sponse (Lee et al. 2008). Moreover, it has been shown that TLR2 signaling protects tight junctions against stress-induced damage in vivo (Cario et al. 2007). Fig. 2 Activation of complement by MBL. MBL (and ficolins) recognizes carbohydrate residues on pathogen cell surfaces. Through interaction with MASPs, C4 and C2 are cleaved and the C3 con- vertase is formed. sMAP is also part of the MBL/MASP complex. C3a and C4a are anaphylatox- ins, while the function of C2b is unknown. Subsequent to the primary cleavage of C3 by the C3 con- 2.2 Main features of the adaptive immune defense vertase, factor I (together with co-factors) cleaves C3b in two steps, thus producing iC3b and C3d. The mammalian adaptive immune system traces its origin to the jawed ver- The downstream effector functions of complement are opsonisation, cell lysis, and inflammation. tebrates about 500 million years ago. The emergence of this more specific, di- verse, and specialized defense system is conventionally linked to the emergence of the recombination activating genes (RAG), which encode recombinase pro- teins essential for somatic recombination (Friec and Kemper 2009). However, 2.1.4 Resistance against innate immunity and mainte- albeit with unknown function, there is a gene cluster with similarity to the ja- wed vertebrate RAG 1 and 2 genes present in the sea urchin genome (Hibino nance of homeostasis et al. 2006), and the sea lamprey possesses variable lymphocyte receptors, with Many pathogens have developed strategies to avoid or exploit the innate somatically rearranged extracellular domains (Pancer et al. 2004). Thus, the ad- immune response to promote their pathogenesis. The pathogenicity of Yersinia aptive immune system seen in today’s mammals may not be the only existing pestis relies on several mechanisms. It is known that this bacterium is capable of form of adaptive immunity (Huang et al. 2008). producing an altered lipid A structure, which does not fully stimulate TLR4, when grown under conditions resembling those in the human host. Moreover, The mammalian adaptive immune defense has a cell-mediated and a hu- it can also stimulate production of IL-10, which is the main inhibitor of acti- moral branch, based on the properties of B and T cells and on antibodies, re- vated macrophages and dendritic cells (DCs), in a TLR2-dependent manner spectively. Somatic recombination accounts, to a large extent, for the diversity and cause depletion of NK cells. Similarly, Mycobacterium tuberculosis is capable of the adaptive immune defense. In this process, the sequences in the immu- of down-regulating IL-12 expression as well as inhibiting macrophage re- noglobulin (Ig) and T cell receptor (TCR) loci in immature B and T cells are γ sponses to IFN- in a TLR2-dependent manner (Portnoy 2005). Moreover, recombined through enzymatic activity, to form functional antigen receptors. Candida albicans, a pathogenic yeast, has got at least two surface proteins which Unlike the innate immune system, the adaptive immune system has got a bind the complement regulators factor H and factor H like protein-1; thus, memory, which enables a quick response to a specific pathogen on the second Candida disguises itself and avoids attention from the complement system (Luo encounter. This property of the adaptive immune defense is the basis for vacci- et al. 2009). nation. Moreover, self-limitation and non-reactivity to self ensures that host injury is prevented (Abbas and Janeway 2000). Since many microorganisms are able to avoid and/or exploit the mammalian innate immune system, it has been suggested that one of its main functions is

22 23 2.3 Bridges between the innate and adaptive immune systems The innate and adaptive immune systems are interconnected mainly thro- ugh antigen presentation by DCs and the complement system (Fig. 3).

Fig. 4 Antigen presentation. Upon ligand binding/uptake and maturation, DCs present antigen to Fig. 3 Schematic representation of the mammalian immune defense. The innate and adaptive im- antigen-specific naïve T cells. Three signals from the DC – antigen presentation (signal 1), expres- mune systems are interconnected mainly through DCs (and IL-12, produced by DCs) and the com- sion of co-stimulatory molecules (signal 2), and secretion of pro-inflammatory cytokines (signal 3) – plement system. The innate immune system is always active, but has limited specificity and no are necessary for T cell activation. memory. The adaptive immune system needs to reach full capacity, but is highly specific and has got a memory. It is believed that direct interaction with pathogens may not be a prerequi- site for DC maturation; signaling via receptors for pro-inflammatory cytokines Antigen presentation is accomplished through uptake and/or processing of released by leukocytes, including other DCs, seems to be enough to elicit the antigen in an antigen presenting cell (APC), followed by exposure of the anti- mature DC phenotype. However, DC maturation is not equivalent to DC im- gen-derived product on the cell surface together with MHC I and II proteins munogenicity, and it seems that DCs which have matured solely under the in- and co-stimulatory molecules, notably CD80 and CD86 (Fig. 4). DCs are the fluence of cytokines are unable to bring about T cell activation, since these DCs most specialized antigen presenting cells. Immature DCs express TLRs and fail to produce pro-inflammatory cytokines. Interestingly, DCs residing in the other PRRs and are present under epithelia and in most organs. When acti- gut, when conditioned by neighboring intestinal epithelial cells, can contribute vated, they mature and migrate to the lymph nodes, where they present antigen to tolerance to commensal microorganisms by releasing IL-10 and IL-6, but to antigen-specific naïve T cells. Pathogen-activated DCs are also a main not IL-12. Similar mechanisms may also apply in the lung (Joffre et al. 2009). source of IL-12, which, in turn, stimulates the differentiation of CD4+ T cells and enhances the cytolytic function of CD8+ T cells (Joffre et al. 2009). The complement system modulates both B and T cell responses. In B cells, this is achieved through CD21 and C3aR. CD21 is part of a co-receptor com- plex expressed on B cells which, when binding C3d deposited on a pathogenic surface, facilitates the activation of the cell in co-operation with the B cell re- ceptor (BCR). Through this co-operation, the threshold value of pathogens needed to achive B cell activation is lowered, since it replaces the otherwise necessary crossbinding of two BCRs (Mongini et al. 1997). On the other hand, B cells also express C3aR, which, when activated, suppresses the polyclonal

24 25 2.3 Bridges between the innate and adaptive immune systems The innate and adaptive immune systems are interconnected mainly thro- ugh antigen presentation by DCs and the complement system (Fig. 3).

Fig. 4 Antigen presentation. Upon ligand binding/uptake and maturation, DCs present antigen to Fig. 3 Schematic representation of the mammalian immune defense. The innate and adaptive im- antigen-specific naïve T cells. Three signals from the DC – antigen presentation (signal 1), expres- mune systems are interconnected mainly through DCs (and IL-12, produced by DCs) and the com- sion of co-stimulatory molecules (signal 2), and secretion of pro-inflammatory cytokines (signal 3) – plement system. The innate immune system is always active, but has limited specificity and no are necessary for T cell activation. memory. The adaptive immune system needs to reach full capacity, but is highly specific and has got a memory. It is believed that direct interaction with pathogens may not be a prerequi- site for DC maturation; signaling via receptors for pro-inflammatory cytokines Antigen presentation is accomplished through uptake and/or processing of released by leukocytes, including other DCs, seems to be enough to elicit the antigen in an antigen presenting cell (APC), followed by exposure of the anti- mature DC phenotype. However, DC maturation is not equivalent to DC im- gen-derived product on the cell surface together with MHC I and II proteins munogenicity, and it seems that DCs which have matured solely under the in- and co-stimulatory molecules, notably CD80 and CD86 (Fig. 4). DCs are the fluence of cytokines are unable to bring about T cell activation, since these DCs most specialized antigen presenting cells. Immature DCs express TLRs and fail to produce pro-inflammatory cytokines. Interestingly, DCs residing in the other PRRs and are present under epithelia and in most organs. When acti- gut, when conditioned by neighboring intestinal epithelial cells, can contribute vated, they mature and migrate to the lymph nodes, where they present antigen to tolerance to commensal microorganisms by releasing IL-10 and IL-6, but to antigen-specific naïve T cells. Pathogen-activated DCs are also a main not IL-12. Similar mechanisms may also apply in the lung (Joffre et al. 2009). source of IL-12, which, in turn, stimulates the differentiation of CD4+ T cells and enhances the cytolytic function of CD8+ T cells (Joffre et al. 2009). The complement system modulates both B and T cell responses. In B cells, this is achieved through CD21 and C3aR. CD21 is part of a co-receptor com- plex expressed on B cells which, when binding C3d deposited on a pathogenic surface, facilitates the activation of the cell in co-operation with the B cell re- ceptor (BCR). Through this co-operation, the threshold value of pathogens needed to achive B cell activation is lowered, since it replaces the otherwise necessary crossbinding of two BCRs (Mongini et al. 1997). On the other hand, B cells also express C3aR, which, when activated, suppresses the polyclonal

24 25 immune response and limits the production of IL-6 and TNF-α. Complement modulation of T cell responses is achieved through receptors on the T cell itself and also via APCs, which express a wide repertoire of complement receptors. PART II: PATTERN RECOG- In the absence of C3, antigen uptake by APCs is suboptimal, which leads to a suboptimal APC maturation and thus suboptimal T cell activation. Comple- NITION RECEPTORS ment also modulates IL-12 production by APCs (Friec and Kemper 2009).

The innate immune defense uses several classes of PRRs for detection of microbial presence and altered self. Some are present in the circulation, while others are found in the cytoplasm or on the cell surface.

3. Selected pattern recognition receptors and their func- tions

C reactive protein (CRP), the first PRR identified, as well as serum amyloid P component (SAP) and pentraxin 3 (PTX3) are secreted PRRs belonging to the protein family of pentraxins (Deban et al. 2009). CRP and SAP are opson- ins and activators of the classical complement pathway (Manfredi et al. 2008; Medzhitov 2001). PTX3 binds to and is non-redundant in its protective func- tion against certain pathogens, notably influenza viruses. Furthermore, it is re- quired for female fertility and behaves as an inhibitor of the fibroblast growth factor-2 (Manfredi et al. 2009).

Protein kinase R (PKR) and nucleotide-binding oligomerization domain- containing protein (NOD) 1 and 2 are intracellular PRRs. PKR interacts with viral RNA present in the cytoplasm and phosphorylates eukaryotic initiation α factor 2 ; this way, translation initiation is inhibited and pathogen replication is prevented (Ranjan et al. 2009). NOD1 and 2 both recognize derivatives from peptidoglycan; subsequent to ligand binding, they undergo conformation κ κ changes and self-oligomerization and induce nuclear factor B (NF B) activa- tion independent of TLRs (Shaw et al. 2008).

Macrophage and macrophage scavenger receptor are cell surface PRRs, while CD14 exists in a cell-bound and a soluble version. Macro- phage mannose receptor and macrophage scavenger receptor function in the phagocytosis of macrophages (Medzhitov 2001). CD14 binds LPS and chap- erones its transfer to the TLR4 receptor complex. Also, signaling through the TLR2/TLR6 heterodimer is partly dependent on CD14: elimination of CD14 impairs sensing of all TLR2/TLR6 ligands. It has been shown that CD14 in- teracts directly with both TLR2 and TLR4. Moreover, CD14 binds double- stranded RNA and interacts physically with TLR3 within intracellular com- partments, thus aiding in NFκB activation (Akashi-Takamura and Miyake K 2008; Lee et al. 2006).

26 27 immune response and limits the production of IL-6 and TNF-α. Complement modulation of T cell responses is achieved through receptors on the T cell itself and also via APCs, which express a wide repertoire of complement receptors. PART II: PATTERN RECOG- In the absence of C3, antigen uptake by APCs is suboptimal, which leads to a suboptimal APC maturation and thus suboptimal T cell activation. Comple- NITION RECEPTORS ment also modulates IL-12 production by APCs (Friec and Kemper 2009).

The innate immune defense uses several classes of PRRs for detection of microbial presence and altered self. Some are present in the circulation, while others are found in the cytoplasm or on the cell surface.

3. Selected pattern recognition receptors and their func- tions

C reactive protein (CRP), the first PRR identified, as well as serum amyloid P component (SAP) and pentraxin 3 (PTX3) are secreted PRRs belonging to the protein family of pentraxins (Deban et al. 2009). CRP and SAP are opson- ins and activators of the classical complement pathway (Manfredi et al. 2008; Medzhitov 2001). PTX3 binds to and is non-redundant in its protective func- tion against certain pathogens, notably influenza viruses. Furthermore, it is re- quired for female fertility and behaves as an inhibitor of the fibroblast growth factor-2 (Manfredi et al. 2009).

Protein kinase R (PKR) and nucleotide-binding oligomerization domain- containing protein (NOD) 1 and 2 are intracellular PRRs. PKR interacts with viral RNA present in the cytoplasm and phosphorylates eukaryotic initiation α factor 2 ; this way, translation initiation is inhibited and pathogen replication is prevented (Ranjan et al. 2009). NOD1 and 2 both recognize derivatives from peptidoglycan; subsequent to ligand binding, they undergo conformation κ κ changes and self-oligomerization and induce nuclear factor B (NF B) activa- tion independent of TLRs (Shaw et al. 2008).

Macrophage mannose receptor and macrophage scavenger receptor are cell surface PRRs, while CD14 exists in a cell-bound and a soluble version. Macro- phage mannose receptor and macrophage scavenger receptor function in the phagocytosis of macrophages (Medzhitov 2001). CD14 binds LPS and chap- erones its transfer to the TLR4 receptor complex. Also, signaling through the TLR2/TLR6 heterodimer is partly dependent on CD14: elimination of CD14 impairs sensing of all TLR2/TLR6 ligands. It has been shown that CD14 in- teracts directly with both TLR2 and TLR4. Moreover, CD14 binds double- stranded RNA and interacts physically with TLR3 within intracellular com- partments, thus aiding in NFκB activation (Akashi-Takamura and Miyake K 2008; Lee et al. 2006).

26 27 4. The Toll-like receptors The first TLR was cloned and characterized in 1997 by Ruslan Medzhitov and co-workers: under the hypothesis that the expression of co-stimulatory molecules and cytokines by APCs was induced by a non-clonal component of immunity, human TLR4 was discovered. It was also found that this molecule could induce activation of the NFκB pathway and expression of the co- stimulatory molecule B7.1, which is necessary for activation of naïve T cells (Medzhitov et al. 1997). TLRs have previously been considered solely dedi- cated to host defense (Imler and Zheng 2004). However, evidence is emerging that these receptors may also be involved in central nervous system develop- ment and maintenance (Larsen et al. 2007; Ma et al. 2006).

Fig. 5 Evolution of TLRs. TLR4 is the eldest TLR gene, dating its origin before vertebrate life. 4.1 Evolutionary perspective TLR10 arose about 300 million years ago while the TLR1/TLR6 split occurred about 130 million TLR4 is the most ancient TLR, dating its origin before vertebrate life (Fig. years ago. 5, Beutler and Rehli 2002). Today’s mammalian TLR family can be subdivided into two groups, which are the result of a gene duplication event prior to the In mice and humans, remnants of a second, disrupted TLR2-like gene is lo- divergence of invertebrates and vertebrates: the TLR1 family, consisting of cated in tandem with the functional version of TLR2. This indicates the occur- TLRs 1, 2, 6, and 10, and a second group, comprised of all other TLRs. The rence of a gene duplication after which one of the gene copies has developed TLR1, TLR6, and TLR10 genes are present only in mammals, suggesting du- into a pseudogene. Since there are two TLR2 genes in chicken, the duplication plication within the mammalian lineage (Hughes and Piontkivska 2008). must have occurred before the divergence of mammals and birds (Boyd et al. TLR10 separated from the TLR1/6 precursor about 300 million years ago, 2007). roughly at the time of divergence from birds, while TLR1 and TLR6 appeared approximately 170 million years later (Beutler and Rehli 2002). Following du- The nematode Caenorhabditis elegans possesses a single TLR (Fig. 6, Tenor plication, gene conversion events between TLR1 and TLR6, limiting the diver- and Aballay 2008). On the other hand, the genome of the purple sea urchin gence between these two genes, have been suggested (Kruithof et al. 2007). encodes 222 TLRs (Hibino et al. 2006), indicating that during evolution, ex- pansion of the TLR repertoire has provided a survival advantage in this organ- ism. Further, the genome of the Florida lancelet is thought to harbor roughly 50 TLRs (Huang et al. 2008), while at least 20 are known in the South African clawed frog (Ishii et al. 2007). In the chicken genome, 10 TLR genes have been identified, out of which 2 – chTLR15 and chTLR21 – are specific for chicken (Boyd et al. 2007). Ten functional TLRs are known in humans and pigs (Uenishi and Shinkai 2008). In mice, TLR11, TLR12, and TLR13 have been identified (Roach et al. 2005), but mouse TLR10 is nonfunctional, due to a retroviral insertion (Hasan et al. 2005). TLR11 is present as a pseudogene in humans (Takeda and Akira 2005).

28 29 4. The Toll-like receptors The first TLR was cloned and characterized in 1997 by Ruslan Medzhitov and co-workers: under the hypothesis that the expression of co-stimulatory molecules and cytokines by APCs was induced by a non-clonal component of immunity, human TLR4 was discovered. It was also found that this molecule could induce activation of the NFκB pathway and expression of the co- stimulatory molecule B7.1, which is necessary for activation of naïve T cells (Medzhitov et al. 1997). TLRs have previously been considered solely dedi- cated to host defense (Imler and Zheng 2004). However, evidence is emerging that these receptors may also be involved in central nervous system develop- ment and maintenance (Larsen et al. 2007; Ma et al. 2006).

Fig. 5 Evolution of TLRs. TLR4 is the eldest TLR gene, dating its origin before vertebrate life. 4.1 Evolutionary perspective TLR10 arose about 300 million years ago while the TLR1/TLR6 split occurred about 130 million TLR4 is the most ancient TLR, dating its origin before vertebrate life (Fig. years ago. 5, Beutler and Rehli 2002). Today’s mammalian TLR family can be subdivided into two groups, which are the result of a gene duplication event prior to the In mice and humans, remnants of a second, disrupted TLR2-like gene is lo- divergence of invertebrates and vertebrates: the TLR1 family, consisting of cated in tandem with the functional version of TLR2. This indicates the occur- TLRs 1, 2, 6, and 10, and a second group, comprised of all other TLRs. The rence of a gene duplication after which one of the gene copies has developed TLR1, TLR6, and TLR10 genes are present only in mammals, suggesting du- into a pseudogene. Since there are two TLR2 genes in chicken, the duplication plication within the mammalian lineage (Hughes and Piontkivska 2008). must have occurred before the divergence of mammals and birds (Boyd et al. TLR10 separated from the TLR1/6 precursor about 300 million years ago, 2007). roughly at the time of divergence from birds, while TLR1 and TLR6 appeared approximately 170 million years later (Beutler and Rehli 2002). Following du- The nematode Caenorhabditis elegans possesses a single TLR (Fig. 6, Tenor plication, gene conversion events between TLR1 and TLR6, limiting the diver- and Aballay 2008). On the other hand, the genome of the purple sea urchin gence between these two genes, have been suggested (Kruithof et al. 2007). encodes 222 TLRs (Hibino et al. 2006), indicating that during evolution, ex- pansion of the TLR repertoire has provided a survival advantage in this organ- ism. Further, the genome of the Florida lancelet is thought to harbor roughly 50 TLRs (Huang et al. 2008), while at least 20 are known in the South African clawed frog (Ishii et al. 2007). In the chicken genome, 10 TLR genes have been identified, out of which 2 – chTLR15 and chTLR21 – are specific for chicken (Boyd et al. 2007). Ten functional TLRs are known in humans and pigs (Uenishi and Shinkai 2008). In mice, TLR11, TLR12, and TLR13 have been identified (Roach et al. 2005), but mouse TLR10 is nonfunctional, due to a retroviral insertion (Hasan et al. 2005). TLR11 is present as a pseudogene in humans (Takeda and Akira 2005).

28 29 global structure, with a conserved patch believed to be important for the inter- action with the down-stream adapter molecule MyD88 (Xu et al. 2000).

TLRs form homodimers and heterodimers, which, in turn, are part of re- ceptor complexes including non-TLR co-receptors (Table 1, Fig. 7). TLRs 1, 2, 4, 5, 6, and 10 are expressed on the cell surface and interact with bacterial ligands as well as with certain MAMPs derived from fungi and parasites, while TLRs 3, 7, 8, and 9 are located in intracellular compartments and recognize vi- ral nucleic acids (Akira et al. 2006). For TLR10, no ligand has been identified (Hasan et al. 2005), but evolutionary analysis indicates that a lipopeptide is probable (Hughes and Piontkivska 2008). Concerning human TLRs, it is be- lieved that those expressed on the cell surface are more redundant than the nu- cleic acid sensors in intracellular compartments (Barreiro et al. 2009).

Table 1 Heterodimeric partner(s), expression site, main ligand(s), and adapter molecules of TLRs 1-10.

a TLR Hd Expression Main ligand (s) Adapter molecule(s)

1 2, 10 cell membrane lipopeptides MyD88 Fig. 6 Schematic and simplified overview of the development of the immune system. Numbers of c 2 1, 6, 10 cell membrane lipopeptides, PG MyD88, Mal TLRs vary between species and tend to decrease during the course of evolution. Proteins belonging to b d the complement system are present in the purple sea urchin. The South African clawed frog possesses 3 int comp dsRNA TRIF e an adaptive immune system based on B and T lymphocytes. 4 cell membrane LPS MyD88, Mal, TRIF, TRAM 5 cell membrane flagelline MyD88 6 2 cell membrane lipopeptides MyD88 b f 4.2 The TLR molecule 7 int comp ssRNA MyD88 2 b f 8 inh int comp ssRNA MyD88 TLRs are type one transmembrane glycoproteins. Their extracellular do- 3 b,1 g mains are composed of 19-25 leucine-rich repeats (LRRs), often flanked by 9 inh int comp CpG DNA MyD88 cystein-rich motifs which represent capping structures at both ends of the su- 10 1, 2 cell membrane MyD88 perhelix formed by the LRRs. The LRRs are 24-29 amino acids in length and a b c contain the motif XLXXLXLXX. The intracellular TLR domain holds a sig- heterodimeric partner(s); intracellular compartments; peptidoglycan from Gram+ d e f naling entity which is homologous to that of the IL-1 receptor and thus termed bacteria; double-stranded RNA; lipopolysaccharide from Gram- bacteria; single- g Toll/IL-1R (TIR) homology domain (Akira et al. 2006; Gay and Gangloff stranded RNA; unmethylated DNA from bacteria 1 2007). In follicular-associated epithelial cells in mouse intestine, TLR9 is expressed on the apical and basolateral surfaces of the cell (Abreau 2010). 2 The crystal structure determined for the TLR1-TLR2-lipopeptide complex Inhibitory interactions with TLR7 and TLR9 are known (Wang et al. 2006). 3 reveals that the extracellular TLR domains adopt horse-shoe-like structures Inhibitory interactions with TLR7 are known (Wang et al. 2006). and that the ligand binding site of this TLR heterodimer is located on the con- vex surface of the horse-shoes (Jin et al. 2007). The latter, however, is not nec- essarily a general rule (Akira et al. 2006). The TIR domain is a more or less

30 31 global structure, with a conserved patch believed to be important for the inter- action with the down-stream adapter molecule MyD88 (Xu et al. 2000).

TLRs form homodimers and heterodimers, which, in turn, are part of re- ceptor complexes including non-TLR co-receptors (Table 1, Fig. 7). TLRs 1, 2, 4, 5, 6, and 10 are expressed on the cell surface and interact with bacterial ligands as well as with certain MAMPs derived from fungi and parasites, while TLRs 3, 7, 8, and 9 are located in intracellular compartments and recognize vi- ral nucleic acids (Akira et al. 2006). For TLR10, no ligand has been identified (Hasan et al. 2005), but evolutionary analysis indicates that a lipopeptide is probable (Hughes and Piontkivska 2008). Concerning human TLRs, it is be- lieved that those expressed on the cell surface are more redundant than the nu- cleic acid sensors in intracellular compartments (Barreiro et al. 2009).

Table 1 Heterodimeric partner(s), expression site, main ligand(s), and adapter molecules of TLRs 1-10.

a TLR Hd Expression Main ligand (s) Adapter molecule(s)

1 2, 10 cell membrane lipopeptides MyD88 Fig. 6 Schematic and simplified overview of the development of the immune system. Numbers of c 2 1, 6, 10 cell membrane lipopeptides, PG MyD88, Mal TLRs vary between species and tend to decrease during the course of evolution. Proteins belonging to b d the complement system are present in the purple sea urchin. The South African clawed frog possesses 3 int comp dsRNA TRIF e an adaptive immune system based on B and T lymphocytes. 4 cell membrane LPS MyD88, Mal, TRIF, TRAM 5 cell membrane flagelline MyD88 6 2 cell membrane lipopeptides MyD88 b f 4.2 The TLR molecule 7 int comp ssRNA MyD88 2 b f 8 inh int comp ssRNA MyD88 TLRs are type one transmembrane glycoproteins. Their extracellular do- 3 b,1 g mains are composed of 19-25 leucine-rich repeats (LRRs), often flanked by 9 inh int comp CpG DNA MyD88 cystein-rich motifs which represent capping structures at both ends of the su- 10 1, 2 cell membrane MyD88 perhelix formed by the LRRs. The LRRs are 24-29 amino acids in length and a b c contain the motif XLXXLXLXX. The intracellular TLR domain holds a sig- heterodimeric partner(s); intracellular compartments; peptidoglycan from Gram+ d e f naling entity which is homologous to that of the IL-1 receptor and thus termed bacteria; double-stranded RNA; lipopolysaccharide from Gram- bacteria; single- g Toll/IL-1R (TIR) homology domain (Akira et al. 2006; Gay and Gangloff stranded RNA; unmethylated DNA from bacteria 1 2007). In follicular-associated epithelial cells in mouse intestine, TLR9 is expressed on the apical and basolateral surfaces of the cell (Abreau 2010). 2 The crystal structure determined for the TLR1-TLR2-lipopeptide complex Inhibitory interactions with TLR7 and TLR9 are known (Wang et al. 2006). 3 reveals that the extracellular TLR domains adopt horse-shoe-like structures Inhibitory interactions with TLR7 are known (Wang et al. 2006). and that the ligand binding site of this TLR heterodimer is located on the con- vex surface of the horse-shoes (Jin et al. 2007). The latter, however, is not nec- essarily a general rule (Akira et al. 2006). The TIR domain is a more or less

30 31 Of these, MyD88 (1 above) and TRIF (3 above) are signaling molecules; hence, the signaling pathways are categorized as MyD88-dependent and TRIF-dependent, respectively. Mal (2 above) and TRAM (4 above) are neces- sary for the recruitment of MyD88 and TRIF, respectively, to TLRs, and are therefore called bridging adapters. SARM (5 above) has got an inhibitory role in TRIF-dependent signaling. MyD88 is used by all TLRs except TLR3, which uses TRIF and no bridging adapter. TLR2 uses Mal, while TLR4 uses all four adapter molecules (Brikos and O’Neill 2008).

Downstream of the TIR domain-containing adapter proteins, signaling ty- pically passes the interleukin-1 receptor-associated kinase 4 and 1 and TNF re- ceptor associated factor 6 complexes. The ultimate result of TLR signaling is activation of transcription factors: either NFκB (subsequent to phosphoryla- tion/degradation of inhibitor of κB), activator protein 1 (via c-Jun N-terminal kinase or protein 38), or one of the interferon regulatory factors. This, in turn, leads to production of pro-inflammatory cytokines, chemokines, and co- stimulatory molecules (Brikos and O’Neill 2008).

Fig. 7 Dimers formed by the TLR1 receptor family, expressed on the cell surface. The TLR2/1 and 4.4 Comparisons to Drosophila Toll and plant TIR TLR2/6 receptor complexes recognize lipopeptides. For TLR10, no ligand has been identified. Seen from an evolutionary point of view and based on the detailed structure of the extracellular domain, TLRs can be divided into vertebrate-like TLRs (V-TLRs) and protostome-like TLRs (P-TLRs) (Hibino et al. 2006). Apply- ing this approach, the Drosophila Toll family includes eight P-TLRs (Toll in- 4.3 TLR signaling cluded) and one V-TLR (Toll-9) (Huang et al. 2008). Moreover, phylogenetic analysis of the TIR domains in the Drosophila Toll family reveals that these – Ligand binding induces dimerisation of the TLR molecules. As a result, Toll-9 being an exception – are more closely related to one another than to conformational changes occur which bring the TIR domains of the two dimeric TLRs in mammals (Bilak et al. 2003). Thus, it is likely that the two groups of partners in close proximity, thus forming a TIR-TIR structure. This structure, receptors – Toll and the Drosophila Toll-related receptors except Toll-9 on the in turn, is the basis for associations with the TIR domain-containing adapter one hand and mammalian TLRs on the other – have evolved independently molecules which link TLRs to downstream intracellular proteins (Brikos and while carrying out different main functions (Bilak et al. 2003). O’Neill 2008). FLS2 in Arabidopsis thaliana, thale cress, is a transmembrane protein in- Five TIR domain-containing adapter molecules are known to interact with volved – like TLR5 – in the recognition of flagelline. FLS2 is equipped with an TLR TIR-TIR structures: extracellular LRR domain. The low degree of sequence similarity between 1. myeloid differentiation primary response protein 88 (MyD88), FLS2 and human TLR5 suggests that these proteins arose independently, as a 2. MyD88-adapter like (Mal), also denoted TIR domain-containing adapter result of convergent evolution. Furthermore, intracellular R proteins containing (TIRAP) β LRR domains fused to TIR domains are present in plants, and a protein re- 3. TIR domain-containing adapter inducing interferon- (TRIF), also denoted sembling CD14, a co-receptor in the TLR4 receptor complex (see above), is TIR-containing adapter molecule-1 (TICAM-1) also known (Nürnberger et al. 2004). In contrast to animal TIRs, plant TIR 4. TRIF-related adapter molecule (TRAM), also denoted TICAM-2 α domains may interact directly with pathogen effectors, while a role in signaling 5. sterile and HEAT-Armadillo motifs (SARM) (Table 1). for plant TIR domains remains to be established. However, some TIR-domain containing R proteins are found in the nucleus and may have roles as transcrip-

32 33 Of these, MyD88 (1 above) and TRIF (3 above) are signaling molecules; hence, the signaling pathways are categorized as MyD88-dependent and TRIF-dependent, respectively. Mal (2 above) and TRAM (4 above) are neces- sary for the recruitment of MyD88 and TRIF, respectively, to TLRs, and are therefore called bridging adapters. SARM (5 above) has got an inhibitory role in TRIF-dependent signaling. MyD88 is used by all TLRs except TLR3, which uses TRIF and no bridging adapter. TLR2 uses Mal, while TLR4 uses all four adapter molecules (Brikos and O’Neill 2008).

Downstream of the TIR domain-containing adapter proteins, signaling ty- pically passes the interleukin-1 receptor-associated kinase 4 and 1 and TNF re- ceptor associated factor 6 complexes. The ultimate result of TLR signaling is activation of transcription factors: either NFκB (subsequent to phosphoryla- tion/degradation of inhibitor of κB), activator protein 1 (via c-Jun N-terminal kinase or protein 38), or one of the interferon regulatory factors. This, in turn, leads to production of pro-inflammatory cytokines, chemokines, and co- stimulatory molecules (Brikos and O’Neill 2008).

Fig. 7 Dimers formed by the TLR1 receptor family, expressed on the cell surface. The TLR2/1 and 4.4 Comparisons to Drosophila Toll and plant TIR TLR2/6 receptor complexes recognize lipopeptides. For TLR10, no ligand has been identified. Seen from an evolutionary point of view and based on the detailed structure of the extracellular domain, TLRs can be divided into vertebrate-like TLRs (V-TLRs) and protostome-like TLRs (P-TLRs) (Hibino et al. 2006). Apply- ing this approach, the Drosophila Toll family includes eight P-TLRs (Toll in- 4.3 TLR signaling cluded) and one V-TLR (Toll-9) (Huang et al. 2008). Moreover, phylogenetic analysis of the TIR domains in the Drosophila Toll family reveals that these – Ligand binding induces dimerisation of the TLR molecules. As a result, Toll-9 being an exception – are more closely related to one another than to conformational changes occur which bring the TIR domains of the two dimeric TLRs in mammals (Bilak et al. 2003). Thus, it is likely that the two groups of partners in close proximity, thus forming a TIR-TIR structure. This structure, receptors – Toll and the Drosophila Toll-related receptors except Toll-9 on the in turn, is the basis for associations with the TIR domain-containing adapter one hand and mammalian TLRs on the other – have evolved independently molecules which link TLRs to downstream intracellular proteins (Brikos and while carrying out different main functions (Bilak et al. 2003). O’Neill 2008). FLS2 in Arabidopsis thaliana, thale cress, is a transmembrane protein in- Five TIR domain-containing adapter molecules are known to interact with volved – like TLR5 – in the recognition of flagelline. FLS2 is equipped with an TLR TIR-TIR structures: extracellular LRR domain. The low degree of sequence similarity between 1. myeloid differentiation primary response protein 88 (MyD88), FLS2 and human TLR5 suggests that these proteins arose independently, as a 2. MyD88-adapter like (Mal), also denoted TIR domain-containing adapter result of convergent evolution. Furthermore, intracellular R proteins containing (TIRAP) β LRR domains fused to TIR domains are present in plants, and a protein re- 3. TIR domain-containing adapter inducing interferon- (TRIF), also denoted sembling CD14, a co-receptor in the TLR4 receptor complex (see above), is TIR-containing adapter molecule-1 (TICAM-1) also known (Nürnberger et al. 2004). In contrast to animal TIRs, plant TIR 4. TRIF-related adapter molecule (TRAM), also denoted TICAM-2 α domains may interact directly with pathogen effectors, while a role in signaling 5. sterile and HEAT-Armadillo motifs (SARM) (Table 1). for plant TIR domains remains to be established. However, some TIR-domain containing R proteins are found in the nucleus and may have roles as transcrip-

32 33 tion factors. Again, these differences in the use of the TIR domain may imply senting any amino acid), leading to the formation of a triple helix structure convergent evolution rather than conservation from an early ancestor of both (Fujita et al. 2004). If a glycine residue is replaced by a more bulky amino acid, animals and plants (Burch-Smith and Dinesh-Kumar 2007). triple helix formation is disturbed. This, in turn, impairs the functionality of the molecule, since this is dependent on oligomerization of the basic trimers (Garred et al. 2003; Garred et al. 2006). 5. Mannan-binding lectin MBL belongs to the protein family of collagenous , which includes multifunctional proteins with defense and house keeping roles (Garred et al. 2009; Lillie et al. 2005). MBL and ficolins activate the lectin pathway, which is the most recently described of the main complement activation pathways, but possibly the most ancient: a glucose-specific lectin which interacts with MASPs has been purified from the sea pineapple, and an ortholog of the mammalian MBLs is present in lamprey (Fujita et al. 2004; Ikeda et al. 1987; Takahashi et al. 2006).

In pigs, as in most mammals, two forms of MBL – MBL-A and MBL-C – are present in the circulation. These are encoded by two separate genes, de- noted MBL1 and MBL2, respectively. Duplication of the MBL gene probably occurred after the divergence of birds and mammals, since only one MBL gene is present in chicken. In humans, MBL1 is non-functional, due to stop codon mutations (Lillie et al. 2005). The protein-coding regions of porcine MBL1 and human MBL2 consist of four exons interrupted by three introns (Garred 2008; Lillie et al. 2006b). Upstream both genes, there are additional exons de- noted exon 0 (Garred 2008; Lillie et al. 2007). In humans, some transcription Fig. 8 The porcine MBL1 gene. Porcine MBL1 (and human MBL2) is comprised of four exons, is initiated from exon 0, but the polypeptides produced are identical to those interrupted by three introns. The cystein-rich domain (C) is encoded by exon 1, the collagen-like translated from transcripts from exon 1 (Garred 2008). domain (CLD) by exons 1 and 2, the neck region (N) by exon 3, and the carbohydrate recognition domain (CRD) by exon 4. The basic trimers oligomerize to form a functional MBL molecule, dis- Single nucleotide polymorphisms (SNPs) affecting MBL concentrations in playing a sertiform (‘bunch of tulips’) formation. serum are known in human MBL2 as well as in porcine MBL1 (Garred 2009; Lillie et al. 2006b). Low-producing MBL alleles are present at varying frequen- The forms and functions of the collagenous lectins vary considerably be- cies in both human (Garred et al. 2009) and porcine (Lillie et al. 2007) popula- tween species. In MBL and ficolins, the CRDs are aggregated at one pole, tions, and concerning human populations, balanced (Bernig et al. 2004) as well which creates their sertiform (‘bunch of tulips’) formation. This facilitates mul- as neutral (Verdu et al. 2006) evolution has been suggested. In vitro experi- tivalent binding to monosaccarides on microbial surfaces. In contrast, in the ments have shown that low levels of MBL in human serum, leading to poor ac- bovine-specific conglutinin and -46, the CRDs are distributed at dia- tivation of the lectin pathway, can be compensated for by the classical pathway metric poles, leading to a cruciform which facilitates agglutination of targets. (Roos et al. 2004). Differences like these may reflect co-evolution with relevant pathogens (Lillie et al. 2005). The basic subunit of MBL is a trimer composed of three identical mono- mers, each consisting of four distinct domains: a cystein-rich N-terminal do- main (C, encoded by exon 1), a collagen-like domain (CLD, encoded by exons 1 and 2), a neck region (N, encoded by exon 3), and a C-terminal carbohydrate recognition domain (CRD, encoded by exon 4) (Fig. 8, Garred 2008). The col- lagen-like domain is composed of numerous G-X-Y motifs (X and Y repre-

34 35 tion factors. Again, these differences in the use of the TIR domain may imply senting any amino acid), leading to the formation of a triple helix structure convergent evolution rather than conservation from an early ancestor of both (Fujita et al. 2004). If a glycine residue is replaced by a more bulky amino acid, animals and plants (Burch-Smith and Dinesh-Kumar 2007). triple helix formation is disturbed. This, in turn, impairs the functionality of the molecule, since this is dependent on oligomerization of the basic trimers (Garred et al. 2003; Garred et al. 2006). 5. Mannan-binding lectin MBL belongs to the protein family of collagenous lectins, which includes multifunctional proteins with defense and house keeping roles (Garred et al. 2009; Lillie et al. 2005). MBL and ficolins activate the lectin pathway, which is the most recently described of the main complement activation pathways, but possibly the most ancient: a glucose-specific lectin which interacts with MASPs has been purified from the sea pineapple, and an ortholog of the mammalian MBLs is present in lamprey (Fujita et al. 2004; Ikeda et al. 1987; Takahashi et al. 2006).

In pigs, as in most mammals, two forms of MBL – MBL-A and MBL-C – are present in the circulation. These are encoded by two separate genes, de- noted MBL1 and MBL2, respectively. Duplication of the MBL gene probably occurred after the divergence of birds and mammals, since only one MBL gene is present in chicken. In humans, MBL1 is non-functional, due to stop codon mutations (Lillie et al. 2005). The protein-coding regions of porcine MBL1 and human MBL2 consist of four exons interrupted by three introns (Garred 2008; Lillie et al. 2006b). Upstream both genes, there are additional exons de- noted exon 0 (Garred 2008; Lillie et al. 2007). In humans, some transcription Fig. 8 The porcine MBL1 gene. Porcine MBL1 (and human MBL2) is comprised of four exons, is initiated from exon 0, but the polypeptides produced are identical to those interrupted by three introns. The cystein-rich domain (C) is encoded by exon 1, the collagen-like translated from transcripts from exon 1 (Garred 2008). domain (CLD) by exons 1 and 2, the neck region (N) by exon 3, and the carbohydrate recognition domain (CRD) by exon 4. The basic trimers oligomerize to form a functional MBL molecule, dis- Single nucleotide polymorphisms (SNPs) affecting MBL concentrations in playing a sertiform (‘bunch of tulips’) formation. serum are known in human MBL2 as well as in porcine MBL1 (Garred 2009; Lillie et al. 2006b). Low-producing MBL alleles are present at varying frequen- The forms and functions of the collagenous lectins vary considerably be- cies in both human (Garred et al. 2009) and porcine (Lillie et al. 2007) popula- tween species. In MBL and ficolins, the CRDs are aggregated at one pole, tions, and concerning human populations, balanced (Bernig et al. 2004) as well which creates their sertiform (‘bunch of tulips’) formation. This facilitates mul- as neutral (Verdu et al. 2006) evolution has been suggested. In vitro experi- tivalent binding to monosaccarides on microbial surfaces. In contrast, in the ments have shown that low levels of MBL in human serum, leading to poor ac- bovine-specific conglutinin and collectin-46, the CRDs are distributed at dia- tivation of the lectin pathway, can be compensated for by the classical pathway metric poles, leading to a cruciform which facilitates agglutination of targets. (Roos et al. 2004). Differences like these may reflect co-evolution with relevant pathogens (Lillie et al. 2005). The basic subunit of MBL is a trimer composed of three identical mono- mers, each consisting of four distinct domains: a cystein-rich N-terminal do- main (C, encoded by exon 1), a collagen-like domain (CLD, encoded by exons 1 and 2), a neck region (N, encoded by exon 3), and a C-terminal carbohydrate recognition domain (CRD, encoded by exon 4) (Fig. 8, Garred 2008). The col- lagen-like domain is composed of numerous G-X-Y motifs (X and Y repre-

34 35

PART III: QUANTITATIVE

TRAITS AND QTL ANALYSIS

Traits influenced by several genes vary continuously and are referred to as quantitative traits. Immune capacity is a typical quantitative trait. The same is true for many traits of primary interest in animal breeding and agriculture, e.g. milk production in cattle and grain yield in crops.

6. Quantitative traits

If a trait is encoded by a single gene occurring as two alleles, of which one is

dominant and one is recessive, and two homozygous individuals are crossed,

the following genotypes will be present in the parental (P), first filial (F1), and

second filial (F2) generations (Fig. 9):

Fig. 9 Mendelian inheritance. The individuals in the P generation are homozygous for the trait of

interest. In the F1 generation, all individuals will be heterozygous. When F1 individuals are inter- crossed, three genotypes – AA, Aa/aA, and aa – will be produced. The AA and Aa/ aA individuals will all exhibit the same (dominant) phenotype.

AA and Aa/aA individuals will all exhibit the same phenotype, namely the

one encoded by the dominant allele; thus, three fourths of the F2 generation

will exhibit the dominant phenotype, while one fourth will exhibit the recessive

phenotype. This is what Gregor Mendel found in his experiments with garden

peas (Mendel s a).

However, if the trait is encoded by two genes, each occurring as a dominant

and a recessive allele, there will be nine genotypes and five phenotypes (Fig. 10):

36 37

PART III: QUANTITATIVE

TRAITS AND QTL ANALYSIS

Traits influenced by several genes vary continuously and are referred to as quantitative traits. Immune capacity is a typical quantitative trait. The same is true for many traits of primary interest in animal breeding and agriculture, e.g. milk production in cattle and grain yield in crops.

6. Quantitative traits

If a trait is encoded by a single gene occurring as two alleles, of which one is

dominant and one is recessive, and two homozygous individuals are crossed,

the following genotypes will be present in the parental (P), first filial (F1), and

second filial (F2) generations (Fig. 9):

Fig. 9 Mendelian inheritance. The individuals in the P generation are homozygous for the trait of

interest. In the F1 generation, all individuals will be heterozygous. When F1 individuals are inter- crossed, three genotypes – AA, Aa/aA, and aa – will be produced. The AA and Aa/ aA individuals will all exhibit the same (dominant) phenotype.

AA and Aa/aA individuals will all exhibit the same phenotype, namely the

one encoded by the dominant allele; thus, three fourths of the F2 generation

will exhibit the dominant phenotype, while one fourth will exhibit the recessive

phenotype. This is what Gregor Mendel found in his experiments with garden

peas (Mendel s a).

However, if the trait is encoded by two genes, each occurring as a dominant

and a recessive allele, there will be nine genotypes and five phenotypes (Fig. 10):

36 37 to a particular quantitative trait of interest can be studied as a separate unit, but when conclusions are drawn, e.g. concerning breeding, the complex interplay between different factors contributing to phenotypes seen in a population must be kept in mind.

In order to determine which genes should be focused on when studying a quantitative trait of interest, a QTL analysis can be applied. QTLs are mapped through their genetic linkage to marker loci (Mackay et al. 2009). Therefore, the prerequisites for a QTL analysis are a linkage map, records of genotypes in the marker loci, and records of phenotypes.

During the first stage (prophase I) of the first meiotic division, replicated homologous chromosomes form pairs. These units, referred to as tetrads, thus consist of two pairs of sister chromatids. In the tetrad, there will be one or more areas – chiasmata – where non-sister chromatids get intertwined. Here, genetic Fig. 10 Nine genotypes and five phenotypes. If a trait is encoded by two genes, each occurring as a material is exchanged, a process referred to as a crossing-over (Fig. 11). dominant and a recessive allele, there will be nine genotypes and five phenotypes. Multiple additive alleles each contribute a portion to the joint phenotype exhibited by the individual.

Here, multiple additive alleles each contribute a portion to the joint pheno- type exhibited by the individual. Many genes influencing a single trait to differ- ent extents, in combination with environmental factors, will bring forth con- tinuous variation in the population. However, the alleles of each of the contrib- uting genes still segregate according to Mendel’s laws. The picture described above is what typically is found concerning immune-related traits. At least some of the traits Gregor Mendel studied – e. g. stem length – are not really encoded by single genes. Since his results still were as if stem length was a monogenic trait, one of his inbred lines must have been homozygous for a re- Fig. 11 Crossing-over. During the first meiotic division, replicated homologous chromosomes form cessive allele with major effect on this trait. tetrads. In the tetrad, there will be one or more areas where non-sister chromatids get intertwined. Here, genetic material is exchanged. A phenotype is influenced not only by genetic factors, but also by the envi- ronment, including conditions intentionally and unintentionally created by In order to construct a linkage map, a family material and a number of po- man. The proportion of genetic versus environmental influence varies between lymorphic marker loci is needed. Two generations of animals are typed for the traits. Heritability estimates are estimates concerning to what degree a certain marker loci, and the result is used to calculate the number of crossing-over trait is influenced by genetic factors. If heritability estimates are low, the re- events that have occurred between these markers during the meiotic division. sponse to selective breeding will be slow. In the opposite case, the status of the This, in turn, will depend on the physical distance between the markers: the trait in a population can be improved by appropriate selection (Falconer and more far apart they are located, the more likely it is that crossing-over events Mackay 1996). will occur (Mackay et al. 2009).

The linkage map will reveal recombination distances between and orders of loci in a chromosomal region of interest. In the next step, this information is 7. QTL analysis used, together with records of genotypes at marker loci and records of pheno- A complex picture emerges, where traits are influenced by multiple genes to types, in a detection and localization process. To detect a QTL, large allelic ef- different extents and by multiple environmental factors. Each gene contributing fects and allele frequencies close to 0.5 will help. To determine its localization

38 39 to a particular quantitative trait of interest can be studied as a separate unit, but when conclusions are drawn, e.g. concerning breeding, the complex interplay between different factors contributing to phenotypes seen in a population must be kept in mind.

In order to determine which genes should be focused on when studying a quantitative trait of interest, a QTL analysis can be applied. QTLs are mapped through their genetic linkage to marker loci (Mackay et al. 2009). Therefore, the prerequisites for a QTL analysis are a linkage map, records of genotypes in the marker loci, and records of phenotypes.

During the first stage (prophase I) of the first meiotic division, replicated homologous chromosomes form pairs. These units, referred to as tetrads, thus consist of two pairs of sister chromatids. In the tetrad, there will be one or more areas – chiasmata – where non-sister chromatids get intertwined. Here, genetic Fig. 10 Nine genotypes and five phenotypes. If a trait is encoded by two genes, each occurring as a material is exchanged, a process referred to as a crossing-over (Fig. 11). dominant and a recessive allele, there will be nine genotypes and five phenotypes. Multiple additive alleles each contribute a portion to the joint phenotype exhibited by the individual.

Here, multiple additive alleles each contribute a portion to the joint pheno- type exhibited by the individual. Many genes influencing a single trait to differ- ent extents, in combination with environmental factors, will bring forth con- tinuous variation in the population. However, the alleles of each of the contrib- uting genes still segregate according to Mendel’s laws. The picture described above is what typically is found concerning immune-related traits. At least some of the traits Gregor Mendel studied – e. g. stem length – are not really encoded by single genes. Since his results still were as if stem length was a monogenic trait, one of his inbred lines must have been homozygous for a re- Fig. 11 Crossing-over. During the first meiotic division, replicated homologous chromosomes form cessive allele with major effect on this trait. tetrads. In the tetrad, there will be one or more areas where non-sister chromatids get intertwined. Here, genetic material is exchanged. A phenotype is influenced not only by genetic factors, but also by the envi- ronment, including conditions intentionally and unintentionally created by In order to construct a linkage map, a family material and a number of po- man. The proportion of genetic versus environmental influence varies between lymorphic marker loci is needed. Two generations of animals are typed for the traits. Heritability estimates are estimates concerning to what degree a certain marker loci, and the result is used to calculate the number of crossing-over trait is influenced by genetic factors. If heritability estimates are low, the re- events that have occurred between these markers during the meiotic division. sponse to selective breeding will be slow. In the opposite case, the status of the This, in turn, will depend on the physical distance between the markers: the trait in a population can be improved by appropriate selection (Falconer and more far apart they are located, the more likely it is that crossing-over events Mackay 1996). will occur (Mackay et al. 2009).

The linkage map will reveal recombination distances between and orders of loci in a chromosomal region of interest. In the next step, this information is 7. QTL analysis used, together with records of genotypes at marker loci and records of pheno- A complex picture emerges, where traits are influenced by multiple genes to types, in a detection and localization process. To detect a QTL, large allelic ef- different extents and by multiple environmental factors. Each gene contributing fects and allele frequencies close to 0.5 will help. To determine its localization

38 39 to the narrowest possible chromosomal region, many individuals and many marker loci are needed (Mackay et al. 2009). PART IV: MUTATIONS, SE-

LECTION, AND EVOLUTION

The word ‘polymorphism’ – literally ‘many forms’ – refers to variations within a DNA sequence and the phenotypic manifestations of these variations. Polymorphisms become evident when sequences or observations of phenotypes from different individuals are compared. Different versions of a gene or a mi- crosatellite (ms) as well as alternative nucleotides at a particular position in the genome are referred to as alleles.

8. Mutations

Mutations, together with phenomena like gene duplication and whole ge-

nome duplication, are the ultimate source of variation in nature and the fuel for

evolution. There are several kinds of mutations, causing different types of

changes within the DNA sequence. Recombination refers to the replacement

of a sequence by another, while insertions and deletions implicate that one or

more nucleotides have been added or removed from the sequence, respectively.

The term inversion signifies the rotation by 180° of a double-stranded segment

of at least two base pairs.

Point mutations, or substitution mutations, refers to the replacement of a

single nucleotide by another. Point mutations give rise to SNPs. Point muta-

tions occur at different frequencies in different kinds of sequences: they are

more common in pseudogenes and introns and less common in coding se-

quences and promoter regions. Also, transitions (i.e. when a purine is ex-

changed for a purine or a pyrimidine for a pyrimidine) are more common than

transversions (i.e. when a purine is exchanged for a pyrimidine, or vice versa).

Point mutations in coding sequences are said to be synonymous when they

leave the amino acid sequence unaffected; otherwise, they are described as non-

synonymous (Graur and Li 2000).

As mentioned previously, there are other phenomena besides mutations

which fuel evolution. In connection with the development of the immune de-

fense, transposition and whole genome duplication have played important roles. As stated above, the origin of the mammalian adaptive immune system is con- ventionally linked to the emergence of the RAG genes. The evolutionary origin of these genes is thought to be a transposon – a mobile genetic element – which invaded an exon belonging to the immunoglobulin superfamily. Furthermore, it

40 41 to the narrowest possible chromosomal region, many individuals and many marker loci are needed (Mackay et al. 2009). PART IV: MUTATIONS, SE-

LECTION, AND EVOLUTION

The word ‘polymorphism’ – literally ‘many forms’ – refers to variations within a DNA sequence and the phenotypic manifestations of these variations. Polymorphisms become evident when sequences or observations of phenotypes from different individuals are compared. Different versions of a gene or a mi- crosatellite (ms) as well as alternative nucleotides at a particular position in the genome are referred to as alleles.

8. Mutations

Mutations, together with phenomena like gene duplication and whole ge-

nome duplication, are the ultimate source of variation in nature and the fuel for

evolution. There are several kinds of mutations, causing different types of

changes within the DNA sequence. Recombination refers to the replacement

of a sequence by another, while insertions and deletions implicate that one or

more nucleotides have been added or removed from the sequence, respectively.

The term inversion signifies the rotation by 180° of a double-stranded segment

of at least two base pairs.

Point mutations, or substitution mutations, refers to the replacement of a

single nucleotide by another. Point mutations give rise to SNPs. Point muta-

tions occur at different frequencies in different kinds of sequences: they are

more common in pseudogenes and introns and less common in coding se-

quences and promoter regions. Also, transitions (i.e. when a purine is ex-

changed for a purine or a pyrimidine for a pyrimidine) are more common than

transversions (i.e. when a purine is exchanged for a pyrimidine, or vice versa).

Point mutations in coding sequences are said to be synonymous when they

leave the amino acid sequence unaffected; otherwise, they are described as non-

synonymous (Graur and Li 2000).

As mentioned previously, there are other phenomena besides mutations

which fuel evolution. In connection with the development of the immune de-

fense, transposition and whole genome duplication have played important roles. As stated above, the origin of the mammalian adaptive immune system is con- ventionally linked to the emergence of the RAG genes. The evolutionary origin of these genes is thought to be a transposon – a mobile genetic element – which invaded an exon belonging to the immunoglobulin superfamily. Furthermore, it

40 41 is believed that two rounds of whole genome duplication were crucial prerequi- Purifying (or negative) selection protects the function of a protein. When sites for the emergence of the adaptive immune system as seen in today’s deleterious alleles are removed, linked neutral alleles may also be removed, par- mammals (Flajnik and Kasahara 2010). ticularly from genomic regions where recombination is infrequent. This phe- nomenon is known as background selection and reduces diversity (Barreiro and Quitana-Murci 2010). 9. Selection and evolution Positive (or directional) selection makes use of mutations which improve Natural selection is defined as the differential reproduction of genetically the fitness of its carrier. It seems that there are a number of genes encoding de- distinct individuals or genotypes within a population. Fitness is a measure of an fense and/or immunity proteins which have evolved rapidly since the diver- individual’s ability to survive and reproduce. The size of a population is limited gence of humans and chimpanzees (Barreiro and Quitana-Murci 2010). This by the carrying capacity of its environment; therefore, the fitness of a certain indicates positive selection and is not surprising, considering the arms’ race type individual is determined in comparison to the other individuals in the popula- of relationship between defense/immunity proteins and microorganisms. Fur- tion. Most new mutations that arise in coding sequences are deleterious, i.e. re- thermore, immune-related genes are overrepresented when recently occurring duce the relative reproductive success of the carrier. These mutations will even- selective events in human populations are considered (Barreiro and Quitana- tually be removed from the population by purifying selection. Rarely, a muta- Murci 2010). This suggests that the human immune defense has been particu- tion arises which is advantageous, i.e. increases the relative reproductive success larly challenged during our most recent past, which, in turn, could be linked to of the carrier; such a mutation may be passed on to the next generation and the massive increase of human population sizes when farming was invented eventually bring about changes in allele frequencies in the population. The fre- and, possibly, exposure to new zoonoses when animals were domesticated. quencies of mutations which are neutral from the point of view of relative re- When a positively selected allele increases in frequency in the population, productive success will not be affected by natural selection (Graur and Li linked neutral alleles will be dragged along with it. Thus, variation that is not 2000). Evolution was originally defined as descent with modification (Darwin associated with the positively selected allele will disappear and the overall ge- 2000), and is currently defined as changes in allele frequencies over time (Fig. netic variation in the population will decrease (Barreiro and Quitana-Murci 12, Graur and Li 2000). 2010).

Balancing selection maintains diversity in a population by two mechanisms: heterozygote advantage and frequency-dependent selection. Heterozygote ad- vantage refers to a situation when the heterozygous genotype confers a greater fitness to its carrier than either of the two homozygous genotypes. The most convincing example found so far is connected with the beta globin gene (HBB): the HbS allele, which causes sickle cell anemia, also protects against malaria, which explains the presence of this allele at high frequencies in Africa. Fre- quency-dependent selection occurs when the fitness of a certain phenotype is dependent on its frequency compared to the frequency of other phenotypes in a population (Barreiro and Quitana-Murci 2010).

It has been shown that TLR2 evolves under purifying as well as positive se- lection (Werling et al. 2009). This is not surprising: a balance between both types of selection would be optimal in order to insure protection of function and at the same time make use of the dynamics brought about by mutations in order to cope with novel virulence potential continuously achieved by microor- Fig. 12 Evolution. Evolution is defined as change in allele frequencies over time. As time passes, ganisms. In consistency with this, the extracellular domain of TLR2, which in- some alleles may disappear and new ones may arise. teracts with ligands, evolves under positive selection, while the TIR domain is subjected to purifying selection (Werling et al. 2009). Moreover, there is accu- mulating evidence of a species-specific component in TLR function, which, in

42 43 is believed that two rounds of whole genome duplication were crucial prerequi- Purifying (or negative) selection protects the function of a protein. When sites for the emergence of the adaptive immune system as seen in today’s deleterious alleles are removed, linked neutral alleles may also be removed, par- mammals (Flajnik and Kasahara 2010). ticularly from genomic regions where recombination is infrequent. This phe- nomenon is known as background selection and reduces diversity (Barreiro and Quitana-Murci 2010). 9. Selection and evolution Positive (or directional) selection makes use of mutations which improve Natural selection is defined as the differential reproduction of genetically the fitness of its carrier. It seems that there are a number of genes encoding de- distinct individuals or genotypes within a population. Fitness is a measure of an fense and/or immunity proteins which have evolved rapidly since the diver- individual’s ability to survive and reproduce. The size of a population is limited gence of humans and chimpanzees (Barreiro and Quitana-Murci 2010). This by the carrying capacity of its environment; therefore, the fitness of a certain indicates positive selection and is not surprising, considering the arms’ race type individual is determined in comparison to the other individuals in the popula- of relationship between defense/immunity proteins and microorganisms. Fur- tion. Most new mutations that arise in coding sequences are deleterious, i.e. re- thermore, immune-related genes are overrepresented when recently occurring duce the relative reproductive success of the carrier. These mutations will even- selective events in human populations are considered (Barreiro and Quitana- tually be removed from the population by purifying selection. Rarely, a muta- Murci 2010). This suggests that the human immune defense has been particu- tion arises which is advantageous, i.e. increases the relative reproductive success larly challenged during our most recent past, which, in turn, could be linked to of the carrier; such a mutation may be passed on to the next generation and the massive increase of human population sizes when farming was invented eventually bring about changes in allele frequencies in the population. The fre- and, possibly, exposure to new zoonoses when animals were domesticated. quencies of mutations which are neutral from the point of view of relative re- When a positively selected allele increases in frequency in the population, productive success will not be affected by natural selection (Graur and Li linked neutral alleles will be dragged along with it. Thus, variation that is not 2000). Evolution was originally defined as descent with modification (Darwin associated with the positively selected allele will disappear and the overall ge- 2000), and is currently defined as changes in allele frequencies over time (Fig. netic variation in the population will decrease (Barreiro and Quitana-Murci 12, Graur and Li 2000). 2010).

Balancing selection maintains diversity in a population by two mechanisms: heterozygote advantage and frequency-dependent selection. Heterozygote ad- vantage refers to a situation when the heterozygous genotype confers a greater fitness to its carrier than either of the two homozygous genotypes. The most convincing example found so far is connected with the beta globin gene (HBB): the HbS allele, which causes sickle cell anemia, also protects against malaria, which explains the presence of this allele at high frequencies in Africa. Fre- quency-dependent selection occurs when the fitness of a certain phenotype is dependent on its frequency compared to the frequency of other phenotypes in a population (Barreiro and Quitana-Murci 2010).

It has been shown that TLR2 evolves under purifying as well as positive se- lection (Werling et al. 2009). This is not surprising: a balance between both types of selection would be optimal in order to insure protection of function and at the same time make use of the dynamics brought about by mutations in order to cope with novel virulence potential continuously achieved by microor- Fig. 12 Evolution. Evolution is defined as change in allele frequencies over time. As time passes, ganisms. In consistency with this, the extracellular domain of TLR2, which in- some alleles may disappear and new ones may arise. teracts with ligands, evolves under positive selection, while the TIR domain is subjected to purifying selection (Werling et al. 2009). Moreover, there is accu- mulating evidence of a species-specific component in TLR function, which, in

42 43 turn, might reflect co-evolution with relevant microorganisms (Werling et al. 2009). Concerning TLR10 and MBL, neutral evolution – guided solely by mi- gration and genetic drift and suggesting redundancy – has been suggested (Bar- PART V: MATERIALS AND reiro et al. 2009; Lillie et al. 2005; Verdu et al. 2006). METHODS Breeding is artificial selection carried out by man. Concerning pigs, it is likely that selection for behavioral traits, like non-aggressiveness, took place during the early stages of domestication. Moreover, it has been established that the variable coat colors exhibited by domestic pigs are the results of man’s in- 11. Pigs in genetic analyses tentional selection, as opposed to relaxed purifying selection or mutations af- Pigs belong to the order of artiodactyls, which also includes cattle, sheep, fecting both coat color and other traits subjected to strong selection (Fang et al. and goats (Graves 1984). It has been established that the wild boar lineages in 2009). Also, selection for growth rate, i.e. meat production, is standard proce- Europe and Asia separated about 900,000 years ago and that domestication dure in farming. Research efforts aiming to promote selective breeding for im- took place independently in Europe and Asia approximately 9,000 years ago mune response is carried out within the disciplines of animal genetics and vet- (Kijas and Andersson 2001). Moreover, microsatellite data as well as analysis of erinary medicine (Nicholas 2005). the mitochondrial genome indicate that Asian domestic pigs trace their origin

to a more diverse wild boar population than their European counterparts (Lar- son et al. 2005; Megens et al. 2008); this, in turn, indicates a bottle neck in the 10. Wild boars and domestic pigs European wild boar population before domestication started, resulting in a Due to differences in selective (natural and artificial) pressure, differences in lower degree of diversity in European pigs than in Asian pigs in wild and do- polymorphic patterns between wild boars on the one hand and domestic pigs mestic populations alike. on the other are to be expected. For some genes studied, wild boars and domes- In large parts of Europe, the wild boar disappeared by the end of the 17th tic pigs differ with regard to allele frequencies (Giuffra et al. 2000; Iso-Touru th et al. 2009), while for genes exposed to strong selection, such as MC1R (Fang century, but their numbers increased again during the latter half of the 19 et al. 2009), wild and domestic populations exhibit completely different alleles. century (Ruiz-Fons et al. 2008). The present Swedish wild boar population Differences in population sizes and local bottle necks will also affect diversity. was founded by a very limited number of animals which escaped from enclo- sures in the 1970s and -80s (personal communication, Alf Lettesjö, Henrik Ba- rnekow, Svenska Jägareförbundet, Öster Malma SE-611 91 Nyköping, Swe- den). Thus, traces of a founder effect in these animals are to be expected. The Austrian wild boars sampled in papers II and III originate from a large enclo- sure (approximately 1,200 ha), established in the 1750s and inhabited by roughly 800 animals (personal communication, Dr. Sabine Hammer, Univer- sity of Veterinary Medicine Vienna, Austria). In contrast, the Czech wild boars, all from the surroundings of Brno, have no history of captivity. On rare occasions, hunters in central Moravia come across wild boars exhibiting a spot- ted coat color phenotype, intermediate between the wild boar and domestic pig phenotypes. These animals are assumed to be the result of crossbreeding be- tween wild boars and freely roaming domestic pigs. However, none of the ani- mals sampled in papers II and III belong to this intermediate phenotype (per- sonal communication, Prof. Aleš Knoll, Mendel University of Agriculture and Forestry Brno, Czech Republic).

It is known that the Large White breed historically has been subjected to strong influence from Asian phenotypes (Giuffra et al. 2000). In paper II, this

44 45 turn, might reflect co-evolution with relevant microorganisms (Werling et al. 2009). Concerning TLR10 and MBL, neutral evolution – guided solely by mi- gration and genetic drift and suggesting redundancy – has been suggested (Bar- PART V: MATERIALS AND reiro et al. 2009; Lillie et al. 2005; Verdu et al. 2006). METHODS Breeding is artificial selection carried out by man. Concerning pigs, it is likely that selection for behavioral traits, like non-aggressiveness, took place during the early stages of domestication. Moreover, it has been established that the variable coat colors exhibited by domestic pigs are the results of man’s in- 11. Pigs in genetic analyses tentional selection, as opposed to relaxed purifying selection or mutations af- Pigs belong to the order of artiodactyls, which also includes cattle, sheep, fecting both coat color and other traits subjected to strong selection (Fang et al. and goats (Graves 1984). It has been established that the wild boar lineages in 2009). Also, selection for growth rate, i.e. meat production, is standard proce- Europe and Asia separated about 900,000 years ago and that domestication dure in farming. Research efforts aiming to promote selective breeding for im- took place independently in Europe and Asia approximately 9,000 years ago mune response is carried out within the disciplines of animal genetics and vet- (Kijas and Andersson 2001). Moreover, microsatellite data as well as analysis of erinary medicine (Nicholas 2005). the mitochondrial genome indicate that Asian domestic pigs trace their origin

to a more diverse wild boar population than their European counterparts (Lar- son et al. 2005; Megens et al. 2008); this, in turn, indicates a bottle neck in the 10. Wild boars and domestic pigs European wild boar population before domestication started, resulting in a Due to differences in selective (natural and artificial) pressure, differences in lower degree of diversity in European pigs than in Asian pigs in wild and do- polymorphic patterns between wild boars on the one hand and domestic pigs mestic populations alike. on the other are to be expected. For some genes studied, wild boars and domes- In large parts of Europe, the wild boar disappeared by the end of the 17th tic pigs differ with regard to allele frequencies (Giuffra et al. 2000; Iso-Touru th et al. 2009), while for genes exposed to strong selection, such as MC1R (Fang century, but their numbers increased again during the latter half of the 19 et al. 2009), wild and domestic populations exhibit completely different alleles. century (Ruiz-Fons et al. 2008). The present Swedish wild boar population Differences in population sizes and local bottle necks will also affect diversity. was founded by a very limited number of animals which escaped from enclo- sures in the 1970s and -80s (personal communication, Alf Lettesjö, Henrik Ba- rnekow, Svenska Jägareförbundet, Öster Malma SE-611 91 Nyköping, Swe- den). Thus, traces of a founder effect in these animals are to be expected. The Austrian wild boars sampled in papers II and III originate from a large enclo- sure (approximately 1,200 ha), established in the 1750s and inhabited by roughly 800 animals (personal communication, Dr. Sabine Hammer, Univer- sity of Veterinary Medicine Vienna, Austria). In contrast, the Czech wild boars, all from the surroundings of Brno, have no history of captivity. On rare occasions, hunters in central Moravia come across wild boars exhibiting a spot- ted coat color phenotype, intermediate between the wild boar and domestic pig phenotypes. These animals are assumed to be the result of crossbreeding be- tween wild boars and freely roaming domestic pigs. However, none of the ani- mals sampled in papers II and III belong to this intermediate phenotype (per- sonal communication, Prof. Aleš Knoll, Mendel University of Agriculture and Forestry Brno, Czech Republic).

It is known that the Large White breed historically has been subjected to strong influence from Asian phenotypes (Giuffra et al. 2000). In paper II, this

44 45 is demonstrated by the haplotype analyses concerning TLR2 and TLR6, where drawbacks are compensated for by their multiallelic nature, as opposed to the haplotypes previously found in Asian breeds were detected (Shinkai et al. biallelic nature of SNPs. 2006b). Furthermore, developed mainly in the United States – possibly from a comparatively narrow breeding basis – the Hampshire breed can be expected to show signs of an origin differing from that of the European Large White and 12.1.3 Interpretation of signals Landrace pig lines. This is demonstrated by the haplotype analysis concerning the TLR6-1-10 gene cluster in paper III. Background noise and replication slippage can be an issue when interpret- ing signals from microsatellites. However, by comparing signals from the off- Domestic animals add extra dimensions to genetic studies. First, they en- spring to those from the P generation, this difficulty can be overcome. able comparisons of populations exposed to different selective criteria and envi- ronmental challenges. Second, by making experimental crosses and using the backcross approach, it is possible to deduce the origin of a certain allele. Third, 12.1.4 CRI-MAP in the case of the pig, comparisons to the wild ancestor – i.e. the wild boar – Based on recombination (crossing-over) events during meiosis and the seg- can shed light on the evolutionary process. Moreover, when constructing family regation of alleles between two generations of animals, CRI-MAP builds link- materials comprising several generations, the short generation time of the pig is age maps along chromosomal regions of interest, expressing distances between convenient, and large litters make it probable that most of the possible genetic markers in cM. LOD scores indicate the probability that two markers are combinations from the P generation will be seen in the offspring. In addition, linked in a particular pedigree (Green et al. 1990). pigs are large animal models for humans.

12.1.5 QTL Express 12. Methodological considerations QTL Express is developed for detection of QTLs in halfsib outbred popu- lations and in F2 populations derived from crosses between either inbred or 12.1 Paper I outbred lines. Three input files, containing a marker map, marker genotypes, and trait data, are used. The mapping procedure has two steps: the identity-by- descent (IBD) probabilities at specific chromosomal locations are determined, 12.1.1 Duplex PCR reactions using a touch-down pro- and a statistical model is fitted to the IBD coefficients and phenotypic data. cedure Permutation tests, to set significance thresholds, and a boot-strapping proce- Duplex – or even triplex – polymerase chain reactions (PCRs) are com- dure, to estimate the confidence interval of a QTL location, are also performed monly used in connection with ms analysis, in order to achieve effectiveness. By (Seaton et al. 2002). means of a touch-down procedure, it is possible to utilize primer pairs with slightly different melting temperatures in the same reaction. Prior to the ampli- fying PCR module, the sample is run through a pre-module, where the anneal- 12.1.6 Limitations in connection with QTL analysis ing temperature is lowered by each cycle. Applying the highest annealing tem- QTL mapping is a slow and labor-intense procedure. Chromosomal regions perature to the first cycle is a measure to ensure specificity. of interest will be revealed, but pinpointing of specific causal genes is not possi-

ble (Mackay et al. 2009). Instead, in following steps, multiple candidate genes have to be evaluated. Moreover, different linkage maps will give slightly differ- 12.1.2 Microsatellites as markers ent results. Furthermore, variation in transcript abundance and other interme- Ms, short sequences comprised of repeated nucleotide motifs, are com- diate molecular phenotypes are not accounted for in a traditional QTL analysis monly used as markers when linkage maps are constructed. Ms are not as fre- (Mackay et al. 2009). quent in the genome as SNPs, and, at least in humans and mice, they are more common towards the ends of chromosomes (Ellegren 2004). However, these

46 47 is demonstrated by the haplotype analyses concerning TLR2 and TLR6, where drawbacks are compensated for by their multiallelic nature, as opposed to the haplotypes previously found in Asian breeds were detected (Shinkai et al. biallelic nature of SNPs. 2006b). Furthermore, developed mainly in the United States – possibly from a comparatively narrow breeding basis – the Hampshire breed can be expected to show signs of an origin differing from that of the European Large White and 12.1.3 Interpretation of signals Landrace pig lines. This is demonstrated by the haplotype analysis concerning the TLR6-1-10 gene cluster in paper III. Background noise and replication slippage can be an issue when interpret- ing signals from microsatellites. However, by comparing signals from the off- Domestic animals add extra dimensions to genetic studies. First, they en- spring to those from the P generation, this difficulty can be overcome. able comparisons of populations exposed to different selective criteria and envi- ronmental challenges. Second, by making experimental crosses and using the backcross approach, it is possible to deduce the origin of a certain allele. Third, 12.1.4 CRI-MAP in the case of the pig, comparisons to the wild ancestor – i.e. the wild boar – Based on recombination (crossing-over) events during meiosis and the seg- can shed light on the evolutionary process. Moreover, when constructing family regation of alleles between two generations of animals, CRI-MAP builds link- materials comprising several generations, the short generation time of the pig is age maps along chromosomal regions of interest, expressing distances between convenient, and large litters make it probable that most of the possible genetic markers in cM. LOD scores indicate the probability that two markers are combinations from the P generation will be seen in the offspring. In addition, linked in a particular pedigree (Green et al. 1990). pigs are large animal models for humans.

12.1.5 QTL Express 12. Methodological considerations QTL Express is developed for detection of QTLs in halfsib outbred popu- lations and in F2 populations derived from crosses between either inbred or 12.1 Paper I outbred lines. Three input files, containing a marker map, marker genotypes, and trait data, are used. The mapping procedure has two steps: the identity-by- descent (IBD) probabilities at specific chromosomal locations are determined, 12.1.1 Duplex PCR reactions using a touch-down pro- and a statistical model is fitted to the IBD coefficients and phenotypic data. cedure Permutation tests, to set significance thresholds, and a boot-strapping proce- Duplex – or even triplex – polymerase chain reactions (PCRs) are com- dure, to estimate the confidence interval of a QTL location, are also performed monly used in connection with ms analysis, in order to achieve effectiveness. By (Seaton et al. 2002). means of a touch-down procedure, it is possible to utilize primer pairs with slightly different melting temperatures in the same reaction. Prior to the ampli- fying PCR module, the sample is run through a pre-module, where the anneal- 12.1.6 Limitations in connection with QTL analysis ing temperature is lowered by each cycle. Applying the highest annealing tem- QTL mapping is a slow and labor-intense procedure. Chromosomal regions perature to the first cycle is a measure to ensure specificity. of interest will be revealed, but pinpointing of specific causal genes is not possi-

ble (Mackay et al. 2009). Instead, in following steps, multiple candidate genes have to be evaluated. Moreover, different linkage maps will give slightly differ- 12.1.2 Microsatellites as markers ent results. Furthermore, variation in transcript abundance and other interme- Ms, short sequences comprised of repeated nucleotide motifs, are com- diate molecular phenotypes are not accounted for in a traditional QTL analysis monly used as markers when linkage maps are constructed. Ms are not as fre- (Mackay et al. 2009). quent in the genome as SNPs, and, at least in humans and mice, they are more common towards the ends of chromosomes (Ellegren 2004). However, these

46 47 12.2 Papers II and III relatively insensitive to the selection coefficient (Kryazhimskiy and Plotkin 2008). Defined as fertile offspring produced by interbreeding, wild boars and domestic pigs are presently a single species. Moreover, to what extent the 12.2.1 Direct sequencing populations are separated is not known, and genotypes are not fixed. Therefore, For papers II and III, a stepwise sequencing strategy was chosen. Using in the context described in this thesis, the dN/dS statistic must be interpreted primer pairs designed to produce overlapping fragments, it is possible to obtain with great care. an unambiguous sequence covering an extended chromosomal region of inter- est, e.g. a gene. 12.3 Paper IV The main issue connected with direct sequencing is impurities in the PCR product, which will create noise. This can be overcome by careful optimization and, if necessary, by applying a touch-down procedure to the PCR reaction. 12.3.1 Sequencing subsequent to cloning When interpreting sequence files, especially at positions displaying mixed se- By ligating a PCR fragment into a vector and transforming competent bac- quence, it is advisable that a single, well experienced person is responsible, in teria, it is possible to obtain effective amplification of the PCR fragment of in- order to avoid arbitrary results. Sequencing of multiple PCR products also en- terest. However, when this approach is applied in connection with diploid or- hances the quality of the interpretation. ganisms, several clones must be sequenced in order to determine the genotype. Furthermore, a certain risk for mutations during the cloning procedure must be taken into account. However, by using reference sequences during the typing 12.2.2 Protein prediction procedure, this problem can be minimized.

The simple modular architecture research tool, SMART, can be used to predict protein domains from an amino acid sequence. Even though it was re- cently expanded, in terms of the number of protein domain models included 12.3.2 ELISA and size of the underlying database (Letunic et al. 2008), SMART does not Detection with a monoclonal antibody is a measure to avoid false positives. necessarily recognise all domains of a certain type (Matsushima et al. 2007), ir- By the usage of a biotinylated antibody in combination with streptavidin-HRP, regularities being an issue. Therefore, results from SMART predictions should detection limits are lowered. be interpreted with care.

12.2.3 Mann-Whitney U-test

In order to determine whether the mean values of samples from two popu- lations differ from one another, a t-test can be applied. However, a t-test pre- supposes normal distribution of the samples from each population. When nor- mal distribution can not be assumed, a Mann-Whitney U-test can be used in- stead.

12.2.4 The dN/dS statistic For quantification of evolutionary pressure, the ratio of substitution rates at non-synonymous and synonymous polymorphic positions is commonly used. However, the dN/dS statistic was originally developed for analysis of sequences from divergent species. When applied to samples from a single population, it is

48 49 12.2 Papers II and III relatively insensitive to the selection coefficient (Kryazhimskiy and Plotkin 2008). Defined as fertile offspring produced by interbreeding, wild boars and domestic pigs are presently a single species. Moreover, to what extent the 12.2.1 Direct sequencing populations are separated is not known, and genotypes are not fixed. Therefore, For papers II and III, a stepwise sequencing strategy was chosen. Using in the context described in this thesis, the dN/dS statistic must be interpreted primer pairs designed to produce overlapping fragments, it is possible to obtain with great care. an unambiguous sequence covering an extended chromosomal region of inter- est, e.g. a gene. 12.3 Paper IV The main issue connected with direct sequencing is impurities in the PCR product, which will create noise. This can be overcome by careful optimization and, if necessary, by applying a touch-down procedure to the PCR reaction. 12.3.1 Sequencing subsequent to cloning When interpreting sequence files, especially at positions displaying mixed se- By ligating a PCR fragment into a vector and transforming competent bac- quence, it is advisable that a single, well experienced person is responsible, in teria, it is possible to obtain effective amplification of the PCR fragment of in- order to avoid arbitrary results. Sequencing of multiple PCR products also en- terest. However, when this approach is applied in connection with diploid or- hances the quality of the interpretation. ganisms, several clones must be sequenced in order to determine the genotype. Furthermore, a certain risk for mutations during the cloning procedure must be taken into account. However, by using reference sequences during the typing 12.2.2 Protein prediction procedure, this problem can be minimized.

The simple modular architecture research tool, SMART, can be used to predict protein domains from an amino acid sequence. Even though it was re- cently expanded, in terms of the number of protein domain models included 12.3.2 ELISA and size of the underlying database (Letunic et al. 2008), SMART does not Detection with a monoclonal antibody is a measure to avoid false positives. necessarily recognise all domains of a certain type (Matsushima et al. 2007), ir- By the usage of a biotinylated antibody in combination with streptavidin-HRP, regularities being an issue. Therefore, results from SMART predictions should detection limits are lowered. be interpreted with care.

12.2.3 Mann-Whitney U-test

In order to determine whether the mean values of samples from two popu- lations differ from one another, a t-test can be applied. However, a t-test pre- supposes normal distribution of the samples from each population. When nor- mal distribution can not be assumed, a Mann-Whitney U-test can be used in- stead.

12.2.4 The dN/dS statistic For quantification of evolutionary pressure, the ratio of substitution rates at non-synonymous and synonymous polymorphic positions is commonly used. However, the dN/dS statistic was originally developed for analysis of sequences from divergent species. When applied to samples from a single population, it is

48 49

AIMS

The overall aim of this thesis has been to explore polymorphic patterns in PRR

genes in pigs. More specifically, the objectives were

 to further develop an existing linkage map and refine previously identi-

fied QTLs for immune-related traits on pig chromosome 8

 to analyze polymorphic patterns in the TLR1 gene family on pig

chromosome 8

 to compare polymorphic patterns in European wild boars and domes-

tic pigs

 to analyze MBL-A serum concentrations and MBL1 genotypes in

European wild boars, Large White pigs, and wild boar/Large White

crossbreds.

50 51

AIMS

The overall aim of this thesis has been to explore polymorphic patterns in PRR

genes in pigs. More specifically, the objectives were

 to further develop an existing linkage map and refine previously identi-

fied QTLs for immune-related traits on pig chromosome 8

 to analyze polymorphic patterns in the TLR1 gene family on pig

chromosome 8

 to compare polymorphic patterns in European wild boars and domes-

tic pigs

 to analyze MBL-A serum concentrations and MBL1 genotypes in

European wild boars, Large White pigs, and wild boar/Large White

crossbreds.

50 51

PART VI: PAPERS, CONCLU-

SIONS, FUTURE

13. Paper I

13.1 Objectives The objective of the paper I study was to refine QTLs for immune-related traits previously identified on pig chromosome 8 (SSC8), close to the mast/stem cell growth factor receptor (KIT) gene and the ms S0225, respec- tively. The study was based on 47 BC2 animals, originating from 2 wild boar x Large White intercrosses, and 191 Landrace/Large White crossbreds, born in a specific pathogen-free herd (Wattrang et al. 2005). The variation in origin of the animals used can be expected to answer to variation in immune-related genes, thus making these pigs suitable for a QTL analysis. For the present analysis, the ms SJA07, SJA10, and SJA11, located within the TLR6-1-10 gene cluster, and HIK10, residing 99 bp upstreams of TLR2, were added to the link- age map.

13.2 Results and discussion

The results reported in paper I confirm the QTLs previously identified (Wattrang et al. 2005), and points at additional regions in their vicinity, possi- bly influencing the traits studied.

As for IgM+ cells, the previously identified QTL was split into two, one lo- cated at the position of KIT and the other residing slightly upstream of this lo- cation (Fig. 13). A potential QTL for ConA induced proliferation and a QTL for PWM induced IL-2 production were both confirmed, close to S0069 and S0225, respectively. Furthermore, for PWM induced proliferation, a significant QTL was detected slightly upstream of KIT. The QTL for this trait detected close to S0225 in the previous study (Wattrang et al. 2005) was not seen in the present analysis, but the F-value was increasing at the corresponding location. Moreover, the increasing F-values at S0225 for hemoglobin and hematocrit levels, as well as for MHCII+ cells, were confirmed. Finally, a suggestive QTL for numbers of segmented neutrophils, not found in the previous study, was lo- cated close to the position of TLR2.

52 53

PART VI: PAPERS, CONCLU-

SIONS, FUTURE

13. Paper I

13.1 Objectives The objective of the paper I study was to refine QTLs for immune-related traits previously identified on pig chromosome 8 (SSC8), close to the mast/stem cell growth factor receptor (KIT) gene and the ms S0225, respec- tively. The study was based on 47 BC2 animals, originating from 2 wild boar x Large White intercrosses, and 191 Landrace/Large White crossbreds, born in a specific pathogen-free herd (Wattrang et al. 2005). The variation in origin of the animals used can be expected to answer to variation in immune-related genes, thus making these pigs suitable for a QTL analysis. For the present analysis, the ms SJA07, SJA10, and SJA11, located within the TLR6-1-10 gene cluster, and HIK10, residing 99 bp upstreams of TLR2, were added to the link- age map.

13.2 Results and discussion

The results reported in paper I confirm the QTLs previously identified (Wattrang et al. 2005), and points at additional regions in their vicinity, possi- bly influencing the traits studied.

As for IgM+ cells, the previously identified QTL was split into two, one lo- cated at the position of KIT and the other residing slightly upstream of this lo- cation (Fig. 13). A potential QTL for ConA induced proliferation and a QTL for PWM induced IL-2 production were both confirmed, close to S0069 and S0225, respectively. Furthermore, for PWM induced proliferation, a significant QTL was detected slightly upstream of KIT. The QTL for this trait detected close to S0225 in the previous study (Wattrang et al. 2005) was not seen in the present analysis, but the F-value was increasing at the corresponding location. Moreover, the increasing F-values at S0225 for hemoglobin and hematocrit levels, as well as for MHCII+ cells, were confirmed. Finally, a suggestive QTL for numbers of segmented neutrophils, not found in the previous study, was lo- cated close to the position of TLR2.

52 53 but a corresponding distribution of allelic effects has not yet been refuted. Studies in humans indicate that quantitative trait nucleotides (QTNs) typically lie outside genes, in regions of the genome likely to have regulatory functions. Since many regulatory regions are not conserved between species, it is unlikely that the causal nucleotides are the same. This, however, does not rule out that quantitative trait genes (QTGs) may coincide, but surprisingly few examples of common QTLs have yet been reported (Flint and Mackay 2009). As for pigs, it follows that the composition of QTLs is likely to be the same that has been found in Drosophila, mice, and – to some extent – in humans, but when it comes to QTGs and QTNs, predictions concerning similarities and differences are difficult to make.

A QTL analysis is a mere first step on the road to relate variation in pheno- types to genetic polymorphism, the ultimate goal being to pinpoint QTNs. However, by comparing the results from a QTL analysis with what is known concerning the functional properties of proteins encoded in the same chromo- somal region, it is possible to draw conclusions concerning candidate genes. From the paper I study, TLR1, TLR2, and TLR6 are among the loci calling for further attention.

14. Papers II and III Fig. 13 F-values indicating chromosome specific QTLs for immune related traits on SSC8. Markers and distances in cM on the x-axis, F-values on the y-axis. Significance levels indicated by horizon- tal lines. ConA Concanavalin A induced proliferation; PWM pokeweed mitogen induced prolifera- 14.1 Objectives tion; IL2 pwm pokeweed mitogen induced IL-2 production; IL2 pha phytohaemagglutinin induced The objective of the studies reported in papers II and III was to explore IL-2 production; Hb hemoglobin; Hem hematocrit; IgM+ IgM positive peripheral blood leucocytes; polymorphisms in TLR1, TLR2, TLR6, and TLR10 and compare their occur- MHCII+ MHCII positive peripheral blood leucocytes; Neu segmented neutrophils. rence in European wild boars and domestic pigs. These particular TLR genes, constituting the porcine TLR1 gene family, were chosen because of their local- When comparing the results reported in paper I with those from the previ- ization, close to QTLs for immune-related traits (see above), and their evolu- ous QTL analysis, the picture has become more complex. This is not surpris- tionary and functional relationships (see above). ing: more markers give more information to a QTL analysis. From the point of view of this thesis, the suggestive QTL for numbers of segmented neutrophils, Porcine TLR6, TLR1, and TLR10 are located in a cluster on the p arm of located close to the position of TLR2, is particularly interesting. In order to de- SSC8 (Shinkai et al. 2006a), while TLR2 is found on the q arm of the same termine whether this QTL really exists there are several options, among which chromosome (Jann et al. 2009). For a number of species, it has been shown that are QTL analysis using a different animal material and evaluations of polymor- the TLR6, TLR1, and TLR10 genes constitute a syntenic block (Jann et al. phisms in candidate genes (Mackay et al. 2009; Ron and Weller 2007). Owing 2009; Roach et al. 2005). Furthermore, phylogenetic analysis has shown that to its localization, HIK10, the ms representing TLR2 in the paper I study, can the extracellular regions of TLR1, TLR6, and TLR10 cluster according to not be excluded as a potential sourse of the QTL effect. gene, while the cytoplasmic regions form one cluster for TLR10 but separate clusters for respective species for TLR1 and TLR6 (Hughes and Piontkivska In Drosophila and mice, it has been found that QTLs are composed of a few 2008; Shinkai et al. 2006a). loci with moderate to large effects and a large number of loci with small effects, allelic effects following an exponential distribution. Comparisons to human QTL studies are not straightforward, due to differences in experimental design,

54 55 but a corresponding distribution of allelic effects has not yet been refuted. Studies in humans indicate that quantitative trait nucleotides (QTNs) typically lie outside genes, in regions of the genome likely to have regulatory functions. Since many regulatory regions are not conserved between species, it is unlikely that the causal nucleotides are the same. This, however, does not rule out that quantitative trait genes (QTGs) may coincide, but surprisingly few examples of common QTLs have yet been reported (Flint and Mackay 2009). As for pigs, it follows that the composition of QTLs is likely to be the same that has been found in Drosophila, mice, and – to some extent – in humans, but when it comes to QTGs and QTNs, predictions concerning similarities and differences are difficult to make.

A QTL analysis is a mere first step on the road to relate variation in pheno- types to genetic polymorphism, the ultimate goal being to pinpoint QTNs. However, by comparing the results from a QTL analysis with what is known concerning the functional properties of proteins encoded in the same chromo- somal region, it is possible to draw conclusions concerning candidate genes. From the paper I study, TLR1, TLR2, and TLR6 are among the loci calling for further attention.

14. Papers II and III Fig. 13 F-values indicating chromosome specific QTLs for immune related traits on SSC8. Markers and distances in cM on the x-axis, F-values on the y-axis. Significance levels indicated by horizon- tal lines. ConA Concanavalin A induced proliferation; PWM pokeweed mitogen induced prolifera- 14.1 Objectives tion; IL2 pwm pokeweed mitogen induced IL-2 production; IL2 pha phytohaemagglutinin induced The objective of the studies reported in papers II and III was to explore IL-2 production; Hb hemoglobin; Hem hematocrit; IgM+ IgM positive peripheral blood leucocytes; polymorphisms in TLR1, TLR2, TLR6, and TLR10 and compare their occur- MHCII+ MHCII positive peripheral blood leucocytes; Neu segmented neutrophils. rence in European wild boars and domestic pigs. These particular TLR genes, constituting the porcine TLR1 gene family, were chosen because of their local- When comparing the results reported in paper I with those from the previ- ization, close to QTLs for immune-related traits (see above), and their evolu- ous QTL analysis, the picture has become more complex. This is not surpris- tionary and functional relationships (see above). ing: more markers give more information to a QTL analysis. From the point of view of this thesis, the suggestive QTL for numbers of segmented neutrophils, Porcine TLR6, TLR1, and TLR10 are located in a cluster on the p arm of located close to the position of TLR2, is particularly interesting. In order to de- SSC8 (Shinkai et al. 2006a), while TLR2 is found on the q arm of the same termine whether this QTL really exists there are several options, among which chromosome (Jann et al. 2009). For a number of species, it has been shown that are QTL analysis using a different animal material and evaluations of polymor- the TLR6, TLR1, and TLR10 genes constitute a syntenic block (Jann et al. phisms in candidate genes (Mackay et al. 2009; Ron and Weller 2007). Owing 2009; Roach et al. 2005). Furthermore, phylogenetic analysis has shown that to its localization, HIK10, the ms representing TLR2 in the paper I study, can the extracellular regions of TLR1, TLR6, and TLR10 cluster according to not be excluded as a potential sourse of the QTL effect. gene, while the cytoplasmic regions form one cluster for TLR10 but separate clusters for respective species for TLR1 and TLR6 (Hughes and Piontkivska In Drosophila and mice, it has been found that QTLs are composed of a few 2008; Shinkai et al. 2006a). loci with moderate to large effects and a large number of loci with small effects, allelic effects following an exponential distribution. Comparisons to human QTL studies are not straightforward, due to differences in experimental design,

54 55 14.2 Results and discussion A total number of 20, 27, 26, and 33 SNPs were detected in TLR1, TLR2, TLR6, and TLR10, respectively (Fig. 14, Fig. 15). Of these, 7 (35%), 12 (44%), 15 (58%), and 20 (61%) were non-synonymous. In TLR1 and TLR2, more SNPs were present in the domestic pigs than in the wild boars. In TLR6, SNP numbers were similar in both animal groups, but the level of heterozygos- ity was higher in the domestic pigs than in the wild boars. In TLR10, again, more SNPs were present in the domestic pigs, and a higher number of non- synonymous SNPs were detected compared to the other genes.

Fig. 14 SNPs in TLR1, TLR2, and TLR6 in 25 European wild boars and 15 unrelated domestic pigs. A total number of 20, 27, and 26 SNPs were detected in TLR1, TLR2, and TLR6, respecti- vely. Of these, 7, 12, and 15 were non-synonymous. Amino acid substitutions are indicated in brackets.

56 57 14.2 Results and discussion A total number of 20, 27, 26, and 33 SNPs were detected in TLR1, TLR2, TLR6, and TLR10, respectively (Fig. 14, Fig. 15). Of these, 7 (35%), 12 (44%), 15 (58%), and 20 (61%) were non-synonymous. In TLR1 and TLR2, more SNPs were present in the domestic pigs than in the wild boars. In TLR6, SNP numbers were similar in both animal groups, but the level of heterozygos- ity was higher in the domestic pigs than in the wild boars. In TLR10, again, more SNPs were present in the domestic pigs, and a higher number of non- synonymous SNPs were detected compared to the other genes.

Fig. 14 SNPs in TLR1, TLR2, and TLR6 in 25 European wild boars and 15 unrelated domestic pigs. A total number of 20, 27, and 26 SNPs were detected in TLR1, TLR2, and TLR6, respecti- vely. Of these, 7, 12, and 15 were non-synonymous. Amino acid substitutions are indicated in brackets.

56 57 The haplotype analyses carried out for each of the four genes separately and for the TLR6-1-10 gene cluster further emphasized the differences between wild boars and domestic pigs. In TLR1, TLR2, and TLR6, both animal groups had group specific haplotypes: 6, 5, and 14 haplotypes were found only in wild boars or domestic pigs, respectively. Furthermore, in the wild boars, one major haplotype, present at a high frequency, was found in all three genes, while the same pattern was exhibited only by TLR2 in the domestic pigs. In TLR10, three of five resolved haplotypes were specific for the domestic pigs. For the TLR6-1-10 gene cluster, 13 haplotypes were resolved. In the wild boars, six haplotypes, out of which four were group-specific, were present, while nine haplotypes, out of which seven were group-specific, were represented in the domestic pigs.

Owing to the domestication process and differences in selective pressure, differences in polymorphic patterns between wild boars and domestic pigs are not surprising. However, since breeding means choosing among genotypes, the opposite pattern – more SNPs in wild boars than in domestic pigs – would have been expected. This expected pattern has been found regarding the MHC complex in commercial chicken breeds (Schou et al. 2007). However, analyses of the mitochondrial as well as the nuclear genome indicate that the European wild boar population went through a bottle neck before domestication started (Larson et al. 2005; Megens et al. 2008), resulting in a low level of diversity. Furthermore, the pattern described in papers II and III is a confirmation of what previously has been found concerning the IGF2 gene (Ojeda et al. 2008). The higher level of diversity seen in the domestic pigs may be caused by relaxed purifying selection during the domestication process.

In paper II, it is further reported that the mean heterozygosity at polymor- phic positions was found to be markedly lower in TLR2 than in TLR1 and TLR6. This might be due to functional differences within the TLR receptor complexes. Using mouse models, it has been shown that TLR1 and TLR6 confer the fine-tuning of the MAMP discriminating capability to the TLR2/1 and TLR2/6 heterodimers (Takeuchi et al. 2001; Takeuchi et al. 2002). Thus, Fig. 15 SNPs in TLR10 in 15 European wild boars and 15 unrelated domestic pigs. A total num- a high level of variation in porcine TLR1 and TLR6 might add more to the ef- ber of 33 SNPs were detected in TLR10. Of these, 20 were non-synonymous. Amino acid substitu- ficiency of the dimer’s antigen recognition, while TLR2 may not need – or may tions are indicated in brackets. be harmed by – the same level of variability. Moreover, homology modeling of human TLR6 based on the TLR1 crystal structure suggests that the structure of TLR6 resembles that of TLR1, except that a narrow hydrophobic channel, constituting part of the ligand binding site in TLR1, is absent (Jin et al. 2007). A similar difference between TLR1 and TLR6 in mice might explain the dif- ferences found (Takeuchi et al. 2001; Takeuchi et al. 2002) regarding their ligand specificity.

58 59 The haplotype analyses carried out for each of the four genes separately and for the TLR6-1-10 gene cluster further emphasized the differences between wild boars and domestic pigs. In TLR1, TLR2, and TLR6, both animal groups had group specific haplotypes: 6, 5, and 14 haplotypes were found only in wild boars or domestic pigs, respectively. Furthermore, in the wild boars, one major haplotype, present at a high frequency, was found in all three genes, while the same pattern was exhibited only by TLR2 in the domestic pigs. In TLR10, three of five resolved haplotypes were specific for the domestic pigs. For the TLR6-1-10 gene cluster, 13 haplotypes were resolved. In the wild boars, six haplotypes, out of which four were group-specific, were present, while nine haplotypes, out of which seven were group-specific, were represented in the domestic pigs.

Owing to the domestication process and differences in selective pressure, differences in polymorphic patterns between wild boars and domestic pigs are not surprising. However, since breeding means choosing among genotypes, the opposite pattern – more SNPs in wild boars than in domestic pigs – would have been expected. This expected pattern has been found regarding the MHC complex in commercial chicken breeds (Schou et al. 2007). However, analyses of the mitochondrial as well as the nuclear genome indicate that the European wild boar population went through a bottle neck before domestication started (Larson et al. 2005; Megens et al. 2008), resulting in a low level of diversity. Furthermore, the pattern described in papers II and III is a confirmation of what previously has been found concerning the IGF2 gene (Ojeda et al. 2008). The higher level of diversity seen in the domestic pigs may be caused by relaxed purifying selection during the domestication process.

In paper II, it is further reported that the mean heterozygosity at polymor- phic positions was found to be markedly lower in TLR2 than in TLR1 and TLR6. This might be due to functional differences within the TLR receptor complexes. Using mouse models, it has been shown that TLR1 and TLR6 confer the fine-tuning of the MAMP discriminating capability to the TLR2/1 and TLR2/6 heterodimers (Takeuchi et al. 2001; Takeuchi et al. 2002). Thus, Fig. 15 SNPs in TLR10 in 15 European wild boars and 15 unrelated domestic pigs. A total num- a high level of variation in porcine TLR1 and TLR6 might add more to the ef- ber of 33 SNPs were detected in TLR10. Of these, 20 were non-synonymous. Amino acid substitu- ficiency of the dimer’s antigen recognition, while TLR2 may not need – or may tions are indicated in brackets. be harmed by – the same level of variability. Moreover, homology modeling of human TLR6 based on the TLR1 crystal structure suggests that the structure of TLR6 resembles that of TLR1, except that a narrow hydrophobic channel, constituting part of the ligand binding site in TLR1, is absent (Jin et al. 2007). A similar difference between TLR1 and TLR6 in mice might explain the dif- ferences found (Takeuchi et al. 2001; Takeuchi et al. 2002) regarding their ligand specificity.

58 59 In vitro experiments suggest that over-expression of human TLRs may 15. Paper IV cause constitutive activation of NFκB, leading to excessive inflammation and tissue damage. All ten human TLRs have been identified based on sequence features, as opposed to identification on the protein level, based on their ability 15.1 Objective to recognize MAMPs (Zhong et al. 2005). Moreover, it has been shown that The study reported in paper IV had two objectives. The first was to analyze surface expression of human TLR1 and TLR4 in monocytes and immature MBL-A serum concentrations in 45 wild boar/Large White crossbred pigs be- DCs is very low (Visintin et al. 2001). Together, this indicates that TLR sig- fore and after immunization with an E-coli vaccine. These pigs constitute the nalling is powerful and needs to be subjected to tight regulation. Protein ex- BC2 generation in a family material originating from two wild boar x Large pression is regulated at transcriptional as well as translational level and through White intercrosses. The second objective was to study genotypes in the MBL1 post-translational modifications. However, codon usage provides an alternative gene in this family material and in a group of wild boars from the south of mechanism for regulation. It has been shown that in most human TLR genes, Sweden. MBL is interesting owing to its disputed biological significance in the codon usage pattern deviates from the general one in humans (Zhong et al. humans and its suggested future role as a breeding parameter in pigs (Juul- 2005). In papers II and III, a similar difference between the general codon us- Madsen 2006; Verdu 2006). age pattern in pigs and that in TLR1, TLR2, TLR6, and TLR10 was detected. It is likely that this specific TLR codon usage pattern – mirroring the condi- Pigs, as most mammals, have two forms of MBL, MBL-A and MBL-C. tions in humans – results in protein expression at a low level, which could be These are encoded by two separate genes, denoted MBL1 and MBL2, respec- interpreted as a safety measure, in order to avoid tissue damage. tively. Porcine MBL1 is located on chromosome 14, in a cluster known as the collectin locus, together with Sp-A and Sp-D (Lillie et al. 2006a). MBL2 has Concerning human TLRs, it is believed that those located on the cell sur- also been mapped to chromosome 14, and a QTL for C3c (a marker for com- face, sensing bacterial products, are more redundant than the nucleic acid sen- plement activation in vivo) serum concentrations has been identified close to sors in intracellular compartments. In a large-scale sequencing project involving the position of the MBL genes (Phatsara et al. 2007). all 10 human TLR genes in African, European, and East-Asian populations, TLR10 was found to be by far the most diverse (Barreiro et al. 2009). As to pigs, it has recently been reported that non-synonymous SNPs are more rare in TLR3, TLR7, and TLR8 than in cell surface TLR genes (Morozumi and Uen- 15.2 Results and discussion ishi 2009; Shinkai et al. 2006b). This suggests a pattern similar to that in hu- MBL-A concentrations in serum in the BC2 animals showed a bimodal dis- mans. Moreover, non-synonymous SNPs in the coding regions of cattle TLR3, tribution (Fig. 16). Fifteen pigs, i.e. one third of the animals sampled, were TLR7, and TLR8 also seem to be few (Cargill and Womack 2007). Thus, it registered with serum concentrations between 0.7 and 1.6 ug/ml. In this group, seems likely that non-synonymous SNPs are tolerated less often in TLR genes concentrations were not affected by the vaccination. The remaining 30 pigs encoding nucleic acid sensors than in those encoding cell surface TLRs. were registered with MBL-A concentrations between 4.5 and 52.0 ug/ml, with Whether originating from microorganisms or the host, there is little change in a significant increase in mean level after vaccination. the structure of nucleic acids over time; consequently, there is no need for posi- tive selection to adapt to new MAMP structures. There is, however, need to protect function through purifying selection. In the studies reported in papers II and III, the highest number of non-synonymous SNPs was detected in TLR10. Possibly, this suggests redundancy for TLR10 in pigs.

60 61 In vitro experiments suggest that over-expression of human TLRs may 15. Paper IV cause constitutive activation of NFκB, leading to excessive inflammation and tissue damage. All ten human TLRs have been identified based on sequence features, as opposed to identification on the protein level, based on their ability 15.1 Objective to recognize MAMPs (Zhong et al. 2005). Moreover, it has been shown that The study reported in paper IV had two objectives. The first was to analyze surface expression of human TLR1 and TLR4 in monocytes and immature MBL-A serum concentrations in 45 wild boar/Large White crossbred pigs be- DCs is very low (Visintin et al. 2001). Together, this indicates that TLR sig- fore and after immunization with an E-coli vaccine. These pigs constitute the nalling is powerful and needs to be subjected to tight regulation. Protein ex- BC2 generation in a family material originating from two wild boar x Large pression is regulated at transcriptional as well as translational level and through White intercrosses. The second objective was to study genotypes in the MBL1 post-translational modifications. However, codon usage provides an alternative gene in this family material and in a group of wild boars from the south of mechanism for regulation. It has been shown that in most human TLR genes, Sweden. MBL is interesting owing to its disputed biological significance in the codon usage pattern deviates from the general one in humans (Zhong et al. humans and its suggested future role as a breeding parameter in pigs (Juul- 2005). In papers II and III, a similar difference between the general codon us- Madsen 2006; Verdu 2006). age pattern in pigs and that in TLR1, TLR2, TLR6, and TLR10 was detected. It is likely that this specific TLR codon usage pattern – mirroring the condi- Pigs, as most mammals, have two forms of MBL, MBL-A and MBL-C. tions in humans – results in protein expression at a low level, which could be These are encoded by two separate genes, denoted MBL1 and MBL2, respec- interpreted as a safety measure, in order to avoid tissue damage. tively. Porcine MBL1 is located on chromosome 14, in a cluster known as the collectin locus, together with Sp-A and Sp-D (Lillie et al. 2006a). MBL2 has Concerning human TLRs, it is believed that those located on the cell sur- also been mapped to chromosome 14, and a QTL for C3c (a marker for com- face, sensing bacterial products, are more redundant than the nucleic acid sen- plement activation in vivo) serum concentrations has been identified close to sors in intracellular compartments. In a large-scale sequencing project involving the position of the MBL genes (Phatsara et al. 2007). all 10 human TLR genes in African, European, and East-Asian populations, TLR10 was found to be by far the most diverse (Barreiro et al. 2009). As to pigs, it has recently been reported that non-synonymous SNPs are more rare in TLR3, TLR7, and TLR8 than in cell surface TLR genes (Morozumi and Uen- 15.2 Results and discussion ishi 2009; Shinkai et al. 2006b). This suggests a pattern similar to that in hu- MBL-A concentrations in serum in the BC2 animals showed a bimodal dis- mans. Moreover, non-synonymous SNPs in the coding regions of cattle TLR3, tribution (Fig. 16). Fifteen pigs, i.e. one third of the animals sampled, were TLR7, and TLR8 also seem to be few (Cargill and Womack 2007). Thus, it registered with serum concentrations between 0.7 and 1.6 ug/ml. In this group, seems likely that non-synonymous SNPs are tolerated less often in TLR genes concentrations were not affected by the vaccination. The remaining 30 pigs encoding nucleic acid sensors than in those encoding cell surface TLRs. were registered with MBL-A concentrations between 4.5 and 52.0 ug/ml, with Whether originating from microorganisms or the host, there is little change in a significant increase in mean level after vaccination. the structure of nucleic acids over time; consequently, there is no need for posi- tive selection to adapt to new MAMP structures. There is, however, need to protect function through purifying selection. In the studies reported in papers II and III, the highest number of non-synonymous SNPs was detected in TLR10. Possibly, this suggests redundancy for TLR10 in pigs.

60 61 2010). The data for the smaller group of animals presented in paper IV are in agreement with this, but the pigs in the paper IV study were clinically healthy, and no difference in total and differential leukocyte counts or hemoglobin and hematocrit levels were found between the two groups of animals.

Sequencing of part of the MBL1 gene revealed a total of seven SNPs in ex- ons 1 and 2 (Fig. 17). One of these – G949T – has previously been detected in domestic pigs and is assumed to affect MBL-A concentration in serum (Lillie et al. 2006b; Juul-Madsen et al 2006). The study reported in paper IV shows that this SNP is also present in wild boars. Furthermore, a novel SNP with po- tential to be functionally important – G131A – was also detected.

Fig. 16 MBL-A concentration in serum in wild boar/Large White crossbreds. Fifteen samples, cor- responding to one third of the animals, had MBL-A concentrations between 0.7 and 1.6 ug/ml, be- fore and after vaccination alike. The other 30 samples had concentrations between 4.5 and 52.0 ug/ml, with a significant increase in mean level after vaccination (13.1±8.3 ug/ml before and 16.4±6.2 ug/ml after vaccination, P<0.05).

Since MBL is an acute phase reactant (Thiel et al. 1992) and the second sampling was carried out three weeks after the first, it is not possible to draw any unambiguous conclusions concerning the impact of the vaccination. How- ever, it is clear that the higher mean MBL-A level after vaccination is not asso- ciated with age. Recently, a time study, monitoring several immune parameters, has shown that at around the time of weaning, MBL-A concentrations rise in pigs expressing normal amounts of the protein at birth. However, this elevated level of MBL-A is temporary: at 12 weeks of age, when the first sampling of the BC2 animals was carried out, MBL-A concentration in serum is lower than at birth and does not change up to 15 weeks of age (Juul-Madsen et al. 2010). Thus, it is tempting to speculate that our data show a connection between vac- cination and MBL-A levels, such that the presence of antigens elicits aug- mented levels of MBL-A expression. In a diseased pig, registered with low Fig. 17 SNPs in MBL1 in eight European wild boars, five Large White pigs, and nine wild MBL-A concentration at birth and excluded from the statistical analysis in the boar/Large White crossbreds. Seven SNPs were detected in exons 1 and 2. Of these, five were non- time study, concentrations barely changed up to week 15 (Juul-Madsen et al. synonymous. Amino acid substitutions are indicated in brackets.

62 63 2010). The data for the smaller group of animals presented in paper IV are in agreement with this, but the pigs in the paper IV study were clinically healthy, and no difference in total and differential leukocyte counts or hemoglobin and hematocrit levels were found between the two groups of animals.

Sequencing of part of the MBL1 gene revealed a total of seven SNPs in ex- ons 1 and 2 (Fig. 17). One of these – G949T – has previously been detected in domestic pigs and is assumed to affect MBL-A concentration in serum (Lillie et al. 2006b; Juul-Madsen et al 2006). The study reported in paper IV shows that this SNP is also present in wild boars. Furthermore, a novel SNP with po- tential to be functionally important – G131A – was also detected.

Fig. 16 MBL-A concentration in serum in wild boar/Large White crossbreds. Fifteen samples, cor- responding to one third of the animals, had MBL-A concentrations between 0.7 and 1.6 ug/ml, be- fore and after vaccination alike. The other 30 samples had concentrations between 4.5 and 52.0 ug/ml, with a significant increase in mean level after vaccination (13.1±8.3 ug/ml before and 16.4±6.2 ug/ml after vaccination, P<0.05).

Since MBL is an acute phase reactant (Thiel et al. 1992) and the second sampling was carried out three weeks after the first, it is not possible to draw any unambiguous conclusions concerning the impact of the vaccination. How- ever, it is clear that the higher mean MBL-A level after vaccination is not asso- ciated with age. Recently, a time study, monitoring several immune parameters, has shown that at around the time of weaning, MBL-A concentrations rise in pigs expressing normal amounts of the protein at birth. However, this elevated level of MBL-A is temporary: at 12 weeks of age, when the first sampling of the BC2 animals was carried out, MBL-A concentration in serum is lower than at birth and does not change up to 15 weeks of age (Juul-Madsen et al. 2010). Thus, it is tempting to speculate that our data show a connection between vac- cination and MBL-A levels, such that the presence of antigens elicits aug- mented levels of MBL-A expression. In a diseased pig, registered with low Fig. 17 SNPs in MBL1 in eight European wild boars, five Large White pigs, and nine wild MBL-A concentration at birth and excluded from the statistical analysis in the boar/Large White crossbreds. Seven SNPs were detected in exons 1 and 2. Of these, five were non- time study, concentrations barely changed up to week 15 (Juul-Madsen et al. synonymous. Amino acid substitutions are indicated in brackets.

62 63 The basic subunit of the MBL molecule is a trimer composed of three iden- play different polymorphic patterns in the TLR1 gene family, the domestic tical monomers, each consisting of four distinct domains: a cystein-rich N- pigs being more diverse. A comparison of polymorphic patterns in TLR3, terminal domain, a collagen-like domain, a neck region, and a C-terminal car- TLR7, and TLR8 in wild boars and domestic pigs would reveal whether this is bohydrate recognition domain. The collagen-like domain consists of numerous the case in these genes as well. If so, the wild boar genotypes are likely to be G-X-Y motifs (X and Y representing any amino acid), leading to the formation close to fixation. of a triple helix structure (Fujita et al. 2004). In the coding sequence of human MBL2, three non-synonymous SNPs, all with decreasing effects on MBL con- Concerning MBL-A, it has been shown that serum concentrations and centrations, have been identified. Two of these implicate the replacement of a heritability estimates differ between Landrace and Duroc pigs, and the usage of glycine residue in a G-X-Y motif with a more bulky amino acid (Garred et al. MBL-A as a breeding parameter has been suggested (Juul-Madsen et al. 2006). 2006). As for pigs, the T allele at position 949 implicates the substitution of the The MBL study reported in paper IV confirmed the presence of a previously glycine residue in G-X-Y motif 16 for a cysteine (Lillie et al. 2006b). Similarly, reported SNP (Lillie et al. 2006b), assumed to affect MBL-A concentrations, the A allele at position 131 substitutes the glycine residue in G-X-Y motif 1 for in purebred Large White pigs. Moreover, a novel SNP, with potential to be an aspartic acid. The variant human MBL alleles are thought to disrupt the tri- functionally important, was detected in the same breed. A comparison between ple helix formation and cause dysfunctional MBL molecules (Garred et al. MBL-A concentrations and heritability estimates in Large White pigs on the 2003; Garred et al. 2006). It is likely that the T allele of G949T and the A al- one hand and Landrace and Duroc pigs on the other might bring forth more lele of G131A interfere with porcine MBL-A triple helix formation in a similar relevant information concerning the usage of MBL-A as a breeding parameter. manner. In two animals registered with low MBL-A concentrations in the pa- Moreover, an analysis focusing on the relationship between the novel SNP and per IV study, the T allele was present at position G949T, probably causing a MBL-A concentrations in serum would reveal to what degree serum concen- reduction of functional MBL-A molecules. trations are affected by this particular SNP.

16. Conclusions

From the work described in papers I-IV, it can be concluded that

 the QTLs for immune-related traits identified on SSC8 are more complex than previously anticipated  European wild boars and domestic pigs display different polymorphic patterns in the TLR1, TLR2, TLR6, and TLR10 genes, the domestic pigs being more diverse than the wild boars,  the TLR10 gene is more diverse than the TLR1, TLR2, and TLR6 genes and may be redundant in pigs,  the SNP G949T, previously detected in domestic pigs and assumed to affect MBL-A concentration in serum, is also present in European wild boars,  the novel SNP G131A is likely to be functionally important.

17. Future It has recently been shown that non-synonymous SNPs are more rare in TLR3, TLR7, and TLR8 than in cell surface TLR genes in various pig breeds (Morozumi and Uenishi 2009; Shinkai et al. 2006b). The studies reported in papers II and III have shown that European wild boars and domestic pigs dis-

64 65 The basic subunit of the MBL molecule is a trimer composed of three iden- play different polymorphic patterns in the TLR1 gene family, the domestic tical monomers, each consisting of four distinct domains: a cystein-rich N- pigs being more diverse. A comparison of polymorphic patterns in TLR3, terminal domain, a collagen-like domain, a neck region, and a C-terminal car- TLR7, and TLR8 in wild boars and domestic pigs would reveal whether this is bohydrate recognition domain. The collagen-like domain consists of numerous the case in these genes as well. If so, the wild boar genotypes are likely to be G-X-Y motifs (X and Y representing any amino acid), leading to the formation close to fixation. of a triple helix structure (Fujita et al. 2004). In the coding sequence of human MBL2, three non-synonymous SNPs, all with decreasing effects on MBL con- Concerning MBL-A, it has been shown that serum concentrations and centrations, have been identified. Two of these implicate the replacement of a heritability estimates differ between Landrace and Duroc pigs, and the usage of glycine residue in a G-X-Y motif with a more bulky amino acid (Garred et al. MBL-A as a breeding parameter has been suggested (Juul-Madsen et al. 2006). 2006). As for pigs, the T allele at position 949 implicates the substitution of the The MBL study reported in paper IV confirmed the presence of a previously glycine residue in G-X-Y motif 16 for a cysteine (Lillie et al. 2006b). Similarly, reported SNP (Lillie et al. 2006b), assumed to affect MBL-A concentrations, the A allele at position 131 substitutes the glycine residue in G-X-Y motif 1 for in purebred Large White pigs. Moreover, a novel SNP, with potential to be an aspartic acid. The variant human MBL alleles are thought to disrupt the tri- functionally important, was detected in the same breed. A comparison between ple helix formation and cause dysfunctional MBL molecules (Garred et al. MBL-A concentrations and heritability estimates in Large White pigs on the 2003; Garred et al. 2006). It is likely that the T allele of G949T and the A al- one hand and Landrace and Duroc pigs on the other might bring forth more lele of G131A interfere with porcine MBL-A triple helix formation in a similar relevant information concerning the usage of MBL-A as a breeding parameter. manner. In two animals registered with low MBL-A concentrations in the pa- Moreover, an analysis focusing on the relationship between the novel SNP and per IV study, the T allele was present at position G949T, probably causing a MBL-A concentrations in serum would reveal to what degree serum concen- reduction of functional MBL-A molecules. trations are affected by this particular SNP.

16. Conclusions

From the work described in papers I-IV, it can be concluded that

 the QTLs for immune-related traits identified on SSC8 are more complex than previously anticipated  European wild boars and domestic pigs display different polymorphic patterns in the TLR1, TLR2, TLR6, and TLR10 genes, the domestic pigs being more diverse than the wild boars,  the TLR10 gene is more diverse than the TLR1, TLR2, and TLR6 genes and may be redundant in pigs,  the SNP G949T, previously detected in domestic pigs and assumed to affect MBL-A concentration in serum, is also present in European wild boars,  the novel SNP G131A is likely to be functionally important.

17. Future It has recently been shown that non-synonymous SNPs are more rare in TLR3, TLR7, and TLR8 than in cell surface TLR genes in various pig breeds (Morozumi and Uenishi 2009; Shinkai et al. 2006b). The studies reported in papers II and III have shown that European wild boars and domestic pigs dis-

64 65

ACKNOWLEDGEMENTS

This thesis was brewed (!), together with large quantities of green (occa-

sionally red) tea, during January and February 2010 in my “cloister cell” at

Norrgård and at home. Martin, you were even more patient than usual during

this period, and I’m forever grateful.

Obviously, many people have contributed to making me who I am today.

Most of you will know who you are, and if you don’t, our Lord certainly does,

so long lists of names are not necessary. However, there are a few people I

would like to mention in particular.

First, Assoc. prof. Inger Edfors, my main supervisor, who has been very

brave, coping with all my phases, and without whom my work had been impos-

sible,

Prof. Kristina Nilsson Ekdahl, my co-supervisor, who has been my never

failing safety net, tremendously supportive and kind,

PhD Kjell Edman, my second co-supervisor, who is the man who speaks

with computers,

Prof. Sven Tågerud, who is loyal to PhD students,

MSc Marlene Norrby, my fellow PhD student and the best room mate ever,

who has been providing “therapy” when necessary and who always takes enter-

tainment value into account,

MSc Maria Gullberg, MSc Stina Israelsson, and MSc Nina Jonsson, my lab

mates, who possess a lot of girl power and have supported me without end,

Seniorforsker Helle R. Juul-Madsen and her group at Foulum, who made

the guest from Sweden feel very welcome,

and, most importantly, my husband Martin. Your steadfast love means eve-

rything to me.

Finally …

Filip, 5 år: ”Vad forskar Ingrid-Maria om?” Mamma Nina: ”Hon tittar på skillnaden mellan vildsvin och vanliga grisar.” Filip, 5 år: ”Det kan inte ta så lång tid. Vildsvin är svarta och bor i skogen och har betar och grisar bor i lagården.”

S.D.G.

66 67

ACKNOWLEDGEMENTS

This thesis was brewed (!), together with large quantities of green (occa-

sionally red) tea, during January and February 2010 in my “cloister cell” at

Norrgård and at home. Martin, you were even more patient than usual during

this period, and I’m forever grateful.

Obviously, many people have contributed to making me who I am today.

Most of you will know who you are, and if you don’t, our Lord certainly does,

so long lists of names are not necessary. However, there are a few people I

would like to mention in particular.

First, Assoc. prof. Inger Edfors, my main supervisor, who has been very

brave, coping with all my phases, and without whom my work had been impos-

sible,

Prof. Kristina Nilsson Ekdahl, my co-supervisor, who has been my never

failing safety net, tremendously supportive and kind,

PhD Kjell Edman, my second co-supervisor, who is the man who speaks

with computers,

Prof. Sven Tågerud, who is loyal to PhD students,

MSc Marlene Norrby, my fellow PhD student and the best room mate ever,

who has been providing “therapy” when necessary and who always takes enter-

tainment value into account,

MSc Maria Gullberg, MSc Stina Israelsson, and MSc Nina Jonsson, my lab

mates, who possess a lot of girl power and have supported me without end,

Seniorforsker Helle R. Juul-Madsen and her group at Foulum, who made

the guest from Sweden feel very welcome,

and, most importantly, my husband Martin. Your steadfast love means eve-

rything to me.

Finally …

Filip, 5 år: ”Vad forskar Ingrid-Maria om?” Mamma Nina: ”Hon tittar på skillnaden mellan vildsvin och vanliga grisar.” Filip, 5 år: ”Det kan inte ta så lång tid. Vildsvin är svarta och bor i skogen och har betar och grisar bor i lagården.”

S.D.G.

66 67

REFERENCES

Abbas AK, Janeway CA Jr (2000) Immunology: improving on nature in the

twenty-first century. Cell 100:129-138.

Abreu MT (2010) Toll-like receptor signalling in the intestinal epithelium:

how bacterial recognition shapes intestinal function. Nat Rev Immunol 10:131-

144.

Akashi-Takamura S, Miyake K (2008) TLR accessory molecules. Curr Opin

Immunol 20:420-425.

Akira S, Uematsu S, Takeuchi O (2006) Pathogen recognition and innate im-

munity. Cell 124:783-801.

Anderson KV, Bokla L, Nüsslein-Volhard C (1985) Establishment of dorsal-

ventral polarity in the drosophila embryo: The induction of polarity by the Toll

gene product. Cell 42:791-798.

Barreiro LB, Ben-Ali M, Quach H, Laval G, Patin E, Pickrell JK, Bouchier C,

Tichit M, Neyrolles O, Gicquel B, Kidd JR, Kidd KK, Alcaïs A, Ragimbeau J,

Pellegrini S, Abel L, Casanova JL, Quintana-Murci L (2009) Evolutionary dy-

namics of human Toll-like receptors and their different contributions to host

defense. PLoS Genet 5(7):e1000562.

Barreiro LB, Quintana-Murci L (2010) From evolutionary genetics to human

immunology: how selection shapes host defence genes. Nat Rev Genet 11:17-

30.

Bernig T, Taylor JG, Foster CB, Staats B, Yeager M, Chanock SJ (2004) Se-

quence analysis of the mannose-binding lectin (MBL2) gene reveals a high de-

gree of heterozygosity with evidence of selection. Genes Immun 5:461-476.

Beutler B, Rehli M (2002) Evolution of the TIR, Tolls and TLRs: functional

inferences from computational biology. Curr Top Microbiol Immunol 270:1-

21.

Bilak H, Tauszig-Delamasure S, Imler JL (2003) Toll and Toll-like receptors

in Drosophila. Biochem Soc Trans 31(Pt 3):648-651.

68 69

REFERENCES

Abbas AK, Janeway CA Jr (2000) Immunology: improving on nature in the

twenty-first century. Cell 100:129-138.

Abreu MT (2010) Toll-like receptor signalling in the intestinal epithelium:

how bacterial recognition shapes intestinal function. Nat Rev Immunol 10:131-

144.

Akashi-Takamura S, Miyake K (2008) TLR accessory molecules. Curr Opin

Immunol 20:420-425.

Akira S, Uematsu S, Takeuchi O (2006) Pathogen recognition and innate im-

munity. Cell 124:783-801.

Anderson KV, Bokla L, Nüsslein-Volhard C (1985) Establishment of dorsal-

ventral polarity in the drosophila embryo: The induction of polarity by the Toll

gene product. Cell 42:791-798.

Barreiro LB, Ben-Ali M, Quach H, Laval G, Patin E, Pickrell JK, Bouchier C,

Tichit M, Neyrolles O, Gicquel B, Kidd JR, Kidd KK, Alcaïs A, Ragimbeau J,

Pellegrini S, Abel L, Casanova JL, Quintana-Murci L (2009) Evolutionary dy-

namics of human Toll-like receptors and their different contributions to host

defense. PLoS Genet 5(7):e1000562.

Barreiro LB, Quintana-Murci L (2010) From evolutionary genetics to human

immunology: how selection shapes host defence genes. Nat Rev Genet 11:17-

30.

Bernig T, Taylor JG, Foster CB, Staats B, Yeager M, Chanock SJ (2004) Se-

quence analysis of the mannose-binding lectin (MBL2) gene reveals a high de-

gree of heterozygosity with evidence of selection. Genes Immun 5:461-476.

Beutler B, Rehli M (2002) Evolution of the TIR, Tolls and TLRs: functional

inferences from computational biology. Curr Top Microbiol Immunol 270:1-

21.

Bilak H, Tauszig-Delamasure S, Imler JL (2003) Toll and Toll-like receptors

in Drosophila. Biochem Soc Trans 31(Pt 3):648-651.

68 69 Boyd A, Philbin VJ, Smith AL (2007) Conserved and distinct aspects of the Friec GL, Kemper C (2009) Complement: coming full circle. Arch Immunol avian Toll-like receptor (TLR) system: implications for transmission and con- Ther Exp (Warsz) 57:393-407. trol of bird-borne zoonoses. Biochem Soc Trans 35(Pt 6):1504-1507. Fujita T, Matsushita M, Endo Y (2004) The lectin-complement pathway-its Brikos C, O’Neill AJO (2008) Signalling of Toll-like receptors. in: Bauer S, role in innate immunity and evolution. Immunol Rev 198:185-202. Hartmann G (eds): Toll-like receptors (TLRs) and innate immunity. Handb Exp Pharmacol vol 183, pp 21-50, Springer-Verlag, Berlin Heidelberg, Ger- Garlanda C, Maina V, Cotena A, Moallia F (2009) The soluble pattern recog- many nition receptor pentraxin-3 in innate immunity, inflammation and fertility. J Reprod Immunol 83:128-133. Burch-Smith TM, Dinesh-Kumar SP (2007) The functions of plant TIR do- mains. Sci STKE 2007(401):pe46. Garred P, Larsen F, Madsen HO, Koch C (2003) Mannose-binding lectin de- ficiency-revisited. Mol Immunol 40:73-84. Cargill EJ, Womack JE (2007) Detection of polymorphisms in bovine toll-like receptors 3, 7, 8, and 9. Genomics 89:745-755. Garred P, Larsen F, Seyfarth J, Fujita R, Madsen HO (2006) Mannose- binding lectin and its genetic variants. Genes Immun 7:85-94. Cario E, Gerken G, Podolsky DK (2007) Toll-like receptor 2 controls mucosal inflammation by regulating epithelial barrier function. Gastroenterology Garred P (2008) Mannose-binding lectin genetics: from A to Z. Biochem Soc 132:1359-1374. Trans 36(Pt 6):1461-1466.

Darwin C (2000) On the origin of species, a facsimile of the first edition, Har- Garred P, Honoré C, Ma YJ, Munthe-Fog L, Hummelshøj T (2009) MBL2, vard University Press, Cambridge (MA), USA FCN1, FCN2 and FCN3 - the genes behind the initiation of the lectin path- way of complement. Mol Immunol 46:2737-2744. Deban L, Bottazzi B, Garlanda C, de la Torre YM, Mantovani A (2009) Pentraxins: multifunctional proteins at the interface of innate immunity and in- Gay NJ, Gangloff M (2007) Structure and function of Toll receptors and their flammation. Biofactors 35:138-145. ligands. Annu Rev Biochem 76:141-165. de Wit, PJ (2007) How plants recognize pathogens and defend themselves. Giuffra E, Kijas JM, Amarger V, Carlborg Ö, Jeon JT, Andersson L (2000) Cell Mol Life Sci 64:2726-2732. The origin of the domestic pig: independent domestication and subsequent in- trogression. Genetics 154:1785-1791. Ellegren H (2004) Microsatellites: simple sequences with complex evolution. Nat Rev Genet 5:435-445. Graur D, Li WH (2000) Fundamentals of molecular evolution, 2 ed, Sinauer Associates, Sunderland (MA), USA Falconer DS, Mackay TFC (1996) Introduction to quantitative genetics, 4 ed, Longman, Harlow, Essex, UK Graves HB (1984) Behavior and Ecology of Wild and Feral Swine (Sus Scrofa). J Anim Sci 1984 58:482-492. Fang M, Larson G, Ribeiro HS, Li N, Andersson L (2009) Contrasting mode of evolution at a coat color locus in wild and domestic pigs. PLoS Genetics Green P, Falls K, Crooks S (1990) Documentation for CRI-MAP, version 2.4. 5:1-6. Washington University School of Medicine St. Louis (MO), USA

Flajnik MF, Kasahara M (2010) Origin and evolution of the adaptive immune Hajishengallis G, Lambris JD (2010) Crosstalk pathways between Toll-like re- system: genetic events and selective pressures. Nat Rev Genet 11:47-59. ceptors and the complement system. Trends Immunol Feb 10. (Epub ahead of print) Flint J, Mackay TF (2009) Genetic architecture of quantitative traits in mice, flies, and humans. Genome Res 19:723-733.

70 71 Boyd A, Philbin VJ, Smith AL (2007) Conserved and distinct aspects of the Friec GL, Kemper C (2009) Complement: coming full circle. Arch Immunol avian Toll-like receptor (TLR) system: implications for transmission and con- Ther Exp (Warsz) 57:393-407. trol of bird-borne zoonoses. Biochem Soc Trans 35(Pt 6):1504-1507. Fujita T, Matsushita M, Endo Y (2004) The lectin-complement pathway-its Brikos C, O’Neill AJO (2008) Signalling of Toll-like receptors. in: Bauer S, role in innate immunity and evolution. Immunol Rev 198:185-202. Hartmann G (eds): Toll-like receptors (TLRs) and innate immunity. Handb Exp Pharmacol vol 183, pp 21-50, Springer-Verlag, Berlin Heidelberg, Ger- Garlanda C, Maina V, Cotena A, Moallia F (2009) The soluble pattern recog- many nition receptor pentraxin-3 in innate immunity, inflammation and fertility. J Reprod Immunol 83:128-133. Burch-Smith TM, Dinesh-Kumar SP (2007) The functions of plant TIR do- mains. Sci STKE 2007(401):pe46. Garred P, Larsen F, Madsen HO, Koch C (2003) Mannose-binding lectin de- ficiency-revisited. Mol Immunol 40:73-84. Cargill EJ, Womack JE (2007) Detection of polymorphisms in bovine toll-like receptors 3, 7, 8, and 9. Genomics 89:745-755. Garred P, Larsen F, Seyfarth J, Fujita R, Madsen HO (2006) Mannose- binding lectin and its genetic variants. Genes Immun 7:85-94. Cario E, Gerken G, Podolsky DK (2007) Toll-like receptor 2 controls mucosal inflammation by regulating epithelial barrier function. Gastroenterology Garred P (2008) Mannose-binding lectin genetics: from A to Z. Biochem Soc 132:1359-1374. Trans 36(Pt 6):1461-1466.

Darwin C (2000) On the origin of species, a facsimile of the first edition, Har- Garred P, Honoré C, Ma YJ, Munthe-Fog L, Hummelshøj T (2009) MBL2, vard University Press, Cambridge (MA), USA FCN1, FCN2 and FCN3 - the genes behind the initiation of the lectin path- way of complement. Mol Immunol 46:2737-2744. Deban L, Bottazzi B, Garlanda C, de la Torre YM, Mantovani A (2009) Pentraxins: multifunctional proteins at the interface of innate immunity and in- Gay NJ, Gangloff M (2007) Structure and function of Toll receptors and their flammation. Biofactors 35:138-145. ligands. Annu Rev Biochem 76:141-165. de Wit, PJ (2007) How plants recognize pathogens and defend themselves. Giuffra E, Kijas JM, Amarger V, Carlborg Ö, Jeon JT, Andersson L (2000) Cell Mol Life Sci 64:2726-2732. The origin of the domestic pig: independent domestication and subsequent in- trogression. Genetics 154:1785-1791. Ellegren H (2004) Microsatellites: simple sequences with complex evolution. Nat Rev Genet 5:435-445. Graur D, Li WH (2000) Fundamentals of molecular evolution, 2 ed, Sinauer Associates, Sunderland (MA), USA Falconer DS, Mackay TFC (1996) Introduction to quantitative genetics, 4 ed, Longman, Harlow, Essex, UK Graves HB (1984) Behavior and Ecology of Wild and Feral Swine (Sus Scrofa). J Anim Sci 1984 58:482-492. Fang M, Larson G, Ribeiro HS, Li N, Andersson L (2009) Contrasting mode of evolution at a coat color locus in wild and domestic pigs. PLoS Genetics Green P, Falls K, Crooks S (1990) Documentation for CRI-MAP, version 2.4. 5:1-6. Washington University School of Medicine St. Louis (MO), USA

Flajnik MF, Kasahara M (2010) Origin and evolution of the adaptive immune Hajishengallis G, Lambris JD (2010) Crosstalk pathways between Toll-like re- system: genetic events and selective pressures. Nat Rev Genet 11:47-59. ceptors and the complement system. Trends Immunol Feb 10. (Epub ahead of print) Flint J, Mackay TF (2009) Genetic architecture of quantitative traits in mice, flies, and humans. Genome Res 19:723-733.

70 71 Hasan U, Chaffois C, Gaillard C, Saulnier V, Merck E, Tancredi S, Guiet C, Jann OC, King A, Corrales NL, Anderson SI, Jensen K, Ait-Ali T, Tang H, Brière F, Vlach J, Lebecque GT, Bates EEM (2005) Human TLR10 is a func- Wu C, Cockett NE, Archibald AL, Glass EJ (2009) Comparative genomics of tional receptor, expressed by B cells and plasmacytoid dendritic cells, which ac- Toll-like receptor signalling in five species. BMC Genomics 10:216. tivates gene transcription through MyD88. J Immunol 174:2942-2950. Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik SG, Lee H, Lee JO Hawlisch H, Köhl J (2006) Complement and Toll-like receptors: key regula- (2007) Crystal structure of the TLR1-TLR2 heterodimer induced by binding tors of adaptive immune responses. Mol Immunol 43:13-21. of a tri-acylated lipopeptide. Cell 130:1071-1082.

Hibino T, Loza-Coll M, Messier C, Majeske AJ, Cohen AH, Terwilliger DP, Joffre O, Nolte MA, Spörri R, Reis e Sousa C (2009) Inflammatory signals in Buckley KM, Brockton V, Nair SV, Berney K, Fugmann SD, Anderson MK, dendritic cell activation and the induction of adaptive immunity. Immunol Rev Pancer Z, Cameron RA, Smith LC, Rast JP (2006) The immune gene reper- 227:234-247. toire encoded in the purple sea urchin genome. Dev Biol 300:349-365. Juul-Madsen HR, Krogh-Meibom T, Henryon M, Palaniyar N, Heegaard Holmskov U, Thiel S, Jensenius JC (2003) Collections and ficolins: humoral PM, Purup S, Willis AC, Tornøe I, Ingvartsen KL, Hansen S, Holmskov U lectins of the innate immune defense. Annu Rev Immunol 21:547-578. (2006) Identification and characterization of porcine mannan-binding lectin A (pMBL-A), and determination of serum concentration heritability. Immuno- Huang S, Yuan S, Guo L, Yu Y, Li J, Wu T, Liu T, Yang M, Wu K, Liu H, genetics 58:129-137. Ge J, Yu Y, Huang H, Dong M, Yu C, Chen S, Xu A (2008) Genomic analy- sis of the immune gene repertoire of amphioxus reveals extraordinary innate Juul-Madsen HR, Jensen KH, Nielsen J, Damgaard BM (2010) Ontogeny and complexity and diversity. Genome Res 18:1112-1126. characterization of blood leukocyte subsets and serum proteins in piglets before and after weaning. Vet Immunol Immunopathol 133:95-108. Hughes AL, Piontkivska H (2008) Functional diversification of the toll-like receptor gene family. Immunogenetics 60:249-256. Kabelitz D, Medzhitov R (2007) Innate immunity--cross-talk with adaptive immunity through pattern recognition receptors and cytokines. Curr Opin Hummelshoj T, Fog LM, Madsen HO, Sim RB, Garred P (2008) Compara- Immunol 19:1-3. tive study of the human ficolins reveals unique features of Ficolin-3 (Hakata antigen). Mol Immunol 45:1623-1632. Kemper C, Atkinson JP, Hourcade DE (2009) Properdin: Emerging Roles of a Pattern-Recognition Molecule. Annu Rev Immunol Nov 30. (Epub ahead of Ikeda K, Sannoh T, Kawasaki N, Kawasaki T, Yamashina I (1987) Serum print) lectin with known structure activates complement through the classical path- way. J Biol Chem 262:7451-7454. Kijas JMH, Andersson L (2001) A phylogenetic study of the origin of the do- mestic pig estimated from the near-complete mtDNA genome. J Mol Evol Imler JL, Zheng L (2004) Biology of Toll receptors: lessons from insects and 52:302-308. mammals. J Leukoc Biol. 75:18-26. Kimura A, Sakaguchi E, Nonaka M (2009) Multi-component complement Ishii A, Kawasaki M, Matsumoto M, Tochinai S, Seya T (2007) Phylogenetic system of Cnidaria: C3, Bf, and MASP genes expressed in the endodermal tis- and expression analysis of amphibian Xenopus Toll-like receptors. Immunoge- sues of a sea anemone, Nematostella vectensis. Immunobiology 214:165-178. netics 59:281-293. Kruithof EKO, Satta N, Liu JW, Dunoyer-Geindre S, Fish RJ (2007) Gene Iso-Touru T, Kantanen J, Li MH, Gizejewski Z, Vilkki J (2009) Divergent conversion limits divergence of mammalian TLR1 and TLR6. BMC Evol Biol evolution in the cytoplasmic domains of PRLR and GHR genes in Artiodac- 7:148. tyla. BMC Evol Biol 9:172. Kryazhimskiy S, Plotkin JB (2008) The population genetics of dN/dS. PLoS Genetics 4(12):e1000304.

72 73 Hasan U, Chaffois C, Gaillard C, Saulnier V, Merck E, Tancredi S, Guiet C, Jann OC, King A, Corrales NL, Anderson SI, Jensen K, Ait-Ali T, Tang H, Brière F, Vlach J, Lebecque GT, Bates EEM (2005) Human TLR10 is a func- Wu C, Cockett NE, Archibald AL, Glass EJ (2009) Comparative genomics of tional receptor, expressed by B cells and plasmacytoid dendritic cells, which ac- Toll-like receptor signalling in five species. BMC Genomics 10:216. tivates gene transcription through MyD88. J Immunol 174:2942-2950. Jin MS, Kim SE, Heo JY, Lee ME, Kim HM, Paik SG, Lee H, Lee JO Hawlisch H, Köhl J (2006) Complement and Toll-like receptors: key regula- (2007) Crystal structure of the TLR1-TLR2 heterodimer induced by binding tors of adaptive immune responses. Mol Immunol 43:13-21. of a tri-acylated lipopeptide. Cell 130:1071-1082.

Hibino T, Loza-Coll M, Messier C, Majeske AJ, Cohen AH, Terwilliger DP, Joffre O, Nolte MA, Spörri R, Reis e Sousa C (2009) Inflammatory signals in Buckley KM, Brockton V, Nair SV, Berney K, Fugmann SD, Anderson MK, dendritic cell activation and the induction of adaptive immunity. Immunol Rev Pancer Z, Cameron RA, Smith LC, Rast JP (2006) The immune gene reper- 227:234-247. toire encoded in the purple sea urchin genome. Dev Biol 300:349-365. Juul-Madsen HR, Krogh-Meibom T, Henryon M, Palaniyar N, Heegaard Holmskov U, Thiel S, Jensenius JC (2003) Collections and ficolins: humoral PM, Purup S, Willis AC, Tornøe I, Ingvartsen KL, Hansen S, Holmskov U lectins of the innate immune defense. Annu Rev Immunol 21:547-578. (2006) Identification and characterization of porcine mannan-binding lectin A (pMBL-A), and determination of serum concentration heritability. Immuno- Huang S, Yuan S, Guo L, Yu Y, Li J, Wu T, Liu T, Yang M, Wu K, Liu H, genetics 58:129-137. Ge J, Yu Y, Huang H, Dong M, Yu C, Chen S, Xu A (2008) Genomic analy- sis of the immune gene repertoire of amphioxus reveals extraordinary innate Juul-Madsen HR, Jensen KH, Nielsen J, Damgaard BM (2010) Ontogeny and complexity and diversity. Genome Res 18:1112-1126. characterization of blood leukocyte subsets and serum proteins in piglets before and after weaning. Vet Immunol Immunopathol 133:95-108. Hughes AL, Piontkivska H (2008) Functional diversification of the toll-like receptor gene family. Immunogenetics 60:249-256. Kabelitz D, Medzhitov R (2007) Innate immunity--cross-talk with adaptive immunity through pattern recognition receptors and cytokines. Curr Opin Hummelshoj T, Fog LM, Madsen HO, Sim RB, Garred P (2008) Compara- Immunol 19:1-3. tive study of the human ficolins reveals unique features of Ficolin-3 (Hakata antigen). Mol Immunol 45:1623-1632. Kemper C, Atkinson JP, Hourcade DE (2009) Properdin: Emerging Roles of a Pattern-Recognition Molecule. Annu Rev Immunol Nov 30. (Epub ahead of Ikeda K, Sannoh T, Kawasaki N, Kawasaki T, Yamashina I (1987) Serum print) lectin with known structure activates complement through the classical path- way. J Biol Chem 262:7451-7454. Kijas JMH, Andersson L (2001) A phylogenetic study of the origin of the do- mestic pig estimated from the near-complete mtDNA genome. J Mol Evol Imler JL, Zheng L (2004) Biology of Toll receptors: lessons from insects and 52:302-308. mammals. J Leukoc Biol. 75:18-26. Kimura A, Sakaguchi E, Nonaka M (2009) Multi-component complement Ishii A, Kawasaki M, Matsumoto M, Tochinai S, Seya T (2007) Phylogenetic system of Cnidaria: C3, Bf, and MASP genes expressed in the endodermal tis- and expression analysis of amphibian Xenopus Toll-like receptors. Immunoge- sues of a sea anemone, Nematostella vectensis. Immunobiology 214:165-178. netics 59:281-293. Kruithof EKO, Satta N, Liu JW, Dunoyer-Geindre S, Fish RJ (2007) Gene Iso-Touru T, Kantanen J, Li MH, Gizejewski Z, Vilkki J (2009) Divergent conversion limits divergence of mammalian TLR1 and TLR6. BMC Evol Biol evolution in the cytoplasmic domains of PRLR and GHR genes in Artiodac- 7:148. tyla. BMC Evol Biol 9:172. Kryazhimskiy S, Plotkin JB (2008) The population genetics of dN/dS. PLoS Genetics 4(12):e1000304.

72 73 Larsen PH, Holm TH, Owens T (2007) Toll-like receptors in brain develop- Ma Y, Li J, Chiu I, Wang Y, Sloane JA, Lü J, Kosaras B, Sidman RL, Volpe ment and homeostasis. Sci. STKE pe47. JJ, Vartanian T (2006) Toll-like receptor 8 functions as a negative regulator of neurite outgrowth and inducer of neuronal apoptosis. J Cell Biol 175:209-215. Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, Robins J, Lowden S, Finlayson H, Brand T, Willerslev E, Rowley-Conwy P, Andersson Mackay TF, Stone EA, Ayroles JF (2009) The genetics of quantitative traits: L, Cooper A (2005) Worldwide phylogeography of wild boar reveals multiple challenges and prospects. Nat Rev Genet 10:565-577. centers of pig domestication. Science 307:1618-1621. Manfredi AA, Rovere-Querini P, Bottazzi B, Garlanda C, Mantovani A Lee HK, Dunzendorfer S, Soldau K, Tobias PS (2006) Double-stranded (2008) Pentraxins, humoral innate immunity and tissue injury. Curr Opin Im- RNA-mediated TLR3 activation is enhanced by CD14. Immunity 24:153- munol 20:538-544. 163. Matsushima N, Tanaka T, Enkhbayar P, Mikami T, Taga M, Yamada K, Ku- Lee J, Gonzales-Navajas JM, Raz E (2008) The "polarizing-tolerizing" mecha- roki Y (2007) Comparative sequence analysis of leucine-rich repeats (LRRs) nism of intestinal epithelium: its relevance to colonic homeostasis. Semin Im- within vertebrate toll-like receptors. BMC Genomics 8:124. munopathol 30:3-9. Medzhitov R, Preston-Hurlburt P, Janeway Jr CA (1997) A human homo- Lemaitre B (2004) The road to Toll. Nat Rev Immunol 4:521-527. logue of the Drosophila Toll protein signals activation of adaptive immunity. Nature 388:394-397. Lemaitre B, Hoffmann J (2007) The host defense of Drosophila melanogaster. Annu Rev Immunol 25:697-743. Medzhitov R (2001) Toll-like receptors and innate immunity. Nat Rev Immu- nol 1:135-145. Letunic I, Doerks T, Bork P (2008) SMART 6: recent updates and new devel- opments. Nucleic Acids Res 37:D229-232. Medzhitov R, Janeway Jr CA (2002) Decoding the patterns of self and nonself by the innate immune system. Science 296:298-300. Lillie BN, Brooks AS, Keirstead ND, Hayes MA (2005) Comparative genetics and innate immune functions of collagenous lectins in animals. Vet Immunol Megens HJ, Crooijmans RP, San Cristobal M, Hui X, Li N, Groenen MA Immunopathol 108:97-110. (2008) Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas Lillie BN, Hammermueller JD, Macinnes JI, Jacques M, Hayes MA (2006a) of domestication. Genet Sel Evol 40:103-128. Porcine mannan-binding lectin A binds to Actinobacillus suis and Haemophi- lus parasuis. Dev Comp Immunol 30:954-965. Mendel G (s a) Experiments in plant hybridization, Harvard University Press, Cambridge (MA), USA Lillie BN, Keirstead ND, Squires EJ, Hayes MA (2006b) Single-nucleotide polymorphisms in porcine mannan-binding lectin A. Immunogenetics 58:983- Mongini PK, Vilensky MA, Highet PF, Inman JK (1997) The affinity thresh- 993. old for human B cell activation via the antigen receptor complex is reduced upon co-ligation of the antigen receptor with CD21 (CR2). J Immunol Lillie BN, Keirstead ND, Squires EJ, Hayes MA (2007) Gene polymorphisms 159:3782-3791. associated with reduced hepatic expression of porcine mannan-binding lectin C. Dev Comp Immunol 31:830-846. Morozumi T, Uenishi H (2009) Polymorphism distribution and structural conservation in RNA-sensing Toll-like receptors 3, 7, and 8 in pigs. Biochim Luo S, Poltermann S, Kunert A, Rupp S, Zipfel PF (2009) Immune evasion of Biophys Acta 1790:267-274. the human pathogenic yeast Candida albicans: Pra1 is a Factor H, FHL-1 and plasminogen binding surface protein. Mol Immunol 47:541-550. Nicholas FW (2005) Animal breeding and disease. Philos Trans R Soc Lond B Biol Sci 360:1529-1536.

74 75 Larsen PH, Holm TH, Owens T (2007) Toll-like receptors in brain develop- Ma Y, Li J, Chiu I, Wang Y, Sloane JA, Lü J, Kosaras B, Sidman RL, Volpe ment and homeostasis. Sci. STKE pe47. JJ, Vartanian T (2006) Toll-like receptor 8 functions as a negative regulator of neurite outgrowth and inducer of neuronal apoptosis. J Cell Biol 175:209-215. Larson G, Dobney K, Albarella U, Fang M, Matisoo-Smith E, Robins J, Lowden S, Finlayson H, Brand T, Willerslev E, Rowley-Conwy P, Andersson Mackay TF, Stone EA, Ayroles JF (2009) The genetics of quantitative traits: L, Cooper A (2005) Worldwide phylogeography of wild boar reveals multiple challenges and prospects. Nat Rev Genet 10:565-577. centers of pig domestication. Science 307:1618-1621. Manfredi AA, Rovere-Querini P, Bottazzi B, Garlanda C, Mantovani A Lee HK, Dunzendorfer S, Soldau K, Tobias PS (2006) Double-stranded (2008) Pentraxins, humoral innate immunity and tissue injury. Curr Opin Im- RNA-mediated TLR3 activation is enhanced by CD14. Immunity 24:153- munol 20:538-544. 163. Matsushima N, Tanaka T, Enkhbayar P, Mikami T, Taga M, Yamada K, Ku- Lee J, Gonzales-Navajas JM, Raz E (2008) The "polarizing-tolerizing" mecha- roki Y (2007) Comparative sequence analysis of leucine-rich repeats (LRRs) nism of intestinal epithelium: its relevance to colonic homeostasis. Semin Im- within vertebrate toll-like receptors. BMC Genomics 8:124. munopathol 30:3-9. Medzhitov R, Preston-Hurlburt P, Janeway Jr CA (1997) A human homo- Lemaitre B (2004) The road to Toll. Nat Rev Immunol 4:521-527. logue of the Drosophila Toll protein signals activation of adaptive immunity. Nature 388:394-397. Lemaitre B, Hoffmann J (2007) The host defense of Drosophila melanogaster. Annu Rev Immunol 25:697-743. Medzhitov R (2001) Toll-like receptors and innate immunity. Nat Rev Immu- nol 1:135-145. Letunic I, Doerks T, Bork P (2008) SMART 6: recent updates and new devel- opments. Nucleic Acids Res 37:D229-232. Medzhitov R, Janeway Jr CA (2002) Decoding the patterns of self and nonself by the innate immune system. Science 296:298-300. Lillie BN, Brooks AS, Keirstead ND, Hayes MA (2005) Comparative genetics and innate immune functions of collagenous lectins in animals. Vet Immunol Megens HJ, Crooijmans RP, San Cristobal M, Hui X, Li N, Groenen MA Immunopathol 108:97-110. (2008) Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas Lillie BN, Hammermueller JD, Macinnes JI, Jacques M, Hayes MA (2006a) of domestication. Genet Sel Evol 40:103-128. Porcine mannan-binding lectin A binds to Actinobacillus suis and Haemophi- lus parasuis. Dev Comp Immunol 30:954-965. Mendel G (s a) Experiments in plant hybridization, Harvard University Press, Cambridge (MA), USA Lillie BN, Keirstead ND, Squires EJ, Hayes MA (2006b) Single-nucleotide polymorphisms in porcine mannan-binding lectin A. Immunogenetics 58:983- Mongini PK, Vilensky MA, Highet PF, Inman JK (1997) The affinity thresh- 993. old for human B cell activation via the antigen receptor complex is reduced upon co-ligation of the antigen receptor with CD21 (CR2). J Immunol Lillie BN, Keirstead ND, Squires EJ, Hayes MA (2007) Gene polymorphisms 159:3782-3791. associated with reduced hepatic expression of porcine mannan-binding lectin C. Dev Comp Immunol 31:830-846. Morozumi T, Uenishi H (2009) Polymorphism distribution and structural conservation in RNA-sensing Toll-like receptors 3, 7, and 8 in pigs. Biochim Luo S, Poltermann S, Kunert A, Rupp S, Zipfel PF (2009) Immune evasion of Biophys Acta 1790:267-274. the human pathogenic yeast Candida albicans: Pra1 is a Factor H, FHL-1 and plasminogen binding surface protein. Mol Immunol 47:541-550. Nicholas FW (2005) Animal breeding and disease. Philos Trans R Soc Lond B Biol Sci 360:1529-1536.

74 75 Nürnberger T, Brunner F, Kemmerling B, Piater L (2004) Innate immunity in Ruiz-Fons F, Segalés J, Gortázar C (2008) A review of viral diseases of the plants and animals: striking similarities and obvious differences. Immunol Rev European wild boar: effects of population dynamics and reservoir rôle. Vet J. 198:249-266. 2008 176:158-169.

Ojeda A, Huang LS, Ren J, Angiolillo A, Cho IC, Soto H, Lemús-Flores C, Schou TW, Permin A, Juul-Madsen HR, Sørensen P, Labouriau R, Nguyên Makuza SM, Folch JM, Pérez-Enciso M (2008) Selection in the making: a TL, Fink M, Pham SL (2007) Gastrointestinal helminths in indigenous and worldwide survey of haplotypic diversity around a causative mutation in porcine exotic chickens in Vietnam: association of the intensity of infection with the IGF2. Genetics 178:1639-1652. Major Histocompatibility Complex. Parasitology 134(Pt 4):561-573.

Palm NW, Medzhitov R (2009) Pattern recognition receptors and control of Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM (2002) QTL Ex- adaptive immunity. Immunol Rev 227:221-233. press: mapping quantitative trait loci in simple and complex pedigrees. Bioin- formatics 18:339-340. Pancer Z, Amemiya CT, Ehrhardt GR, Ceitlin J, Gartland GL, Cooper MD (2004) Somatic diversification of variable lymphocyte receptors in the agnathan Shaw MH, Reimer T, Kim YG, Nuñez G (2008) NOD-like receptors sea lamprey. Nature 430:174-180. (NLRs): bona fide intracellular microbial sensors. Curr Opin Immunol 20:377- 382. Phatsara C, Jennen DG, Ponsuksili S, Murani E, Tesfaye D, Schellander K, Wimmers K (2007) Molecular genetic analysis of porcine mannose-binding Shinkai H, Muneta Y, Suzuki K, Eguchi-Ogawa T, Awata T, Uenishi H lectin genes, MBL1 and MBL2, and their association with complement activ- (2006a) Porcine Toll-like receptor 1, 6, and 10 genes: Complete sequencing of ity. Int J Immunogenet 34:55-63. genomic region and expression analysis. Mol Immunol 43:1474-1480.

Portnoy DA (2005) Manipulation of innate immunity by bacterial pathogens. Shinkai H, Tanaka M, Morozumi T, Eguchi-Ogawa T, Okumura N, Muneta Curr Opin Immunol 17:25-28. Y, Awata T, Uenishi H (2006b) Biased distribution of single nucleotide poly- morphisms (SNPs) in porcine Toll-like receptor 1 (TLR1), TLR2, TLR4, Ranjan P, Bowzard JB, Schwermann JW, Jeisy-Scott V, Fujita T, Sambhara S TLR5, and TLR6 genes. Immunogenetics 58:324-330. (2009) Cytoplasmic nucleic acid sensors in antiviral immunity. Trends Mol Med 15:359-368. Skjoedt MO, Hummelshoj T, Palarasah Y, Honore C, Koch C, Skjodt K, Garred P (2010) A novel MBL/ficolin associated protein is highly expressed in Roach JC, Glusman G, Rowen L, Kaur A, Purcell MK, Smith KD, Hood LE, heart and skeletal muscle tissues and inhibits complement activation. J Biol Aderem A (2005) The evolution of vertebrate Toll-like receptors. Pros Natl Chem 285:8234-8243. Acad Sci U S A 102:9577-9582. Takahashi M, Iwaki D, Matsushita A, Nakata M, Matsushita M, Endo Y, Fu- Ron M, Weller JI (2007) From QTL to QTN identification in livestock-- jita T (2006) Cloning and characterization of mannose-binding lectin from winning by points rather than knock-out: a review. Anim Genet 38:429-439. lamprey (Agnathans). J Immunol 176:4861-4868.

Roos A, Garred P, Wildenberg ME, Lynch NJ, Munoz JR, Zuiverloon TC, Takahashi M, Ishida Y, Iwaki D, Kanno K, Suzuki T, Endo Y, Homma Y, Bouwman LH, Schlagwein N, Fallaux van den Houten FC, Faber-Krol MC, Fujita T (2010) Essential role of mannose-binding lectin-associated serine pro- Madsen HO, Schwaeble WJ, Matsushita M, Fujita T, Daha MR (2004) Anti- tease-1 in activation of the complement factor D. J Exp Med 207:29-37. body-mediated activation of the classical pathway of complement may compen- sate for mannose-binding lectin deficiency. Eur J Immunol 34:2589-2598. Takeda K, Akira S (2005) Toll-like receptors in innate immunity. Int Immunol 17:1-14. Rose D, Zhu X, Kose H, Hoang B, Cho J, Chiba A (1997) Toll, a muscle cell surface molecule, locally inhibits synaptic initiation of the RP3 motoneuron growth cone in Drosophila. Development 124:1561-1571.

76 77 Nürnberger T, Brunner F, Kemmerling B, Piater L (2004) Innate immunity in Ruiz-Fons F, Segalés J, Gortázar C (2008) A review of viral diseases of the plants and animals: striking similarities and obvious differences. Immunol Rev European wild boar: effects of population dynamics and reservoir rôle. Vet J. 198:249-266. 2008 176:158-169.

Ojeda A, Huang LS, Ren J, Angiolillo A, Cho IC, Soto H, Lemús-Flores C, Schou TW, Permin A, Juul-Madsen HR, Sørensen P, Labouriau R, Nguyên Makuza SM, Folch JM, Pérez-Enciso M (2008) Selection in the making: a TL, Fink M, Pham SL (2007) Gastrointestinal helminths in indigenous and worldwide survey of haplotypic diversity around a causative mutation in porcine exotic chickens in Vietnam: association of the intensity of infection with the IGF2. Genetics 178:1639-1652. Major Histocompatibility Complex. Parasitology 134(Pt 4):561-573.

Palm NW, Medzhitov R (2009) Pattern recognition receptors and control of Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM (2002) QTL Ex- adaptive immunity. Immunol Rev 227:221-233. press: mapping quantitative trait loci in simple and complex pedigrees. Bioin- formatics 18:339-340. Pancer Z, Amemiya CT, Ehrhardt GR, Ceitlin J, Gartland GL, Cooper MD (2004) Somatic diversification of variable lymphocyte receptors in the agnathan Shaw MH, Reimer T, Kim YG, Nuñez G (2008) NOD-like receptors sea lamprey. Nature 430:174-180. (NLRs): bona fide intracellular microbial sensors. Curr Opin Immunol 20:377- 382. Phatsara C, Jennen DG, Ponsuksili S, Murani E, Tesfaye D, Schellander K, Wimmers K (2007) Molecular genetic analysis of porcine mannose-binding Shinkai H, Muneta Y, Suzuki K, Eguchi-Ogawa T, Awata T, Uenishi H lectin genes, MBL1 and MBL2, and their association with complement activ- (2006a) Porcine Toll-like receptor 1, 6, and 10 genes: Complete sequencing of ity. Int J Immunogenet 34:55-63. genomic region and expression analysis. Mol Immunol 43:1474-1480.

Portnoy DA (2005) Manipulation of innate immunity by bacterial pathogens. Shinkai H, Tanaka M, Morozumi T, Eguchi-Ogawa T, Okumura N, Muneta Curr Opin Immunol 17:25-28. Y, Awata T, Uenishi H (2006b) Biased distribution of single nucleotide poly- morphisms (SNPs) in porcine Toll-like receptor 1 (TLR1), TLR2, TLR4, Ranjan P, Bowzard JB, Schwermann JW, Jeisy-Scott V, Fujita T, Sambhara S TLR5, and TLR6 genes. Immunogenetics 58:324-330. (2009) Cytoplasmic nucleic acid sensors in antiviral immunity. Trends Mol Med 15:359-368. Skjoedt MO, Hummelshoj T, Palarasah Y, Honore C, Koch C, Skjodt K, Garred P (2010) A novel MBL/ficolin associated protein is highly expressed in Roach JC, Glusman G, Rowen L, Kaur A, Purcell MK, Smith KD, Hood LE, heart and skeletal muscle tissues and inhibits complement activation. J Biol Aderem A (2005) The evolution of vertebrate Toll-like receptors. Pros Natl Chem 285:8234-8243. Acad Sci U S A 102:9577-9582. Takahashi M, Iwaki D, Matsushita A, Nakata M, Matsushita M, Endo Y, Fu- Ron M, Weller JI (2007) From QTL to QTN identification in livestock-- jita T (2006) Cloning and characterization of mannose-binding lectin from winning by points rather than knock-out: a review. Anim Genet 38:429-439. lamprey (Agnathans). J Immunol 176:4861-4868.

Roos A, Garred P, Wildenberg ME, Lynch NJ, Munoz JR, Zuiverloon TC, Takahashi M, Ishida Y, Iwaki D, Kanno K, Suzuki T, Endo Y, Homma Y, Bouwman LH, Schlagwein N, Fallaux van den Houten FC, Faber-Krol MC, Fujita T (2010) Essential role of mannose-binding lectin-associated serine pro- Madsen HO, Schwaeble WJ, Matsushita M, Fujita T, Daha MR (2004) Anti- tease-1 in activation of the complement factor D. J Exp Med 207:29-37. body-mediated activation of the classical pathway of complement may compen- sate for mannose-binding lectin deficiency. Eur J Immunol 34:2589-2598. Takeda K, Akira S (2005) Toll-like receptors in innate immunity. Int Immunol 17:1-14. Rose D, Zhu X, Kose H, Hoang B, Cho J, Chiba A (1997) Toll, a muscle cell surface molecule, locally inhibits synaptic initiation of the RP3 motoneuron growth cone in Drosophila. Development 124:1561-1571.

76 77 Takeuchi O, Kawai T, Muhlradt PF, Morr M, Radolf JD, Zychlinsky A, Ta- Zhang X, Kimura Y, Fang C, Zhou L, Sfyroera G, Lambris JD, Wetsel RA, keda K, Akira S (2001) Discrimination of bacterial lipoproteins by Toll-like re- Miwa T, Song WC (2007) Regulation of Toll-like receptor-mediated inflam- ceptor 6. Int Immunol 13:933-940. matory response by complement in vivo. Blood 110:228-236.

Takeuchi O, Sato S, Horiuchi T, Hoshino K, Takeda K, Dong Z, Modlin RL, Zhong F, Cao W, Chan E, Tay PN, Cahya FF, Zhang H, Lu, J (2005) Devia- Akira S (2002) Cutting edge: role of Toll-like receptor 1 in mediating immune tion from major codons in the Toll-like receptor genes is associated with low response to microbial lipoproteins. J Immunol 169:10-4. Toll-like receptor expression. Immunology 114:83-93.

Tenor JL, Aballay A (2008) A conserved Toll-like receptor is required for Zänker KS (2008) General introduction to innate immunity: Dr. Jekyl/Mr. Caenorhabditis elegans innate immunity. EMBO Rep 9:103-109. Hyde quality of the innate immune system. in: Egesten A, Schmidt A, Her- wald H (eds): Trends in Innate Immunity. Contrib Microbiol vol 15, pp 12-20, Thiel S, Holmskov U, Hviid L, Laursen SB, Jensenius JC (1992) The concen- Karger, Basel, Switzerland tration of the C-type lectin, mannan-binding protein, in human plasma in- creases during an acute phase response. Clin Exp Immunol 90:31-35.

Uenishi H, Shinkai H (2008) Porcine Toll-like receptors: the front line of pathogen monitoring and possible implications for disease resistance. Dev Comp Immunol 33:353-361.

Verdu P, Barreiro LB, Patin E, Gessain A, Cassar O, Kidd JR, Kidd KK, Be- har DM, Froment A, Heyer E, Sica L, Casanova JL, Abel L, Quintana-Murci L (2006) Evolutionary insights into the high worldwide prevalence of MBL2 deficiency alleles. Hum Mol Genet 15:2650-2658.

Visintin A, Mazzoni A, Spitzer JH, Wyllie DH, Dower SK, Segal DM (2001) Regulation of Toll-like receptors in human monocytes and dendritic cells. J Immunol 166:249-255.

Wang J, Shao Y, Bennett TA, Shankar RA, Wightman PD, Reddy LG (2006) The functional effects of physical interactions among Toll-like receptors 7, 8, and 9. J Biol Chem 281:37427-37434.

Wattrang E, Almqvist M, Johansson A, Fossum C, Wallgren P, Pielberg G, Andersson L, Edfors-Lilja I (2005) Confirmation of QTL on porcine chromo- somes 1 and 8 influencing leukocyte numbers, haematological parameters and leukocyte function. Anim Genet 36:337-345.

Werling D, Jann OC, Offord V, Glass EJ, Coffey TJ (2009) Variation matters: TLR structure and species-specific pathogen recognition. Trends Immunol 30:124-130.

Xu Y, Tao X, Shen B, Horng T, Medzhitov R, Manley JL, Tong L (2000) Structural basis for signal transduction by the Toll/interleukin-1 receptor do- mains. Nature 408:111-115. Front cover photograph by Inger Edfors (private).

78 79 Takeuchi O, Kawai T, Muhlradt PF, Morr M, Radolf JD, Zychlinsky A, Ta- Zhang X, Kimura Y, Fang C, Zhou L, Sfyroera G, Lambris JD, Wetsel RA, keda K, Akira S (2001) Discrimination of bacterial lipoproteins by Toll-like re- Miwa T, Song WC (2007) Regulation of Toll-like receptor-mediated inflam- ceptor 6. Int Immunol 13:933-940. matory response by complement in vivo. Blood 110:228-236.

Takeuchi O, Sato S, Horiuchi T, Hoshino K, Takeda K, Dong Z, Modlin RL, Zhong F, Cao W, Chan E, Tay PN, Cahya FF, Zhang H, Lu, J (2005) Devia- Akira S (2002) Cutting edge: role of Toll-like receptor 1 in mediating immune tion from major codons in the Toll-like receptor genes is associated with low response to microbial lipoproteins. J Immunol 169:10-4. Toll-like receptor expression. Immunology 114:83-93.

Tenor JL, Aballay A (2008) A conserved Toll-like receptor is required for Zänker KS (2008) General introduction to innate immunity: Dr. Jekyl/Mr. Caenorhabditis elegans innate immunity. EMBO Rep 9:103-109. Hyde quality of the innate immune system. in: Egesten A, Schmidt A, Her- wald H (eds): Trends in Innate Immunity. Contrib Microbiol vol 15, pp 12-20, Thiel S, Holmskov U, Hviid L, Laursen SB, Jensenius JC (1992) The concen- Karger, Basel, Switzerland tration of the C-type lectin, mannan-binding protein, in human plasma in- creases during an acute phase response. Clin Exp Immunol 90:31-35.

Uenishi H, Shinkai H (2008) Porcine Toll-like receptors: the front line of pathogen monitoring and possible implications for disease resistance. Dev Comp Immunol 33:353-361.

Verdu P, Barreiro LB, Patin E, Gessain A, Cassar O, Kidd JR, Kidd KK, Be- har DM, Froment A, Heyer E, Sica L, Casanova JL, Abel L, Quintana-Murci L (2006) Evolutionary insights into the high worldwide prevalence of MBL2 deficiency alleles. Hum Mol Genet 15:2650-2658.

Visintin A, Mazzoni A, Spitzer JH, Wyllie DH, Dower SK, Segal DM (2001) Regulation of Toll-like receptors in human monocytes and dendritic cells. J Immunol 166:249-255.

Wang J, Shao Y, Bennett TA, Shankar RA, Wightman PD, Reddy LG (2006) The functional effects of physical interactions among Toll-like receptors 7, 8, and 9. J Biol Chem 281:37427-37434.

Wattrang E, Almqvist M, Johansson A, Fossum C, Wallgren P, Pielberg G, Andersson L, Edfors-Lilja I (2005) Confirmation of QTL on porcine chromo- somes 1 and 8 influencing leukocyte numbers, haematological parameters and leukocyte function. Anim Genet 36:337-345.

Werling D, Jann OC, Offord V, Glass EJ, Coffey TJ (2009) Variation matters: TLR structure and species-specific pathogen recognition. Trends Immunol 30:124-130.

Xu Y, Tao X, Shen B, Horng T, Medzhitov R, Manley JL, Tong L (2000) Structural basis for signal transduction by the Toll/interleukin-1 receptor do- mains. Nature 408:111-115. Front cover photograph by Inger Edfors (private).

78 79