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.     

An investigation of genetic variation within a series of congenic strains of mice

2 congenic strains of mice, B6N.AKN-Ahk and D2N.B6N-Ahb, imported from the USA, were found to be either segregating or fixed for an incorrect allele at a number of biochemical loci. B6N.AKN-Ahk, supposedly congenic with C57BL/6N, had the wrong genotype at 6 out of 12 biochemical loci; D2N.B6N-Ahb, supposedly congenic with DBA/2N, was segregating at 3 out of 9 loci. There was genetic variation in mandible shape within the 2 strains but no abnormal coat colours were found and no hybrid vigour in breeding performance was detected. Analyses in the USA confirmed these results and showed that 2 other congenic strains, C3N.D2N-Ahd and AKN.B6J-Ahb, were also segregating at a number of loci. Some of the alleles found in the C3N.D2N-Ahd mice must be the result of a genetic contamination. The simplest explanation for this breakdown in the backcrossing programme is genetic contamination with other congenic strains or recombinant inbred lines under development in the same laboratory. These findings emphasize the importance of continual genetic monitoring of all genetic stocks at regular intervals and in particular during the development of congenic and recombinant lines.     

Genetic analysis of the diabetes-prone C57BLKS/J mouse strain reveals genetic contribution from multiple strains

The C57BLKS/J (BKS) inbred mouse strain is a widely used animal model of type 2 diabetes. In the presence of the diabetes (db) mutation, obese BKS-db mice develop severe diabetes. Genetic studies of diabetes-susceptibility in this strain are facilitated by the fact that BKS is a genetic composite between the diabetes-resistant C57BL/6J (B6) and susceptible DBA/2J (DBA) strains. On this basis, it has been hypothesized that diabetes-susceptibility in BKS is conferred by DBA-derived alleles. However, recent studies revealed non-B6/non-DBA genetic material in BKS. To identify the origin of this genetic component, we generated a genomic map of BKS using 537 microsatellite markers. Our results demonstrate that, in addition to B6 and DBA, BKS contains alleles from at least three other strains, including 129, C57BL/10 and an unidentified mouse strain. We also analyzed two congenic strains, B6-db and BKS-db, which are widely used for the genetic mapping of diabetes-susceptibility loci. We identified several donor-derived genomic regions introduced during the generation of these congenic strains. In summary, our study reveals novel aspects of the genetic fine-structure of BKS and related strains and facilitates the identification of diabetes-susceptibility loci in this mouse model.     

Segregation of a QTL cluster for home-cage activity using a new mapping method based on regression analysis of congenic mouse strains

Recent genetic studies have shown that genetic loci with significant effects in whole-genome quantitative trait loci (QTL) analyses were lost or weakened in congenic strains. Characterisation of the genetic basis of this attenuated QTL effect is important to our understanding of the genetic mechanisms of complex traits. We previously found that a consomic strain, B6-Chr6C^MSM, which carries chromosome 6 of a wild-derived strain MSM/Ms on the genetic background of C57BL/6J, exhibited lower home-cage activity than C57BL/6J. In the present study, we conducted a composite interval QTL analysis using the F2 mice derived from a cross between C57BL/6J and B6-Chr6C^MSM. We found one QTL peak that spans 17.6 Mbp of chromosome 6. A subconsomic strain that covers the entire QTL region also showed lower home-cage activity at the same level as the consomic strain. We developed 15 congenic strains, each of which carries a shorter MSM/Ms-derived chromosomal segment from the subconsomic strain. Given that the results of home-cage activity tests on the congenic strains cannot be explained by a simple single-gene model, we applied regression analysis to segregate the multiple genetic loci. The results revealed three loci (loci 1–3) that have the effect of reducing home-cage activity and one locus (locus 4) that increases activity. We also found that the combination of loci 3 and 4 cancels out the effects of the congenic strains, which indicates the existence of a genetic mechanism related to the loss of QTLs.     

A genome-wide set of congenic mouse strains derived from DBA/2J on a C57BL/6J background

