Quantitative trait loci and transcriptome signatures associated with avian heritable resistance to Campylobacter.

Campylobacter is the leading cause of bacterial foodborne gastroenteritis worldwide. Handling or consumption of contaminated poultry meat is a key risk factor for human campylobacteriosis. One potential control strategy is to select poultry with increased resistance to Campylobacter. We associated high-density genome-wide genotypes (600K single nucleotide polymorphisms) of 3000 commercial broilers with Campylobacter load in their caeca. Trait heritability was modest but significant (h2 = 0.11 ± 0.03). Results confirmed quantitative trait loci (QTL) on chromosomes 14 and 16 previously identified in inbred chicken lines, and detected two additional QTLs on chromosomes 19 and 26. RNA-Seq analysis of broilers at the extremes of colonisation phenotype identified differentially transcribed genes within the QTL on chromosome 16 and proximal to the major histocompatibility complex (MHC) locus. We identified strong cis-QTLs located within MHC suggesting the presence of cis-acting variation in MHC class I and II and BG genes. Pathway and network analyses implicated cooperative functional pathways and networks in colonisation, including those related to antigen presentation, innate and adaptive immune responses, calcium, and renin-angiotensin signalling. While co-selection for enhanced resistance and other breeding goals is feasible, the frequency of resistance-associated alleles was high in the population studied and non-genetic factors significantly influenced Campylobacter colonisation.

Assembly of hundreds of novel bacterial genomes from the chicken caecum.

BACKGROUND: Chickens are a highly important source of protein for a large proportion of the human population. The caecal microbiota plays a crucial role in chicken nutrition through the production of short-chain fatty acids, nitrogen recycling, and amino acid production. In this study, we sequence DNA from caecal content samples taken from 24 chickens belonging to either a fast or a slower growing breed consuming either a vegetable-only diet or a diet containing fish meal. RESULTS: We utilise 1.6 T of Illumina data to construct 469 draft metagenome-assembled bacterial genomes, including 460 novel strains, 283 novel species, and 42 novel genera. We compare our genomes to data from 9 European Union countries and show that these genomes are abundant within European chicken flocks. We also compare the abundance of our genomes, and the carbohydrate active enzymes they produce, between our chicken groups and demonstrate that there are both breed- and diet-specific microbiomes, as well as an overlapping core microbiome. CONCLUSIONS: This data will form the basis for future studies examining the composition and function of the chicken caecal microbiota.

Highly multiplexed quantitative PCR-based platform for evaluation of chicken immune responses.

To address the need for sensitive high-throughput assays to analyse avian innate and adaptive immune responses, we developed and validated a highly multiplexed qPCR 96.96 Fluidigm Dynamic Array to analyse the transcription of chicken immune-related genes. This microfluidic system permits the simultaneous analysis of expression of 96 transcripts in 96 samples in 6 nanolitre reactions and the 9,216 reactions are ready for interpretation immediately. A panel of 89 genes was selected from an RNA-seq analysis of the transcriptional response of chicken macrophages, dendritic cells and heterophils to agonists of innate immunity and from published transcriptome data. Assays were confirmed to be highly specific by amplicon sequencing and melting curve analysis and the reverse transcription and preamplification steps were optimised. The array was applied to RNA of various tissues from a commercial line of broiler chickens housed at two different levels of biosecurity. Gut-associated lymphoid tissues, bursa, spleen and peripheral blood leukocytes were isolated and transcript levels for immune-related genes were defined. The results identified blood cells as a potentially reliable indicator of immune responses among all the tissues tested with the highest number of genes significantly differentially transcribed between birds housed under varying biosecurity levels. Conventional qPCR analysis of three differentially transcribed genes confirmed the results from the multiplex qPCR array. A highly multiplexed qPCR-based platform for evaluation of chicken immune responses has been optimised and validated using samples from commercial chickens. Apart from applications in selective breeding programmes, the array could be used to analyse the complex interplay between the avian immune system and pathogens by including pathogen-specific probes, to screen vaccine responses, and as a predictive tool for immune robustness.

Bi-Functional Chicken Immunoglobulin-Like Receptors With a Single Extracellular Domain (ChIR-AB1): Potential Framework Genes Among a Relatively Stable Number of Genes Per Haplotype.

