Skip to main content

Whole-genome haplotyping approaches and genomic medicine

Abstract

Genomic information reported as haplotypes rather than genotypes will be increasingly important for personalized medicine. Current technologies generate diploid sequence data that is rarely resolved into its constituent haplotypes. Furthermore, paradigms for thinking about genomic information are based on interpreting genotypes rather than haplotypes. Nevertheless, haplotypes have historically been useful in contexts ranging from population genetics to disease-gene mapping efforts. The main approaches for phasing genomic sequence data are molecular haplotyping, genetic haplotyping, and population-based inference. Long-read sequencing technologies are enabling longer molecular haplotypes, and decreases in the cost of whole-genome sequencing are enabling the sequencing of whole-chromosome genetic haplotypes. Hybrid approaches combining high-throughput short-read assembly with strategic approaches that enable physical or virtual binning of reads into haplotypes are enabling multi-gene haplotypes to be generated from single individuals. These techniques can be further combined with genetic and population approaches. Here, we review advances in whole-genome haplotyping approaches and discuss the importance of haplotypes for genomic medicine. Clinical applications include diagnosis by recognition of compound heterozygosity and by phasing regulatory variation to coding variation. Haplotypes, which are more specific than less complex variants such as single nucleotide variants, also have applications in prognostics and diagnostics, in the analysis of tumors, and in typing tissue for transplantation. Future advances will include technological innovations, the application of standard metrics for evaluating haplotype quality, and the development of databases that link haplotypes to disease.

Introduction

Technological progress has enabled the routine resequencing of human genomes. These genomes include rare variants at high frequency [1],[2] that are the result of exponential human population growth over the past hundred generations [3]. These variants can affect single nucleotides or larger genomic ranges by substitution, insertion, deletion, or change in copy number. Combinations of variants are present in cis on the same physical molecule or in trans on homologous chromosomes. This set of cis and trans relationships between the variants, known as the phase of the variants, affects the interpretation and implications of the relationships between genotypes and phenotypes, including disease phenotypes [4]-[6]. To simplify the discussion, we define a haplotype in general terms as a contiguous subset of the information contained in a molecule of DNA (Box 1). An example of a haplotype by this definition, grounded on molecular observation, is the actual sequence inherited from one parent and spanning one or more genes in a specific genomic region of interest. A corollary of this definition is that the longest possible haplotype, the `chromosome haplotype, is the full sequence of a chromosome that an individual inherited from one parent. Haplotypes have a number of important roles and applications that are listed in Box 2.

Analysis of haplotypes falls generally into three categories: 1) elucidating the `haplotype block structure of the genome, 2) employing haplotypes as genetic markers, and 3) understanding haplotypes as functional units [7]. As the number of observable genetic variants increases, so does the number of observable haplotypes. This increase in observed variants is largely the result of rare variants assayed by whole-genome sequencing (WGS). As a result, there are now many observed haplotypes with frequencies too small to estimate using population-based inference methods. Recent technological advances have enabled the determination of haplotypes by direct observation of molecular data and from genetic information, and with decreased reliance on population-based statistical estimation.

Historically, it has been difficult to distinguish between the homologous haplotypes of autosomes inherited from each of the parents. For that reason, allele information from each pair of autosomal chromosomes is typically comingled into one sequence of information: the unphased genotype sequence. Phasing (or haplotyping) describes the process of determining haplotypes from the genotype data. Until recently, cost, lack of data, and computational intractability has limited the availability of phased whole-genome haplotypes [4].

There are three basic methods for phasing: molecular, genetic, and population analysis. Molecular haplotyping is rooted in the analysis of single molecules (Figures1 and 2, Table1). If the molecule haplotyped is shorter than a chromosome, molecular haplotyping can be followed by haplotype assembly. Increasingly clever methods are being deployed to exploit high-throughput parallelization, combining data from measurements of many single molecules to build longer haplotypes. Genetic haplotyping requires pedigrees and can yield chromosome-length haplotypes [8]. Population haplotyping requires models of the population structure of haplotypes and can only phase common variations. These three approaches can also be combined to create hybrid strategies. As a general rule, these methods can be used to phase any combination of single nucleotide variants (SNVs; commonly called single nucleotide polymorphisms (SNPs) when they are frequent enough in the population), insertions and deletions (indels), and copy number variants (CNVs). SNVs and short indels are typically easier to phase because they can be observed within individual sequence reads. Larger variants, such as CNVs, are typically assessed using genotyping arrays [9],[10].

Figure 1
figure1

Overview of the three methods for phasing whole-genome sequence data. Phasing is achieved by (1) molecular methods, (2) genetic analysis, or (3) population inference. Molecular methods focus on individual samples and involve either (a) processing genomic DNA prior to sequencing or (b) leveraging the single-molecule characteristic of each physical read. Genetic analysis and population inference process sequenced genomes from pedigrees and population cohorts, respectively.

Figure 2
figure2

The properties of phasing methods. The level of confidence in phasing and the achievable range of phased sequence length both vary depending on the method used. Molecular methods provide direct observations from single molecules and therefore the level of confidence in the results is high. The phased sequence length that can be achieved by these methods has a wide range, which depends on the method employed. Molecular observations can be assembled into haplotypes (dashed arrows), adding moderately to the range of phased sequence length, and potentially introducing inference error. Genetic analyses infer phase by leveraging the property of Mendelian segregation and can phase entire chromosomes. Population inference methods are probabilistic and limited to the generation of short-range haplotype blocks.

Table 1 Overview of whole-genome haplotyping methods

In this review, we describe in detail the three main methodologies for phasing variants and their integration into combination strategies, and we provide quality metrics (Box 3). Finally, we provide an overview of the applications of whole-genome haplotyping in genomic medicine.

Molecular haplotyping

Molecular haplotyping involves the direct observation of alleles on a single molecule. These molecules are often single sequence reads, ranging in size from tens of bases to thousands of bases. When two variants are observed in the same physical read, or in paired reads derived from the same molecule, they can be directly phased. Therefore, same-read molecular haplotyping gains power with sequencing read length, and the major source of error is the sequencing error rate. Although often overlooked because of its simplicity, sequencing is the most common form of molecular haplotyping. Other forms of molecular haplotyping include restriction fragment analysis, optical mapping, and coded fluorescence hybridization approaches [11].

Long-range binning can be achieved by physical separation of the two haploid genomes prior to sequencing. Binning methods are able to resolve private and rare haplotypes and can be used to generate personalized genome-resolved haplotypes. Sequencing isolated sperm cell genomes [12] is one simple approach, but applicable only to males. Chromosome isolation methods do not require sequencing coverage to the depth needed to resolve possibly heterozygous positions [10]. Whole-chromosome sorting methods include microdissection, fluorescence-activated cell sorting (FACS) and microfluidics. Chromosomes are individually tagged or separated into pools that tend to contain at most one copy of a chromosome. These are genotyped or sequenced to generate whole-chromosome haplotypes. Microdissection involves arresting cells in metaphase and spreading the chromosomes to isolate them [13]. FACS separates individual chromosomes, which are then amplified and tagged before sequencing [14]. The `direct deterministic phasing method uses microfluidic devices to capture a single metaphase cell and partition each intact chromosome [10].

Semiconductor-based nanotechnologies are being applied to assay single DNA molecules, deriving very long-range haplotype information. NabSys (Providence, RI, USA) tags DNA molecules with probes that are specific to particular chromosomal locations and passes single molecules of DNA with bound tags through solid-state nanodetectors to identify the locations of bound tags [15]. BioNano Genomics (San Diego, CA, USA) labels DNA using nicking endonucleases and fluorescently labeled nucleotides, then visualizes single molecules in linearized nanochannels [16]. Both technologies yield de novo genome-wide maps, informing structural variation and haplotypes. Nevertheless, no current technology captures all variants in the genome; for example, most are unable to assay trinucleotide repeats. These technologies are changing rapidly, and Buermans and den Dunnen [17] have provided a recent review of the types of variants assayed by some of these technologies. The principles and methods of haplotyping described here will apply even as methods change. In some cases, combining a technology that assays large variants (for example, BioNano) with one that assays SNPs (for example, pairwise end sequencing) may best address a particular need.

Haplotype assembly

Any set of two or more overlapping haplotypes can be assembled into a single haplotype. Typically, after generation of many individual molecular haplotypes, sequence assembly is used to identify overlapping sequences and thus to infer a longer haplotype [18]-[20]. The haplotypes being assembled may be derived from heterogeneous data sources but haplotype assembly is most commonly based on a set of molecular haplotypes [21]-[23], so we will discuss this prior to the discussion of genetic and population-inferred haplotypes.

Assembly of molecular sequences from fragments predates the ability to sequence DNA. Assembly was originally employed for determining the sequence of proteins [24]. Before the Human Genome Project (HGP), genome maps were assembled from restriction-fragment haplotypes [25]. During the HGP, haplotype reconstruction relied on the assembly of matched-end sequences of clones. As the HGP wound down, for economy of scale, there was a general shift away from long-read towards short-read sequencing. This shift increased the difficulty of haplotype assembly directly from shotgun reads, and resulted in a revival in algorithms for haplotype assembly. Lancia et al. [26] describe the `SNPs haplotyping problem by looking at the fundamental constraint shared by the group of algorithms that solve this problem: that all sequence reads must be partitioned into exactly two bins. These algorithms generally allow for sets of raw reads to be constrained to co-occur in the same bin. Such constraints arise either from paired end data or from pooling strategies. Clever experimental designs have maximized the utility of these constraints, particularly those that use statistical or molecular techniques to bin reads from a particular haplotype.

