Immune response and resistance to infectious bursal disease virus of chicken lines selected for high or low antibody response to Escherichia coli

Genetic resistance to avian pathogenic Escherichia coli (APEC): current status and opportunities

Infections with avian pathogenic Escherichia coli (APEC) can be extremely detrimental to poultry health and production. Investigating host genetic variation could identify the biological mechanisms that control resistance to this pathogen and allow selection for improved resistance in experimental and commercial poultry populations. In this review, the current knowledge of how host genetics contributes to APEC resistance and future opportunities that would benefit the understanding or application of genetic resistance are discussed. Phenotypes, such as antibody responses, lesion scores, and mortality, revealed that genetic background impacts APEC resistance and interacts with other factors including the environment and challenge conditions. Experiments have used divergent selection for APEC-specific antibody levels to facilitate genetic studies, estimated heritabilities in relevant traits, detected quantitative trait loci using microsatellites, and made associations with sequence variation in the major histocompatibility complex, which collectively suggest that improving APEC resistance through selection is feasible, although genetic control is partial, complex, and highly polygenic. Additionally, functional genomics techniques have identified antimicrobial responses, toll-like receptor and cytokine signalling, and the cell cycle as central pathways in the host response to APEC challenge. Opportunities for future research are discussed, including the expansion of existing lines of research and the application of new technologies that are relevant to the study of host genetics and APEC. This review closes with prospective strategies for improvement of host genetic resistance to APEC.

Genome-wide association studies of immune, disease and production traits in indigenous chicken ecotypes

Background

The majority of chickens in sub-Saharan Africa are indigenous ecotypes, well adapted to the local environment and raised in scavenging production systems. Although they are generally resilient to disease challenge, routine vaccination and biosecurity measures are rarely applied and infectious diseases remain a major cause of mortality and reduced productivity. Management and genetic improvement programmes are hampered by lack of routine data recording. Selective breeding based on genomic technologies may provide the means to enhance sustainability. In this study, we investigated the genetic architecture of antibody response to four major infectious diseases [infectious bursal disease (IBDV), Marek’s disease (MDV), fowl typhoid (SG), fowl cholera (PM)] and resistance to Eimeria and cestode parasitism, along with two production traits [body weight and body condition score (BCS)] in two distinct indigenous Ethiopian chicken ecotypes. We conducted variance component analyses, genome-wide association studies, and pathway and selective sweep analyses.

Results

The large majority of birds was found to have antibody titres for all pathogens and were infected with both parasites, suggesting almost universal exposure. We derived significant moderate to high heritabilities for IBDV, MDV and PM antibody titres, cestodes infestation, body weight and BCS. We identified single nucleotide polymorphisms (SNPs) with genome-wide significance for each trait. Based on these associations, we identified for each trait, pathways, networks and functional gene clusters that include plausible candidate genes. Selective sweep analyses revealed a locus on chromosome 18 associated with viral antibody titres and resistance to Eimeria parasitism that is within a positive selection signal. We found no significant genetic correlations between production, immune and disease traits, implying that selection for altered antibody response and/or disease resistance will not affect production.

Conclusions

We confirmed the presence of genetic variability and identified SNPs significantly associated with immune, disease and production traits in indigenous village chickens. Results underpin the feasibility of concomitant genetic improvement for enhanced antibody response, resistance to parasitism and productivity within and across indigenous chicken ecotypes.

Electronic supplementary material

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

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 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.

Differences in genetic background influence the induction of innate and acquired immune responses in chickens depending on the virulence of the infecting infectious bursal disease virus (IBDV) strain

