Development and Evaluation of SoySNP50K,a High-Density Genotyping Array for Soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35–45bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina’s manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.     

High-throughput genotyping in citrus accessions using an SNP genotyping array

We developed a 384 multiplexed SNP array, named CitSGA-1, for the genotyping of Citrus cultivars, and evaluated the performance and reliability of the genotyping. SNPs were surveyed by direct sequence comparison of the sequence tagged site (STS) fragment amplified from genomic DNA of cultivars representing the genetic diversity of citrus breeding in Japan. Among 1497 SNPs candidates, 384 SNPs for a high-throughput genotyping array were selected based on physical parameters of Illumina’s bead array criteria. The assay using CitSGA-1 was applied to a hybrid population of 88 progeny and 103 citrus accessions for breeding in Japan, which resulted in 73,726 SNP calls. A total of 351 SNPs (91 %) could call different genotypes among the DNA samples, resulting in a success rate for the assay comparable to previously reported rates for other plant species. To confirm the reliability of SNP genotype calls, parentage analysis was applied, and it indicated that the number of reliable SNPs and corresponding STSs were 276 and 213, respectively. The multiplexed SNP genotyping array reported here will be useful for the efficient construction of linkage map, for the detection of markers for marker-assisted breeding, and for the identification of cultivars.     

Identification of Capsicum species using SNP markers based on high resolution melting analysis

Single nucleotide polymorphisms (SNPs) derived from both nuclear and cytoplasmic DNA sequences were developed to identify distinct species of Capsicum. Species identification was achieved by detecting allelic variations of these type of markers via high resolution melting analysis (HRM). We used the HRM polymorphisms of COSII markers and the Waxy gene from the nuclear sequence, in addition to the intergenic spacer between trnL and trnF from cytoplasmic DNA as our SNP markers. A total of 31 accessions of Capsicum, representing six species, were analyzed using this method. As single markers were insufficient for identifying Capsicum species, combinations of all markers unambiguously identified all six. A phylogeny based on the SNP markers was consistent with the current taxonomy of Capsicum species. These observations demonstrate that the markers developed in this study are useful for rapid identification of new germplasm for management of Capsicum species.     

High resolution melting analysis of almond SNPs derived from ESTs

High resolution melting curve (HRM) is a recent advance for the detection of SNPs. The technique measures temperature induced strand separation of short PCR amplicons, and is able to detect variation as small as one base difference between samples. It has been applied to the analysis and scan of mutations in the genes causing human diseases. In plant species, the use of this approach is limited. We applied HRM analysis to almond SNP discovery and genotyping based on the predicted SNP information derived from the almond and peach EST database. Putative SNPs were screened from almond and peach EST contigs by HRM analysis against 25 almond cultivars. All 4 classes of SNPs, INDELs and microsatellites were discriminated, and the HRM profiles of 17 amplicons were established. The PCR amplicons containing single, double and multiple SNPs produced distinctive HRM profiles. Additionally, different genotypes of INDEL and microsatellite variations were also characterised by HRM analysis. By sequencing the PCR products, 100 SNPs were validated/revealed in the HRM amplicons and their flanking regions. The results showed that the average frequency of SNPs was 1:114 bp in the genic regions, and transition to transversion ratio was 1.16:1. Rare allele frequencies of the SNPs varied from 0.02 to 0.5, and the polymorphic information contents of the SNPs were from 0.04 to 0.53 at an average of 0.31. HRM has been demonstrated to be a fast, low cost, and efficient approach for SNP discovery and genotyping, in particular, for species without much genomic information such as almond.     

In silico single nucleotide polymorphism discovery and application to marker-assisted selection in soybean

The general approach to discovering single nucleotide polymorphisms (SNPs) requires locus-specific PCR amplification. To enhance the efficiency of SNP discovery in soybean, we used in silico analysis prior to re-sequencing as it is both rapid and inexpensive. In silico analysis was performed to detect putative SNPs in expressed sequence tag (EST) contigs assembled using publicly available ESTs from 18 different soybean genotypes. SNP validation by direct sequencing of six soybean cultivars and a wild soybean genotype was performed with PCR primers designed from EST contigs aligned with at least 5 out of 18 soybean genotypes. The efficiency of SNP discovery among the confirmation genotypes was 81.2%. Furthermore, the efficiency of SNP discovery between Pureunkong and Jinpumkong 2 genotypes was 47.4%, a great improvement on our previous finding based on direct sequencing (22.3%). Using SNPs between Pureunkong and Jinpumkong 2 in EST contigs, which were linked to target traits, we were able to genotype 90 recombinant inbred lines by high-resolution melting (HRM) analysis. These SNPs were mapped onto the expected locations near quantitative trait loci for water-logging tolerance and seed pectin concentration. Thus, our protocol for HRM analysis can be applied successfully not only to genetic diversity studies, but also to marker-assisted selection (MAS). Our study suggests that a combination of in silico analysis and HRM can reduce the cost and labor involved in developing SNP markers and genotyping SNPs. The markers developed in this study can also easily be applied to MAS if the markers are associated with the target traits.     

