Direct determination of molecular haplotypes by chromosome microdissection.

Direct observation of haplotypes is still technically challenging. Here we report a method for the determination of haplotypes through chromosome microdissection. We determined human long-range chromosomal haplotypes with more than 98.85% accuracy at 24,245 genome-wide heterozygous single-nucleotide polymorphism (SNP) loci.

Haplotype refers to a group of alleles inherited on a single chromosome. Haplotype analysis plays an important role for mapping disease genes, elucidating population histories, studying the evolutionary genetics, and exploring cis-interactions in the regulation of gene expression 1,2. Statistical and computational methods have been developed to reconstruct haplotypes from conventional genotype data. Although statistical and computational methods can provide highly accurate estimates of haplotypes in most situations, even the best method can have a large error rate of more than 35% for some real datasets 3. The accuracy of statistical methods for haplotype reconstruction can be greatly improved when they are applied to trio datasets using information from relatives 2; however, the parental samples are not always available to make this feasible. Experimental approaches have been developed to determine haplotypes 4,5 with limited capability of haplotyping SNPs across a large chromosomal region spanning megabase distances, except the conversion technology, which can determine the long-range haplotypes by creating human-mouse hybrid cells that contain a subset of human chromosomes 6,7, and the polony (polymerase colony) technology, which can determine the long-range haplotype of a pair of distant SNP loci on the same immobilized colony within an acrylamide gel 8,9. Conversion technology involves considerable cost and time for large-scale applications, the polony technology is limited by the small number of SNPs that can be haplotyped in genome-wide studies 2,10.

Here we report , 7 (poly) Dimensional DNA (7DDNA), a method built upon the concept that a somatic cell of a diploid organism contains two and only two sets of homologous chromosomes. If we collect ~ 23 chromosomes from one cell instead of its entire 46 chromosomes, the collection may contain only one copy of some chromosomes, it may also contain no copy or both copies of other chromosomes (Supplementary Fig. 1).We genotype the 7DDNA sample and an entire genomic sample in parallel; the latter allows us to determine which SNPs are heterozygous or homozygous in each sample. . If a 7DDNA sample shows “homozygous” genotypes at all polymorphic loci along a chromosome known to be heterozygous, this 7DDNA sample contains one single copy of this chromosome (Fig. 1, Supplementary Fig. 1 and Supplementary Table 1), which will enable us to read one haplotype of this chromosome directly from its conventional high-throughput genotyping output. The unobserved dropout alleles will constitute the other haplotype of this chromosome. If a 7DDNA sample shows heterozygous genotypes at loci known to be polymorphic, this harvest contains both chromosome copies. If a 7DDNA sample fails to call any SNPs along a chromosome, this chromosome is absent in this sample.

Figure 1

The principle and results on chromosome 5 of the 7DDNA haplotyping method. (a) 7DDNA haplotyping procedure. (b) haplotyped chromosome 5, on which 2,514 SNPs were resolved. Each dot represents a SNP. Scale bars, 10 μM.

We have performed a proof-of-principle study with a HapMap sample 11, GM10847, since the haplotypes can be computationally reconstructed with the parental information2,12. These computationally configured haplotypes should be accurate on those non-triple-heterozygote SNP loci unless there is mis-inheritance, genotyping errors, SNP recurrence, or an overlapping of a SNP with an insertion-deletion polymorphism (indel), a repeat, or a copy number variation (CNV), or segmental duplications.

In our experiment, we microdissected three pieces from the chromosome spread of an individual (Online Methods). The microdissection harvests were subjected to whole genome amplification (WGA), which provided us with about 5-8 μg DNA from each harvest. After genotyping with an Illumina BeadChip (CNV370), we collectively haplotyped 15 chromosomes in these three harvested samples, among which eight chromosomes were covered by one harvest, five chromosomes were covered by two harvests, and one chromosome was covered by all three harvests (Supplementary Table 2). Totally 24,245 heterozygous SNPs on these 15 chromosomes were haplotyped (Supplementary Table 2), among which 2,089 SNPs were phased by more than one microdissection harvests, 98.85% of these SNPs (2,065 SNPs) showed consistent results between different microdissection harvests (Table 1). Our data showed a faithful amplification of our microdissection harvests, in which the error rate was 0.73% (Supplementary Table 3). We noticed that the no-call rate in the Illumina genotyping on haploid chromosomes of the microdissection harvests was as high as ~70% (Supplementary Table 4). The possible reasons may be that our WGA method was insufficient to amplify every locus when the quantity of input DNA was as small as one single copy, or this genotyping platform was insensitive to the DNA template amplified from one single molecule of template. We anticipate that this issue may be solved by a new robust WGA method for amplifying single-copy template or other genotyping platforms such as next-generation sequencers.