In the analysis of complex traits, congenic strains are powerful tools because they allow characterization of a single locus in the absence of genetic variation throughout the remainder of the genome. Here, we report the construction and initial characterization of a genome-wide panel of congenic strains derived from the donor strain DBA/2J on the background strain C57BL/6J. For many strains, we have carried out high-density SNP genotyping to precisely map the congenic interval and to identify any contaminating regions. Certain strains exhibit striking variation in litter size and in the ratio of females to males. We illustrate the utility of the set by "Mendelizing" the complex trait of myocardial calcification. These 65 strains cover more than 95% of the autosomal genome and should facilitate the analysis of the many genetic trait differences that have been reported between these parental strains.     

The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits

Genome-wide association studies have identified hundreds of loci for type 2 diabetes, coronary artery disease and myocardial infarction, as well as for related traits such as body mass index, glucose and insulin levels, lipid levels, and blood pressure. These studies also have pointed to thousands of loci with promising but not yet compelling association evidence. To establish association at additional loci and to characterize the genome-wide significant loci by fine-mapping, we designed the "Metabochip," a custom genotyping array that assays nearly 200,000 SNP markers. Here, we describe the Metabochip and its component SNP sets, evaluate its performance in capturing variation across the allele-frequency spectrum, describe solutions to methodological challenges commonly encountered in its analysis, and evaluate its performance as a platform for genotype imputation. The metabochip achieves dramatic cost efficiencies compared to designing single-trait follow-up reagents, and provides the opportunity to compare results across a range of related traits. The metabochip and similar custom genotyping arrays offer a powerful and cost-effective approach to follow-up large-scale genotyping and sequencing studies and advance our understanding of the genetic basis of complex human diseases and traits.     

High-resolution melting analysis for SNP genotyping and mapping in tetraploid alfalfa (Medicago sativa L.)

Single nucleotide polymorphisms (SNPs) represent the most abundant type of genetic polymorphism in plant genomes. SNP markers are valuable tools for genetic analysis of complex traits of agronomic importance, linkage and association mapping, genome-wide selection, map-based cloning, and marker-assisted selection. Current challenges for SNP genotyping in polyploid outcrossing species include multiple alleles per loci and lack of high-throughput methods suitable for variant detection. In this study, we report on a high-resolution melting (HRM) analysis system for SNP genotyping and mapping in outcrossing tetraploid genotypes. The sensitivity and utility of this technology is demonstrated by identification of the parental genotypes and segregating progeny in six alfalfa populations based on unique melting curve profiles due to differences in allelic composition at one or multiple loci. HRM using a 384-well format is a fast, consistent, and efficient approach for SNP discovery and genotyping, useful in polyploid species with uncharacterized genomes. Possible applications of this method include variation discovery, analysis of candidate genes, genotyping for comparative and association mapping, and integration of genome-wide selection in breeding programs.     

Identification of Mom12 and Mom13, two novel modifier loci of ApcMin-mediated intestinal tumorigenesis

Colorectal cancer is a heterogeneous disease resulting from a combination of genetic and environmental factors. The C57BL/6J (B6) Apc^Min/+ mouse develops polyps throughout the gastrointestinal tract and has been a valuable model for understanding the genetic basis of intestinal tumorigenesis. Apc^Min/+ mice have been used to study known oncogenes and tumor suppressor genes on a controlled genetic background. These studies often utilize congenic knockout alleles, which can carry an unknown amount of residual donor DNA. The Apc^Min model has also been used to identify modifer loci, known as Modifier of Min (Mom) loci, which alter Apc^Min-mediated intestinal tumorigenesis. B6 mice carrying a knockout allele generated in WW6 embryonic stem cells were crossed to B6 Apc^Min/+ mice to determine the effect on polyp multiplicity. The newly generated colony developed significantly more intestinal polyps than Apc^Min/+ controls. Polyp multiplicity did not correlate with inheritance of the knockout allele, suggesting the presence of one or more modifier loci segregating in the colony. Genotyping of simple sequence length polymorphism (SSLP) markers revealed residual 129X1/SvJ genomic DNA within the congenic region of the parental knockout line. An analysis of polyp multiplicity data and SSLP genotyping indicated the presence of two Mom loci in the colony: (1) Mom12, a dominant modifier linked to the congenic region on chromosome 6 and (2) Mom13, which is unlinked to the congenic region and whose effect is masked by Mom12. The identification of Mom12 and Mom13 demonstrates the potential problems resulting from residual heterozygosity present in congenic lines.Key words: adenomatous polyposis coli, modifier of min, congenic mice, caveolin-1, cancer susceptibility     

