Informatics resources for the Collaborative Cross and related mouse populations.

With the completion of the mouse genome sequence, it is possible to define the amount, type, and organization of the genetic variation in this species. Recent reports have provided an overview of the structure of genetic variation among classical laboratory mice. On the other hand, little is known about the structure of genetic variation among wild-derived strains with the exception of the presence of higher levels of diversity. We have estimated the sequence diversity due to substitutions and insertions/deletions among 20 inbred strains of Mus musculus, chosen to enable interpretation of the molecular variation within a clear evolutionary framework. Here, we show that the level of sequence diversity present among these strains is one to two orders of magnitude higher than the level of sequence diversity observed in the human population, and only a minor fraction of the sequence differences observed is found among classical laboratory strains. Our analyses also demonstrate that deletions are significantly more frequent than insertions. We estimate that 50% of the total variation identified in M. musculus may be recovered in intrasubspecific crosses. Alleles at variants positions can be classified into 164 strain distribution patterns, a number exceeding those reported and predicted in panels of classical inbred strains. The number of strains, the analysis of multiple loci scattered across the genome, and the mosaic nature of the genome in hybrid and classical strains contribute to the observed diversity of strain distribution patterns. However, phylogenetic analyses demonstrate that ancient polymorphisms that segregate across species and subspecies play a major role in the generation of strain distribution patterns.

The genomes of classical inbred mouse strains include genes derived from all three major subspecies of the house mouse, Mus musculus. We recently posited that genetic diversity in the immunoglobulin heavy chain (IGH) gene loci of C57BL/6 and BALB/c mice reflects differences in subspecies origin. To investigate this hypothesis, we conducted high-throughput sequencing of IGH gene rearrangements to document IGH variable (IGHV), joining (IGHJ) and diversity (IGHD) genes in four inbred wild-derived mouse strains (CAST/EiJ, LEWES/EiJ, MSM/MsJ and PWD/PhJ) and a single disease model strain (NOD/ShiLtJ), collectively representing genetic backgrounds of several major mouse subspecies. A total of 341 germline IGHV sequences were inferred in the wild-derived strains, including 247 not curated in the international ImMunoGeneTics information system. By contrast, 83/84 inferred NOD IGHV genes had previously been observed in C57BL/6 mice. Variability among the strains examined was observed for only a single IGHJ gene, involving a description of a novel allele. By contrast, unexpected variation was found in the IGHD gene loci, with four previously unreported IGHD gene sequences being documented. Very few IGHV sequences of C57BL/6 and BALB/c mice were shared with strains representing major subspecies, suggesting that their IGH loci may be complex mosaics of genes of disparate origins. This suggests a similar level of diversity is likely present in the IGH loci of other classical inbred strains. This must now be documented if we are to properly understand interstrain variation in models of antibody-mediated disease.

The probabilities that two individuals share 0, 1, or 2 alleles identical by descent (IBD) at a given genotyped marker locus are quantities of fundamental importance for disease gene and quantitative trait mapping and in family-based tests of association. Until recently, genotyped markers were sufficiently sparse that founder haplotypes could be modelled as having been drawn from a population in linkage equilibrium for the purpose of estimating IBD probabilities. However, with the advent of high-throughput single nucleotide polymorphism genotyping assays, this is no longer a reasonable assumption. Indeed, the imminent arrival of individual sequencing will enable high-density single nucleotide polymorphism genotyping on a scale for which current algorithms are not equipped. In this paper, we present a simple new model in which founder haplotypes are modelled as a Markov chain. Another important innovation is that genotyping errors are explicitly incorporated into the model. We compare results obtained using the new model to those obtained using the popular genetic linkage analysis package Merlin, with and without using the cluster model of linkage disequilibrium that is incorporated into that program. We find that the new model results in accuracy approaching that of Merlin with haplotype blocks, but achieves this with orders of magnitude faster run times. Moreover, the new algorithm scales linearly with number of markers, irrespective of density, whereas Merlin scales supralinearly. We also confirm a previous finding that ignoring linkage disequilibrium in founder haplotypes can cause errors in the calculation of IBD probabilities.

