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.     

High‐precision genetic mapping of behavioral traits in the diversity outbred mouse population

Historically our ability to identify genetic variants underlying complex behavioral traits in mice has been limited by low mapping resolution of conventional mouse crosses. The newly developed Diversity Outbred (DO) population promises to deliver improved resolution that will circumvent costly fine‐mapping studies. The DO is derived from the same founder strains as the Collaborative Cross (CC), including three wild‐derived strains. Thus the DO provides more allelic diversity and greater potential for discovery compared to crosses involving standard mouse strains. We have characterized 283 male and female DO mice using open‐field, light–dark box, tail‐suspension and visual‐cliff avoidance tests to generate 38 behavioral measures. We identified several quantitative trait loci (QTL) for these traits with support intervals ranging from 1 to 3 Mb in size. These intervals contain relatively few genes (ranging from 5 to 96). For a majority of QTL, using the founder allelic effects together with whole genome sequence data, we could further narrow the positional candidates. Several QTL replicate previously published loci. Novel loci were also identified for anxiety‐ and activity‐related traits. Half of the QTLs are associated with wild‐derived alleles, confirming the value to behavioral genetics of added genetic diversity in the DO. In the presence of wild‐alleles we sometimes observe behaviors that are qualitatively different from the expected response. Our results demonstrate that high‐precision mapping of behavioral traits can be achieved with moderate numbers of DO animals, representing a significant advance in our ability to leverage the mouse as a tool for behavioral genetics .     

Architecture of energy balance traits in emerging lines of the Collaborative Cross

The potential utility of the Collaborative Cross (CC) mouse resource was evaluated to better understand complex traits related to energy balance. A primary focus was to examine if genetic diversity in emerging CC lines (pre-CC) would translate into equivalent phenotypic diversity. Second, we mapped quantitative trait loci (QTL) for 15 metabolism- and exercise-related phenotypes in this population. We evaluated metabolic and voluntary exercise traits in 176 pre-CC lines, revealing phenotypic variation often exceeding that seen across the eight founder strains from which the pre-CC was derived. Many phenotypic correlations existing within the founder strains were no longer significant in the pre-CC population, potentially representing reduced linkage disequilibrium (LD) of regions harboring multiple genes with effects on energy balance or disruption of genetic structure of extant inbred strains with substantial shared ancestry. QTL mapping revealed five significant and eight suggestive QTL for body weight (Chr 4, 7.54 Mb; CI 3.32-10.34 Mb; Bwq14), body composition, wheel running (Chr 16, 33.2 Mb; CI 32.5-38.3 Mb), body weight change in response to exercise (1: Chr 6, 77.7Mb; CI 72.2-83.4 Mb and 2: Chr 6, 42.8 Mb; CI 39.4-48.1 Mb), and food intake during exercise (Chr 12, 85.1 Mb; CI 82.9-89.0 Mb). Some QTL overlapped with previously mapped QTL for similar traits, whereas other QTL appear to represent novel loci. These results suggest that the CC will be a powerful, high-precision tool for examining the genetic architecture of complex traits such as those involved in regulation of energy balance.     

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.     

The Collaborative Cross, a groundbreaking tool to tackle complex traits

Complex traits, like the susceptibility to common diseases, are controlled by numerous genomic regions which individual effect is generally weak. These observations led geneticists to develop an experimental system to dissect the genetic of complex traits in the mouse. The Collaborative Cross (CC) is a genetic reference population of over 300 inbred lines derived from eight inbred strains of three Mus musculus sub-species that captures 90% of the genetic variation known in the mouse genome. We present here the generation and the characteristics of the CC and we report the results of the first experiments with partially inbred CC lines.     

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.     

Lz-0 × Berkeley: a new Arabidopsis recombinant inbred line population for the mapping of complex traits

