An efficient method for multi-locus molecular haplotyping

Many methods exist for genotyping—revealing which alleles an individual carries at different genetic loci. A harder problem is haplotyping—determining which alleles lie on each of the two homologous chromosomes in a diploid individual. Conventional approaches to haplotyping require the use of several generations to reconstruct haplotypes within a pedigree, or use statistical methods to estimate the prevalence of different haplotypes in a population. Several molecular haplotyping methods have been proposed, but have been limited to small numbers of loci, usually over short distances. Here we demonstrate a method which allows rapid molecular haplotyping of many loci over long distances. The method requires no more genotypings than pedigree methods, but requires no family material. It relies on a procedure to identify and genotype single DNA molecules, and reconstruction of long haplotypes by a ‘tiling’ approach. We demonstrate this by resolving haplotypes in two regions of the human genome, harbouring 20 and 105 single-nucleotide polymorphisms, respectively. The method can be extended to reconstruct haplotypes of arbitrary complexity and length, and can make use of a variety of genotyping platforms. We also argue that this method is applicable in situations which are intractable to conventional approaches.     

Molecular haplotyping at high throughput

Reconstruction of haplotypes, or the allelic phase, of single nucleotide polymorphisms (SNPs) is a key component of studies aimed at the identification and dissection of genetic factors involved in complex genetic traits. In humans, this often involves investigation of SNPs in case/control or other cohorts in which the haplotypes can only be partially inferred from genotypes by statistical approaches with resulting loss of power. Moreover, alternative statistical methodologies can lead to different evaluations of the most probable haplotypes present, and different haplotype frequency estimates when data are ambiguous. Given the cost and complexity of SNP studies, a robust and easy-to-use molecular technique that allows haplotypes to be determined directly from individual DNA samples would have wide applicability. Here, we present a reliable, automated and high-throughput method for molecular haplotyping in 2 kb, and potentially longer, sequence segments that is based on the physical determination of the phase of SNP alleles on either of the individual paternal haploids. We demonstrate that molecular haplotyping with this technique is not more complicated than SNP genotyping when implemented by matrix-assisted laser desorption/ionisation mass spectrometry, and we also show that the method can be applied using other DNA variation detection platforms. Molecular haplotyping is illustrated on the well-described β₂-adrenergic receptor gene.     

Assaying chromosomal inversions by single-molecule haplotyping

Inversions are an important form of structural variation, but they are difficult to characterize, as their breakpoints often fall within inverted repeats. We have developed a method called 'haplotype fusion' in which an inversion breakpoint is genotyped by performing fusion PCR on single molecules of human genomic DNA. Fusing single-copy sequences bracketing an inversion breakpoint generates orientation-specific PCR products, exemplified by a genotyping assay for the int22 hemophilia A inversion on Xq28. Furthermore, we demonstrated that inversion events with breakpoints embedded within long (>100 kb) inverted repeats can be genotyped by haplotype-fusion PCR followed by bead-based single-molecule haplotyping on repeat-specific markers bracketing the inversion breakpoint. We illustrate this method by genotyping a Yp paracentric inversion sponsored by >300-kb-long inverted repeats. The generality of our methods to survey for, and genotype chromosomal inversions should help our understanding of the contribution of inversions to genomic variation, inherited diseases and cancer.     

Haplotyping using a combination of polymerase chain reaction–single-strand conformational polymorphism analysis and haplotype-specific PCR amplification

A single nucleotide polymorphism (SNP) may have an impact on phenotype, but it may also be influenced by multiple SNPs within a gene; hence, the haplotype or phase of multiple SNPs needs to be known. Various methods for haplotyping SNPs have been proposed, but a simple and cost-effective method is currently unavailable. Here we describe a haplotyping approach using two simple techniques: polymerase chain reaction–single-strand conformational polymorphism (PCR–SSCP) and haplotype-specific PCR. In this approach, individual regions of a gene are analyzed by PCR–SSCP to identify variation that defines sub-haplotypes, and then extended haplotypes are assembled from the sub-haplotypes either directly or with the additional use of haplotype-specific PCR amplification. We demonstrate the utility of this approach by haplotyping ovine FABP4 across two variable regions that contain seven SNPs and one indel. The simplicity of this approach makes it suitable for large-scale studies and/or diagnostic screening.     

A locus-wide approach to assessing variation in the avian MHC: the B-locus of the wild turkey

