Imputation of single-nucleotide polymorphisms in inbred mice using local phylogeny

We present full-genome genotype imputations for 100 classical laboratory mouse strains, using a novel method. Using genotypes at 549,683 SNP loci obtained with the Mouse Diversity Array, we partitioned the genome of 100 mouse strains into 40,647 intervals that exhibit no evidence of historical recombination. For each of these intervals we inferred a local phylogenetic tree. We combined these data with 12 million loci with sequence variations recently discovered by whole-genome sequencing in a common subset of 12 classical laboratory strains. For each phylogenetic tree we identified strains sharing a leaf node with one or more of the sequenced strains. We then imputed high- and medium-confidence genotypes for each of 88 nonsequenced genomes. Among inbred strains, we imputed 92% of SNPs genome-wide, with 71% in high-confidence regions. Our method produced 977 million new genotypes with an estimated per-SNP error rate of 0.083% in high-confidence regions and 0.37% genome-wide. Our analysis identified which of the 88 nonsequenced strains would be the most informative for improving full-genome imputation, as well as which additional strain sequences will reveal more new genetic variants. Imputed sequences and quality scores can be downloaded and visualized online.     

SNP genotyping for the genetic monitoring of laboratory mice by using a microarray-based method with dualcolour fluorescence hybridisation

Ensuring the genetic homogeneity of the mice used in laboratory experiments contributes to the Reduction aspect of the Three Rs, by maximising the quality of the data obtained from any animals that are used for these purposes, and ultimately reducing the numbers of animals used. Single nucleotide polymorphism (SNP) genotyping is especially suitable for use in the analysis of the genetic purity of model organisms such as the mouse, because bi-allelic markers remain fully informative when used to characterise crosses between inbred strains. Here, we attempted to apply a microarray-based method for a SNP marker to monitor the genetic quality of inbred mouse strains, so as to validate the reliability, stability and applicability of this SNP genotyping panel. The amplified PCR products containing four different SNP loci from four inbred mouse strains were spotted and immobilised onto amino-modified glass slides to generate a microarray. This was then interrogated through hybridisation with dual-colour probes, to determine the SNP genotypes of each sample. The results indicated that this microarray-based method could effectively determine the genotypes of the four selected SNPs with a high degree of accuracy. We have developed a new SNP genotyping technique for effective use in the genetic monitoring of inbred mouse strains.     

Development of a SNP genotyping panel for genetic monitoring of the laboratory mouse

We have developed a genotyping system for detecting genetic contamination in the laboratory mouse based on assaying single-nucleotide polymorphism (SNP) markers positioned on all autosomes and the X chromosome. This system provides a fast, reliable, and cost-effective way for genetic monitoring, while maintaining a very high degree of confidence. We describe the allelic distribution of 235 SNPs in 48 mouse strains, thereby creating a database of polymorphisms useful for genotyping purposes. The SNP markers used in this study were chosen from publicly available SNP databases. Four genotyping methods were evaluated, and dynamic two-tube allele-specific PCR assays were developed for each marker and tested on a set of 48 inbred mouse strains. The minimal number of assays sufficient to distinguish groups consisting of different numbers of mouse strains was estimated, and a panel of 28 SNPs sufficient to distinguish virtually all of the inbred strains tested was selected. Amplifluor SNP detection assays were developed for these markers and tested on an extended list of 96 strains. This panel was used as a genetic quality control approach to monitor the genotypes of nearly 300 inbred, wild-derived, congenic, consomic, and recombinant inbred strains maintained at The Jackson Laboratory. We have concluded that this marker panel is sufficient for genetic contamination monitoring in colonies containing a large number of genetically diverse mouse strains and that reduced versions of the panel could be implemented in facilities housing a lower number of strains.     

Genomic prediction using imputed whole-genome sequence data in Holstein Friesian cattle

Background

In contrast to currently used single nucleotide polymorphism (SNP) panels, the use of whole-genome sequence data is expected to enable the direct estimation of the effects of causal mutations on a given trait. This could lead to higher reliabilities of genomic predictions compared to those based on SNP genotypes. Also, at each generation of selection, recombination events between a SNP and a mutation can cause decay in reliability of genomic predictions based on markers rather than on the causal variants. Our objective was to investigate the use of imputed whole-genome sequence genotypes versus high-density SNP genotypes on (the persistency of) the reliability of genomic predictions using real cattle data.

Methods

