Quantitative trait locus analysis using recombinant inbred intercrosses: theoretical and empirical considerations

We describe a new approach, called recombinant inbred intercross (RIX) mapping, that extends the power of recombinant inbred (RI) lines to provide sensitive detection of quantitative trait loci (QTL) responsible for complex genetic and nongenetic interactions. RIXs are generated by producing F1 hybrids between all or a subset of parental RI lines. By dramatically extending the number of unique, reproducible genomes, RIXs share some of the best properties of both the parental RI and F2 mapping panels. These attributes make the RIX method ideally suited for experiments requiring analysis of multiple parameters, under different environmental conditions and/or temporal sampling. However, since any pair of RIX genomes shares either one or no parental RIs, this cross introduces an unusual population structure requiring special computational approaches for analysis. Herein, we propose an efficient statistical procedure for QTL mapping with RIXs and describe a novel empirical permutation procedure to assess genome-wide significance. This procedure will also be applicable to diallel crosses. Extensive simulations using strain distribution patterns from CXB, AXB/BXA, and BXD mouse RI lines show the theoretical power of the RIX approach and the analysis of CXB RIXs demonstrates the limitations of this procedure when using small RI panels.     

Genetic characterization of a new set of recombinant inbred lines (LGXSM) formed from the intercross of SM/J and LG/J inbred mouse strains

A new set of LGXSM recombinant inbred (RI) strains is presented. The RI strain panel consists of 18 remaining strains of the original 55 founding strains. Strain characterization is based on 506 polymorphic microsatellites and 4289 single nucleotide polymorphisms (SNPs) distributed across the genome. Average microsatellite intermarker distance is 4.80 ± 4.84 Mb or 2.91 ± 3.21 F₂ cM. SNPs are more densely spaced at 0.57 ± 1.27 Mb. Ninety-five percent of all microsatellite intermarker intervals are separated by less than 15.00 Mb or 8.50 F₂ cM, while 95% of the SNPs are less than 0.95 Mb apart. Strains show expected low levels of nonsyntenic association among loci and complete genomic independence. During inbreeding, the RI strains went through strong natural selection on the agouti locus on Chromosome 2, especially when the epistatically interacting tyrosinase locus on Chromosome 7 carried the wild-type allele. The LG/J and SM/J strains differ in a large number of biomedically important traits, and they and their intercross progeny have been used in multiple mapping studies. The LG×SM RI strain panel provides a powerful new resource for mapping the genetic bases of complex traits and should prove to be of great biomedical utility in modeling complex human diseases such as obesity and diabetes. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users.     

Bayesian multiple quantitative trait loci mapping for recombinant inbred intercrosses

The Collaborative Cross (CC) is a renewable mouse resource that mimics the genetic diversity in humans. The recombinant inbred intercrosses (RIX) generated from CC recombinant inbred (RI) lines share similar genetic structures to those of F(2) individuals. In contrast to F(2) mice, genotypes of RIX can be inferred from the genotypes of their RI parents and can be produced repeatedly. Also, RIX mice do not typically share the same degree of relatedness. This unbalanced genetic relatedness requires careful statistical modeling to avoid a large number of false positive findings. For complex traits, mapping multiple genes simultaneously is arguably more powerful than mapping one gene at a time. In this article, we describe how we have developed a Bayesian quantitative trait locus (QTL) mapping method that simultaneously deals with the special genetic architecture of RIX and maps multiple genes. The performance of the proposed method is evaluated by extensive simulations. In addition, for a given set of RI lines, there are numerous ways to generate RIX samples. To provide a general guideline on future RIX studies, we compare several RIX designs through simulations.     

Polymorphisms revealed by PCR with single,short-sized,arbitrary primers are reliable markers for mouse and rat gene mapping

Ten single, arbitrarily designed oligodeoxynucleotide primers, with 50–70% (G+C) content, were used to amplify by polymerase chain reaction (PCR) sequences with DNA templates from several mouse species (Mus spretus, Mus musculus musculus, and Mus musculus domesticus), as well as DNA from the laboratory rat (Rattus norvegicus). Eight of these ten primers, used either individually or associated in pairs, generated a total of 13 polymorphic products which were used as genetic markers. All of these polymorphic sequences but one were mapped to a particular mouse chromosome, by use of DNA panels prepared either from interspecific backcross progeny of the type (C57BL/6 x Mus spretus)F₁ x C57BL/6 or DNA samples prepared from two sets of recombinant inbred (RI) strains (AKXL and BXD). Six rat-specific DNA segments were also assigned to a particular chromosome with DNA panels prepared from 18 rat/mouse somatic cell hybrids segregating rat chromosomes. From these experiments we conclude that, under precisely standardized PCR conditions, the DNA molecules amplified with these arbitrarily designed primers are useful and reliable markers for genetic mapping in both mouse and rat.     