We studied the patterns of linkage disequilibrium (LD) in the human genome among three populations: African Americans, Caucasians and Ashkenazi Jews. These three populations represent admixed, outbred and isolated populations, respectively. The study examined defined chromosomal regions across the whole genome. We found that SNP allele frequencies are highly correlated between Ashkenazi Jews and Caucasians and somewhat distinct in African Americans. In addition, Ashkenazi Jews have a modest increase in LD compared with Caucasians, and both have greater LD than African Americans. The three populations differed more significantly with regard to haplotype heterogeneity. We found, as expected, that Ashkenazi Jews display the greatest extent of homogeneity and African Americans the greatest extent of heterogeneity. We found that most of the variance in LD can be attributed to the difference between regions and markers rather than to that between different population types. The average recombination rates estimated by low-resolution genetic maps can only explain a small fraction of the variance between regions. We found that LD (in terms of r(2)) decreases as a function of distance even within the so-called 'haplotype blocks'. This has significant consequences when using LD mapping for the genetic dissection of complex traits, as higher density SNP maps will be required to scan the genome. PMID: 12651872

Microsatellite loci are uniformly distributed at approximately 100-kbp intervals on all chromosomes except the chromosome Y, and genetic information about more than 9000 loci and high-throughput polymorphism analysis are now available. Taking advantage of these properties, we carried out whole-genome scanning using eight common inbred strains (CIS) of laboratory mice, including A/J, C57BL/6J, CBA/J, DBA/2J, SM/J, SWR/J, NC/Nga, and 129/SvJ, and eight wild-derived inbred strains (WIS), BGL2/Ms, CAST/Ei, JF1/Ms, MSM/Ms, NJL/Ms, PGN2/Ms, SK/CamEi, and SWN/Ms. We selected and located 1226 informative loci at 1.2-cM average intervals on all of the chromosomes of the 16 strains and compared the polymorphisms of the eight CIS with those from the eight WIS as subspecies representatives. More than 50% of the loci can be identified as WIS (therefore, subspecies-specific) alleles in the CIS genomes. We also discovered that the CIS chromosomes form a mosaic structure with an average ratio of domesticus to non-domesticus alleles of 3:1. Furthermore, the domesticus alleles were present much more frequently on the CIS chromosome X than on their autosomes, suggesting that successive backcrossing of non-domesticus stocks to domesticus stocks had been undergone at the beginning of CIS history.

Estimation of linkage disequilibrium levels and haplotype block structure in Chinese Simmental and Wagyu beef cattle using high-density genotypes

Telomere length and telomerase activity directly affect the replicative capacity of primary human cells. Some have suggested that telomere length influences organismal lifespan. We compared telomere length distributions in a number of inbred and outbred established mouse strains with those of strains recently derived from wild mice. Telomere length was considerably shorter in wild-derived strains than in the established strains. We found no correlation of telomere length with lifespan, even among closely related inbred mouse strains. Thus, while telomere length plays a role in cellular lifespan in cultured human cells, it is not a major factor in determining organismal lifespan.

BackgroundGenetic isolates such as the Ashkenazi Jews (AJ) potentially offer advantages in mapping novel loci in whole genome disease association studies. To analyze patterns of genetic variation in AJ, genotypes of 101 healthy individuals were determined using the Affymetrix EAv3 500 K SNP array and compared to 60 CEPH-derived HapMap (CEU) individuals. 435,632 SNPs overlapped and met annotation criteria in the two groups.ResultsA small but significant global difference in allele frequencies between AJ and CEU was demonstrated by a mean FST of 0.009 (P < 0.001); large regions that differed were found on chromosomes 2 and 6. Haplotype blocks inferred from pairwise linkage disequilibrium (LD) statistics (Haploview) as well as by expectation-maximization haplotype phase inference (HAP) showed a greater number of haplotype blocks in AJ compared to CEU by Haploview (50,397 vs. 44,169) or by HAP (59,269 vs. 54,457). Average haplotype blocks were smaller in AJ compared to CEU (e.g., 36.8 kb vs. 40.5 kb HAP). Analysis of global patterns of local LD decay for closely-spaced SNPs in CEU demonstrated more LD, while for SNPs further apart, LD was slightly greater in the AJ. A likelihood ratio approach showed that runs of homozygous SNPs were approximately 20% longer in AJ. A principal components analysis was sufficient to completely resolve the CEU from the AJ.ConclusionLD in the AJ versus was lower than expected by some measures and higher by others. Any putative advantage in whole genome association mapping using the AJ population will be highly dependent on regional LD structure.