Highly accurate phenotypes based on daughter performance and Illumina BovineHD Beadchip genotypes were available for 5503 Holstein Friesian bulls. The BovineHD genotypes (631,428 SNPs) of each bull were used to impute whole-genome sequence genotypes (12,590,056 SNPs) using the Beagle software. Imputation was done using a multi-breed reference panel of 429 sequenced individuals. Genomic estimated breeding values for three traits were predicted using a Bayesian stochastic search variable selection (BSSVS) model and a genome-enabled best linear unbiased prediction model (GBLUP). Reliabilities of predictions were based on 2087 validation bulls, while the other 3416 bulls were used for training.

Results

Prediction reliabilities ranged from 0.37 to 0.52. BSSVS performed better than GBLUP in all cases. Reliabilities of genomic predictions were slightly lower with imputed sequence data than with BovineHD chip data. Also, the reliabilities tended to be lower for both sequence data and BovineHD chip data when relationships between training animals were low. No increase in persistency of prediction reliability using imputed sequence data was observed.

Conclusions

Compared to BovineHD genotype data, using imputed sequence data for genomic prediction produced no advantage. To investigate the putative advantage of genomic prediction using (imputed) sequence data, a training set with a larger number of individuals that are distantly related to each other and genomic prediction models that incorporate biological information on the SNPs or that apply stricter SNP pre-selection should be considered.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0149-x) contains supplementary material, which is available to authorized users.     

A comprehensive mouse IBD database for the efficient localization of quantitative trait loci

Traditional fine-mapping approaches in mouse genetics that go from a linkage region to a candidate gene are very costly and time consuming. Shared ancestry regions, along with the combination of genetics and genomics approaches, provide a powerful tool to shorten the time and effort required to identify a causative gene. In this article we present a novel methodology that predicts IBD (identical by descent) regions between pairs of inbred strains using single nucleotide polymorphism (SNP) maps. We have validated this approach by comparing the IBD regions, estimated using different algorithms, to the results derived using the sequence information in the strains present in the Celera Mouse Database. We showed that based on the current publicly available SNP genotypes, large IBD regions (>1 Mb) can be identified successfully. By assembling a list of 21,514 SNPs in 61 common inbred strains, we inferred IBD regions between all pairs of strains and confirmed, for the first time, that existing quantitative trait genes (QTG) and susceptibility genes all lie outside of IBD regions. We also illustrated how knowledge of IBD structures can be applied to strain selection for future crosses. We have made our results available for data mining and download through a public website ( ). Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users.     

Evaluation of measures of correctness of genotype imputation in the context of genomic prediction: a review of livestock applications

In livestock, many studies have reported the results of imputation to 50k single nucleotide polymorphism (SNP) genotypes for animals that are genotyped with low-density SNP panels. The objective of this paper is to review different measures of correctness of imputation, and to evaluate their utility depending on the purpose of the imputed genotypes. Across studies, imputation accuracy, computed as the correlation between true and imputed genotypes, and imputation error rates, that counts the number of incorrectly imputed alleles, are commonly used measures of imputation correctness. Based on the nature of both measures and results reported in the literature, imputation accuracy appears to be a more useful measure of the correctness of imputation than imputation error rates, because imputation accuracy does not depend on minor allele frequency (MAF), whereas imputation error rate depends on MAF. Therefore imputation accuracy can be better compared across loci with different MAF. Imputation accuracy depends on the ability of identifying the correct haplotype of a SNP, but many other factors have been identified as well, including the number of genotyped immediate ancestors, the number of animals with genotypes at the high-density panel, the SNP density on the low- and high-density panel, the MAF of the imputed SNP and whether imputed SNP are located at the end of a chromosome or not. Some of these factors directly contribute to the linkage disequilibrium between imputed SNP and SNP on the low-density panel. When imputation accuracy is assessed as a predictor for the accuracy of subsequent genomic prediction, we recommend that: (1) individual-specific imputation accuracies should be used that are computed after centring and scaling both true and imputed genotypes; and (2) imputation of gene dosage is preferred over imputation of the most likely genotype, as this increases accuracy and reduces bias of the imputed genotypes and the subsequent genomic predictions.     

Quantitative trait locus and haplotype mapping in closely related inbred strains identifies a locus for open field behavior

