Hierarchical high-throughput SNP genotyping of the human Y chromosome using MALDI-TOF mass spectrometry

We have established the use of a primer extension/mass spectrometry method (the PinPoint assay) for high-throughput SNP genotyping of the human Y chromosome. 118 markers were used to define 116 haplogroups and typing was organised in a hierarchical fashion. Twenty multiplex PCR/primer extension reactions were set up and each sample could be assigned to a haplogroup with only two to five of these multiplex analyses. A single aliquot of one enzyme was found to be sufficient for both PCR and primer extension. We observed 100% accuracy in blind validation tests. The technique thus provides a reliable, cost-effective and automated method for Y genotyping, and the advantages of using a hierarchical strategy can be applied to any DNA segment lacking recombination.     

A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome

SNP genotyping arrays have been useful for many applications that require a large number of molecular markers such as high-density genetic mapping, genome-wide association studies (GWAS), and genomic selection. We report the establishment of a large maize SNP array and its use for diversity analysis and high density linkage mapping. The markers, taken from more than 800,000 SNPs, were selected to be preferentially located in genes and evenly distributed across the genome. The array was tested with a set of maize germplasm including North American and European inbred lines, parent/F1 combinations, and distantly related teosinte material. A total of 49,585 markers, including 33,417 within 17,520 different genes and 16,168 outside genes, were of good quality for genotyping, with an average failure rate of 4% and rates up to 8% in specific germplasm. To demonstrate this array's use in genetic mapping and for the independent validation of the B73 sequence assembly, two intermated maize recombinant inbred line populations - IBM (B73×Mo17) and LHRF (F2×F252) - were genotyped to establish two high density linkage maps with 20,913 and 14,524 markers respectively. 172 mapped markers were absent in the current B73 assembly and their placement can be used for future improvements of the B73 reference sequence. Colinearity of the genetic and physical maps was mostly conserved with some exceptions that suggest errors in the B73 assembly. Five major regions containing non-colinearities were identified on chromosomes 2, 3, 6, 7 and 9, and are supported by both independent genetic maps. Four additional non-colinear regions were found on the LHRF map only; they may be due to a lower density of IBM markers in those regions or to true structural rearrangements between lines. Given the array's high quality, it will be a valuable resource for maize genetics and many aspects of maize breeding.     

Mutagenic primer design for mismatch PCR-RFLP SNP genotyping using a genetic algorithm

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) is useful in small-scale basic research studies of complex genetic diseases that are associated with single nucleotide polymorphism (SNP). Designing a feasible primer pair is an important work before performing PCR-RFLP for SNP genotyping. However, in many cases, restriction enzymes to discriminate the target SNP resulting in the primer design is not applicable. A mutagenic primer is introduced to solve this problem. GA-based Mismatch PCR-RFLP Primers Design (GAMPD) provides a method that uses a genetic algorithm to search for optimal mutagenic primers and available restriction enzymes from REBASE. In order to improve the efficiency of the proposed method, a mutagenic matrix is employed to judge whether a hypothetical mutagenic primer can discriminate the target SNP by digestion with available restriction enzymes. The available restriction enzymes for the target SNP are mined by the updated core of SNP-RFLPing. GAMPD has been used to simulate the SNPs in the human SLC6A4 gene under different parameter settings and compared with SNP Cutter for mismatch PCR-RFLP primer design. The in silico simulation of the proposed GAMPD program showed that it designs mismatch PCR-RFLP primers. The GAMPD program is implemented in JAVA and is freely available at http://bio.kuas.edu.tw/gampd/.     

High-throughput SNP genotyping with the GoldenGate assay in maize

