Genome-wide association mapping of fruit-quality traits using genotyping-by-sequencing approach in citrus landraces,modern cultivars,and breeding lines in Japan

Association mapping is an attractive method to identify QTLs in perennial horticultural crops such as citrus, as it does not need a designed cross between parental genotypes and can save time and labor to construct a segregating population. It usually requires more genetic markers than linkage-based QTL mapping owing to a lower degree of linkage disequilibrium (LD). However, recent advances in next-generation sequencing offer high-throughput, cost-effective methods, including genotyping-by-sequencing (GBS), for genotyping massive amounts of single nucleotide polymorphisms (SNPs). In this study, we performed a genome-wide association study (GWAS) of fruit-quality traits in citrus using SNPs obtained by GBS. We evaluated 110 citrus accessions, including landraces, modern cultivars, and breeding lines, for eight fruit-quality traits (fruit weight, fruit skin color, fruit surface texture, peelability, pulp firmness, segment firmness, sugar content, and acid content) during 2005 to 2012 (except 2007). GBS found 2309 SNPs, which we anchored to the clementine reference genome. We evaluated LD in the 110 accessions and confirmed that GBS gave enough SNPs to conduct GWAS. We identified seven QTLs, including four novel ones, comprising four significant QTLs for fruit weight and one QTL each for fruit skin color, pulp firmness, and segment firmness. These QTLs offer promise for use in citrus crossbreeding.     

High-throughput genotyping and mapping of single nucleotide polymorphisms in loblolly pine (Pinus taeda L.)

The development and application of genomic tools to loblolly pine (Pinus taeda L.) offer promising insights into the organization and structure of conifer genomes. The application of a high-throughput genotyping assay across diverse forest tree species, however, is currently limited taxonomically. This is despite the ongoing development of genome-scale projects aiming at the construction of expressed sequence tag (EST) libraries and the resequencing of EST-derived unigenes for a diverse array of forest tree species. In this paper, we report on the application of Illumina’s high-throughput GoldenGate™ SNP genotyping assay to a loblolly pine mapping population. Single nucleotide polymorphisms (SNPs) were identified through resequencing of previously identified wood quality, drought tolerance, and disease resistance candidate genes prior to genotyping. From that effort, a 384 multiplexed SNP assay was developed for high-throughput genotyping. Approximately 67% of the 384 SNPs queried converted into high-quality genotypes for the 48 progeny samples. Of those 257 successfully genotyped SNPs, 70 were segregating within the mapping population. A total of 27 candidate genes were subsequently mapped onto the existing loblolly pine consensus map, which consists of 12 linkage groups spanning a total map distance of 1,227.6 cM. The ability of SNPs to be mapped to the same position as fragment-based markers previously developed within the same candidate genes, as well as the pivotal role that SNPs currently play in the dissection of complex phenotypic traits, illustrate the usefulness of high-throughput SNP genotyping technologies to the continued development of pine genomics. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Identification of Polymorphic Markers by High-Resolution Melting (HRM) Assay for High-Throughput SNP Genotyping in Maize

The development of next generation sequencing (NGS) and high throughput genotyping are important techniques for the QTL mapping and genetic analysis of different crops. High-resolution melting (HRM) is an emerging technology used for detecting single-nucleotide polymorphisms (SNPs) in various species. However, its use is still limited in maize. The HRM analysis was integrated with SNPs to identify three types of populations (NIL population, RIL population and natural population), and the useful tags were screened. The patterns of temperature-shifted melting curves were investigated from the HRM analysis, and compared these with the kit. Among all 48 pairs of primers, 10 pairs of them were selected: six pairs of primers for the NIL population, three pairs of primers for the RIL population, and one pair of primer for the natural population. The marker for the natural population was developed with a matching rate of 80% for the plant height trait, based on the data of the phenotypic characteristics measured in the field. This study provides an effective method for maize genotyping in the classification of maize germplasm resources, which can be applied to other plants for high-throughput SNP genotyping or further mapping.     