Quantitative trait locus (QTL) mapping in the mouse typically utilizes inbred strains that exhibit significant genetic and phenotypic diversity. The development of dense SNP panels in a large number of inbred strains has eliminated the need to maximize genetic diversity in QTL studies as plenty of SNP markers are now available for almost any combination of strains. We conducted a QTL mapping experiment using both a backcross (N₂) and an intercross (F₂) between two genetically similar inbred mouse strains: C57BL/6J (B6) and C57L/J (C57). A set of additive QTLs for activity behaviors was identified on Chrs 1, 9, 13, and 15. We also identified additive QTLs for anxiety-related behaviors on Chrs 7, 9, and 16. A QTL on Chr 11 is sex-specific, and we revealed pairwise interactions between QTLs on Chrs 1 and 13 and Chrs 10 and 18. The Chr 9 activity QTL accounts for the largest amount of phenotypic variance and was not present in our recent analysis of a B6 × C58/J (C58) intercross (Bailey et al. in Genes Brain Behav 7:761–769, 2008). To narrow this QTL interval, we used a dense SNP haplotype map with over 7 million real and imputed SNP markers across 74 inbred mouse strains (Szatkiewicz et al. in Mamm Genome 19(3):199–208, 2008). Evaluation of shared and divergent haplotype blocks among B6, C57, and C58 strains narrowed the Chr 9 QTL interval considerably and highlights the utility of QTL mapping in closely related inbred strains.     

Fine Mapping in 94 Inbred Mouse Strains Using a High-Density Haplotype Resource

