Detection of different quantitative trait loci for antibody responses to keyhole lympet hemocyanin and Mycobacterium butyricum in two unrelated populations of laying hens

SNP prioritization in targeted sequencing data associated with humoral immune responses in chicken

Our study aimed to identify single nucleotide polymorphisms (SNPs) with a significant impact on the innate immunity represented by antibody response against lipopolysaccharide (LPS) and lipoteichoid acid (LTA) and the adaptive immune response represented toward keyhole limpet hemocyanin (KLH) using the SNP prioritization method. Data set consisted of 288 F2 experimental individuals, created by crossing Green-legged Partridgelike and White Leghorn. The analyzed SNPs were located within 24 short genomic regions of GGA1, GGA2, GGA3, GGA4, GGA9, GGA10, GGA14, GGA18, and GGZ, pre-targeted based on literature references and database information. For the specific antibody response toward KLH at d 0 the most highly prioritized SNP for additive and dominance effects were located on GGA2 in the 3’UTR of MYD88. For the response at d 7, the most highly prioritized SNP pointed at the 3’UTR of MYD88, but potential causal additive variants were located within ADIPOQ and one in PROCR. The highest priority for additive and dominance effects in the antibody response toward lipoteichoic acid at d 0 was attributed to the same SNP, located on GGA2 in the 3’UTR region of MYD88. Two SNPs among the top-10 for additive effect were located in the exon of NOCT. SNPs selected for their additive effect on antibody response toward lipopolysaccharide at d 0 marked 3 genes – NOCT, MYD88, and SNX8, while SNPs selected for their dominance effect marked – NOCT, ADIPOQ, and MYD88. The top-10 variants identified in our study were located in different functional parts of the genome. In the context of causality three groups can be distinguished: variants located in exons of protein coding genes (ADIPOQ, NOCT, PROCR, SNX8), variants within exons of non-coding transcripts, and variants located in genes’ UTR regions. Variants from the first group influence protein structure and variants from both latter groups’ exhibit regulatory roles on DNA (UTR) or RNA (lncRNA).

Sensitivity analysis based on the random forest machine learning algorithm identifies candidate genes for regulation of innate and adaptive immune response of chicken

Two categories of immune responses—innate and adaptive immunity—have both polygenic backgrounds and a significant environmental component. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chickens using machine learning–based sensitivity analysis for single-nucleotide polymorphisms (SNPs) located in candidate genes defined in quantitative trait loci regions. Here the adaptive immunity is represented by the specific antibody response toward keyhole limpet hemocyanin (KLH), whereas the innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). The analysis consisted of 3 basic steps: an identification of candidate SNPs via feature selection, an optimisation of the feature set using recursive feature elimination, and finally a gene-level sensitivity analysis for final selection of models. The predictive model based on 5 genes (MAPK8IP3 CRLF3, UNC13D, ILR9, and PRCKB) explains 14.9% of variance for KLH adaptive response. The models obtained for LTA and LPS use more genes and have lower predictive power, explaining respectively 7.8 and 4.5% of total variance. In comparison, the linear models built on genes identified by a standard statistical analysis explain 1.5, 0.5, and 0.3% of variance for KLH, LTA, and LPS response, respectively. The present study shows that machine learning methods applied to systems with a complex interaction network can discover phenotype-genotype associations with much higher sensitivity than traditional statistical models. It adds contribution to evidence suggesting a role of MAPK8IP3 in the adaptive immune response. It also indicates that CRLF3 is involved in this process as well. Both findings need additional verification.

The growing feather as a dermal test site: Comparison of leukocyte profiles during the response to Mycobacterium butyricum in growing feathers, wattles, and wing webs