Studies of major histocompatibility complex (MHC) diversity in non-model vertebrates typically focus on structure and sequence variation in the antigen-presenting loci: the highly variable and polymorphic class I and class IIB genes. Although these studies provide estimates of the number of genes and alleles/locus, they often overlook variation in functionally related and co-inherited genes important in the immune response. This study utilizes the sequence of the MHC B-locus derived from a commercial turkey to investigate MHC variation in wild birds. Sequences were obtained for nine interspersed MHC amplicons (non-class I/II) from each of 40 birds representing 3 subspecies of wild turkey (Meleagris gallopavo). Analysis of aligned sequences identified 238 single-nucleotide variants approximately one-third of which had minor allele frequencies >0.2 in the sampled birds. PHASE analysis identified 70 prospective MHC haplotypes in the wild turkeys, whereas a combined analysis with commercial birds identified almost 100 haplotypes in the species. Denaturing gradient gel electrophoresis (DGGE) of the class IIB loci was used to test the efficacy of single-nucleotide polymorphism (SNP) haplotyping to capture locus-wide variation. Diversity in SNP haplotypes and haplotype sharing among individuals was directly reflected in the DGGE patterns. Utilization of a reference haplotype to sequence interspersed regions of the MHC has significant advantages over other methods of surveying diversity while identifying high-frequency SNPs for genotyping. SNP haplotyping provides a means to identify both divergent haplotypes and homozygous individuals for assessment of immunological variation in wild and domestic populations.     

Linkage disequilibrium-based quality control for large-scale genetic studies