Genetic variations in the wild-derived inbred mouse strains are more diverse than thatof classical laboratory inbred mouse strains, including C57BL/6J (B6). The sleep/wake andmonoamine properties of six wild-derived inbred mouse strains (PGN2, NJL, BLG2, KJR, MSM,HMI) were characterized and compared with those of B6 mice. All examined mice werenocturnal and had a polyphasic sleep pattern with a “main sleep period” identified duringthe light period. However, there were three sleep/wake phenotypic differences between thewild-derived mouse strains and B6 strain. First, the amount of sleep during the dark phasewas comparable with that of B6 mice. However, the amount of sleep during the light phasewas more varied among strains, in particular, NJL and HMI had significantly less sleepcompared with that of B6 mice. Second, PGN2, NJL, BLG2, and KJR mice showed a “highlyawake period” (in which the hourly total sleep time was <10%) immediately after theonset of the dark period, which was not seen in B6 mice. Third, relative to that of B6mice, PGN2 and KJR mice showed longer duration of wakefulness episodes during the 12-hdark phase. Differences in whole brain noradrenaline, dopamine, and 5-hydroxy-tryptaminecontents between the wild-derived mouse strains and B6 strain were also found. Theseidentified phenotypes might be potentially under strong genetic control. Hence,wild-derived inbred mice could be useful for identifying the genetic factors underlyingthe regulation of sleep and wakefulness.

We compared four inbred mouse strains in their physical performance, measured as a maximal treadmill running time, characteristics of soleus muscle, anatomic character, and growth. The strains used were Mus musculus domesticus [C57BL/6 (B6) and BALB/c], Mus musculus molossinus (MSM/Ms), and Mus spretus. Maximal running time was significantly different among these four mouse strains. Running time until exhaustion was highest in MSM/Ms and lowest in M. spretus. Maximal times for the laboratory mouse strains were nearly identical. Soleus muscle fiber type and cross-sectional area also differed significantly among the species. In particular, M. spretus was significantly different from the other inbred mouse strains. Growth in the wild-derived inbred mice appeared to be complete earlier than in the laboratory mice, and the body size of the wild strains was about half that of the laboratory strains. From these results, we propose that wild-derived inbred mouse strains are useful models for enhancing phenotypic variation in physical performance and adaptability.

Fragile foal syndrome (FFS) is a disease caused by a recessive lethal missense mutation in the PLOD1 gene located on ECA2. Despite its harmful effect, a relatively high frequency of FFS carriers was observed in Warmblood breeds spanning from 7.4% in a random sample of Swedish Warmblood breed to 17% in the Hanoverian and Danish Warmblood, indicating potential heterozygous advantage. Balancing selection can be further studied based on haplotype blocks and via detection of heterozygosity‐rich region (ROHet) around the target of selection. In this study we evaluated the presence of haplotype blocks and ROHet on ECA2 in 380 Swedish Warmblood horses. We compared the results of ROHet with the rest of the genome. On average, 11.7 heterozygosity rich regions were identified per horse on ECA2, with no significant difference in numbers and length compared to what was found in other chromosomes. A unique haplotype block containing 28 markers was found in the FFS haplotype, while there were several haplotype blocks in the non‐carrier haplotype. This unique haplotype block mostly spanned the region upstream of the PLOD1 gene and included the MFN2 gene. The presence of this extended haplotype, shared by multiple individuals and including both the FFS variant and the MFN2 gene, suggests that this region may be under selection. While we did not find a clear heterozygosity‐rich region around the FFS variant, the extended haplotype may reflect either a signature of balancing selection or linkage disequilibrium with a positively selected variant in MFN2, PLOD1, or nearby loci.

The triangular correlation heatmap aiming to visualize the linkage disequilibrium (LD) pattern and haplotype block structure of SNPs is ubiquitous component of population-based genetic studies. However, current tools suffered from the problem of time and memory consuming. Here, we developed LDBlockShow, an open source software, for visualizing LD and haplotype blocks from variant call format files. It is time and memory saving. In a test dataset with 100 SNPs from 60000 subjects, it was at least 10.60 times faster and used only 0.03-13.33% of physical memory as compared with other tools. In addition, it could generate figures that simultaneously display additional statistical context (e.g. association P-values) and genomic region annotations. It can also compress the SVG files with a large number of SNPs and support subgroup analysis. This fast and convenient tool will facilitate the visualization of LD and haplotype blocks for geneticists.

