Selection for high and low antibody responses to sheep red blood cells influences cytokine and chemokine expression in chicken peripheral blood leukocytes and splenic tissue

White Leghorn chickens from a common founder population have been divergently selected for high (HAS) or low (LAS) antibody responses to sheep red blood cells (SRBC) for 49 generations resulting in 2 diverse lines for this trait. Much has been studied in these two lines; however, the impact of these selection pressures on cytokine and chemokine expression is not fully understood. The purpose of this study is to determine if selection for antibody response to SRBC impacts cytokine and chemokine expression in peripheral blood leukocytes (PBL) and spleen from HAS and LAS chickens. Total RNA was isolated from PBL and spleen after which mRNA expression of cytokines (IL4, IL6, IL10, TGF-β4) and chemokines (CXCL8, CCL4) were determined by quantitative real-time RT-PCR (qRT-PCR). The data were analyzed using Student's t test comparing HAS and LAS (P < 0.05) and are reported as corrected 40-CT. PBL and spleen samples were analyzed separately. With respect to PBL, expression of IL6 was higher (P < 0.05) in PBL isolated from LAS chickens compared to those from the HAS line whereas there were no differences (P > 0.05) in IL4, IL10, CXCL8, CCL4, or TGF-β4. The cytokine and chemokine mRNA expression profiles were different in the spleen between the two lines. IL4 and CXCL8 expression were higher (P < 0.05) in spleen samples from HAS chickens than LAS. The expression of IL6, IL10, CCL4, or TGF-β4 in the spleens did not differ (P > 0.05) between the lines. The data indicate that selection for specific antibody responses to SRBC impacts the cytokine and chemokine expression profile in PBL and spleens but in different ways in HAS and LAS. These studies provide insight into the influence that selection pressures for antibody responses have on different immune response components, specifically cytokines and chemokines typically involved in the innate response.

Genetic Variation in Natural and Induced Antibody Responses in Layer Chickens

Selection of livestock for disease resistance is challenging due to the difficulty in obtaining reliable phenotypes. Antibodies are immunological molecules that provide direct and indirect defenses against infection and link the activities of both the innate and adaptive compartments of the immune system. As a result, antibodies have been used as a trait in selection for immune defense. The goal of this study was to identify genomic regions associated with natural and induced antibodies in chickens using low-pass sequencing. Enzyme-linked immunosorbent assays were used to quantify innate (natural) antibodies binding KLH, OVA, and PHA and induced (adaptive) antibodies binding IBD, IBV, NDV, and REO. We collected plasma from four White Leghorn (WL), two White Plymouth Rock (WPR), and two Rhode Island Red (RIR) lines. Samples numbers ranged between 198 and 785 per breed. GWAS was performed within breed on data pre-adjusted for Line-Hatch-Sex effects using GCTA. A threshold of p = 10-6 was used to select genes for downstream annotation and enrichment analysis with SNPEff and Panther. Significant enrichment was found for the defense/immunity protein, immunoglobulin receptor superfamily, and the antimicrobial response protein in RIR; and the immunoglobulin receptor superfamily, defense/immunity protein, and protein modifying enzyme in WL. However, none were present in WPR, but some of the selected SNP were annotated in immune pathways. This study provides new insights regarding the genetics of the antibody response in layer chickens.

Differences in the microbiome of the small intestine of Leghorn lines divergently selected for antibody titer to sheep erythrocytes suggest roles for commensals in host humoral response

For forty generations, two lines of White Leghorn chickens have been selected for high (HAS) or low (LAS) antibody response to a low dose injection of sheep red blood cells (SRBCs). Their gut is home to billons of microorganisms and the largest number of immune cells in the body; therefore, the objective of this experiment was to gain understanding of the ways the microbiome may influence the differential antibody response observed in these lines. We achieved this by characterizing the small intestinal microbiome of HAS and LAS chickens, determining their functional microbiome profiles, and by using machine learning to identify microbes which best differentiate HAS from LAS and associating the abundance of those microbes with host gene expression. Microbiome sequencing revealed greater diversity in LAS but statistically higher abundance of several strains, particularly those of Lactobacillus, in HAS. Enrichment of microbial metabolites implicated in immune response such as lactic acid, short chain fatty acids, amino acids, and vitamins were different between HAS and LAS. The abundance of several microbial strains corresponds to enriched host gene expression pathways related to immune response. These data provide a compelling argument that the microbiome is both likely affected by host divergent genetic selection and that it exerts influence on host antibody response by various mechanisms.