The leukocyte receptor complex (LRC) in humans encodes many receptors with immunoglobulin-like (Ig-like) extracellular domains, including the killer Ig-like receptors (KIRs) expressed on natural killer (NK) cells among others, the leukocyte Ig-like receptors (LILRs) expressed on myeloid and B cells, and an Fc receptor (FcR), all of which have important roles in the immune response. These highly-related genes encode activating receptors with positively-charged residues in the transmembrane region, inhibitory receptors with immuno-tyrosine based motifs (ITIMs) in the cytoplasmic tail, and bi-functional receptors with both. The related chicken Ig-like receptors (ChIRs) are almost all found together on a microchromosome, with over 100 activating (A), inhibitory (B), and bi-functional (AB) genes, bearing either one or two extracellular Ig-like domains, interspersed over 500-1,000 kB in the genome of an individual chicken. Sequencing studies have suggested rapid divergence and little overlap between ChIR haplotypes, so we wished to begin to understand their genetics. We chose to use a hybridization technique, reference strand-mediated conformational analysis (RSCA), to examine the ChIR-AB1 family, with a moderate number of genes dispersed across the microchromosome. Using fluorescently-labeled references (FLR), we found that RSCA and sequencing of ChIR-AB1 extracellular exon gave two groups of peaks with mobility correlated with sequence relationship to the FLR. We used this system to examine widely-used and well-characterized experimental chicken lines, finding only one or a few simple ChIR haplotypes for each line, with similar numbers of peaks overall. We found much more complicated patterns from a broiler line from a commercial breeder and a flock of red junglefowl, but trios of parents and offspring from another commercial chicken line show that the complicated patterns are due to heterozygosity, indicating a relatively stable number of peaks within haplotypes of these birds. Some ChIR-AB1 peaks were found in all individuals from the commercial lines, and some of these were shared with red junglefowl and the experimental lines derived originally from egg-laying chickens. Overall, this analysis suggests that there are some simple features underlying the apparent complexity of the ChIR locus.

Development of the duodenal, ileal, jejunal and caecal microbiota in chickens.

BACKGROUND: The chicken intestinal microbiota plays a large role in chicken health and productivity and a greater understanding of its development may lead to interventions to improve chicken nutrition, disease resistance and welfare. RESULTS: In this study we examine the duodenal, jejunal, ileal and caecal microbiota of chickens from day of hatch to 5 weeks of age (day 1, 3, 7, 14 and week 5). DNA was extracted from intestinal content samples and the V4 region of the 16S rRNA gene was amplified and sequenced. We identified significant differences in microbial community composition, diversity and richness between samples taken from different locations within the chicken intestinal tract. We also characterised the development of the microbiota at each intestinal site over time. CONCLUSIONS: Our study builds upon existing literature to further characterise the development of the chicken intestinal microbiota.

Phenotypic and genetic variation in the response of chickens to Eimeria tenella induced coccidiosis.

BACKGROUND: Coccidiosis is a major contributor to losses in poultry production. With emerging constraints on the use of in-feed prophylactic anticoccidial drugs and the relatively high costs of effective vaccines, there are commercial incentives to breed chickens with greater resistance to this important production disease. To identify phenotypic biomarkers that are associated with the production impacts of coccidiosis, and to assess their covariance and heritability, 942 Cobb500 commercial broilers were subjected to a defined challenge with Eimeria tenella (Houghton). Three traits were measured: weight gain (WG) during the period of infection, caecal lesion score (CLS) post mortem, and the level of a serum biomarker of intestinal inflammation, i.e. circulating interleukin 10 (IL-10), measured at the height of the infection. RESULTS: Phenotypic analysis of the challenged chicken cohort revealed a significant positive correlation between CLS and IL-10, with significant negative correlations of both these traits with WG. Eigenanalysis of phenotypic covariances between measured traits revealed three distinct eigenvectors. Trait weightings of the first eigenvector, (EV1, eigenvalue = 59%), were biologically interpreted as representing a response of birds that were susceptible to infection, with low WG, high CLS and high IL-10. Similarly, the second eigenvector represented infection resilience/resistance (EV2, 22%; high WG, low CLS and high IL-10), and the third eigenvector tolerance (EV3, 19%; high WG, high CLS and low IL-10), respectively. Genome-wide association studies (GWAS) identified two SNPs that were associated with WG at the suggestive level. CONCLUSIONS: Eigenanalysis separated the phenotypic impact of a defined challenge with E. tenella on WG, caecal inflammation/pathology, and production of IL-10 into three major eigenvectors, indicating that the susceptibility-resistance axis is not a single continuous quantitative trait. The SNPs identified by the GWAS for body weight were located in close proximity to two genes that are involved in innate immunity (FAM96B and RRAD).

