QTL mapping in outbred populations: successes and challenges.

Quantitative trait locus (QTL) mapping in animal populations has been a successful strategy for identifying genomic regions that play a role in complex diseases and traits. When conducted in an F2 intercross or backcross population, the resulting QTL is frequently large, often encompassing 30 Mb or more and containing hundreds of genes. To narrow the locus and identify candidate genes, additional strategies are needed. Congenic strains have proven useful but work less well when there are multiple tightly linked loci, frequently resulting in loss of phenotype. As an alternative, we discuss the use of highly recombinant outbred models for directly fine-mapping QTL to only a few megabases. We discuss the use of several currently available models such as the advanced intercross (AI), heterogeneous stocks (HS), the diversity outbred (DO), and commercially available outbred stocks (CO). Once a QTL has been fine-mapped, founder sequence and expression QTL mapping can be used to identify candidate genes. In this regard, the large number of alleles found in outbred stocks can be leveraged to identify causative genes and variants. We end this review by discussing some important statistical considerations when analyzing outbred populations. Fine-resolution mapping in outbred models, coupled with full genome sequence, has already led to the identification of several underlying causative genes for many complex traits and diseases. These resources will likely lead to additional successes in the coming years.

Protease inhibitor 15, a candidate gene for abdominal aortic internal elastic lamina ruptures in the rat.

The inbred Brown Norway (BN) rat develops spontaneous ruptures of the internal elastic lamina (RIEL) of the abdominal aorta (AA) and iliac arteries. Prior studies with crosses of the BN/Orl RJ (susceptible) and LOU/M (resistant) showed the presence of a significant QTL on chromosome 5 and the production of congenic rats proved the involvement of this locus. In this study, we further dissected the above-mentioned QTL by creating a new panel of LOU.BN(chr5) congenic and subcongenic lines and reduced the locus to 5.2 Mb. Then we studied 1,002 heterogeneous stock (HS) rats, whose phenotyping revealed a low prevalence and high variability for RIEL. High-resolution mapping in the HS panel detected the major locus on chromosome 5 (log P > 35) and refined it to 1.4 Mb. Subsequently, RNA-seq analysis on AA of BN, congenics, and LOU revealed expression differences for only protease inhibitor 15 (Pi15) gene and a putative long intergenic noncoding RNA (lincRNA) within the linkage region. The high abundance of lincRNA with respect to reduced Pi15 expression, in conjunction with exertion of longitudinal strain, may be related to RIEL, indicating the potential importance of proteases in biological processes related to defective aortic internal elastic lamina structure. Similar mechanisms may be involved in aneurysm initiation in the human AA.

A Cost-effective, High-throughput, Highly Accurate Genotyping Method for Outbred Populations.

Affordable sequencing and genotyping methods are essential for large scale genome-wide association studies. While genotyping microarrays and reference panels for imputation are available for human subjects, non-human model systems often lack such options. Our lab previously demonstrated an efficient and cost-effective method to genotype heterogeneous stock rats using double-digest genotyping-by-sequencing. However, low-coverage whole-genome sequencing offers an alternative method that has several advantages. Here, we describe a cost-effective, high-throughput, high-accuracy genotyping method for N/NIH heterogeneous stock rats that can use a combination of sequencing data previously generated by double-digest genotyping-by-sequencing and more recently generated by low-coverage whole-genome-sequencing data. Using double-digest genotyping-by-sequencing data from 5,745 heterogeneous stock rats (mean 0.21x coverage) and low-coverage whole-genome-sequencing data from 8,760 heterogeneous stock rats (mean 0.27x coverage), we can impute 7.32 million bi-allelic single-nucleotide polymorphisms with a concordance rate >99.76% compared to high-coverage (mean 33.26x coverage) whole-genome sequencing data for a subset of the same individuals. Our results demonstrate the feasibility of using sequencing data from double-digest genotyping-by-sequencing or low-coverage whole-genome-sequencing for accurate genotyping, and demonstrate techniques that may also be useful for other genetic studies in non-human subjects.

Genome-wide association study of delay discounting in Heterogeneous Stock rats.

