Detection of a Cis Eqtl Controlling Bmco1 Gene Expression Leads to the Identification of a Qtg for Chicken Breast Meat Color

Classical quantitative trait loci (QTL) analysis and gene expression QTL (eQTL) were combined to identify the causal gene (or QTG) underlying a highly significant QTL controlling the variation of breast meat color in a F2 cross between divergent high-growth (HG) and low-growth (LG) chicken lines. Within this meat quality QTL, BCMO1 (Accession number GenBank: AJ271386), encoding the b-carotene 15, 159-monooxygenase, a key enzyme in the conversion of b-carotene into colorless retinal, was a good functional candidate. Analysis of the abundance of BCMO1 mRNA in breast muscle of the HG x LG F2 population allowed for the identification of a strong cis eQTL. Moreover, reevaluation of the color QTL taking BCMO1 mRNA levels as a covariate indicated that BCMO1 mRNA levels entirely explained the variations in meat color. Two fully-linked single nucleotide polymorphisms (SNP) located within the proximal promoter of BCMO1 gene were identified. Haplotype substitution resulted in a marked difference in BCMO1 promoter activity in vitro. The association study in the F2 population revealed a threefold difference in BCMO1 expression leading to a difference of 1 standard deviation in yellow color between the homozygous birds at this haplotype. This difference in meat yellow color was fully consistent with the difference in carotenoid content (i.e. lutein and zeaxanthin) evidenced between the two alternative haplotypes. A significant association between the haplotype, the level of BCMO1 expression and the yellow color of the meat was also recovered in an unrelated commercial broiler population. The mutation could be of economic importance for poultry production by making possible a gene-assisted selection for color, a determining aspect of meat quality. Moreover, this natural genetic diversity constitutes a new model for the study of b-carotene metabolism which may act upon diverse biological processes as precursor of the vitamin A. Copyright: ß 2011 Le Bihan-Duval et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: The creation of the original F2 resource population from the HG x LG intercross and the subsequent quantitative trait loci (QTL) analysis were funded by grants from the United States Department of Agriculture-Initiative for Future Agriculture and Food Systems Animal Genome Program [Grant # 00-52100-9614 to LAC], AGENAE (France) and INRA Departments of Animal Genetics and Animal Physiology (France). The expression QTL study was funded by a …


Introduction
For more than half of a century, commercial poultry breeding programs have focused mainly on improvements of two major production traits, growth rate and feed efficiency, in meat-type (broiler) chickens. Furthermore, different experimental lines of chickens have been created to increase our understanding of genetic control over other important production traits, including meat quality. Our unique model is a population of meat-type chickens that was divergently selected for high (HG) or low growth (LG) rate, based on a difference in body weight (BW) at both 8 and 36 weeks of age [1]. A genetic analysis of the highly heritable growth curve from this experimental selection has been described in detail [2,3]. The HG and LG broiler lines have been extensively studied to understand the physiological and genetic basis of marked differences in growth rate and skeletal muscle development [4,5]. An increase in fiber diameter and at a less extent in the total number of muscle fibers accounts for the greater breast and leg muscle weights of the HG birds [6]. Recently, we found that the HG chickens exhibit a paler meat characterized by higher lightness (BCo-L), lower redness (BCo-R) and yellowness (BCo-Y) than that of LG birds. Several QTL for meat quality were detected in the F2 resource population created from the HG x LG intercross, among these was a strong QTL on GGA11 that affects color of breast meat (i.e., yellowness and redness) [7]. Here, we describe the identification through a cis eQTL analysis of the causal gene and regulatory mutation underlying this meat quality QTL.

Refinement of the Meat Color QTL Region
At the step of initial mapping, the position of the meat color QTL was imprecise, due to the lack of microsatellite markers on the distal end of GGA11. To refine the QTL location, three additional markers were developed and the QTL re-analyzed. As shown in Figure 1, the QTL for the breast muscle yellowness (BCo-Y) trait was detected at about the same position as with the initial mapping (at 68 cM vs. 69 cM), although the significance level was higher (F statistic of 95 vs. 63). In addition, the confidence interval of the QTL was reduced from 35 (13.3-21.9 Mb) to 17 cM . The origin of the high allele for BCo-Y was traced back to the LG line, which was consistent with the more intense yellow color of breast meat in this genotype. The QTL on GGA11 explained about 14% of the phenotypic variance in the BCo-Y trait in the HG x LG F2 population. Half-sib analysis showed that all F1 sires were heterozygous at this QTL, which suggested fixation of alternative alleles in the founder lines originally created from a mixing of several commercial breeds. Similar improvements in refinement of the QTL were obtained for the breast muscle redness (BCo-R), although the level of significance for this trait remained lower. The QTL region for the breast muscle redness trait (BCo-R) overlapped that for yellowness trait (BCo-Y). These observations suggested a strong impact of the GGA11 QTL on the variation of meat color observed between the HG and LG lines.

Identification through eQTL Analysis of the Causal Gene
Refinement of this meat quality QTL region reduced the number of positional candidate genes from about 140 in the initial mapping to 30 genes. The most obvious candidate gene was BCMO1, being located on GGA11 between 17.085 and 17.105 Mbp, the approximate location of the observed BCo-Y QTL ( Figure 1). Furthermore, BCMO1 is a good functional candidate because it encodes b-carotene 15, 159-monooxygenase, an enzyme responsible for the conversion of b-carotene (a yellow pigment) into two colorless retinal (pro-vitamin A) molecules [8]. We first compared BCMO1 mRNA levels in the breast muscle of HG and LG birds across six ages (1-11 wk). As reported in Figure 2, the level of BCMO1 mRNA was consistently higher in the Pectoralis major muscle of HG chickens compared to LG chickens, regardless of age. This large difference in abundance of BCMO1 transcripts between LG and HG birds was evident with or without normalization to 18S ribosomal RNA levels. We then examined the relationship between variations of BCMO1 mRNA levels and the yellowness of breast meat (BCo-Y) in the segregating HG x LG F2 population. Breast muscle BCMO1 mRNA levels were quantified in one of the five F1 half-sib families (n = 134). A significant negative correlation (r = 20.47) between BCMO1 mRNA levels and the BCo-Y measurements was observed. Subsequently, a QTL analysis was performed on GGA11 using BCMO1 mRNA levels as a quantitative trait. As shown in Figure 3a, a strong cis expression QTL (eQTL) was detected that was even more significant than the original BCo-Y QTL ( Fig. 3b; F value of 79 vs. 14). The high allele for the eQTL (BCMO1 mRNA) was traced back to the HG line, indicating that higher expression of the BCMO1 gene is linked to lower meat yellowness in HG birds. The QTL analysis for the BCo-Y trait was also completed using breast muscle BCMO1 mRNA levels as a covariate. As illustrated by Figure 3b, the QTL for BCo-Y previously detected on GGA11 disappeared after this correction. Collectively, these observations clearly show that the variation in the BCo-Y trait is linked to the variation of BCMO1 mRNA levels in breast muscle. Gene Sequencing and Genotyping of Putative Causative SNPs in the F2 Population Next, we sequenced the promoter region and 11 exons of the BCMO1 gene to identify mutations underlying variations in BCMO1 transcript levels. The causal mutation was expected to be heterozygous in the five F1 sires, which were all heterozygous for BCo-Y QTL. Amplified fragments of the BCMO1 gene were sequenced in the five F1 sires from the HG x LG intercross, where a total of 23 SNP and 4 insertion/deletion polymorphisms were identified in the BCMO1 gene and its promoter. Six  . Twenty-one new polymorphisms were observed (accession numbers ss325995734 to ss325995754). Six of the 23 SNP were heterozygous for all F1 sires: one semiconservative SNP (Glu.Lys, AA residue 524) was within exon 11, without any supposed effect on BCMO1 mRNA level, two SNP were within intron 8, one within intron 9 and, interestingly, two SNPs within the promoter region. Since the QTL controlled the BCMO1 gene expression and since the promoter is known as a region regulating expression, the last two SNPs positioned at 2678 (SNP1) and 2621 (SNP2) upstream of ATG ( Figure S1) were considered as the most likely candidate mutations. A genotyping assay, based on high-resolution melting (HRM) curve analysis, was developed to test the association between the SNP1 and SNP2 mutations and variations in BCMO1 mRNA levels and meat color traits (BCo-Y and BCo-R) in the HG x LG F2 population. The two SNPs in the BCMO1 promoter were fully linked defining two AN 57 A and GN 57 G haplotypes. No recombinant AN 57 G or GN 57 A haplotype was observed in any of the F0, F1 or F2 birds. Genotyping of 35 F0 birds showed a difference in allele frequencies between the two lines. The GN 57 G haplotype was fixed in F0 LG birds (N = 20) while the AN 57 A haplotype was predominant in F0 HG birds (9 of 13 birds tested had the AN 57 A type). As shown in Figure 4a

Impact of the haplotype on breast muscle carotenoid contents and relationship with breast meat color
The impact of the haplotype on breast muscle carotenoid content was evaluated. Only lutein and zeaxanthin pigments were detected in breast muscle, while b-carotene was undetectable. Carotenoid content was markedly higher in GN 57 G (n = 11) than in AN 57 A (n = 9) with a ratio of 1.9 and 1.8 for lutein and zeaxanthin, respectively ( Figure 5). Across all samples (n = 20) a significant positive relationship was observed between the total content in carotenoids (lutein + zeaxanthin) and the yellow color of breast meat (Pearson correlation, R = +0.6).

Screening and effect of the mutation in a commercial chicken population
Five out of the ten sires from a commercial population of slowgrowing broilers used for ''Label Rouge'' production were found to be heterozygous for the BCMO1 haplotype. As for the experimental population, the association study performed in the commercial population showed a significant effect of the haplotype on the color of the breast meat, with differences of 0.5 unit of b* (p = 0.007) and 0.4 unit of a* (p = 0.006) between the homozygous birds (LSM of 8.5 and 8.0 for the b* color of Figure 2. Relative BCMO1 mRNA levels in the HG and LG lines according to age. The relative BCMO1 mRNA levels (2 (dCT BCMO1-dCT 18S) ) have been measured in the Pectoralis major muscle of HG and LG chickens slaughtered at 1, 3, 5, 7, 9, and 11 weeks of age (n = 6 by line/age). Age and line effects were tested by two-way analysis of variance. A highly significant effect of the line (P,0.0001) was evidenced, without any interaction with age. doi:10.1371/journal.pone.0014825.g002 the homozygous GN 57 G and AN 57 A respectively; LSM of 21.3 and 21.7 for the a* color of the homozygous GN 57 G and AN 57 A respectively). The effect of the haplotype on BCMO1 expression was evaluated in a sample of 11 AN 57 A and 12 GN 57 G homozygous birds. BCMO1 mRNA levels were threefold higher while BCo-Y values were significantly lower in AN 57 A compared to GN 57 G ( Figure 6).

The Two Fully Linked SNPs have a Functional Effect on the Promoter Activity
The functional effect of the wild type (AN 57 A and GN 57 G) and recombinant (GN 57 A and AN 57 G) haplotypes on the promoter activity was further investigated using a gene reporter strategy. As the AN 57 A, GN 57 A and AN 57 G haplotypes were obtained by mutagenesis of the GN 57 G haplotype, the gene expression variation can directly be associated to the two targeted SNPs. Chicken hepatoma cells (LMH) were transfected with a luciferase reporter construction containing a 2734 to +97 fragment of the chicken BCMO1 promoter, corresponding to a unique sequence except for the two candidate SNPs. The AN 57 A construction was associated with a 3.7 fold higher luciferase activity than the GN 57 G construction (Figure 7). Consistent results were obtained using human hepatoma cell lines (data not shown). For both GN 57 A and AN 57 G constructions, the luciferase activity was similar to the AN 57 A haplotype but significantly higher than the one observed with the GN 57 G haplotype (Figure 7). These results were confirmed in three independent experiments and demonstrated that both SNP are required to affect the activity of BCMO1 promoter and consequently the level of BCMO1 transcripts.

Discussion
While several hundreds QTL have been reported in livestock species including chicken (http://www.animalgenome.org/ QTLdb), the identification of the underlying genes and mutations is still exceptional [9,10]. This study is a good example of the strength of combining positional and functional data since the eQTL approach validated the choice of a candidate gene and rapidly gave strong indications about the underlying mutation(s). Our study supports the conclusion that BCMO1 is the gene underlying the chicken meat color QTL on GGA11, through a cis eQTL controlling its expression. Two fully linked SNPs in the BCMO1 promoter defined two haplotypes (AN 57 A and GN 57 G) with respectively high and low promoter activity as proven in vivo and in vitro. Genomatix was used to check if the mutations could alter a transcription factor binding site. If no results were obtained for the first mutation, the second mutation alters a Vitamin D receptor (VDR), a DR5 binding site detected only with the A allele. However, we clearly demonstrated that both mutations are needed to modify the gene expression. This could be due either to a transcription factor binding site undetected in the first mutation region or most likely to an interaction between the two mutations, one hypothesis being that the first mutation alters the genome 3D structure and affects the fixation of transcription factors in this region. These variations are likely to impact on the biological activity of the BCMO1 enzyme as they lead to respectively low or high carotenoid content in breast meat. Of the major food quality attributes, color is the first critical one since it is the main visual factor involved in the consumer's choice [11]. The identified BCMO1 polymorphism mainly affected the yellow color of breast meat, by altering the content of the yellow pigments lutein and zeaxanthin. These two carotenoids, although not considered as substrates of the BCMO1 enzyme [12], could be converted into b-carotene before symmetric cleavage by BCMO1, as already described for related xanthophylls (canthaxanthin and astaxanthin) in fish [13] and in rat [14]. The effect of BCMO1 polymorphism on red color was less marked, likely because redness may be influenced by other mechanisms such as muscle blood content at death. This last parameter depends on several factors such as genotype [15], bird's activity on the shackle line [16], stunning and bleeding efficiency [17]. The economic value of the haplotype presently identified for meat-type chicken selection will have to be studied more extensively by estimating the frequency of alternative alleles in commercial populations and their effect on color and other important economic traits. Any significant effect of the haplotype on body weight as well as breast or abdominal fatness percentage was observed in the F2 experimental cross which suggested that the effect of the mutation was specific to meat color. Fast genetic tests for the two SNPs by High Resolution Melting method have been developed (for an illustration of the results obtained for SNP2 see Figure S2) and make the association studies possible in various populations. Interestingly, the genotyping of a first commercial population showed that this polymorphism was not specific to the divergent HG and LG lines. The effect of the haplotype on the color of the meat was also significant in this commercial line.
Several mutations have been found in the coding region of BCMO1 gene in humans. One, most likely very rare, leads to a drastic enzyme loss-of-function (90%) and hypercarotenemia and hypovitaminosis A [18]. Others, more frequently distributed also reduce enzyme activity but to a lower extent [19]. BCMO1 knockout in mice had no evident effects on growth but revealed, as expected, a dramatic decrease of vitamin A level in several tissues [20]. Liver steatosis developed without an increase in liver fatty acid synthase, independently of dietary vitamin A levels. Serum free fatty acid levels increased in BCMO1 2/2 mice, whereas no changes in serum glucose level were noticed. In response to dietary oil supplementation (30%), BCMO1 2/2 mice became heavier. In our chicken model, in addition to be of lower weight, LG chickens carrying the ''deficient'' BCMO1 promoter gene are leaner and exhibit higher plasma glucose but normal free fatty acid levels as compared to HG chickens [21]. As body weight and metabolic parameters are under the control of many loci that are segregating between the two divergent HG and LG lines [21], further investigations using BCMO1 mutant chickens with a common genetic background will be needed to precise the specific effects of that gene and to know whether chicken and mice respond differently to variations in BCMO1 activity.
The color of chicken skin, easily visible at the level of the shank, is also dependent upon dietary carotenoid levels and carotenoid metabolism. Interestingly, a recent study [22] showed that the yellow skin phenotype is caused by mutation(s) in the tissue-specific regulatory elements of BCMO2 gene (formerly called BCDO2) which encodes b-carotene 99, 109-monooxygenase. These mutations prevent expression of the gene in the skin but not in liver. While BCMO1 and BCMO2 have close functions (the first

Animals and Color Measurements
The HG and LG chicken lines have been divergently selected for BW for more than 20 generations [1][2][3]. Characteristics of the F2 experimental cross (N = 698) used for QTL analysis have been Figure 5. Effect of the BCMO1 haplotype on carotenoid pigment content of Pectoralis major muscle. ***, P,0.001 Carotenoid concentration and composition were determined through extraction, purification and chromatographic analysis using UPLC in Pectoralis major muscle samples from populations of chicken with the two linked SNP (n = 11 for GN 57 G and n = 9 for AN 57 A). Lutein and zeaxanthin were the sole carotenoids detected in the samples. Their concentrations were 1.9 and 1.8-fold higher for lutein and for zeaxanthin, respectively, in GN 57 G than in AN 57 A chicken (P,0.001). doi:10.1371/journal.pone.0014825.g005 Figure 6. Effect of the BCMO1 haplotype in an unrelated commercial population. Breast meat yellowness (BCo-Y) and relative BCMO1 mRNA level (2 (dCT BCMO1-dCT 18S) ) were measured in the Pectoralis major muscle of birds homozygous at the BCMO1 locus (n = 11 for AN 57 G and n = 12 for GN 57 G). BCo-Y value was significantly lower (***, P,0.001) while BCMO1 mRNA level was higher (***, P,0.001) in AN 57 G than in GN 57 G. doi:10.1371/journal.pone.0014825.g006 described in detail previously [7]. Breast meat color of 9 week old birds was measured 24 h post-slaughter on the upper ventral side of the P. major muscle using a Miniscan Spectrocolorimeter. Color was measured by the CIELAB system, where L* represents lightness, a* redness and b* yellowness. Higher L*, a* and b* values correspond to paler, redder and more yellow meat, respectively.

Genetic Markers
Genomic DNA was extracted from whole blood and QTL mapping on GGA11 was originally developed with three polymorphic microsatellite markers (MCW097, ADL210 and ADL308) [7]. A re-analysis of the meat color QTL on GGA11 was performed by adding three additional markers (SE-QALL0528, SEQALL0529 and MCW230). These markers were already available [23] and were selected for their informativeness among microsatellite markers developed from the chicken genome assembly [21]. Fluorescently-labeled microsatellite markers were analyzed on an ABI 3700 DNA Analyzer (Applied Biosystems, Foster City, CA USA) and genotypes were called using GeneMapper software (Applied Biosystems, Foster City, CA USA). All F0, F1 and F2 (males and females) animals were genotyped with the six microsatellite markers.

QTL Mapping
Prior to QTL detection, the data were corrected for sex and hatch effects as estimated by PEST software [24]. QTL detection, with both F2 and half sib analyses, were conducted by QTL Express, a web based software program [25]. Permutations were done with 10,000 iterations to obtain an empirical threshold at the chromosome level [26]. The 95% confidence interval for a QTL position was estimated by the two LOD (logarithm of the odds) drop-off method [27].

Measurement of BCMO1 mRNA Levels
The expression of BCMO1 gene was first measured in the Pectoralis major muscle of six HG and LG individuals slaughtered at 1, 3, 5, 7, 9 and 11 weeks of age. Then to test the association between SNP mutations and BCMO1 expression, BCMO1 mRNA levels were measured on 373 F2 birds from the HG x LG intercross, including all 134 offspring of the F1 sire family used for eQTL detection and additional 239 birds issued from the other F1 sire families. Methods used for preparation of RNA samples and quantitative real-time RT-PCR assay were described earlier [28]. The chicken BCMO1 sequence (143 bp spanning exon 1 and 2) was amplified using 59-GCCAAGCCATCAAACCAGTG-39 and 59-AACAAAGAAGAGCATCCAGAGCC-39 as reverse and forward primers, respectively. The mRNA levels were estimated based on the cycle threshold (CT) deviation of an unknown sample versus a control cDNA (consisting of a pool of chicken muscle cDNA) according to the equation proposed by Pfaffl [29] and the data were expressed as dCT BCMO1 or dCT BCMO1-dCT 18S.

DNA Sequencing
PCR primers were designed to amplify the 11 BCMO1 exons and the putative promoter region from the chicken genome assembly (http://genome.ucsc.edu/index.html) using Lightcycler Probe Design 2 software (Roche) (See Table S1). PCR amplifications were carried out for each marker and each F1 sire in a 15 ml reaction containing 25 ng DNA, 0.4 mM primers, 0.25 units Taq polymerase (Go Taq, Promega), 1X buffer (Promega), 2 mM MgCl 2 and 0.2 mM dNTP on a GeneAmpH PCR System 9700 thermocycler (Applied Biosystems). The first 5-min denaturation was followed by 38 cycles of: denaturation at 94uC for 30 sec, annealing for 30 sec and elongation at 72uC for 30 sec, with a final elongation of 10 min at 72uC. After Exo/SAP-IT (USB Europe GmbH, Staufen, Germany) purification, PCR products were sequenced on both strands using the same primers and a Big Dye Terminator V3.1 kit (Applied Biosystems) on an ABI 3100 DNA Analyzer (Applied Biosystems). All new polymorphisms have been deposited in dbSNP (http://www.ncbi.nlm.nih. gov/projects/SNP/).

Genotyping
Candidate mutations were screened by high-resolution melting curve (HRM) analysis (see Table S2 for primers and probes). For SNP1 analysis, PCR amplifications were carried out in a 15 ml reaction containing 25 ng DNA, 0.5 mM primer 1, 0.1 mM primer 2, 0.25 units Taq polymerase (Go Taq, Promega), 1X buffer, 2 mM MgCl 2 and 0.2 mM dNTP on a GeneAmpH PCR System 9700 thermocycler (Applied Biosystems). The first 5-min denaturation was followed by 45 cycles, each of denaturation at 95uC for 20 sec, annealing at 58uC for 20 sec and elongation at 72uC for 20 sec, with a final elongation for 10 min at 72uC. After addition of 0.4 mM of allele-specific probe and 1X LC Green (Idaho Technology) to 8 mL of PCR products, the samples were heated to 95uC for 10 min and then progressively cooled to 60uC by 5uC increments and then to 40uC by 10uC increments, 10 sec for each step during, with a ramp rate of 2.2uC/sec, on a LightCycler 480 (Roche). Then the mix was heated to 95uC (0.02uC/sec). The fluorescence data were subsequently analyzed by LightCycler 480 Gene scanning software. SNP2 analysis was performed similarly, except probe and LightCycler Green were added to the PCR reaction mix, and the amplification and HRM analyses were realized on the LightCycler 480: the first 10-min denaturation was followed by 50 cycles, each of denaturation at 95uC for 1 s, annealing at 58uC for 1 sec and elongation at 72uC for 10 sec; then the PCR products were heated to 95uC for 1 sec, and the dissociation curve was obtained by cooling to 40uC and heating to 95uC at a ramp rate of 0.02uC/sec.

Muscle carotenoid measurements
Muscle carotenoid contents were compared between homozygous birds from the F2 population (N = 9 and 11 for AN 57 A and GN 57 G, respectively). Birds from the two genotypes were chosen within the same sire family and exhibited significant differences in BCMO1 expression (1.7 vs -1.3 on the dCT scale for AN 57 A and GN 57 G, respectively) and yellow breast meat color (9.5 vs 12.7 for AN 57 A and GN 57 G, respectively). Carotenoid concentrations were determined in Pectoralis major muscle by using Ultra Performance Liquid Chromatography (UPLC) equipped with a 15062.1 mm Acquity UPLC HSS T3, 1.8 mm column and a 2996 PDA detector (Waters, Saint-Quentin-en-Yvelines, France) according to an initial procedure used for plasma by Chauveau-Duriot et al. [30] adapted to muscle samples. Carotenoids were extracted by hexane (264 ml) from around 0.8 g of muscle sample after ethyl alcohol deproteinisation. Sample treatments were performed at room temperature under yellow light, with echinenone as the internal standard. Lutein, zeaxanthin and all-trans b-carotene were detected at 450 nm and identified by comparison of retention times and spectral analyses with those of pure standards (.95%) purchased from Carotenature (Lupsingen, Switzerland). Concentrations of carotenoids were calculated by using an external standard curve and then were adjusted by percent recovery of the added internal standard.

Study of an unrelated commercial population
Birds originated from a slow-growing commercial line selected by a breeding company since 1994. A large pedigree population was constituted with 868 male and female birds issued from 10 sires and 100 dams. Birds were reared in three successive batches under similar conditions of free-range chickens with access to outdoor area after 6 weeks of age. At 12 weeks of age, birds were slaughtered after a 7 hours feed withdrawal. As for the experimental population, color was measured on the upper ventral side of the P. major muscle by using a Miniscan spectrocolorimeter. It was measured in the CIELAB trichromatic system by Lightness (L*), redness (a*) and yellowness (b*). All the sires and dams were genotyped as described previously. In a second step, genotyping was performed on the offspring of the heterozygous sires (N = 470). BCMO1 gene expression was measured in a sub-population of 11 AN 57 A and 12 GN 57 G homozygous male birds chosen within one batch of hatching.

Plasmid Constructions
Chicken BCMO1 promoter with the GN 57 G haplotype was amplified by PCR from bases -734 to + 97 relative to the transcription start site from exon 1, using the following primers: BCMOF: 59-AGCCTGTGATTTCCTCTGC-39 and BCMO1R: 59-TCGCTCTCCTGTGTCCTCTA-39. First, this fragment was subcloned into PCR3.1 plasmid (Invitrogen, Carlsbad, CA) by TA cloning procedure. The promoter insert of PCR3.1 plasmid was then cloned into pGL2 basic plasmid (Promega, Madison, WI), upstream of the luciferase reporter gene, to create pGL2-BCMO1-Luc plasmid. Site-directed mutagenesis were performed with the QuikChange II Site-Directed Mutagenesis Kit (Stratagene) using the following mutagenic primers (only forward primers are shown), 59-AATCCCCCTTCTTCTTAATTATTCCTACCATACTTT-CGCAGAGG-39 and 59-GTGGTGGAATCAGATCACACG-TAGGAGAATGTAGAATAGTGCC-39 (the mutated bases are in bold) to create the AN 57 A, AN 57 G and GN 57 A haplotypes.
Cell culture and Transactivation Assay. The chicken hepatocellular carcinoma cell line (LMH) was obtained from ATCC (cat # CRL-2117). LMH cells were grown in Waymouth's medium (Invitrogen, Carlsbad, CA) supplemented with 10% foetal calf serum (FCS, Invitrogen, Carlsbad, CA), 100 IU/ml penicillin (Invitrogen, Carlsbad, CA), and 100 mg/ml streptomycin (Invitrogen, Carlsbad, CA). The human hepatocellular carcinoma cells (HuH-7) were provided by the unit 6061-INSERM-France. HuH-7 cells were grown in William's medium (Invitrogen, Carlsbad, CA), 10% FCS, 1.6 ml BSA 30%, 500 nM Hydroxy-Cholesterol (GE Healthcare, Piscataway, NJ), 1 mg/mL Insuline (SIGMA, St. Louis, MO). All cultures were incubated at 37uC in a humidified atmosphere containing 5% CO 2 . For transactivation assay, transient transfections of LMH cells were performed in triplicate in 6 well plates using the FuGENE lipofection protocol (Roche, Neuilly sur Seine, France). Each well was transfected with 200 ng of pGL2-BCMO1-Luc reporter plasmid and 20 ng of pCMV-bgal reference plasmid containing a bacterial b-galactosidase gene. After 24 h following transfection, the cells were washed once with phosphate buffered saline, and incubated with fresh medium. After 48 h following transfection, the cells were lysed (Lysis buffer -Promega, Madison, WI). Cell lysates were stored at -80uC before analyses. After centrifugation 5 min at 12000 g, supernatants (10 ml) were incubated 5 s in the presence of 40 ml of luciferase assay buffer (Promega, Madison, WI) including luciferin (470 mM). Luciferase activity was determined using a luminometer (Veritas Turner Biosystems, Sunnyvale, XA). The b-galactosidase activity was measured by hydrolysation of Ortho Nitro Phenyl Galactopyranoside (ONPG, SIGMA, St. Louis, MO) and used to normalize the luciferase activity. Figure S1 Polymorphisms within the BCMO1 promoter sequence (-817 to +41 bp). All the SNPs found in at least one HG x LG F1 sire are italicized (light red font), while the two linked candidate SNP are in red bold font and underlined (SNP1: CTG/ ACC, SNP2: GTG/AGG). The ATG codon is shown as a shaded box. Found at: doi:10.1371/journal.pone.0014825.s001 (0.33 MB TIF) Figure S2 Melting curves obtained for SNP2 (GTG/AGG). The homozygous AN57A are characterized by a fluorescence peak at 63uC (blue line), the homozygous GN57G by a fluorescence peak at 68uC (red line) and the heterozygous by two fluorescence peaks at 63uC and 68uC (green line). Found at: doi:10.1371/journal.pone.0014825.s002 (0.32 MB TIF)