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.     

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.     

Quantitative trait mapping in a diallel cross of recombinant inbred lines

A recombinant inbred intercross (RIX) is created by generating diallel F₁ progeny from one or more panels of recombinant inbred (RI) strains. This design was originally introduced to extend the power of small RI panels for the confirmation of quantitative trait loci (QTL) provisionally detected in a parental RI set. For example, the set of 13 C × B (C57BL/6ByJ × BALB/cByJ) RI strains can, in principle, be supplemented with 156 isogenic F₁s. We describe and test a method of analysis, based on a linear mixed model, that accounts for the correlation structure of RIX populations. This model suggests a novel permutation algorithm that is needed to obtain appropriate threshold values for genome-wide scans of an RIX population. Despite the combinational multiplication of unique genotypes that can be generated using an RIX design, the effective sample size of the RIX population is limited by the number of progenitor RI genomes that are combined. When using small RI panels such as the C × B there appears to be only modest advantage of the RIX design when compared with the original RI panel for detecting QTLs with additive effects. The RIX, however, does have an inherent ability to detect dominance effects, and, unlike RI strains, the RIX progeny are genetically reproducible but are not fully inbred, providing somewhat more natural genetic context. We suggest a breeding strategy, the balanced partial RIX, that balances the advantage of RI and RIX designs. This involves the use of a partial RIX population derived from a large RI panel in which the available information is maximized by minimizing correlations among RIX progeny.     

Heritability and coefficient of genetic variation analyses of phenotypic traits provide strong basis for high-resolution QTL mapping in the Collaborative Cross mouse genetic reference population

Most biological traits of human importance are complex in nature; their manifestation controlled by the cumulative effect of many genetic factors interacting with one another and with the individual’s life history. Because of this, mouse genetic reference populations (GRPs) consisting of collections of inbred lines or recombinant inbred lines (RIL) derived from crosses between inbred lines are of particular value in analysis of complex traits, since massive amounts of data can be accumulated on the individual lines. However, existing mouse GRPs are derived from inbred lines that share a common history, resulting in limited genetic diversity, and reduced mapping precision due to long-range gametic disequilibrium. To overcome these limitations, the Collaborative Cross (CC) a genetically highly diverse collection of mouse RIL was established. The CC, now in advanced stages of development, will eventually consist of about 500 RIL derived from reciprocal crosses of eight divergent founder strains of mice, including three wild subspecies. Previous studies have shown that the CC indeed contains enormous diversity at the DNA level, that founder haplotypes are inherited in expected frequency, and that long-range gametic disequilibrium is not present. We here present data, primarily from our own laboratory, documenting extensive genetic variation among CC lines as expressed in broad-sense heritability (H²) and by the well-known “coefficient of genetic variation,” demonstrating the ability of the CC resource to provide unprecedented mapping precision leading to identification of strong candidate genes.     

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     

Development of pre-isogenic lines for rice blast-resistance by marker-aided selection from a recombinant inbred population

To increase the available set of near-isogenic lines (NILs) for blast-resistance in rice, we have developed a general method for establishing NILs from populations of fixed recombinants that have been used for gene mapping. We demonstrated the application of this method by the selection of lines carrying genes from the rice cultivar Moroberekan. Moroberekan is a West African japonica cultivar that is considered to have durable resistance to rice blast. Multiple genes from Moroberekan conferring complete and partial resistance to blast have previously been mapped using a recombinant inbred (RI) population derived from a cross between Moroberekan and the highly and broadly susceptible indica cultivar CO39. To analyze individual blast-resistance genes, it is desirable to transfer them individually into a susceptible genetic background. This RI population, and the associated data sets on blast reaction and restriction fragment length polymorphism (RFLP) genotypes, were used for selection of lines likely to carry individual blast-resistance genes and a minimum number of chromosomal segments from Moroberekan. Because skewed segregation in the RI population favored CO39 (indica) alleles, resistant lines carrying 8.7–17.5% of Moroberekan alleles (the proportion expected after two or three backcrosses) could be selected. We chose three RI lines carrying different complete resistance genes to blast and two RI lines carrying partial resistance genes to blast as potential parents for the development of NILs. These lines were subjected to genetic analysis, which allowed clarification of some issues that could not be resolved during the initial gene-mapping study.     

