亲缘群体单体型推断的研究进展

随着高通量的基因分型技术的出现,越来越多的单核苷酸多态性( single nucleotide polymorphism,SNP)可被检测。基因组研究中,单体型分析已被运用到SNP位点与复杂性状的关联分析中。现有的单体型推断方法包括基于非亲缘群体的方法,DNA pool样本方法和系谱资料方法。不过家庭及亲缘数据可为单体型推断提供很重要的信息,利用这些信息进行单体型推断会增加准确性。如今对亲缘群体单体型推断的算法和软件已有很多报道。该文的目的是回顾这些算法、软件以及它们的优缺点和适用的群体类型,同时探讨仍然存在的问题和面临的挑战。     

Noninvasive prenatal diagnosis by genome-wide haplotyping of cell-free plasma DNA

PurposeWhereas noninvasive prenatal screening for aneuploidies is widely implemented, there is an increasing need for universal approaches for noninvasive prenatal screening for monogenic diseases. Here, we present a cost-effective, generic cell-free fetal DNA (cffDNA) haplotyping approach to scan the fetal genome for the presence of inherited monogenic diseases.MethodsFamilies participating in the preimplantation genetic testing for monogenic disorders (PGT-M) program were recruited for this study. Two hundred fifty thousand single-nucleotide polymorphisms (SNPs) captured from maternal plasma DNA along with genomic DNA from family members were massively parallel sequenced. Parental genotypes were phased via an available genotype from a close relative, and the fetal genome-wide haplotype and copy number were determined using cffDNA haplotyping analysis based on estimation and segmentation of fetal allele presence in the maternal plasma.ResultsIn all families tested, mutational profiles from cffDNA haplotyping are consistent with embryo biopsy profiles. Genome-wide fetal haplotypes are on average 97% concordant with the newborn haplotypes and embryo haplotypes.ConclusionWe demonstrate that genome-wide targeted capture and sequencing of polymorphic SNPs from maternal plasma cell-free DNA (cfDNA) allows haplotyping and copy-number profiling of the fetal genome during pregnancy. The method enables the accurate reconstruction of the fetal haplotypes and can be easily implemented in clinical practice.     

Single-molecule analysis for molecular haplotyping

In the genome era, there is great hope that genetic approaches such as linkage equilibrium mapping can be used to study common human disorders using a case-control population association study design. Ideally, the parental chromosomes are marked so that chromosomal regions in the form of haplotypes are compared in these studies to increase the power of association. Determining the haplotypes in a diploid individual is a major technical challenge in genetic studies of complex traits. A molecular approach to haplotyping is therefore highly desirable. Recent advances in DNA preparation, separation, labeling, and image analysis provide hope that a strategy of using a three-dye system coupled with DNA distance measurements between alleles will yield haplotype information of sufficiently high quality for genetic studies. In this work, we present the outline of the major challenges one must meet in developing a robust strategy for SNP detection and molecular haplotyping using single molecule analysis.     

A gel-free SNP genotyping method: bioluminometric assay coupled with modified primer extension reactions (BAMPER) directly from double-stranded PCR products

Inexpensive, high-throughput genotyping methods are needed for analyzing human genetic variations. We have successfully applied the regular bioluminometric assay coupled with modified primer extension reactions (BAMPER) method to single-nucleotide polymorphism (SNP) typing as well as the allele frequency determination for various SNPs. This method includes the production of single-strand target DNA from a genome and a primer extension reaction coupled with inorganic pyrophosphate (PPi) detection by a bioluminometric assay. It is an efficient way to get accurate allele frequencies for various SNPs, while single-strand DNA preparation is labor intensive. The procedure can be simplified in the typing of SNPs. We demonstrate that a modified BAMPER method in which we need not prepare a single-strand DNA can be carried out in one tube. A PCR product is directly used as a template for SNP typing in the new BAMPER method. Generally, tremendous amounts of PPi are produced in a PCR process, as well as many residual dNTPs, and residual PCR primers remain in the PCR products, which cause a large background signal in a bioluminometric assay. Here, shrimp alkaline phosphatase (SAP) and E. coli exonuclease I were used to degrade these components prior to BAMPER detection. The specific primer extension reactions in BAMPER were carried out under thermocycle conditions. The primers were extended to produce large amounts of PPi only when their bases at 3'-termini were complementary to the target. The extension products, PPis, were converted to ATP to be analyzed using the luciferin-luciferase detection system. We successfully demonstrated that PCR products can be directly genotyped by BAMPER in one tube for SNPs with various GC contents. As all reactions can be carried out in a single tube, the method will be useful for realizing a fully automated genotyping system.     