Production of congenic mouse strains carrying genomic intervals containing SLE-susceptibility genes derived from the SLE-prone NZM2410 strain

Systemic lupus erythematosus is inherited as a complex polygenic trait. Four genomic intervals containing major SLE-susceptibility loci were previously identified by interval mapping in the NZM2410 mouse model. In this paper, we utilized a marker-assisted selection protocol to produce four congenic mouse strains, each carrying an NZM2410-derived SLE-susceptibility interval on a C57BL/6-resistant background. Each strain carries only one susceptibility allele derived from this polygenic model and consequently can be used to characterize the specific component phenotypes contributed by individual SLE-susceptibility genes. We illustrate the efficacy of this approach with phenotypic data for one of our congenic strains, B6.NZMH2 ^z . Our results indicate that this single genomic interval from Chromosome (Chr) 17 of NZM2410 can mediate increased levels of IgG autoantibodies specific for chromatin and that, similar to results obtained in our original genetic cross, B6.NZMH2 ^z/b heterozygotes are more prone than B6.NZMH2 ^z homozygotes to the development of humoral autoimmunity to nuclear antigens. These results illustrate the feasibility of using congenic strains to dissect the complex pathogenic mechanisms that mediate polygenic SLE. These congenic strains will be valuable tools in the genetic analysis of SLE susceptibility. In future studies, these congenic strains will be interbred to produce bi- and tri-congenic strains in order to assess the role of genetic interactions in the expression of specific components of SLE pathogenesis. They will also be instrumental to the positional cloning and identification of the genes responsible for SLE susceptibility, via the production of congenic recombinants. Received: 1 September 1995 / Accepted: 20 December 1995     

Identification of Mom12 and Mom13, two novel modifier loci of ApcMin-mediated intestinal tumorigenesis

Colorectal cancer is a heterogeneous disease resulting from a combination of genetic and environmental factors. The C57BL/6J (B6) Apc^Min/+ mouse develops polyps throughout the gastrointestinal tract and has been a valuable model for understanding the genetic basis of intestinal tumorigenesis. Apc^Min/+ mice have been used to study known oncogenes and tumor suppressor genes on a controlled genetic background. These studies often utilize congenic knockout alleles, which can carry an unknown amount of residual donor DNA. The Apc^Min model has also been used to identify modifer loci, known as Modifier of Min (Mom) loci, which alter Apc^Min-mediated intestinal tumorigenesis. B6 mice carrying a knockout allele generated in WW6 embryonic stem cells were crossed to B6 Apc^Min/+ mice to determine the effect on polyp multiplicity. The newly generated colony developed significantly more intestinal polyps than Apc^Min/+ controls. Polyp multiplicity did not correlate with inheritance of the knockout allele, suggesting the presence of one or more modifier loci segregating in the colony. Genotyping of simple sequence length polymorphism (SSLP) markers revealed residual 129X1/SvJ genomic DNA within the congenic region of the parental knockout line. An analysis of polyp multiplicity data and SSLP genotyping indicated the presence of two Mom loci in the colony: 1) Mom12, a dominant modifier linked to the congenic region on chromosome 6, and 2) Mom13, which is unlinked to the congenic region and whose effect is masked by Mom12. The identification of Mom12 and Mom13 demonstrates the potential problems resulting from residual heterozygosity present in congenic lines.     

Immunodominance in the T cell response to multiple non-H-2 histocompatibility antigens. IV. Partial tissue distribution and mapping of immunodominant antigens