Combining supervised machine learning with statistics reveals differential gene expression patterns related to energy metabolism in the jejuna of chickens divergently selected for antibody response to sheep red blood cells

Since the 1970s, 2 lines of White Leghorn chickens, HAS and LAS, have been continuously divergently selected for 5-day postinjection antibody titer to injection with sheep red blood cells (SRBC). Antibody response is a complex genetic trait and characterizing differences in gene expression could facilitate better understanding of physiological changes due to selection and antigen exposure. At 41 d of age, randomly selected HAS and LAS chickens, which had been coraised from hatch, were either injected with SRBC (HASI and LASI) or kept as the noninjected cohort (HASN and LASN). Five days later, all were euthanized, and samples collected from the jejunum for RNA isolation and sequencing. Resulting gene expression data were analyzed combining traditional statistics with machine learning to obtain signature gene lists for functional analysis. Differences in ATP production and cellular processes were observed in the jejunum between lines and following SRBC injection. HASN vs. LASN exhibited upregulation of ATP production, immune cell motility, and inflammation. LASI exhibits upregulation of ATP production and protein synthesis vs. LASN, reflective of what was observed in HASN vs. LASN. In contrast, no corresponding upregulation of ATP production was observed in HASI vs. HASN, and most other cellular processes appear inhibited. Without exposure to SRBC, gene expression in the jejunum indicates HAS generates more ATP than LAS, suggesting HAS maintains a "primed" system; and gene expression of HASI vs. HASN further suggests this basal ATP production is sufficient for robust antibody responses. Conversely, LASI vs. LASN jejunal gene expression implies a physiological need for increased ATP production with only minimal correlating antibody production. The results of this experiment provide insight into energetic resource needs and allocations in the jejunum in response to genetic selection and antigen exposure in HAS and LAS which may help explain phenotypic differences observed in antibody response.

Genome-Wide Association Study of Growth Performance and Immune Response to Newcastle Disease Virus of Indigenous Chicken in Rwanda

A chicken genome has several regions with quantitative trait loci (QTLs). However, replication and confirmation of QTL effects are required particularly in African chicken populations. This study identified single nucleotide polymorphisms (SNPs) and putative genes responsible for body weight (BW) and antibody response (AbR) to Newcastle disease (ND) in Rwanda indigenous chicken (IC) using genome-wide association studies (GWAS). Multiple testing was corrected using chromosomal false detection rates of 5 and 10% for significant and suggestive thresholds, respectively. BioMart data mining and variant effect predictor tools were used to annotate SNPs and candidate genes, respectively. A total of four significant SNPs (rs74098018, rs13792572, rs314702374, and rs14123335) significantly (p ≤ 7.6E-5) associated with BW were identified on chromosomes (CHRs) 8, 11, and 19. In the vicinity of these SNPs, four genes such as pre-B-cell leukaemia homeobox 1 (PBX1), GPATCH1, MPHOSPH6, and MRM1 were identified. Four other significant SNPs (rs314787954, rs13623466, rs13910430, and rs737507850) all located on chromosome 1 were strongly (p ≤ 7.6E-5) associated with chicken antibody response to ND. The closest genes to these four SNPs were cell division cycle 16 (CDC16), zinc finger, BED-type containing 1 (ZBED1), myxovirus (influenza virus) resistance 1 (MX1), and growth factor receptor bound protein 2 (GRB2) related adaptor protein 2 (GRAP2). Besides, other SNPs and genes suggestively (p ≤ 1.5E-5) associated with BW and antibody response to ND were reported. This work offers a useful entry point for the discovery of causative genes accountable for essential QTLs regulating BW and antibody response to ND traits. Results provide auspicious genes and SNP-based markers that can be used in the improvement of growth performance and ND resistance in IC populations based on gene-based and/or marker-assisted breeding selection.

Major histocompatibility complex recombinant R13 antibody response against bovine red blood cells