Linear allele-specific long-range amplification: a novel method of long-range molecular haplotyping

Haplotypes have been repeatedly shown to be more powerful than collections of single-locus markers in gene-mapping studies. Various haplotyping methods including statistical estimation are employed but molecular haplotyping, the acquisition of information directly on physical DNA sequences, has been in demand for its accuracy and independence from family pedigrees. We investigated the allelic specificity of long-range PCR, which was successful for long-range haplotyping in recent reports, and found problems of initial mispriming and crossover amplification significantly confounded its application. Based on these observations, we designed a novel method based on linear amplification of a hemizygous DNA segment with a single phosphorothioate-modified oligonucleotide. Our results revealed, with a single nucleotide polymorphism as the discriminative marker, downstream haplotypes of 14-15 kb DNA segment could be confidently scored. With two rounds of the method and five single nucleotide polymorphisms, molecular haplotypes of 29.3 kb spanning the HCR and CDSN genes, two genes associated with the susceptibility of psoriasis, of 11 members, belonging to a CEPH family, were revealed. Clear Mendelian segregation of 35 highly heterozygous SNPs confirmed the accuracy of the method. Problems of low specificity associated with long-range PCR were not observed. The simplicity, along with long-sequence accessibility and feasibility of a single nucleotide difference as the discriminative marker indicated our method holds promise for future gene-mapping studies.     

Haplotype block partitioning and tag SNP selection using genotype data and their applications to association studies

Recent studies have revealed that linkage disequilibrium (LD) patterns vary across the human genome with some regions of high LD interspersed by regions of low LD. A small fraction of SNPs (tag SNPs) is sufficient to capture most of the haplotype structure of the human genome. In this paper, we develop a method to partition haplotypes into blocks and to identify tag SNPs based on genotype data by combining a dynamic programming algorithm for haplotype block partitioning and tag SNP selection based on haplotype data with a variation of the expectation maximization (EM) algorithm for haplotype inference. We assess the effects of using either haplotype or genotype data in haplotype block identification and tag SNP selection as a function of several factors, including sample size, density or number of SNPs studied, allele frequencies, fraction of missing data, and genotyping error rate, using extensive simulations. We find that a modest number of haplotype or genotype samples will result in consistent block partitions and tag SNP selection. The power of association studies based on tag SNPs using genotype data is similar to that using haplotype data.     

正常汉族人群FcαR1基因结构研究

目的：构建我国汉族人群FcαR1 基因单倍型及单倍型域。方法：提取100例正常广东汉族人基因组，分段扩增FcαR1基因并直接测序鉴定SNP基因型。选取MAF>0.10的SNP构建单倍型及单倍型域。结果：汉族人群FcαR1基因共有3个单倍型域，每个单倍型域中包含的单倍型数3至7个，平均为5.33个，单倍型域中常见单倍型数2至3个，平均2.3个。结论：FcαR1 基因存在单倍型及单倍型域，本研究结果将为进一步以该基因为候选基因进行疾病基因定位奠定基础。     

Exo-proofreading,a versatile SNP scoring technology

We report the validation of a new assay for typing single nucleotide polymorphisms (SNPs) that takes advantage of the 3'-to-5' exonuclease proofreading activity of many DNA polymerases. The assay uses one or more primers labeled on the 3' nucleotide base, and can be implemented in a variety of formats including a one-step PCR reaction that allows SNP typing directly from genomic DNA samples. The detection of genotypes can be accomplished by means of fluorescence detection on assays that have been purified to remove excess primer, or by means of fluorescence polarization without any additional cleanup. We also demonstrate that the Exo-Proofreading SNP assay can be used on pooled samples to obtain allele frequency data.     

Pyrosequencing-based SNP allele frequency estimation in DNA pools

Association screening involving numerous genetic markers is facilitated by the analysis of pooled DNA samples rather than individual samples. Several genotyping methods have shown high accuracy and precision of allele frequency estimation in pools. Here, we expand the validation of SNP allele frequency estimation in DNA pools using Pyrosequencing by analyzing 186 pools for three SNPs representing complex sequencing cases. The correlation coefficient between estimated and true allele frequencies ranged between 0.979 and 0.996 and tended to increase with pool size, whereas the difference between estimated and true allele frequencies was 2.37+/-0.11%, in post-PCR pools. The precision was 1.73%. Pool size had no significant effect on accuracy and precision. A comparison between post-PCR and pre-PCR pools showed that for pre-PCR pooling efforts to accurately quantify the genomic DNA samples to be pooled and subsequently amplified are critical. To conclude, Pyrosequencing can be used for allele frequency estimation in DNA pools of SNPs with complex sequencing scenarios with accuracy and precision values in ranges comparable with those of other SNP typing techniques. Considering the ease of use, short run and analysis times, and little instrument maintenance requirements, Pyrosequencing may even be a preferred option.     