Delay discounting refers to the behavioral tendency to devalue rewards as a function of their delay in receipt. Heightened delay discounting has been associated with substance use disorders and multiple co-occurring psychopathologies. Human and animal genetic studies have established that delay discounting is heritable, but only a few associated genes have been identified. We aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogeneous Stock (HS) rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female HS rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at various delays. Preference switch points were calculated and both exponential and hyperbolic functions were fitted to these indifference points. Area under the curve (AUC) and the discounting parameter k of both functions were used as delay discounting measures. GWAS for AUC, exponential k, and one indifference point identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family, was the sole gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression might be responsible for the association with behavior. Adgrl3, which encodes a latrophilin subfamily G-protein coupled receptor, was the sole gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

Y and Mitochondrial Chromosomes in the Heterogeneous Stock Rat Population.

Genome-wide association studies typically evaluate the autosomes and sometimes the X Chromosome, but seldom consider the Y or mitochondrial Chromosomes. We genotyped the Y and mitochondrial Chromosomes in heterogeneous stock rats (Rattus norvegicus), an outbred population created from eight inbred strains. We identified 8 distinct Y and 4 distinct mitochondrial Chromosomes among the 8 founders. However, only two types of each nonrecombinant chromosome were observed in our modern heterogeneous stock rat population (generations 81-97). Despite the relatively large sample size, there were virtually no significant associations for behavioral, physiological, metabolome, or microbiome traits after correcting for multiple comparisons. However, both Y and mitochondrial Chromosomes were strongly associated with expression of a few genes located on those chromosomes, which provided a positive control. Our results suggest that within modern heterogeneous stock rats there are no Y and mitochondrial Chromosomes differences that strongly influence behavioral or physiological traits. These results do not address other ancestral Y and mitochondrial Chromosomes that do not appear in modern heterogeneous stock rats, nor do they address effects that may exist in other rat populations, or in other species.

Association of tameness and sociability but no sign of domestication syndrome in mice selectively bred for active tameness.

Domesticated animals have been developed by selecting desirable traits following the initial unconscious selection stage, and now exhibit phenotypes desired by humans. Tameness is a common behavioural trait found in all domesticated animals. At the same time, these domesticated animals exhibit a variety of morphological, behavioural, and physiological traits that differ from their wild counterparts of their ancestral species. These traits are collectively referred to as domestication syndrome. However, whether this phenomenon exists is debatable. Previously, selective breeding has been used to enhance active tameness, a motivation to interact with humans, in wild heterogeneous stock mice derived from eight wild inbred strains. In the current study, we used tame mice to study how selective breeding for active tameness affects behavioural and morphological traits. A series of behavioural and morphological analyses on mice showed an increased preference for social stimuli and a longer duration of engagement in non-aggressive behaviour. However, no differences were observed in exploratory or anxiety-related behaviours. Similarly, selection for tameness did not affect ultrasonic vocalisations in mice, and no changes were observed in known morphological traits associated with domestication syndrome. These results suggest that there may be a link between active tameness and sociability and provide insights into the relationship between tameness and other behaviours in the context of domestication.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Radiation cataract in Heterogeneous Stock mice after γ-ray or HZE ion exposure.

Health effects of space radiation are a serious concern for astronauts on long-duration missions. The lens of the eye is one of the most radiosensitive tissues in the body and, therefore, ocular health risks for astronauts is a significant concern. Studies in humans and animals indicate that ionizing radiation exposure to the eye produces characteristic lens changes, termed "radiation cataract," that can affect visual function. Animal models of radiation cataractogenesis have previously utilized inbred mouse or rat strains. These studies were essential for determining morphological changes and dose-response relationships between radiation exposure and cataract. However, the relevance of these studies to human radiosensitivity is limited by the narrow phenotypic range of genetically homogeneous animal models. To model radiation cataract in genetically diverse populations, longitudinal cataract phenotyping was nested within a lifetime carcinogenesis study in male and female heterogeneous stock (HS/Npt) mice exposed to 0.4 Gy HZE ions (n = 609) or 3.0 Gy γ-rays (n = 602) and in unirradiated controls (n = 603). Cataractous change was quantified in each eye for up to 2 years using Merriam-Focht grading criteria by dilated slit lamp examination. Virtual Optomotry™ measurement of visual acuity and contrast sensitivity was utilized to assess visual function in a subgroup of mice. Prevalence and severity of posterior lens opacifications were 2.6-fold higher in HZE ion and 2.3-fold higher in γ-ray irradiated mice compared to unirradiated controls. Male mice were at greater risk for spontaneous and radiation associated cataracts. Risk for cataractogenesis was associated with family structure, demonstrating that HS/Npt mice are well-suited to evaluate genetic determinants of ocular radiosensitivity. Last, mice were extensively evaluated for cataract and tumor formation, which revealed an overlap between individual susceptibility to both cancer and cataract.

