Identification of parental line specific effects of MLF2 on resistance to coccidiosis in chickens

Background

Avian coccidiosis impairs the growth and feed efficiency of infected chickens [1]. There were evidences that resistance to avian coccidiosis is associated with inheritance and coccidiosis resistant line could be established by selection [2]. The resistance to avian coccidiosis QTL has been identified near two microsatellite markers on chromosome 1 [3,4]. Associations between parameters of resistance to coccidiosis and single nucleotide polymorphisms (SNPs) in 3 candidate genes located around QTL on chromosome 1 (Zyxin, TCR-b, MLF2) were determined [5,6]. These studies showed the SNPs in MLF2 are the most probable locus associated with coccidiosis resistance in chickens. The current analysis was conducted to identify the parental line specific haplotype of MLF2. Although MLF2 explained variation of body weight affected by coccidiosis, the association of causative mutation was not obviously identified [6]. Thus, this study was suggested to identify segregation of MLF2 allele associated with coccidiosis resistance for practical application and future study for identifying causative mutation. Since oocyst shedding significantly deviated from normal distribution, transformation method may affect association test. In previous studies, oocyst shedding was transformed using log transformation, but normality was fully satisfied (p ~ 0.01) by log transformation compared with Box-Cox transformation (p > 0.05). Box-Cox transformation of oocyst shedding was used to evaluate association between oocyst shedding and SNPs, including candidate genes analyzed in previous studies.

Methods

Chickens from coccidiosis resistant and susceptible commercial broiler lines were crossed to produce F₁ chickens [3]. Twelve pairs of F₁ individuals were mated to produce 290 F2 offspring from four hatches. F₂ chickens were orally infected with 1.0 x 10⁴ of sporulated E. maxima oocysts at 4 weeks of age [3]. Body weight was measured on days 3, 6, and 9 post-infection (PI) as previously described [7,8]. The number of oocysts was evaluated using the fecal sample collected between Days 5 to 9 PI by method of Lillehoj and Ruff (1987).

Genomic DNA was extracted from erythrocytes using the GenElute Blood Genomic DNA kit (Sigma, St. Louis, MO). PCR product was sequenced with internal primers using an ABI 3730 DNA analyzer (Applied Biosystem, Foster City, CA). SNP was genotyped by genomic DNA sequencing with the same gene-specific internal primers [5,6]. Then, haplotype phase was decided based on pedigree information. Associations between SNP genotypes or haplotypes and disease phenotypes were evaluated using linear model including family, sex, and hatch effects [5]. The number of oocyst was transformed by log transformation, and it was also transformed using an empirical approach suggested by Box-Cox [9]. These methods use the profile likelihood function for the largest linear model to be considered as a guide in choosing a value for parameter [10].

Results

Ten SNPs were identified in the zyxin, 12 in the TCR-b and 4 in the MLF2 gene [5,6]. At various times following experimental infection of the F₂ generation with Eimeria maxima, body weights, fecal oocyst shedding, and biochemical parameters were measured as parameters of coccidiosis resistance [3]. Single marker and haplotype-based tests were applied to determine the associations between SNPs and the parameters of coccidiosis resistance including body weight and Box-Cox transformed oocyst shedding. The maximum additive genetic effect on disease resistance of a SNP in MLF2 was explained by body weight (p = 0.0002) and this SNP of MLF2 also significantly associated with body weight was also associated with fecal oocyst shedding [6]. This was confirmed using Box-Cox transformed oocyst shedding data in this study. Box-Cox transformed fecal oocyst shedding was close to normal distribution (Kolmogorov-Smirnov normality test, p > 0.1) compared with log transformed data (p > 0.01). Associations between Box-Cox transformed oocyst shedding number and SNP was summarized in Table 1. This analysis support evidence of association between markers in MLF2 gene and oocyst shedding. Association between transformed oocyst shedding and the other genes including zyxin /TCR-b were also detected (p < 0.05) using BoxCox transformed oocyst shedding, which was not significant using log transformed data. The significance level of association between the other markers and oocyst shedding was not changed obviously at various thresholds (p = 0.05, 0.01 or 0.001) compared with analysis using log transformed oocyst shedding.

Identification of parental line specific effects of MLF2 on resistance to coccidiosis in chickens - Image 1

Four haplotypes accounted for 98% of all observed MLF2 SNP (Table 2). The haplotype association test was used to determine the relationship between the oocyst shedding number and the haplotypes of the MLF2. The allele substitution effect of MLF2 haplotype 4 versus haplotypes 1 and 3 was significantly associated with increased body weight and oocyst shedding number (Table 2). The MLF2 homozygous haplotypes 2/2 and 4/ 4 were found in 2 and 4 animals, respectively (Tables 3, 4). However, homozygous haplotypes 1 and 3 were not observed in 24 parent chickens, which implies that only haplotypes 1 (n=4) and 3 (n=10) were possibly originated from a specific parental chicken line. The haplotype substitution effect of haplotype 3 was higher than allele substitution effect of haplotypes 1 and 2 (Tables 3, 4). Haplotypes 2 and 4 were found in both parental line, but haplotypes 1 and 3 were likely to be originated from specific parental line. However, haplotypes 1 and 3 were likely to be originated from different parental chicken lines considering high heterozygous frequency of haplotypes 1 and 3 (Table 3).

Identification of parental line specific effects of MLF2 on resistance to coccidiosis in chickens - Image 2

Identification of parental line specific effects of MLF2 on resistance to coccidiosis in chickens - Image 3

Identification of parental line specific effects of MLF2 on resistance to coccidiosis in chickens - Image 4

Conclusions

Single marker association analysis for 16 traits of coccidiosis resistance showed SNPs of TCR-b and MLF2 were associated with oocyst shedding and body weights in previous study [6], and it was confirmed by additional association tests in current study. In this study allele substitution effect analysis and classification of MLF2 haplotype elucidated parental origin of haplotype associated with coccidiosis resistance. The fourth [A/G] SNP in MLF2 was major determinant of haplotype 3 vs the other haplotypes, but allele A of the first SNP [A/G] was common in haplotypes 1, 2 and 3. Haplotype association test was applied to estimate effect of each haplotype which was not revealed clearly by single marker association tests. Difference between additive effect (body weight and oocyst) of haplotypes 1 and 3 was smaller than that of haplotypes 2 and 3. However, haplotypes 1 and 3 are likely to be originated from the different parental lines. In genome-wide [3] and fine mapping [4] of coccidiosis resistance, QTL was not detected in some chicken families. These previous linkage QTL mapping studies could be affected by heterogeneity in parental lines. The haplotype (ATTG) of MLF2 originated from a specific parental line was associated with elevated body weight and oocyst shedding number. Identifying haplotype diversity of MLF2 in other chicken families will provide useful information for experimental designs in the future studies.

This article was originally published in BMC Proceedings 2011, 5(Suppl 4):S21 http://www.biomedcentral.com/1753-6561/5/S4/S21. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0).