Rapid Identification of Major-Effect Genes Using the Collaborative Cross

The Collaborative Cross (CC) was designed to facilitate rapid gene mapping and consists of hundreds of recombinant inbred lines descended from eight diverse inbred founder strains. A decade in production, it can now be applied to mapping projects. Here, we provide a proof of principle for rapid identification of major-effect genes using the CC. To do so, we chose coat color traits since the location and identity of many relevant genes are known. We ascertained in 110 CC lines six different coat phenotypes: albino, agouti, black, cinnamon, and chocolate coat colors and the white-belly trait. We developed a pipeline employing modifications of existing mapping tools suitable for analyzing the complex genetic architecture of the CC. Together with analysis of the founders’ genome sequences, mapping was successfully achieved with sufficient resolution to identify the causative genes for five traits. Anticipating the application of the CC to complex traits, we also developed strategies to detect interacting genes, testing joint effects of three loci. Our results illustrate the power of the CC and provide confidence that this resource can be applied to complex traits for detection of both qualitative and quantitative trait loci.     

Chromosomal locations and gonadal dependence of genes that mediate resistance to ectromelia (mousepox) virus-induced mortality.

Four genetic loci were tested for linkage with loci that control genetic resistance to lethal ectromelia virus infection in mice. Three of the loci were selected because of concordance with genotypes assigned to recombinant inbred (RI) strains of mice derived from resistant C57BL/6 and susceptible DBA/2 (BXD) mice on the basis of their responses to challenge infection. Thirty-six of 167 male (C57BL/6 x DBA/2)F1 x DBA/2 backcross (BC) mice died (22%), of which 27 (75%) were homozygous for DBA/2 alleles at Hc and H-2D. Twenty-eight percent of sham-castrated and 6% of sham-ovariectomized BC mice were susceptible to lethal mousepox, whereas 50% of gonadectomized mice were susceptible. There was no linkage evident between Hc or H-2D and loci that controlled resistance to lethal ectromelia virus infection in 44 castrated BC mice. Mortality among female mice of BXD RI strains with susceptible or intermediate male phenotypes was strongly correlated (r = 0.834) with male mortality. Gonadectomized C57BL/6 mice were as resistant as intact mice to lethal ectromelia virus infection. These results indicate that two gonad-dependent genes on chromosomes 2 and 17 and one gonad-independent gene control resistance to mousepox virus infection, that males and females share gonad-dependent genes, and that the gonad-independent gene is fully protective.     

The genome architecture of the Collaborative Cross mouse genetic reference population

The Collaborative Cross Consortium reports here on the development of a unique genetic resource population. The Collaborative Cross (CC) is a multiparental recombinant inbred panel derived from eight laboratory mouse inbred strains. Breeding of the CC lines was initiated at multiple international sites using mice from The Jackson Laboratory. Currently, this innovative project is breeding independent CC lines at the University of North Carolina (UNC), at Tel Aviv University (TAU), and at Geniad in Western Australia (GND). These institutions aim to make publicly available the completed CC lines and their genotypes and sequence information. We genotyped, and report here, results from 458 extant lines from UNC, TAU, and GND using a custom genotyping array with 7500 SNPs designed to be maximally informative in the CC and used a novel algorithm to infer inherited haplotypes directly from hybridization intensity patterns. We identified lines with breeding errors and cousin lines generated by splitting incipient lines into two or more cousin lines at early generations of inbreeding. We then characterized the genome architecture of 350 genetically independent CC lines. Results showed that founder haplotypes are inherited at the expected frequency, although we also consistently observed highly significant transmission ratio distortion at specific loci across all three populations. On chromosome 2, there is significant overrepresentation of WSB/EiJ alleles, and on chromosome X, there is a large deficit of CC lines with CAST/EiJ alleles. Linkage disequilibrium decays as expected and we saw no evidence of gametic disequilibrium in the CC population as a whole or in random subsets of the population. Gametic equilibrium in the CC population is in marked contrast to the gametic disequilibrium present in a large panel of classical inbred strains. Finally, we discuss access to the CC population and to the associated raw data describing the genetic structure of individual lines. Integration of rich phenotypic and genomic data over time and across a wide variety of fields will be vital to delivering on one of the key attributes of the CC, a common genetic reference platform for identifying causative variants and genetic networks determining traits in mammals.     