Geographic distribution and recombination of genomic fragments on the short arm of chromosome 2 of Arabidopsis thaliana

Abstract: Range expansion from Pleistocene refugia and anthropogenic influences contribute to the present distribution pattern of Arabidopsis thaliana. We scored a genome-wide set of CAPSs and found two markers with an east-west geographic distribution across the Eurasian range of the species. Regions around the two SNPs were sequenced in 98 accessions, including newly collected plants from Middle Asia and Western Siberia. These regions correspond to a gene (∼ 1500 bp) and a non-coding region (∼ 500 bp) 300 kbp apart on chromosome 2. Nucleotide diversities, π, of the two sequenced fragments were 0.0032 and 0.0130. The haplotypes of both sequences belonged to one of two groups: a rather uniform "Asian" and a more variable "European" haplotype group, on the basis of non-disjunct clusters of SNPs. Recombination between "Asian" and "European" haplotypes occurs where they meet. Especially in the "European" haplotype, many rare SNP variants representing independent mutations are scattered among the shared haplotype-specific SNPs. This agrees with previous suggestions of two large haplotype groups in A. thaliana and the post-glacial colonization of central Europe from the east and the west. A clear correlation between climatic factors and the haplotype distribution may reflect the dispersal history rather than local climate adaptation. The pattern of SNP variation within the contiguous sequences explains why only a minority of SNPs selected across the genome show evidence of this geographic pattern.     

Rapid genotyping of soybean cultivars using high throughput sequencing

Soybean (Glycine max) breeding involves improving commercially grown varieties by introgressing important agronomic traits from poor yielding accessions and/or wild relatives of soybean while minimizing the associated yield drag. Molecular markers associated with these traits are instrumental in increasing the efficiency of producing such crosses and Single Nucleotide Polymorphisms (SNPs) are particularly well suited for this task, owing to high density in the non-genic regions and thus increased likelihood of finding a tightly linked marker to a given trait. A rapid method to develop SNP markers that can differentiate specific loci between any two parents in soybean is thus highly desirable. In this study we investigate such a protocol for developing SNP markers between multiple soybean accessions and the reference Williams 82 genome. To restrict sampling frequency reduced representation libraries (RRLs) of genomic DNA were generated by restriction digestion followed by library construction. We chose to sequence four accessions Dowling (PI 548663), Dwight (PI 597386), Komata (PI200492) and PI 594538A for their agronomic importance as well as Williams 82 as a control.MseI was chosen to digest genomic DNA based on predictions that it will cut sparingly in the mathematically defined high-copy-number regions of the genome. All RRLs were sequenced on the Illumina genome analyzer. Reads were aligned to the Glyma1 reference assembly and SNP calls made from the alignments. We identified from 4294 to 14550 SNPs between the four accessions and the Williams 82 reference. In addition a small number of SNPs (1142) were found by aligning Williams 82 reads to the reference assembly (Glyma1) suggesting limited genetic variation within the Williams 82 line. The SNP data allowed us to estimate genetic diversity between the four lines and Williams 82. Restriction digestion of soybean genomic DNA with MseI followed by high throughput sequencing provides a rapid and reproducible method for generating SNP markers.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

Application of denaturing high-performance liquid chromatography for mapping of single nucleotide polymorphisms in barley (Hordeum vulgare L.).

Recent advances in DNA sequence analysis and the establishment of high-throughput assays have provided the framework for large-scale discovery and analysis of DNA sequence variation. In this context, single nucleotide polymorphisms (SNPs) are of particular interest. To initiate a systematic approach to develop an SNP map of barley (Hordeum vulgare L.), we have employed denaturing high-performance liquid chromatography (DHPLC) to analyse segregating SNP patterns in a doubled-haploid (DH) mapping population. To this end, SNPs between the parental genotypes were identified using a direct sequencing approach. Once a SNP was established between the parents, the optimal melting temperature of the PCR fragment containing the SNP was predicted for its analysis by DHPLC. Following the detection of the optimal temperature, the DH lines were analysed for the presence of either of the alleles. To test the utility of the analysis, data from previously mapped RFLP markers from which these SNPs were derived were compared. Results from these experiments indicate that DHPLC can be efficiently employed in analysing SNPs on a high-throughput scale.     