Using the response to Mycobacterium butyricum as the test-immune response, the main goal of this study was to demonstrate the suitability of the growing feather (GF) as a dermal test site and window into in vivo cellular/tissue responses (US-Patent 8,216,551). Using M. butyricum immunized chickens, the specific objectives were to: 1) compare the leukocyte infiltration response to intra-dermally injected M. butyricum in GF, wattles, and wing webs; 2) use GF as the test site to monitor leukocyte response profiles to recall antigen in the same individuals; and 3) gain new knowledge regarding the local response to M. butyricum in chickens. For objective 1, chickens were euthanized for tissue collection at 4 to 6, 24, 48, and 72 h after intra-dermal antigen injection. Leukocyte infiltration profiles were determined using immunochemical and conventional histology. Data from this study established the similarities between the cellular response in GF, wattles, and wing webs and uncovered many advantages of working with GF. For objective 2, antigen was injected into multiple GF per individual. GF were collected before and at 0.25, 1, 2, 3, and 7 d post injection and processed for cell population analysis by flow cytometry. Advantages of the approach used in objective 2 included a technically easier, more comprehensive, and more objective leukocyte profile analysis; same-day data acquisition; and, most importantly, easy, minimally invasive sample collection from the same individual throughout the study. Both studies contributed new knowledge regarding the local cutaneous response to M. butyricum in M. butyricum immunized chickens and confirmed the cell-mediated nature of the immune response to M. butyricum (e.g., elevated levels [P < 0.05] of T cells [CD4+ and CD8+], macrophages and MHC class II+-cells on days one to 3 post injection in M. butyricum- compared to PBS-injected tissues). The use of GF as an "in vivo test tube" to monitor local innate and adaptive immune activities will find direct application in vaccine development, as well as in the assessment and optimization of immune system development and function in poultry.

The identification of loci for immune traits in chickens using a genome-wide association study

The genetic improvement of disease resistance in poultry continues to be a challenge. To identify candidate genes and loci responsible for these traits, genome-wide association studies using the chicken 60k high density single nucleotide polymorphism (SNP) array for six immune traits, total serum immunoglobulin Y (IgY) level, numbers of, and the ratio of heterophils and lymphocytes, and antibody responses against Avian Influenza Virus (AIV) and Sheep Red Blood Cell (SRBC), were performed. RT-qPCR was used to quantify the relative expression of the identified candidate genes. Nine significantly associated SNPs (P < 2.81E-06) and 30 SNPs reaching the suggestively significant level (P < 5.62E-05) were identified. Five of the 10 SNPs that were suggestively associated with the antibody response to SRBC were located within or close to previously reported QTL regions. Fifteen SNPs reached a suggestive significance level for AIV antibody titer and seven were found on the sex chromosome Z. Seven suggestive markers involving five different SNPs were identified for the numbers of heterophils and lymphocytes, and the heterophil/lymphocyte ratio. Nine significant SNPs, all on chromosome 16, were significantly associated with serum total IgY concentration, and the five most significant were located within a narrow region spanning 6.4kb to 253.4kb (P = 1.20E-14 to 5.33E-08). After testing expression of five candidate genes (IL4I1, CD1b, GNB2L1, TRIM27 and ZNF692) located in this region, changes in IL4I1, CD1b transcripts were consistent with the concentrations of IgY, while abundances of TRIM27 and ZNF692 showed reciprocal changes to those of IgY concentrations. This study has revealed 39 SNPs associated with six immune traits (total serum IgY level, numbers of, and the ratio of heterophils and lymphocytes, and antibody responses against AIV and SRBC) in Beijing-You chickens. The narrow region spanning 247kb on chromosome 16 is an important QTL for serum total IgY concentration. Five candidate genes related to IgY level validated here are novel and may play critical roles in the modulation of immune responses. Potentially useful candidate SNPs for marker-assisted selection for disease resistance are identified. It is highly likely that these candidate genes play roles in various aspects of the immune response in chickens.

Identification of candidate genes and mutations in QTL regions for immune responses in chicken

There are two categories of immune responses – innate and adaptive immunity – both having polygenic backgrounds and a significant environmental component. In our study, adaptive immunity was represented by the specific antibody response toward keyhole limpet hemocyanin (KLH); innate immunity was represented by natural antibodies toward lipopolysaccharide (LPS) and lipoteichoic acid (LTA). Defining genetic bases of immune responses leads from defining quantitative trait loci (QTL) toward a single mutation responsible for variation in the phenotypic trait. The goal of the reported study was to define candidate genes and mutations for the immune traits of interest in chicken by performing an association study of SNPs located in candidate genes defined in QTL regions. Candidate genes and SNPs in QTL regions were selected in silico. SNP association was based on a custom SNP panel, GoldenGate genotyping assay (Illumina) and two statistical models: random mixed model and CAR score. The most significant SNP for immune response toward KLH was located in the JMJD6 gene located on GGA18. Four SNPs in candidate genes FOXJ1 (GGA18), EPHB1 (GGA9), PTGER4 (GGAZ) and PRKCB (GGA14) showed association with natural antibodies for LPS. A single SNP in ITGB4 (GGA18) was associated with natural antibodies for LTA. All associated SNPs mentioned above showed additive effects.