Gene expression of AvBD6-10 in broiler chickens is independent of AvBD6, 9, and 10 peptide potency.

The Avian ß-defensin (AvBD) gene cluster contains fourteen genes; within this, two groups (AvBD6/7 and AvBD8 -10) encode charged peptides of >+5 (AvBD6/7), indicative of potent microbial killing activities, and ≤+4 (AvBD8-10), suggestive of reduced antimicrobial activities. Chicken broiler gut tissues are constantly exposed to microbes in the form of commensal bacteria. This study examined whether tissue expression patterns of AvBD6-10 reflected microbial exposure and the encoded peptides a functional antimicrobial hierarchy. Gut AvBD6-10 gene expression was observed in hatch to day 21 birds, although the AvBD8-10 profiles were eclipsed by those detected in the liver and kidney tissues. In vitro challenges of chicken CHCC-OU2 cells using the gut commensal Lactobacillus johnsonii (104 CFU) did not significantly affect AvBD8-10 gene expression patterns, although upregulation (P < 0.05) of IL-Iß gene expression was observed. Similarly, in response to Bacteriodes doreii, IL-Iß and IL-6 gene upregulation were detected (P < 0.05), but AvBD10 gene expression remained unaffected. These data suggested that AvBD8-10 gene expression was not induced by commensal gut bacteria. Bacterial time-kill assays employing recombinant (r)AvBD6, 9 and 10 peptides (0.5µM - 12µM), indicated an antimicrobial hierarchy, linked to charge, of AvBD6 > AvBD9 > AvBD10 against Escherichia coli, but AvBD10 > AvBD9 > AvBD6 using Enterococcus faecalis. rAvBD10, selected due to its reduced cationic charge was, using CHCC-OU2 cells, investigated for cell proliferation and wound healing properties, but none were observed. These data suggest that in healthy broiler chicken tissues AvBD6/7 and AvBD8-10 gene expression profiles are independent of the in vitro antimicrobial hierarchies of the encoded AvBD6, 9 and 10 peptides.

Colonization of a commercial broiler line by Campylobacter is under limited genetic control and does not significantly impair performance or intestinal health.

Campylobacter is the leading bacterial cause of foodborne diarrheal illness in humans and source attribution studies unequivocally identify handling or consumption of poultry meat as a key risk factor. Campylobacter colonizes the avian intestines in high numbers and rapidly spreads within flocks. A need therefore exists to devise strategies to reduce Campylobacter populations in poultry flocks. There has been a great deal of research aiming to understand the epidemiology and transmission characteristics of Campylobacter in poultry as a means to reduce carriage rates in poultry and reduce infection in humans. One potential strategy for control is the genetic selection of poultry for increased resistance to colonization by Campylobacter. The potential for genetic control of colonization has been demonstrated in inbred populations following experimental challenge with Campylobacter where quantitative trait loci associated with resistance have been identified. Currently in the literature there is no information of the genetic basis of Campylobacter colonization in commercial broiler lines and it is unknown whether these QTL are found in commercial broiler lines. The aim of this study was to estimate genetic parameters associated with Campylobacter load and genetic correlations with gut health and production traits following natural exposure of broiler chickens to Campylobacter.The results from the analysis show a low but significant heritability estimate (0.095 ± 0.037) for Campylobacter load which indicates a limited genetic basis and that non-genetic factors have a greater influence on the level of Campylobacter found in the broiler chicken.Furthermore, through examination of macroscopic intestinal health and absorptive capacity, our study indicated that Campylobacter has no detrimental effects on intestinal health and bird growth following natural exposure in the broiler line under study. These data indicate that whilst there is a genetic component to Campylobacter colonization worthy of further investigation, there is a large proportion of phenotypic variance under the influence of non-genetic effects. As such the control of Campylobacter will require understanding and manipulation of non-genetic host and environmental factors.

AvBD1 nucleotide polymorphisms, peptide antimicrobial activities and microbial colonisation of the broiler chicken gut.