Single nucleotide polymorphisms (SNPs) are abundant and evenly distributed throughout the genomes of most plant species. They have become an ideal marker system for genetic research in many crops. Several high throughput platforms have been developed that allow rapid and simultaneous genotyping of up to a million SNP markers. In this study, a custom GoldenGate assay containing 1,536 SNPs was developed based on public SNP information for maize and used to genotype two recombinant inbred line (RIL) populations (Zong3 x 87-1, and B73 x By804) and a panel of 154 diverse inbred lines. Over 90% of the SNPs were successfully scored in the diversity panel and the two RIL populations, with a genotyping error rate of less than 2%. A total of 975 SNP markers detected polymorphism in at least one of the two mapping populations, with a polymorphic rate of 38.5% in Zong3 x 87-1 and 52.6% in B73 x By804. The polymorphic SNPs in B73 x By804 have been integrated with previously mapped simple sequence repeat markers to construct a high-density linkage map containing 662 markers with a total length of 1,673.7 cM and an average of 2.53 cM between two markers. The minor allelic frequency (MAF) was distributed evenly across 10 continued classes from 0.05 to 0.5, and about 16% of the SNP markers had a MAF below 10% in the diversity panel. Polymorphism rates for individual SNP markers in pair-wise comparisons of genotypes tested ranged from 0.3 to 63.8% with an average of 36.3%. Most SNPs used in this GoldenGate assay appear to be equally useful for diversity analysis, marker-trait association studies, and marker-aided breeding.     

Construction of a genetic map based on high-throughput SNP genotyping and genetic mapping of a TuMV resistance locus in Brassica rapa

Brassica rapa is a member of the Brassicaceae family and includes vegetables and oil crops that are cultivated worldwide. The introduction of durable resistance against turnip mosaic virus (TuMV) into agronomically important cultivars has been a significant challenge for genetic and horticultural breeding studies of B. rapa. Based on our previous genome-wide analysis of DNA polymorphisms between the TuMV-resistant doubled haploid (DH) line VC40 and the TuMV-susceptible DH line SR5, we constructed a core genetic map of the VCS-13M DH population, which is composed of 83 individuals derived from microspore cultures of a F₁ cross between VC40 and SR5, by analyzing the segregation of 314 sequence-characterized genetic markers. The genetic markers correspond to 221 SNPs and 31 InDels of genes as well as 62 SSRs, covering 1,115.9 cM with an average distance of 3.6 cM between the adjacent marker loci. The alignment and orientation of the constructed map showed good agreement with the draft genome sequence of Chiifu, thus providing an efficient strategy to map genic sequences. Using the genetic map, a novel dominant TuMV resistance locus (TuMV-R) in the VCS-13M DH population was identified as a 0.34 Mb region in the short arm of chromosome A6 in which four CC–NBS–LRR resistance genes and two pathogenesis-related-1 genes reside. The genetic map developed in this study can play an important role in the genetic study of TuMV resistance and the molecular breeding of B. rapa.     

Genotools SNP manager: a new software for automated high-throughput MALDI-TOF mass spectrometry SNP genotyping

Analysis of single nucleotide polymorphisms (SNPs) is a rapidly growing field of research that provides insights into the most common type of differences between individual genomes. The resulting information has a strong impact in the fields of pharmacogenomics, drug development, forensic medicine, and diagnostics of specific disease markers. The technique of matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS) has been shown to be a highly suitable tool for the analysis of DNA. It supplies a very versatile method for addressing a high-throughput SNP genotyping approach. Here, we present the Bruker genotools SNP MANAGER, a new software tool suitable for highly automated MALDI-TOF MS SNP genotyping. The genotools SNP MANAGER administers the sample preparation data, calculates masses of allele-specific primer extension products, performs genotyping analysis, and displays the results. In the current study, we have used the genotools SNP MANAGER to perform an automated duplex SNP analysis of two biallelic markers from the promoter of the gene encoding the inflammatory mediator interleukin-6.     

Linkage disequilibrium in European elite maize germplasm investigated with SSRs

Information about the extent and genomic distribution of linkage disequilibrium (LD) is of fundamental importance for association mapping. The main objectives of this study were to (1) investigate genetic diversity within germplasm groups of elite European maize (Zea mays L.) inbred lines, (2) examine the population structure of elite European maize germplasm, and (3) determine the extent and genomic distribution of LD between pairs of simple sequence repeat (SSR) markers. We examined genetic diversity and LD in a cross section of European and US elite breeding material comprising 147 inbred lines genotyped with 100 SSR markers. For gene diversity within each group, significant (P<0.05) differences existed among the groups. The LD was significant (P<0.05) for 49% of the SSR marker pairs in the 80 flint lines and for 56% of the SSR marker pairs in the 57 dent lines. The ratio of linked to unlinked loci in LD was 1.1 for both germplasm groups. The high incidence of LD suggests that the extent of LD between SSR markers should allow the detection of marker-phenotype associations in a genome scan. However, our results also indicate that a high proportion of the observed LD is generated by forces, such as relatedness, population stratification, and genetic drift, which cause a high risk of detecting false positives in association mapping.     

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.     

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.     

