Document Detail

Detection of two QTL on chicken chromosome 14 for keyhole lymphet haemocyanin.
Jump to Full Text
MedLine Citation:
PMID:  22048895     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
A keyhole lymphet haemocyanin is an antigen which triggers Th1 type of immune response. A QTL for a primary immune response towards keyhole lymphet haemocyanin has been detected on chicken chromosome 14 in three populations. The results from the most recent population were inconsistent and varied depending on the applied QTL detection model. The major goal of the current study was the reanalysis of this data using a 2 QTL model. Additionally, in order to provide more accurate estimates of QTL effects and positions, epistasis between the QTL was considered as a potential important contributor to quantitative traits. Four statistical models were assumed: M1: A model assuming marginal additive effects of two QTL; M2: A model assuming marginal and epistatic additive effects of two QTL; M3: A model assuming marginal additive and dominance effects of two QTL; M4: A model assuming marginal additive and dominance effects of two QTL and all possible pairwise epistases. Two QTL with significant additive and dominance effects were detected on chicken chromosome 14 using model M3. One QTL was detected at 63 cM between MCW0123 and ROS0005, another at 76 cM between ROS0005 and MCW0225/NTN2Lsts1 (FDR = 0.0051). Modelling only additive effects resulted in a significantly worse fit. On the other hand, including epistatic effects did not improve fit significantly. The current study confirms previous reports of the QTL location on GGA14. A notable finding of this study is recognition of two closely related QTL for a keyhole lymphet haemocyanin response at the distal part of chicken chromosome 14.
Authors:
Maria Siwek; Joanna Szyda; Anna Sławińska; Marek Bednarczyk
Related Documents :
20552275 - Modeling the complex dynamics of derecruitment in the lung.
12897805 - Adaptation in the respiratory control system.
11247925 - A method of reconstruction of clinical gas-analyzer signals corrupted by positive-press...
15587645 - Mutual information based ct registration of the lung at exhale and inhale breathing sta...
7471385 - New approach to cystic fibrosis diagnosis by use of an improved sweat-induction/collect...
24860035 - Design, development and clinical evaluation of the electronic mobility cane for vision ...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2011-11-03
Journal Detail:
Title:  Journal of applied genetics     Volume:  53     ISSN:  2190-3883     ISO Abbreviation:  J. Appl. Genet.     Publication Date:  2012 Feb 
Date Detail:
Created Date:  2012-01-25     Completed Date:  2012-05-14     Revised Date:  2013-11-14    
Medline Journal Info:
Nlm Unique ID:  9514582     Medline TA:  J Appl Genet     Country:  England    
Other Details:
Languages:  eng     Pagination:  115-9     Citation Subset:  IM    
Affiliation:
Department of Animal Biotechnology, University of Technology and Life Sciences, Mazowiecka 28, Bydgoszcz, Poland. siwek@utp.edu.pl
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Animals
Antibody Formation
Chickens / genetics*,  immunology
Chromosomes / genetics*
Crosses, Genetic
DNA / analysis,  genetics
Epistasis, Genetic
Female
Hemocyanin / genetics*,  immunology
Linear Models
Male
Microsatellite Repeats
Phenotype
Quantitative Trait Loci*
Chemical
Reg. No./Substance:
9007-49-2/DNA; 9013-72-3/Hemocyanin; FV4Y0JO2CX/keyhole-limpet hemocyanin
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): J Appl Genet
ISSN: 1234-1983
ISSN: 2190-3883
Publisher: Springer-Verlag, Berlin/Heidelberg
Article Information
Download PDF
© The Author(s) 2011
Received Day: 11 Month: 7 Year: 2011
Revision Received Day: 10 Month: 10 Year: 2011
Accepted Day: 10 Month: 10 Year: 2011
Electronic publication date: Day: 3 Month: 11 Year: 2011
pmc-release publication date: Day: 3 Month: 11 Year: 2011
Print publication date: Month: 2 Year: 2012
Volume: 53 Issue: 1
First Page: 115 Last Page: 119
ID: 3265721
PubMed Id: 22048895
Publisher Id: 74
DOI: 10.1007/s13353-011-0074-9

Detection of two QTL on chicken chromosome 14 for keyhole lymphet heamocyanin
Maria Siwek1 Address: +48-523749767 +48-523228158 siwek@utp.edu.pl
Joanna Szyda2
Anna Sławińska1
Marek Bednarczyk1
1Department of Animal Biotechnology, University of Technology and Life Sciences, Mazowiecka 28, 85-225 Bydgoszcz, Poland
2Department of Animal Genetics, Wrocław University of Environmental and Life Sciences, Kożuchowska 7, 51-631 Wroclaw, Poland

