Kellie Watson

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

BackgroundHigh 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.ResultsWe 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.ConclusionsThis 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.

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.

Accuracy of genomic prediction using an evenly spaced, low-density single nucleotide polymorphism panel in broiler chickens

One approach for cost-effective implementation of genomic selection is to genotype training individuals with a high-density (HD) panel and selection candidates with an evenly spaced, low-density (ELD) panel. The purpose of this study was to evaluate the extent to which the ELD approach reduces the accuracy of genomic estimated breeding values (GEBV) in a broiler line, in which 1,091 breeders from 3 generations were used for training and 160 progeny of the third generation for validation. All birds were genotyped with an Illumina Infinium platform HD panel that included 20,541 segregating markers. Two subsets of HD markers, with 377 (ELD-1) or 766 (ELD-2) markers, were selected as ELD panels. The ELD-1 panel was genotyped using KBiosciences KASPar SNP genotyping chemistry, whereas the ELD-2 panel was simulated by adding markers from the HD panel to the ELD-1 panel. The training data set was used for 2 traits: BW at 35 d on both sexes and hen house production (HHP) between wk 28 and 54. Methods Bayes-A, -B, -C and genomic best linear unbiased prediction were used to estimate HD-marker effects. Two scenarios were used: (1) the 160 progeny were ELD-genotyped, and (2) the 160 progeny and their dams (117 birds) were ELD-genotyped. The missing HD genotypes in ELD-genotyped birds were imputed by a Gibbs sampler, capitalizing on linkage within families. In scenario (1), the correlation of GEBV for BW (HHP) of the 160 progeny based on observed HD versus imputed genotypes was greater than 0.94 (0.98) with the ELD-1 panel and greater than 0.97 (0.99) with the ELD-2 panel. In scenario (2), the correlation of GEBV for BW (HHP) was greater than 0.92 (0.96) with the ELD-1 panel and greater than 0.95 (0.98) with the ELD-2 panel. Hence, in a pedigreed population, genomic selection can be implemented by genotyping selection candidates with about 400 ELD markers with less than 6% loss in accuracy. This leads to substantial savings in genotyping costs, with little sacrifice in accuracy.

Identification of stable reference genes for quantitative PCR in cells derived from chicken lymphoid organs.

Quantitative polymerase chain reaction (qPCR) is a powerful technique for quantification of gene expression, especially genes involved in immune responses. Although qPCR is a very efficient and sensitive tool, variations in the enzymatic efficiency, quality of RNA and the presence of inhibitors can lead to errors. Therefore, qPCR needs to be normalised to obtain reliable results and allow comparison. The most common approach is to use reference genes as internal controls in qPCR analyses. In this study, expression of seven genes, including β-actin (ACTB), β-2-microglobulin (B2M), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), β-glucuronidase (GUSB), TATA box binding protein (TBP), α-tubulin (TUBAT) and 28S ribosomal RNA (r28S), was determined in cells isolated from chicken lymphoid tissues and stimulated with three different mitogens. The stability of the genes was measured using geNorm, NormFinder and BestKeeper software. The results from both geNorm and NormFinder were that the three most stably expressed genes in this panel were TBP, GAPDH and r28S. BestKeeper did not generate clear answers because of the highly heterogeneous sample set. Based on these data we will include TBP in future qPCR normalisation. The study shows the importance of appropriate reference gene normalisation in other tissues before qPCR analysis.

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

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

Use of high-density SNP data to identify patterns of diversity and signatures of selection in broiler chickens

Summary The development of broiler chickens over the last 70 years has been accompanied by large phenotypic changes, so that the resulting genomic signatures of selection should be detectable by current statistical techniques with sufficiently dense genetic markers. Using two approaches, this study analysed high‐density SNP data from a broiler chicken line to detect low‐diversity genomic regions characteristic of past selection. Seven regions with zero diversity were identified across the genome. Most of these were very small and did not contain many genes. In addition, fifteen regions were identified with diversity increasing asymptotically from a low level. These regions were larger and thus generally included more genes. Several candidate genes for broiler traits were found within these ‘regression regions’, including IGF1,GPD2 and MTNR1AI. The results suggest that the identification of zero‐diversity regions is too restrictive for characterizing regions under selection, but that regions showing patterns of diversity along the chromosome that are consistent with selective sweeps contain a number of genes that are functional candidates for involvement in broiler development. Many regions identified in this study overlap or are close to regions identified in layer chicken populations, possibly due to their shared precommercialization history or to shared selection pressures between broilers and layers.

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.

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

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

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.