An exponential increase in QTL detection with an increased sample size.

Power analyses are often used to determine the number of animals required for a genome-wide association study (GWAS). These analyses are typically intended to estimate the sample size needed for at least 1 locus to exceed a genome-wide significance threshold. A related question that is less commonly considered is the number of significant loci that will be discovered with a given sample size. We used simulations based on a real data set that consisted of 3,173 male and female adult N/NIH heterogeneous stock rats to explore the relationship between sample size and the number of significant loci discovered. Our simulations examined the number of loci identified in subsamples of the full data set. The subsampling analysis was conducted for 4 traits with low (0.15 ± 0.03), medium (0.31 ± 0.03 and 0.36 ± 0.03), and high (0.46 ± 0.03) SNP-based heritabilities. For each trait, we subsampled the data 100 times at different sample sizes (500, 1,000, 1,500, 2,000, and 2,500). We observed an exponential increase in the number of significant loci with larger sample sizes. Our results are consistent with similar observations in human GWAS and imply that future rodent GWAS should use sample sizes that are significantly larger than those needed to obtain a single significant result.

Y and Mitochondrial Chromosomes in the Heterogeneous Stock Rat Population.

Genome-wide association studies typically evaluate the autosomes and sometimes the X Chromosome, but seldom consider the Y or mitochondrial Chromosomes. We genotyped the Y and mitochondrial chromosomes in heterogeneous stock rats (Rattus norvegicus), which were created in 1984 by intercrossing eight inbred strains and have subsequently been maintained as an outbred population for 100 generations. As the Y and mitochondrial Chromosomes do not recombine, we determined which founder had contributed these chromosomes for each rat, and then performed association analysis for all complex traits (n=12,055; intersection of 12,116 phenotyped and 15,042 haplotyped rats). We found the eight founders had 8 distinct Y and 4 distinct mitochondrial Chromosomes, however only two of each were observed in our modern heterogeneous stock rat population (Generations 81-97). Despite the unusually large sample size, the p-value distribution did not deviate from expectations; there were no significant associations for behavioral, physiological, metabolome, or microbiome traits after correcting for multiple comparisons. However, both Y and mitochondrial Chromosomes were strongly associated with expression of a few genes located on those chromosomes, which provided a positive control. Our results suggest that within modern heterogeneous stock rats there are no Y and mitochondrial Chromosomes differences that strongly influence behavioral or physiological traits. These results do not address other ancestral Y and mitochondrial Chromosomes that do not appear in modern heterogeneous stock rats, nor do they address effects that may exist in other rat populations, or in other species.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks.

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Abstinence from Escalation of Cocaine Intake Changes the microRNA Landscape in the Cortico-Accumbal Pathway.

Cocaine administration alters the microRNA (miRNA) landscape in the cortico-accumbal pathway. These changes in miRNA can play a major role in the posttranscriptional regulation of gene expression during withdrawal. This study aimed to investigate the changes in microRNA expression in the cortico-accumbal pathway during acute withdrawal and protracted abstinence following escalated cocaine intake. Small RNA sequencing (sRNA-seq) was used to profile miRNA transcriptomic changes in the cortico-accumbal pathway [infralimbic- and prelimbic-prefrontal cortex (IL and PL) and nucleus accumbens (NAc)] of rats with extended access to cocaine self-administration followed by an 18-h withdrawal or a 4-week abstinence. An 18-h withdrawal led to differential expression (fold-change > 1.5 and p < 0.05) of 21 miRNAs in the IL, 18 miRNAs in the PL, and two miRNAs in the NAc. The mRNAs potentially targeted by these miRNAs were enriched in the following pathways: gap junctions, neurotrophin signaling, MAPK signaling, and cocaine addiction. Moreover, a 4-week abstinence led to differential expression (fold-change > 1.5 and p < 0.05) of 23 miRNAs in the IL, seven in the PL, and five miRNAs in the NAc. The mRNAs potentially targeted by these miRNAs were enriched in pathways including gap junctions, cocaine addiction, MAPK signaling, glutamatergic synapse, morphine addiction, and amphetamine addiction. Additionally, the expression levels of several miRNAs differentially expressed in either the IL or the NAc were significantly correlated with addiction behaviors. Our findings highlight the impact of acute and protracted abstinence from escalated cocaine intake on miRNA expression in the cortico-accumbal pathway, a key circuit in addiction, and suggest developing novel biomarkers and therapeutic approaches to prevent relapse by targeting abstinence-associated miRNAs and their regulated mRNAs.