Introduction

The identification of quantitative trait loci (QTL) usually follows the same pathway, which includes QTL detection in the whole genome scan, fine mapping and narrowing down the QTL region, confirmation of QTL location in the independent population and finally functional and positional candidate gene approach. In chicken, there is a long list of traits for which QTL have been detected (Abasht et al. 2006). Among them 83 are related to immune traits. So far, only a few of these QTL were confirmed in an independent population or fine mapped (Kim et al. 2006; Calenge et al. 2009). The current study is a follow up analysis of the QTL validated on a GGA14 for an immune response towards a non pathogenic antigen (Siwek et al.2010).

In the current study we use a novel protein named keyhole lymphet heamocyanin (KLH). KLH is a copper-containing high-molecular-weight protein. This protein is involved in the Th2 immune pathways in mammals (Bliss et al. 1996). This antigen has been used in several studies of chicken immunity (Minozzi et al. 2008; Minozzi et al.2007; Hangalapura et al.2004) and in QTL studies based on microsatellite markers (Siwek et al.2003). This non pathogenic antigen is a potent immunostimulator used as hapten carrier to assess the magnitude of immune responses (Marcia et al. 2006). A QTL for a primary antibody response towards KLH has been first detected on GGA14 in two independent experimental populations (Siwek et al.2003). The first confirmation of that QTL was performed on a feather pecking cross (Rodenburg and Koene 2001; Buitenhuis et al.2003; Siwek et al.2003). In the current experimental population, which is a White Leghorn and Green legged partidgelike cross (WL/ZK), a QTL on GGA14 was detected, but its most probable location varied depending on the applied QTL detection model (LDLA or Half sib; Siwek et al. 2010).

Because of the inconsistency of the WL/ZK results we were interested in checking whether a single or multiple QTL reside on GGA14. Consequently the major goal of the current study was the reanalysis of the data using a 2 QTL model. Additionally, in order to provide more accurate estimates of QTL effects (Zeng et al.2005) and positions, epistasis between the QTL was considered as a potentially important contributor to quantitative traits (Mackay 2001).


Materials and methods
Experimental population

The experimental population was created from a cross between two breeds of chickens: White Leghorn (WL), which is a commercial layer and Green legged partidgelike (ZK), which is a Polish native chicken breed. Animals were kept on the farm at the University of Life Sciences in Lublin. We used floor system to keep the birds. Entire experimental population was vaccinated against Salmonella, Gumboro disease, bronchitis, Bourse Fabricius disease and encephalomyelitis. The reproduction of the population was as follows: five F0 roosters (WL) were mated to six – eight hens (ZK). From that cross of the parental breeds we obtained a F1 generation consisting of 115 individuals. Subsequently six F1 roosters were mated to six – nine hens each. In the F2 generation we obtained 506 individuals from 35 full sib families. Every time fertilized eggs were collected and incubated in an automatic incubator. The selection of the F1 parents was based on the KLH antibody titre in such a way that a high responding rooster was mated to low responding hens or the other way around.

Genotypic data

Genomic DNA was isolated from blood using a commercial kit MasterPure DNA Purification Kit for Blood (Epicenter ® biotechnologies, MG 711100). PCR reaction was used to amplify microsatellite markers. Seven microsatellite markers located along 77 cM of GGA14 were used for genotyping: MCW0296 (at 0 cM), ADL0200 (at 16 cM), ADL0118 (at 20 cM), MCW0136 (at 20 cM), MCW0123 (at 45 cM), ROS0005 (at 71 cM), and MCW0225/NTN2Lsts1 (at 77 cM). The reaction was carried out in a volume of 10 microlitres. Each reaction consisted of a mixture containing 75 mM Tris-HCl (pH 8.8), 20 mM (NH4)2SO4, 0.01% Tween 20, 5.0 nmol TMACl, 2.0 mM MgCl2, 20 ng DNA matrix, 0.25 units of recombinant Taq polymerase (Fermentas, Vilnius, Lithuania), and a pair of specific oligonucleotides for each marker, in an amount from 0.2 to 0.75 pmol each. After PCR reaction samples were loaded on a 6% denaturating polyacrylamide gel and subjected to electrophoresis. Fragment length analysis was based on a laser fluorescence detection, using an automatic plate sequenator Li-Cor 4200, together with software designed for genotyping of microsatellite markers Saga GT.

