Using transcriptome profiling to characterize QTL regions on chicken chromosome 5
© Le Mignon et al; licensee BioMed Central Ltd. 2009
Received: 1 August 2009
Accepted: 2 December 2009
Published: 2 December 2009
Although many QTL for various traits have been mapped in livestock, location confidence intervals remain wide that makes difficult the identification of causative mutations. The aim of this study was to test the contribution of microarray data to QTL detection in livestock species. Three different but complementary approaches are proposed to improve characterization of a chicken QTL region for abdominal fatness (AF) previously detected on chromosome 5 (GGA5).
Hepatic transcriptome profiles for 45 offspring of a sire known to be heterozygous for the distal GGA5 AF QTL were obtained using a 20 K chicken oligochip. mRNA levels of 660 genes were correlated with the AF trait. The first approach was to dissect the AF phenotype by identifying animal subgroups according to their 660 transcript profiles. Linkage analysis using some of these subgroups revealed another QTL in the middle of GGA5 and increased the significance of the distal GGA5 AF QTL, thereby refining its localization. The second approach targeted the genes correlated with the AF trait and regulated by the GGA5 AF QTL region. Five of the 660 genes were considered as being controlled either by the AF QTL mutation itself or by a mutation close to it; one having a function related to lipid metabolism (HMGCS1). In addition, a QTL analysis with a multiple trait model combining this 5 gene-set and AF allowed us to refine the QTL region. The third approach was to use these 5 transcriptome profiles to predict the paternal Q versus q AF QTL mutation for each recombinant offspring and then refine the localization of the QTL from 31 cM (100 genes) at a most probable location confidence interval of 7 cM (12 genes) after determining the recombination breakpoints, an interval consistent with the reductions obtained by the two other approaches.
The results showed the feasibility and efficacy of the three strategies used, the first revealing a QTL undetected using the whole population, the second providing functional information about a QTL region through genes related to the trait and controlled by this region (HMGCS1), the third could drastically refine a QTL region.
In spite of success in QTL research for complex traits in livestock species in the last twenty years, location confidence intervals of many QTL are wide, possibly harboring hundreds of genes. This is the major obstacle to finding causative mutations underlying any QTL identified. In addition, fine mapping techniques and positional cloning to reduce the location confidence interval of the initial QTL are time-consuming, especially for livestock species compared to plant and animal models. This is mainly due to a lack of inbred lines, long generation intervals, the cost of maintaining each animal and also the difficulty of producing transgenic or "knock-out" individuals to confirm the causative nature of the mutation of the trait of interest. Few mutations underlying QTL have therefore been identified in livestock (e.g., the DGAT1 gene in dairy cattle , IGF2 gene in swine [2, 3], GDF8 gene in sheep etc, for review see Ron & Weller  and Georges ). Several groups have proposed combining QTL detection programs and high throughput transcriptome data to elucidate biological pathways associated with complex traits and their underlying genetic determinants. [7–14]. This new integrative approach, known as "Genetical Genomics (GG)" or "Integrative Genomics", treats the expression level of each gene present on a microarray as a quantitative trait and use genetic markers to identify genomic regions that regulate gene expression phenotypes. Such regions are named eQTL (expression Quantitative Trait Loci). Independently of the context of QTL identification for a complex trait, the eQTL identification approach was first applied in 2002 by Brem et al.  in order to understand the genetic architecture of natural variations in gene expression in yeast. This approach was soon extended to eukaryotes [10, 11, 15–20]. An eQTL region close to the physical location of a gene controlled by this region is referred to as a cis-eQTL [10, 15]. In such a case, a mutation in the gene itself might be responsible for variability in its own expression at the mRNA level. When an eQTL region for a given gene maps to a location on the genome other than the localization of this gene, it is referred to as a trans e-QTL. Very little is known of the molecular nature of cis-acting and (even more so) trans-acting eQTL regions.
In the context of QTL identification for a complex trait, GG studies have mainly been undertaken on plant or animal models such as flies , mice , rats , eucalyptus , Arabidopsis . GG is not yet usually used because it requires skills in both genetic and genomic fields and the cost of microarray is high, which can be a real limitation for GG in which several animals have to be analyzed. The present study aims at testing the microarray contribution to a QTL research program in livestock species in which the population structure and marker density are less favorable compared to animal models to QTL localization (no consanguineous lines, low marker density...). Our question was therefore whether GG could be transposed in livestock species to reduce a QTL region of interest and to provide new functional information about the causative mutation. To answer this question, this study focused on chicken species, with abdominal fatness as the complex trait. Although various QTL for the fatness trait have been reported in this species , these QTL regions remain wide and no causative mutation has been clearly detected. We chose to apply three complementary strategies integrating transcriptome data in order to improve characterization of a QTL for abdominal fatness (AF) previously detected on the chicken chromosome GGA5 (p < 0.01) with an effect of 1 phenotypic standard deviation . Because of the cost of microarray, the experimental design in the present study included 45 birds. Preliminary linkage analysis on this design revealed the expected AF QTL on the distal GGA5 (p < 0.05) showing that, despite the fairly small size of the experimental design, it is possible to detect QTL with an effect of ~1 phenotypic standard deviation, thus justifying continuing the study. One of the strategies used was based on dissection of the complex trait using the elementary gene expression profiles, as first performed in 2003 by Schadt et al. . To the best to our knowledge, no study using this approach has been published since this first report. The second strategy commonly used by authors working in this context was based on the identification of genes with eQTL co-localizing with the QTL responsible for the complex trait of interest. The function of such genes can provide new functional information about the candidate positional and functional gene sought in the QTL region as causative to the trait of interest. Only one study conducted in livestock species has been reported using this approach . In the present study, we used this strategy in this way in order to characterize the QTL functionally; we also used this strategy in another way using a multiple trait model for QTL analysis in order to refine this QTL region. A third approach was to use hepatic transcriptome profiles to predict for each recombinant offspring the Q versus q GGA5 AF QTL allele (at the mutation looked for) inherited from its sire and then to refine the localization of the GGA5 AF QTL after determining the recombination breakpoints.
Animal design and microarray setup
Previous studies using a three-generation design performed by inter-crossing two experimental chicken lines divergently selected on abdominal fatness have revealed 6 QTL for abdominal fatness (AF) on GGA1, 3, 5 and 7 in male meat-type chickens, [22, 24]. Different recombinant backcross 1 (BC1) and BC2 males were then produced to refine AF QTL on chromosome 5 by crossing F1 sires heterozygous for these QTL with lean line dams [22, 24]. The present study focused on a BC1 chicken sire (and its 71 BC2 male offspring) known to be heterozygous for the AF QTL on the chicken distal GGA5 chromosome, the other QTL on GGA1, 3 and 7 were not detected in this family. Our aim was to integrate hepatic transcriptome data to refine the location confidence interval of this AF QTL on GGA5 and to characterize it functionally. To reduce the cost of microarray experiments, only 46 animals randomly selected for AF values from the 71 BC2 male offspring were studied. Using a 20 K chicken oligo array (Ark-genomics), 20461 gene expression measurements were obtained from the livers of these animals. One microarray was discarded from the 46 because of lower hybridization signals, confirming the good quality of the technical procedures (see Materials and Methods), and leaving a total of 45 birds for further study. Fifty-five percent of the 20461 genes (11213) were selected as expressed in the liver (see Materials and Methods), among which a human ortholog with a HUGO symbol for 4002 genes was determined. The raw and normalized microarray data were deposited in the Gene Expression Omnibus (GEO) public repository . The accession number for the series is GSE12319 and the sample series can be retrieved with accession numbers GSM309564 to GSM309609.
Selection of 660 genes "correlated" with the AF trait
Six hundred sixty genes were found to be "associated" with the AF trait by analysis of their correlation with AF (P < 0.05 at the gene level) or their differential expression between the 10 fattest and 10 leanest birds out of the 45 (P < 0.05 at the gene level). Despite the absence of correction for multiple tests, principal component analysis (PCA) generated with the 660 transcript levels showed appropriate separation between fat and lean chickens on the two principal components, explaining 30% of the data variance (Additional file 1).
Approach 1: Refining GGA5 AF QTL by dissection of the AF complex trait using the 660 gene-set transcriptome profiles
In order to improve the understanding of the characteristics of the subgroups obtained by HCA, we calculated (on the basis of the marker information only) the probability of each offspring receiving from its sire the Q or q haplotype for the two QTL at 102 cM and 175 cM; we considered only birds having a probability > 0.95 (43 and 35 birds for the two QTL at 102 cM and 175 cM, respectively). For the two AF QTL detected, the animal labels on Figure 2A show the distribution of birds for which the AF value was in disagreement with the paternal Q/q haplotype. A large proportion of these birds were included in subgroups whose suppression allowed an increase in power of QTL detection (subgroup 1 for the QTL at 102 cM and subgroups 2 and 3 for QTL at 175 cM). These results may explain the better power of linkage analysis using this "subgroup" strategy.
Approach 2: refining distal GGA5 AF QTL by eQTL mapping
Selection of genes correlated with the AF trait and having an eQTL colocalizing with GGA5 AF QTL
Out of the previously selected 660 genes, we identified 46 genes (6.9%) that had an eQTL (p-value < 0.1) that colocalized with the location confidence interval of the GGA5 AF QTL (156-187 cM). Among the 11213 genes analyzed, 285 genes (2.3%) had an eQTL that mapped in the GGA5 AF QTL. Therefore, using the 660 genes correlated with the AF trait led to a 3-fold increase in genes (6.9% against 2.3%) having an eQTL colocalizing with the GGA5 AF QTL. The same analysis was performed on the other chromosomes (GGA1, GGA3 and GGA7), and no enrichment of genes with an eQTL in a particular region of these chromosomes was observed in this 660 gene-set, probably because of the absence of other AF QTL regions on these chromosomes. As previously observed with the AF QTL analysis, the design of 45 animals allowed us to detect significant eQTL regions with an effect on an expression trait of about 1 within-family residual standard deviation, with probably some false positives (no correction for multiple tests). Maximum likelihood ratio test (LRT) locations of eQTL for the 46 previously selected genes were evenly distributed over the 156-187 cM location confidence interval of the GGA5 AF QTL. Correlations between gene expressions were variable (from ~0 up to 0.7 and -0.64), suggesting different independent gene networks. These results strongly suggested that this 46 gene-set reflected at the mRNA level the impact of different linked mutations in the location confidence interval of the GGA5 AF QTL.
Selection of the genes regulated by the GGA5 AF QTL mutation or a mutation close to it
To select the closest eQTL mutations to the location of the GGA5 AF QTL mutation, we performed 46 linkage analyses on the residual AF value corrected for each transcript. Genes could be sorted by this step according to the degree of correlation between their transcript levels and AF trait and/or the high proximity of the maximum LRT location for their eQTL with one of the AF QTL. As a result, 12 genes out of the 46 genes effectively corrected the AF QTL (P > 0.1).
Validation by RT-PCR
Because of the high probability of false positives (no correction for multiple tests) we needed to validate these results using another method of mRNA quantification. Expression of the 12 genes was quantified by qRT-PCR. The previous results were confirmed for 5 genes. This 5 gene-set was considered to be the best gene-set giving new information about the position of the causative mutation underlying AF QTL and possibly about the impact at the mRNA level of this mutation for some of them. We therefore used these 5 genes to refine confidence interval of AF QTL location and also analyzed their functions to identify those likely to be regulated by the mutation itself.
Functional analysis of the 5 gene-set
According to the V3.2 annotation of the 20 K chicken oligochip (see Materials and Methods), 3 of the 5 genes of interest had a precise gene name: HGMCS1, TCF3 and SALL4. RIGG10516 and RIGG19646 genes both encode a hypothetical protein. A gene-by-gene bibliography indicated that TCF3 and SALL4 are both transcriptional factors involved in large molecular processes such as B-cell development and pluripotent stem cell generation, respectively. Neither seems to be directly involved in lipid metabolism. However, HGMCS1 (3-Hydroxy-3-MethylGlutaryl-CoA Synthase 1) is clearly involved in lipid metabolism, more precisely in cholesterol synthesis .
Refining GGA5 AF QTL by mapping of the 5 gene-set combined variable
Summary of reduction of GGA5 AF QTL using the 3 approaches
Confidence interval (CI) of GGA5 AF QTL according to the strategy used
Initial AF QTL CI
52.2 - 57 (4.8)
54 - 57 (3)
54.16 - 55.1 (1.04)
Overlapped CI taking into account all approaches
The next step was to use these 5 genes to refine the location of this QTL with another approach (Figure 1).
Approach 3: Refining distal GGA5 AF QTL by prediction of the paternal Q/q GGA5 AF QTL mutation, using the 5 gene-set for the recombinant birds
The first approach consisted of "dissection" of the complex AF trait by grouping the offspring in accordance with their 660 gene transcriptome profiles. Schadt et al.  were the first to use this strategy and they improved the significance of a fatness QTL previously detected on a chromosome and even detected a new QTL on another. We found similar results in the present study. First, we found an increase in the LRT of GGA5 AF QTL by using subgroups 1 and 4. Second, despite the small size of our experimental design (45 animals), we detected using subgroups 2 and 3 an AF QTL at 102 cM on GGA5, as previously detected using an experimental design with 1300 birds . These results show the power of the approach. In view of the polygenic influence on the complex traits, the variations in abdominal fatness are probably due to variations in several biological pathways, impacted by multiple mutations acting separately or in interaction. Transcriptome data offer the possibility of dissecting such a complex trait in more elementary phenotypes (gene expressions correlated with the trait) and therefore make it possible to separate the population into genetically homogenous subgroups using transcriptome profiling. Some combinations of subgroups possibly reflect the effects of a precise mutation whereas others reflect the signature of other mutations. In this study, we showed the over representation in some subgroups of birds for which the AF value was in disagreement with the paternal Q/q haplotype of the two GGA5 AF QTL, thus increasing the QTL detection power when they were removed. These results indicate that the QTL each significantly affect only a subset of the population analyzed. This heterogeneity observed between offspring clearly demonstrates the complexity underlying traits such as fatness. In summary, the 660 genes correlated with the AF allowed classification of offspring in a relevant way to dissect AF QTL on the GGA5 chromosome, despite the small number of animals analyzed.
To the best of our knowledge, no study using this approach has been published since the first publication in 2003 . Our results obtained with a livestock design and those of Schadt et al. obtained with mice indicate that the identification of subgroups in a population on the basis of transcriptome profiles would be an effective way to improve the power of QTL detection by linkage analysis.
The second strategy aimed at improving characterization of GGA5 AF QTL by eQTL mapping. Such a strategy is widely used in the context of QTL detection of complex traits using transcriptome data [10, 11, 19, 23, 30–37]. The principle is to identify genes correlated with the complex trait that have an eQTL co-localizing with the QTL of interest. Most of these authors then focused on genes with a function related to the complex trait and having a cis-eQTL, allowing them to hypothesize that the mutation responsible for the complex trait is in the cis-eQTL gene [10, 30, 32, 34–37]. In our study, no cis-eQTL was detected among the 46 genes correlated with the AF trait and having an eQTL in the GGA5 AF QTL region. Correlation pairs between these genes strongly suggested that these 46 trans-eQTL gene expressions were probably controlled by different mutations in the location confidence interval of the GGA5 AF QTL, as previously commented by Georges  and Schadt et al. . Because the CI of the QTL was large (31 cM), we therefore selected 5 genes that effectively corrected the AF QTL (P > 0.1), and therefore likely to be controlled by the mutation sought or by a mutation close to it. Functional analysis of these 5 genes was still limited by the partial functional annotation of genes. However, we identified one gene related to lipid metabolism that could be affected by the mutation sought. This gene encodes HGMCS1, known to be involved in cholesterol metabolism. We have recently shown that its regulation in response to fasting is different in chickens compared to mammals . Further experiments will be necessary to clarify its role in fatty acid metabolism and its regulation in chickens in order to target a potential regulatory gene in the distal GGA5 AF QTL.
Moreover, the whole 5 gene-set considered as a signature of the mutation underlying the QTL of interest or a mutation close to it may be useful to refine this QTL using a multivariate model that takes advantage of the correlation between these 5 expression traits and AF. Multivariate analysis combining the CV5 variable and the AF trait led to a significant increase in maximum LRT compared to the AF trait. This result supported the hypotheses of the existence of QTL affecting both AF and CV5 at the same position or the existence of different close mutations in linkage disequilibrium. We were unable to reach a conclusion with such a small number of animals analyzed. However, this result makes it possible to reduce the location confidence interval of GGA5 AF QTL from 156-187 cM (31 cM) to 166-184 cM (18 cM).
Finally, an original approach to refine a QTL region was proposed in this study. We used the same 5 gene-set to find the best gene expression combination discriminating the paternal Q from q haplotypes (corresponding to the whole confidence interval of the GGA5 AF QTL) and used it to predict the Q versus q mutation received by the recombinant animals. Genotyping these birds with additional markers drastically reduced the region to 166-173 cM (7 cM). Contrary to conventional approaches used to refine a QTL, such a strategy avoided generating new offspring to test the QTL genotype of the recombinant birds and saved on high levels of genotyping, thus gaining time and saving money. However, it is important to remember that it was based on the relevance of the gene-set considered as the signature of the QTL mutation or mutations close to it.
The gains in power and precision of QTL detection offered by approaches 1 and 2 were probably limited by the low density of markers and size of the experimental design used in this study (45 birds). Nevertheless, these approaches allowed substantial reduction of the GGA5 AF QTL region (20% up to 50%). Approach 3 was more effective (80% reduction), depending on the recombination breakpoints in the recombinant birds. This third approach allowed us to refine the GGA5 AF QTL from 156-187 cM (31 cM) to a most probable location confidence interval of 166-173 cM (7 cM) (Table 1). It can be seen that this reduction of GGA5 AF QTL region was consistent with the other more limited reductions obtained by the other two approaches. Unfortunately, a gene by gene bibliography analysis did not allow us to propose a good functional and positional gene candidate as regulator of HMGCS1 and the AF trait.
Our results showed the value of using "Genetical Genomics (GG)" to characterize QTL responsible for complex trait variability in livestock. The originality of this study was to propose complementary approaches allowing a reduction of a QTL region and also providing functional information on it. The most common way to use GG in the QTL detection field is to identify an eQTL region colocalizing with a QTL of interest, making it possible to propose candidate genes possibly regulated by the QTL mutation. In the present study, we identified HGMCS1 that could be affected by the GGA5 AF QTL. In addition, as previously reported by Schadt et al.  using a different animal design, we showed that the identification of animal subgroups on the basis of their transcriptome profiles is an effective way to partially eclipse the polygenic effects which interfere in classical QTL analyses. Such an approach improves the precision of previously detected QTL and also localizes new ones. Finally, the original procedure, consisting of predicting QTL mutation allele for recombinant animals using a haplotype signature based on transcriptome profiles, may lead to drastic reduction of the QTL region. Detection of causative gene mutations underlying the GGA5 AF QTL will be further studied by increasing marker density in the QTL region and the number of animals analyzed, the number of recombinants being crucial in QTL detection analyses. Because of the availability in the future of high marker density combined with the drastic price reduction of microarrays, larger eQTL experimental designs are expected in livestock and this should accelerate identification of causative genes responsible for economic trait variability.
Animals and experimental procedures
All animals were bred at INRA, avian experimental unit UE1295, F-37380 Nouzilly, in accordance with European Union guidelines for animal care and under the authorization 6290 delivered to Nadine Sellier by the French Ministry of Agriculture. A three-generation design was generated by-intercrossing two experimental meat-type chicken lines that were divergently selected on fatness , referred to as fat (FL) and lean (LL) lines. First generation F1 birds were generated from 14 FL males mated with 24 LL females and 4 LL males mated with 8 FL females. A full genome scan in one F2 experimental design bird allowed us to identify several QTL for AF and to refine one of them located on chromosome 5, and 81 backcross (BC1) males were produced from an F1 male (known to be heterozygous for the GGA5 AF QTL) mated with 10 LL females. One of the BC1 animals, recombinant for the paternal chromosome in the QTL region, was mated with 8 LL females to generate 71 male BC2 animals obtained in four hatches. Progeny testing demonstrated that this male was heterozygous for the GGA5 AF QTL and was therefore chosen for transcriptome analysis. Forty-six animals were randomly chosen from the 71 birds to be analyzed by microarray procedures. BC2 chickens were fed ad libitum using a conventional starter diet (0-3 weeks: 12.8 MJ of metabolizable energy) and then a growing broiler diet (4-9 weeks: 13.0 MJ of metabolizable energy). Light/dark periods were 24 h light for the first 2 days, then 14 h light/10 h night up to slaughtering. At 4 weeks of age, blood samples were collected for DNA extraction and genotyping. At 9 weeks of age, the birds were fed ad libitum for a minimum of 4 hours after overnight fasting and then weighed and sacrificed by electrical stunning in the experimental processing plant. Following sacrifice, livers were collected, quickly frozen in liquid nitrogen and stored at -80°C until RNA extraction for transcriptome analyses. After evisceration, carcasses were stored overnight at 4°C before dissection and weighing of abdominal fat.
DNA extraction and marker genotyping
Genomic DNA was extracted from blood samples (100 μl) from the 71 male BC2 according to the modified phenol/chloroform method . DNA was quantified by the saran method  or optical density reading and diluted to 10 ng/μL. Genotyping for chromosomes 1, 3, 5 and 7 was performed for 11, 7, 10 and 5 markers, respectively (Additional Files 4 and 5). Fluorescence-labelled microsatellites were analyzed on an ABI 3100 DNA sequencer (Applied Biosystems, Foster City, CA, USA). Genotypes were interpreted using both the GenoMapper™ software (Applied Biosystems, Foster City, CA, USA) and the GEMMA database . The Kosambi genetic distances in centi-Morgans (cM) were newly estimated with all animals using the "build" option in the CriMap linkage program . Marker order was explored using the FLIPS command until the marker order that maximized the likelihood was obtained.
Total RNA was extracted with TRIzol® reagent (Invitrogen, Cergy Pontoise, France) according to the manufacturer's instructions. Quality and concentration of extracted RNA were assessed using a 2100 Bioanalyzer (Agilent Technologies, Massy, France).
The 20 K chicken array printed in singlets was produced by ARK-Genomics (Roslin Institute - UK: http://www.ark-genomics.org), and the array design was published in the ArrayExpress  repository with Accession N° A-MEXP-820 ArrayExpress 2003 and in the Gene Expression Omnibus with the name GPL5480 . The functional annotations used in the present study (version V3.2) are available on the web site: http://www.sigenae.org. They were obtained by a bioinformatics procedure developed by SIGENAE (INRA-France) .
mRNA labeling, hybridization and data acquisition
All these procedures were as previously described by Desert et al. . Briefly, 5 μg of each mRNA sample were reverse-transcribed and Cy5 fluorescence-labeled. Each Cy5-labeled mRNA sample was hybridized to the microarray with the same Cy3-labeling reference probe according to the procedure of Transcriptome-Biochip Laboratory of Genopole "Toulouse Midi-Pyrénées" (France). The reference RNA pool was made up from equal amounts of RNA derived from all liver samples. The fluorescence ratio for each gene reflected the relative abundance of the mRNA of interest in each experimental sample compared with the same reference mRNA. The reference thus made it possible to take into account any eventual "spot × array" interaction.
Fluorescence signals were detected with a laser scanner (GenePix 4000A from Axon Instrument, CA) keeping a constant PMT gain for each channel. The images were then analyzed with GenepixPro 4.0 software (Axon instruments, Inc., Union City, CA). The raw files were stored in genepix files compatible with the LIMMA library of the R-project statistical and Bioconductor environment  which was used for the normalization and analysis of the data.
The first step was to select the genes considered as expressed in the liver (roughly 60% expected) showing a good contrast between spot and background intensities (SNR (Signal to Noise Ratio) defined as greater than 2). Among the 20461 genes on the microarray, we selected 11590 genes (57%) and eliminated one microarray for which the average SNR was equal to 1.5. We then removed 377 genes for which a minimum of 20% of spots among the 45 microarrays did not respect the 2 criteria: i) absence of the genepix flag (automatically performed by GenepixPro 4.0 ) and ii) symmetry of the spot. A total of 11213 genes were finally selected. Finally we targeted the isolated "bad" spots that did not conform to the two criteria in order to avoid taking them into account during the normalization procedure (2.8% of spots).
The procedure was previously described by Desert et al. . The Cy5/Cy3 ratio used was expressed as the Log2 of the ratio of median pixel intensity of the red and green spots. Median Log2 ratio values were then normalized (ratio centered on zero) according to the hypothesis that the majority of gene expressions did not differ between two samples. The normalization was performed by a nonlinear regression method (Lowess fitness)  to take into account the intensity dependence of the fluorescence bias.
Raw data analyses were performed using a code written in R and softwares from the open-source Bioconductor Project . The animal labels were defined as follows: F1 to F20 for the 20 fattest animals, L1 to L20 for the 20 leanest animals and I for the 5 intermediates. Analysis of variance between the 10 fattest (F1-F10) and 10 leanest chicken (L1-L10) groups for the AF trait value were performed with "aov" function. Two-way Hierarchical Cluster Analysis (HCA) was performed using the "hclust" function with "1-cor" as distance function and "ward" as aggregation criterion; the "heatmap" function was used to generate images. Principal Component Analysis (PCA) was performed with the "pca" function of the FactoMiner library. The predictions of Q versus q allele at the QTL mutation for the recombinants were estimated by discriminant analysis or a logistic regression model using the "lda" and "glm" functions of R, respectively.
Real time quantitative RT-PCR (qRT-PCR) assay
Reverse transcription (RT) was carried out using the high-capacity cDNA archive kit (Applied Biosystems, Foster City, CA) according to the manufacturer's protocol. Briefly, 200 μL of each reaction mixture containing 20 μL of 10× RT buffer, 8 μL of 25× dNTPs, 20 μL of 10× random primers, 10 μL of MultiScribe Reverse Transcriptase (50 U/μL), and total RNA (10 μg) was incubated for 10 min at 25°C followed by 2 h at 37°C. A 1/10 or 1/20 dilution (depending on the gene) of each RT reaction was further used for real time quantitative PCR (qPCR). cDNA samples were mixed with 20 μl ABsolute SYBR Green Mix (Abgene, UK) and 300, 450 or 600 nM (according to the gene) of specific reverse and forward primers. Reaction mixtures were incubated in an iCycler iQ Multicolour Real-Time PCR Detector (Bio-Rad, Marne la Coquette, France) programed to perform one cycle (95°C for 15 min) and 40 cycles (95°C for 15 s and 59°C for 45 s). A melting curve program was then performed for each gene to check the presence of a single product with a specific melting temperature. For each sample and each gene, PCR runs were performed in duplicate. For each gene, serial PCR reactions constructed with 2-fold serial dilutions from a pool of the cDNA samples were systematically added on each microplate for the calibration curve and determination of the amplification rate (R) of the Taq polymerase. For all genes, the amplification rates were within the range 99% to 100% and could be considered as equal to 1. Thus, for the same sample, the gene expression level could be normalized relative to the B-actin expression level.
QTL and eQTL mapping
Before QTL analyses, the AF trait values of the sire family (71 birds) were adjusted for hatch and dam effects by two-way variance analysis, including body weight at slaughter as a covariate (SAS GLM procedure). For the eQTL analyses, no adjustment of the gene variables was performed for hatch and dam effects because of the small size of the population studied (45 birds). QTLMAP software based on an interval mapping method described by Elsen et al. , was used to detect QTL (or eQTL) affecting the AF trait (or a gene expression phenotype). Gene expression values were obtained by microarray or qRT-PCR procedure. The statistical variable for testing the presence of one QTL (or eQTL) versus no QTL (or no eQTL) at one location was an approximate likelihood ratio test (LRT) . Significance thresholds were empirically determined for AF QTL and transcript level eQTL from 2000 simulations. The widely used "one LOD drop-off method" was applied to obtain 95% confidence intervals of the QTL location . QTLMAP software was also used to perform multivariate QTL analysis. As only two traits (CV5 and AF) was studied, we were able to apply a multivariate model with a multinormal penetrance distribution, which is the most powerful and accurate method, even though it is time consuming .
Gene Expression Omnibus
Gallus Gallus, GO: Gene Ontology
Hierarchical Cluster Analysis
Kyoto Encyclopedia of Genes and Genomes
3-Hydroxy-3-MethylGlutaryl-CoA Synthase 1
Human Genome Organization
Likelihood Ratio Test
Principal Component Analysis
quantitative Reverse Transcription-Polymerase Chain Reaction
Expression Quantitative Trait Locus.
GLM is a Ph.D. fellow supported by the French Technical Institute for Poultry (ITAVI). The research program was supported by grants from a French society for genomics in poultry (AGENAVI), INRA and the Agence Nationale de la Recherche (Grant N°0426). Genotyping was performed at Toulouse-Midi-Pyrénées Génopole® (France). The authors are grateful to André Neau (INRA, Jouy en Josas) for archiving of phenotype and genotyping data in the INRA Animal Genetics Department database, Nicole Millet and Nadine Sellier (UE 1295, INRA, Nouzilly) for supervising the maintenance of the LL and FL flocks, Cécile Duby and Frédéric Heraut for RNA extraction, Richard Talbot (ARK Genomics) for providing microarrays and Vincent Le Guyader for his help in programming in R language. The authors thank Doreen Raine for editing the English.
- Grisart B, Farnir F, Karim L, Cambisano N, Kim JJ, Kvasz A, Mni M, Simon P, Frere JM, Coppieters W, et al: Genetic and functional confirmation of the causality of the DGAT1 K232A quantitative trait nucleotide in affecting milk yield and composition. Proc Natl Acad Sci USA. 2004, 101 (8): 2398-2403. 10.1073/pnas.0308518100.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Laere AS, Nguyen M, Braunschweig M, Nezer C, Collette C, Moreau L, Archibald AL, Haley CS, Buys N, Tally M, et al: A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature. 2003, 425 (6960): 832-836. 10.1038/nature02064.View ArticlePubMedGoogle Scholar
- Georges M, Andersson L: Positional identification of structural and regulatory quantitative trait nucleotides in domestic animal species. Cold Spring Harb Symp Quant Biol. 2003, 68: 179-187. 10.1101/sqb.2003.68.179.View ArticlePubMedGoogle Scholar
- Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibe B, Bouix J, Caiment F, Elsen JM, Eychenne F, et al: A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Genet. 2006, 38 (7): 813-818. 10.1038/ng1810.View ArticlePubMedGoogle Scholar
- Ron M, Weller JI: From QTL to QTN identification in livestock--winning by points rather than knock-out: a review. Anim Genet. 2007, 38 (5): 429-439. 10.1111/j.1365-2052.2007.01640.x.View ArticlePubMedGoogle Scholar
- Georges M: Mapping, fine mapping, and molecular dissection of quantitative trait Loci in domestic animals. Annu Rev Genomics Hum Genet. 2007, 8: 131-162. 10.1146/annurev.genom.8.080706.092408.View ArticlePubMedGoogle Scholar
- Mootha VK, Lepage P, Miller K, Bunkenborg J, Reich M, Hjerrild M, Delmonte T, Villeneuve A, Sladek R, Xu F, et al: Identification of a gene causing human cytochrome c oxidase deficiency by integrative genomics. Proc Natl Acad Sci USA. 2003, 100 (2): 605-610. 10.1073/pnas.242716699.PubMed CentralView ArticlePubMedGoogle Scholar
- Schadt EE, Lamb J, Yang X, Zhu J, Edwards S, Guhathakurta D, Sieberts SK, Monks S, Reitman M, Zhang C, et al: An integrative genomics approach to infer causal associations between gene expression and disease. Nat Genet. 2005, 37 (7): 710-717. 10.1038/ng1589.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghazalpour A, Wang X, Lusis AJ, Mehrabian M: Complex inheritance of the 5-lipoxygenase locus influencing atherosclerosis in mice. Genetics. 2006, 173 (2): 943-951. 10.1534/genetics.106.057455.PubMed CentralView ArticlePubMedGoogle Scholar
- Schadt EE, Monks SA, Drake TA, Lusis AJ, Che N, Colinayo V, Ruff TG, Milligan SB, Lamb JR, Cavet G, et al: Genetics of gene expression surveyed in maize, mouse and man. Nature. 2003, 422 (6929): 297-302. 10.1038/nature01434.View ArticlePubMedGoogle Scholar
- Hubner N, Wallace CA, Zimdahl H, Petretto E, Schulz H, Maciver F, Mueller M, Hummel O, Monti J, Zidek V, et al: Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease. Nat Genet. 2005, 37 (3): 243-253. 10.1038/ng1522.View ArticlePubMedGoogle Scholar
- Wayne ML, McIntyre LM: Combining mapping and arraying: An approach to candidate gene identification. Proc Natl Acad Sci USA. 2002, 99 (23): 14903-14906. 10.1073/pnas.222549199.PubMed CentralView ArticlePubMedGoogle Scholar
- Kirst M, Myburg AA, De Leon JP, Kirst ME, Scott J, Sederoff R: Coordinated genetic regulation of growth and lignin revealed by quantitative trait locus analysis of cDNA microarray data in an interspecific backcross of eucalyptus. Plant Physiol. 2004, 135 (4): 2368-2378. 10.1104/pp.103.037960.PubMed CentralView ArticlePubMedGoogle Scholar
- DeCook R, Lall S, Nettleton D, Howell SH: Genetic regulation of gene expression during shoot development in Arabidopsis. Genetics. 2006, 172 (2): 1155-1164. 10.1534/genetics.105.042275.PubMed CentralView ArticlePubMedGoogle Scholar
- Brem RB, Yvert G, Clinton R, Kruglyak L: Genetic dissection of transcriptional regulation in budding yeast. Science. 2002, 296 (5568): 752-755. 10.1126/science.1069516.View ArticlePubMedGoogle Scholar
- Monks SA, Leonardson A, Zhu H, Cundiff P, Pietrusiak P, Edwards S, Phillips JW, Sachs A, Schadt EE: Genetic inheritance of gene expression in human cell lines. Am J Hum Genet. 2004, 75 (6): 1094-1105. 10.1086/426461.PubMed CentralView ArticlePubMedGoogle Scholar
- Yvert G, Brem RB, Whittle J, Akey JM, Foss E, Smith EN, Mackelprang R, Kruglyak L: Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors. Nat Genet. 2003, 35 (1): 57-64. 10.1038/ng1222.View ArticlePubMedGoogle Scholar
- Kirst M, Basten CJ, Myburg AA, Zeng ZB, Sederoff RR: Genetic architecture of transcript-level variation in differentiating xylem of a eucalyptus hybrid. Genetics. 2005, 169 (4): 2295-2303. 10.1534/genetics.104.039198.PubMed CentralView ArticlePubMedGoogle Scholar
- Mehrabian M, Allayee H, Stockton J, Lum PY, Drake TA, Castellani LW, Suh M, Armour C, Edwards S, Lamb J, et al: Integrating genotypic and expression data in a segregating mouse population to identify 5-lipoxygenase as a susceptibility gene for obesity and bone traits. Nat Genet. 2005, 37 (11): 1224-1233. 10.1038/ng1619.View ArticlePubMedGoogle Scholar
- Bystrykh L, Weersing E, Dontje B, Sutton S, Pletcher MT, Wiltshire T, Su AI, Vellenga E, Wang J, Manly KF, et al: Uncovering regulatory pathways that affect hematopoietic stem cell function using 'genetical genomics'. Nat Genet. 2005, 37 (3): 225-232. 10.1038/ng1497.View ArticlePubMedGoogle Scholar
- Abasht B, Dekkers JC, Lamont SJ: Review of quantitative trait loci identified in the chicken. Poult Sci. 2006, 85 (12): 2079-2096.View ArticlePubMedGoogle Scholar
- Lagarrigue S, Pitel F, Carre W, Abasht B, Le Roy P, Neau A, Amigues Y, Sourdioux M, Simon J, Cogburn L, et al: Mapping quantitative trait loci affecting fatness and breast muscle weight in meat-type chicken lines divergently selected on abdominal fatness. Genet Sel Evol. 2006, 38 (1): 85-97. 10.1186/1297-9686-38-1-85.PubMed CentralView ArticlePubMedGoogle Scholar
- Ponsuksili S, Jonas E, Murani E, Phatsara C, Srikanchai T, Walz C, Schwerin M, Schellander K, Wimmers K: Trait correlated expression combined with expression QTL analysis reveals biological pathways and candidate genes affecting water holding capacity of muscle. BMC Genomics. 2008, 9: 367-10.1186/1471-2164-9-367.PubMed CentralView ArticlePubMedGoogle Scholar
- Abasht B, Pitel F, Lagarrigue S, Le Bihan-Duval E, Le Roy P, Demeure O, Vignoles F, Simon J, Cogburn L, Aggrey S, et al: Fatness QTL on chicken chromosome 5 and interaction with sex. Genet Sel Evol. 2006, 38 (3): 297-311. 10.1186/1297-9686-38-3-297.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. 2007, D760-765. 10.1093/nar/gkl887. 35 Database
- Le Mignon G, Pitel F, Gilbert H, Le Bihan-Duval E, Vignoles F, Demeure O, Lagarrigue S, Simon J, Cogburn LA, Aggrey SE, et al: A comprehensive analysis of QTL for abdominal fat and breast muscle weights on chicken chromosome 5 using a multivariate approach. Anim Genet. 2009, 40 (2): 157-164. 10.1111/j.1365-2052.2008.01817.x.View ArticlePubMedGoogle Scholar
- Kattar-Cooley PA, Wang HH, Mende-Mueller LM, Miziorko HM: Avian liver 3-hydroxy-3-methylglutaryl-CoA synthase: distinct genes encode the cholesterogenic and ketogenic isozymes. Arch Biochem Biophys. 1990, 283 (2): 523-529. 10.1016/0003-9861(90)90677-Q.View ArticlePubMedGoogle Scholar
- Gilbert H, Le Roy P: Comparison of three multitrait methods for QTL detection. Genet Sel Evol. 2003, 35 (3): 281-304. 10.1186/1297-9686-35-3-281.PubMed CentralView ArticlePubMedGoogle Scholar
- Lan H, Stoehr JP, Nadler ST, Schueler KL, Yandell BS, Attie AD: Dimension reduction for mapping mRNA abundance as quantitative traits. Genetics. 2003, 164 (4): 1607-1614.PubMed CentralPubMedGoogle Scholar
- Bao L, Peirce JL, Zhou M, Li H, Goldowitz D, Williams RW, Lu L, Cui Y: An integrative genomics strategy for systematic characterization of genetic loci modulating phenotypes. Hum Mol Genet. 2007, 16 (11): 1381-1390. 10.1093/hmg/ddm089.View ArticlePubMedGoogle Scholar
- Farber CR, Aten JE, Farber EA, de Vera V, Gularte R, Islas-Trejo A, Wen P, Horvath S, Lucero M, Lusis AJ, et al: Genetic dissection of a major mouse obesity QTL (Carfhg2): integration of gene expression and causality modeling. Physiol Genomics. 2009, 37 (3): 294-302. 10.1152/physiolgenomics.90245.2008.PubMed CentralView ArticlePubMedGoogle Scholar
- Holdt LM, Thiery J, Breslow JL, Teupser D: Increased ADAM17 mRNA expression and activity is associated with atherosclerosis resistance in LDL-receptor deficient mice. Arterioscler Thromb Vasc Biol. 2008, 28 (6): 1097-1103. 10.1161/ATVBAHA.108.165654.PubMed CentralView ArticlePubMedGoogle Scholar
- Papeta N, Chan KT, Prakash S, Martino J, Kiryluk K, Ballard D, Bruggeman LA, Frankel R, Zheng Z, Klotman PE, et al: Susceptibility loci for murine HIV-associated nephropathy encode trans-regulators of podocyte gene expression. J Clin Invest. 2009, 119 (5): 1178-1188. 10.1172/JCI37131.PubMed CentralView ArticlePubMedGoogle Scholar
- Doss S, Schadt EE, Drake TA, Lusis AJ: Cis-acting expression quantitative trait loci in mice. Genome Res. 2005, 15 (5): 681-691. 10.1101/gr.3216905.PubMed CentralView ArticlePubMedGoogle Scholar
- GuhaThakurta D, Xie T, Anand M, Edwards SW, Li G, Wang SS, Schadt EE: Cis-regulatory variations: a study of SNPs around genes showing cis-linkage in segregating mouse populations. BMC Genomics. 2006, 7: 235-10.1186/1471-2164-7-235.PubMed CentralView ArticlePubMedGoogle Scholar
- Yaguchi H, Togawa K, Moritani M, Itakura M: Identification of candidate genes in the type 2 diabetes modifier locus using expression QTL. Genomics. 2005, 85 (5): 591-599. 10.1016/j.ygeno.2005.01.006.View ArticlePubMedGoogle Scholar
- Yamashita S, Wakazono K, Nomoto T, Tsujino Y, Kuramoto T, Ushijima T: Expression quantitative trait loci analysis of 13 genes in the rat prostate. Genetics. 2005, 171 (3): 1231-1238. 10.1534/genetics.104.038174.PubMed CentralView ArticlePubMedGoogle Scholar
- Desert C, Duclos MJ, Blavy P, Lecerf F, Moreews F, Klopp C, Aubry M, Herault F, Le Roy P, Berri C, et al: Transcriptome profiling of the feeding-to-fasting transition in chicken liver. BMC Genomics. 2008, 9 (1): 611-10.1186/1471-2164-9-611.PubMed CentralView ArticlePubMedGoogle Scholar
- Leclecrq B, Blum JC, Boyer JP: Selecting broilers for low or high abdominal fat: initial observations. British Poultry Sciences. 1980, 107-113.Google Scholar
- Assaf S, Lagarrigue S, Daval S, Sansom M, Leclercq B, Michel J, Pitel F, Alizadeh M, Vignal A, Douaire M: Genetic linkage and expression analysis of SREBP and lipogenic genes in fat and lean chicken. Comp Biochem Physiol B Biochem Mol Biol. 2004, 137 (4): 433-441. 10.1016/j.cbpc.2004.02.005.View ArticlePubMedGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory cloning: a laboratory manual. New York. 1989Google Scholar
- Ianucelli N, Wolosyn N, Arhainx J, Gellin J, Milan D: GEMMA: A database to automate microsatellite genotyping. Animals Genetics. 1996, 27 (suppl 2): 55-Google Scholar
- Green P: Documentation for CRI-MAP version 2.4. Edited by: Falls K. 1990, CS. St Louis, MO: Washington University School of MedecineGoogle Scholar
- Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, Holloway E, Kolesnykov N, Lilja P, Lukk M: ArrayExpress--a public database of microarray experiments and gene expression profiles. Nucleic Acids Res. 2007, D747-750. 10.1093/nar/gkl995. 35 Database
- Casel P, Moreews F, Lagarrigue S, Klopp C: sigReannot: an oligo-set re-annotation pipeline based on similarities with the Ensembl transcripts and Unigene clusters. BMC Proc. 2009, 3 (Suppl 4): S3-10.1186/1753-6561-3-s4-s3.PubMed CentralView ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.PubMed CentralView ArticlePubMedGoogle Scholar
- Genepix User's guide and tutorial. [http://microarray.okstate.edu/files/GenePix_Pro.pdf/]
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30 (4): e15-10.1093/nar/30.4.e15.PubMed CentralView ArticlePubMedGoogle Scholar
- Elsen JM, Mangin B, Goffinet B, Boichard D, Le Roy P: Alternatives models for QTL detection in livestock.I.General introduction. Genetic Selection Evolution. 1999, 31: 213-224. 10.1186/1297-9686-31-3-213.View ArticleGoogle Scholar
- Le Roy P, Elsen JM, Boichard D, Mangin M, Bidanel JP, Goffinet B: An algorithm for QTL detection in mixture of full and half sib families. 6th World Congress of Genetic Applied to Livestock Production: 1998. 1998, University of Nex England, Armidale, 257-260.Google Scholar
- Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121 (1): 185-199.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.