Quality control (QC) is a critical step in large-scale studies of genetic variation. While, on average, high-throughput single nucleotide polymorphism (SNP) genotyping assays are now very accurate, the errors that remain tend to cluster into a small percentage of "problem" SNPs, which exhibit unusually high error rates. Because most large-scale studies of genetic variation are searching for phenomena that are rare (e.g., SNPs associated with a phenotype), even this small percentage of problem SNPs can cause important practical problems. Here we describe and illustrate how patterns of linkage disequilibrium (LD) can be used to improve QC in large-scale, population-based studies. This approach has the advantage over existing filters (e.g., HWE or call rate) that it can actually reduce genotyping error rates by automatically correcting some genotyping errors. Applying this LD-based QC procedure to data from The International HapMap Project, we identify over 1,500 SNPs that likely have high error rates in the CHB and JPT samples and estimate corrected genotypes. Our method is implemented in the software package fastPHASE, available from the Stephens Lab website (http://stephenslab.uchicago.edu/software.html).     

Direct molecular haplotyping by melting curve analysis of hybridization probes: beta 2-adrenergic receptor haplotypes as an example

Direct determination of the association of multiple genetic polymorphisms, or haplotyping, in individual samples is challenging because of chromosome diploidy. Here, we describe the ability of hybridization probes, commonly used as genotyping tools, to establish single nucleotide polymorphism (SNP) haplotypes in a single step. Three haplotypes found in the beta 2-adrenergic receptor (β2AR) gene and characterized by three different SNPs combinations are presented as examples. Each combination of SNPs has a unique stability, recorded by its melting temperature, even when intervening sequences from the template must loop out during probe hybridization. In the course of this study, two haplotypes in β2AR not described previously were discovered. This approach provides a tool for molecular haplotyping that should prove useful in clinical molecular genetics diagnostics and pharmacogenetic research where methods for direct haplotyping are needed.     

Molecular haplotyping of tandem single nucleotide polymorphisms by allele-specific PCR

BARD1 Val507Met (1592A>G) is an interesting marker for association studies on cancer risk. However, studies are scarce in the literature, probably reflecting the methodological problem imposed by the fact that next to the 1592A>G stands the 1591C>T single nucleotide polymorphism (SNP). We have designed an allele-specific PCR method capable of molecular haplotyping tandem SNPs. In the tandem SNPs haplotyping assay (tSNPh), four reverse primers are designed to be perfect matches of each potential haplotype. The forward primer is labeled with a fluorochrome. PCR products are analyzed by capillary electrophoresis. Haplotyping is performed by size calling. To ascertain the accuracy and reproducibility of the assay, we measured the level of concordance with sequencing data in 124 samples. In vitro-generated templates have been used for further testing. We developed a novel and reliable assay that permits typing two SNPs directly adjacent to each other, avoiding mutual interferences. The method is amenable to automation and high throughput. We expect that this assay will contribute to clarifying the role of BARD1 in cancer susceptibility. In addition, we suggest that tandem SNPs are potentially interesting polymorphic markers in which molecular haplotyping can be performed easily.     

Selection of genetic markers for association analyses,using linkage disequilibrium and haplotypes

The genotyping of closely spaced single-nucleotide polymorphism (SNP) markers frequently yields highly correlated data, owing to extensive linkage disequilibrium (LD) between markers. The extent of LD varies widely across the genome and drives the number of frequent haplotypes observed in small regions. Several studies have illustrated the possibility that LD or haplotype data could be used to select a subset of SNPs that optimize the information retained in a genomic region while reducing the genotyping effort and simplifying the analysis. We propose a method based on the spectral decomposition of the matrices of pairwise LD between markers, and we select markers on the basis of their contributions to the total genetic variation. We also modify Clayton's "haplotype tagging SNP" selection method, which utilizes haplotype information. For both methods, we propose sliding window-based algorithms that allow the methods to be applied to large chromosomal regions. Our procedures require genotype information about a small number of individuals for an initial set of SNPs and selection of an optimum subset of SNPs that could be efficiently genotyped on larger numbers of samples while retaining most of the genetic variation in samples. We identify suitable parameter combinations for the procedures, and we show that a sample size of 50-100 individuals achieves consistent results in studies of simulated data sets in linkage equilibrium and LD. When applied to experimental data sets, both procedures were similarly effective at reducing the genotyping requirement while maintaining the genetic information content throughout the regions. We also show that haplotype-association results that Hosking et al. obtained near CYP2D6 were almost identical before and after marker selection.     

Improved allelic differentiation using sequence-specific oligonucleotide hybridization incorporating an additional base-analogue mismatch

Sequence-specific oligonucleotide hybridization (SSOH, 'dot-blotting') is a widely employed method of typing single nucleotide polymorphisms (SNPs), but it is often compromised by lack of allelic differentiation. We describe a novel improvement to SSOH that incorporates an additional mismatch into the oligonucleotide probe using the universal base analogue 3-nitropyrrole. This method greatly increases allelic differentiation compared to standard SSOH where oligonucleotides contain only SNP-defining base changes. Moreover, stringency of the hybridisation is predictably maintained over a wide range of temperatures, which can be calculated empirically, thus facilitating the genotyping of multiple SNPs using similar conditions. This improved method increases the usefulness of hybridisation-based methods of rapid genotyping of SNPs and may have implications for array methodologies.     

Improved Allelic Differentiation Using Sequence‐Specific Oligonucleotide Hybridization Incorporating an Additional Base‐Analogue Mismatch

Sequence‐specific oligonucleotide hybridization (SSOH, ‘dot‐blotting’) is a widely employed method of typing single nucleotide polymorphisms (SNPs), but it is often compromised by lack of allelic differentiation. We describe a novel improvement to SSOH that incorporates an additional mismatch into the oligonucleotide probe using the universal base analogue 3‐nitropyrrole. This method greatly increases allelic differentiation compared to standard SSOH where oligonucleotides contain only SNP‐defining base changes. Moreover, stringency of the hybridisation is predictably maintained over a wide range of temperatures, which can be calculated empirically, thus facilitating the genotyping of multiple SNPs using similar conditions. This improved method increases the usefulness of hybridisation‐based methods of rapid genotyping of SNPs and may have implications for array methodologies.     

A pipeline for high throughput detection and mapping of SNPs from EST databases

Single nucleotide polymorphisms (SNPs) represent the most abundant type of genetic variation that can be used as molecular markers. The SNPs that are hidden in sequence databases can be unlocked using bioinformatic tools. For efficient application of these SNPs, the sequence set should be error-free as much as possible, targeting single loci and suitable for the SNP scoring platform of choice. We have developed a pipeline to effectively mine SNPs from public EST databases with or without quality information using QualitySNP software, select reliable SNP and prepare the loci for analysis on the Illumina GoldenGate genotyping platform. The applicability of the pipeline was demonstrated using publicly available potato EST data, genotyping individuals from two diploid mapping populations and subsequently mapping the SNP markers (putative genes) in both populations. Over 7000 reliable SNPs were identified that met the criteria for genotyping on the GoldenGate platform. Of the 384 SNPs on the SNP array approximately 12% dropped out. For the two potato mapping populations 165 and 185 SNPs segregating SNP loci could be mapped on the respective genetic maps, illustrating the effectiveness of our pipeline for SNP selection and validation.     

Direct micro-haplotyping by multiple double PCR amplifications of specific alleles (MD-PASA)

Analysis of haplotypes is an important tool in population genetics, familial heredity and gene mapping. Determination of haplotypes of multiple single nucleotide polymorphisms (SNPs) or other simple mutations is time consuming and expensive when analyzing large populations, and often requires the help of computational and statistical procedures. Based on double PCR amplification of specific alleles, described previously, we have developed a simple, rapid and low-cost method for direct haplotyping of multiple SNPs and simple mutations found within relatively short specific regions or genes (micro-haplotypes). Using this method, it is possible to directly determine the physical linkage of multiple heterozygous alleles, by conducting a series of double allele-specific PCR amplification sets with simple analysis by gel electrophoresis. Application of the method requires prior information as to the sequence of the segment to be haplotyped, including the polymorphic sites. We applied the method to haplotyping of nine sites in the chicken HSP108 gene. One of the haplotypes in the population apparently arose by recombination between two existing haplotypes, and we were able to locate the point of recombination within a segment of 19 bp. We anticipate rapidly growing needs for SNP haplotyping in human (medical and pharmacogenetics), animal and plant genetics; in this context, the multiple double PCR amplifications of specific alleles (MD-PASA) method offers a useful haplotyping tool.     

SNiPORK - a microarray of SNPs in candidate genes potentially associated with pork yield and quality - development and validation in commercial breeds

SNiPORK is an oligonucleotide microarray based on the arrayed primer extension (APEX) technique, allowing genotyping of single nucleotide polymorphisms (SNPs) in genes of interest for pork yield and quality traits. APEX consists of a sequencing reaction primed by an oligonucleotide anchored with its 5' end to a glass slide and terminating one nucleotide before the polymorphic site. Extension with one fluorescently labeled dideoxynucleotide complementary to the template reveals the polymorphism. Ninety SNPs were selected from those associated directly or potentially with pork traits. Of the 90 SNPs, 5 did not produce a positive signal. For 85 SNPs, 100% repeatiblity was proved by double genotyping of 13 randomly chosen boars. In addition, the accuracy of genotyping was verified in 2 sib-families by a Mendelian inheritance of 49-50 homozygous genotypes from sire to sons. Three genotype discrepancies were found (97% accuracy rate). All inaccurities were confirmed by an alternative method (sequencing and PCR-RFLP assays). Moreover, the exclusion power of the chip was evalueted by an SNP inheritance analysis of unrelated boars within each sib-family. In the validation step, 88 boars (13 Pietrain, 31 Landrace, 16 Large White, 8 Duroc, 7 Hampshire x Pietrain crosses, and 13 other hybrid lines) were screened to validate SNPs. Among the 85 selected SNPs, 12 were found to be monoallelic, the rest showing at least two genotypes for the entire population under study. The primary application of the SNiPORK chip is the simultaneous genotyping of dozens of SNPs to study gene interaction and consequently better understand the genetic background of pork yield and quality. The chip may prospectively be used for evolutionary studies, evaluation of genetic distances between wild and domestic pig breeds, traceability tests, as well as the starting point for developing a platform for identification and paternity analysis.     

High-throughput genotyping of single nucleotide polymorphisms using new biplex invader technology

The feasibility of large-scale genome-wide association studies of complex human disorders depends on the availability of accurate and efficient genotyping methods for single nucleotide polymorphisms (SNPs). We describe a new platform of the invader assay, a biplex assay, where both alleles are interrogated in a single reaction tube. The assay was evaluated on over 50 different SNPs, with over 20 SNPs genotyped in study cohorts of over 1500 individuals. We assessed the usefulness of the new platform in high-throughput genotyping and compared its accuracy to genotyping results obtained by the traditional monoplex invader assay, TaqMan genotyping and sequencing data. We present representative data for two SNPs in different genes (CD36 and protein tyrosine phosphatase 1β) from a study cohort comprising over 1500 individuals with high or low-normal blood pressure. In this high-throughput application, the biplex invader assay is very accurate, with an error rate of <0.3% and a failure rate of 1.64%. The set-up of the assay is highly automated, facilitating the processing of large numbers of samples simultaneously. We present new analysis tools for the assignment of genotypes that further improve genotyping success. The biplex invader assay with its automated set-up and analysis offers a new efficient high-throughput genotyping platform that is suitable for association studies in large study cohorts.     

SNP discovery and population structure analysis in Lassi and Marecha camel breeds using a genotyping by sequencing method

Pakistani camels have been classified socio-geographically into 20 breeds, but they have not yet been subjected to substantial selective pressures and the genetic basis for these breeds is not understood. However, it should be possible to distinguish them by use of molecular data. This study investigated the genetic diversity and population structure within and between two major Pakistani camel breeds, Marecha and Lassi. As no SNP array is currently available, we first identified 63 619 SNPs using a genotyping by sequencing approach. After quality control, a panel of 36 926 SNPs was used in the analysis. Population structure was investigated with a principal coordinate analysis as well as a cluster analysis using NetView , and multilocus heterozygosity analysis to explore between- and within-breed genetic variation. In addition, between-breed variation was explored using the fixation index, F_ST. We also compared relationship matrices computed using the VanRaden SNP-based method and a method developed specifically for genotyping by sequencing data. Among the two camel breeds, Lassi showed a lower level of genetic diversity whereas Marecha showed a higher level. As a genotyping platform has not yet been developed for the camel, the SNPs discovered in this study will be useful in future genetic studies in camels.     

Linkage disequilibrium over short physical distances measured in sheep using a high‐density SNP chip

The extent of linkage disequilibrium (LD) between genetic loci has implications for both association studies and the accuracy of genomic prediction. To characterise the persistence of LD in diverse sheep breeds, two SNP genotyping platforms were used. First, existing SNP genotypes from 63 breeds obtained using the ovine SNP50 BeadChip (49 034 loci) were used to estimate LD decay in populations with contrasting levels of genetic diversity. Given the paucity of marker pairs separated by short physical distances on the SNP50 BeadChip, genotyping was subsequently performed for four breeds using the recently developed ovine HD BeadChip that assays approximately 600 000 SNPs with an average genomic spacing of 5 kb. This facilitated a highly accurate estimate of LD over short genomic distances (<30 kb) and revealed LD varies considerably between sheep breeds. Further, sheep appear to contain generally lower levels of LD than do other domestic species, likely a reflection of aspects of their past population history.     

Evidence of admixture from haplotyping in an epidemiological study of UK Caucasian males: implications for association analyses

OBJECTIVE: Cohort and case-control genetic association studies offer the greatest power to detect small genotypic influences on disease phenotypes, relative to family-based designs. However, genetic subdivisions could confound studies involving unrelated individuals, but the topic has been little investigated. We examined geographical and interallelic association of SNP and microsatellite haplotypes of the Y chromosome, of regions of chromosome 11, and of autosomal SNP genotypes relevant to cardiovascular risk traits in a UK-wide epidemiological survey. RESULTS: We show evidence (p = 0.00001) of the Danelaw history of the UK, marked by a two-fold excess of a Viking Y haplotype in central England. We also found evidence for a (different) single-centre geographical over-representation of one haplotype, both for APOC3-A4-A5 and for IGF2. The basis of this remains obscure but neither reflect genotyping error nor correlate with the phenotypic associations by centre of these markers. A panel of SNPs relevant to cardiovascular risks traits showed neither association with geographical location nor with Y haplotypes. CONCLUSION: Combinations of Y haplotyping, autosomal haplotyping, and genome-wide SNP typing, taken together with phenotypic2 associations, should improve epidemiological recognition and interpretation of possible confounding by genetic subdivision.     

MilkProtChip--a microarray of SNPs in candidate genes associated with milk protein biosynthesis--development and validation

MilkProtChip is an oligonucleotide microarray based on the arrayed primer extension (APEX) technique, allowing genotyping of single nucleotide polymorphisms (SNPs) in genes of interest for bovine milk protein biosynthesis. APEX consists of a sequencing reaction primed by an oligonucleotide anchored with its 5'end to a glass slide and terminating one nucleotide before the polymorphic site. The extension with one fluorescently labeled dideoxy nucleotide complementary to the template reveals the polymorphism. A total of 75 SNPs were selected among those associated directly or potentially with milk protein content. Among the 75 SNPs, 4 did not produce a positive signal. Most of the remaining SNPs produced a signal for both strands, except for 4 (one strand). In the validation step, 12 Polish Holstein bulls, 1 Polish Red bull, 1 bison (Bison bonasus), 11 Jersey cows and 25 Polish Holstein cows were screened to validate SNPs. Among the 71 selected SNPs--26 were found monoallelic, the rest showing at least two genotypes for the entire population under study. All the animals were earlier genotyped for 2-5 SNPs by PCR-RFLP and PCR sequencing and all showed complete concordance with APEX genotyping. APEX reactions showed relatively high signal frequencies: more than 0.9, 0.9-0.8 and below 0.8, for 65, 4 and 2 DNA samples, respectively. The primary application of the MilkProtChip is the simultaneous genotyping of dozens of SNPs to reveal and clarify the genetic background of milk protein biosynthesis. The chip may possibly be used for dairy cattle identification and paternity analysis, evolutionary studies, the evaluation of genetic distances between wild and domestic cattle breeds and the domestication history of bovine species.