Now that sequencing of the mouse genome has been completed, the function of each gene remains to be elucidated through phenotypic analysis. The "genetic background" (in which each gene functions) is defined as the genotype of all other related genes that may interact with the gene of interest, and therefore potentially influences the specific phenotype. To understand the nature and importance of genetic background on phenotypic expression of specific genes, it is necessary to know the origin and evolutionary history of the laboratory mouse genome. Molecular analysis has indicated that the fancy mice of Japan and Europe contributed significantly to the origin of today's laboratory mice. The genetic background of present-day laboratory mice varies by mouse strain, but is mainly derived from the European domesticus subspecies group and to a lesser degree from Asian mice, probably Japanese fancy mice, which belong to the musculus subspecies group. Inbred laboratory mouse strains are genetically uniform due to extensive inbreeding, and they have greatly contributed to the genetic analysis of many Mendelian traits. Meanwhile, for a variety of practical reasons, many transgenic and targeted mutant mice have been created in mice of mixed genetic backgrounds to elucidate the function of the genes, although efforts have been made to create inbred transgenic mice and targeted mutant mice with coisogenic embryonic stem cell lines. Inbred mouse strains have provided uniform genetic background for accurate evaluation of specific genes phenotypes, thus eliminating the phenotypic variations caused by mixed genetic backgrounds. However, the process of inbreeding and selection of various inbred strain characteristics has resulted in inadvertent selection of other undesirable genetic characteristics and mutations that may influence the genotype and preclude effective phenotypic analysis. Because many of the common inbred mouse stains have been established from relatively small gene pools, common inbred strains have limitations in their genetic polymorphisms and phenotypic variations. Wild-derived mouse strains can complement deficiencies of common inbred mouse strains, providing novel allelic variants and phenotypes. Although wild-derived strains are not as tame as the common laboratory strains, their genetic characteristics are attractive for the future study of gene function.

Genome-wide association is the most widely used methodology for mapping of disease genes for complex traits. Among many limitations of this approach is the requirement for investigating large disease populations and use of sizeable number of genetic markers, which necessitates correction for multiple comparisons and is a reason for its reduced power for signal detection. These limitations have prompted alternative approaches to identify novel disease loci. One such approach is the admixture mapping, which provides a more powered alternative for gene discovery in admixed populations for a trait that is differentially distributed in the ancestral populations. One such example is the African American population, which is at increased risk for chronic and end-stage renal disease, HIV-associated nephropathy, and focal segmental glomerulosclerosis, which are believed to originate from African ancestry. See Article by Shendre et al In the classical sense, admixture is the result of breeding between individuals from ≥2 isolated populations. The consequence of the gene flow is temporary generation of long haplotype blocks that contains genetic variants of one or another population. The approach in admixture mapping is based on the assumption that in the admixed populations the disease-causing alleles will be more frequent on chromosome segments derived from the ancestral population with the higher prevalence for the disease. The admixture mapping associates ancestry of haplotype block with a specific trait using differences in allele frequencies of the ancestral population. The haplotype blocks become shorter with increasing age of the admixed population through higher recombination events during each meiosis. Shorter admixed linkage disequilibrium (LD) blocks require a higher density of markers to differentiate chromosome ancestry transition. In admixed populations such as African Americans or Hispanics, the gene flow has originated within the last several hundred years, resulting in linked alleles that show extended LD relative to the ancestral populations. …

BackgroundHigh-density single-nucleotide polymorphism (SNP) genotyping arrays are a powerful tool for genome-wide association studies and can give valuable insight into patterns of population structure and linkage disequilibrium (LD). In this study we used the Brassica 60kSNP Illumina consortium genotyping array to assess the influence of selection and breeding for important agronomic traits on LD and haplotype structure in a diverse panel of 203 Chinese semi-winter rapeseed (Brassica napus) breeding lines.ResultsPopulation structure and principal coordinate analysis, using a subset of the SNPs, revealed diversification into three subpopulations and one mixed population, reflecting targeted introgressions from external gene pools during breeding. Pairwise LD analysis within the A- and C-subgenomes of allopolyploid B. napus revealed that mean LD, at a threshold of r2 = 0.1, decayed on average around ten times more rapidly in the A-subgenome (0.25-0.30 Mb) than in the C-subgenome (2.00-2.50 Mb). A total of 3,097 conserved haplotype blocks were detected over a total length of 182.49 Mb (15.17% of the genome). The mean size of haplotype blocks was considerably longer in the C-subgenome (102.85 Kb) than in the A-subgenome (33.51 Kb), and extremely large conserved haplotype blocks were found on a number of C-genome chromosomes. Comparative sequence analysis revealed conserved blocks containing homoloeogous quantitative trait loci (QTL) for seed erucic acid and glucosinolate content, two key seed quality traits under strong agronomic selection. Interestingly, C-subgenome QTL were associated with considerably greater conservation of LD than their corresponding A-subgenome homoeologues.ConclusionsThe data we present in this paper provide evidence for strong selection of large chromosome regions associated with important rapeseed seed quality traits conferred by C-subgenome QTL. This implies that an increase in genetic diversity and recombination within the C-genome is particularly important for breeding. The resolution of genome-wide association studies is also expected to vary greatly across different genome regions.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-1170) contains supplementary material, which is available to authorized users.