This study describes the generation and test of a genetic resource suited to identify determinants of cell biological traits in plants. The use of quantitative trait loci (QTL) mapping for a better genetic understanding of cell biological traits is still at an early stage, even for biotechnologically important cell properties, such as the dimensions of fiber cells. A common strategy, the mapping of QTLs in recombinant inbred line (RIL) populations, is limited by the fact that the existing RIL populations exploit only a small fraction of the existing natural variation. Here, we report the mapping of QTLs impacting on the length of fiber cells in Arabidopsis inflorescence stems in a newly generated RIL population derived from a cross between the accessions Berkeley and the little known Lz-0. Through inbreeding of individual F₂ plants, a total of 159 new F₈ lines were produced and genotyped with a set of 49 single nucleotide polymorphism markers. The population was successfully used not only for the mapping of three QTLs controlling fiber length, but also to map five QTL controlling flowering time under short and long-day conditions. Our study demonstrates the usefulness of this new genetic resource by mapping in it QTLs underlying a poorly explored cellular trait as well as an already better explored regulatory pathway. The new RIL population and an online platform for the continuous supplementation of genetic markers will be generally available to substantially broaden the genetic diversity through which loci with impact on plant quantitative traits can be identified.     

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     

The evolution of genetic architecture. I. Diversification of genetic backgrounds by genetic drift

The genetic architecture of a phenotype plays a critical role in determining phenotypic evolution through its effects on patterns of genetic variation. Genetic architecture is often considered to be constant in evolutionary quantitative genetic models. However, genetic architecture may be variable and itself evolve when there are dominance and epistatic interactions among alleles at the same and different loci, respectively. The evolution of genetic architecture by genetic drift is examined here by testing the breeding value of four standard inbred mouse strains mated across a set of 26 related recombinant quasi-inbred (RqI) lines generated from the intercross of the Large (LG/J) and Small (SM/J) inbred mouse strains. Phenotypes of interest include age-specific body weights, growth, and adult body composition. If the genetic architecture of these traits has differentiated by genetic drift during the production of the RqI strains, we should observe interactions between tester strain and RqI strain. The breeding values of the tester strains will change relative to one another depending on which RqI strain they are crossed to. The study included an average of 15.1 offspring per cross, over a total of 100 different crosses. Multivariate and univariate analyses of variance indicate that there is strongly significant interaction for all traits. Interaction is more pronounced in males than in females and accounted for an average of about 40% of the explained variation in males and 30% in females. These results indicate that the genetic architecture of these traits has differentiated by genetic drift in the RqI strains since their isolation from a common founder population. Further analysis indicates that this differentiation results in changes in the order of tester strain effects so that common patterns of selection in these differentiated populations could result in the fixation of different alleles.     

The genetic architecture of amino acids dissection by association and linkage analysis in maize

Amino acids are both constituents of proteins, providing the essential nutrition for humans and animals, and signalling molecules regulating the growth and development of plants. Most cultivars of maize are deficient in essential amino acids such as lysine and tryptophan. Here, we measured the levels of 17 different total amino acids, and created 48 derived traits in mature kernels from a maize diversity inbred collection and three recombinant inbred line (RIL) populations. By GWAS, 247 and 281 significant loci were identified in two different environments, 5.1 and 4.4 loci for each trait, explaining 7.44% and 7.90% phenotypic variation for each locus in average, respectively. By linkage mapping, 89, 150 and 165 QTLs were identified in B73/By804, Kui3/B77 and Zong3/Yu87‐1 RIL populations, 2.0, 2.7 and 2.8 QTLs for each trait, explaining 13.6%, 16.4% and 21.4% phenotypic variation for each QTL in average, respectively. It implies that the genetic architecture of amino acids is relative simple and controlled by limited loci. About 43.2% of the loci identified by GWAS were verified by expression QTL, and 17 loci overlapped with mapped QTLs in the three RIL populations. GRMZM2G015534, GRMZM2G143008 and one QTL were further validated using molecular approaches. The amino acid biosynthetic and catabolic pathways were reconstructed on the basis of candidate genes proposed in this study. Our results provide insights into the genetic basis of amino acid biosynthesis in maize kernels and may facilitate marker‐based breeding for quality protein maize.     

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.     

The relationship between parental genetic or phenotypic divergence and progeny variation in the maize nested association mapping population