Phenotypic data

Total antibody responses to KLH were measured in individual plasma samples obtained at 7 days after s.c. immunization with 1 mg KLH in 1 ml PBS (pH 7.2) at 12 wk of age. Antibody titres to KLH of all birds were measured by an indirect ELISA as described by Sijben et al. (2000). Titres were expressed as the log2 values of the highest dilution giving a positive reaction.

Statistical analysis

Four mixed linear models with 2 QTL effects were applied to the data. All the models were evaluated for putative QTL positions every 1 cM in a 2-dimentional grid.

M1: A model assuming marginal additive effects of two QTL:

[Formula ID: Equa]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbf{y}} = {\mathbf{Z\alpha }} + {{\mathbf{X}}_{{{a_1}}}}{{\mathbf{q}}_{{{a_1}}}} + {{\mathbf{X}}_{{{a_2}}}}{{\mathbf{q}}_{{{a_2}}}} + {\mathbf{e}}, $$\end{document}]
where y contains trait values of F2 individuals, [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbf{\alpha }}\sim N\left( {0,{\mathbf{A}}\sigma_{\alpha }^2} \right) $$\end{document}] is a vector of random polygenic effects of F2 animals with A representing an additive polygenic relationship matrix among individuals, and [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sigma_{\alpha }^2 $$\end{document}] being a component of the total additive genetic variance attributed to polygenes, [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{q}}_{{{a_i}}}} $$\end{document}] represents a fixed additive effect of ith QTL, and [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbf{e}}\sim N\left( {0,{\mathbf{I}}\sigma_e^2} \right) $$\end{document}] is a residual term, [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{X}}_{{{a_i}}}} $$\end{document}] is the matrix with coefficients for an additive effect of ith QTL calculated based on individuals' marker genotypes for a putative position of ith QTL using the QTL Express software (Seaton et al.2002), Z is a design matrix for α, which is independent of the putative QTL positions considered.

M2: A model assuming marginal and epistatic additive effects of two QTL:

[Formula ID: Equb]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbf{y}} = {\mathbf{Z\alpha }} + {{\mathbf{X}}_{{{a_1}}}}{{\mathbf{q}}_{{{a_1}}}} + {{\mathbf{X}}_{{{a_2}}}}{{\mathbf{q}}_{{{a_2}}}} + {{\mathbf{X}}_{{{a_1}{a_2}}}}{{\mathbf{q}}_{{{a_1}{a_2}}}} + {\mathbf{e}}, $$\end{document}]
where [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{q}}_{{{a_1}{a_2}}}} $$\end{document}] represents an additive-by-additive epistatic effect between both QTL with the corresponding coefficient matrix [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{X}}_{{{a_1}{a_2}}}} $$\end{document}] .

M3: A model assuming marginal additive and dominance effects of two QTL:

[Formula ID: Equc]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbf{y}} = {\mathbf{Z\alpha }} + {{\mathbf{X}}_{{{a_1}}}}{{\mathbf{q}}_{{{a_1}}}} + {{\mathbf{X}}_{{{a_2}}}}{{\mathbf{q}}_{{{a_2}}}} + {{\mathbf{X}}_{{{d_1}}}}{{\mathbf{q}}_{{{d_1}}}} + {{\mathbf{X}}_{{{d_2}}}}{{\mathbf{q}}_{{{d_2}}}} + {\mathbf{e}}, $$\end{document}]
where [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{q}}_{{{d_i}}}} $$\end{document}] represents a vector of fixed dominance effect of ith QTL with the corresponding coefficient matrix [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{X}}_{{{d_i}}}} $$\end{document}] .

M4: A model assuming marginal additive and dominance effects of two QTL and all possible pairwise epistases:

