Microsatellite-dense genetic map: towards genome coverage in a tropical maize (Zea mays L.) population

Dense molecular genetic maps are used for an efficient quantitative trait loci (QTL) mapping and in the marker-assisted selection programs. A dense genetic map was generated with 139 microsatellite markers using 256 F2 plants generated by the crossing of two tropical maize inbred lines (L-02-03D and L-20-01F). This map presented 1,858.61 cM in length, where 10 linkage groups were found spanned, with an average interval of 13.47 cM between adjacent markers. Seventy seven percent of the maize genetic mapping bins were covered, which means an increase of 14% coverage in relation to the previous tropical maize maps. The results provide a more detailed and informative genetic map in a tropical maize population representing the first step to make possible the studies of genetic architecture to identify and map QTL and estimate their effects on the variation of quantitative traits, thus allowing the manipulation and use in tropical maize breeding programs.

Évaluation neurobiologique des souris spontanément hypertendues : Du vieillissement à la génomique comparative

Le but de cette thèse est premièrement d’évaluer l’effet du vieillissement sur les fonctions psychomotrices des souches de souris sélectionnées génétiquement en fonction de leur tension artérielle (TA); deuxièmement, de localiser les déterminants génétiques des phénotypes psychophysiologiques à partir de souches recombinantes congéniques (RCS). Ces travaux ont mené à la publication de 4 articles. Le premier article décrit l’évaluation des fonctions psychomotrices des souches avec une tension artérielle élevée (HBP), basse (LBP) et normale (NBP). La performance aux épreuves d’exploration, d’habiletés motrices et d’apprentissage spatial, a été mesurée sur deux cohortes âgées respectivement de 12 mois et de trois mois. Indépendamment de l’âge, les HBPs sont hyperactives dans l’open-field (OF), mais pas dans le test d’exploration de trous. Inversement, les LBP explorent moins d’espaces que les NBP et, à trois mois seulement, sont hypoactives dans l’OF. Par ailleurs, les HBPs et les LBP présentent des déficits précoces de coordination motrice et des fonctions visuo-motrices. Le second article concerne l’évaluation longitudinale de la coordination motrice, de l’anxiété et de l’apprentissage spatial des souches HBP, LBP et NBP, à l’âge de deux mois et de 12 mois. Le vieillissement accentue l’hyperactivité des HBPs dans l’OF. Par contre, l’hypoactivité des souris LBP est détectable seulement à l’âge de deux mois. Indépendamment de l’âge, les souris HBP et LBP montrent une perception réduite du danger dans l’épreuve d’anxiété et des dysfonctions visuo-motrices au labyrinthe aquatique. Enfin, des déficits précoces de coordination motrice se manifestent seulement chez les HBPs. Il reste à déterminer si les déficits observés sont liés à des déterminants génétiques indépendants ou secondaires aux altérations de la tension artérielle. Le troisième article présente la comparaison entre les souches consanguines A/J et C57Bl/6J (B6) aux épreuves de l’OF, de la planche à trous, du labyrinthe aquatique et du cintre (coordination motrice). Les B6 explore d’avantage l’OF et la planche à trous. Les B6 sont moins rapides sur le cintre, mais supérieurs aux A/J dans le labyrinthe aquatique, avec une plate-forme invisible ou visible. Ces résultats démontrent l’implication de déterminants génétiques. Cette thèse se termine par un quatrième article sur la localisation des déterminants génétiques de la susceptibilité au stress dans les RCS, dérivées de A/J et B6, et présentant un agencement spécifique de 12.5% du génome. La réactivité émotionnelle est évaluée dans l’OF et le plus-maze; la réponse de stress est mesurée par radio télémétrie de la température interne pendant le stress d’immobilisation (SI) sous diète régulière et riche en sel; l’excrétion des électrolytes urinaires est dosée après 24 heures de diète salée. Les loci les plus significatifs sont situés dans les régions suivantes: de l’émotionalité dans l’OF (Emo1) sur le chr. 1 (LOD=4.6) correspondant à la région homologue impliquée dans la cohorte d’hypertension familiale du Saguenay; de la dopa décarboxylase (ddc) sur le chr. 11 pour l’émergence du plus-maze (LOD=4.7); de la protéine liant l’endotoxine (lbp) sur le chr. 2 pour l’hypothermie initiale en réponse au SI (LOD=4); et de HSP90 sur le chr. 12 pour l’excrétion de Ca++ (LOD=4.6). Des banques de données sont ensuite interrogées pour recenser les polymorphismes des régions régulatrices ou codantes des gènes candidats chez les souches ancestrales A/J et B6, dont les séquences sont disponibles pour le génome entier. Des utilitaires web permettent de dévoiler les changements dans la structure secondaire de l’ARNm, l’interférence avec des microARN ou avec d’autres motifs de liaison. Plusieurs SNPs fonctionnels ont été identifiés pour le QTL du chr. 1, particulièrement dans les éléments de régulation; ceux-ci impliquant des gènes reliés avec les réponses inflammatoire/immunitaire ou avec le système cardiovasculaire. La quantification par la PCR confirme une régulation à la baisse d’atp1a2 dans le cœur et le cerveau des souches susceptibles à l’anxiété. Ces résultats confirment l’intrication des altérations de la susceptibilité au stress et de la régulation de la TA.