Identification of Hepatocarcinogen-Resistance Genes in Dba/2 Mice

Male DBA/2J mice are ~20-fold more susceptible than male C57BL/6J mice to hepatocarcinogenesis induced by perinatal treatment with N,N-diethylnitrosamine (DEN). In order to elucidate the genetic control of hepatocarcinogenesis in DBA/2J mice, male BXD recombinant inbred, D2B6F(1) X B6 backcross, and D2B6F(2) intercross mice were treated at 12 days of age with DEN and liver tumors were enumerated at 32 weeks. Interestingly, the distribution of mean tumor multiplicities among BXD recombinant inbred strains indicated that hepatocarcinogen-sensitive DBA/2 mice carry multiple genes with opposing effects on the susceptibility to liver tumor induction. By analyzing D2B6F(1) X B6 backcross and D2B6F(2) intercross mice for their liver tumor multiplicity phenotypes and for their genotypes at simple sequence repeat marker loci, we mapped two resistance genes carried by DBA/2J mice, designated Hcr1 and -2, to chromosomes 4 and 10, respectively. Hcr1 and Hcr2 resolved the genetic variance in the backcross population well, indicating that these resistance loci are the major determinants of the variance in the backcross population. Although our collection of 100 simple sequence repeat markers allowed linkage analysis for ~95% of the genome, we failed to map any sensitivity alleles for DBA/2J mice. Thus, it is likely that the susceptibility of DBA/2J mice is the consequence of the combined effects of multiple sensitivity loci.     

The Power of QTL Mapping with RILs

QTL (quantitative trait loci) mapping is commonly used to identify genetic regions responsible to important phenotype variation. A common strategy of QTL mapping is to use recombinant inbred lines (RILs), which are usually established by several generations of inbreeding of an F1 population (usually up to F6 or F7 populations). As this inbreeding process involves a large amount of labor, we are particularly interested in the effect of the number of inbreeding generations on the power of QTL mapping; a part of the labor could be saved if a smaller number of inbreeding provides sufficient power. By using simulations, we investigated the performance of QTL mapping with recombinant inbred lines (RILs). As expected, we found that the power of F4 population could be almost comparable to that of F6 and F7 populations. A potential problem in using F4 population is that a large proportion of RILs are heterozygotes. We here introduced a new method to partly relax this problem. The performance of this method was verified by simulations with a wide range of parameters including the size of the segregation population, recombination rate, genome size and the density of markers. We found our method works better than the commonly used standard method especially when there are a number of heterozygous markers. Our results imply that in most cases, QTL mapping does not necessarily require RILs at F6 or F7 generations; rather, F4 (or even F3) populations would be almost as useful as F6 or F7 populations. Because the cost to establish a number of RILs for many generations is enormous, this finding will cause a reduction in the cost of QTL mapping, thereby accelerating gene mapping in many species.     

Quantitative trait loci for thermotolerance phenotypes in Drosophila melanogaster

For insects, temperature is a major environmental variable that can influence an individual's behavioral activities and fitness. Drosophila melanogaster is a cosmopolitan species that has had great success in adapting to and colonizing diverse thermal niches. This adaptation and colonization has resulted in complex patterns of genetic variation in thermotolerance phenotypes in nature. Although extensive work has been conducted documenting patterns of genetic variation, substantially less is known about the genomic regions or genes that underlie this ecologically and evolutionarily important genetic variation. To begin to understand and identify the genes controlling thermotolerance phenotypes, we have used a mapping population of recombinant inbred (RI) lines to map quantitative trait loci (QTL) that affect variation in both heat- and cold-stress resistance. The mapping population was derived from a cross between two lines of D. melanogaster (Oregon-R and 2b) that were not selected for thermotolerance phenotypes, but exhibit significant genetic divergence for both phenotypes. Using a design in which each RI line was backcrossed to both parental lines, we mapped seven QTL affecting thermotolerance on the second and third chromosomes. Three of the QTL influence cold-stress resistance and four affect heat-stress resistance. Most of the QTL were trait or sex specific, suggesting that overlapping but generally unique genetic architectures underlie resistance to low- and high-temperature extremes. Each QTL explained between 5 and 14% of the genetic variance among lines, and degrees of dominance ranged from completely additive to partial dominance. Potential thermotolerance candidate loci contained within our QTL regions are identified and discussed.     