A genome-wide association study identifies major loci affecting the immune response against infectious bronchitis virus in chicken

•

The genetic basis of host responses to infectious bronchitis virus is unclear.

•

We detected 20 significant markers for the antibody response to infectious bronchitis virus in chicken.

•

Loci on chromosomes 1 and 5 explained 12% and 13% of phenotypic variation.

•

The host immune response cluster had 13 beta-defensin and interleukin-17F genes.

•

Our results will contribute to the control of outbreaks of infectious bronchitis.

Coronaviruses are a hot research topic because they can cause severe diseases in humans and animals. Infectious bronchitis virus (IBV), belonging to gamma-coronaviruses, causes a highly infectious respiratory viral disease and can result in catastrophic economic losses to the poultry industry worldwide. Unfortunately, the genetic basis of the host immune responses against IBV is poorly understood. In the present study, the antibody levels against IBV post-immunization were measured by an enzyme-linked immunosorbent assay in the serum of 511 individuals from a commercial chicken (Gallus gallus) population. A genome-wide association study using 43,211 single nucleotide polymorphism markers was performed to identify the major loci affecting the immune response against IBV. This study detected 20 significant (P < 1.16 × 10⁻⁶) effect single nucleotide polymorphisms for the antibody level against IBV. These single nucleotide polymorphisms were distributed on five chicken chromosomes (GGA), involving GGA1, GGA3, GGA5, GGA8, and GGA9. The genes in the 1-Mb windows surrounding each single nucleotide polymorphism with significant effect for the antibody level against IBV were associated with many biological processes or pathways related to immunity, such as the defense response and mTOR signaling pathway. A genomic region containing a cluster of 13 beta-defensin (GAL1–13) and interleukin-17F genes on GGA3 probably plays an important role in the immune response against IBV. In addition, the major loci significantly associated with the antibody level against IBV on GGA1 and GGA5 could explain about 12% and 13% of the phenotypic variation, respectively. This study suggested that the chicken genome has several important loci affecting the immune response against IBV, and increases our knowledge of how to control outbreaks of infectious bronchitis.

Meta - and combined - QTL analysis of different experiments on immune traits in chickens

Meta and/or combined QTL analysis from multiple studies can improve quantitative trait loci (QTL) position estimates compared to the individual experiments. Hereby we present results of a meta-analysis of QTL on chicken chromosome 9, 14 and 18 using data from three separate experiments and joint QTL analysis for chromosome 14 and 18. Meta QTL analysis uses information from multiple QTLs studies. Joint QTL analysis is based on combining raw data from different QTL experimental populations. QTLs under the study were related to specific antibody response to keyhole lymphet hemocyanin (KLH), and natural antibodies to environmental antigens, lipopolisaccharide (LPS) and lipoteichoic acid (LTA). Meta QTL analysis resulted in narrowing down the confidence interval for two QTLs on GGA14. The first one for natural antibodies against LTA and the second one for specific antibody response toward KLH. Also, a confidence interval of a QTL for natural antibodies against LPS located on GGA18 was narrowed down. Combined QTL analysis was successful for two QTLs: for specific antibody response toward KLH on GGA14, and for natural antibodies against LPS on GGA18. The greatest statistical power for QTL detection in joint analysis was achieved when raw data from segregating half–sib families from different populations under the study was used.

Detection of two QTL on chicken chromosome 14 for keyhole lymphet haemocyanin

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.

Quantitative trait loci associated with the humoral innate immune response in chickens were confirmed in a cross between Green-Legged Partridgelike and White Leghorn