Genome-wide association study finds multiple loci associated with intraocular pressure in HS rats.

Elevated intraocular pressure (IOP) is influenced by environmental and genetic factors. Increased IOP is a major risk factor for most types of glaucoma, including primary open angle glaucoma (POAG). Investigating the genetic basis of IOP may lead to a better understanding of the molecular mechanisms of POAG. The goal of this study was to identify genetic loci involved in regulating IOP using outbred heterogeneous stock (HS) rats. HS rats are a multigenerational outbred population derived from eight inbred strains that have been fully sequenced. This population is ideal for a genome-wide association study (GWAS) owing to the accumulated recombinations among well-defined haplotypes, the relatively high allele frequencies, the accessibility to a large collection of tissue samples, and the large allelic effect size compared to human studies. Both male and female HS rats (N = 1,812) were used in the study. Genotyping-by-sequencing was used to obtain ∼3.5 million single nucleotide polymorphisms (SNP) from each individual. SNP heritability for IOP in HS rats was 0.32, which agrees with other studies. We performed a GWAS for the IOP phenotype using a linear mixed model and used permutation to determine a genome-wide significance threshold. We identified three genome-wide significant loci for IOP on chromosomes 1, 5, and 16. Next, we sequenced the mRNA of 51 whole eye samples to find cis-eQTLs to aid in identification of candidate genes. We report 5 candidate genes within those loci: Tyr, Ctsc, Plekhf2, Ndufaf6 and Angpt2. Tyr, Ndufaf6 and Angpt2 genes have been previously implicated by human GWAS of IOP-related conditions. Ctsc and Plekhf2 genes represent novel findings that may provide new insight into the molecular basis of IOP. This study highlights the efficacy of HS rats for investigating the genetics of elevated IOP and identifying potential candidate genes for future functional testing.

The utility of a closed breeding colony of Peromyscus leucopus for dissecting complex traits.

Deermice of the genus Peromyscus are well suited for addressing several questions of biologist interest, including the genetic bases of longevity, behavior, physiology, adaptation, and their ability to serve as disease vectors. Here, we explore a diversity outbred approach for dissecting complex traits in Peromyscus leucopus, a nontraditional genetic model system. We take advantage of a closed colony of deer-mice founded from 38 individuals and subsequently maintained for ∼40-60 generations. From 405 low-pass short-read sequenced deermice we accurate impute genotypes at 16 million single nucleotide polymorphisms. Conditional on observed genotypes simulations were conducted in which three different sized quantitative trait loci contribute to a complex trait under three different genetic models. Using a stringent significance threshold power was modest, largely a function of the percent variation attributable to the simulated quantitative trait loci, with the underlying genetic model having only a subtle impact. We additionally simulated 2,000 pseudo-individuals, whose genotypes were consistent with those observed in the genotyped cohort and carried out additional power simulations. In experiments employing more than 1,000 mice power is high to detect quantitative trait loci contributing greater than 2.5% to a complex trait, with a localization ability of ∼100 kb. We finally carried out a Genome-Wide Association Study on two demonstration traits, bleeding time and body weight, and uncovered one significant region. Our work suggests that complex traits can be dissected in founders-unknown P. leucopus colony mice and similar colonies in other systems using easily obtained genotypes from low-pass sequencing.

Genome-Wide Association Study on Three Behaviors Tested in an Open Field in Heterogeneous Stock Rats Identifies Multiple Loci Implicated in Psychiatric Disorders.

Many personality traits are influenced by genetic factors. Rodents models provide an efficient system for analyzing genetic contribution to these traits. Using 1,246 adolescent heterogeneous stock (HS) male and female rats, we conducted a genome-wide association study (GWAS) of behaviors measured in an open field, including locomotion, novel object interaction, and social interaction. We identified 30 genome-wide significant quantitative trait loci (QTL). Using multiple criteria, including the presence of high impact genomic variants and co-localization of cis-eQTL, we identified 17 candidate genes (Adarb2, Ankrd26, Cacna1c, Cacng4, Clock, Ctu2, Cyp26b1, Dnah9, Gda, Grxcr1, Eva1a, Fam114a1, Kcnj9, Mlf2, Rab27b, Sec11a, and Ube2h) for these traits. Many of these genes have been implicated by human GWAS of various psychiatric or drug abuse related traits. In addition, there are other candidate genes that likely represent novel findings that can be the catalyst for future molecular and genetic insights into human psychiatric diseases. Together, these findings provide strong support for the use of the HS population to study psychiatric disorders.