Genetic Variation in the Shape of the Mouse Mandible and Its Relationship to Glucocorticoid-Induced Cleft Palate Analyzed by Using Recombinant Inbred Lines

Variation in mandible shape has been investigated in a set of recombinant inbred (RI) lines of mice, the C57BL/6J X A/J (BXA;AXB) RI lines. Considerable genetic variation was detected between the RI lines, but most lines were intermediate in shape when compared with the parent lines. Variation in mandible shape could not be explained by any single gene differences known between the parent lines including the H-2 locus. Some RI lines had mandible shapes unlike either parent, and one in particular, line BXA1, had an unusual shape with a pronounced condyloid process. It was concluded that mandible shape has a complex inheritance involving a number of genes, each with small effects. In some cases, recombination of the genes can produce bone shapes quite different from those of the original parent line.--There was no evidence that the variability in steroid-induced cleft palate incidence in the BXA;AXB RI lines is related to the variation in adult mandible shape as detected in this study.     

Gene and QTL detection in a three-way barley cross under selection by a mixed model with kinship information using SNPs

Quantitative trait locus (QTL) detection is commonly performed by analysis of designed segregating populations derived from two inbred parental lines, where absence of selection, mutation and genetic drift is assumed. Even for designed populations, selection cannot always be avoided, with as consequence varying correlation between genotypes instead of uniform correlation. Akin to linkage disequilibrium mapping, ignoring this type of genetic relatedness will increase the rate of false-positives. In this paper, we advocate using mixed models including genetic relatedness, or ‘kinship’ information for QTL detection in populations where selection forces operated. We demonstrate our case with a three-way barley cross, designed to segregate for dwarfing, vernalization and spike morphology genes, in which selection occurred. The population of 161 inbred lines was screened with 1,536 single nucleotide polymorphisms (SNPs), and used for gene and QTL detection. The coefficient of coancestry matrix was estimated based on the SNPs and imposed to structure the distribution of random genotypic effects. The model incorporating kinship, coancestry, information was consistently superior to the one without kinship (according to the Akaike information criterion). We show, for three traits, that ignoring the coancestry information results in an unrealistically high number of marker–trait associations, without providing clear conclusions about QTL locations. We used a number of widely recognized dwarfing and vernalization genes known to segregate in the studied population as landmarks or references to assess the agreement of the mapping results with a priori candidate gene expectations. Additional QTLs to the major genes were detected for all traits as well.     

Genetic analysis of susceptibility to isoniazid-induced seizures in mice

Four sets of recombinant inbred lines of mice have been used to analyze genetic differences in acute toxicity of the drug, isonicotinic acid hydrazide. Standard inbred strains, their F₁ hybrids and recombinant inbred strains were all challenged with a single dose of the drug. The percent mortality of the different groups was analyzed to estimate heritability and the number of genes affecting resistance. The data indicated that resistance factors were dominant, heritability was moderate (.25-.37), and more than one gene was involved in each of four different sets of recombinant inbred lines. Possible approaches for identifying and mapping individual genes affecting resistance are discussed.     

Genetic regulation of gene-specific mRNA by ethanolin vivo and its possible role in ethanol preference in a cross with RI lines in mice