Recombination within the chicken major histocompatibility complex (MHC) has enabled more precise identification of genes controlling immune responses. Chicken MHC genes include BF, MHC class I; BL, MHC class II; and BG, MHC class IV that are closely linked on chromosome 16. A new recombination occurred during the 10th backcross generation to develop congenic lines on the inbred Line UCD 003 (B17B17) background. Recombinant R13 (BF17-BG23) was found in a single male chick from the Line 003.R1 (BF24-BG23) backcross. An additional backcross of this male to Line UCD 003 females increased the number of R13 individuals. Two trials tested this new recombinant for antibody production against the T cell-dependent antigen, bovine red blood cells. Fifty-one progeny segregating for R13R13 (n = 10), R13B17 (n = 26), and B17B17 (n = 15) genotypes were produced by a single R13B17 male mated to 5 R13B17 dams. One milliliter of 2.5% bovine red blood cell was injected intravenously into all genotypes at 4 and 11 wk of age to stimulate primary and secondary immune responses, respectively. Blood samples were collected 7 d after injection. Serum total and mercaptoethanol-resistant antibodies against bovine red blood cell were measured by microtiter methods. The least squares ANOVA used to evaluate all antibody titers included trial and B genotype as main effects. Significant means were separated by Fisher's protected least significant difference at P < 0.05. R13R13 chickens had significantly lower primary total and mercaptoethanol-resistant antibodies than did the R13B17 and B17B17 genotypes. Secondary total and mercaptoethanol-resistant antibodies were significantly lower in R13R13 chickens than in R13B17 but not B17B17 chickens. Gene differences generated through recombination impacted the antibody response of R13 compared with B17. Secondary antibody titers were not substantially higher than the primary titers suggesting that the memory response had waned in the 7-wk interval between injections. Overall, the results suggest that the lower antibody response in R13R13 homozygotes may be caused by recombination affecting a region that contributes to higher antibody response.

Lack of statistical power as a major limitation in understanding MHC-mediated immunocompetence in wild vertebrate populations

Disentangling the sources of variation in developing an effective immune response against pathogens is of major interest to immunoecology and evolutionary biology. To date, the link between immunocompetence and genetic variation at the major histocompatibility complex (MHC) has received little attention in wild animals, despite the key role of MHC genes in activating the adaptive immune system. Although several studies point to a link between MHC and immunocompetence, negative findings have also been reported. Such disparate findings suggest that limited statistical power might be affecting studies on this topic, owing to insufficient sample sizes and/or a generally small effect of MHC on the immunocompetence of wild vertebrates. To clarify this issue, we investigated the link between MHC variation and seven immunocompetence proxies in a large sample of barn owls and estimated the effect sizes and statistical power of this and published studies on this topic. We found that MHC poorly explained variation in immunocompetence of barn owls, with small-to-moderate associations between MHC and immunocompetence in owls (effect size: .1 ≥ r ≤ .3) similar to other vertebrates studied to date. Such small-to-moderate effects were largely associated with insufficient power, which was only sufficient (>0.8) to detect moderate-to-large effect sizes (r ≥ .3). Thus, studies linking MHC variation with immunocompetence in wild populations are underpowered to detect MHC effects, which are likely to be of generally small magnitude. Larger sample sizes (>200) will be required to achieve sufficient power in future studies aiming to robustly test for a link between MHC variation and immunocompetence.

Genomic loci associated with antibody-mediated immune responses in an F2 chicken population

Immunity-related traits are heritable in chicken, therefore, it is possible to improve the inherent immunity by breeding programs. In this study using the Illumina chicken 60K single nucleotide polymorphisms (SNPs) chip, we performed a set of genome-wide association studies to determine candidate genes and loci responsible for primary and secondary antibody-mediated responses against sheep red blood cell. A F2 population descended from a commercial meat-type breed and an Iranian indigenous chicken was used for this study. Statistical analysis was based on a mixed linear model utilizing genomic relationship matrix to prevent spurious associations. Correction for multiple testing was done by applying 5% and 10% chromosomal false discovery rates (FDRs) for significant and suggestive thresholds, respectively. Nine significant and 17 suggestive associated SNPs were identified. Most of the SNPs that were suggestively associated with the primary response of total plasma immunoglobulins were also significantly associated with this trait in secondary response. Three SNPs were located within a narrow region of 23 kb on chromosome 16. Pathway analysis for the genes surrounding the associated SNPs showed that they are involve in antigen processing and presentation, primary immunodeficiency, vitamin digestion and absorption, cell adhesion molecules, phagosome, influenza A, folding, assembly and peptide loading of class I major histocompatibility complex, lipid digestion, mobilization, and transport (FDR &lt; 0.1). Interestingly, there were common regains associated with multiple immune-related traits.

Evaluation of genetic resistance to Salmonella Pullorum in three chicken lines