A generic hidden Markov model for multiparent populations.

A common step in the analysis of multiparent populations (MPPs) is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of MPPs. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly available data on Diversity Outbred and Collaborative Cross mice.

Bisphenol F Exposure in Adolescent Heterogeneous Stock Rats Affects Growth and Adiposity.

Bisphenol F (BPF) is increasingly substituting bisphenol A in manufacturing polycarbonates and consumer products. The cardiometabolic effects of BPF in either humans or model organisms are not clear, and no studies to date have investigated the role of genetic background on susceptibility to BPF-induced cardiometabolic traits. The primary goal of this project was to determine if BPF exposure influences growth and adiposity in male N:NIH heterogeneous stock (HS) rats, a genetically heterogeneous population. Littermate pairs of male HS rats were randomly exposed to either vehicle (0.1% ethanol) or 1.125 µg/ml BPF in 0.1% ethanol for 5 weeks in drinking water starting at 3 weeks-of-age. Water consumption and body weight was measured weekly, body composition was determined using nuclear magnetic resonance, urine and feces were collected in metabolic cages, and blood and tissues were collected at the end of the study. BPF-exposed rats showed significantly increased body growth and abdominal adiposity, risk factors for cardiometabolic disease. Urine output was increased in BPF-exposed rats, driving a trend in increased creatinine clearance. We also report the first relationship between a bisphenol metabolizing enzyme and a bisphenol-induced phenotype. Preliminary heritability estimates of significant phenotypes suggest that BPF exposure may alter trait variation. These findings support BPF exposure as a cardiometabolic disease risk factor and indicate that the HS rat will be a useful model for dissecting gene by BPF interactions on metabolic health.

Combined change of behavioral traits for domestication and gene-networks in mice selectively bred for active tameness.

Tameness is a major element of animal domestication and involves two components: motivation to approach humans (active tameness) and reluctance to avoid humans (passive tameness). To understand the behavioral and genetic mechanisms of active tameness in mice, we had previously conducted selective breeding for long durations of contact and heading toward human hands in an active tameness test using a wild-derived heterogeneous stock. Although the study showed a significant increase in contacting and heading with the 12th generation of breeding, the effect on other behavioral indices related to tameness and change of gene expression levels underlying selective breeding was unclear. Here, we analyzed nine tameness-related traits at a later stage of selective breeding and analyzed how gene expression levels were changed by the selective breeding. We found that five traits, including contacting and heading, showed behavioral change in the selective groups comparing to the control through the generations. Furthermore, we conducted cluster analyses to evaluate the relationships among the nine traits and found that contacting and heading combined in an independent cluster in the selected groups, but not in the control groups. RNA-Seq of hippocampal tissue revealed differential expression of 136 genes between the selection and control groups, while the pathway analysis identified the networks associated with these genes. These results suggest that active tameness was hidden in the control groups but became apparent in the selected populations by selective breeding, potentially driven by changes in gene expression networks.

Adapting Genotyping-by-Sequencing and Variant Calling for Heterogeneous Stock Rats.

The heterogeneous stock (HS) is an outbred rat population derived from eight inbred rat strains. HS rats are ideally suited for genome wide association studies; however, only a few genotyping microarrays have ever been designed for rats and none of them are currently in production. To address the need for an efficient and cost effective method of genotyping HS rats, we have adapted genotype-by-sequencing (GBS) to obtain genotype information at large numbers of single nucleotide polymorphisms (SNPs). In this paper, we have outlined the laboratory and computational steps we took to optimize double digest genotype-by-sequencing (ddGBS) for use in rats. We evaluated multiple existing computational tools and explain the workflow we have used to call and impute over 3.7 million SNPs. We have also compared various rat genetic maps, which are necessary for imputation, including a recently developed map specific to the HS. Using our approach, we obtained concordance rates of 99% with data obtained using data from a genotyping array. The principles and computational pipeline that we describe could easily be adapted for use in other species for which reliable reference genome sets are available.