Previous studies and field observations have suggested that genetic background influences infectious bursal disease virus (IBDV) pathogenesis. However, the influence of the virulence of the infecting IBDV strain and the mechanisms underlying the differences in susceptibility are not known. In the present study IBDV pathogenesis was compared between specific-pathogen-free layer-type (LT) chickens, which are the most susceptible chicken for IBDV and have been used as the model for pathogenesis studies, and broiler-type (BT) chickens, which are known to be less susceptible to clinical infectious bursal disease (IBD). The innate and acquired immune responses were investigated after inoculation of an intermediate (i), virulent (v) or very virulent (vv) strain of IBDV. IBDV pathogenesis was comparable among genetic backgrounds after infection with iIBDV. After infection with vIBDV and vvIBDV, LT birds showed severe clinical disease and mortality, higher bursal lesion scores and IBDV-antigen load relative to BT birds. Circulating cytokine induction varied significantly in both timing and quantity between LT and BT birds and among virus strains (P<0.05). Evaluation of different immune cell populations by flow-cytometric analysis in the bursa of Fabricius provided circumstantial evidence of a stronger local T cell response in BT birds vs. LT birds after infection with the virulent strain. On the other hand, LT birds showed a more significant increase in circulating macrophage-derived immune mediators such as total interferon (IFN) and serum nitrite than BT birds on days 2 and 3 post-vIBDV infection (P<0.05). Stronger stimulation of innate immune reactions especially after vIBDV infection in the early phase may lead to faster and more severe lesion development accompanied by clinical disease and death in LT chickens relative to BT chickens. Interestingly, no significant differences were seen between genetic backgrounds in induction of the IBDV-specific humoral response: timing of IBDV-antibody induction and antibody levels were comparable between BT and LT birds. This study clearly demonstrates a significant influence of chickens' genetic background on disease outcome. The difference between backgrounds in IBDV susceptibility is further influenced by the virulence of the infecting virus strain.

Bilateral effects of vaccination against infectious bursal disease and Newcastle disease in specific-pathogen-free layers and commercial broiler chickens

Different infectious bursal disease virus (IBDV) live vaccines (intermediate, intermediate plus) were compared for their immunosuppressive abilities in specific-pathogen-free (SPF) layer-type chickens or commercial broilers. The Newcastle disease virus (NDV) vaccination model was applied to determine not only IBDV-induced immunosuppression but also bilateral effects between IBDV and NDV. None of the IBDV vaccines abrogated NDV vaccine-induced protection. All NDV-vaccinated SPF layers and broilers were protected against NDV challenge independent of circulating NDV antibody levels. Sustained suppression of NDV antibody development was observed in SPF layers, which had received the intermediate plus IBDV vaccine. We observed a temporary suppression of NDV antibody development in broilers vaccinated with one of the intermediate, as well as the intermediate plus, IBDV vaccines. Different genetic backgrounds, ages, and residual maternal antibodies might have influenced the pathogenesis of IBDV in the different types of chickens. Temporary suppression of NDV antibody response in broilers was only seen if the NDV vaccine was administered before and not, as it was speculated previously, at the time the peak of IBDV-induced bursa lesions was detected. For the first time, we have demonstrated that the NDV vaccine had an interfering effect with the pathogenesis of the intermediate as well as the intermediate plus IBDV vaccine. NDV vaccination enhanced the incidence of IBDV bursa lesions and IBDV antibody development. This observation indicates that this bilateral effect of an IBDV and NDV vaccination should be considered in the field and could have consequences for the performance of broiler flocks.

A comparison of immune responses to infection with virulent infectious bursal disease virus (IBDV) between specific-pathogen-free chickens infected at 12 and 28 days of age

Previous studies indicated that the pathogenesis of infectious bursal disease virus (IBDV) may vary depending on the age of the infected chicken. The reason for this difference is not clear, immunopathogenesis studies comparing different age groups are lacking. In the present study the IBDV-immunopathogenesis was compared in 12-day-old and 28-day-old specific-pathogen-free layer-type chickens. Birds were inoculated with the classic virulent strain IBDV-IM. Interestingly, the up-regulation of interferon (IFN)-gamma and interleukin-1beta expression varied between the two IBDV-infected age groups. The cytokine expression pattern was different between IBDV-infected age groups in magnitude and timing not only in bursal tissue and caecal tonsils but also in the spleen. No variations between the different age groups were observed in IBDV-induced bursa lesions, accumulation of intrabursal T cells and macrophages, and systemic stimulation of the release of nitric oxide inducing factors. Furthermore, virus-neutralizing antibody levels were comparable between age groups. Under these experimental conditions, age-related differences were found only for the cytokine patterns but interestingly this did not influence the outcome of the disease with IBDV-IM.