Resistance to diseases varies considerably among populations of the same species and can be ascribed to both genetic and environmental factors. Salmonella Pullorum (SP) is responsible for significant losses in the poultry industry, especially in developing countries. To better understand SP resistance in chicken populations with different genetic backgrounds, we orally challenged 3 chicken lines with SP-a highly selected commercial breed (Rhode Island Red, RIR), a local Chinese chicken (Beijing You, BY), and a synthetic layer line (dwarf, DW)-at 4 d of age. Two traits related to SP resistance, survival, and bacterial carriage in the spleen were evaluated after infection. Survival rates were recorded up to 40 d of age when all chickens still alive were killed to verify the presence of SP in the spleen to determine carrier state. Mortalities for RIR, BY, and DW chicks were 25.1%, 8.3%, and 22.7%, respectively, and the corresponding carrier-states in the spleens were 17.9%, 0.6%, and 15.8%. Survival and carrier-state heritabilities were estimated using an animal threshold model. Survival heritability was 0.197, 0.091, and 0.167 in RIR, BY, and DW populations, respectively, and the heritabilities of carrier state for DW and RIR were 0.32 and 0.16, respectively. This is the first time that the heritability of the SP carrier state has been evaluated in chickens. Our study provides experimental evidence that chickens with various genetic background exhibited significantly different SP-resistant activities and heritabilities. These results may be useful for selecting lines with better disease resistance.

Gut Microbiota Co-microevolution with Selection for Host Humoral Immunity

To explore coevolution between the gut microbiota and the humoral immune system of the host, we used chickens as the model organism. The host populations were two lines (HAS and LAS) developed from a common founder that had undergone 40 generations of divergent selection for antibody titers to sheep red blood cells (SRBC) and two relaxed sublines (HAR and LAR). Analysis revealed that microevolution of host humoral immunity contributed to the composition of gut microbiota at the taxa level. Relaxing selection enriched some microorganisms whose functions were opposite to host immunity. Particularly, Ruminococcaceae and Oscillospira enriched in high antibody relaxed (HAR) and contributed to reduction in antibody response, while Lactobacillus increased in low antibody relaxed (LAR) and elevated the antibody response. Microbial functional analysis showed that alterations were involved in pathways relating to the immune system and infectious diseases. Our findings demonstrated co-microevolution relationships of host-microbiota and that gut microorganisms influenced host immunity.

Genome-wide association study of cell-mediated immune response in chicken

Cell-mediated immunity (CMI) causes the intracellular destruction of the antigen or elimination of the host cell to make animals resistant against exogenous antigens and cancers. In this study, a genome-wide association study (GWAS) was carried out to identify genomic regions associated with CMI in chicken using chicken 60k high-density single nucleotide polymorphism (SNP) array. Genomic relationships were taken into account to adjust for population structure. In order to account for multiple testing, chromosome-wise false discovery rate was controlled at 5% and 10% levels. Moreover, a comparison of the power of fixed and mixed linear models based on genomic inflation factor was carried out. Mixed linear model (MLM) had better inflation rate, and therefore the results from MLM were used for subsequent analysis. Three significantly associated SNPs (FDR < 0.05) on chromosome 24 and linkage group E22C19W28_E50C23, and three suggestively associated SNPs (FDR < 0.1) on chromosome 1, 5 and 16 were identified. Pathway analysis showed that two biological pathways, which are related to immune response, were strongly associated with the candidate genes surrounding identified SNPs, and their influences were mostly on antigen processing and presentation, and cellular structure.

Genome-wide standing variation facilitates long-term response to bidirectional selection for antibody response in chickens

Background

Long-term selection experiments provide a powerful approach to gain empirical insights into adaptation, allowing researchers to uncover the targets of selection and infer their contributions to the mode and tempo of adaptation. Here we implement a pooled genome re-sequencing approach to investigate the consequences of 39 generations of bidirectional selection in White Leghorn chickens on a humoral immune trait: antibody response to sheep red blood cells.

Results

We observed wide genome involvement in response to this selection regime. Many genomic regions were highly differentiated resulting from this experimental selection regime, an involvement of up to 20% of the chicken genome (208.8 Mb). While genetic drift has certainly contributed to this, we implement gene ontology, association analysis and population simulations to increase our confidence in candidate selective sweeps. Three strong candidate genes, MHC, SEMA5A and TGFBR2, are also presented.

