Low-coverage whole-genome sequencing facilitates accurate and cost-effective haplotype reconstruction in complex mouse crosses

The search for the underlying genetic contributions to complex traits and diseases relies on accurate genetic data from populations of interest. Outbred populations, like the Diversity Outbred (DO), are commonly genotyped using commercial SNP arrays, such as the Giga Mouse Universal Genotyping Array (GigaMUGA). However, array genotypes are expensive to collect, subject to significant ascertainment bias, and too sparse to capture the genetic structure of highly recombined mouse crosses. We investigated the efficacy of sequencing-based genotyping by comparing genotyping results between the GigaMUGA, double-digest restriction-site associated DNA sequencing (ddRADseq), and low-coverage whole-genome sequencing (lcWGS). We aligned reads at [~]1X coverage and imputed segregating SNPs from the eight DO founder strains onto 48 DO genomes and reconstructed their haplotypes using R/qtl2. Haplotype reconstructions derived from all three methods were highly concordant. However, lcWGS more faithfully recapitulated crossover counts and identified more small (< 1 Mb) haplotype blocks at as low as 0.1X coverage. Over 90% of local expression quantitative trait loci identified in a set of 183 DO-derived embryoid bodies using the GigaMUGA were recalled by lcWGS at coverages as low as 0.1X. We recommend that lcWGS be adopted as the primary method of genotyping complex crosses, and cell-based resources derived from them because they are as accurate as array-based reconstructions, robust to ultra-low sequencing depths, may more accurately model haplotypes of the mouse genome that are difficult to resolve with dense reference data, and cost-effective.