Differential expression of inducible nitric oxide synthase and cytokine mRNA in chicken lines divergent for cutaneous hypersensitivity response

Phytohemagglutinin (PHA)-induced delayed-type hypersensitivity is an immunocompetent trait considered an indicator of cell-mediated immune or T-cell responses. Divergent selection was performed to generate high and low lines for response to PHA-P. Extreme-responder birds of the F2 generation in each line were used to study possible differences in macrophage activity and the associated functional genes. To evaluate macrophage activity, nitric oxide (NO) was estimated both systemically in serum and in in vitro monocyte culture. Semi-quantitative RT-PCR was used to detect the differential mRNA expression patterns of iNOS and MIP-1beta in monocyte culture, whereas T(H)1 cytokines (IL-2 and IFN-gamma) were studied in peripheral blood mononuclear cells (PBMC) at different time intervals after lipopolysaccharide (LPS) induction. The high line showed strong systemic, as well as in vitro NO production, compared to the low line, upon stimulation with NDV and LPS, similar to early and high iNOS mRNA expression. Following the pattern of iNOS gene expression, an early strong expression of cytokines with powerful iNOS-inducing action, such as IFN-gamma and the chemokine MIP-1beta, was observed in the high line. In contrast, for response to PHA-P, low expression of IL-2 was observed in the high compared to the low line. In conclusion, the study revealed that divergent selection for response to PHA-P resulted in a divergent effect on T(H)1 cell activity, resulting in altered macrophage function in chickens. Selection, based on response to PHA-P, could lead to more resistant birds or birds with an enhanced immune response.

Behavioral and physiological features of chickens diversely selected for resistance to Avian Disease. 1. Selected inbred lines differ in behavioral and physical responses to social stress

To test the hypothesis that genetic variations in response to social stress modulate susceptibility to disease in poultry, aggressive behaviors induced by social stress were measured in chickens of different inbred lines selected for disease resistance (line 63) or susceptibility (lines 72 and 15I5), as well as 2 recombinant congenic strains (B and X). At 15 wk of age, roosters from each genetic line or strain were randomly assigned to pairs for intraline male-male aggression tests (n = 8 per line). Based on the results of the intraline aggression tests, the roosters were divided into 2 groups, winners and losers. At 16 wk of age, the roosters were randomly paired as winners vs. winners and losers vs. losers for interline aggression tests, i.e., line 63 vs. 72 and 15I5; line 73 vs. line 15I5; and strain X vs. strain B. Similarly, at 17 wk of age, line 63 vs. strains X and B, and line 72 vs. strains X and B were tested. The tests were conducted in a novel cage that was similar to their home cages, to provide a neutral space for both roosters being tested. Each pair was videotaped for 15 min. Male-male interaction-induced aggressive behaviors were markedly different among the genetic lines. Compared with roosters of lines 15I5 and 72, line 63 roosters generally showed fewer aggressive behaviors, including aggressive pecks and fights, as well as durations (P < 0.05). Roosters of the recombinant congenic strains X and B, each possessing a unique random 87.5% genome of line 63, exhibited low aggressive behaviors, which were similar or equal to the level of line 63 in both intraline and interline aggression tests (P = 0.05). These results may indicate that some of the gene(s) commonly carried between strains X and B as well as line 63 likely played an important role in governing their lower levels of aggression. The present chicken lines may be used as animal models for investigation of the cellular mechanisms of genetic-environmental interactions on disease resistance and stress responses.

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