[Formula ID: Equd]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\mathbf{y}} = {\mathbf{Z\alpha }} + {{\mathbf{X}}_{{{a_1}}}}{{\mathbf{q}}_{{{a_1}}}} + {{\mathbf{X}}_{{{a_2}}}}{{\mathbf{q}}_{{{a_2}}}} + {{\mathbf{X}}_{{{d_1}}}}{{\mathbf{q}}_{{{d_1}}}} + {{\mathbf{X}}_{{{d_2}}}}{{\mathbf{q}}_{{{d_2}}}} + {{\mathbf{X}}_{{{a_1}{a_2}}}}{{\mathbf{q}}_{{{a_1}{a_2}}}} + {{\mathbf{X}}_{{{a_1}{d_2}}}}{{\mathbf{q}}_{{{a_1}{d_2}}}} + {{\mathbf{X}}_{{{d_1}{a_2}}}}{{\mathbf{q}}_{{{d_1}{a_2}}}} + {{\mathbf{X}}_{{{d_1}{d_2}}}}{{\mathbf{q}}_{{{d_1}{d_2}}}} + {\mathbf{e}}, $$\end{document}]
where [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{q}}_{{{a_1}{d_2}}}} $$\end{document}] , [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{q}}_{{{d_1}{a_2}}}} $$\end{document}] , and [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{q}}_{{{d_1}{d_2}}}} $$\end{document}] represent respectively an additive-by-dominance, dominance-by-additive, and dominance-by-dominance epistatic effects between both QTL with the corresponding coefficient matrices [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{X}}_{{{a_1}{d_2}}}} $$\end{document}] , [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{X}}_{{{d_1}{a_2}}}} $$\end{document}] , and [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {{\mathbf{X}}_{{{d_1}{d_2}}}} $$\end{document}] . Note that all the coefficient matrices for epistatic effects result from the multiplication of the corresponding coefficient matrices for marginal QTL effects.

Model parameters were estimated using Maximum Likelihood, under the assumption that the variance component underlying a random polygenic effect are known (i.e. this variance component was not estimated). Models were compared using the likelihood ratio test [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \lambda = - 2\left[ {\ln L\left( {{{\hat{\beta }}_0}} \right) - \ln L\left( {{{\hat{\beta }}_1}} \right)} \right] $$\end{document}] , where [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ L\left( {{{\hat{\beta }}_0}} \right) $$\end{document}] and [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ L\left( {{{\hat{\beta }}_1}} \right) $$\end{document}] represent values of the likelihood function for a more parsimonious and a less parsimonious model respectively. The nominal P values of the test were calculated based on the χ2 distribution with the number of degrees of freedom corresponding to a difference in the number of parameters between the two models. Since 4466 pairwise location for the two QTL were tested, the false discovery rate (FDR; Benjamini and Hochberg 1995) was calculated, based on the nominal P values, in order to circumvent the multiple testing problem. The FDR adjusted nominal P values (PFDR) were calculated as follows: [\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ P_{{4466}}^{{FDR}} = {P_{{4466}}},\,\forall i \in \left\{ {2,. \cdots, 4465} \right\}{:}P_{{4466 - i}}^{{FDR}} = \min \left( {P_{{4466 - i + 1}}^{{FDR}},\frac{{4466}}{{4466 - i}}{P_{{4466 - i}}}} \right), P_1^{{FDR}} = \min \left( {P_2^{{FDR}},4466{P_1}} \right) $$\end{document}] , where index i represents the order of nominal P values from the smallest (i = 1) to the largest (i = 4 466). Estimation of model parameters and P value adjustment were calculated using the SAS software (SAS Institute Inc.). PROC MIXED procedure was used for the estimation of model parameters and PROC MULTTEST procedure was used for P value adjustment.


Results and discussion

Two QTL with significant additive and dominance effects at 63 cM between MCW0123 and ROS0005 and 76 cM between ROS0005 and MCW0225/NTN2Lsts1 (FDR = 0.0051; Fig. 1) were detected on chicken chromosome 14 (GGA14). Modelling only additive effects resulted in a significantly worse fit. On the other hand, including epistatic effect did not improve fit significantly.

Non-random LD in chickens, which extends over short distances reported by Andreescu et al. 2007, justifies the choice of linkage mapping over the association analysis as it was done in the presented manuscript.