Mitochondrial DNA haplotyping revealed the presence of mixed up benign and neoplastic tissue sections from two individuals on the same prostatic biopsy slide

DNA typing was requested to investigate a presumptive cancer diagnosis error by confirming whether benign and cancerous prostatic tissue in the same presurgical haematoxylin and eosin stained slide belonged to the same person. After independent histological re-examination of the slide by a pathologist, manual slide dissection was used to guarantee independent and high recovery DNA isolation from each tissue section, avoiding carryover and background contamination. Nuclear DNA quantification performed by real time polymerase chain reaction (PCR) revealed the absence of human DNA for short tandem repeat (STR) typing. Mitochondrial DNA was only obtained by performing PCR of very short fragments ( approximately 100 bp), indicating high DNA degradation. Different low frequency hypervariable region I haplotypes were obtained from each tissue section (normal tissue section haplotype: 16224C, 16234T, 16311C, 16356C; cancer tissue section haplotype: 16256T, 16270T, 16293G). Only the normal tissue section haplotype matched that obtained from the patient's blood sample, indicating that the cancer tissue section originated from an unknown patient. These results supported the hypothesis of sample mix up during block processing or slide preparation by a carryover mechanism. Mitochondrial genetic typing is recommended to exclude the possibility of carryover artefacts when low DNA content and high degradation compromise conventional STR typing.     

Determination of haplotypes from single DNA molecules: a method for single-molecule barcoding

Determining the haplotypes in a diploid individual is a major technical challenge in genetic studies of human complex traits. Here we report a method of molecular haplotyping by directly imaging multiple polymorphic sites on individual DNA molecules simultaneously. DNA fragments amplified by long-range PCR were labeled with fluorescent dyes at each polymorphic site using a modified gap-filled padlock probe ligation approach. The labeled DNA molecules were then stretched into linear form on a functionalized glass surface and imaged with multicolor total internal reflection fluorescence (TIRF) microscopy. By determining the colors and positions of the fluorescent labels with respect to the backbone at polymorphic sites, the haplotype can be inferred accurately, in a manner similar to reading a barcode, even when the DNA fragments are not fully labeled. The feasibility of this technology is demonstrated by the determination of the haplotype of a 9.3-kbp DNA fragment containing four SNPs.     

A novel method for SNP detection using a new duplex-specific nuclease from crab hepatopancreas

We have characterized a novel nuclease from the Kamchatka crab, designated duplex-specific nuclease (DSN). DSN displays a strong preference for cleaving double-stranded DNA and DNA in DNA-RNA hybrid duplexes, compared to single-stranded DNA. Moreover, the cleavage rate of short, perfectly matched DNA duplexes by this enzyme is essentially higher than that for nonperfectly matched duplexes of the same length. Thus, DSN differentiates between one-nucleotide variations in DNA. We developed a novel assay for single nucleotide polymorphism (SNP) detection based on this unique property, termed "duplex-specific nuclease preference" (DSNP). In this innovative assay, the DNA region containing the SNP site is amplified and the PCR product mixed with signal probes (FRET-labeled short sequence-specific oligonucleotides) and DSN. During incubation, only perfectly matched duplexes between the DNA template and signal probe are cleaved by DSN to generate sequence-specific fluorescence. The use of FRET-labeled signal probes coupled with the specificity of DSN presents a simple and efficient method for detecting SNPs. We have employed the DSNP assay for the typing of SNPs in methyltetrahydrofolate reductase, prothrombin and p53 genes on homozygous and heterozygous genomic DNA.     

High-throughput SNP allele-frequency determination in pooled DNA samples by kinetic PCR

We have developed an accurate, yet inexpensive and high-throughput, method for determining the allele frequency of biallelic polymorphisms in pools of DNA samples. The assay combines kinetic (real-time quantitative) PCR with allele-specific amplification and requires no post-PCR processing. The relative amounts of each allele in a sample are quantified. This is performed by dividing equal aliquots of the pooled DNA between two separate PCR reactions, each of which contains a primer pair specific to one or the other allelic SNP variant. For pools with equal amounts of the two alleles, the two amplifications should reach a detectable level of fluorescence at the same cycle number. For pools that contain unequal ratios of the two alleles, the difference in cycle number between the two amplification reactions can be used to calculate the relative allele amounts. We demonstrate the accuracy and reliability of the assay on samples with known predetermined SNP allele frequencies from 5% to 95%, including pools of both human and mouse DNAs using eight different SNPs altogether. The accuracy of measuring known allele frequencies is very high, with the strength of correlation between measured and known frequencies having an r(2) = 0.997. The loss of sensitivity as a result of measurement error is typically minimal, compared with that due to sampling error alone, for population samples up to 1000. We believe that by providing a means for SNP genotyping up to thousands of samples simultaneously, inexpensively, and reproducibly, this method is a powerful strategy for detecting meaningful polymorphic differences in candidate gene association studies and genome-wide linkage disequilibrium scans.     