Modularisation épistatique des loci à trait quantitatif associés à la pression artérielle et identification de gènes candidats pour l’hypertension

Problématique: L’hypertension artérielle essentielle, facteur de risque majeur dans le développement des maladies cardiovasculaires, est un trait multigénique complexe dont les connaissances sur le déterminisme génétique nécessitent d’être approfondies. De nombreux loci à trait quantitatif (QTLs); soit des gènes responsables de faire varier la pression artérielle (PA), ont été identifiés chez l’humain et le modèle animal. Cependant, le mystère plane encore sur la façon dont ces gènes fonctionnent ensemble pour réguler la PA. Hypothèse et objectif: Plutôt qu’une addition de QTLs ayant chacun une action infinitésimale sur la PA, une interaction épistatique entre les gènes serait responsable du phénotype hypertendu. Ainsi, l’étude de cette épistasie entre les gènes impliqués, directement ou indirectement, dans l’homéostasie de la PA nous permettrait d’explorer de nouvelles voies de régulation moléculaire en cause dans cette maladie. Méthodes: Via la réalisation de souches congéniques de rats, où un segment chromosomique provenant d’une souche receveuse hypertendue (Dahl Salt Sensitive, SS/Jr) est remplacé par son homologue provenant d’une souche donneuse normotendue (Lewis, LEW), des QTLs peuvent être mis en évidence. Dans ce contexte, la combinaison de QTLs via la création de doubles ou multiples congéniques constitue la première démonstration fonctionnelle des interactions intergéniques. Résultats: Vingt-sept combinaisons au total nous ont menés à l’appréciation d’une modularisation des QTLs. Ces derniers ont été catégorisés selon deux principaux modules épistatiques (EMs) où les QTLs appartenant à un même EM sont épistatiques entre eux et participent à une même voie régulatrice. Les EMs/cascades agissent alors en parallèle pour réguler la PA. Grâce à l’existence de QTLs ayant des effets opposés sur la PA, nous avons pu établir l’ordre hiérarchique entre trois paires de QTLs. Cependant, lorsque cette suite régulatrice ne peut être déterminée, d’autres approches sont nécessaires. Nos travaux nous ont mené à l’identification d’un QTL situé sur le chromosome 16 du rat (C16QTL), appartenant au EM1 et qui révélerait une nouvelle voie de l’homéostasie de la PA. Le gène retinoblastoma-associated protein 140 (Rap140)/family with sequence similarity 208 member A (Fam208a), présentant une mutation non synonyme entre SS/Jr et LEW est le gène candidat le plus plausible pour représenter C16QTL. Celui-ci code pour un facteur de transcription et semblerait influencer l’expression de Solute carrier family 7 (cationic amino acid transporter, y+ system) member 12 (Slc7a12), spécifiquement et significativement sous exprimé dans les reins de la souche congénique portant C16QTL par rapport à la souche SS/Jr. Rap140/Fam208a agirait comme un inhibiteur de la transcription de Slc7a12 menant à une diminution de la pression chez Lewis. Conclusions: L’architecture complexe de la régulation de la PA se dévoile mettant en scène de nouveaux acteurs, pour la plupart inconnus pour leur implication dans la PA. L’étude de la nouvelle voie de signalisation Rap140/Fam208a - Slc7a12 nous permettra d’approfondir nos connaissances quant à l’homéostasie de la pression artérielle et de l’hypertension chez SS/Jr. À long terme, de nouveaux traitements anti-hypertenseurs, ciblant plus d’une voie de régulation à la fois, pourraient voir le jour.

Domestication of plants in the americas: Insights from mendelian and molecular genetics