BACKGROUND: The importance of poultry as a global source of protein underpins the chicken genome and associated SNP data as key tools in selecting and breeding healthy robust birds with improved disease resistance. SNPs affecting host peptides involved in the innate defences tend to be rare, but three non-synonymous SNPs in the avian ß-defensin (AvBD1) gene encoding the variant peptides NYH, SSY and NYY were identified that segregated specifically to three lines of commercial broiler chickens Line X (LX), Line Y(LY) and Line Z. The impacts of such amino acid changes on peptide antimicrobial properties were analysed in vitro and described in relation to the caecal microbiota and gut health of LX and LY birds. RESULTS: Time-kill and radial immune diffusion assays indicated all three peptides to have antimicrobial properties against gram negative and positive bacteria with a hierarchy of NYH > SSY > NYY. Calcein leakage assays supported AvBD1 NYH as the most potent membrane permeabilising agent although no significant differences in secondary structure were identified to explain this. However, distinct claw regions, identified by 3D modelling and proposed to play a key role in microbial membrane attachment, and permeation, were more distinct in the NYH model. In vivo AvBD1 synthesis was detected in the bird gut epithelia. Analyses of the caecal gut microbiota of young day 4 birds suggested trends in Lactobacilli sp. colonisation at days 4 (9% LX vs × 30% LY) and 28 (20% LX vs 12% LY) respectively, but these were not statistically significant (P > 0.05). CONCLUSION: Amino acid changes altering the killing capacity of the AvBD1 peptide were associated with two different bird lines, but such changes did not impact significantly on caecal gut microbiota.

Effects of rearing environment on the gut antimicrobial responses of two broiler chicken lines.

To reduce the risk of enteric disease in poultry, knowledge of how bird gut innate defences mature with age while also responding to different rearing environments is necessary. In this study the gut innate responses of two phylogenetically distinct lines of poultry raised from hatch to 35days, in conditions mimicing high hygiene (HH) and low hygiene (LH) rearing environments, were compared. Analyses focussed on the proximal gut antimicrobial activities and the duodenal and caecal AvBD1, 4 and 10 defensin profiles. Variability in microbial killing was observed between individual birds in each of the two lines at all ages, but samples from day 0 birds (hatch) of both lines exhibited marked killing properties, Line X: 19±11% (SEM) and Line Y: 8.5±12% (SEM). By day 7 a relaxation in killing was observed with bacterial survival increased from 3 (Line Y (LY)) to 11 (Line X (LX)) fold in birds reared in the HH environment. A less marked response was observed in the LH environment and delayed until day 14. At day 35 the gut antimicrobial properties of the two lines were comparable. The AvBD 1, 4 and 10 data relating to the duodenal and caecal tissues of day 0, 7 and 35 birds LX and LY birds revealed gene expression trends specific to each line and to the different rearing environments although the data were confounded by inter-individual variability. In summary elevated AvBD1 duodenal expression was detected in day 0 and day 7 LX, but not LY birds, maintained in LH environments; Line X and Y duodenal AvBD4 profiles were detected in day 7 birds reared in both environments although duodenal AvBD10 expression was less sensitive to bird age and rearing background. Caecal AvBD1 expression was particularly evident in newly hatched birds. These data suggest that proximal gut antimicrobial activity is related to the bird rearing environments although the roles of the AvBDs in such activities require further investigation.

The genetic basis of pectoralis major myopathies in modern broiler chicken lines.

This is the first report providing estimates of the genetic basis of breast muscle myopathies (BMM) and their relationship with growth and yield in broiler chickens. In addition, this paper addresses the hypothesis that genetic selection for increase breast yield has contributed to the onset of BMM. Data were analyzed from ongoing recording of BMM within the Aviagen breeding program. This study focused on three BMM: deep pectoral myopathy (DPM; binary trait), white striping (WS; 4 categories) and wooden breast (WB; 3 categories). Data from two purebred commercial broiler lines (A and B) were utilized providing greater than 40,000 meat quality records per line. The difference in selection history between these two lines has resulted in contrasting breast yield (BY): 29% for Line A and 21% for Line B. Data were analyzed to estimate genetic parameters using a multivariate animal model including six traits: body weight (BW), processing body weight (PW), BY, DPM, WB, and WS, in addition to the appropriate fixed effects and permanent environmental effect of the dam. Results indicate similar patterns of heritability and genetic correlations for the two lines. Heritabilities (h2) of BW, PW and BY ranged from 0.271-0.418; for DPM and WB h2<0.1; and for WS h2≤0.338. Genetic correlations between the BMM and BW, PW, or BY were ≤0.132 in Line A and ≤0.248 in Line B. This paper demonstrates the polygenic nature of these traits and the low genetic relationships with BW, PW, and BY, which facilitates genetic improvement across all traits in a balanced breeding program. It also highlights the importance of understanding the environmental and/or management factors that contribute greater than 65% of the variance in the incidence of white striping of breast muscle and more than 90% of the variance of the incidence of wooden breast and deep pectoral myopathy in broiler chickens.