The genetics of phenotypic variation in inbred mice has for nearly a century provided a primary weapon in the medical research arsenal. A catalog of the genetic variation among inbred mouse strains, however, is required to enable powerful positional cloning and association techniques. A recent whole-genome resequencing study of 15 inbred mouse strains captured a significant fraction of the genetic variation among a limited number of strains, yet the common use of hundreds of inbred strains in medical research motivates the need for a high-density variation map of a larger set of strains. Here we report a dense set of genotypes from 94 inbred mouse strains containing 10.77 million genotypes over 121,433 single nucleotide polymorphisms (SNPs), dispersed at 20-kb intervals on average across the genome, with an average concordance of 99.94% with previous SNP sets. Through pairwise comparisons of the strains, we identified an average of 4.70 distinct segments over 73 classical inbred strains in each region of the genome, suggesting limited genetic diversity between the strains. Combining these data with genotypes of 7570 gap-filling SNPs, we further imputed the untyped or missing genotypes of 94 strains over 8.27 million Perlegen SNPs. The imputation accuracy among classical inbred strains is estimated at 99.7% for the genotypes imputed with high confidence. We demonstrated the utility of these data in high-resolution linkage mapping through power simulations and statistical power analysis and provide guidelines for developing such studies. We also provide a resource of in silico association mapping between the complex traits deposited in the Mouse Phenome Database with our genotypes. We expect that these resources will facilitate effective designs of both human and mouse studies for dissecting the genetic basis of complex traits.PHENOTYPIC variation among inbred mouse strains exposed to a disease-causing agent (be it genetic, infectious, or environmental) provides potential insight into human disease processes that often cannot be practically achieved through direct human studies. Indeed, hundreds of phenotype measurements related to human diseases are available for dozens of inbred strains in common use over the past 50–100 years (Bogue et al. 2007; Grubb et al. 2009). As with the direct study of chronic disease in humans, key steps toward determining the genetic underpinnings of this phenotypic variation are to develop a catalog of the genetic variation among inbred mouse strains and to interpret the structure of variation patterns across the strains. Recent advances in high-throughput genotyping and DNA resequencing technologies are making it possible to rapidly uncover the genetic variation maps of many model organisms (Lindblad-Toh et al. 2005; Mackay and Anholt 2006; Borevitz et al. 2007; Frazer et al. 2007; International Hapmap Consortium 2007; Star Consortium 2008). A recent whole-genome resequencing study of 15 inbred mouse strains captured a significant fraction of the genetic variation among a limited number of strains, allowing researchers to infer patterns of genetic variation and to identify the ancestral origin of the genetic variation (Frazer et al. 2007; Yang et al. 2007). Yet the availability and common experimental employment of hundreds of inbred strains, including >190 stocks available from the Jackson Laboratory, motivates the need for a high-density variation map for a larger set of strains. We have assembled the Mouse HapMap, a resource consisting of a dense set of genotypes for a total of 138,980 unique biallelic single nucleotide polymorphisms (SNPs) in 94 inbred mouse strains at an average spacing of 20 kb on chromosomes 1–19 and X.This resource is ideal for performing high-resolution mapping studies under QTL peaks. We evaluate the feasibility and effectiveness of such studies by examining a typical study from the Mouse Phenome Database (MPD) (Bogue et al. 2007; Grubb et al. 2009) (http://www.jax.org/phenome) and measure the statistical power to detect genetic associations in regions of various sizes. We provide several resources to the mouse genetics community for supporting such studies and a webserver that can estimate the significance threshold, compute the statistical power of a proposed study, and perform in the fine mapping of measured phenotypes. In addition, we provide a database of associations for all phenotypes contained in the MPD. The web resources are available at http://mouse.cs.ucla.edu/.     

Imputation of genotypes from low- to high-density genotyping platforms and implications for genomic selection

The objective of this study was to quantify the accuracy achievable from imputing genotypes from a commercially available low-density marker panel (2730 single nucleotide polymorphisms (SNPs) following edits) to a commercially available higher density marker panel (51 602 SNPs following edits) in Holstein-Friesian cattle using Beagle, a freely available software package. A population of 764 Holstein-Friesian animals born since 2006 were used as the test group to quantify the accuracy of imputation, all of which had genotypes for the high-density panel; only SNPs on the low-density panel were retained with the remaining SNPs to be imputed. The reference population for imputation consisted of 4732 animals born before 2006 also with genotypes on the higher density marker panel. The concordance between the actual and imputed genotypes in the test group of animals did not vary across chromosomes and was on average 95%; the concordance between actual and imputed alleles was, on average, 97% across all SNPs. Genomic predictions were undertaken across a range of production and functional traits for the 764 test group animals using either their real or imputed genotypes. Little or no mean difference in the genomic predictions was evident when comparing direct genomic values (DGVs) using real or imputed genotypes. The average correlation between the DGVs estimated using the real or imputed genotypes for the 15 traits included in the Irish total merit index was 0.97 (range of 0.92 to 0.99), indicating good concordance between proofs from real or imputed genotypes. Results show that a commercially available high-density marker panel can be imputed from a commercially available lower density marker panel, which will also have a lower cost, thereby facilitating a reduction in the cost of genomic selection. Increased available numbers of genotyped and phenotyped animals also has implications for increasing the accuracy of genomic prediction in the entire population and thus genetic gain using genomic selection.     

Use of a dense single nucleotide polymorphism map for in silico mapping in the mouse

Rapid expansion of available data, both phenotypic and genotypic, for multiple strains of mice has enabled the development of new methods to interrogate the mouse genome for functional genetic perturbations. In silico mapping provides an expedient way to associate the natural diversity of phenotypic traits with ancestrally inherited polymorphisms for the purpose of dissecting genetic traits. In mouse, the current single nucleotide polymorphism (SNP) data have lacked the density across the genome and coverage of enough strains to properly achieve this goal. To remedy this, 470,407 allele calls were produced for 10,990 evenly spaced SNP loci across 48 inbred mouse strains. Use of the SNP set with statistical models that considered unique patterns within blocks of three SNPs as an inferred haplotype could successfully map known single gene traits and a cloned quantitative trait gene. Application of this method to high-density lipoprotein and gallstone phenotypes reproduced previously characterized quantitative trait loci (QTL). The inferred haplotype data also facilitates the refinement of QTL regions such that candidate genes can be more easily identified and characterized as shown for adenylate cyclase 7.     

Haplotype-based genetics in mice and rats

Haplotype blocks are conceptually defined as genomic segments harbouring sets of coupled polymorphisms that reflect a common ancestral origin. Experimentally, however, haplotype blocks are characterized using computational algorithms based on incomplete inventories of polymorphisms. Haplotype blocks and their deduced strain-distribution patterns are considered to be extremely powerful for use in genetic association and mapping experiments in laboratory mice and rats. However, recent high-density SNP screening in commonly used mouse inbred strains reveals a complex pattern, suggesting that the current expectations for the use of haplotype blocks in genetic mapping will have to be revisited.     

Genomic copy number variation in Mus musculus

Background

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

Results

We found 9634 putative autosomal CNVs across the samples affecting 6.87 % of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

Conclusion

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1713-z) contains supplementary material, which is available to authorized users.     

Genomic landscape of rat strain and substrain variation

Background