The immunization of C57BL/6 responder mice with spleen cells from H-2-matched BALB.B donors, which differ by multiple non-H-2 histocompatibility (H) antigens, results in the generation of cytotoxic T lymphocytes (CTL) that are specific for only a limited number of immunodominant antigens. Previous analysis of the genes encoding these dominant antigens has not mapped these genes to any of the non-H-2 H loci defined by congenic strains. It would have been expected that the histogenetic techniques employed for congenic strain selection would have preferentially identified the "strongest" H antigens. Therefore, we have investigated the possibility that immunodominant antigens do not belong to the class of non-H-2 H antigens encoded by genes mapping to H loci defined and mapped by congenic strains. The first experiments were aimed at identifying antigens that were expressed by independently derived inbred strains and were cross-reactive with the immunodominant cytotoxic T cell target (CTT-1) antigen of BALB.B. Strong cross-reaction with the C3H.SW (H-2b) strain was observed; the C3H gene encoding this antigen was mapped with BXH recombinant inbred strains. Contrary to the mapping of the CTT-1 gene to chromosome 1 in BALB.B, the C3H gene was shown to map to either chromosome 4 or chromosome 7. This result indicates that identical, or at least extensively cross-reactive, non-H-2 antigens may be encoded by genes mapping to independently segregating loci in different inbred strains. The tissue distribution of immunodominant antigens was approached by determining the reactivity of CTL specific for these antigens with either lymphoid-derived or fibroblast-derived targets. These CTL effectively lysed lymphoblast and lymphoid tumor targets but did not lyse an SV40-transformed fibroblast line that was shown to be efficiently lysed by CTL specific for non-H-2 H antigens defined by congenic strains. Therefore, it was concluded that immunodominant antigens detected by B6 anti-BALB.B CTL have a restricted tissue distribution in comparison to non-H-2 H antigens defined by congenic strains. The implications of these results for our understanding of the origin and heterogeneity of non-H-2 cell-surface antigen recognized by effector T cells are discussed.     

Single nucleotide polymorphism genotyping by mini-primer allele-specific amplification with universal reporter primers for identification of degraded DNA

Single nucleotide polymorphism (SNP) is informative for human identification, and much shorter regions are targeted in analysis of biallelic SNP compared with highly polymorphic short tandem repeat (STR). Therefore, SNP genotyping is expected to be more sensitive than STR genotyping of degraded human DNA. To achieve simple, economical, and sensitive SNP genotyping for identification of degraded human DNA, we developed 18 loci for a SNP genotyping technique based on the mini-primer allele-specific amplification (ASA) combined with universal reporter primers (URP). The URP/ASA-based genotyping consisted of two amplifications followed by detection using capillary electrophoresis. The sizes of the target genome fragments ranged from 40 to 67 bp in length. In the Japanese population, the frequencies of minor alleles of 18 SNPs ranged from 0.36 to 0.50, and these SNPs are informative for identification. The success rate of SNP genotyping was much higher than that of STR genotyping of artificially degraded DNA. Moreover, we applied this genotyping method to case samples and showed successful SNP genotyping of severely degraded DNA from a 4-year buffered formalin-fixed tissue sample for human identification.     

Single nucleotide polymorphism genotyping in polyploid wheat with the Illumina GoldenGate assay

Single nucleotide polymorphisms (SNPs) are indispensable in such applications as association mapping and construction of high-density genetic maps. These applications usually require genotyping of thousands of SNPs in a large number of individuals. Although a number of SNP genotyping assays are available, most of them are designed for SNP genotyping in diploid individuals. Here, we demonstrate that the Illumina GoldenGate assay could be used for SNP genotyping of homozygous tetraploid and hexaploid wheat lines. Genotyping reactions could be carried out directly on genomic DNA without the necessity of preliminary PCR amplification. A total of 53 tetraploid and 38 hexaploid homozygous wheat lines were genotyped at 96 SNP loci. The genotyping error rate estimated after removal of low-quality data was 0 and 1% for tetraploid and hexaploid wheat, respectively. Developed SNP genotyping assays were shown to be useful for genotyping wheat cultivars. This study demonstrated that the GoldenGate assay is a very efficient tool for high-throughput genotyping of polyploid wheat, opening new possibilities for the analysis of genetic variation in wheat and dissection of genetic basis of complex traits using association mapping approach. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A comprehensive SNP-based genetic analysis of inbred mouse strains