Estimating autozygosity from high-throughput information: effects of SNP density and genotyping errors

Background

Runs of homozygosity are long, uninterrupted stretches of homozygous genotypes that enable reliable estimation of levels of inbreeding (i.e., autozygosity) based on high-throughput, chip-based single nucleotide polymorphism (SNP) genotypes. While the theoretical definition of runs of homozygosity is straightforward, their empirical identification depends on the type of SNP chip used to obtain the data and on a number of factors, including the number of heterozygous calls allowed to account for genotyping errors. We analyzed how SNP chip density and genotyping errors affect estimates of autozygosity based on runs of homozygosity in three cattle populations, using genotype data from an SNP chip with 777 972 SNPs and a 50 k chip.

Results

Data from the 50 k chip led to overestimation of the number of runs of homozygosity that are shorter than 4 Mb, since the analysis could not identify heterozygous SNPs that were present on the denser chip. Conversely, data from the denser chip led to underestimation of the number of runs of homozygosity that were longer than 8 Mb, unless the presence of a small number of heterozygous SNP genotypes was allowed within a run of homozygosity.

Conclusions

We have shown that SNP chip density and genotyping errors introduce patterns of bias in the estimation of autozygosity based on runs of homozygosity. SNP chips with 50 000 to 60 000 markers are frequently available for livestock species and their information leads to a conservative prediction of autozygosity from runs of homozygosity longer than 4 Mb. Not allowing heterozygous SNP genotypes to be present in a homozygosity run, as has been advocated for human populations, is not adequate for livestock populations because they have much higher levels of autozygosity and therefore longer runs of homozygosity. When allowing a small number of heterozygous calls, current software does not differentiate between situations where these calls are adjacent and therefore indicative of an actual break of the run versus those where they are scattered across the length of the homozygous segment. Simple graphical tests that are used in this paper are a current, yet tedious solution.     

SNP genotyping using single-tube fluorescent bidirectional PCR

SNP genotyping is a well-populatedfield with a large number of assay formats offering accurate allelic discrimination. However, there remains a discord between the ultimate goal of rapid, inexpensive assays that do not require complex design considerations and involved optimization strategies. We describe the first integration of bidirectional allele-specific amplification, SYBR Green I, and rapid-cycle PCR to provide a homogeneous SNP-typing assay. Wild-type, mutant, and heterozygous alleles were easily discriminated in a single tube using melt curve profiling of PCR products alone. We demonstrate the effectiveness and reliability of this assay with a blinded trial using clinical samples from individuals with sickle cell anemia, sickle cell trait, or unaffected individuals. The tests were completed in less than 30 min without expensive fluorogenic probes, prohibiting design rules, or lengthy downstream processing for product analysis.     

Automated SNP genotype clustering algorithm to improve data completeness in high-throughput SNP genotyping datasets from custom arrays

High-throughput SNP genotyping platforms use automated genotype calling algorithms to assign genotypes. While these algorithms work efficiently for individual platforms, they are not compatible with other platforms, and have individual biases that result in missed genotype calls. Here we present data on the use of a second complementary SNP genotype clustering algorithm. The algorithm was originally designed for individual fluorescent SNP genotyping assays, and has been optimized to permit the clustering of large datasets generated from custom-designed Affymetrix SNP panels. In an analysis of data from a 3K array genotyped on 1,560 samples, the additional analysis increased the overall number of genotypes by over 45,000, significantly improving the completeness of the experimental data. This analysis suggests that the use of multiple genotype calling algorithms may be advisable in high-throughput SNP genotyping experiments. The software is written in Perl and is available from the corresponding author.     

Identification of functional genetic variations underlying drought tolerance in maize using SNP markers