Since the completion of the rat reference genome in 2003, whole-genome sequencing data from more than 40 rat strains have become available. These data represent the broad range of strains that are used in rat research including commonly used substrains. Currently, this wealth of information cannot be used to its full extent, because the variety of different variant calling algorithms employed by different groups impairs comparison between strains. In addition, all rat whole genome sequencing studies to date used an outdated reference genome for analysis (RGSC3.4 released in 2004).

Results

Here we present a comprehensive, multi-sample and uniformly called set of genetic variants in 40 rat strains, including 19 substrains. We reanalyzed all primary data using a recent version of the rat reference assembly (RGSC5.0 released in 2012) and identified over 12 million genomic variants (SNVs, indels and structural variants) among the 40 strains. 28,318 SNVs are specific to individual substrains, which may be explained by introgression from other unsequenced strains and ongoing evolution by genetic drift. Substrain SNVs may have a larger predicted functional impact compared to older shared SNVs.

Conclusions

In summary we present a comprehensive catalog of uniformly analyzed genetic variants among 40 widely used rat inbred strains based on the RGSC5.0 assembly. This represents a valuable resource, which will facilitate rat functional genomic research. In line with previous observations, our genome-wide analyses do not show evidence for contribution of multiple ancestral founder rat subspecies to the currently used rat inbred strains, as is the case for mouse. In addition, we find that the degree of substrain variation is highly variable between strains, which is of importance for the correct interpretation of experimental data from different labs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1594-1) contains supplementary material, which is available to authorized users.     

Integrating QTL and high-density SNP analyses in mice to identify Insig2 as a susceptibility gene for plasma cholesterol levels

The use of inbred strains of mice to dissect the genetic complexity of common diseases offers a viable alternative to human studies, given the control over experimental parameters that can be exercised. Central to efforts to map susceptibility loci for common diseases in mice is a comprehensive map of DNA variation among the common inbred strains of mice. Here we present one of the most comprehensive high-density, single nucleotide polymorphism (SNP) maps of mice constructed to date. This map consists of 10,350 SNPs genotyped in 62 strains of inbred mice. We demonstrate the utility of these data via a novel integrative genomics approach to mapping susceptibility loci for complex traits. By integrating in silico quantitative trait locus (QTL) mapping with progressive QTL mapping strategies in segregating mouse populations that leverage large-scale mapping of the genetic determinants of gene expression traits, we not only facilitate identification of candidate quantitative trait genes, but also protect against spurious associations that can arise in genetic association studies due to allelic association among unlinked markers. Application of this approach to our high-density SNP map and two previously described F2 crosses between strains C57BL/6J (B6) and DBA/2J and between B6 ApoE(-/-) and C3H/HeJ ApoE(-/-) results in the identification of Insig2 as a strong candidate susceptibility gene for total plasma cholesterol levels.     

Mouse inbred strain sequence information and yin-yang crosses for quantitative trait locus fine mapping

The shared ancestry of mouse inbred strains, together with the availability of sequence and phenotype information, is a resource that can be used to map quantitative trait loci (QTL). The difficulty in using only sequence information lies in the fact that in most instances the allelic state of the QTL cannot be unambiguously determined in a given strain. To overcome this difficulty, the performance of multiple crosses between various inbred strains has been proposed. Here we suggest and evaluate a general approach, which consists of crossing the two strains used initially to map the QTL and any new strain. We have termed these crosses "yin-yang," because they are complementary in nature as shown by the fact that the QTL will necessarily segregate in only one of the crosses. We used the publicly available SNP database of chromosome 16 to evaluate the mapping resolution achievable through this approach. Although on average the improvement of mapping resolution using only four inbred strains was relatively small (i.e., reduction of the QTL-containing interval by half at most), we found a great degree of variability among different regions of chromosome 16 with regard to mapping resolution. This suggests that with a large number of strains in hand, selecting a small number of strains may provide a significant contribution to the fine mapping of QTL.     

单管双向等位基因专一性扩增方法在近交系小鼠遗传检测中的应用

目的将新近建立的单管双向等位基因专一性扩增(single-tube bi-directional allele specific amplification,SB-ASA)方法用于分析近交系小鼠基因组中的单核苷酸多态性(SNP)。方法以5个近交系小鼠为研究对象,采用SB-ASA方法对其16个SNP位点进行检测,并通过双盲实验和测序验证该方法的可靠性;且考察了该方法中PCR反应各成分及扩增条件对结果的影响。结果16个SNP位点,SB-ASA都成功地对5个品系小鼠进行了分型,与测序结果完全一致;双盲实验结果显示通过3个SNP位点即可鉴别5个品系。结论SB-ASA方法可用于近交系小鼠SNP的遗传检测,可望作为一种新的分子生物学遗传检测方法推广应用。     