Background Plant domestication occurred independently in four different regions of the Americas. In general, different species were domesticated in each area, though a few species were domesticated independently in more than one area. The changes resulting from human selection conform to the familiar domestication syndrome, though different traits making up this syndrome, for example loss of dispersal, are achieved by different routes in crops belonging to different families. Genetic and Molecular Analyses of Domestication Understanding of the genetic control of elements of the domestication syndrome is improving as a result of the development of saturated linkage maps for major crops, identification and mapping of quantitative trait loci, cloning and sequencing of genes or parts of genes, and discoveries of widespread orthologies in genes and linkage groups within and between families. As the modes of action of the genes involved in domestication and the metabolic pathways leading to particular phenotypes become better understood, it should be possible to determine whether similar phenotypes have similar underlying genetic controls, or whether human selection in genetically related but independently domesticated taxa has fixed different mutants with similar phenotypic effects. Conclusions Such studies will permit more critical analysis of possible examples of multiple domestications and of the origin(s) and spread of distinctive variants within crops. They also offer the possibility of improving existing crops, not only major food staples but also minor crops that are potential export crops for developing countries or alternative crops for marginal areas.

Quantitative and qualitative differences in morphological traits revealed between diploid Fragaria species

Background and Aims: The aims of this investigation were to highlight the qualitative and quantitative diversity apparent between nine diploid Fragaria species and produce interspecific populations segregating for a large number of morphological characters suitable for quantitative trait loci analysis. Methods: A qualitative comparison of eight described diploid Fragaria species was performed and measurements were taken of 23 morphological traits from 19 accessions including eight described species and one previously undescribed species. A principal components analysis was performed on 14 mathematically unrelated traits from these accessions, which partitioned the species accessions into distinct morphological groups. Interspecific crosses were performed with accessions of species that displayed significant quantitative divergence and, from these, populations that should segregate for a range of quantitative traits were raised. Key Results: Significant differences between species were observed for all 23 morphological traits quantified and three distinct groups of species accessions were observed after the principal components analysis. Interspecific crosses were performed between these groups, and F2 and backcross populations were raised that should segregate for a range of morphological characters. In addition, the study highlighted a number of distinctive morphological characters in many of the species studied. Conclusions: Diploid Fragaria species are morphologically diverse, yet remain highly interfertile, making the group an ideal model for the study of the genetic basis of phenotypic differences between species through map-based investigation using quantitative trait loci. The segregating interspecific populations raised will be ideal for such investigations and could also provide insights into the nature and extent of genome evolution within this group.

Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea

Calcium (Ca) and magnesium (Mg) are the most abundant group II elements in both plants and animals. Genetic variation in shoot Ca and shoot Mg concentration (shoot Ca and Mg) in plants can be exploited to biofortify food crops and thereby increase dietary Ca and Mg intake for humans and livestock. We present a comprehensive analysis of within-species genetic variation for shoot Ca and Mg, demonstrating that shoot mineral concentration differs significantly between subtaxa (varietas). We established a structured diversity foundation set of 376 accessions to capture a high proportion of species-wide allelic diversity within domesticated Brassica oleracea, including representation of wild relatives (C genome, 1n = 9) from natural populations. These accessions and 74 modern F-1 hybrid cultivars were grown in glasshouse and field environments. Shoot Ca and Mg varied 2- and 2.3-fold, respectively, and was typically not inversely correlated with shoot biomass, within most subtaxa. The closely related capitata (cabbage) and sabauda (Savoy cabbage) subtaxa consistently had the highest mean shoot Ca and Mg. Shoot Ca and Mg in glasshouse-grown plants was highly correlated with data from the field. To understand and dissect the genetic basis of variation in shoot Ca and Mg, we studied homozygous lines from a segregating B. oleracea mapping population. Shoot Ca and Mg was highly heritable (up to 40). Quantitative trait loci (QTL) for shoot Ca and Mg were detected on chromosomes C2, C6, C7, C8, and, in particular, C9, where QTL accounted for 14 to 55 of the total genetic variance. The presence of QTL on C9 was substantiated by scoring recurrent backcross substitution lines, derived from the same parents. This also greatly increased the map resolution, with strong evidence that a 4-cM region on C9 influences shoot Ca. This region corresponds to a 0.41-Mb region on Arabidopsis (Arabidopsis thaliana) chromosome 5 that includes 106 genes. There is also evidence that pleiotropic loci on C8 and C9 affect shoot Ca and Mg. Map-based cloning of these loci will reveal how shoot-level phenotypes relate to Ca 21 and Mg 21 uptake and homeostasis at the molecular level.

Shoot zinc (Zn) concentration varies widely within Brassica oleracea L. and is affected by soil Zn and phosphorus (P) levels