Appropriate selection of parents for the development of mapping populations is pivotal to maximizing the power of quantitative trait loci detection. Trait genotypic variation within a family is indicative of the family's informativeness for genetic studies. Accurate prediction of the most useful parental combinations within a species would help guide quantitative genetics studies. We tested the reliability of genotypic and phenotypic distance estimators between pairs of maize inbred lines to predict genotypic variation for quantitative traits within families derived from biparental crosses. We developed 25 families composed of ~200 random recombinant inbred lines each from crosses between a common reference parent inbred, B73, and 25 diverse maize inbreds. Parents and families were evaluated for 19 quantitative traits across up to 11 environments. Genetic distances (GDs) among parents were estimated with 44 simple sequence repeat and 2303 single-nucleotide polymorphism markers. GDs among parents had no predictive value for progeny variation, which is most likely due to the choice of neutral markers. In contrast, we observed for about half of the traits measured a positive correlation between phenotypic parental distances and within-family genetic variance estimates. Consequently, the choice of promising segregating populations can be based on selecting phenotypically diverse parents. These results are congruent with models of genetic architecture that posit numerous genes affecting quantitative traits, each segregating for allelic series, with dispersal of allelic effects across diverse genetic material. This architecture, common to many quantitative traits in maize, limits the predictive value of parental genotypic or phenotypic values on progeny variance.     

Patterns of variation in desiccation resistance in a set of recombinant inbred lines in Drosophila melanogaster

Desiccation, resulting from extremely dry environmental conditions, is a serious obstacle to the survival of organisms. Water is vital for the maintenance of intracellular structure and prevents the irreversible formation of aggregates, an occurrence leading to loss of cellular function. To characterize genetic variation in desiccation stress resistance (DSR) in Drosophila melanogaster Meigen, an intercontinental set of recombinant inbred lines (RIL) is used. Flies are exposed to a low humidity environment (<10% relative humidity) at a constant temperature of 25 °C. Desiccation stress resistance is higher in RIL derived from a backcross to the parental stock sensitive to heat stress (from Denmark) than in RIL derived from the reciprocal backcross to the heat‐stress resistant stock (from Australia). Composite interval mapping reveals significant quantitative trail loci (QTL) for DSR in the set of RIL. Both major and minor effects QTL are detected, suggesting a complex genetic architecture. When compared with a previous investigation performed on the same set of RIL, the present study indicates that not all traits of resistance to environmental stressors are affected in the same direction by segregating co‐localized QTL.     

The Collaborative Cross,developing a resource for mammalian systems genetics: A status report of the Wellcome Trust cohort

We report on the progress of a project funded by the Wellcome Trust to produce over 100 recombinant inbred mouse lines as part of the Collaborative Cross (CC) genetic reference panel. These new strains of mice are being derived from a set of eight genetically diverse founders. The genomes of the finished strains will be mosaics of the founder strains’ genomes with a high density of independent recombination breakpoints. The CC mice will be available for distribution free of any intellectual property constraints to serve as a community resource for systems genetics studies.     

Linkage analysis of quantitative trait loci in multiple line crosses

Simple line crosses, for example, backcross and F₂, are commonly used in mapping quantitative trait loci (QTL). However, these simple crosses are rarely used alone in commercial plant breeding; rather, crosses involving multiple inbred lines or several simple crosses but connected by shared inbred lines may be common in plant breeding. Mapping QTL using crosses of multiple lines is more relevant to plant breeding. Unfortunately, current statistical methods and computer programs of QTL mapping are all designed for simple line crosses or multiple line crosses but under a regular mating system. It is not straightforward to extend the existing methods to handle multiple line crosses under irregular and complicated mating designs. The major hurdle comes from irregular inbreeding, multiple generations, and multiple alleles. In this study, we develop a Bayesian method implemented via the Markov chain Monte Carlo (MCMC) algorithm for mapping QTL using complicated multiple line crosses. With the MCMC algorithm, we are able to draw a complete path of the gene flow from founder alleles to their descendents via a recursive process. This has greatly simplified the problem caused by irregular mating and inbreeding in the mapping population. Adopting the reversible jump MCMC algorithm, we are able to simultaneously search for multiple QTL along the genome. We can even infer the posterior distribution of the number of QTL, one of the most important parameters in QTL study. Application of the new MCMC based QTL mapping procedure is demonstrated using two different mating designs. Design I involves two inbred lines and their derived F₁, F₂, and BC populations. Design II is a half-diallel cross involving three inbred lines. The two designs appear different, but can be handled with the same robust computer program.     

The genetic architecture of Drosophila sensory bristle number