Natural antibodies (NA) create a crucial barrier at the initial steps of the innate humoral immune response. The main role of NA in the defense system is to bind the pathogens at early stages of infection. Different pathogens are recognized by the presence of highly conserved antigen determinant [e.g., lipopolysaccharide (LPS) in gram-negative bacteria or lipoteichoic acid (LTA) in gram-positive bacteria]. In chickens, a different genetic background of NA binds LPS and LTA antigens, encoded by different QTL. The main objective of this work was to confirm known QTL associated with LPS and LTA NA. For this purpose a chicken reference population was created by crossing 2 breeds: a commercial layer, White Leghorn, and a Polish indigenous chicken, Green-Legged Partridgelike. The chromosomal regions analyzed harbored to GGA3, GGA5, GGA6, GGA8, GGA9, GGA10, GGA14, GGA15, GGA18, and GGAZ. The data collected consisted of the NA titers binding LPS and LTA (determined by ELISA at 12 wk of age) as well as the genotypes (30 short tandem repeat markers; average of 3 markers/chromosome, collected for generations F(0), F(1), and F(2)). The analyses were performed with 3 statistical models (paternal and maternal half-sib, line cross, and linkage analysis and linkage disequilibrium) implemented in GridQTL software (http://www.gridqtl.org.uk/). The QTL study of humoral innate immune response traits resulted in the confirmation of 3 QTL associated with NA titers binding LPS (located on GGA9, GGA18, and GGAZ) and 2 QTL associated with NA titers binding LTA (located on GGA5 and GGA14). A set of candidate genes within the regions of the validated QTL has been proposed.

New QTL for resistance to Salmonella carrier-state identified on fowl microchromosomes

Chicken's ability to carry Salmonella without displaying disease symptoms leads to an invisible propagation of Salmonella in poultry stocks. Using chicken lines more resistant to carrier state could improve both animal health and food safety. Previous studies identified several QTL for resistance to carrier state. To improve genome coverage and QTL detection power we produced a new set of 480 informative SNP markers and genotyped a larger number of animals. Ten additional microchromosomes could be covered when compared with previous studies. These new data led to the identification of 18 QTL significant at the chromosome-wide level. The only QTL significant at the genome-wide level were identified on microchromosomes 14 and 22 and have never been identified previously. Using a higher number of animals improved the power and the precision of QTL detection. Some of the QTL newly identified are located close to candidate genes or microsatellite markers previously identified for their involvement in the genetic control of resistance to Salmonella, which confirms their interest for selection purposes.

A quantitative trait locus for a primary antibody response to keyhole limpet hemocyanin on chicken chromosome 14--confirmation and candidate gene approach

A QTL involved in the primary antibody response toward keyhole limpet hemocyanin (KLH) was detected on chicken chromosome 14 in the experimental population, which was created by crossing commercial White Leghorn and a Polish native chicken breed (green-legged partridgelike). The current QTL location is a validation of previous experiments pointing to the same genomic location for the QTL linked to a primary antibody response to KLH. An experimental population was typed with microsatellite markers distributed over the chicken chromosome 14. Titers of antibodies binding KLH were measured for all individuals by ELISA. Statistical models applied in the Grid QTL Web-based software were used to analyze the data: a half-sib model, a line-cross model, and combined analysis in a linkage disequilibrium and linkage analysis model. Candidate genes that have been proposed were genotyped with SNP located in genes exons. Statistical analyses of single SNP associations were performed pointing out 2 SNP of an axis inhibitor protein (AXIN1) gene as significantly associated with the trait of an interest.

Across-line SNP association study of innate and adaptive immune response in laying hens

The aim of the present study was to detect quantitative trait loci (QTL) for innate and adaptive immunity in laying hens. For this purpose, the associations between 1022 single nucleotide polymorphism (SNP) markers and immune traits were studied in 583 hens from nine different layer lines. Immune traits were natural antibodies for keyhole limpet haemocyanin (KLH) and lipopolysaccharide (LPS) at 20, 40 and 65 weeks, acquired antibodies to the vaccinal virus of Newcastle disease at 20 weeks, and complement activity measured on sheep and bovine red blood cells at 20, 40 and 65 weeks. We adopted a novel approach based on across-line analysis and testing of the SNP-by-line interaction. Among lines, linkage disequilibrium is conserved at shorter distances than in individual lines; therefore, SNPs significantly associated with immune traits across lines are expected to be near the functional mutations. In the analysis, the SNPs that had a significant across-line effect but did not show significant SNP-by-line interaction were identified to test whether the association was consistent in the individual lines. Ultimately, 59 significant associations between SNPs and immune traits were detected. Our results confirmed some previously identified QTL and identified new QTL potentially involved in the immune function. We found evidence for a role of IL17A (chromosome 3) in natural and acquired antibody titres and in the classical and alternative pathways of complement activation. The major histocompatibility genes on chromosome 16 showed significant association with natural and acquired antibody titres and classical complement activity. The IL12B gene on chromosome 13 was associated with natural antibody titres.

Correlated effects of selection for immunity in White Leghorn chicken lines on natural antibodies and specific antibody responses to KLH and M. butyricum

Background

The effect of selection for three general immune response traits on primary antibody responses (Ab) to Mycobacterium butyricum or keyhole limpet hemocyanin (KLH) was studied in four experimental lines of White Leghorn chicken. Birds underwent 12 generations of selection for one of three different general immune criteria; high antibody response to Newcastle disease virus 3 weeks after vaccination (ND3), high cell-mediated immune response, using the wing web response to phytohemglutinin (PHA) and high phagocytic activity, measured as carbone clearance (CC). Line ND3-L was selected on ND3, line PHA-L was selected on PHA, and line CC-L on CC, but all lines were measured for all three traits. The fourth line was a contemporary random bred control maintained throughout the selection experiment. Principal component analysis was used to distinguish clusters based on the overall set of immune measures.

Results

In the KLH immunised group, no differences were present between lines for natural antibodies binding to KLH and LPS, and, lines ND3-L and PHA-L had higher titers to LTA and anti-Gal titers measured before the immunisation protocol. The measure of ND3 was correlated positively with LPS titers measured post KLH immunisation and with the difference between LPS titers measured at day 0 and 7 post immunisation. In the M. butyricum immunised group, Line ND3-L showed significantly higher specific antibody response to M. butyricum, and this result agrees well with the hypothesis that the Th-1 pathway was expected to be selected for in this line.

Conclusion

This study has shown that the two different antigens KLH and M. butyricum gave rise to different responses in the set of selected lines, and that the response was only enhanced for the antigen associated with the same response mechanism as that for the trait (ND3, PHA or CC) for which the line was selected. Interactions between innate and acquired immunity have been observed mainly for the high antibody selected trait, indicating there was a specific interaction due to the selection criterion. Furthermore, the results confirmed the independence between the three selected traits. Finally, principal component analysis contributed to visually discriminate high and low responders to the two new antigens in the four lines.

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 Traits

Planned crosses were designed to produce an F(2) and 2 backcross populations from 2 lines of White Leghorn chickens previously selected over 10 generations for 2 different in vivo immune responses. The selection criteria applied on the 2 grandparental lines were as follows: high antibody response to Newcastle disease virus vaccine 3 wk after vaccination (ND3) and high cell-mediated immune response [response to phytohemagglutinin]. Furthermore a control line was kept by random breeding. The objective of the study was to estimate if the 2 selection criteria applied on the pure lines had changed the level of and type of immune (humoral) response to a new antigen, keyhole limpet hemocyanin (KLH), in the various second-generation progeny groups. In addition, correlations between parameters of acquired and innate immunity were tested. Primary total (IgT) and isotype-specific (IgG and IgM) antibody response to KLH 1 wk after immunization and levels of natural antibodies (NAB) binding to Salmonella enteriditis-derived lipopolysaccharide (LPS) were measured. Although no differences were present between IgM and IgG antibodies to KLH and the phytohemagglutinin skin-swelling response, significant differences were present between all the progeny groups for IgT to KLH and ND3 and NAB binding to LPS. The mean values for IgT to ND3 and KLH were significantly different between the crosses using the selected lines compared with the control line, indicating a contribution of the previous selection. In addition, a sex effect was found for IgM to KLH and NAB to LPS, for which females had a higher response than males in both cases. No interaction between progeny type and sex was found. Furthermore, significant positive correlations were found between NAB to LPS and specific antibody titers to KLH. Finally, the results of the present study demonstrated an interaction between innate and acquired immunity under this strategy of selection and crossbreeding and confirmed the effect of selection on general immune response to a new antigen in second-generation crosses.

Review of Quantitative Trait Loci Identified in the Chicken

Methods for mapping QTL are actively used in the chicken to identify chromosomal regions contributing to variation in traits related to growth, disease resistance, egg production, behavior, and metabolic parameters. However, higher-resolution mapping and better knowledge of the genetic architecture underlying QTL are needed for successful application of this information into breeding programs. Therefore, this paper summarizes and integrates original, primary QTL studies in the chicken to identify basic information on the genetic architecture of quantitative traits in chickens. The results of this review show several instances of consensus of QTL locations for similar traits from independent studies. Furthermore, the consensus of QTL location for different traits and evidence for QTL with parent-of-origin effect, transgressive alleles, epistatic QTL, and QTL x sex interaction in chicken are presented and discussed. This information can be helpful in identifying genes or mutations underlying the QTL and in the application of genomic information in marker-assisted breeding programs.

Detection of QTL for innate: non-specific antibody levels binding LPS and LTA in two independent populations of laying hens

In the current study results are presented of an experiment dealing with the Natural antibodies which are measured by level of homotopes LPS and LTA. Two independent populations were examined (F2 population descendant from a cross between chickens divergently selected for either High or Low specific Ab responses to SRBC (HL) and F2 cross descendant from lines expressed different behavior concerning feather pecking (FP)). In total 12 QTL were detected to non-specific antibody titers directed to LTA and LPS and at two ages after applying two statistical models in an F2 HL population. Similarly in an FP cross overall seven QTL were detected. Based on obtained results it might be concluded that different QTL regions are associated with immune responses to homotopes LPS and LTA in poultry.

The genetic dissection of immune response using gene-expression studies and genome mapping

Functional genomics has been applied to the genetic dissection of immune response in different ways: (1) experimental crosses between lines that differ in their (non-) specific immune response have been used to detect quantitative trait loci (QTL) underlying these differences. (2) The measurement of gene expression levels for thousands of genes using microarrays or oligonucleotide chips to identify differential expression with regard to antigen challenge: (a) before and after infection, (b) resistant versus susceptible lines, or (c) combinations of both. Interpretation of QTL results is hampered by the fact that confidence regions of the QTL are large and can contain hundreds of potential candidate genes for the QTL. At the same time, the microarray experiments tend to show large numbers of differentially expressed genes without identifying the relationships between these genes. In the recently proposed 'genetical genomics' framework, members of a segregating population are characterised for genome-wide molecular markers and for gene expression levels. This facilitates the mapping of expression-QTL (eQTL): loci in the genome that control the expression of genes. Initial applications of this approach are critically reviewed and potential applications of this approach with regard to immune response are presented.

Genetic and Phenotypic Correlations Between Feather Pecking Behavior, Stress Response, Immune Response, and Egg Quality Traits in Laying Hens

The objective of the current study was to estimate genetic and phenotypic correlations among feather pecking (FP) behavior and stress response, immune response, and egg quality parameters. These traits have been measured in an F2 cross, coming from a cross between a high and a low FP line of laying hens. Heritabilities (h2) of stress response (32 wk), primary immune response to keyhole limpet hemocyanin (KLH) (36 wk) and Mycobacterium butyricum (39 wk), and egg quality parameters (35, 44, and 50 wk of age) were calculated. The h2 was 0.05 +/- 0.05 (SE) for stress response, 0.15 +/- 0.07 for antibody response to KLH, and 0.08 +/- 0.06 for antibody response to M. butyricum. The h2 for egg quality traits were in the range of 0.12 to 0.30. Significant phenotypic correlations were found between gentle FP in adult hens and stress response, egg weight at 44 and 50 wk, and egg deformation at 50 wk. Significant additive genetic correlations were found between severe FP in adult hens and antibody response to KLH (0.79 +/- 0.35), and between ground pecking in adult hens and egg deformation at 50 wk (0.63 +/- 0.26), and between ground pecking and eggshell strength at 35, 44, and 50 wk of age (-0.86 +/- 0.29, -0.81 +/- 0.20, -0.76 +/- 0.24, respectively).