1 + 1 = 3: Development and Validation of a SNP-Based Algorithm to Identify Genetic Contributions from Three Distinct Inbred Mouse Strains

State-of-the-art, genome-wide assessment of mouse genetic background uses single nucleotide polymorphism (SNP) PCR. As SNP analysis can use multiplex testing, it is amenable to high-throughput analysis and is the preferred method for shared resource facilities that offer genetic background assessment of mouse genomes. However, a typical individual SNP query yields only two alleles (A vs. B), limiting the application of this methodology to distinguishing contributions from no more than two inbred mouse strains. By contrast, simple sequence length polymorphism (SSLP) analysis yields multiple alleles but is not amenable to high-throughput testing. We sought to devise a SNP-based technique to identify donor strain origins when three distinct mouse strains potentially contribute to the genetic makeup of an individual mouse. A computational approach was used to devise a three-strain analysis (3SA) algorithm that would permit identification of three genetic backgrounds while still using a binary-output SNP platform. A panel of 15 mosaic mice with contributions from BALB/c, C57Bl/6, and DBA/2 genetic backgrounds was bred and analyzed using a genome-wide SNP panel using 1449 markers. The 3SA algorithm was applied and then validated using SSLP. The 3SA algorithm assigned 85% of 1449 SNPs as informative for the C57Bl/6, BALB/c, or DBA/2 backgrounds, respectively. Testing the panel of 15 F2 mice, the 3SA algorithm predicted donor strain origins genome-wide. Donor strain origins predicted by the 3SA algorithm correlated perfectly with results from individual SSLP markers located on five different chromosomes (n=70 tests). We have established and validated an analysis algorithm based on binary SNP data that can successfully identify the donor strain origins of chromosomal regions in mice that are bred from three distinct inbred mouse strains.     

Efficient control of population structure in model organism association mapping

Genomewide association mapping in model organisms such as inbred mouse strains is a promising approach for the identification of risk factors related to human diseases. However, genetic association studies in inbred model organisms are confronted by the problem of complex population structure among strains. This induces inflated false positive rates, which cannot be corrected using standard approaches applied in human association studies such as genomic control or structured association. Recent studies demonstrated that mixed models successfully correct for the genetic relatedness in association mapping in maize and Arabidopsis panel data sets. However, the currently available mixed-model methods suffer from computational inefficiency. In this article, we propose a new method, efficient mixed-model association (EMMA), which corrects for population structure and genetic relatedness in model organism association mapping. Our method takes advantage of the specific nature of the optimization problem in applying mixed models for association mapping, which allows us to substantially increase the computational speed and reliability of the results. We applied EMMA to in silico whole-genome association mapping of inbred mouse strains involving hundreds of thousands of SNPs, in addition to Arabidopsis and maize data sets. We also performed extensive simulation studies to estimate the statistical power of EMMA under various SNP effects, varying degrees of population structure, and differing numbers of multiple measurements per strain. Despite the limited power of inbred mouse association mapping due to the limited number of available inbred strains, we are able to identify significantly associated SNPs, which fall into known QTL or genes identified through previous studies while avoiding an inflation of false positives. An R package implementation and webserver of our EMMA method are publicly available.     

Identification of a genetic contamination in a commercial mouse strain using two panels of polymorphic markers

Rapid detection of genetic contamination is critical in mouse studies involving inbred strains. During a Quantitative Trait Locus (QTL) study using simple sequence length polymorphism (SSLP) markers, we noticed heterozygosity at some loci of a commercially available inbred C57BL/6N mouse strain, suggesting a contamination by another mouse strain. A panel of 100 single-nucleotide polymorphism (SNP) markers was used to confirm and specify the genetic contamination suspected. Retrospective analyses demonstrated that the contamination took place as early as autumn 2003 and has persisted ever since at a fairly constant level. Contaminating alleles most probably originated from a DBA strain. Our data demonstrate the suitability of SNP markers for rapid detection and identification of the source of genetic contamination. Further, our results show the importance of a state-of-the-art genetic monitoring of the authenticity of murine inbred strains.