Quantitative trait loci involved in the primary antibody response to keyhole lympet hemocyanin (KLH) and Mycobacterium butyricum were detected in two independent populations of laying hens. The first population was an F2 cross (H/L) of lines divergently selected for either high or low primary antibody responses to SRBC, and the second population was an F2 cross between 2 commercial layer lines displaying differences in feather pecking behavior (FP). Both populations were typed with microsatellite markers widely distributed over the genome with similar intervals between markers. Titers of antibodies binding KLH and M. butyricum were measured for all individuals by ELISA. Two genetic models were applied to detect QTL involved in the humoral immune response: a half-sib model and a line-cross model, both using the regression interval method. In the half-sib analysis, 2 QTL (on GGA14 and GGA27) were detected for the antibody response to KLH for the H/L population, and 2 QTL (on GGA14 and GGA18) were detected for the FP population. Only 1 QTL was detected for M. butyricum on GGA14 in the FP population using the half-sib analysis model. Two QTL were detected for the FP population on GGA2 and GGA3 using the line-cross analysis model. A QTL for the primary antibody response to KLH detected on GGA14 was validated in both populations under the half-sib analysis model. The present data suggest differences in the genetic regulation of antibody responses to two different T-cell dependent antigens.

Reproductive effort reduces long-term immune function in breeding tree swallows (Tachycineta bicolor)

We examined whether strategies of reproductive allocation may reduce long-term immunocompetence through the effects of manipulated effort on secondary or acquired immunity. We tested whether increased reproductive effort leads to reduced immune function and survival by manipulating brood size in tree swallows (Tachycineta bicolor) and exposing breeding females to a primary and secondary exposure of sheep red blood cells to elicit a humoral immune response. Females raising enlarged broods produced fewer secondary antibodies than did females raising control or reduced broods. Most importantly, individuals with high secondary responses were more likely to survive to breed 3 years after brood manipulations, suggesting that differences in disease susceptibility may be caused by trade-offs in reproductive allocation. We also found that individual quality, measured by clutch initiation date, mediated the effects of brood manipulations, with higher-quality birds showing a greater ability to deal with increases in effort.

Associations of six candidate genes with antibody response kinetics in hens

The chicken B-cell marker (ChB6), caspase-1, inhibitor of apoptosis protein-1 (IAP-1), interleukin-15 receptor alpha-chain (IL-15Ralpha), interleukin-2 receptor gamma-chain (IL-2Rgamma), and immunoglobulin supfamily gene (ZOV3), as physiological candidate genes for chicken immune response, were selected to investigate associations with antibody kinetics to SRBC and killed B. abortus. An F2 population was derived from mating highly inbred (>99%) males of two MHC-congenic Fayoumi lines (named M5.1 and M15.2) with G-B1 Leghorn hens. Antibody response to SRBC and B. abortus after immunization at 19 and 22 wk were measured. Secondary phase parameters of maximum titers (Ymax) and time required to achieve Ymax (Tmax) were estimated from postsecondary titers by using a nonlinear regression model. The DNA polymorphisms of six genes were identified, and associations of single nucleotide polymorphisms (SNP) in the six genes with antibody response parameters were analyzed. Significant main effects of the gene polymorphisms were mostly found in the lineage of the M5.1 grandsire and primarily on antibody response to B. abortus. There was general agreement of allelic effect within antibody parameters among genes. These results suggest that the SNP characterized in the study may serve as markers for genetic enhancement of humoral immune capacity in the chicken.

Microsatellite markers associated with quantitative trait loci controlling antibody response to Escherichia coli and Salmonella enteritidis in young broilers