Table 1

Replication of the 7DDNA haplotyping results.

Chromosome	Repeated SNPs*	Consistent SNPs**	Consistency rate %
Chr3	517	516	99.81
Chr4	565	554	98.05
Chr16	189	186	98.41
Chr17	217	215	99.08
Chr18	212	208	98.11
Chr20	389	386	99.23
Total	2,089	2,065	98.85

Repeated SNPs are those SNPs that were haplotyped in 2 or 3 microdissection harvests in our study.

Consistent SNPs are those repeated SNPs that obtained the identical results on haplotypes from different microdissection harvests.

Compared with computationally reconstructed haplotypes of this person (GM10847) on those non-triple-heterozygous SNP loci, our 7DDNA haplotypes reported discordant allele phases on only 0.65% of SNPs (Supplementary Table 2). Comparing our experimentally determined haplotypes of GM10847 with her parents’ haplotypes downloaded from the HapMap database that were computationally configured by IMPUTE++ 13, the difference was only 2.15%, among which 58.7% of those discrepancies occurred on triple-heterozygote loci. This data not only demonstrated a high accuracy of haplotypes configured with trio information 2, it also validated our 7DDNA haplotype results.

About 10.8% of heterozygous SNPs on the diploid chromosomes in our microdissection harvests remained to be heterozygous (Supplementary Table 5); the other heterozygous SNPs on those diploid chromosomes received a homozygous genotype call due to allele dropout. However, on the haploid chromosomes in our microdissection harvests, only 0.08% of heterozygous SNPs appeared to be heterozygous (Supplementary Table 5). This substantial disparity (10.8% vs. 0.08%) enabled us to distinguish those haploid chromosomes from diploid chromosomes in a microdissection harvest (Supplementary Fig. 2 and Supplementary Fig. 3).

To haplotype all chromosomes of each individual, multiple microdissection harvests should be collected to cover all of a person's chromosomes. Our statistical analysis (Supplementary Note-1) showed that 4-8 microdissection samplings of each individual will provide a 93% probability to obtain a person's whole genome haplotype (Supplementary Table 6 and Supplementary Fig. 4). In other words, 52% of individuals may have completed their genome haplotyping with no more than five microdissections, 93% of individuals will need no more than eight microdissections to complete their genome haplotyping (Supplementary Table 6 and Supplementary Fig. 4). When we reach 12 microdissection samplings, 99.6% of all individuals will have their whole genome haplotyped. In addition to chromosome coverage, this cumulative strategy will also effectively increase the SNP coverage on resolved haplotypes (Supplementary Table 7). We have estimated the cost of this strategy to be $0.0023 per SNP per sample (Supplementary Note-2). The entire 7DDNA procedure will take three days, with a ~10-min work on the first day, no work on the second day, and a 6-8 hour work on the third day. Based on this schedule, we may do five experiments each week, and handle 24 individuals in each experiment. A 1,000-sample project may need only 2-6 months.

Because many researchers may be interested in particular chromosomes rather than the entire genome, such as in human leukocyte antigen (HLA) studies, we have explored an alternative approach in which a single chromosome is collected in a microdissection with an Objective HCX 150X/0.90 lens. Chromosome recognition can be achieved by karyotype morphology in combination with PCR validation.

It is possible that a chromosome is broken by the microdissection procedure. When only a partial chromosome is harvested without accompaniment of its sister chromosome, there will be failed allele calls at those SNP loci located on the absent chromosomal segment. If a partial chromosome is harvested accompanied by its full-length sister chromosome in the microdissection, the heterozygous SNPs located on this partial chromosomal segment will receive heterozygous genotypes in the 7DDNA genotyping results, but those heterozygous SNPs located outside this partial chromosome piece will appear as “homozygous” loci in the 7DDNA genotyping results (Supplementary Fig. 2). The segmental boundaries will be deciphered from the genotyping results and the haplotypes of entire chromosome will be assembled from broken pieces in different microdissection harvests by their overlapping portion.

The 7DDNA technology provides the first high-throughput experimental method for determination of long-range haplotypes for large-scale and genome-wide studies. Compared with conversion technology 6,7,14 and polony haplotyping 8,9, 7DDNA exhibits advantages on labor, time, throughput and the number of SNPs. Chromosomal haplotypes will be important for functional interpretations of SNPs in cis-interactions and the integration of genetic and epigenetic datasets along the chromosomes.

AOP:

Comparing the haplotypes of a few randomly microdissected chromosomes to a full genome wide haplotype allows one to determine the long-range haplotype of chromosomes for which only one copy was captured.

Issue

Comparing the haplotypes of a few randomly microdissected chromosomes to a full genome wide haplotype allows one to determine the long-range haplotype of chromosomes for which only one copy was captured.