Recent patents on high-throughput single nucleotide polymorphism (SNP) genotyping methods

Single nucleotide polymorphisms (SNPs) are single-base inheritable variations in a given and defined genetic location that occur in at least 1% of the population. SNPs are useful markers for genetic association studies in disease susceptibility or adverse drug reactions, in evolutionary studies and forensic science. Given the potential impact of SNPs, the biotechnology industry has focused on the development of high-throughput methods for SNP genotyping. Many highthroughput SNP genotyping technologies are currently available and many others are being patented recently. Each offers a unique combination of scale, accuracy, throughput and cost. In this review, we described some of the most important recent SNP genotyping methods and also recent patents associated with it.     

Development of a 135K SNP genotyping array for Actinidia arguta and its applications for genetic mapping and QTL analysis in kiwifruit

Kiwifruit (Actinidia spp) is a woody, perennial and deciduous vine. In this genus, there are multiple ploidy levels but the main cultivated cultivars are polyploid. Despite the availability of many genomic resources in kiwifruit, SNP genotyping is still a challenge given these different levels of polyploidy. Recent advances in SNP array technologies have offered a high-throughput genotyping platform for genome-wide DNA polymorphisms. In this study, we developed a high-density SNP genotyping array to facilitate genetic studies and breeding applications in kiwifruit. SNP discovery was performed by genome-wide DNA sequencing of 40 kiwifruit genotypes. The identified SNPs were stringently filtered for sequence quality, predicted conversion performance and distribution over the available Actinidia chinensis genome. A total of 134 729 unique SNPs were put on the array. The array was evaluated by genotyping 400 kiwifruit individuals. We performed a multidimensional scaling analysis to assess the diversity of kiwifruit germplasm, showing that the array was effective to distinguish kiwifruit accessions. Using a tetraploid F1 population, we constructed an integrated linkage map covering 3060.9 cM across 29 linkage groups and performed QTL analysis for the sex locus that has been identified on Linkage Group 3 (LG3) in Actinidia arguta. Finally, our dataset presented evidence of tetrasomic inheritance with partial preferential pairing in A. arguta. In conclusion, we developed and evaluated a 135K SNP genotyping array for kiwifruit. It has the advantage of a comprehensive design that can be an effective tool in genetic studies and breeding applications in this high-value crop.     

High density SNP mapping and QTL analysis for fruit quality characteristics in peach (Prunus persica L.)

Single nucleotide polymorphisms (SNPs) were used to construct an integrated SNP linkage map of peach (Prunus persica (L.) Batsch). A set of 1,536 SNPs were evaluated with the GoldenGate® Genotyping assay in two mapping populations, Pop-DF, and Pop-DG. After genotyping and filtering, a final set of 1,400 high quality SNPs in Pop-DF and 962 in Pop-DG with full map coverage were selected and used to construct two linkage maps with JoinMap®4.0. The Pop-DF map covered 422 cM of the peach genome and included 1,037 SNP markers, and Pop-DG map covered 369 cM and included 738 SNPs. A consensus map was constructed with 588 SNP markers placed in eight linkage groups (n?=?8 for peach), with map coverage of 454 cM and an average distance of 0.81 cM/marker site. Placements of SNPs on the “peach v1.0” physical map were compared to placement on the linkage maps and several differences were observed. Using the SNP linkage map of Pop-DG and phenotypic data collected for three harvest seasons, a QTL analysis for fruit quality traits and chilling injury symptoms was carried out with the mapped SNPs. Significant QTL effects were detected for mealiness (M) and flesh bleeding (FBL) QTLs on linkage group 4 and flesh browning (FBr) on linkage group 5. This study represents one of the first examples of QTL detection for quality traits and chilling injury symptoms using a high-density SNP map in a single peach F1 family.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

Highly efficient genotyping of rice biparental populations by GoldenGate assays based on parental resequencing