应用复合诱导突变分离PCR法检测mtDNA的单核苷酸多态性

目的应用复合诱导突变分离PCR(multiplexed mutagenically separated PCR,MS-PCR)技术、银染分型,建立线粒体DNA(mitochondrial DNA,mtDNA)编码区单核苷酸多态(single nucleotide polymorphism,SNP)分型系统,探讨其应用价值。并调查了成都汉族群体mtDNA编码区4个SNP基因座等位基因频率和单倍型分布情况。方法根据SNP基因座(C12705T、A8701G、G8584A、C10400T)设计两条片段相差4个碱基的等位基因特异性引物和一条公共引物,4个SNP基因座复合扩增,PCR产物经聚丙烯酰胺凝胶电泳、银染显带后确定样本的基因型。结果不同SNP基因座为长度不同的单一谱带,其分型结果与直接测序一致。在成都汉族160名无关个体中,4个SNP基因座C12705T、A8701G、G8584A、C10400T等位基因频率分别为0.3813/0.6187、0.4813/0·5187、0.8250/0.1750、0.4938/0.5062;共检出6种单倍型,单倍型的基因多样性为0.7137。结论建立的MMS-PCR银染分型系统是一种简单、快速、准确、有效的SNP分型方法,对建立mtDNA编码区SNP数据库,研究群体遗传学、进化学和进行法医学个人识别和亲子鉴定有重要意义。     

Identification of the sources of error in allele frequency estimations from pooled DNA indicates an optimal experimental design

Genotyping costs still preclude analysis of a comprehensive SNP map in thousands of individual subjects in the search for disease susceptibility loci. Allele frequency estimation in DNA pools from cases and controls offers a partial solution, but variance in these estimates will result in some loss of statistical power. However, there has been no systematic attempt to quantify the several sources of error in previous studies. We report an analysis of the magnitude of variance components of each experimental stage in DNA pooling studies, and find that a design based on the formation of numerous small pools of approximately 50 individuals is superior to the formation of fewer, larger pools and the replication of any of the experimental stages. We conclude that this approach may retain an effective sample size greater than 68% of the true sample size, whilst offering a 60-fold reduction in DNA usage and a greater than 30-fold saving in cost, compared to individual genotyping. The possibility of combining pooling with informed selection of haplotype tag SNPs is also considered. In this way further savings in efficiency may be possible by using pooled allele frequency estimates to infer haplotype frequencies and hence, allele frequencies at untyped markers.     

Identification of the sources of error in allele frequency estimations from pooled DNA indicates an optimal experimental design

Genotyping costs still preclude analysis of a comprehensive SNP map in thousands of individual subjects in the search for disease susceptibility loci. Allele frequency estimation in DNA pools from cases and controls offers a partial solution, but variance in these estimates will result in some loss of statistical power. However, there has been no systematic attempt to quantify the several sources of error in previous studies. We report an analysis of the magnitude of variance components of each experimental stage in DNA pooling studies, and find that a design based on the formation of numerous small pools of approximately 50 individuals is superior to the formation of fewer, larger pools and the replication of any of the experimental stages. We conclude that this approach may retain an effective sample size greater than 68% of the true sample size, whilst offering a 60-fold reduction in DNA usage and a greater than 30-fold saving in cost, compared to individual genotyping. The possibility of combining pooling with informed selection of haplotype tag SNPs is also considered. In this way further savings in efficiency may be possible by using pooled allele frequency estimates to infer haplotype frequencies and hence, allele frequencies at untyped markers.     

An SNP resource for rice genetics and breeding based on subspecies indica and japonica genome alignments