The use of replicated progenies in marker-based mapping of QTL's

Summary Seven types of progeny are described which can be used in detection of linkage between marker loci and quantitative trait loci (QTL) in a cross between two inbred lines. Three types of progeny: recombinant inbred lines (RI); doubled haploid lines (DH); and S₁ lines can be used to detect linked main effects, d. DH and RI lines can be used to detect smaller effects than S₁ lines. However, S₁ lines can also be used to detect within-population dominance effects, h. The smallest d detectible is in the range of 1/2 to 1/12 the size of the corresponding LSD_(0.05) for the quantitative trait, using 100 lines and 6 replicates. The smallest h detectible is 3–4 times this size. Four types of progeny can be used to detect differences in the dominance behavior of alleles within the population relative to an allele in another inbred line (P₄: DH lines x P₄; RI lines x P₄; either F₂ x P₄ or S₁ lines x P₄; and progeny generated by crossing (F₁ x P₃) x P₄. Dominance differences in the range of 1 1/4 to 1/6 the size of the corresponding LSD_(0.05) are routinely detectible using 100 lines and 6 replicates. Increasing the numbers of progeny evaluated or the number of replicates allows for the detection of relatively smaller linked effects.Contribution of United AgriSeeds, Inc.     

Use of AFLP Markers for Gene Mapping and QTL Detection in the Rat

The AFLP technique is a new DNA marker technology based on the selective amplification of restriction fragments. Multiple polymorphic markers are simultaneously produced and can be tested in one PCR. No prior information on genomic DNA sequences is needed. In the current study, we contribute 18 AFLP markers to the linkage map of the rat. Seven AFLP markers were assigned to specific chromosomes by analysis of a (BN × ACI)F1 × ACI backcross progeny. Another 11 AFLP markers were mapped by using a panel of the H × B/B × H recombinant inbred (RI) strains. Genotypes of these AFLP markers were also tested for correlations with some blood pressure phenotypes in the RI strains. Suggestive correlation was found between the mean arterial pressure and two closely linked AFLP markers located on chromosome 20. The current study illustrates the value of AFLP markers for the construction of linkage maps and the detection of quantitative trait loci.     

Genetic loci for diet-induced atherosclerotic lesions and plasma lipids in mice

Genetic factors independent of those affecting plasma lipid levels are a major contributor to risk for atherosclerosis in humans, yet the basis for these is poorly understood. This study examined plasma lipids and diet-induced atherosclerosis in 16-month-old female mice of strains C56BL/6J and DBA/2J. Mice of the parental strains, from recombinant inbred strains derived from these (BXD RI), and F₂ progeny were fed an atherogenic diet for 16 weeks, beginning at 1 year of age. This induced atherosclerotic lesion formation in both parental strains, accompanied by increased plasma LDL levels. However, individual BXD RI strains and the BXD F₂ mice demonstrated a range of atherosclerotic lesion formation that was not or at best weakly correlated with plasma lipid levels. Quantitative trait locus (QTL) analysis of the BXD F₂ mice identified a locus with significant linkage (lod 4.5) for aortic lesion size on Chromosome (Chr) 10 that was independent of plasma lipids. Other loci with suggestive or significant linkage for various plasma lipid measures were identified on Chr 2, 3, 4, 5, 6, 7, 11, and 17. In this intercross, the genes primarily influencing atherosclerosis are distinct from those controlling plasma lipid levels.     

High genetic susceptibility to ethanol withdrawal predicts low ethanol consumption

C57BL/6J (B6) inbred mice are well known to drink large amounts of alcohol (ethanol) voluntarily and to have only modest ethanol-induced withdrawal under fixed dose conditions. In contrast, DBA/2J (D2) mice are ``teetotallers' and exhibit severe ethanol withdrawal. Speculation that an inverse genetic relationship existed between these two traits was substantiated by meta-analysis of existing data collected in multiple genetic models, including large panels of standard and recombinant inbred strains, their crosses, and selectively bred mouse lines. Despite methodological differences among laboratories in measurement of both preference drinking and withdrawal, a nearly universal finding was that genotypes consuming large amounts of 10% ethanol (calculated as g/kg/day) during two-bottle choice preference drinking were genetically predisposed to low withdrawal scores in independent studies after either acute or chronic ethanol treatment. Conversely, low-drinking genotypes had higher withdrawal severity scores. The genetic relationship appears to be strongest in populations derived from B6 and D2, where data from more genotypes (BXD RIs, B6D2F₂s, BXD RI F₁s, and B6D2F₂-derived selectively bred lines) were available for analysis. Gene mapping studies in these populations identified four chromosome regions [on Chromosomes (Chrs) 1, 2, 4, and 15] where genes might potentially influence both traits. Among genotypes with greater genetic diversity (for example, a panel of standard inbred strains or selectively bred lines), the relationship was less pronounced. Thus, reduced susceptibility to the development of high alcohol use may be supported by increased genetic susceptibility to ethanol withdrawal symptoms. Received: 15 September 1998 / Accepted: 8 October 1998     

Mapping of theLyb-4 gene to different chromosomes in DBA/2J and C3H/HeJ mice

Linkage has been established between the Lyb-4 alloantigen locus and the chromosome 4 markersLyb- 2 andMup- 1 using recombinant inbred (RI) strains. Only 2 of 24 BXD RI strains possess recombinant genotypes with respect to the B cell alloantigen lociLyb- 4 andLyb- 2, for an estimated recombination frequency of 0.024 ±0.019. One additional BXD RI strain was a recombinant with respect toLyb- 4 andMup- 1 (major urinary protein locus) for an estimated recombination frequency of 0.039 ± 0.026. These linkages were confirmed and further quantitated in a (C57BL/6J × DBA/2J)F₁ × C57BL/6J backcross population, in which the recombination frequency betweenLyb- 4 andMup- 1 was 0.049 ± 0.019. No recombination between the expression of Lyb-4.1 antigen and the ability of anti-Lyb-4.1 serum to suppress MLC reactivity was found, indicating that the genes controlling the antigenic determinant which is recognized with cytotoxic antibodies in anti-Lyb-4.1 serum is the same as, or is very closely linked to, the gene which is responsible for augmentation of the MLC response. In contrast, no linkage was observed between the gene controlling the Lyb-4.1 determinant andMup- 1 in RI strain and backcross mice derived from the cross of C3H/HeJ and C57BL/6J. Again, there was complete concordance between the serologically recognized determinant and the ability of anti-Lyb-4.1 serum to suppress the MLC response. Absorption of anti-Lyb-4.1 serum with C3H/HeJ, DBA/2J, and C57BL/6J lymphocytes, followed by the cytotoxic assay of the absorbed sera on lymphocytes of each of these three strains showed that serologically the Lyb-4.1 antigenic determinant on DBA/2 mice was indistinguishable from that on C3H/HeJ mice. Thus, both traits appear to be under the control of single genes in both DBA/2J and C3H/HeJ, but the C3H/HeJ gene appears to be nonallelic and unlinked to the DBA/2J gene.Abbreviations used in this paper LAD lymphocyte activating determinants - LPS lipopolysaccharide - MLC mixed lymphocyte culture - RI recombinant inbred     

Genetic structure of the LXS panel of recombinant inbred mouse strains: a powerful resource for complex trait analysis

The set of LXS recombinant inbred (RI) strains is a new and exceptionally large mapping panel that is suitable for the analysis of complex traits with comparatively high power. This panel consists of 77 strains—more than twice the size of other RI sets— and will typically provide sufficient statistical power (=0.8) to map quantitative trait loci (QTLs) that account for 25% of genetic variance with a genomewide p < 0.05. To characterize the genetic architecture of this new set of RI strains, we genotyped 330 MIT microsatellite markers distributed on all autosomes and the X Chromosome and assembled error-checked meiotic recombination maps that have an average F₂-adjusted marker spacing of 4 cM. The LXS panel has a genetic structure consistent with random segregation and subsequent fixation of alleles, the expected 3–4 × map expansion, a low level of nonsyntenic association among loci, and complete independence among all 77 strains. Although the parental inbred strains—Inbred Long-Sleep (ILS) and Inbred Short-Sleep (ISS)—were derived originally by selection from an 8-way heterogeneous stock selected for differential sensitivity to sedative effects of ethanol, the LXS panel is also segregating for many other traits. Thus, the LXS panel provides a powerful new resource for mapping complex traits across many systems and disciplines and should prove to be of great utility in modeling the genetics of complex diseases in human populations.(Robert W. Williams and Beth Bennett)These authors contributed equally to this work.     

A molecular marker that segregates with sorghum leaf blight resistance in one cross is maternally inherited in another

Leaf blight-resistant sorghum accession SC326-6 was crossed to the susceptible cultivar BTx623 to analyze the genetic basis for resistance. Field scoring of inoculated F₂ progeny revealed that resistance was transmitted as a dominant single-gene trait. By combining the random amplified polymorphic DNA (RAPD) technique with bulked-segregant analysis, it was possible to identify PCR amplification products that␣segregated with disease response. Primer OPD12 amplified a 323-bp band (D12R) that segregated with resistance. Creation of longer primers, or SCARs (sequence characterized amplified regions) for D12R resulted in the amplification of a single major band of the predicted size from all the resistant F₂ progeny and the resistant parent SC326-6, but not from BTx623 or 24 of 29 susceptible F₂ progeny. The SCAR primers also amplified a single band with DNA from IS3620C, the female parent in a cross with BTx623 that has been used to produce a recombinant inbred population for RFLP mapping. An equivalent band was amplified from all 137 recombinant inbred progeny, indicating that organelle DNA is the amplification target in this cross. Received: 31 July 1998 / Accepted: 23 November 1998     

Genotyping new BXD recombinant inbred mouse strains and comparison of BXD and consensus maps

Nine additional BXD recombinant inbred (RI) strains have been developed from the F₂ cross of C57BL/6J and DBA/2J mouse strains. A tenth line stopped breeding in the F₁₂ generation. F₂₀ generation breeding pairs from the nine surviving strains and an F₁₂ pair from the extinct line were genotyped at 319 genetic markers (primarily microsatellites) spanning most of the genome. Where typing data were lacking, the established set of 26 BXD strains also were genotyped at these same loci. The availability of these additional nine strains enhances the value of the BXD RI set for analysis of complex phenotypic traits. The proportion of loci still segregating at the F₂₀ generation was found to closely approximate expectation, suggesting that selection favoring the retention of heterozygosity is not a strong factor. However, the number of crossovers between adjacent markers was frequently less than predicted from consensus map distances. A significant deficiency of recombinants was observed on Chrs 3, 4, 14, and X. On Chr 14, the estimated cumulative BXD map distance between the most proximal and distal markers was only 30.2 cM, compared with a distance of 60.0 cM in the consensus map. On the X Chr, the estimated and predicted cumulative distances were 38.8 and 69.5 cM, respectively. Over all chromosomes, the BXD RI map is 14.5% shorter than predicted from the consensus map. It is suggested that distances in some of the consensus maps are inflated. Alternatively, recombinant genotypes could be selected against during inbreeding owing to allelic interactions affecting fitness. The latter interpretation implies that relatively strong intrachromosomal epistasis is common. Received: 2 October 1998 / Accepted: 15 December 1998     

Simulating the collaborative cross: power of quantitative trait loci detection and mapping resolution in large sets of recombinant inbred strains of mice

It has been suggested that the collaborative cross, a large set of recombinant inbred strains derived from eight inbred mouse strains, would be a powerful resource for the dissection of complex phenotypes. Here we use simulation to investigate the power of the collaborative cross to detect and map small genetic effects. We show that for a fixed population of 1000 individuals, 500 RI lines bred using a modified version of the collaborative cross design are adequate to map a single additive locus that accounts for 5% of the phenotypic variation to within 0.96 cM. In the presence of strong epistasis more strains can improve detection, but 500 lines still provide sufficient resolution to meet most goals of the collaborative cross. However, even with a very large panel of RILs, mapping resolution may not be sufficient to identify single genes unambiguously. Our results are generally applicable to the design of RILs in other species.     

The locus Om,responsible for the DDK syndrome,maps close to Sigje on mouse Chromosome 11

The DDK inbred strain of mouse has a striking particularity: when DDK females are crossed to males of other strains they exhibit a reduced fertility, whereas the reciprocal crosses (non-DDK females x DDK males) are fertile (Wakasugi et al. 1967; Wakasugi 1973). The low fertility results from an early embryonic lethality, the F₁ embryos dying near the late morula-early blastocyst stage. Genetic analyses (Wakasugi 1974) and nuclear and cytoplasmic transfers (Renard and Babinet 1986; Babinet et al. 1990; Mann 1986), have shown that the failure of the embryos to develop is due to an incompatibility between a DDK maternally encoded cytoplasmic product and the non-DDK paternal genome. In order to elucidate the genetic determinism of this embryonic lethality, we have analyzed the fertility of male progeny from a backcross BALB/c females x (BALB/c x DDK)F₁ males and that of males from a set of recombinant inbred (RI) strains, established from DDK and BALB/c progenitors, when mated with DDK females. Our results indicate that a single locus, Om, is responsible for the DDK syndrome and is located on Chromosome (Chr) 11, very close to the Sigje locus.     

Blood pressure in 15 inbred mouse strains and its lack of relation with obesity and insulin resistance in the progeny of an NZO/HILtJ × C3H/HeJ intercross

We characterized the systolic and diastolic blood pressures of 10-week-old males from 15 inbred mouse strains and found that blood pressures among strains were continuously distributed and that strain C3H/HeJ had the lowest mean systolic and diastolic pressure (100.5 ± 3.2 and 66.8 ± 3.5 mmHg), and a strain with obesity and diabetes, NZO/HILtJ, had the highest (132.4 ± 3.1 and 86.6 ± 6.9 mmHg). To understand the relationship of blood pressure with insulin resistance and obesity, we produced F₁ and F₂ progeny from reciprocal crosses of NZO, the strain with obesity, diabetes, and high blood pressure, and the strain with the lowest blood pressures, C3H/HeJ. Mean systolic pressures of 10-week-old (NZO × C3H)F₁ and (C3H × NZO)F₁ males were similar to each other (114.9 ± 3.8 and 117.2 ± 5.0 mmHg) and were intermediate to those of the parental strains. Systolic pressure of F₂ males (n = 223) was distributed normally about the mean, suggesting that blood pressure is a polygenic trait. The body mass index (BMI) and plasma insulin levels of F₂ progeny correlated significantly and positively with plasma leptin levels, suggesting that obesity is associated with insulin resistance. In contrast, systolic pressure did not correlate with BMI, plasma leptin levels, and plasma insulin levels, suggesting that genes underlying the development of hypertension in this intercross are not associated with the development of obesity and insulin resistance. Our results demonstrate that the progeny of NZO and C3H intercrosses are a practical and powerful tool for identifying blood pressure genes and for understanding human polygenic hypertension.(Fumihiro Sugiyama) These authors contributed equally to this study.     

Genetics and ontogeny of alcohol dehydrogenase isozymes in the mouse: Evidence for a cis-acting regulator gene (Adt-1) controlling C2 isozyme expression in reproductive tissues and close linkage of Adh-3 and Adt-1 on chromsome 3

An electrophoretic variant previously reported for the stomach isozyme of alcohol dehydrogenase (ADH-C2) in inbred strains of Mus musculus (Holmes, 1977) has been used to localize the gene encoding this enzyme (Adh-3) on chromosome 3 near Va (varitint) (9.6 ± 3.6% recombinants). Genetic variation of ADH-C2 activity in male and female reproductive tissues among inbred strains and Harwell linkage testing stocks was also observed. Reproductive tissue ADH-C2 phenotypes were inherited in a normal Mendelian fashion among F₂ progeny of an F₁ (LII × C57BL/Go) × C57BL/Go backcross as though controlled by a single cis-acting regulator locus (designated Adt-1) with two alleles: Adt-1 ^a (presence of ADH-C2) and Adt-1 ^b (absence or low activity of ADH-C2). No recombinants were observed among 73 progeny or among 13 inbred strains and six Harwell linkage testing stocks of mice, indicating that Adh-3 and Adt-1 are closely linked or identical genes. A single recombinant phenotype was observed in Peru-Coppock mice, suggesting that they are separate genes. Ontogenetic analyses demonstrated that ADH-B₂ is present throughout development from late fetal stages in stomach, liver, and kidney; similar results were found for ADH-C2 in developing kidney and stomach extracts, whereas ADH-A2 exhibited high activity in liver extracts after 3 weeks of age in both sexes and in male kidney extracts after 6 weeks.     

Basal serum immunoglobulin levels

Two congenic strains of mice were identified that differ in their serum immunoglobulin levels. The strains were crossed, the F₁ progeny were intercrossed, and the serum immunoglobulin levels of the F₁ and F₂ progeny were analyzed. The F₁ mice have serum immunoglobulin levels like that of the high parent, and the low-immunoglobulin phenotype segregates in the F₂ population. Six other inbred strains of mice were also characterized for basal serum levels of five classes of immunoglobulin.     

Mapping quantitative trait loci using molecular marker linkage maps

Summary High-density restriction fragment length polymorphism (RFLP) and allozyme linkage maps have been developed in several plant species. These maps make it technically feasible to map quantitative trait loci (QTL) using methods based on flanking marker genetic models. In this paper, we describe flanking marker models for doubled haploid (DH), recombinant inbred (RI), backcross (BC), F₁ testcross (F₁TC), DH testcross (DHTC), recombinant inbred testcross (RITC), F₂, and F₃ progeny. These models are functions of the means of quantitative trait locus genotypes and recombination frequencies between marker and quantitative trait loci. In addition to the genetic models, we describe maximum likelihood methods for estimating these parameters using linear, nonlinear, and univariate or multivariate normal distribution mixture models. We defined recombination frequency estimators for backcross and F₂ progeny group genetic models using the parameters of linear models. In addition, we found a genetically unbiased estimator of the QTL heterozygote mean using a linear function of marker means. In nonlinear models, recombination frequencies are estimated less efficiently than the means of quantitative trait locus genotypes. Recombination frequency estimation efficiency decreases as the distance between markers decreases, because the number of progeny in recombinant marker classes decreases. Mean estimation efficiency is nearly equal for these methods.     

Fine mapping of a major quantitative trait loci, qSSP7, controlling the number of spikelets per panicle as a single Mendelian factor in rice

In our previous studies, one putative QTL affecting number of spikelets per panicle (SPP) was identified in the pericentromeric region of rice chromosome 7 using a recombinant inbred population. In order to define the QTL (qSPP7), RI50, a recombinant inbred line with 70% of genetic background same as the female parent of Zhenshan 97, was selected to produce near-isogenic lines for the target region in the present study. In a BC₂F₂ population consisting of 190 plants, the frequency distribution of SPP was shown to be discontinuous and followed the expected Mendelian ratios (1:2:1 by progeny test) for single locus segregation. qSPP7 was mapped to a 0.4 cM region between SSR marker RM3859 and RFLP marker C39 based on tests of the BC₂F₂ population and its progeny. Its additive and dominant effects on SPP were 51.1 and 24.9 spikelets, respectively. Of great interest, the QTL region also had effects on grain yield per plant (YD), 1,000 grain weight (GW), tillers per plant (TPP) and seed setting ratio (SR). Significant correlations were observed between SPP and YD (r = 0.66) and between SPP and SR (r = −0.29) in the progeny test. 1082 extremely small panicle plants of a BC₃F₂ population containing 8,400 individuals were further used to fine map the QTL. It turns out that qSPP7 co-segregated with two markers, RM5436 and RM5499 spanning a physical distance of 912.4 kb. Overall results suggested that recombination suppression occurred in the region and positional cloning strategy is infeasible for qSPP7 isolation. The higher grain yield of Minghui 63 homozygote as compared to the heterozygote suggested that Minghui 63 homozygote at qSPP7 in hybrid rice could further improve its yield. Y. Z. Xing and W. J. Tang contributed equally to this work.     

Varying coefficient models for mapping quantitative trait loci using recombinant inbred intercrosses

There has been a great deal of interest in the development of methodologies to map quantitative trait loci (QTL) using experimental crosses in the last 2 decades. Experimental crosses in animal and plant sciences provide important data sources for mapping QTL through linkage analysis. The Collaborative Cross (CC) is a renewable mouse resource that is generated from eight genetically diverse founder strains to mimic the genetic diversity in humans. The recombinant inbred intercrosses (RIX) generated from CC recombinant inbred (RI) lines share similar genetic structures of F(2) individuals but with up to eight alleles segregating at any one locus. In contrast to F(2) mice, genotypes of RIX can be inferred from the genotypes of their RI parents and can be produced repeatedly. Also, RIX mice typically do not share the same degree of relatedness. This unbalanced genetic relatedness requires careful statistical modeling to avoid false-positive findings. Many quantitative traits are inherently complex with genetic effects varying with other covariates, such as age. For such complex traits, if phenotype data can be collected over a wide range of ages across study subjects, their dynamic genetic patterns can be investigated. Parametric functions, such as sigmoidal or logistic functions, have been used for such purpose. In this article, we propose a flexible nonparametric time-varying coefficient QTL mapping method for RIX data. Our method allows the QTL effects to evolve with time and naturally extends classical parametric QTL mapping methods. We model the varying genetic effects nonparametrically with the B-spline bases. Our model investigates gene-by-time interactions for RIX data in a very flexible nonparametric fashion. Simulation results indicate that the varying coefficient QTL mapping has higher power and mapping precision compared to parametric models when the assumption of constant genetic effects fails. We also apply a modified permutation procedure to control overall significance level.