Volume 5 Supplement 4
The three-way relationship of polymorphisms of porcine genes encoding terminal complement components, their differential expression, and health-related phenotypes
© Wimmers et al; licensee BioMed Central Ltd. 2011
Published: 3 June 2011
The complement system is an evolutionary ancient mechanism that plays an essential role in innate immunity and contributes to the acquired immune response. Three modes of activation, known as classical, alternative and lectin pathway, lead to the initiation of a common terminal lytic pathway. The terminal complement components (TCCs: C6, C7, C8A, C8B, and C9) are encoded by the genes C6, C7, C8A, C8B, C8G, and C9. We aimed at experimentally testing the porcine genes encoding TCCs as candidate genes for immune competence and disease resistance by addressing the three-way relationship of genotype, health related phenotype, and mRNA expression.
Comparative sequencing of cDNAs of animals of the breeds German Landrace, Piétrain, Hampshire, Duroc, Vietnamese Potbelly Pig, and Berlin Miniature Pig (BMP) revealed 30 SNPs (21 in protein domains, 12 with AA exchange). The promoter regions (each ~1.5 kb upstream the transcription start sites) of C6, C7, C8A, C8G, and C9 exhibited 29 SNPs. Significant effects of the TCC encoding genes on hemolytic complement activity were shown in a cross of Duroc and BMP after vaccination against Mycoplasma hyopneumoniae, Aujeszky disease virus and PRRSV by analysis of variance using repeated measures mixed models. Family based association tests (FBAT) confirmed the associations. The promoter SNPs were associated with the relative abundance of TCC transcripts obtained by real time RT-PCR of 311 liver samples of commercial slaughter pigs. Complement gene expression showed significant relationship with the prevalence of acute and chronic lung lesions.
The analyses point to considerable variation of the porcine TCC genes and promote the genes as candidate genes for disease resistance.
Infectious diseases in farm animals and their control by therapeutic, prophylactic and metaphylactic treatments and veterinary management schemes cause immense costs and animal welfare concerns and increasingly burden man and environment with chemotherapeutics. Genetically improved health and resistance status of farm animals is associated with higher product quality, reduced danger of transmission of infectious agents or traces of antibiotics together with its positive impact on production costs and animal welfare. Breeding for specific and/or global disease resistance is hampered by the difficulty to discriminate the phenotypes `resistant' and `susceptible' either after experimental or natural infections. Alternatively, it has been proposed to measure immune response traits as indicator traits of improved defense power. However, determination of high and low immune responsiveness is still very inefficient and too complex for most farm animal breeding programs. Consequently, the availability of DNA-markers for disease resistance would be particularly useful in farm animal breeding. Since infectious factorial diseases occur in all porcine age groups and are especially responsible for economic losses and high mortality-rates in weaners and young feeder pigs affected by respiratory and enteritic diseases we focus on broad resistance against various pathogens.
The complement system is a mechanism of the innate immune system capable of direct killing of microorganisms and modulating phagocytosis, inflammation, humoral and cellular immune response. The complement system is a highly regulated and complex set of interacting proteins in blood plasma and on cells surfaces. There are three pathways of complement activation, the classical, the alternative and the lectin pathway. Activation of the complement system through the classical pathway depends on the presence of immune complexes and links the innate immune system to the adaptive. The alternative and lectin pathway represent major means of natural resistance to microorganisms in the absence of specific antibodies. All the three pathways lead to the initiation of a common terminal lytic pathway and the formation of the membrane attack complex (MAC), which causes lysis of the invading pathogens/microbes and is an assembly of the homologous components C5b, C6, C7, C8, and C9. The terminal complement components (TCCs) are encoded by the genes C6, C7, C8A, C8B, C8G, and C9 originating from an ancestral gene. We have previously shown association of complement genes of the initiation pathways and other immune relevant genes with complement activity and acute phase response [1–4]. Here we aimed at structural and functional characterization of the porcine genes encoding TCCs and experimentally qualifying TCCs as candidate genes for immune competence and disease resistance by addressing the three-way relationship of genotype, health related phenotype, and mRNA expression.
Experimental population: Animals (n=457) of an experimental F2 population based on reciprocal crossbreeding of Duroc and Berlin Miniature Pig (DUMI)  were immunized with Mycoplasma hyopneumoniae (Mh), Aujeszky disease (ADV), and PRRS (PRRSV) vaccines at 6, 14 and 16 weeks of age, respectively. Blood samples were taken immediately prior to immunisation (day 0) and at day 4 and 10 after Mh and ADV vaccination, but only at day 10 days after PRRSV vaccination. Hemolytic complement activities in the classical (CH50) and alternative pathways (AH50), C3c (C3c) and haptoglobin (HP) acute phase protein levels were measured . The DUMI animals were genotyped at 88 loci (72 microsatellites, 16 biallelic markers) covering the porcine autosomes. Linkage analysis was performed using CRI MAP . The QTL analysis was done using QTLexpress  under the line cross and the half sib model and after adjustment of phenotypic data for systematic effects by mixed models with and without repeated measures statement. Association and linkage of SNPs detected in the cDNA sequence of C6, C7, C8A, C8B, and C9 were tested by means of analysis of variance using repeated measures mixed models and family based association test (FBAT) testing the null hypothesis of no association in the presence of linkage (option `eÂ´, computing the test statistic using the empirical variance ).
Commercial herds: Animals (n=311) of commercial herds of the breeds German Landrace, Large White and their crossbreeds (F1) and Piétrain × F1 were kept and performance tested at the Pig performance test station Jürgenstorf and slaughtered at the experimental abattoir of the Leibniz Institute of Farm Animals Biology (FBN). Liver tissue was immediately sampled after exsanguination for RNA isolation using Tri-Reagent (Sigma-Aldrich, Taufkirchen, Germany) and NucleoSpin RNA II kit (Macherey-Nagel, Düren, Germany) including DNase treatment. Lung lesions were recorded during veterinary inspection and relative abundance of transcripts of C6, C7, C8A, C8B, C8G, and C9 was determined in liver tissue. Therefore, real time PCR was performed using the LightCycler 480 system (Roche, Mannheim, Germany) with cDNA synthesized from 1µg of total RNA using random and oligo d(T) 13VN primers, Superscript III reverse transcriptase (Invitrogen, Karlsruhe, Germany), and the LightCycler 480 SYBR Green I Master kit (Roche). Standard curves were generated by amplifying serial dilutions of specific PCR products. Normalization of variation in RT-PCR efficiency and initial RNA input was performed using the RPL32 and ALB genes as internal standards by dividing the calculated mRNA copy numbers by a mean normalization factor derived from the expression of the reference genes (relative abundance). Alternatively, ΔCT was calculated by subtracting the mean ct-values of the reference genes from the ct-values obtained for single TCC.
Association of SNPs detected in the promoters or cDNAs of C6, C7, C8A, and C9 with relative transcript abundance were analyzed by mixed model analysis of variance (SAS, Proc mixed) taking into account the fixed effects of breed, farm where the piglets were born and raised until seven weeks of age when they were transported to the performance test station, genotype of the respective gene and the random effect of father. The relationship between the respective genotypes and the prevalence of lung lesions was analyzed by case control analysis using JMP genomics. The relationship between the transcript abundance and lung lesions was evaluated by analysis of variance taking breed and lung lesions (affected, non-affected) as fixed effects on the dependent variable of relative transcript abundance.
Results and discussion
For genetic mapping of the TCCs and of QTL for complement activity and acute phase response as well as studying association and linkage of the complement genes animals of the experimental population DUMI were examined where multiple vaccinations and samplings were made. The QTL analysis for hemolytic complement activity in the classical and alternative pathway and acute phase response, i.e. C3c and HP serum concentrations, revealed 67 QTL with p < .05 genome-wide significance. The proximal region of SSC2 (orthologous to HSA11 0-70 Mb), the distal region of SSC4 (HSA1 95-155 Mb), and the intermediate region of SSC16 (HSA5 0-73 Mb & 150-174 Mb) showed several QTL clustered together . In these regions a single true QTL might underlie the various QTL positions estimated for related traits. With regard to the number and the magnitude of their impact, QTL for humoral innate immune traits behave like those for other quantitative traits. Radiation hybrid and genetic mapping of complement genes (C8G: SSC1q2.13; C8A and C8b: SSC6q3.1-q3.5; C6, C7, C9: SSC 16q14) highlights them as positional candidate genes for the QTL .
In order to analyze association and linkage of the TCCs with complement activity the cDNA sequences were obtained and SNPs were detected by comparative sequencing of animals of the breeds German Landrace, Piétrain, Hampshire, Duroc, Vietnamese Potbelly Pig, and Berlin Miniature Pig (BMP). In total 30 SNPs were detected (21 in protein domains, 12 nonsynonymous) in the cDNA sequences of the genes C6 (Acc. # DQ333199), C7 (Acc. # AF162274; ), C8A (Acc. # DQ333200), C8B (Acc. # DQ333201), C8G (Acc. # DQ333202) and C9 (Acc # DQ333198). Analysis of variance revealed association of C7 (c.881A>G) and C9 (c.407C>T) with complement activity at p<.05 and for C8A (c.1544C>T) at p=0.06 . Family based association tests (FBAT) confirmed the associations (p<0.05). The results indicate considerable polymorphism of the coding region of porcine TCCs and promote the TCCs as candidate genes for innate humoral defense.
The analyses in the experimental F2 population allowed multiple vaccinations and samplings, mapping of QTL for complement activity and acute phase response as well as detection of association and linking of complement genes. However, in order to evaluate the results and to get closer to commercial conditions in terms of breeds and pedigrees - especially in view of the available SNP chips - and also diseases, further studies were focused on commercial herds, where relative transcript abundance of TCCs and lung lesions were used as traits related to health and the three-way relationship of genotype, health related phenotype, and mRNA expression was explored. Consequently, focus was on promoter regions for detection of further SNPs. The promoter regions (each ~1.5 kb upstream the transcription start sites) of C6, C7, C8A, C8G, and C9 exhibited 29 SNPs. In silico analyses indicated that most of these SNPs were located in conserved regions or even in putative transcription factor binding sites. For association studies SNPs in the promoter of C6 and C9 (C6: g.-80C>G; C9: g.-737C>T) and in the cDNA-sequence of C7 and C8A (C7: c.154A>G; C8A: c.1544C>T), respectively, were genotyped in 311 animals of commercial herds, of which relative abundance of transcripts of TCCs and lung lesion score were available. The minor allele frequencies obtained were ≥ 0.2. For C7 (p=.05) and C9 (p=.0003) association with relative transcript abundance was shown. Moreover, the analysis revealed interaction of polymorphisms and levels of expression among various complement genes. For C9 significant differences of the genotype distribution among non-affected and affected animals was evident (p=.03). Moreover, the level of expression of C7 (p=.02), C8A (p=.02), C8B (p=.09), and C8G (p=.05) differed between non-affected and affected animals .
Due to the lack of trait records related to animal health from the field we proposed to evaluate transcript levels of immune genes including TCCs obtained at slaughter as biomarkers for disease resistance. Relationships between polymorphisms of the complement genes and their expression were apparent for several loci. Here, lung lesions, which were recorded during veterinary inspection, were used as parameters of health in order to get a first insight of any links between the polymorphisms, differential basal expression of TCCs, and disease resistance. For some TCCs this three-way relationship could be shown. Subsequent eQTL analyses and relating of traits relevant to health and performance will facilitate the identification of DNA-markers – which are the ultimate biomarkers for breeding purpose – for disease resistance. The results encourage further studies including functional assays of promoter SNPs, genome wide association studies for lung lesion prevalence and transcript abundance providing eQTL, monitoring the health status and sampling within commercial herds.
This work was supported by the EC-funded FP6 Project SABRE, the BMBF/DLR (Federal Ministry of Education and Research (BMBF)/ International Bureau, Project Management Agency of the German Aerospace Center (DLR); Project-No: VNB 03/B01), the BMBF program FUGATOplus (project GeneDialog, Project-No: 0315130A) and internal funds of the FBN.
This article has been published as part of BMC Proceedings Volume 5 Supplement 4, 2011: Proceedings of the International Symposium on Animal Genomics for Animal Health (AGAH 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/5?issue=S4.
- Kumar KG, Ponsuksili S, Schellander K, Wimmers K: Molecular cloning and sequencing of porcine C5 gene and its association with immunological traits. Immunogenetics. 2004, 55: 811-817. 10.1007/s00251-003-0637-y.View ArticlePubMedGoogle Scholar
- Wimmers K, Mekchay S, Schellander K, Ponsuksili S: Molecular characterization of the pig C3 gene and its association with complement activity. Immunogenetics. 2003, 54: 714-724.PubMedGoogle Scholar
- Wimmers K, Ponsuksili S, Schellander K: BF, HP, DQB and DRB are associated with haemolytic complement activity, acute phase protein reaction and antibody response in the pig. Vet Immunol Immunopathol. 2004, 99: 215-228.View ArticlePubMedGoogle Scholar
- Phatsara C, Jennen DGJ, Ponsuksili S, Murani E, Tesfaye D, Schellander K, Wimmers K: Molecular genetic analysis of porcine mannose-binding lectin genes, MBL1 and MBL2, and their association with complement activity. International Journal of Immunogenetics. 2007, 34: 55-63. 10.1111/j.1744-313X.2007.00656.x.View ArticlePubMedGoogle Scholar
- Hardge T, Koepke K, Reissmann M, Wimmers K: Maternal influences on litter size and growth in reciprocal crossed Miniature Pigs and Durocs. Archiv Tierzucht. 1999, 42: 83-92.Google Scholar
- Wimmers K, Lipperheide C, Ponsuksili S, Schmoll F, Hardge T, Petersen B, Schellander K: Haemolytic complement activity and C3c serum concentration in pigs. Archiv Tierzucht. 1999, 42: 93-102.Google Scholar
- Green P, Falls K, Crooks S: Documentation for CRIMAP, version 2.4. 1990, Washington University, School of Medicine, St Louis, MOGoogle Scholar
- Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in simple and complex pedigrees. Bioinformatics. 2002, 18: 339-340. 10.1093/bioinformatics/18.2.339.View ArticlePubMedGoogle Scholar
- Lake SL, Blacker D, Laird NM: Family-based tests of association in the presence of linkage. American Journal of Human Genetics. 2000, 67: 1515-1525. 10.1086/316895.PubMed CentralView ArticlePubMedGoogle Scholar
- Wimmers K, Murani E, Schellander K, Ponsuksili S: QTL for traits related to humoral immune response estimated from data of a porcine F2 resource population. International Journal of Immunogenetics. 2009, 36: 141-151. 10.1111/j.1744-313X.2009.00838.x.View ArticlePubMedGoogle Scholar
- Do VAK, Ponsuksili S, Murani E, Wimmers K: IMpRH mapping of the porcine complement genes of the terminal lytic pathway. Proceedings of International Workshop on Biotechnology in Agriculture. 2006, Nong Lam University, Ho Chi Minh City, Vietnam, 32-34.Google Scholar
- Agah A, Montalto MC, Kiesecker CL, Morrissey M, Grover M, Whoolery KL, Rother RP, Stahl GL: Isolation, characterization, and cloning of porcine complement component C7. J Immunol. 2000, 165: 1059-1065.View ArticlePubMedGoogle Scholar
- Do VAK, Ponsuksili S, Murani E, Huynh TPL, Wimmers K: Molecular characterization of cDNA sequences of porcine genes encoding complement compoment of the terminal lytic pathway. Proceedings of International Workshop on Biotechnology in Agriculture. 2006, Nong Lam University, Ho Chi Minh City, Vietnam, 29-31.Google Scholar
- Schütze S, Murani E, Ponsuksili S, Wimmers K: Einflüsse von Polymorphismen in 5Â´ regulierenden Abschnitten porciner Komplementgene auf ihre Expression und Merkmale der Tiergesundheit. Schriftenreihe des FBN 16, 11th Day of the Doctoral students. Edited by: Seyfert H-M, Viereck G. 2010, Dummerstorf, Germany, 43-46.Google Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.