In 2007, Levy et al.[27] used single sequence reads together with some mate pairs to build long-range haplotypes for an individual genome, with haplotypes reaching several hundred kilobases. In 2009, McKernan et al. [28] used a ligation-based sequencing assay to phase a single genome physically into blocks averaging several kilobases. In 2011, Kitzman et al. [29] produced 350kb haplotype blocks by subpooling a fosmid library. Suk et al. [30] also used fosmid pool-based sequencing to assemble variants into haplotypes of approximately 1Mb, up to a maximum length of 6.3Mb; fosmids were tiled into contiguous molecular haplotype sequences based on sequence overlaps [31]. In 2013, Kaper et al. [32] also used a dilution, amplification and sequencing approach to compile haplotypes of several hundred kilobases in length.

Extreme dilution of genomic DNA can generate long-range haplotypes without requiring the sorting of metaphase chromosomes or cloning. These methods recreate, with twists, the basic method used to sequence the human genome: local haplotypes (in the order of tens of kilobases) are first carefully sequenced and then strung together by aligning overlaps. Dilution methodologies allow long fragments to be shotgun sequenced with short reads [18]. If these long fragments overlap with a sufficient fingerprint [33], then haplotypes up to 1Mb may be achieved by chromosomal walking [34]. The number of DNA molecules in a pool is small enough that there is little chance that repeated or duplicate sequences will occur within a pool. Therefore, DNA dilution methods simplify both de novo assembly and mapping reads to a reference genome. Nevertheless, these methods can be confounded by the local presence of repetitive sequences. Commercialization of dilution methodologies now includes Complete Genomics `long fragment read (LFR) [35] and Illuminas Moleculo technology [36]. For LFR, long parental DNA fragments are separated into distinct pools and sequenced using pairwise end sequencing. Moleculo implements statistically aided, long-read haplotyping (SLRH) by further phasing initial contigs with population information using the Prism software (Table2).

Table 2 Summary of selected software available for whole-genome haplotyping

Several algorithms exist to assemble reads into haplotypes (Table2). HASH (haplotype assembly for single human) uses a Markov chain Monte Carlo (MCMC) algorithm and graph partitioning approach to assemble haplotypes given a list of heterozygous variants and a set of shotgun sequence reads mapped to a reference genome assembly [21]. HapCut uses the overlapping structure of the fragment matrix and max-cut computations to find the optimum minimum error correction (MEC) solution for haplotype assembly [22]. There are many other sequence assembly algorithms, reviewed elsewhere [54],[55]. Duitama et al. [31] reviewed eight algorithms for the `SNPs haplotyping problem with binned reads as input. They concluded that, under a reasonable compromise between accuracy, completeness, and computational resources, ReFHap (Reliable and Fast Haplotyping) [37] yields the best results for a low-coverage fosmid pooling approach, which they term single individual haplotyping (SIH). More recent algorithms claim improvements on parameters such as speed and accuracy (for example, H-BOP[38]) or focus on improving performance in the presence of high error rates [39],[56]. MixSH shows good performance as evaluated by pair consistency, a version of a metric described in Box 3 [39].

The process of assembly may introduce phase errors at the joins between component haplotypes, and so should best be done when the overlaps between fragments can be inferred with high confidence. Such confidence can be gained both by identification of unique overlapping fingerprints or by physical separation of the original molecules. Haplotype assembly has worked very well when the underlying haplotypes are long, such as those determined by sequencing a clonal source such as a cosmid or bacterial artificial chromosome (BAC) [14],[29]. We therefore expect to see increasing development of technologies that generate sequence reads in the range of many thousands of bases to facilitate haplotype assembly. These long sequences will be generated by strobe sequencing, nanopore sequencing, and perhaps other technologies [57],[58].

The existence of chromosome territories in the nucleus can also be exploited for long-range haplotyping. In an innovative approach, pairs of reads that are likely to come from the same haplotype are generated by cross-linking chromatin sites that are potentially distant along a chromosome but spatially close within the nucleus. This technique is known as `Hi-C, and was simultaneously exploited by three different groups for sequence assembly [59]-[61]. Selvaraj et al. [62] focused on haplotyping using Hi-C (which they term `HaploSeq), and in their initial report using low coverage sequencing they phased approximately 81% of sequenced alleles.

Disparate sources of haplotyping information and markers can also be assembled. For example, the 1000 Genomes Project Consortium recently produced an integrated haplotype map of SNPs, small indels and larger deletions derived from SNP arrays or from exome and whole-genome sequencing [63].

Genetic analysis

The principles of Mendelian segregation of alleles in pedigrees can be used to deduce the phasing of variants observed in ordered genotypes. At the simplest level of a family trio (both parents and one child), very simple rules indicate which alleles in the child were inherited from each parent, thus largely separating the two haplotypes in the child. The remaining (not inherited) parental haplotypes can then be reconstructed using a simple exclusion rule. As the locations of recombinations are not known, the inferred parental haplotypes will have a phase error at each recombination. These low-frequency phase errors (Box 1) will have little effect on short-range haplotypes but will scramble chromosomal haplotypes.

In the context of a family quartet (two full siblings and their parents), whole-genome sequences from high-throughput paired-end short-reads can generate complete chromosomal haplotypes for all family members [8],[64]. The method can be extended to larger pedigrees by tiling or MCMC approaches [65]. Tiling can accumulate small errors with each tile, and so MCMC and similar approaches are likely to be the best methods for pedigrees spanning more than four generations. Haploscribe is a suite of software scripts that phase whole-genome data across entire chromosomes by genetic analysis [8]. Haploscribe implements a parsimony approach to generate meiosis-indicator (inheritance state) vectors and uses a hidden Markov model (HMM) to deduce haplotypes spanning entire chromosomes. These haplotypes are nearly 100% accurate and potentially suitable for medical diagnostics.

The rule-based nature of genetic phasing has a useful property: some family genotypes are not consistent with the expected patterns of inheritance, and are thus highlighted as probable sequencing errors or, rarely, as de novo mutations [64]. Mendelian inheritance errors (MIEs) are sites in which the genotype of a child is inconsistent with inheritance from one or both parents. In state consistency errors (SCEs), the genotype of each child is consistent with both parents but the combination of offspring genotypes is inconsistent with the prevailing inheritance state around that locus, as determined from neighboring sites.

Genetic analysis enables the phasing of rare alleles that cannot otherwise be accomplished by reference to population-based data. Phasing information obtained through the sequencing of the genomes of family members maps recombinations and inheritance states at high resolution, highlighting the regions of the genome where causal variants segregate. The resulting haplotypes are highly accurate and complete. Nevertheless, genetic analysis cannot phase positions in which all family members are heterozygous. Furthermore, it is not always feasible to recruit the required participants for family-based studies. In the absence of a family context, molecular haplotyping is an excellent choice because it does not require DNA samples from other family members. We predict that, in the next decade, molecular haplotyping will largely supplant the need for genetic analysis.