Single nucleotide polymorphism (SNP) is a common form of genetic variation and popularly exists in maize genome. An Illumina GoldenGate assay with 1 536 SNP markers was used to genotype maize inbred lines and identified the functional genetic variations underlying drought tolerance by association analysis. Across 80 lines, 1 006 polymorphic SNPs (65.5% of the total) in the assay with good call quality were used to estimate the pattern of genetic diversity, population structure, and familial relatedness. The analysis showed the best number of fixed subgroups was six, which was consistent with their original sources and results using only simple sequence repeat markers. Pairwise linkage disequilibrium (LD) and association mapping with phenotypic traits investigated under water-stressed and well-watered regimes showed rapid LD decline within 100-500 kb along the physical distance of each chromosome, and that 29 SNPs were associated with at least two phenotypic traits in one or more environments, which were related to drought-tolerant or drought-responsive genes. These drought-tolerant SNPs could be converted into functional markers and then used for maize improvement by marker-assisted selection.     

Objective evaluation measures of genetic marker selection in large-scale SNP genotyping

High-throughput single nucleotide polymorphism (SNP) genotyping systems provide two kinds of fluorescent signals detected from different alleles. In current technologies, the process of genotype discrimination requires subjective judgments by expert operators, even when using clustering algorithms. Here, we propose two evaluation measures to manage fluorescent scatter data with nonclear plot aggregation. The first is the marker ranking measure, which provides a ranking system for the SNP markers based on the distance between the scatter plot distribution and a user-defined ideal distribution. The second measure, called individual genotype membership, uses the membership probability of each genotype related to an individual plot in the scatter data. In verification experiments, the marker ranking measure determined the ranking of SNP markers correlated with the subjective order of SNP markers judged by an expert operator. The experiment using the individual genotype membership measure clarified that the total number of unclassified individuals was remarkably reduced compared to that of manually unclassified ones. These two evaluation measures were implemented as the GTAssist software. GTAssist provides objective standards and avoids subjective biases in SNP genotyping workflows.     

Potential of SNP markers for the characterization of Brazilian cassava germplasm

Key message

High-throughput markers, such as SNPs, along with different methodologies were used to evaluate the applicability of the Bayesian approach and the multivariate analysis in structuring the genetic diversity in cassavas.

Abstract

The objective of the present work was to evaluate the diversity and genetic structure of the largest cassava germplasm bank in Brazil. Complementary methodological approaches such as discriminant analysis of principal components (DAPC), Bayesian analysis and molecular analysis of variance (AMOVA) were used to understand the structure and diversity of 1,280 accessions genotyped using 402 single nucleotide polymorphism markers. The genetic diversity (0.327) and the average observed heterozygosity (0.322) were high considering the bi-allelic markers. In terms of population, the presence of a complex genetic structure was observed indicating the formation of 30 clusters by DAPC and 34 clusters by Bayesian analysis. Both methodologies presented difficulties and controversies in terms of the allocation of some accessions to specific clusters. However, the clusters suggested by the DAPC analysis seemed to be more consistent for presenting higher probability of allocation of the accessions within the clusters. Prior information related to breeding patterns and geographic origins of the accessions were not sufficient for providing clear differentiation between the clusters according to the AMOVA analysis. In contrast, the F _ST was maximized when considering the clusters suggested by the Bayesian and DAPC analyses. The high frequency of germplasm exchange between producers and the subsequent alteration of the name of the same material may be one of the causes of the low association between genetic diversity and geographic origin. The results of this study may benefit cassava germplasm conservation programs, and contribute to the maximization of genetic gains in breeding programs.     

Population genetic analysis and phylogeny reconstruction in Eucalyptus (Myrtaceae) using high-throughput, genome-wide genotyping