Conclusions

The extensive genomic changes highlight the polygenic genetic architecture of antibody response in these chicken populations, which are derived from a common founder population, demonstrating the extent of standing immunogenetic variation available at the onset of selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3414-7) contains supplementary material, which is available to authorized users.

Genetic Determinism of Fearfulness, General Activity and Feeding Behavior in Chickens and Its Relationship with Digestive Efficiency

The genetic relationships between behavior and digestive efficiency were studied in 860 chickens from a cross between two lines divergently selected on digestive efficiency. At 2 weeks of age each chick was video-recorded in the home pen to characterize general activity and feeding behavior. Tonic immobility and open-field tests were also carried out individually to evaluate emotional reactivity (i.e. the propensity to express fear responses). Digestive efficiency was measured at 3 weeks. Genetic parameters of behavior traits were estimated. Birds were genotyped on 3379 SNP markers to detect QTLs. Heritabilities of behavioral traits were low, apart from tonic immobility (0.17-0.18) and maximum meal length (0.14). The genetic correlations indicated that the most efficient birds fed more frequently and were less fearful. We detected 14 QTL (9 for feeding behavior, 3 for tonic immobility, 2 for frequency of lying). Nine of them co-localized with QTL for efficiency, anatomy of the digestive tract, feed intake or microbiota composition. Four genes involved in fear reactions were identified in the QTL for tonic immobility on GGA1.

Screening for genes involved in antibody response to sheep red blood cells in the chicken, Gallus gallus

Antibody response, an important trait in both agriculture and biomedicine, plays a part in protecting animals from infection. Dissecting molecular basis of antibody response may improve artificial selection for natural disease resistance in livestock and poultry. A number of genetic markers associated with antibody response have been identified in the chicken and mouse by linkage-based association studies, which only define genomic regions by genetic markers but do not pinpoint genes for antibody response. In contrast, global expression profiling has been applied to define the molecular bases of a variety of biological traits through identification of differentially expressed genes (DEGs). Here, we employed Affimetrix GeneChip Chicken Genome Arrays to identify differentially expressed genes for antibody response to sheep red blood cells (SRBC) using chickens challenged with and without SRBC or chickens with high and low anti-SRBC titers. The DEGs include those with known (i.e., MHC class I and IgH genes) or unknown function in antibody response. Classification test of these genes suggested that the response of the chicken to intravenous injection of SRBC involved multiple biological processes, including response to stress or other different stimuli, sugar, carbohydrate or protein binding, and cell or soluble fraction, in addition to antibody response. This preliminary study thus provides an insight into molecular basis of antibody response to SRBC in the chicken.

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.

A genome-wide detection of copy number variation using SNP genotyping arrays in Beijing-You chickens

Copy number variation (CNV) has been recently examined in many species and is recognized as being a source of genetic variability, especially for disease-related phenotypes. In this study, the PennCNV software, a genome-wide CNV detection system based on the 60 K SNP BeadChip was used on a total sample size of 1,310 Beijing-You chickens (a Chinese local breed). After quality control, 137 high confidence CNVRs covering 27.31 Mb of the chicken genome and corresponding to 2.61 % of the whole chicken genome. Within these regions, 131 known genes or coding sequences were involved. Q-PCR was applied to verify some of the genes related to disease development. Results showed that copy number of genes such as, phosphatidylinositol-5-phosphate 4-kinase II alpha, PHD finger protein 14, RHACD8 (a CD8α- like messenger RNA), MHC B-G, zinc finger protein, sarcosine dehydrogenase and ficolin 2 varied between individual chickens, which also supports the reliability of chip-detection of the CNVs. As one source of genomic variation, CNVs may provide new insight into the relationship between the genome and phenotypic characteristics.

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.

Genome-wide association study for Marek's disease mortality in layer chickens