Dewey et al. [66] employed family inheritance-state analysis to control sequencing error and inform haplotype phasing to quantify genome-wide compound heterozygosity from high-throughput sequencing data. To define the inheritance states of neighboring SNVs in the family quartet, Dewey et al. first used a heuristic algorithm that binned allele assortments, followed by a HMM in which the hidden states corresponded to the four possible inheritance states in the quartet and the two error states described by Roach et al. [64]. A combination of pedigree data and statistical phasing based on inheritance state analysis was then used to infer phase for the majority of positions in each childs genome. For uniformly heterozygous positions, the minor allele was assigned to the paternal and maternal chromosome scaffolds using pair-wise pre-computed population linkage disequilibrium (LD) data from the SNP Annotation and Proxy Search (SNAP) database [67]. These algorithms successfully determined genome-wide, long-ranging phased haplotypes in the family quartet. Phased variant data were also used to determine parental contribution to each childs disease risk in the context of thrombophilia.

Population inference

Population analysis leverages shared ancestry information to infer the most likely phasing of variants. The reference population can range from the very large (for example, the global human population), to the narrowly defined (for example, an isolated community). Because population relationships may be distant or cryptic, methodologies for population analysis are statistical and not deterministic. Also, because many more meioses separate all of the genomes in a large population, the length of haplotypes determined by population analysis is typically limited to thousands or tens of thousands of bases. Population inference methods work well on genotyping panels, which are compilations of common SNPs. As marker density increases, brute-force algorithms become less tractable, and algorithms such as those based on HMMs are employed [68]. Discerning private and rare haplotypes by population-based methods is highly challenging. Population analysis cannot phase de novo mutations, rare variants or structural variants. If a rare variant is assigned to a haplotype by other methods, however, its presence on a haplotype determined by common SNPs can be probabilistically inferred [69].

Parsimony approaches such as Clarks algorithm [70] attempt to find the minimum number of unique haplotypes in a data set. The accuracy of this method depends on the assumption that markers are tightly linked and largely assignable to common haplotypes. Therefore, such algorithms over-predict common haplotypes. Coalescent-based methods and HMMs are also commonly employed to model population haplotype frequencies. The software programs PHASE[45], fastPHASE[41], MaCH[44], IMPUTE2[43], SHAPEIT[71], and Beagle[40] implement such methods (Table2). These methods estimate parameters iteratively, so they work well with a small number of genetic markers residing on a short haplotype block.

SHAPEIT (segmented haplotype estimation and imputation tool) scales linearly with the number of variants, samples and conditional haplotypes used in each iteration and can be run efficiently on whole chromosomes [46],[71]: it was applied to generate a haplotype map of 38 million SNPs for phase 1 of the 1000 Genomes Project [63],[69]. This program is versatile for population-based studies as it is able to handle data from combinations of unrelated individuals, duos and trios. SHAPEIT2[51] adds the ability to incorporate molecular information from sequence reads, incorporating calls and base-quality scores in a probabilistic model. It works best for high-coverage sequence. OConnell et al. [72] have incorporated SHAPEIT2 into a general haplotyping workflow that can also detect state consistency errors in pedigrees [64].

Wang et al. [48] developed a population imputation pipeline, SNPTools, to phase low-coverage data obtained from phase 1 of the 1000 Genomes Project. SNPTools statistically models sequence data, scores polymorphic loci, and generates genotype likelihoods using a Binary sequence map (BAM)-specific Binomial Mixture Model (BBMM). The genotype likelihoods can then be integrated into SNPTools imputation algorithm or other algorithms such as Beagle to estimate phased genotypes and haplotypes. SNPTools haplotype imputation algorithm employs a four-state constrained HMM sampling scheme that assumes that the individual haplotype is a mosaic of the four parental haplotypes. WinHAP estimates multi-SNV haplotypes from large-scale genotype data [49]. This software simplifies the 2SNP algorithm, using pairs of heterozygous SNVs to generate initial haplotypes and subsequently to construct a linear tree that makes it possible to infer a solution for the haplotype phase [73]. These haplotypes are then improved by applying scalable sliding windows. Last, parsimony is used to iteratively restrict the number of haplotypes.

The accuracy of population haplotyping can be improved by modeling population substructure and detecting cryptic relatedness. Such issues may be overcome by exploiting algorithms originally conceived for identical-by-descent (IBD) detection [74]. Programs such as fastIBD[40] and GERMLINE[18] leverage population level IBD to define haplotypes [75]. The extent of shared IBD between a pair of individuals depends on the number of generations since their last common ancestor as recombination and mutation increase genetic diversity. GERMLINE directly matches portions of haplotypes between individuals from phased genotype data. FastIBD uses a HMM approach for IBD detection of shared haplotypes from unphased genotype data. IBD segments are identified by modeling shared haplotype frequencies that account for background levels of LD.

Most of the available algorithms for population inference of haplotypes from WGS require careful balancing of computational speed and accuracy. They also rely on the availability of well-characterized, population-matched reference datasets [76]; these need to be large enough to sample rare variants. Population-based phasing methods are probabilistic, limited to generation of short haplotype blocks, and will incorrectly phase rare combinations of variants, exactly those combinations most likely to be medically important. Moreover, haplotypes derived from algorithms that include population inference are likely to have an error rate that is unacceptably high for medical purposes.

If an individual is a member of a completely characterized isolated population, the accuracy of population-based haplotypes can be very high. Such haplotyping has been demonstrated by Kong et al. [77] for the Icelandic population. Use of such databases in combination with methods to phase de novo mutations and haplotypes resulting from recent recombinations could both permit increased haplotype quality and reduce the need for genetic and molecular haplotyping in patients from these populations.

Combination strategies

Combinations of molecular, genetic and population-based methods may work better than any single approach, by combining strengths and minimizing weaknesses (Table2). HARSH evolved from Hap-seqX, combining haplotype assembly with a reference population dataset to predict haplotypes from WGS data [50],[78]. Prism, mentioned earlier, is another recent hybrid algorithm [36]. PPHS (Perfect Phylogeny Haplotypes from Sequencing) is another combination approach that combines population and molecular analysis by using raw sequence data to construct a phylogenetic tree for every short region [52]. The phylogeny model assumes that there are no recurrent mutations or recombination within each short sequence region. For each set of SNPs in a given window, the algorithm reconstructs a local phylogenetic tree by expectation maximization and assigns haplotypes to individuals. The results for each window are then stitched together using the output of Beagle as a guide. Combination strategies such as these may increase the accuracy of population inference methods by leveraging the information provided by sequence data, or may supplement genetic analyses with population data, as described by Dewey et al.[66]. A combination of genetic, physical and population-based approaches in a quartet yielded complete genome phasing, including phasing of 99.8% of fully heterozygous variants [53].

Clinical applications of whole-genome haplotyping

Local haplotyping has been and will remain important for genomic diagnostics. The immediate impact of whole-genome haplotyping will be to provide all local haplotypes. Local haplotypes are well known for the major histocompatibility complex (MHC) and several other loci, including the ApoE4 haplotype of the ApoE locus. MHC haplotypes are important for predicting graft compatibility and for prediction of the risks and protectivity of many phenotypes, notably type 1 diabetes [79]. In many cases, the causative variant is not known, and the observed haplotype serves as a proxy for assaying the unknown single variant that lies in or is linked to that haplotype. In other cases, such as ApoE4, multiple coding variants must lie on the same haplotype within a single coding sequence in order to effect a particular phenotypic change. Family-based haplotyping to identify compound heterozygosity as a cause of recessive Mendelian disease is also fairly routine. Fetal and newborn diagnostics will also benefit from haplotyping. Spearheading such an approach in 2012, Kitzman et al. [80] inferred haplotypes of a human fetus using haplotype-resolved genome sequencing of the mother, shotgun sequencing of the father, and deep sequencing of maternal plasma.

Pathogenic rare variants will be a significant source of concern when practicing genomic medicine. Thus, an important clinical application of haplotypes will occur at the largely unseen analysis stage - in improving variant calling and avoiding false alarms. Already, software tools such as Platypus (www.well.ox.ac.uk/platypus) are being developed to produce improved base calling as informed by haplotypes [63].