Key message

A new time- and cost-effective strategy was developed for medium-density SNP genotyping of rice biparental populations, using GoldenGate assays based on parental resequencing.

Abstract

Since the advent of molecular markers, crop researchers and breeders have dedicated huge amounts of effort to detecting quantitative trait loci (QTL) in biparental populations for genetic analysis and marker-assisted selection (MAS). In this study, we developed a new time- and cost-effective strategy for genotyping a population of progeny from a rice cross using medium-density single nucleotide polymorphisms (SNPs). Using this strategy, 728,362 “high quality” SNPs were identified by resequencing Teqing and Lemont, the parents of the population. We selected 384 informative SNPs that were evenly distributed across the genome for genotyping the biparental population using the Illumina GoldenGate assay. 335 (87.2 %) validated SNPs were used for further genetic analyses. After removing segregation distortion markers, 321 SNPs were used for linkage map construction and QTL mapping. This strategy generated SNP markers distributed more evenly across the genome than previous SSR assays. Taking the GW5 gene that controls grain shape as an example, our strategy provided higher accuracy (0.8 Mb) and significance (LOD 5.5 and 10.1) in QTL mapping than SSR analysis. Our study thus provides a rapid and efficient strategy for genetic studies and QTL mapping using SNP genotyping assays.     

Precision-mapping and statistical validation of quantitative trait loci by machine learning

Background

DNA sequence diversity within the human genome may be more greatly affected by copy number variations (CNVs) than single nucleotide polymorphisms (SNPs). Although the importance of CNVs in genome wide association studies (GWAS) is becoming widely accepted, the optimal methods for identifying these variants are still under evaluation. We have previously reported a comprehensive view of CNVs in the HapMap DNA collection using high density 500 K EA (Early Access) SNP genotyping arrays which revealed greater than 1,000 CNVs ranging in size from 1 kb to over 3 Mb. Although the arrays used most commonly for GWAS predominantly interrogate SNPs, CNV identification and detection does not necessarily require the use of DNA probes centered on polymorphic nucleotides and may even be hindered by the dependence on a successful SNP genotyping assay.

Results

In this study, we have designed and evaluated a high density array predicated on the use of non-polymorphic oligonucleotide probes for CNV detection. This approach effectively uncouples copy number detection from SNP genotyping and thus has the potential to significantly improve probe coverage for genome-wide CNV identification. This array, in conjunction with PCR-based, complexity-reduced DNA target, queries over 1.3 M independent NspI restriction enzyme fragments in the 200 bp to 1100 bp size range, which is a several fold increase in marker density as compared to the 500 K EA array. In addition, a novel algorithm was developed and validated to extract CNV regions and boundaries.

Conclusion

Using a well-characterized pair of DNA samples, close to 200 CNVs were identified, of which nearly 50% appear novel yet were independently validated using quantitative PCR. The results indicate that non-polymorphic probes provide a robust approach for CNV identification, and the increasing precision of CNV boundary delineation should allow a more complete analysis of their genomic organization.     

Clinical significance of germline copy number variation in susceptibility of human diseases

Germline copy number variation (CNV) is considered to be an important form of human genetic polymorphisms. Previous studies have identified amounts of CNVs in human genome by advanced technologies, such as comparative genomic hybridization, single nucleotide genotyping, and high-throughput sequencing. CNV is speculated to be derived from multiple mechanisms, such as nonallelic homologous recombination (NAHR) and nonhomologous end-joining (NHEJ). CNVs cover a much larger genome scale than single nucleotide polymorphisms (SNPs), and may alter gene expression levels by means of gene dosage, gene fusion, gene disruption, and long-range regulation effects, thus affecting individual phenotypes and playing crucial roles in human pathogenesis. The number of studies linking CNVs with common complex diseases has increased dramatically in recent years. Here, we provide a comprehensive review of the current understanding of germline CNVs, and summarize the association of germline CNVs with the susceptibility to a wide variety of human diseases that were identified in recent years. We also propose potential issues that should be addressed in future studies.     