Dense genetic maps of mammalian genomes facilitate a variety of biological studies including the mapping of polygenic traits, positional cloning of monogenic traits, mapping of quantitative or qualitative trait loci, marker association, allelic imbalance, speed congenic construction, and evolutionary or phylogenetic comparison. In particular, single nucleotide polymorphisms (SNPs) have proved useful because of their abundance and compatibility with multiple high-throughput technology platforms. SNP genotyping is especially suited for the genetic analysis of model organisms such as the mouse because biallelic markers remain fully informative when used to characterize crosses between inbred strains. Here we report the mapping and genotyping of 673 SNPs (including 519 novel SNPs) in 55 of the most commonly used mouse strains. These data have allowed us to construct a phylogenetic tree that correlates and expands known genealogical relationships and clarifies the origin of strains previously having an uncertain ancestry. All 55 inbred strains are distinguishable genetically using this SNP panel. Our data reveal an uneven SNP distribution consistent with a mosaic pattern of inheritance and provide some insight into the changing dynamics of the physical architecture of the genome. Furthermore, these data represent a valuable resource for the selection of markers and the design of experiments that require the genetic distinction of any pair of mouse inbred strains such as the generation of congenic mice, positional cloning, and the mapping of quantitative or qualitative trait loci.The content of this publication does not necessarily reflect the view or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.     

Rapid and Efficient Zebrafish Genotyping Using PCR with High-resolution Melt Analysis

Zebrafish is a powerful vertebrate model system for studying development, modeling disease, and performing drug screening. Recently a variety of genetic tools have been introduced, including multiple strategies for inducing mutations and generating transgenic lines. However, large-scale screening is limited by traditional genotyping methods, which are time-consuming and labor-intensive. Here we describe a technique to analyze zebrafish genotypes by PCR combined with high-resolution melting analysis (HRMA). This approach is rapid, sensitive, and inexpensive, with lower risk of contamination artifacts. Genotyping by PCR with HRMA can be used for embryos or adult fish, including in high-throughput screening protocols.     

Development and Characterization of a High Density SNP Genotyping Assay for Cattle

The success of genome-wide association (GWA) studies for the detection of sequence variation affecting complex traits in human has spurred interest in the use of large-scale high-density single nucleotide polymorphism (SNP) genotyping for the identification of quantitative trait loci (QTL) and for marker-assisted selection in model and agricultural species. A cost-effective and efficient approach for the development of a custom genotyping assay interrogating 54,001 SNP loci to support GWA applications in cattle is described. A novel algorithm for achieving a compressed inter-marker interval distribution proved remarkably successful, with median interval of 37 kb and maximum predicted gap of <350 kb. The assay was tested on a panel of 576 animals from 21 cattle breeds and six outgroup species and revealed that from 39,765 to 46,492 SNP are polymorphic within individual breeds (average minor allele frequency (MAF) ranging from 0.24 to 0.27). The assay also identified 79 putative copy number variants in cattle. Utility for GWA was demonstrated by localizing known variation for coat color and the presence/absence of horns to their correct genomic locations. The combination of SNP selection and the novel spacing algorithm allows an efficient approach for the development of high-density genotyping platforms in species having full or even moderate quality draft sequence. Aspects of the approach can be exploited in species which lack an available genome sequence. The BovineSNP50 assay described here is commercially available from Illumina and provides a robust platform for mapping disease genes and QTL in cattle.     

Susceptibility to Seizure-Induced Excitotoxic Cell Death Is Regulated by an Epistatic Interaction between Chr 18 (Sicd1) and Chr 15 (Sicd2) Loci in Mice