A high-throughput SNP marker system for parental polymorphism screening, and diversity analysis in common bean (Phaseolus vulgaris L.)

Single nucleotide polymorphism (SNP) detection has become a marker system of choice, because of the high abundance of source polymorphisms and the ease with which allele calls are automated. Various technologies exist for the evaluation of SNP loci and previously we validated two medium throughput technologies. In this study, our goal was to utilize a 768 feature, Illumina GoldenGate assay for common bean (Phaseolus vulgaris L.) developed from conserved legume gene sequences and to use the new technology for (1) the evaluation of parental polymorphisms in a mini-core set of common bean accessions and (2) the analysis of genetic diversity in the crop. A total of 736 SNPs were scored on 236 diverse common bean genotypes with the GoldenGate array. Missing data and heterozygosity levels were low and 94 % of the SNPs were scorable. With the evaluation of the parental polymorphism genotypes, we estimated the utility of the SNP markers in mapping for inter-genepool and intra-genepool populations, the latter being of lower polymorphism than the former. When we performed the diversity analysis with the diverse genotypes, we found Illumina GoldenGate SNPs to provide equivalent evaluations as previous gene-based SNP markers, but less fine-distinctions than with previous microsatellite marker analysis. We did find, however, that the gene-based SNPs in the GoldenGate array had some utility in race structure analysis despite the low polymorphism. Furthermore the SNPs detected high heterozygosity in wild accessions which was probably a reflection of ascertainment bias. The Illumina SNPs were shown to be effective in distinguishing between the genepools, and therefore were most useful in saturation of inter-genepool genetic maps. The implications of these results for breeding in common bean are discussed as well as the advantages and disadvantages of the GoldenGate system for SNP detection.     

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.     

2SNP: scalable phasing based on 2-SNP haplotypes

2SNP software package implements a new very fast scalable algorithm for haplotype inference based on genotype statistics collected only for pairs of SNPs. This software can be used for comparatively accurate phasing of large number of long genome sequences, e.g. obtained from DNA arrays. As an input 2SNP takes genotype matrix and outputs the corresponding haplotype matrix. On datasets across 79 regions from HapMap 2SNP is several orders of magnitude faster than GERBIL and PHASE while matching them in quality measured by the number of correctly phased genotypes, single-site and switching errors. For example, 2SNP requires 41 s on Pentium 4 2 Ghz processor to phase 30 genotypes with 1381 SNPs (ENm010.7p15:2 data from HapMap) versus GERBIL and PHASE requiring more than a week and admitting no less errors than 2SNP.     

A Definitive Haplotype Map as Determined by Genotyping Duplicated Haploid Genomes Finds a Predominant Haplotype Preference at Copy-Number Variation Events

The majority of complete hydatidiform moles (CHMs) harbor duplicated haploid genomes that originate from sperm. This makes CHMs more advantageous than conventional diploid cells for determining haplotypes of SNPs and copy-number variations (CNVs), because all of the genetic variants in a CHM genome are homozygous. Here we report SNP and CNV haplotype structures determined by analysis of 100 CHMs from Japanese subjects via high-density DNA arrays. The obtained haplotype map should be useful as a reference for the haplotype structure of Asian populations. We resolved common CNV regions (merged CNV segments across the examined samples) into CNV events (clusters of CNV segments) on the basis of mutual overlap and found that the haplotype backgrounds of different CNV events within the same CNV region were predominantly similar, perhaps because of inherent structural instability.     

Rapid genotyping of APOA5 -1131T>C polymorphism using high resolution melting analysis with unlabeled probes

The APOA5 -1131 T/C polymorphism (rs662799) exhibits a very strong association with elevated TG levels in different racial groups. High resolution melting (HRM) analysis with the use of unlabeled probes has shown to be a convenient and reliable tool to genotyping, but not yet been used for detecting rs662799 polymorphism. We applied the unlabeled probe HRM analysis and direct DNA sequencing to assay the -1131T>C SNP in 130 cases DNA samples blindly. This HRM analysis can be completed in <3 min for each sample. The two melting peaks were displayed at 66.1±0.4°C for CC homozygote and 68.7±0.2°C for TT homozygote; TC heterozygote showed the both melting peaks. The genotyping results by HRM method were completely concordant with direct DNA sequencing. The distribution of CC, TC, and TT genotypes for the -1131T>C SNP was 9.2, 49.2, and 41.5%, respectively. This assay was sensitive enough to detect C allele down to 20% and 10% for T allele. The limit of detection for C and T allele was 6.2 and 2.5 ng/μL DNA, respectively. The developed unlabeled probe HRM method provides an alternative mean to detect ApoA5 -1131T>C SNP rapidly and accurately.     