Genome-wide SNP identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus): applications to QTL mapping of Vibrio anguillarum disease resistance and comparative genomic analysis

High-resolution genetic maps are essential for fine mapping of complex traits, genome assembly, and comparative genomic analysis. Single-nucleotide polymorphisms (SNPs) are the primary molecular markers used for genetic map construction. In this study, we identified 13,362 SNPs evenly distributed across the Japanese flounder (Paralichthys olivaceus) genome. Of these SNPs, 12,712 high-confidence SNPs were subjected to high-throughput genotyping and assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 3,497.29 cM with an average distance of 0.47 cM between loci, thereby representing the densest genetic map currently reported for Japanese flounder. Nine positive quantitative trait loci (QTLs) forming two main clusters for Vibrio anguillarum disease resistance were detected. All QTLs could explain 5.1–8.38% of the total phenotypic variation. Synteny analysis of the QTL regions on the genome assembly revealed 12 immune-related genes, among them 4 genes strongly associated with V. anguillarum disease resistance. In addition, 246 genome assembly scaffolds with an average size of 21.79 Mb were anchored onto the LGs; these scaffolds, comprising 522.99 Mb, represented 95.78% of assembled genomic sequences. The mapped assembly scaffolds in Japanese flounder were used for genome synteny analyses against zebrafish (Danio rerio) and medaka (Oryzias latipes). Flounder and medaka were found to possess almost one-to-one synteny, whereas flounder and zebrafish exhibited a multi-syntenic correspondence. The newly developed high-resolution genetic map, which will facilitate QTL mapping, scaffold assembly, and genome synteny analysis of Japanese flounder, marks a milestone in the ongoing genome project for this species.     

Large-scale SNP genotyping with canine buccal swab DNA

The dog is an attractive model for genetic studies of complex disease. With drafts of the canine genome complete, a large number of single-nucleotide polymorphisms (SNPs) that are potentially useful for gene-mapping studies and empirical estimations of canine diversity and linkage disequilibrium (LD) are now available. Unfortunately, most canine SNPs remain uncharacterized, and the amount and quality of DNA available from population-based samples are limited. We assessed how these real-world challenges influence automated SNP genotyping methods such as Illumina's GoldenGate assay. We examined 384 SNPs on canine chromosome 9 and successfully genotyped a minimum of 217 and a maximum of 275 SNPs using buccal swab samples for 181 dogs (86 beagles, 76 border collies, and 15 Australian shepherds). Call rates per SNP and sample averaged 97%, with reproducibility within and between analyses averaging 98%. The majority of these SNPs were polymorphic across all 3 breeds. We observed extensive LD, albeit less than reported for surveys using fewer dogs, consistent between breeds. Analyses of population substructure indicated that beagles are distinct from border collies and Australian shepherds. These results demonstrate the suitability of amplified canine buccal samples for high-throughput multiplex genotyping and confirm extensive LD in the dog.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

Applications of single nucleotide polymorphisms in crop genetics

The discovery of single nucleotide polymorphisms (SNPs) and insertions/deletions, which are the basis of most differences between alleles, has been simplified by recent developments in sequencing technology. SNP discovery in many crop species, such as corn and soybean, is relatively straightforward because of their high level of intraspecific nucleotide diversity, and the availability of many gene and expressed sequence tag (EST) sequences. For these species, direct readout of SNP haplotypes is possible. Haplotype-based analysis is more informative than analysis based on individual SNPs, and has more power in analyzing association with phenotypes. The elite germplasm of some crops may have been subjected to bottlenecks relatively recently, increasing the amount of linkage disequilibrium (LD) present and facilitating the association of SNP haplotypes at candidate gene loci with phenotypes. Whole-genome scans may help identify genome regions that are associated with interesting phenotypes if sufficient LD is present. Technological improvements make the use of SNP and indel markers attractive for high-throughput use in marker-assisted breeding, EST mapping and the integration of genetic and physical maps.     