Dense coverage of the rice genome with polymorphic DNA markers is an invaluable tool for DNA marker-assisted breeding, positional cloning, and a wide range of evolutionary studies. We have aligned drafts of two rice subspecies, indica and japonica, and analyzed levels and patterns of genetic diversity. After filtering multiple copy and low quality sequence, 408,898 candidate DNA polymorphisms (SNPs/INDELs) were discerned between the two subspecies. These filters have the consequence that our data set includes only a subset of the available SNPs (in particular excluding large numbers of SNPs that may occur between repetitive DNA alleles) but increase the likelihood that this subset is useful: Direct sequencing suggests that 79.8% +/- 7.5% of the in silico SNPs are real. The SNP sample in our database is not randomly distributed across the genome. In fact, 566 rice genomic regions had unusually high (328 contigs/48.6 Mb/13.6% of genome) or low (237 contigs/64.7 Mb/18.1% of genome) polymorphism rates. Many SNP-poor regions were substantially longer than most SNP-rich regions, covering up to 4 Mb, and possibly reflecting introgression between the respective gene pools that may have occurred hundreds of years ago. Although 46.2% +/- 8.3% of the SNPs differentiate other pairs of japonica and indica genotypes, SNP rates in rice were not predictive of evolutionary rates for corresponding genes in another grass species, sorghum. The data set is freely available at http://www.plantgenome.uga.edu/snp.     

Multiple ant colony algorithm method for selecting tag SNPs

The search for the association between complex disease and single nucleotide polymorphisms (SNPs) or haplotypes has recently received great attention. Finding a set of tag SNPs for haplotyping in a great number of samples is an important step to reduce cost for association study. Therefore, it is essential to select tag SNPs with more efficient algorithms. In this paper, we model problem of selection tag SNPs by MINIMUM TEST SET and use multiple ant colony algorithm (MACA) to search a smaller set of tag SNPs for haplotyping. The various experimental results on various datasets show that the running time of our method is less than GTagger and MLR. And MACA can find the most representative SNPs for haplotyping, so that MACA is more stable and the number of tag SNPs is also smaller than other evolutionary methods (like GTagger and NSGA-II). Our software is available upon request to the corresponding author.     

Polymorphism analysis within the HLA-A locus by universal oligonucleotide array

Human leukocyte antigen (HLA) class I genes present some of the most complex single nucleotide polymorphism (SNP) patterns in the human genome. HLA typing is therefore extremely challenging. In this article, we use the ligation detection reaction (LDR) combined with a universal array (UA) as a robust and efficient method to analyze SNPs within the HLA-A region that includes HLA-A alleles of interest for immunotherapy in tumor diseases. The LDR, combined with a UA platform, has been optimized for the detection of 27 alleles distributed within exons 2 and 3 of HLA-A. The assay involves the amplification by PCR of the HLA-A genomic region (1,900 bp), the cycled ligation reaction, followed by the capture of ligated products through hybridization onto a UA. Each slide was designed to allow the detection of up to eight samples in parallel. The PCR/LDR/UA HLA-A assay was evaluated by analyzing 62 individuals (31 homozygous and 31 heterozygous) previously typed by direct sequencing. We demonstrate that the microarray genotyping procedure described here is a robust and efficient method for unambiguous detection of HLA alleles. HLA genotyping by PCR/LDR/UA is in perfect agreement with typing obtained by direct sequencing. Our results clearly demonstrate that the combination of enzymatic processing (LDR) and a demultiplexing hybridization onto a UA is a robust tool for SNP discrimination within the highly polymorphic HLA region. We demonstrate the specificity and efficiency of such an approach, suggesting the feasibility of a PCR/LDR/UA low resolution HLA typing procedure.     

Efficient approach to unique single-nucleotide polymorphism discovery

Single-nucleotide polymorphisms (SNPs) are the most frequently found DNA sequence variations in the human genome. It has been argued that a dense set of SNP markers can be used to identify genetic factors associated with complex disease traits. Because all high-throughput genotyping methods require precise sequence knowledge of the SNPs, any SNP discovery approach must involve both the determination of DNA sequence and allele frequencies. Furthermore, high-throughput genotyping also requires a genomic DNA amplification step, making it necessary to develop sequence-tagged sites (STSs) that amplify only the DNA fragment containing the SNP and nothing else from the rest of the genome. In this report, we demonstrate the utility of a SNP-screening approach that yields the DNA sequence and allele frequency information while screening out duplications with minimal cost and effort. Our approach is based on the use of a homozygous complete hydatidiform mole (CHM) as the reference. With this homozygous reference, one can identify and estimate the allele frequencies of common SNPs with a pooled DNA-sequencing approach (rather than having to sequence numerous individuals as is commonly done). More importantly, the CHM reference is preferable to a single individual reference because it reveals readily any duplicated regions of the genome amplified by the PCR assay before the duplicated sequences are found in GenBank. This approach reduces the cost of SNP discovery by 60% and eliminates the costly development of SNP markers that cannot be amplified uniquely from the genome.