As knowledge and methods improve, understanding the functional interactions between regulatory elements and coding regions will permit medical decision-making that is based not only on the predicted effects of variants on the function of a protein, but also on combining separate predictions of (a) the functions of the two proteins produced by the two alleles of the encoding gene, and (b) the effects of the two cis-regulators of these two proteins. For example, if one of the cis-regulators markedly increases expression while the other decreases expression, and one protein is defective while the other is normal, then one combination of cis-regulators with the protein-coding alleles will produce wellness whereas the other combination will produce disease [4].

Conclusions and future directions

High-throughput short-read sequencing enabled rapid advances following the HGP. Unfortunately, haplotyping got left by the wayside, as the long reads characteristic of the HGP gave way to cheaper short reads. Now a combination of new strategies and new technologies is enabling the determination of personal haplotypes that will soon be economical for more routine medical use. The new strategies that we have discussed enable the use of cheap short reads for inferring longer haplotypes, typically by physically or computationally placing these reads into haplotype bins. Some new technologies, such as Hi-C, facilitate this binning process, whereas other new technologies will enable the generation of cheaper long reads.

Considering the garbage-in-garbage-out principle, and that most current algorithms perform near perfectly on error-free data, improving sequencing error rate is probably the most critical factor for improving haplotypes [81]. In other words, to improve haplotypes for use in genomic medicine, a focus on phasing algorithms and methodologies is not necessarily the greatest requirement, but rather a focus on improving the input data. More phase errors can arise with whole-genome data than with genotyping chip data. The SNPs included in genotyping chips tend to be selected for Hardy-Weinberg equilibrium, and so any SNP with a heterozygote frequency that is unexpected in relation to the allele frequencies is excluded. Such pre-selection is not done for WGS data. The ability to phase WGS data can be confounded by reference sequence errors, reference gaps and centromeres, and long interspersed nuclear elements (LINEs) [8]. Methods are needed for filtering out these regions or for handling them in a probabilistic framework with appropriate confidence statistics. In conjunction with improvements in sequence quality, the generation of long sequence reads (of 10,000 bases or longer) is another key factor for haplotype improvement [18],[19],[82]. Reduction of sequencing error will have the greatest impact on high-frequency switch error, whereas improvements in read length will have a greater impact on low-frequency switch error (Boxes 1 and 3, and Figure3). Because of the importance of compound heterozygote analysis and within-gene phasing for medical applications, high-frequency switch errors must be minimized. Commoditization of haplotyping will also be necessary, and as this occurs, the costs of the various approaches will become less opaque. Other commodity technologies, such as sequencing, are best performed in high-throughput operations because such facilities offer a concentration of expertise, economy of scale, standard operating procedures, and rigorous quality control. Clinical haplotype databases will need to be developed in parallel with haplotype commoditization, much like the ClinVar database for individual variants associated with disease [83].

Figure 3
figure3

Switch errors and quality metrics. We present four scenarios involving low- and high-frequency switch errors and the resulting values for two quality metrics: total switch rate, and the pairwise metric 2*cse-1, which we term the `phase accuracy (Box 3). The latter was computed on the assumption that each local `back and forth switch error affects 1% of the markers in the chromosome. The true haplotypes (maternal (M) and paternal (P)) are shown at the top. In Example 1, there is a single phase-error switch at `a; variants on opposite sides of the error are incorrectly phased (denoted by the `x on the curved arrow). This situation might arise if a parents haplotype is inferred by subtraction of a childs haplotype, missing a meiotic recombination. Example 1 thus has a low-frequency phase error but no high-frequency errors. In Example 2, there are many high-frequency errors, but no low-frequency errors. Most variant pairs are separated by even numbers of switch errors, and are thus properly phased (curved arrows). This error pattern might arise from a long-read technology such as strobe sequencing. In Example 3, there are three switch errors, one at `a and two at `b. This haplotype could have arisen from a false haplotype assembly join at `a and a sequence error of a single base at `b. Example 3 thus has a mix of low- and high-frequency errors. In terms of the pairwise metric, the single switch in Example 1 most strongly affects the long-range haplotype quality, while the several localized switch errors (each pair affecting just one or a few markers) degrade haplotype quality only modestly.

It is now routine in medical genetics to consider compound heterozygosity in identifying disease risks and causes in patients. Typically, this search is either carried out by genetic haplotyping if the sequences of parents are available or is achieved by considering all possible haplotypes of detrimental variations within a gene. We identified compound heterozygosity causing Miller syndrome in our analysis of the first whole-genome sequenced family [64], but such analyses have been routine for years in the analysis of candidate genes, such as those for cystic fibrosis and breast cancer risk. Numerous examples include the identification of compound heterozygous causes of diseases, including the gene that encodes protein C in cerebral palsy [84], Charcot-Marie-Tooth neuropathy [85], and the gene encoding lysyl-tRNA synthetase in peripheral neuropathy [86]. Currently, much clinical screening for compound heterozygosity is done with exome sequencing, but we predict a shift towards WGS as costs drop. As the understanding and annotation of regulatory variants continues to improve, we will see an increasing number of reports of cis-acting regulatory elements that alter gene expression and cause disease. Examples that have already been reported include a mutation in a RET enhancer that increases Hirschsprung disease risk [87],[88] as well as mutations that affect thalassemia, hemophilia, hypercholesterolemia, and pyruvate kinase deficiency [89]. Increasingly, the phase of these regulatory elements with respect to the coding variants will be part of routine diagnostics. For example, Ludlow et al. [90] described a patient with a mutation in the promoter of one allele of GP1BB (encoding platelet glycoprotein Ib beta) and a deletion of the other allele, which together resulted in Bernard-Soulier syndrome [90].

DNA diagnostics and prognostics that have clinical applications in oncology are expanding rapidly. For example, particular haplotypes of GNMT (encoding glycine N-methyltransferase) differently predispose individuals to prostate cancer [91]. In many oncological applications, genetic and population phasing is of little value because of the large number of somatic mutations that may be present in tumor cells. Molecular phasing will therefore be the primary tool in this area, and algorithms that allow for multiple ploidy states will be important in handling the complexities of tumor genomes; currently most haplotype assembly algorithms assume diploidy. The MHC/HLA (human leukocyte antigen) locus is the most important haplotype influencing disease; and an understanding of the value of MHC haplotypes is therefore nothing new [92],[93]. It has traditionally been difficult to use molecular techniques that avoid low-frequency switch errors between genes of the MHC. Applications of some of the new long-range haplotyping techniques, particularly those capable of de novo assembly of regions of personal genomes within the MHC that are not in the reference genome, are likely to allow better utility of MHC typing for research, prognostics, diagnostics, and tissue transplantation [93]. The genomic medicine of the future will rely on accurately mining patient sequence data to identify disease, wellness and actionable genes [6]. Genomics must move beyond simple single allelic and genotypic tests of association and familial-segregation to explain phenotypes. At the simplest level, whether two detrimental variants that affect `the same gene lie in cis or in trans may spell the difference between a healthy carrier and a diseased compound heterozygote. The paradigm of medical understanding must be shifted from `the function of a gene in an individual to `the functions of each allele of each gene in an individual. To achieve this, we must transform the conceptualization of the genome in the minds of both clinicians and researchers from one that contains 22 autosomes and two sex chromosomes to one that contains 44 autosomes and two sex chromosomes, each with its own haplotype. Individual genome sequencing is being applied at all stages of life, from preimplantation, prenatal and neonatal diagnosis, to `no phenotype personalized genomics [94]. Whole-genome haplotypes will improve the precision of personalized predictive, preventive and participatory medicine.

Abbreviations

BAC:

Bacterial artificial chromosome

BAM:

Binary alignment map

BBMM:

BAM-specific Binomial Mixture Model

CNV:

Copy number variant

FACS:

Fluorescence-activated cell sorting

HASH :

Haplotype assembly for single human

HGP:

Human Genome Project

HLA:

Human leukocyte antigen

HMM:

Hidden Markov model

IBD:

Identical-by-descent

LD:

Linkage disequilibrium

LFR:

Long fragment read

LINE:

Long interspersed nuclear element

MCMC:

Markov chain Monte Carlo

MEC:

Minimum error correction

MHC:

Major histocompatibility complex

MIE:

Mendelian inheritance error

PPHS:

Perfect Phylogeny Haplotypes from Sequencing

SCE:

State consistency error

SIH:

Single individual haplotyping

SLRH:

Statistically aided, long-read haplotyping

SNAP:

SNP Annotation and Proxy Search

SNP:

Single nucleotide polymorphism

SNV:

Single nucleotide variant

WGS:

Whole-genome sequencing

References

  1. 1.

    Tennessen JA, Bigham AW, OConnor TD, Fu W, Kenny EE, Gravel S, McGee S, Do R, Liu X, Jun G, Kang HM, Jordan D, Leal SM, Gabriel S, Rieder MJ, Abecasis G, Altshuler D, Nickerson DA, Boerwinkle E, Sunyaev S, Bustamante CD, Bamshad MJ, Akey JM: Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. 2012, 337: 64-69.

    PubMed  PubMed Central  Google Scholar 

  2. 2.

    Gravel S, Henn BM, Gutenkunst RN, Indap AR, Marth GT, Clark AG, Yu F, Gibbs RA, Bustamante CD: Demographic history and rare allele sharing among human populations. Proc Natl Acad Sci U S A. 2011, 108: 11983-11988.

    PubMed  PubMed Central  Google Scholar 

  3. 3.

    Keinan A, Clark AG: Recent explosive human population growth has resulted in an excess of rare genetic variants. Science. 2012, 336: 740-743.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Tewhey R, Bansal V, Torkamani A, Topol EJ, Schork NJ: The importance of phase information for human genomics. Nat Rev Genet. 2011, 12: 215-223.

    PubMed  PubMed Central  Google Scholar 

  5. 5.

    Muers M: Genomics: no half measures for haplotypes. Nat Rev Genet. 2011, 12: 77-

    PubMed  Google Scholar 

  6. 6.

    Tian Q, Hood L: Systems approaches to biology and disease enable translational systems medicine. Genomics Proteomics Bioinformatics. 2012, 10: 181-185.

    PubMed  Google Scholar 

  7. 7.

    Hoehe MR: Haplotypes and the systematic analysis of genetic variation in genes and genomes. Pharmacogenomics. 2003, 4: 547-570.

    PubMed  Google Scholar 

  8. 8.

    Roach JC, Glusman G, Hubley R, Montsaroff SZ, Holloway AK, Mauldin DE, Srivastava D, Garg V, Pollard KS, Galas DJ, Hood L, Smit AF: Chromosomal haplotypes by genetic phasing of human families. Am J Hum Genet. 2011, 89: 382-397.

    PubMed  PubMed Central  Google Scholar 

  9. 9.

    McCarroll SA, Kuruvilla FG, Korn JM, Cawley S, Nemesh J, Wysoker A, Shapero MH, de Bakker PIW, Maller JB, Kirby A, Elliott AL, Parkin M, Hubbell E, Webster T, Mei R, Veitch J, Collins PJ, Handsaker R, Lincoln S, Nizzari M, Blume J, Jones KW, Rava R, Daly MJ, Gabriel SB, Altshuler D: Integrated detection and population-genetic analysis of SNPs and copy number variation. Nat Genet. 2008, 40: 1166-1174.

    PubMed  Google Scholar 

  10. 10.

    Fan HC, Wang J, Potanina A, Quake SR: Whole-genome molecular haplotyping of single cells. Nat Biotechnol. 2011, 29: 51-57.

    PubMed  PubMed Central  Google Scholar 

  11. 11.

    Geiss GK, Bumgarner RE, Birditt B, Dahl T, Dowidar N, Dunaway DL, Fell HP, Ferree S, George RD, Grogan T, James JJ, Maysuria M, Mitton JD, Oliveri P, Osborn JL, Peng T, Ratcliffe AL, Webster PJ, Davidson EH, Hood L, Dimitrov K: Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat Biotechnol. 2008, 26: 317-325.

    PubMed  Google Scholar 

  12. 12.

    Kirkness EF, Grindberg RV, Yee-Greenbaum J, Marshall CR, Scherer SW, Lasken RS, Venter JC: Sequencing of isolated sperm cells for direct haplotyping of a human genome. Genome Res. 2013, 23: 826-832.

    PubMed  PubMed Central  Google Scholar 

  13. 13.

    Ma L, Xiao Y, Huang H, Wang Q, Rao W, Feng Y, Zhang K, Song Q: Direct determination of molecular haplotypes by chromosome microdissection. Nat Methods. 2010, 7: 299-301.

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Yang H, Chen X, Wong WH: Completely phased genome sequencing through chromosome sorting. Proc Natl Acad Sci U S A. 2011, 108: 12-17.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Thompson JF, Oliver JS: Mapping and sequencing DNA using nanopores and nanodetectors. Electrophoresis. 2012, 33: 3429-3436.

    PubMed  Google Scholar 

  16. 16.

    Lam ET, Hastie A, Lin C, Ehrlich D, Das SK, Austin MD, Deshpande P, Cao H, Nagarajan N, Xiao M, Kwok P-Y: Genome mapping on nanochannel arrays for structural variation analysis and sequence assembly. Nat Biotechnol. 2012, 30: 771-776.

    PubMed  Google Scholar 

  17. 17.

    Buermans HPJ, den Dunnen JT: Next generation sequencing technology: advances and applications.Biochim Biophys Acta 2014,5: 464-473. doi:10.1016/j.bbadis.2014.06.01.

    Google Scholar 

  18. 18.

    Roach JC: Random subcloning. Genome Res. 1995, 5: 464-473.

    PubMed  Google Scholar 

  19. 19.

    Roach JC, Boysen C, Wang K, Hood L: Pairwise end sequencing: a unified approach to genomic mapping and sequencing. Genomics. 1995, 26: 345-353.

    PubMed  Google Scholar 

  20. 20.

    Roach JC, Thorsson V, Siegel AF: Parking strategies for genome sequencing. Genome Res. 2000, 10: 1020-1030.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Bansal V, Halpern AL, Axelrod N, Bafna V: An MCMC algorithm for haplotype assembly from whole-genome sequence data. Genome Res. 2008, 18: 1336-1346.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Bansal V, Bafna V: HapCUT: an efficient and accurate algorithm for the haplotype assembly problem. Bioinformatics. 2008, 24: i153-i159.

    PubMed  Google Scholar 

  23. 23.

    He D, Choi A, Pipatsrisawat K, Darwiche A, Eskin E: Optimal algorithms for haplotype assembly from whole-genome sequence data. Bioinformatics. 2010, 26: i183-i190.

    PubMed  PubMed Central  Google Scholar 

  24. 24.

    Tsugita A, Gish DT, Young J, Fraenkel-Conrat H, Knight CA, Stanley WM: The complete amino acid sequence of the protein of Tobacco Mosaic Virus. Proc Natl Acad Sci U S A. 1960, 46: 1463-1469.

    PubMed  PubMed Central  Google Scholar 

  25. 25.

    Olson MV, Dutchik JE, Graham MY, Brodeur GM, Helms C, Frank M, MacCollin M, Scheinman R, Frank T: Random-clone strategy for genomic restriction mapping in yeast. Proc Natl Acad Sci U S A. 1986, 83: 7826-7830.

    PubMed  PubMed Central  Google Scholar 

  26. 26.

    Lancia G, Bafna V, Istrail S, Lippert R, Schwartz R: SNPs problems, complexity and algorithms. In Algorithms - ESA 2001. Volume 2161. Edited by auf der Heide FM. Berlin, Heidelberg: Springer; 2001:182193.

  27. 27.

    Levy S, Sutton G, Ng PC, Feuk L, Halpern AL, Walenz BP, Axelrod N, Huang J, Kirkness EF, Denisov G, Lin Y, MacDonald JR, Pang AWC, Shago M, Stockwell TB, Tsiamouri A, Bafna V, Bansal V, Kravitz SA, Busam DA, Beeson KY, McIntosh TC, Remington KA, Abril JF, Gill J, Borman J, Rogers Y-H, Frazier ME, Scherer SW, Strausberg RL:et al, The diploid genome sequence of an individual human. PLoS Biol. 2007, 5: e254-

    PubMed  PubMed Central  Google Scholar 

  28. 28.

    McKernan KJ, Peckham HE, Costa GL, McLaughlin SF, Fu Y, Tsung EF, Clouser CR, Duncan C, Ichikawa JK, Lee CC, Zhang Z, Ranade SS, Dimalanta ET, Hyland FC, Sokolsky TD, Zhang L, Sheridan A, Fu H, Hendrickson CL, Li B, Kotler L, Stuart JR, Malek JA, Manning JM, Antipova AA, Perez DS, Moore MP, Hayashibara KC, Lyons MR, Beaudoin RE: et al, Sequence and structural variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding. Genome Res. 2009, 19: 1527-1541.

    PubMed  PubMed Central  Google Scholar 

  29. 29.

    Kitzman JO, Mackenzie AP, Adey A, Hiatt JB, Patwardhan RP, Sudmant PH, Ng SB, Alkan C, Qiu R, Eichler EE, Shendure J: Haplotype-resolved genome sequencing of a Gujarati Indian individual. Nat Biotechnol. 2011, 29: 59-63.

    PubMed  PubMed Central  Google Scholar 

  30. 30.

    Suk E-K, Schulz S, Nowick K, Duitama J, Peckham H, Lee C, McLaughlin S, Schreiber S, Palczewski S, Holloway DT, McEwen GK, Hoehe MR: A comprehensively molecular haplotype-resolved genome of a European individual. Genome Res. 2011, 21: 1672-1685.

    PubMed  PubMed Central  Google Scholar 

  31. 31.

    Duitama J, McEwen GK, Huebsch T, Palczewski S, Schulz S, Verstrepen K, Suk E-K, Hoehe MR: Fosmid-based whole genome haplotyping of a HapMap trio child: evaluation of single individual haplotyping techniques. Nucleic Acids Res. 2012, 40: 2041-2053.

    PubMed  PubMed Central  Google Scholar 

  32. 32.

    Kaper F, Swamy S, Klotzle B, Munchel S, Cottrell J, Bibikova M, Chuang H-Y, Kruglyak S, Ronaghi M, Eberle MA, Fan J-B: Whole-genome haplotyping by dilution, amplification, and sequencing. Proc Natl Acad Sci U S A. 2013, 110: 5552-5557.

    PubMed  PubMed Central  Google Scholar 

  33. 33.

    Siegel AF, Roach JC, van den Engh G: Expectation and variance of true and false fragment matches in DNA restriction mapping. J Comput Biol. 1998, 5: 101-111.

    PubMed  Google Scholar 

  34. 34.

    Yen PH, Davidson N: The gross anatomy of a tRNA gene cluster at region 42A of the D. melanogaster chromosome. Cell. 1980, 22: 137-148.

    PubMed  Google Scholar 

  35. 35.

    Peters BA, Kermani BG, Sparks AB, Alferov O, Hong P, Alexeev A, Jiang Y, Dahl F, Tang YT, Haas J, Robasky K, Zaranek AW, Lee J-H, Ball MP, Peterson JE, Perazich H, Yeung G, Liu J, Chen L, Kennemer MI, Pothuraju K, Konvicka K, Tsoupko-Sitnikov M, Pant KP, Ebert JC, Nilsen GB, Baccash J, Halpern AL, Church GM, Drmanac R: Accurate whole-genome sequencing and haplotyping from 10 to 20 human cells. Nature. 2012, 487: 190-195.

    PubMed  PubMed Central  Google Scholar 

  36. 36.

    Kuleshov V, Xie D, Chen R, Pushkarev D, Ma Z, Blauwkamp T, Kertesz M, Snyder M: Whole-genome haplotyping using long reads and statistical methods. Nat Biotechnol. 2014, 32: 261-266.

    PubMed  PubMed Central  Google Scholar 

  37. 37.

    Duitama J, Huebsch T, Mcewen G, Suk E, Hoehe MR: ReFHap: a reliable and fast algorithm for single individual haplotyping. Proceedings of the First ACM international Conference on Bioinformatics and Computational Biology: August 2-4, 2010. 2010, 160-169. ACM, Niagara Falls, New York

    Google Scholar 

  38. 38.

    Xie M, Wang J, Jiang T: A fast and accurate algorithm for single individual haplotyping. BMC Syst Biol. 2012, 6 (Suppl 2): S8-

    PubMed  PubMed Central  Google Scholar 

  39. 39.

    Matsumoto H, Kiryu H: MixSIH: a mixture model for single individual haplotyping. BMC Genomics. 2013, 14 (Suppl 2): S5-

    PubMed  PubMed Central  Google Scholar 

  40. 40.

    Browning SR, Browning BL: Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007, 81: 1084-1097.

    PubMed  PubMed Central  Google Scholar 

  41. 41.

    Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78: 629-644.

    PubMed  PubMed Central  Google Scholar 

  42. 42.

    Gusev A, Lowe JK, Stoffel M, Daly MJ, Altshuler D, Breslow JL, Friedman JM, Peer I: Whole population, genome-wide mapping of hidden relatedness. Genome Res. 2009, 19: 318-326.

    PubMed  PubMed Central  Google Scholar 

  43. 43.

    Howie BN, Donnelly P, Marchini J: A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009, 5: e1000529-

    PubMed  PubMed Central  Google Scholar 

  44. 44.

    Li Y, Willer CJ, Ding J, Scheet P, Abecasis GR: MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet Epidemiol. 2010, 34: 816-834.

    PubMed  PubMed Central  Google Scholar 

  45. 45.

    Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73: 1162-1169.

    PubMed  PubMed Central  Google Scholar 

  46. 46.

    Delaneau O, Zagury J-F, Marchini J: Improved whole-chromosome phasing for disease and population genetic studies. Nat Methods. 2013, 10: 5-6.

    PubMed  Google Scholar 

  47. 47.

    Delaneau O, Marchini J, Zagury JF: A linear complexity phasing method for thousands of genomes. Nat Methods. 2011, 9: 179-181.

    PubMed  Google Scholar 

  48. 48.

    Wang Y, Lu J, Yu J, Gibbs RA, Yu F: An integrative variant analysis pipeline for accurate genotype/haplotype inference in population NGS data. Genome Res. 2013, 23: 833-842.

    PubMed  PubMed Central  Google Scholar 

  49. 49.

    Cheng W, Zhou F, Nie P, Xu Y: WinHAP: an efficient haplotype phasing algorithm based on scalable sliding windows. PLoS One. 2012, 7: e43163-

    PubMed  PubMed Central  Google Scholar 

  50. 50.

    Yang W-Y, Hormozdiari F, Wang Z, He D, Pasaniuc B, Eskin E: Leveraging reads that span multiple single nucleotide polymorphisms for haplotype inference from sequencing data. Bioinformatics. 2013, 29: 2245-2252.

    PubMed  PubMed Central  Google Scholar 

  51. 51.

    Delaneau O, Howie B, Cox AJ, Zagury J-F, Marchini J: Haplotype estimation using sequencing reads. Am J Hum Genet. 2013, 93: 687-696.

    PubMed  PubMed Central  Google Scholar 

  52. 52.

    Efros A, Halperin E: Haplotype reconstruction using perfect phylogeny and sequence data. BMC Bioinformatics. 2012, 13 (Suppl 6): S3-

    PubMed  PubMed Central  Google Scholar 

  53. 53.

    Lajugie J, Mukhopadhyay R, Schizas M, Lailler N, Fourel N, Bouhassira EE: Complete genome phasing of family quartet by combination of genetic, physical and population-based phasing analysis. PLoS One. 2013, 8: e64571-

    PubMed  PubMed Central  Google Scholar 

  54. 54.

    El-Metwally S, Hamza T, Zakaria M, Helmy M: Next-generation sequence assembly: four stages of data processing and computational challenges. PLoS Comput Biol. 2013, 9: e1003345-

    PubMed  PubMed Central  Google Scholar 

  55. 55.

    Geraci F: A comparison of several algorithms for the single individual SNP haplotyping reconstruction problem. Bioinformatics. 2010, 26: 2217-2225.

    PubMed  PubMed Central  Google Scholar 

  56. 56.

    Deng F, Cui W, Wang L: A highly accurate heuristic algorithm for the haplotype assembly problem. BMC Genomics. 2013, 14 (Suppl 2): S2-

    PubMed  PubMed Central  Google Scholar 

  57. 57.

    Lo C, Bashir A, Bansal V, Bafna V: Strobe sequence design for haplotype assembly. BMC Bioinformatics. 2011, 12 (Suppl 1): S24-

    PubMed  PubMed Central  Google Scholar 

  58. 58.

    Taussig DM, McGinn S, Gut IG: DNA sequencing - spanning the generations. N Biotechnol. 2013, 30: 366-372.

    Google Scholar 

  59. 59.

    Korbel JO, Lee C: Genome assembly and haplotyping with Hi-C. Nat Biotechnol. 2013, 31: 1099-1101.

    PubMed  Google Scholar 

  60. 60.

    Burton JN, Adey A, Patwardhan RP, Qiu R, Kitzman JO, Shendure J: Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013, 31: 1119-1125.

    PubMed  PubMed Central  Google Scholar 

  61. 61.

    Kaplan N, Dekker J: High-throughput genome scaffolding from in vivo DNA interaction frequency. Nat Biotechnol. 2013, 31: 1143-1147.

    PubMed  Google Scholar 

  62. 62.

    Selvaraj SR, Dixon J, Bansal V, Ren B: Whole-genome haplotype reconstruction using proximity-ligation and shotgun sequencing. Nat Biotechnol. 2013, 31: 1111-1118.

    PubMed  PubMed Central  Google Scholar 

  63. 63.

    Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, Kang HM, Marth GT, McVean GA: An integrated map of genetic variation from 1, 092 human genomes. Nature. 2012, 491: 56-65.

    PubMed  Google Scholar 

  64. 64.

    Roach JC, Glusman G, Smit AFA, Huff CD, Hubley R, Shannon PT, Rowen L, Pant KP, Goodman N, Bamshad M, Shendure J, Drmanac R, Jorde LB, Hood L, Galas DJ: Analysis of genetic inheritance in a family quartet by whole-genome sequencing. Science. 2010, 328: 636-639.

    PubMed  PubMed Central  Google Scholar 

  65. 65.

    Wijsman EM, Rothstein JH, Thompson EA: Multipoint linkage analysis with many multiallelic or dense diallelic markers: Markov chain-Monte Carlo provides practical approaches for genome scans on general pedigrees. Am J Hum Genet. 2006, 79: 846-858.

    PubMed  PubMed Central  Google Scholar 

  66. 66.

    Dewey FE, Chen R, Cordero SP, Ormond KE, Caleshu C, Karczewski KJ, Whirl-Carrillo M, Wheeler MT, Dudley JT, Byrnes JK, Cornejo OE, Knowles JW, Woon M, Sangkuhl K, Gong L, Thorn CF, Hebert JM, Capriotti E, David SP, Pavlovic A, West A, Thakuria JV, Ball MP, Zaranek AW, Rehm HL, Church GM, West JS, Bustamante CD, Snyder M, Altman RB: et al, Phased whole-genome genetic risk in a family quartet using a major allele reference sequence. PLoS Genet. 2011, 7: e1002280-

    PubMed  PubMed Central  Google Scholar 

  67. 67.

    Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, ODonnell CJ, de Bakker PIW: SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008, 24: 2938-2939.

    PubMed  PubMed Central  Google Scholar 

  68. 68.

    Browning SR, Browning BL: Haplotype phasing: existing methods and new developments. Nat Rev Genet. 2011, 12: 703-714.

    PubMed  PubMed Central  Google Scholar 

  69. 69.

    Abecasis G, Altschuler D, Auton A, Brooks L, Durbin R, Gibbs R, Hurles M, McVean G: A map of human genome variation from population-scale sequencing. Nature. 2010, 467: 1061-1073.

    PubMed  Google Scholar 

  70. 70.

    Clark AG: Inference of haplotypes from PCR-amplified samples of diploid populations. Mol Biol Evol. 1990, 7: 111-122.

    PubMed  Google Scholar 

  71. 71.

    Zagury J-F, Marchini J, Delaneau O: A linear complexity phasing method for thousands of genomes. Nat Methods. 2011, 9: 179-181.

    PubMed  Google Scholar 

  72. 72.

    OConnell J, Gurdasani D, Delaneau O, Pirastu N, Ulivi S, Cocca M, Traglia M, Huang J, Huffman JE, Rudan I, McQuillan R, Fraser RM, Campbell H, Polasek O, Asiki G, Ekoru K, Hayward C, Wright AF, Vitart V, Navarro P, Zagury J-F, Wilson JF, Toniolo D, Gasparini P, Soranzo N, Sandhu MS, Marchini J: A general approach for haplotype phasing across the full spectrum of relatedness. PLoS Genet. 2014, 10: e1004234-

    Google Scholar 

  73. 73.

    Browning BL, Browning SR: A fast, powerful method for detecting identity by descent. Am J Hum Genet. 2011, 88: 173-182.

    PubMed  PubMed Central  Google Scholar 

  74. 74.

    Brinza D, Zelikovsky A: 2SNP: scalable phasing method for trios and unrelated individuals. IEEE/ACM Trans Comput Biol Bioinform. 2008, 5: 313-318.

    PubMed  Google Scholar 

  75. 75.

    Su S-Y, Kasberger J, Baranzini S, Byerley W, Liao W, Oksenberg J, Sherr E, Jorgenson E: Detection of identity by descent using next-generation whole genome sequencing data. BMC Bioinformatics. 2012, 13: 121-

    PubMed  PubMed Central  Google Scholar 

  76. 76.

    Rao W, Ma Y, Ma L, Zhao J, Li Q, Gu W, Zhang K, Bond VC, Song Q: High-resolution whole-genome haplotyping using limited seed data. Nat Methods. 2013, 10: 6-7.

    PubMed  Google Scholar 

  77. 77.

    Kong A, Masson G, Frigge ML, Gylfason A, Zusmanovich P, Thorleifsson G, Olason PI, Ingason A, Steinberg S, Rafnar T, Sulem P, Mouy M, Jonsson F, Thorsteinsdottir U, Gudbjartsson DF, Stefansson H, Stefansson K: Detection of sharing by descent, long-range phasing and haplotype imputation. Nat Genet. 2008, 40: 1068-1075.

    PubMed  PubMed Central  Google Scholar 

  78. 78.

    He D, Eskin E: Hap-seqX: expedite algorithm for haplotype phasing with imputation using sequence data. Gene. 2013, 518: 2-6.

    PubMed  Google Scholar 

  79. 79.

    Trowsdale J, Knight JC: Major histocompatibility complex genomics and human disease. Annu Rev Genomics Hum Genet. 2013, 14: 301-323.

    PubMed  PubMed Central  Google Scholar 

  80. 80.

    Kitzman JO, Snyder MW, Ventura M, Lewis AP, Qiu R, Simmons LE, Gammill HS, Rubens CE, Santillan DA, Murray JC, Tabor HK, Bamshad MJ, Eichler EE, Shendure J: Noninvasive whole-genome sequencing of a human fetus. Sci Transl Med. 2012, 4: 137ra76-

    PubMed  PubMed Central  Google Scholar 

  81. 81.

    Liu H, Motoda H: On issues of instance selection. Data Min Knowl Discov. 2002, 6: 115-130.

    Google Scholar 

  82. 82.

    Lo C, Liu R, Lee J, Robasky K, Byrne S, Lucchesi C, Aach J, Church G, Bafna V, Zhang K: On the design of clone-based haplotyping. Genome Biol. 2013, 14: R100-

    PubMed  PubMed Central  Google Scholar 

  83. 83.

    Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, Maglott DR: ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014, 42 (Database issue): D980-D985.

    PubMed  PubMed Central  Google Scholar 

  84. 84.

    Fong CYI, Mumford AD, Likeman MJ, Jardine PE: Cerebral palsy in siblings caused by compound heterozygous mutations in the gene encoding protein C. Dev Med Child Neurol. 2010, 52: 489-493.

    PubMed  Google Scholar 

  85. 85.

    Lupski JR, Reid JG, Gonzaga-Jauregui C, Rio Deiros D, Chen DCY, Nazareth L, Bainbridge M, Dinh H, Jing C, Wheeler DA, McGuire AL, Zhang F, Stankiewicz P, Halperin JJ, Yang C, Gehman C, Guo D, Irikat RK, Tom W, Fantin NJ, Muzny DM, Gibbs RA: Whole-genome sequencing in a patient with Charcot-Marie-Tooth neuropathy. N Engl J Med. 2010, 362: 1181-1191.

    PubMed  PubMed Central  Google Scholar 

  86. 86.

    McLaughlin HM, Sakaguchi R, Liu C, Igarashi T, Pehlivan D, Chu K, Iyer R, Cruz P, Cherukuri PF, Hansen NF, Mullikin JC, Biesecker LG, Wilson TE, Ionasescu V, Nicholson G, Searby C, Talbot K, Vance JM, Zchner S, Szigeti K, Lupski JR, Hou Y-M, Green ED, Antonellis A: Compound heterozygosity for loss-of-function lysyl-tRNA synthetase mutations in a patient with peripheral neuropathy. Am J Hum Genet. 2010, 87: 560-566.

    PubMed  PubMed Central  Google Scholar 

  87. 87.

    Emison ES, McCallion AS, Kashuk CS, Bush RT, Grice E, Lin S, Portnoy ME, Cutler DJ, Green ED, Chakravarti A: A common sex-dependent mutation in a RET enhancer underlies Hirschsprung disease risk. Nature. 2005, 434: 857-863.

    PubMed  Google Scholar 

  88. 88.

    Grice EA, Rochelle ES, Green ED, Chakravarti A, McCallion AS: Evaluation of the RET regulatory landscape reveals the biological relevance of a HSCR-implicated enhancer. Hum Mol Genet. 2005, 14: 3837-3845.

    PubMed  Google Scholar 

  89. 89.

    De Vooght KMK, van Wijk R, van Solinge WW: Management of gene promoter mutations in molecular diagnostics. Clin Chem. 2009, 55: 698-708.

    PubMed  Google Scholar 

  90. 90.

    Ludlow LB, Schick BP, Budarf ML, Driscoll DA, Zackai EH, Cohen A, Konkle BA: Identification of a mutation in a GATA binding site of the platelet glycoprotein Ibbeta promoter resulting in the Bernard-Soulier syndrome. J Biol Chem. 1996, 271: 22076-22080.

    PubMed  Google Scholar 

  91. 91.

    Huang Y-C, Lee C-M, Chen M, Chung M-Y, Chang Y-H, Huang WJ-S, Ho DM-T, Pan C-C, Wu TT, Yang S, Lin M-W, Hsieh J-T, Chen Y-MA: Haplotypes, loss of heterozygosity, and expression levels of glycine N-methyltransferase in prostate cancer. Clin Cancer Res. 2007, 13: 1412-1420.

    PubMed  Google Scholar 

  92. 92.

    De Bakker PIW, McVean G, Sabeti PC, Miretti MM, Green T, Marchini J, Ke X, Monsuur AJ, Whittaker P, Delgado M, Morrison J, Richardson A, Walsh EC, Gao X, Galver L, Hart J, Hafler DA, Pericak-Vance M, Todd JA, Daly MJ, Trowsdale J, Wijmenga C, Vyse TJ, Beck S, Murray SS, Carrington M, Gregory S, Deloukas P, Rioux JD: A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC. Nat Genet. 2006, 38: 1166-1172.

    PubMed  PubMed Central  Google Scholar 

  93. 93.

    Petersdorf EW, Malkki M, Gooley TA, Martin PJ, Guo Z: MHC haplotype matching for unrelated hematopoietic cell transplantation. PLoS Med. 2007, 4: e8-

    PubMed  PubMed Central  Google Scholar 

  94. 94.

    Glusman G: Clinical applications of sequencing take center stage. Genome Biol. 2013, 14: 303-

    PubMed  PubMed Central  Google Scholar 

  95. 95.

    Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne EH, McCarroll SA, Gaudet R, Schaffner SF, Lander ES, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H:et al, Genome-wide detection and characterization of positive selection in human populations. Nature. 2007, 449: 913-918.

    PubMed  PubMed Central  Google Scholar 

  96. 96.

    Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz MH-Y, Hansen NF, Durand EY, Malaspinas A-S, Jensen JD, Marques-Bonet T, Alkan C, Prfer K, Meyer M, Burbano HA, Good JM, Schultz R, Aximu-Petri A, Butthof A, Hber B, Hffner B, Siegemund M, Weihmann A, Nusbaum C, Lander ES, Russ C:et al, A draft sequence of the Neandertal genome. Science. 2010, 328: 710-722.

    PubMed  Google Scholar 

  97. 97.

    Lawson DJ, Hellenthal G, Myers S, Falush D: Inference of population structure using dense haplotype data. PLoS Genet. 2012, 8: e1002453-

    PubMed  PubMed Central  Google Scholar 

  98. 98.

    Patterson N, Hattangadi N, Lane B, Lohmueller KE, Hafler DA, Oksenberg JR, Hauser SL, Smith MW, OBrien SJ, Altshuler D, Daly MJ, Reich D: Methods for high-density admixture mapping of disease genes. Am J Hum Genet. 2004, 74: 979-1000.

    PubMed  PubMed Central  Google Scholar 

  99. 99.

    Marchini J, Howie B, Myers S, McVean G, Donnelly P: A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007, 39: 906-913.

    PubMed  Google Scholar 

  100. 100.

    Kang HM, Zaitlen NA, Eskin E: EMINIM: an adaptive and memory-efficient algorithm for genotype imputation. J Comput Biol. 2010, 17: 547-560.

    PubMed  PubMed Central  Google Scholar 

  101. 101.

    Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229.

    PubMed  Google Scholar 

  102. 102.

    Langefeld CD, Fingerlin TE: Association methods in human genetics. Methods Mol Biol. 2007, 404: 431-460.

    PubMed  Google Scholar 

  103. 103.

    Liu N, Zhang K, Zhao H: Haplotype-association analysis. Adv Genet. 2008, 60: 335-405.

    PubMed  Google Scholar 

  104. 104.

    Murdoch JD, Speed WC, Pakstis AJ, Heffelfinger CE, Kidd KK: Worldwide population variation and haplotype analysis at the serotonin transporter gene SLC6A4 and implications for association studies. Biol Psychiatry. 2013, 74: 879-889.

    PubMed  Google Scholar 

  105. 105.

    Kasowski M, Grubert F, Heffelfinger C, Hariharan M, Asabere A, Waszak SM, Habegger L, Rozowsky J, Shi M, Urban AE, Hong M-Y, Karczewski KJ, Huber W, Weissman SM, Gerstein MB, Korbel JO, Snyder M: Variation in transcription factor binding among humans. Science. 2010, 328: 232-235.

    PubMed  PubMed Central  Google Scholar 

  106. 106.

    Tycko B: Allele-specific DNA methylation: beyond imprinting. Hum Mol Genet. 2010, 19: R210-R220.

    PubMed  PubMed Central  Google Scholar 

  107. 107.

    Aretz S, Uhlhaas S, Caspari R, Mangold E, Pagenstecher C, Propping P, Friedl W: Frequency and parental origin of de novo APC mutations in familial adenomatous polyposis. Eur J Hum Genet. 2004, 12: 52-58.

    PubMed  Google Scholar 

  108. 108.

    Olson MV, Green P: Criterion for the completeness of large-scale physical maps of DNA. Cold Spring Harb Symp Quant Biol. 1993, 58: 349-355.

    PubMed  Google Scholar 

  109. 109.

    Ewing B, Green P: Base-calling of automated sequencer traces using phred, II. Error probabilities. Genome Res. 1998, 8: 186-194.

    PubMed  Google Scholar 

  110. 110.

    Lin S, Cutler DJ, Zwick ME, Chakravarti A: Haplotype inference in random population samples. Am J Hum Genet. 2002, 71: 1129-1137.

    PubMed  PubMed Central  Google Scholar 

  111. 111.

    Rastas P, Koivisto M, Mannila H, Ukkonen E: Phasing genotypes using a hidden Markov model. Bioinformics Algorithms :Techniques and Applications. Edited by: Mǎndoiu II, Zelikovsky A. 2008, 355-362. John Wiley & Sons, Inc, Hoboken, NJ

    Google Scholar 

Download references

Acknowledgments

GG, HCC and JCR received support from the University of Luxembourg - Institute for Systems Biology Program. GG and JCR received support from and the National Institute of General Medical Sciences Center for Systems Biology (P50 GM076547). We thank the anonymous reviewers for their contributions.

Author information

Affiliations

Authors

Corresponding author

Correspondence to Jared C Roach.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Glusman, G., Cox, H.C. & Roach, J.C. Whole-genome haplotyping approaches and genomic medicine. Genome Med 6, 73 (2014). https://doi.org/10.1186/s13073-014-0073-7

Download citation

Keywords

  • Major Histocompatibility Complex
  • Human Genome Project
  • Genomic Medicine
  • Compound Heterozygosity
  • Phase Accuracy