We have mapped quantitative trait loci (QTL) for Drosophila mechanosensory bristle number in six recombinant isogenic line (RIL) mapping populations, each of which was derived from an isogenic chromosome extracted from a line selected for high or low, sternopleural or abdominal bristle number and an isogenic wild-type chromosome. All RILs were evaluated as male and female F(1) progeny of crosses to both the selected and the wild-type parental chromosomes at three developmental temperatures (18 degrees, 25 degrees, and 28 degrees ). QTL for bristle number were mapped separately for each chromosome, trait, and environment by linkage to roo transposable element marker loci, using composite interval mapping. A total of 53 QTL were detected, of which 33 affected sternopleural bristle number, 31 affected abdominal bristle number, and 11 affected both traits. The effects of most QTL were conditional on sex (27%), temperature (14%), or both sex and temperature (30%). Epistatic interactions between QTL were also common. While many QTL mapped to the same location as candidate bristle development loci, several QTL regions did not encompass obvious candidate genes. These features are germane to evolutionary models for the maintenance of genetic variation for quantitative traits, but complicate efforts to understand the molecular genetic basis of variation for complex traits.     

High-resolution genetic mapping using the Mouse Diversity outbred population

The JAX Diversity Outbred population is a new mouse resource derived from partially inbred Collaborative Cross strains and maintained by randomized outcrossing. As such, it segregates the same allelic variants as the Collaborative Cross but embeds these in a distinct population architecture in which each animal has a high degree of heterozygosity and carries a unique combination of alleles. Phenotypic diversity is striking and often divergent from phenotypes seen in the founder strains of the Collaborative Cross. Allele frequencies and recombination density in early generations of Diversity Outbred mice are consistent with expectations based on simulations of the mating design. We describe analytical methods for genetic mapping using this resource and demonstrate the power and high mapping resolution achieved with this population by mapping a serum cholesterol trait to a 2-Mb region on chromosome 3 containing only 11 genes. Analysis of the estimated allele effects in conjunction with complete genome sequence data of the founder strains reduced the pool of candidate polymorphisms to seven SNPs, five of which are located in an intergenic region upstream of the Foxo1 gene.     

A large nested association mapping population for breeding and quantitative trait locus mapping in Ethiopian durum wheat

The Ethiopian plateau hosts thousands of durum wheat (Triticum turgidum subsp. durum) farmer varieties (FV) with high adaptability and breeding potential. To harness their unique allelic diversity, we produced a large nested association mapping (NAM) population intercrossing fifty Ethiopian FVs with an international elite durum wheat variety (Asassa). The Ethiopian NAM population (EtNAM) is composed of fifty interconnected bi‐parental families, totalling 6280 recombinant inbred lines (RILs) that represent both a powerful quantitative trait loci (QTL) mapping tool, and a large pre‐breeding panel. Here, we discuss the molecular and phenotypic diversity of the EtNAM founder lines, then we use an array featuring 13 000 single nucleotide polymorphisms (SNPs) to characterize a subset of 1200 EtNAM RILs from 12 families. Finally, we test the usefulness of the population by mapping phenology traits and plant height using a genome wide association (GWA) approach. EtNAM RILs showed high allelic variation and a genetic makeup combining genetic diversity from Ethiopian FVs with the international durum wheat allele pool. EtNAM SNP data were projected on the fully sequenced AB genome of wild emmer wheat, and were used to estimate pairwise linkage disequilibrium (LD) measures that reported an LD decay distance of 7.4 Mb on average, and balanced founder contributions across EtNAM families. GWA analyses identified 11 genomic loci individually affecting up to 3 days in flowering time and more than 1.6 cm in height. We argue that the EtNAM is a powerful tool to support the production of new durum wheat varieties targeting local and global agriculture.     

Out of the bottleneck: the Diversity Outcross and Collaborative Cross mouse populations in behavioral genetics research

The historical origins of classical laboratory mouse strains have led to a relatively limited range of genetic and phenotypic variation, particularly for the study of behavior. Many recent efforts have resulted in improved diversity and precision of mouse genetic resources for behavioral research, including the Collaborative Cross and Diversity Outcross population. These two populations, derived from an eight way cross of common and wild-derived strains, have high precision and allelic diversity. Behavioral variation in the population is expanded, both qualitatively and quantitatively. Variation that had once been canalized among the various inbred lines has been made amenable to genetic dissection. The genetic attributes of these complementary populations, along with advances in genetic and genomic technologies, makes a systems genetic analyses of behavior more readily tractable, enabling discovery of a greater range of neurobiological phenomena underlying behavioral variation.