The low availability of zinc (Zn) in soils and crops affects dietary Zn intake worldwide. This study sought to determine if the natural genetic variation in shoot Zn concentrations (Zn(shoot)) is sufficient to pursue a crop improvement breeding strategy in a leafy vegetable crop. The gene-pool of Brassica oleracea L. was sampled using a large (n = 376) diversity foundation set (DFS), representing almost all species-wide common allelic variation, and 74 commercial varieties (mostly F(1)). The DFS genotypes were grown at low and high soil phosphorus (P) levels under glasshouse and field conditions, and also in a Zn-deficient soil, with or without Zn-fertilisation, in a glasshouse. Despite the large variation in Zn(shoot) among genotypes, environment had a profound effect on Zn(shoot) The heritability of Zn(shoot) was significant, but relatively low, among 90 doubled-haploid (DH) lines from a mapping population. While several quantitative trait loci (QTL) associated with Zn(shoot) occurred on chromosomes C2, C3, C5, C7, and C9, these were generally weak and conditional upon growth conditions. Breeding for Zn(shoot) in B. oleracea is therefore likely to be challenging. Shoot P concentrations increased substantially in all genotypes under low soil Zn conditions. Conversely, only some genotypes had increased Zn(shoot) at low soil P levels. Sufficient natural genetic variation may therefore exist to study some of the interactions between Zn and P nutrition.

Shoot yield drives phosphorus use efficiency in Brassica oleracea and correlates with root architecture traits

The environmental and financial costs of using inorganic phosphate fertilizers to maintain crop yield and quality are high. Breeding crops that acquire and use phosphorus (P) more efficiently could reduce these costs. The variation in shoot P concentration (shoot-P) and various measures of P use efficiency (PUE) were quantified among 355 Brassica oleracea L. accessions, 74 current commercial cultivars, and 90 doubled haploid (DH) mapping lines from a reference genetic mapping population. Accessions were grown at two or more external P concentrations in glasshouse experiments; commercial and DH accessions were also grown in replicated field experiments. Within the substantial species-wide diversity observed for shoot-P and various measures of PUE in B. oleracea, current commercial cultivars have greater PUE than would be expected by chance. This may be a consequence of breeding for increased yield, which is a significant component of most measures of PUE, or early establishment. Root development and architecture correlate with PUE; in particular, lateral root number, length, and growth rate. Significant quantitative trait loci associated with shoot-P and PUE occur on chromosomes C3 and C7. These data provide information to initiate breeding programmes to improve PUE in B. oleracea.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast-and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

Natural genetic variation in caesium (Cs) accumulation by Arabidopsis thaliana

Ingestion of caesium (Cs) radioisotopes poses a health risk to humans. Crop varieties that accumulate less Cs in their edible tissues may provide a useful countermeasure. This study was performed to determine whether quantitative genetics on a model plant (Arabidopsis thaliana) might inform such 'safe'-crop strategies. Arabidopsis accessions and recombinant inbred lines (RILs), from Landsberg erecta (Ler) x Cape Verdi Island (Cvi), Ler x Columbia (Col), and Niederzenz (Nd) x Col mapping populations, were grown on agar supplemented with subtoxic levels of Cs. Shoot Cs concentration varied up to three-fold, and shoot f. wt varied up to 25-fold within populations. The heritability of growth and Cs accumulation traits ranged from 0.06 to 0.28. Four quantitative trait loci (QTL) accounted for > 80 of the genetic contribution to the total phenotypic variation in shoot Cs concentration in the Ler x Col population. QTL identified in this study, in particular, QTL co-localizing to the top and bottom regions of Chromosomes I and V in two different mapping populations, are amenable to positional cloning and, through collinearity, may inform selection or breeding strategies for the development of 'safe' crops.

High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities

Background and Aims: Phosphate (Pi) deficiency in soils is a major limiting factor for crop growth worldwide. Plant growth under low Pi conditions correlates with root architectural traits and it may therefore be possible to select these traits for crop improvement. The aim of this study was to characterize root architectural traits, and to test quantitative trait loci (QTL) associated with these traits, under low Pi (LP) and high Pi (HP) availability in Brassica napus. Methods: Root architectural traits were characterized in seedlings of a double haploid (DH) mapping population (n = 190) of B. napus 'Tapidor' x 'Ningyou 7' (TNDH) using high-throughput phenotyping methods. Primary root length (PRL), lateral root length (LRL), lateral root number (LRN), lateral root density (LRD) and biomass traits were measured 12 d post-germination in agar at LP and HP. Key Results: In general, root and biomass traits were highly correlated under LP and HP conditions. 'Ningyou 7' had greater LRL, LRN and LRD than 'Tapidor', at both LP and HP availability, but smaller PRL. A cluster of highly significant QTL for LRN, LRD and biomass traits at LP availability were identified on chromosome A03; QTL for PRL were identified on chromosomes A07 and C06. Conclusions: High-throughput phenotyping of Brassica can be used to identify root architectural traits which correlate with shoot biomass. It is feasible that these traits could be used in crop improvement strategies. The identification of QTL linked to root traits under LP and HP conditions provides further insights on the genetic basis of plant tolerance to P deficiency, and these QTL warrant further dissection.