According to the Animal QTL data base chicken chromosome 14 harbours altogether 28 QTL for various traits (http://www.animalgenome.org/cgi-bin/QTLdb/GG/index). We were particularly interested in a QTL for a primary antibody response towards KLH. This trait has been studied in three experimental populations: SRBC population - which is a cross between two chicken lines divergently selected for high (H) and low (L) primary antibody response towards sheep red blood cells (SRBC) (Siwek et al. 2003); F/P population - a cross between two chicken lines which differed in their propensity for feather pecking (Rodenburg and Koene 2001) and the White Leghorn/ ZK cross (Siwek et al. 2010). In all three populations (SRBC, F/P, WL/ZK) all methods used for the analysis were the same, individuals received the same treatment which allows comparison of the results from three experiments. The only difference was the time of the KLH injection. The F/P population birds were injected at 36 weeks of age, and at 12 weeks of age in the SRBC and WL/ZK populations. At first a QTL for a primary antibody response towards KLH was detected in the SRBC population in the whole genome scan approach. The evidence of the QTL for a KLH antigen was provided by two models: paternal half – sib model, and a line cross model. Both models pointed towards the same location between microsatellite markers: MCW0123 and MCW0225 (Siwek et al. 2003). First confirmation of this QTL was performed in the F/P population. Paternal half – sib model, which was applied to analyse the dataset, pointed a QTL location towards the proximal part of GGA14 in the marker bracket: MCW0296 - ADL0118. Despite that the time of the KLH stimulation differed in both populations, the confidence interval estimated for a QTL position on GGA14 was large, the marker density in both cases (SRBC population and F/P population) was low which drove us to the conclusion that we dealt with the same QTL in both populations and slightly different peak position is due to the diverse time of KLH injection or differences in the distribution of recombination events along GGA14. The second confirmation of the QTL for a primary antibody response towards KLH on GGA14 was performed in the experimental cross of White Leghorn and Green legged partidgelike (Siwek et al. 2010). In the WL/ZK reference population a QTL for a primary antibody response towards KLH was detected in the marker bracket: MCW0136 – MCW0125 with a maternal half – sib analysis model. In fact there were two QTLs detected on GGA14 with initial half sib linkage analysis models. One QTL in the middle part of GGA14 detected with maternal half sib QTL model, presented in Fig. 1A (Siwek et al. 2010) and a second one at the distal part of GGA14 detected with paternal half sib QTL model. The second QTL was just below the significance threshold (Table 3, Siwek et al. 2010). Another model (LDLA) applied to this data set moved the QTL peak towards the end of the chromosome, to the marker bracket: ROS0005 – MCW0225. Multiple QTL analysis confirmed location of the QTL pointed by LDLA analysis and separated QTL location detected with the maternal – half sib model as a second independent QTL.

Results of the candidate genes approach have been reported elsewhere (Siwek et al. 2010). In short, three candidate genes where selected in the QTL region: AXIN1 (axis inhibition protein 1), TRAP (TNF receptor-associated protein 1), and MAPK8IP3 (mitogen-activated protein kinase 8 interacting protein 3). Given the current findings on two QTL located at the distal part of GGA14 these three genes are distributed over the two QTL marker brackets. AXIN1 and TRAP being located in the first QTL region and MAPK8IP3 being located in the second QTL region. These genes have been genotyped with SNP located in gene exons. Statistical analyses of single SNP associations were performed pointing out two SNPs within the AXIN1 gene as significantly associated with the primary antibody response to KLH (Siwek et al. 2010). Even though the presented candidate gene approach does not provide solid evidence for a second gene related with a trait of interest in the other QTL region, there is a long list of hypothetical proteins and undefined sequences in the distal part of chicken chromosome 14 which gives room for new detections.

In summary, we hereby present results of a multi QTL analysis aimed at the identification of QTL involved in the primary antibody response to KLH. Current study confirms previous reports of the QTL location on GGA14. A notable finding of this study is recognition of two closely related QTL for KLH response at the distal part of chicken chromosome 14.


Acknowledgement

The study was financed by the Polish State Committee for Scientific Research, grants no. N N311 450434 and N N311 558640.

Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.


References
Andreescu C,Avendano S,Brown SR,Hassen A,Lamont SJ,Dekkers JCM. Linkage disequilibrium in related breeding lines of chickensGeneticsYear: 20071772161216910.1534/genetics.107.08220617947400
Abasht B, Dekkers JC, Lamont SJ (2006) Review of quantitative trait loci identified in the chicken. Poult Sci 85:2079–2096. Erratum: Poult. Sci. 86:206
Benjamini Y,Hochberg Y. Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple TestingJournal of the Royal Statistical Society. BYear: 199557289300
Bliss J,Cleave V,Murray K,Wiencis A,Ketchum M,Maylor R,Haire T,Resmini C,Abbas AK,Wolf SF. IL-12, as an adjuvant, promotes a T helper 1 cell, but does not suppress a T helper 2 cell recall responseJ ImmunolYear: 19961568878648558014
Buitenhuis AJ,Rodenburg TB,Hierden Y,Siwek M,Cornelissen SJB,Nieuwland MGB,Crooijmans RPMA,Groenen MAM,Koene P,Korte SM,Bovenhuis H,Poel JJ. Mapping Quantitative Trait Loci affecting feather pecking behaviour and stress response in laying hensPoult SciYear: 2003821215122212943291
Calenge F,Lecerf F,Demars J,Feve K,Vignoles F,Pitel F,Vignal A,Velge P,Sellier N,Beaumont C. QTL for resistance to Salmonella carrier state confirmed in both experimental and commercial chicken linesAnim GenetYear: 20094059059710.1111/j.1365-2052.2009.01884.x19422366
Hangalapura BN,Nieuwland MG,Vries RG,Brand GH,Kemp B,Parmentier HK. Durations of cold stress modulates overall immunity of chicken lines divergently selected for antibody responsesPoult SciYear: 20048376577515141834
Kim ES,Hong YM,Min W,Lillehoj HS. Fine-mapping of coccidia-resistant quantitative trait loci in chickensPoult SciYear: 2006852028203017032840
Marcia L,Delacre M,Abboud G. Impairment of dendritic cell functionality and steady-state number in obese miceJ. ImmunolYear: 20061775997600617056524
Minozzi G,Parmentier HK,Nieuwland MG,Bed'hom B,Minvielle F,Gourichon D,Pinard-van der Laan MH. Antibody responses to keyhole limpet hemocyanin, lipopolysaccharide, and Newcastle Disease virus vaccine in F2 and backcrosses of white Leghorn lines selected for two different immune response traitsPoult SciYear: 2007861316132217575177
Minozzi G, Parmentier HK, Mignon-Grasteau S, Nieuwland MGB, Bed'hom B, Gourichon D, Minvielle F, Pinard-van der Laan MH (2008) Correlated effects of selection for immunity in White Leghorn chicken lines on natural antibodies and specific antibody responses to KLH and M. butyricum. BMC Genetics 9, 5
Rodenburg B, Koene P (2001) Comparison of individual and social feather pecking tests in two lines of White Leghorn Laying hens. In: Proceedings of the Six European Symposium on Poultry Welfare, Zollikofen, Switzerland, pp 250-253
Seaton G,Haley CS,Knott SA,Kearsey M,Visscher PM. QTL Express: mapping quantitative trait loci in simple and complex pedigreesBioinformaticsYear: 20021833934010.1093/bioinformatics/18.2.33911847090
Sijben JW,Groot H,Nieuwland MGB,Schrama JW,Parmentier HK. Dietary linoleic acid divergently affects immune responsiveness of growing layer hensPoult SciYear: 2000791106111510947178
Siwek M,Buitenhuis AJ,Cornelissen SJB,Nieuwland MGB,Bovenhuis H,Crooijmans RPMA,Groenen MAM,Vries-Reilingh G,Parmentier HK,Poel JJ. Detection of different QTL for antibody response to Keyhole Lymphet Haemocynain and Mycobacterium butiricum in two unrelated populations of laying hensPoult SciYear: 2003821845184214717541
Siwek M,Slawinska A,Nieuwland M,Witkowski A,Zięba G,Minozzi G,Knol EF,Bednarczyk M. A quantitative trait locus for a primary antibody response to keyhole limpet hemocyanin on chicken chromosome 14—Confirmation and candidate gene approachPoult SciYear: 2010891850185710.3382/ps.2010-0075520709969
SAS Institute Inc., SAS 9.1.3 Help and Documentation, Cary, NC: SAS Institute Inc., 2000–2004
Zeng ZB,Wang T,Zou W. Modeling Quantitative Trait Loci and Interpretation of ModelsGeneticsYear: 20051691711172510.1534/genetics.104.03585715654105

Figures

[Figure ID: Fig1]
Fig. 1 

Likelihood ratio test profile for a two QTL model. Test corresponds to the hypothesis of additive vs. additive and dominance effects (M1 vs. M3) of the QTL for a primary immune response towards keyhole lymphet heamocyanin. CM1 and CM2 - distance in centiMorgan; L - likelihood ratio test



Article Categories:
  • Animal Genetics • Original Paper

Keywords: Keywords Chicken, Epistatic effect, Immune trait, Multiple quantitative trait loci.

Previous Document:  Climate change experiments in temperate grasslands: synthesis and future directions.
Next Document:  FGFR2 mutation confers a less drastic gain of function in mesenchymal stem cells than in fibroblasts...