Campylobacter colonization and proliferation in the broiler chicken upon natural field challenge is not affected by the bird growth rate or breed.

The zoonotic association between Campylobacter bacteria in poultry and humans has been characterized by decades of research which has attempted to elucidate the epidemiology of this complex relationship and to reduce carriage within poultry. While much work has focused on the mechanisms facilitating its success in contaminating chicken flocks (and other animal hosts), it remains difficult to consistently exclude Campylobacter under field conditions. Within the United Kingdom poultry industry, various bird genotypes with widely varying growth rates are available to meet market needs and consumer preferences. However, little is known about whether any differences in Campylobacter carriage exist across this modern broiler range. The aim of this study was to establish if a relationship exists between growth rate or breed and cecal Campylobacter concentration after natural commercial flock Campylobacter challenge. In one investigation, four pure line genotypes of various growth rates were grown together, while in the second, eight different commercial broiler genotypes were grown individually. In both studies, the Campylobacter concentration was measured in the ceca at 42 days of age, revealing no significant difference in cecal load between birds of different genotypes both in mixed- and single-genotype pens. This is important from a public health perspective and suggests that other underlying reasons beyond genotype are likely to control and affect Campylobacter colonization within chickens. Further studies to gain a better understanding of colonization dynamics and subsequent proliferation are needed, as are novel approaches to reduce the burden in poultry.

Development of a high density 600K SNP genotyping array for chicken.

BACKGROUND: High density (HD) SNP genotyping arrays are an important tool for genetic analyses of animals and plants. Although the chicken is one of the most important farm animals, no HD array is yet available for high resolution genetic analysis of this species. RESULTS: We report here the development of a 600 K Affymetrix® Axiom® HD genotyping array designed using SNPs segregating in a wide variety of chicken populations. In order to generate a large catalogue of segregating SNPs, we re-sequenced 243 chickens from 24 chicken lines derived from diverse sources (experimental, commercial broiler and layer lines) by pooling 10-15 samples within each line. About 139 million (M) putative SNPs were detected by mapping sequence reads to the new reference genome (Gallus_gallus_4.0) of which ~78 M appeared to be segregating in different lines. Using criteria such as high SNP-quality score, acceptable design scores predicting high conversion performance in the final array and uniformity of distribution across the genome, we selected ~1.8 M SNPs for validation through genotyping on an independent set of samples (n = 282). About 64% of the SNPs were polymorphic with high call rates (>98%), good cluster separation and stable Mendelian inheritance. Polymorphic SNPs were further analysed for their population characteristics and genomic effects. SNPs with extreme breach of Hardy-Weinberg equilibrium (P < 0.00001) were excluded from the panel. The final array, designed on the basis of these analyses, consists of 580,954 SNPs and includes 21,534 coding variants. SNPs were selected to achieve an essentially uniform distribution based on genetic map distance for both broiler and layer lines. Due to a lower extent of LD in broilers compared to layers, as reported in previous studies, the ratio of broiler and layer SNPs in the array was kept as 3:2. The final panel was shown to genotype a wide range of samples including broilers and layers with over 100 K to 450 K informative SNPs per line. A principal component analysis was used to demonstrate the ability of the array to detect the expected population structure which is an important pre-investigation step for many genome-wide analyses. CONCLUSIONS: This Affymetrix® Axiom® array is the first SNP genotyping array for chicken that has been made commercially available to the public as a product. This array is expected to find widespread usage both in research and commercial application such as in genomic selection, genome-wide association studies, selection signature analyses, fine mapping of QTLs and detection of copy number variants.