Comparative LD mapping using single SNPs and haplotypes identifies QTL for plant height and biomass as secondary traits of drought tolerance in maize

Drought often delays developmental events so that plant height and above-ground biomass are reduced, resulting in yield loss due to inadequate photosynthate. In this study, plant height and biomass measured by the Normalized Difference Vegetation Index (NDVI) were used as criteria for drought tolerance. A total of 305 lines representing temperate, tropical and subtropical maize germplasm were genotyped using two single nucleotide polymorphism (SNP) chips each containing 1536 markers, from which 2052 informative SNPs and 386 haplotypes each constructed with two or more SNPs were used for linkage disequilibrium (LD) or association mapping. Single SNP- and haplotype-based LD mapping identified two significant SNPs and three haplotype loci [a total of four quantitative trait loci (QTL)] for plant height under well-watered and water-stressed conditions. For biomass, 32 SNPs and 12 haplotype loci (30 QTL) were identified using NDVIs measured at seven stages under the two water regimes. Some significant SNP and haplotype loci for NDVI were shared by different stages. Comparing significant loci identified by single SNP- and haplotype-based LD mapping, we found that six out of the 14 chromosomal regions defined by haplotype loci each included at least one significant SNP for the same trait. Significant SNP haplotype loci explained much higher phenotypic variation than individual SNPs. Moreover, we found that two significant SNPs (two QTL) and one haplotype locus were shared by plant height and NDVI. The results indicate the power of comparative LD mapping using single SNPs and SNP haplotypes with QTL shared by plant height and biomass as secondary traits for drought tolerance in maize.     

Development and evaluation of a 9K SNP array for peach by internationally coordinated SNP detection and validation in breeding germplasm

Although a large number of single nucleotide polymorphism (SNP) markers covering the entire genome are needed to enable molecular breeding efforts such as genome wide association studies, fine mapping, genomic selection and marker-assisted selection in peach [Prunus persica (L.) Batsch] and related Prunus species, only a limited number of genetic markers, including simple sequence repeats (SSRs), have been available to date. To address this need, an international consortium (The International Peach SNP Consortium; IPSC) has pursued a coordinated effort to perform genome-scale SNP discovery in peach using next generation sequencing platforms to develop and characterize a high-throughput Illumina Infinium® SNP genotyping array platform. We performed whole genome re-sequencing of 56 peach breeding accessions using the Illumina and Roche/454 sequencing technologies. Polymorphism detection algorithms identified a total of 1,022,354 SNPs. Validation with the Illumina GoldenGate® assay was performed on a subset of the predicted SNPs, verifying ∼75% of genic (exonic and intronic) SNPs, whereas only about a third of intergenic SNPs were verified. Conservative filtering was applied to arrive at a set of 8,144 SNPs that were included on the IPSC peach SNP array v1, distributed over all eight peach chromosomes with an average spacing of 26.7 kb between SNPs. Use of this platform to screen a total of 709 accessions of peach in two separate evaluation panels identified a total of 6,869 (84.3%) polymorphic SNPs.The almost 7,000 SNPs verified as polymorphic through extensive empirical evaluation represent an excellent source of markers for future studies in genetic relatedness, genetic mapping, and dissecting the genetic architecture of complex agricultural traits. The IPSC peach SNP array v1 is commercially available and we expect that it will be used worldwide for genetic studies in peach and related stone fruit and nut species.     

Entropy-based SNP selection for genetic association studies