Structure of Pseudomonas aeruginosa populations analyzed by single nucleotide polymorphism and pulsed-field gel electrophoresis genotyping

Pseudomonas aeruginosa has a wide ecological distribution that includes natural habitats and clinical settings. To analyze the population structure and distribution of P. aeruginosa, a collection of 111 isolates of diverse habitats and geographical origin, most of which contained a genome with a different SpeI macrorestriction profile, was typed by restriction fragment length polymorphism based on 14 single nucleotide polymorphisms (SNPs) located at seven conserved loci of the core genome (oriC, oprL, fliC, alkB2, citS, oprI, and ampC). The combination of these SNPs plus the type of fliC present (a or b) allowed the assignment of a genetic fingerprint to each strain, thus providing a simple tool for the discrimination of P. aeruginosa strains. Thirteen of the 91 identified SNP genotypes were found in two or more strains. In several cases, strains sharing their SNP genotype had different SpeI macrorestriction profiles. The highly virulent CHA strain shared its SNP genotype with other strains that had different SpeI genotypes and which had been isolated from nonclinical habitats. The reference strain PAO1 also shared its SNP genotype with other strains that had different SpeI genotypes. The P. aeruginosa chromosome contains a conserved core genome and variable amounts of accessory DNA segments (genomic islands and islets) that can be horizontally transferred among strains. The fact that some SNP genotypes were overrepresented in the P. aeruginosa population studied and that several strains sharing an SNP genotype had different SpeI macrorestriction profiles supports the idea that changes occur at a higher rate in the accessory DNA segments than in the conserved core genome.     

Genome-wide haplotype changes produced by artificial selection during modern rice breeding in Japan

During the last 90 years, the breeding of rice has delivered cultivars with improved agronomic and economic characteristics. Crossing of different lines and successive artificial selection of progeny based on their phenotypes have changed the chromosomal constitution of the ancestors of modern rice; however, the nature of these changes is unclear. The recent accumulation of data for genome-wide single-nucleotide polymorphisms (SNPs) in rice has allowed us to investigate the change in haplotype structure and composition. To assess the impact of these changes during modern breeding, we studied 177 Japanese rice accessions, which were categorized into three groups: landraces, improved cultivars developed from 1931 to 1974 (the early breeding phase), and improved cultivars developed from 1975 to 2005 (the late breeding phase). Phylogenetic tree and structure analysis indicated genetic differentiation between non-irrigated (upland) and irrigated (lowland) rice groups as well as genetic structuring within the irrigated rice group that corresponded to the existence of three subgroups. Pedigree analysis revealed that a limited number of landraces and cultivars was used for breeding at the beginning of the period of systematic breeding and that 11 landraces accounted for 70% of the ancestors of the modern improved cultivars. The values for linkage disequilibrium estimated from SNP alleles and the haplotype diversity determined from consecutive alleles in five-SNP windows indicated that haplotype blocks became less diverse over time as a result of the breeding process. A decrease in haplotype diversity, caused by a reduced number of polymorphisms in the haplotype blocks, was observed in several chromosomal regions. However, our results also indicate that new haplotype polymorphisms have been generated across the genome during the breeding process. These findings will facilitate our understanding of the association between particular haplotypes and desirable phenotypes in modern Japanese rice cultivars.     

Discrimination of SNP genotypes associated with complex haplotypes by high resolution melting analysis in almond: implications for improved marker efficiencies

Developed recently, high resolution melting (HRM) analysis is an efficient, accurate and inexpensive method for distinguishing DNA polymorphisms. HRM has been used to identify mutations in human genes, and to detect SNPs, INDELs and microsatellites in plants. However, its capacity to discriminate DNA variants in the context of complex haplotypes involving INDEL as well as SNP variants has not been examined until now. In this study, we genotyped an almond (Prunus dulcis (Mill.) D. A. Webb, syn. Prunus amygdalus Batsch) pseudo-testcross mapping population that showed segregation of complex haplotypes associated with CYP79D16 promoter sequence. The 175 bp region in question included a 7 bp INDEL and 3 SNPs, and manifested as three different haplotypes in the parents. Thus, with one homozygous and one heterozygous parent, two relevant genotypes were identified in the mapping population. Although the population displayed monomorphism with respect to the INDEL and one of the SNPs, HRM was sufficiently sensitive to distinguish genotypes on the basis of the two informative SNPs, and the resulting data were used to map CYP79D16 to linkage group 6 of the almond genome. Thus the capacity of HRM to resolve genotypes arising from complex haplotypes has been demonstrated, and this has important implications for the design of efficient HRM markers for various genetic applications including mapping, population studies and biodiversity analyses.