Seizure-induced cell death is believed to be regulated by multiple genetic components in addition to numerous external factors. We previously defined quantitative trait loci that control susceptibility to seizure-induced cell death in FVB/NJ (susceptible) and C57BL/6J (resistant) mice. Two of these quantitative trait loci assigned to chromosomes 18 (Sicd1) and 15 (Sicd2), control seizure-induced cell death resistance. In this study, through the use of a series of novel congenic strains containing the Sicd1 and Sicd2 congenic strains and different combinations of the Sicd1 or Sicd2 sub region(s), respectively, we defined these genetic interactions. We generated a double congenic strain, which contains the two C57BL/6J differential segments from chromosome 18 and 15, to determine how these two segments interact with one another. Phenotypic comparison between FVB-like littermates and the double congenic FVB.B6-Sicd1/Sicd2 strain identified an additive effect with respect to resistance to seizure-induced excitotoxic cell death. It thus appears that C57BL/6J alleles located on chromosomes 18 and 15 interact epistatically in an additive manner to control the extent of seizure-induced excitotoxic cell death. Three interval-specific congenic lines were developed, in which either segments of C57BL/6J Chr 18 or C57BL/6J Chr 15 were introduced in the FVB/NJ genetic background, and progeny were treated with kainate and examined for the extent of seizure-induced cell death. All of the interval-specific congenic lines exhibited reduced cell death in both area CA3 and the dentate hilus, associated with the C57BL/6J phenotype. These experiments demonstrate functional interactions between Sicd1 and Sicd2 that improve resistance to seizure-induced excitotoxic cell death, validating the critical role played by gene-gene interactions in excitotoxic cell death.     

A large QTL for fear and anxiety mapped using an F2 cross can be dissected into multiple smaller QTLs

Using chromosome substitution strains (CSS), we previously identified a large quantitative trait locus (QTL) for conditioned fear (CF) on mouse chromosome 10. Here, we used an F₂ cross between CSS‐10 and C57BL/6J (B6) to localize that QTL to distal chromosome 10. That QTL accounted for all the difference between CSS‐10 and B6. We then produced congenic strains to fine‐map that interval. We identified two congenic strains that captured some or all the QTL. The larger congenic strain (Line 1: 122.387121–129.068 Mb; build 37) appeared to account for all the difference between CSS‐10 and B6. The smaller congenic strain (Line 2: 127.277–129.068 Mb) was intermediate between CSS‐10 and B6. We used haplotype mapping followed by quantitative polymerase chain reaction to identify one gene that was differentially expressed in both lines relative to B6 (Rnf41) and one that was differentially expressed between only Line 1 and B6 (Shmt2). These cis‐eQTLs may cause the behavioral QTLs; however, further studies are required to validate these candidate genes. More generally, our observation that a large QTL mapped using CSS and F₂ crosses can be dissected into multiple smaller QTLs shows a weaknesses of two‐stage approaches that seek to use coarse mapping to identify large regions followed by fine‐mapping. Indeed, additional dissection of these congenic strains might result in further subdivision of these QTL regions. Despite these limitations, we have successfully fine‐mapped two QTLs to small regions and identified putative candidate genes, showing that the congenic approach can be effective for fine‐mapping QTLs .     

Using Genotyping-By-Sequencing (GBS) for Genomic Discovery in Cultivated Oat

Advances in next-generation sequencing offer high-throughput and cost-effective genotyping alternatives, including genotyping-by-sequencing (GBS). Results have shown that this methodology is efficient for genotyping a variety of species, including those with complex genomes. To assess the utility of GBS in cultivated hexaploid oat (Avena sativa L.), seven bi-parental mapping populations and diverse inbred lines from breeding programs around the world were studied. We examined technical factors that influence GBS SNP calls, established a workflow that combines two bioinformatics pipelines for GBS SNP calling, and provided a nomenclature for oat GBS loci. The high-throughput GBS system enabled us to place 45,117 loci on an oat consensus map, thus establishing a positional reference for further genomic studies. Using the diversity lines, we estimated that a minimum density of one marker per 2 to 2.8 cM would be required for genome-wide association studies (GWAS), and GBS markers met this density requirement in most chromosome regions. We also demonstrated the utility of GBS in additional diagnostic applications related to oat breeding. We conclude that GBS is a powerful and useful approach, which will have many additional applications in oat breeding and genomic studies.