Because of their abundance, density, and ease of practical use, single-nucleotide polymorphisms (SNPs) have become the major source of information for association gene mapping in humans. Sensible strategies for selecting practically useful SNPs are therefore required. Among the factors influencing the mapping utility of a given set of SNPs are (1) their individual diversity, (2) their haplotype structure in the population of interest, and (3) their physical distribution. We propose a strategy integrating these aspects into a single mapping utility measure, which is based upon Shannon entropy, and which maximizes the amount of information extracted from a genomic region under a Malecot model of linkage disequilibrium (LD) decay. The same utility measure has also been used to define a criterion guiding SNP discovery and rational decision-making about the continuation or termination of a mapping study. The proposed strategy performs consistently well in a data set comprising 549 German control individuals, genotyped for 136 SNPs from four genomic regions of different LD structure. Adoption of the method in practice is estimated to save up to 30% of genotyping load when compared with equidistant SNP localization or pair-wise LD minimization alone.     

The identification of SNPs with indeterminate positions using the Equine SNP50 BeadChip

We have used linkage disequilibrium (LD) to identify single nucleotide polymorphisms (SNPs) on the Illumina Equine SNP50 BeadChip, which may be incorrectly positioned on the genome map. A total of 1201 Thoroughbred horses were genotyped using the Illumina Equine SNP50 BeadChip. LD was evaluated in a pairwise fashion between all autosomal SNPs, both within and across chromosomes. Filters were then applied to the data, firstly to identify SNPs that may have been mapped to the wrong chromosome and secondly to identify SNPs that may have been incorrectly positioned within chromosomes. We identified a single SNP on ECA28, which showed low LD with neighbouring SNPs but considerable LD with a group of SNPs on ECA10. Furthermore, a cluster of SNPs on ECA5 showed unusually low LD with surrounding SNPs. A total of 39 SNPs met the criteria for unusual within-chromosome LD. The results of this study indicate that some SNPs may be misplaced. This finding is significant, as misplaced SNPs may lead to difficulties in the application of genomic methods, such as homozygosity mapping, for which SNP order is important.     

Copy number variation in the cattle genome

Copy number variations (CNVs) are gains and losses of genomic sequence greater than 50?bp between two individuals of a species. While single nucleotide polymorphisms (SNPs) are more frequent, CNVs impact a higher percentage of genomic sequence and have potentially greater effects, including the changing of gene structure and dosage, altering gene regulation and exposing recessive alleles. In particular, segmental duplications (SDs) were shown to be one of the catalysts and hotspots for CNV formation. Substantial progress has been made in understanding CNVs in mammals, especially in humans and rodents. CNVs have been shown to be important in both normal phenotypic variability and disease susceptibility. Recently, interest in CNV study has extended into domesticated animals, including cattle. Multiple genome-wide cattle CNV studies have been carried out using both microarray and next generation sequencing technologies. Integration of SD and CNV results with SNP and other datasets are beginning to reveal impacts of CNVs on cattle domestication, health, and production traits.     

MIG-seq is an effective method for high-throughput genotyping in wheat (Triticum spp.)

MIG-seq (Multiplexed inter-simple sequence repeats genotyping by sequencing) has been developed as a low cost genotyping technology, although the number of polymorphisms obtained is assumed to be minimal, resulting in the low application of this technique to analyses of agricultural plants. We applied MIG-seq to 12 plant species that include various crops and investigated the relationship between genome size and the number of bases that can be stably sequenced. The genome size and the number of loci, which can be sequenced by MIG-seq, are positively correlated. This is due to the linkage between genome size and the number of simple sequence repeats (SSRs) through the genome. The applicability of MIG-seq to population structure analysis, linkage mapping, and quantitative trait loci (QTL) analysis in wheat, which has a relatively large genome, was further evaluated. The results of population structure analysis for tetraploid wheat showed the differences among collection sites and subspecies, which agreed with previous findings. Additionally, in wheat biparental mapping populations, over 3,000 SNPs/indels with low deficiency were detected using MIG-seq, and the QTL analysis was able to detect recognized flowering-related genes. These results revealed the effectiveness of MIG-seq for genomic analysis of agricultural plants with large genomes, including wheat.