Genetic analysis of potassium use efficiency in Brassica oleracea

Potassium (K) fertilizers are used in intensive and extensive agricultural systems to maximize production. However, there are both financial and environmental costs to K-fertilization. It is therefore important to optimize the efficiency with which K-fertilizers are used. Cultivating crops that acquire and/or utilize K more effectively can reduce the use of K-fertilizers. The aim of the present study was to determine the genetic factors affecting K utilization efficiency (KUtE), defined as the reciprocal of shoot K concentration (1/K(shoot)), and K acquisition efficiency (KUpE), defined as shoot K content, in Brassica oleracea. Genetic variation in K(shoot) was estimated using a structured diversity foundation set (DFS) of 376 accessions and in 74 commercial genotypes grown in glasshouse and field experiments that included phosphorus (P) supply as a treatment factor. Chromosomal quantitative trait loci (QTL) associated with K(shoot) and KUpE were identified using a genetic mapping population grown in the glasshouse and field. Putative QTL were tested using recurrent backcross substitution lines in the glasshouse. More than two-fold variation in K(shoot) was observed among DFS accessions grown in the glasshouse, a significant proportion of which could be attributed to genetic factors. Several QTL associated with K(shoot) were identified, which, despite a significant correlation in K(shoot) among genotypes grown in the glasshouse and field, differed between these two environments. A QTL associated with K(shoot) in glasshouse-grown plants (chromosome C7 at 62 center dot 2 cM) was confirmed using substitution lines. This QTL corresponds to a segment of arabidopsis chromosome 4 containing genes encoding the K(+) transporters AtKUP9, AtAKT2, AtKAT2 and AtTPK3. There is sufficient genetic variation in B. oleracea to breed for both KUtE and KUpE. However, as QTL associated with these traits differ between glasshouse and field environments, marker-assisted breeding programmes must consider carefully the conditions under which the crop will be grown.

Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

Key message We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species. Abstract Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F5 recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene—a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.

An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation, properties, and validation

MAGIC populations represent one of a new generation of crop genetic mapping resources combining high genetic recombination and diversity. We describe the creation and validation of an eight-parent MAGIC population consisting of 1091 F7 lines of winter-sown wheat (Triticum aestivum L.). Analyses based on genotypes from a 90,000-single nucleotide polymorphism (SNP) array find the population to be well-suited as a platform for fine-mapping quantitative trait loci (QTL) and gene isolation. Patterns of linkage disequilibrium (LD) show the population to be highly recombined; genetic marker diversity among the founders was 74% of that captured in a larger set of 64 wheat varieties, and 54% of SNPs segregating among the 64 lines also segregated among the eight founder lines. In contrast, a commonly used reference bi-parental population had only 54% of the diversity of the 64 varieties with 27% of SNPs segregating. We demonstrate the potential of this MAGIC resource by identifying a highly diagnostic marker for the morphological character "awn presence/absence" and independently validate it in an association-mapping panel. These analyses show this large, diverse, and highly recombined MAGIC population to be a powerful resource for the genetic dissection of target traits in wheat, and it is well-placed to efficiently exploit ongoing advances in phenomics and genomics. Genetic marker and trait data, together with instructions for access to seed, are available at http://www.niab.com/MAGIC/.

Detecting Major Genes Controlling Robustness of Chicken Body Weight Using Double Generalized Linear Models

Detecting both the majors genes that control the phenotypic mean and those controlling phenotypic variance has been raised in quantitative trait loci analysis. In order to mapping both kinds of genes, we applied the idea of the classic Haley-Knott regression to double generalized linear models. We performed both kinds of quantitative trait loci detection for a Red Jungle Fowl x White Leghorn F2 intercross using double generalized linear models. It is shown that double generalized linear model is a proper and efficient approach for localizing variance-controlling genes. We compared two models with or without fixed sex effect and prefer including the sex effect in order to reduce the residual variances. We found that different genes might take effect on the body weight at different time as the chicken grows.