A set of over 8000 Diversity Arrays Technology (DArT) markers was tested for its utility in high-resolution population and phylogenetic studies across a range of Eucalyptus taxa. Small-scale population studies of Eucalyptus camaldulensis, Eucalyptus cladocalyx, Eucalyptus globulus, Eucalyptus grandis, Eucalyptus nitens, Eucalyptus pilularis and Eucalyptus urophylla demonstrated the potential of genome-wide genotyping with DArT markers to differentiate species, to identify interspecific hybrids and to resolve biogeographic disjunctions within species. The population genetic studies resolved geographically partitioned clusters in E. camaldulensis, E. cladocalyx, E. globulus and E. urophylla that were congruent with previous molecular studies. A phylogenetic study of 94 eucalypt species provided results that were largely congruent with traditional taxonomy and ITS-based phylogenies, but provided more resolution within major clades than had been obtained previously. Ascertainment bias (the bias introduced in a phylogeny from using markers developed in a small sample of the taxa that are being studied) was not detected. DArT offers an unprecedented level of resolution for population genetic, phylogenetic and evolutionary studies across the full range of Eucalyptus species.     

Biophysical characterization of double-stranded oligonucleotides using ETBR and isothermal fluorescence spectroscopy: implication for SNP genotyping

The UV-absorption, fluorescence and CD spectra of aps 23 bp oligoduplexes were performed for potential diagnostic purpose. These oligonucleotide sequences were mimicked from natural mutations (mitochondrial genome) of human population (unpublished). This work was designed on the basis of hybridization of non-self complementary oligoduplexes (aps) containing no mismatch, one-mismatch and two-mismatches. Since melting temperature™ is dependent on concentration of the oligoduplex, various concentrations were used in this study protocol. The thermal spectra profiles (UV absorbance and fluorescence) of these oligoduplexes (aps) are different for a particular concentration, and can be implicated for mutations. − dF/dT (or dA/dT) vs T, lnK (or RlnK) vs T_M, ΔG vs T_M, ΔS vs T_M and ΔH vs T_M are also variable for those sequences. All these thermodynamic data were calculated from absorbance (at 260 nm) data. On the contrary to the 23 bp oligoduplexes (aps), the PCR products of 97 bp and 256 bp length were genotyped with ETBR (excitation 530 nm, emission 600 nm) fluorimetrically. But our attempts to genotype these PCR sequences with isothermal UV absorbance spectroscopy were unsuccessful. Isothermal UV absorbance spectra has a limitation of sequence length. However, the structural conformation (all B-type) of the oligoduplexes (aps) was determined using CD. The minor discrepancy in CD spectra of these oligoduplexes are not significant for mutational analysis. 97 bp nested PCR product was an amplicon having either GcT or AcC mutation of mitochondria of normal human population, whereas 256 bp PCR product was an amplicon of human BRCA2 gene (NCBI Accession No. AY151039) of chromosome 13 having either A or G mutation at position − 26.     

High-throughput retrotransposon-based genetic diversity of maize germplasm assessment and analysis

Maize is one of the world’s most important crops and a model for grass genome research. Long terminal repeat (LTR) retrotransposons comprise most of the maize genome; their ability to produce new copies makes them efficient high-throughput genetic markers. Inter-retrotransposon-amplified polymorphisms (IRAPs) were used to study the genetic diversity of maize germplasm. Five LTR retrotransposons (Huck, Tekay, Opie, Ji, and Grande) were chosen, based on their large number of copies in the maize genome, whereas polymerase chain reaction primers were designed based on consensus LTR sequences. The LTR primers showed high quality and reproducible DNA fingerprints, with a total of 677 bands including 392 polymorphic bands showing 58% polymorphism between maize hybrid lines. These markers were used to identify genetic similarities among all lines of maize. Analysis of genetic similarity was carried out based on polymorphic amplicon profiles and genetic similarity phylogeny analysis. This diversity was expected to display ecogeographical patterns of variation and local adaptation. The clustering method showed that the varieties were grouped into three clusters differing in ecogeographical origin. Each of these clusters comprised divergent hybrids with convergent characters. The clusters reflected the differences among maize hybrids and were in accordance with their pedigree. The IRAP technique is an efficient high-throughput genetic marker-generating method.

     

SNP frequency, haplotype structure and linkage disequilibrium in elite maize inbred lines

Background  

Recent studies of ancestral maize populations indicate that linkage disequilibrium tends to dissipate rapidly, sometimes within 100 bp. We set out to examine the linkage disequilibrium and diversity in maize elite inbred lines, which have been subject to population bottlenecks and intense selection by breeders. Such population events are expected to increase the amount of linkage disequilibrium, but reduce diversity. The results of this study will inform the design of genetic association studies.