A genome-wide association study (GWAS) using Bayesian variable selection was performed to determine genomic regions associated with mortality due to Marek's disease virus (MDV) infection in layers. Mortality (%) under experimental disease challenge (500 plaque-forming units of a very virulent plus MDV strain) was recorded for progeny groups (average 15.5 birds; range 3 to 30) of 253 genotyped sires from four generations of a brown-egg layer line. An additional generation of 43 sires with progeny data was used to validate results. Sires were genotyped with a 42K Illumina single-nucleotide polymorphism (SNP) chip. Methods BayesB (pi = 0.995) and BayesCpi, with or without weighting residuals by the size of progeny groups were applied. The proportion of genetic variance contributed by SNPs within each 1-megabase (Mb) genomic region was quantified. Average mortality was 33% but differed significantly between generations. Genetic markers explained about 11% of phenotypic variation in mortality. Correlations between genomic estimated breeding values and percentage of progeny mortality for the validation generation (sons of individuals in training) were 0.12, 0.17, 0.02, and 0.16 for BayesB, weighted BayesB, BayesCpi, and weighted BayesCpi, respectively, when using the whole genome, and 0.03, 0.20, -0.06, and 0.14, when using only SNP from the 10, 1-Mb regions, explaining the largest proportion of genetic variance according to each method. Results suggest that regions on chromosomes 2, 3, 4, 9, 15, 18, and 21 are associated with Marek's disease resistance and can be used for selection and that accounting for the size of progeny groups has a large impact on correct localization of such genomic regions.

Genome-wide association study of antibody response to Newcastle disease virus in chicken

Background

Since the first outbreak in Indonesia in 1926, Newcastle disease has become one of the most common and contagious bird diseases throughout the world. To date, enhancing host antibody response by vaccination remains the most efficient strategy to control outbreaks of Newcastle disease. Antibody response plays an important role in host resistance to Newcastle disease, and selection for antibody response can effectively improve disease resistance in chickens. However, the molecular basis of the variation in antibody response to Newcastle disease virus (NDV) is not clear. The aim of this study was to detect genes modulating antibody response to NDV by a genome-wide association study (GWAS) in chickens.

Results

To identify genes or chromosomal regions associated with antibody response to NDV after immunization, a GWAS was performed using 39,833 SNP markers in a chicken F₂ resource population derived from a cross between two broiler lines that differed in their resistance. Two SNP effects reached 5% Bonferroni genome-wide significance (P<1.26×10^-6). These two SNPs, rs15354805 and rs15355555, were both on chicken (Gallus gallus) chromosome 1 and spanned approximately 600 Kb, from 100.4 Mb to 101.0 Mb. Rs15354805 is in intron 7 of the chicken Roundabout, axon guidance receptor, homolog 2 (ROBO2) gene, and rs15355555 is located about 243 Kb upstream of ROBO2. Rs15354805 explained 5% of the phenotypic variation in antibody response to NDV, post immunization, in chickens. Rs15355555 had a similar effect as rs15354805 because of its linkage disequilibrium with rs15354805 (r²=0.98).

Conclusion

The region at about 100 Mb from the proximal end of chicken chromosome 1, including the ROBO1 and ROBO2 genes, has a strong effect on the antibody response to the NDV in chickens. This study paves the way for further research on the host immune response to NDV.

Phenotypic responses of chickens to long-term selection for high or low antibody titers to sheep red blood cells

A long-term bidirectional selection experiment was conducted to study antibody response to SRBC. Lines, high antibody selection (HAS) and low antibody selection (LAS), originating from the same White Leghorn base population had undergone 37 generations of selection for either high or low antibody response 5 d after a single intravenous injection of 0.1 mL of a 0.25% suspension of SRBC antigen. Subpopulations, where selection was relaxed, were maintained as contemporaries with the selected lines from generations 16 to 24 [high antibody relaxed (HAR) and low antibody relaxed (LAR)] and 24 to 37. Body weights were obtained at 4, 24, and 38 wk of age and at the onset of lay (BW at first egg). Also measured were age in days to first egg, percentages of hen-day ovulations and normal egg production, and percentages of normal and defective eggs from total ovulation (PNE and PDE). Selection lead to a large divergence in antibody titers between the selected lines, with a plateau reached in line LAS. Line HAS and HAR females displayed higher antibody titers, lower BW4, and matured at older ages than those from LAS and LAR (P < 0.05). Correlations between BW at 4 wk and antibody titers were different between the selected lines, being positive in line LAS and negative in line HAS. Quadratic regression models fit well with antibody titers, BW4, and PNE, with limiting values for these traits calculated based on regression curves. For line HAS, plots showed that an increased tendency of antibody titers was followed by decreased BW4 and increased PNE. For line LAS, however, antibody titers and BW4 decreased in parallel while PNE increased. It appears that at the phenotypic level there was a resource balance between immune response, growth, and reproductive traits, which during long-term selection, individuals altered their dynamic of resource allocations to satisfy certain needs.