A unique resource population was produced to facilitate detection of microsatellite markers associated with quantitative trait loci controlling antibody (Ab) response in broiler chickens. Three F1 males were produced by mating two lines divergently selected on Ab response to Escherichia coli vaccination. Each F1 male was mated with females from four genetic backgrounds: F1, high-Ab line (HH), low-Ab line and commercial line, producing three resource families, each with four progeny types. About 1700 chicks were immunized with E. coli and Salmonella enteritidis vaccines. Selective genotyping was conducted on the individuals with highest or lowest average Ab to E. coli and S. enteritidis within each progeny type in each sire family. Twelve markers were significantly associated with Ab to E. coli and six of them were also associated with Ab to S. enteritidis, mostly exhibiting a similar low effect (approximately 0.35 phenotypic SD) in all progeny types. Four markers exhibited a highly significant and much larger effect (approximately 1.7 SD), but only in progeny of females from the HH, suggesting that a backcross to the high parental line should be preferred over the commonly used F2 population. Results from two markers suggested a quantitative trait locus on chromosome 2 around 400 cM. The marker MCW0083, significant in two sire families, is closely linked to the bone morphogenetic protein 2 (BMP2) gene, known to be associated with the control of T-cell transformation in humans.

Genetic and phenotypic correlations between antibody responses to Escherichia coli, infectious bursa disease virus (IBDV), and Newcastle disease virus (NDV), in broiler lines selected on antibody response to Escherichia coli

The genetic control of antibody (Ab) response to Escherichia coli (EC), infectious bursa disease virus, and Newcastle disease virus and the genetic and phenotypic correlation between these Ab responses, were evaluated under farm conditions in which chicks were simultaneously exposed to these antigens. The experimental population comprised five groups: two lines divergently selected for high (HH) or low (LL) Ab response to EC vaccination; a commercial broiler dam-line (CC), from which HH and LL had been derived; and the HH x CC and LL x CC hybrid groups (HC and LC, respectively). Lines LL and HH expressed similar symmetric divergence to all three antigens. The ranking of the LL, LC, CC, HC, and HH genetic groups according to their mean Ab responses and their very high linear correlation with the LL vs. HH genomic scale clearly indicate the additive nature of the genetic divergence between these lines. Several estimates of correlation were calculated between Ab responses of each pair of antigens and between BW and Ab to each antigen. The high correlation between group means, the near-zero within-group correlation, and the low phenotypic correlation indicate the strongly positive genetic correlation between Ab responses and no correlation with BW. The results of this study suggest that overall immunocompetence of commercial broilers can be improved by selection for high Ab response of young chicks to controlled immunization with a single antigen, without counteracting further selection for high BW.

Antibody responses and morbidity following infection with infectious bronchitis virus and challenge with Escherichia coli, in lines divergently selected on antibody response

We evaluated the association between antibody (Ab) production and disease resistance. A controlled-challenge protocol was developed to mimic natural infection and to yield a higher rate of mortality following Escherichia coli (EC) challenge. Chicks were first infected with infectious bronchitis virus (IBV) by injecting a high dose of vaccine (attenuated virus) into their air sacs and then were infected with pathogenic EC introduced intratracheally. The experimental population consisted of lines divergently selected for high (HH) or low (LL) Ab response to EC vaccination, an HH x LL cross (HL), and commercial broilers (CC). When chicks were vaccinated with EC vaccine, mean Ab titer 15 d post-EC challenge was threefold higher in HH than LL lines, but both lines exhibited very low mortality (approximately 2%). When chicks were not vaccinated prior to EC challenge, high mortality (8 to 20%) occurred in the slow-growing HH, LL, and HL lines, and much higher mortality (approximately 40%) occurred among the CC broilers that were 38% heavier than the HH, LL, and HL lines. Mean level of Ab to EC, 7 d after EC challenge, was about twofold higher in HH vs. LL chicks and intermediate in HL and CC chicks. Within each line, Ab levels were higher in chicks exhibiting colibacillosis than in healthy ones, suggesting that these Ab were produced as a result of ongoing infection but were too late to fully prevent morbidity and mortality. These results indicate that rapid growth rate substantially reduces broiler viability, whereas Ab levels produced in response to acute pathogenic challenge without prior vaccination do not contribute to disease resistance. Among the relatively slow-growing lines, mortality was about twofold higher in HH than in LL lines. This finding may confirm previous reports that without prior vaccination, high Ab response to acute challenge increases consequent mortality; alternatively, the LL line may be superior in nonspecific defense mechanisms.