Relative ethanol preference is a well-recognized phenotype in a number of species, including mice, but the molecular basis for this phenotype remains speculative. We generated novel recombinant inbred (RI) mouse lines from C57BL/6J (ethanol preferring) and BALB/c (ethanol avoiding) strains and evaluated the effect of ethanol feeding on the mRNA levels of three genes (Adh-1, Adh-2, andCas-1) of alcohol metabolism. Ethanol feeding affects the mRNA levels of all three genes in both a gene- and a genotype-specific manner. The effect of ethanol feeding onAhd-2 mRNA, in particular, is highly correlated with the relative ethanol acceptance of the genotypes. DNA sequencing of ∼500 bp of the 5′ upstream region of theAhd-2 gene has yielded identical sequence for the two strains and the genetically determined associated factors are hypothesized to be regulatory proteins. Quantitative trait locus analysis on the RI lines should lead to the molecular characterization and mapping of such gene-specific regulatory factors.     

Integration of simple sequence repeat (SSR) markers into a molecular linkage map of common bean (Phaseolus vulgaris L.)

Microsatellite or simple sequence repeat (SSR) markers have been successfully used for genomic mapping, DNA fingerprinting, and marker-assisted selection in many plant species. Here we report the first successful assignment of 15 SSR markers to the Phaseolus vulgaris molecular linkage map. A total of 37 SSR primer pairs were developed and tested for amplification and product-length polymorphism with BAT93 and Jalo EEP558, the parental lines of an F7 recombinant inbred (RI) population previously used for the construction of a common bean molecular linkage map. Sixteen of the SSRs polymorphic to the parental lines were analyzed for segregation and 15 of them were assigned to seven different linkage groups, indicating a widespread distribution throughout the bean genome. Map positions for genes coding for DNAJ-like protein, pathogenesis-related protein 3, plastid-located glutamine synthetase, endochitinase, sn-glycerol-3 phosphate acyltransferase, NADP-dependent malic enzyme, and protein kinase were determined for the first time. Addition of three SSR loci to linkage group B4 brought two separated smaller linkage groups together to form a larger linkage group. Analysis of allele segregation in the F7 RI population revealed that all 16 SSRs segregated in the expected 1:1 ratio. These SSR markers were stable and easy to assay by polymerase chain reaction (PCR). They should be useful markers for genetic mapping, genotype identification, and marker-assisted selection of common beans.     

The collaborative cross: a recombinant inbred mouse population for the systems genetic era

The mouse is the most extensively used mammalian model for biomedical and aging research, and an extensive catalogue of laboratory resources is available to support research using mice: classical inbred lines, genetically modified mice (knockouts, transgenics, and humanized mice), selectively bred lines, consomics, congenics, recombinant inbred panels, outbred and heterogeneous stocks, and an expanding set of wild-derived strains. However, these resources were not designed or intended to model the heterogeneous human population or for a systematic analysis of phenotypic effects due to random combinations of uniformly distributed natural variants. The Collaborative Cross (CC) is a large panel of recently established multiparental recombinant inbred mouse lines specifically designed to overcome the limitations of existing mouse genetic resources for analysis of phenotypes caused by combinatorial allele effects. The CC models the complexity of the human genome and supports analyses of common human diseases with complex etiologies originating through interactions between allele combinations and the environment. The CC is the only mammalian resource that has high and uniform genomewide genetic variation effectively randomized across a large, heterogeneous, and infinitely reproducible population. The CC supports data integration across environmental and biological perturbations and across space (different labs) and time.     

Initial characterization of STS markers in the LSXSS series of recombinant inbred strains

We are mapping quantitative trait loci (QTLs) that influence ethanol-induced anesthesia (sleep time) in the Long-sleep (LS) and Short-sleep (SS) slected lines of mice. Fifty microsatellite-STS markers were initially screened for simple-sequence length polymorphisms between the LS and SS lines. Nineteen markers were polymorphic. Eleven markers unequivocally differentiated the LS and SS lines and were used to establish strain distribution patterns for the LSXSS series of recombinant inbred strains. Five markers each accounted for at least 5% of sleep-time genetic variance among the RI strains. Linkage of provisional QTLs detected among RIs will be confirmed or disproved in a large F₂ population. This ongoing QTL-mapping project eventually will result in a strain distribution pattern for the LSXSS RI series with an average marker spacing of 5 centimorgans.