A model for the clustered distribution of SNPs in the human genome

Motivated by a non-random but clustered distribution of SNPs, we introduce a phenomenological model to account for the clustering properties of SNPs in the human genome. The phenomenological model is based on a preferential mutation to the closer proximity of existing SNPs. With the Hapmap SNP data, we empirically demonstrate that the preferential model is better for illustrating the clustered distribution of SNPs than the random model. Moreover, the model is applicable not only to autosomes but also to the X chromosome, although the X chromosome has different characteristics from autosomes. The analysis of the estimated parameters in the model can explain the pronounced population structure and the low genetic diversity of the X chromosome. In addition, correlation between the parameters reveals the population-wise difference of the mutation probability. These results support the mutational non-independence hypothesis against random mutation.     

Genetic and expression analysis of all non‐synonymous single nucleotide polymorphisms in the human deoxyribonuclease I‐like 1 and 2 genes

Members of the human DNase I family, DNase I‐like 1 and 2 (DNases 1L1 and 1L2), with physiological role(s) other than those of DNase I, possess three and one non‐synonymous SNPs in the genes, respectively. However, only limited population data are available, and the effect of these SNPs on the catalytic activity of the enzyme remains unknown. Genotyping of all the non‐synonymous SNPs was performed in three ethnic groups including six different populations using the PCR‐RFLP method newly developed. Asian and African groups including Japanese, Koreans, Ghanaians and Ovambos were typed as a single genotype at each SNP, but polymorphism at only SNP V122I in DNase 1L1 was found in Caucasian groups including Germans and Turks; thus a Caucasian‐specific allele was identified. The DNase 1L1 and 1L2 genes show relatively low genetic diversity with regard to these non‐synonymous SNPs. The level of activity derived from the V122I, Q170H and D227A substituted DNase 1L1 corresponding to SNPs was similar to that of the wild‐type, whereas replacement of the Asp residue at position 197 in the DNase 1L2 protein with Ala, corresponding to SNP D197A, reduced its activity greatly. Thus, SNP V122I in DNase 1L1 exhibiting polymorphism exerts no effect on the catalytic activity, and furthermore SNP D197A in DNase 1L2, affecting its catalytic activity, shows no polymorphism. These findings permit us to postulate that the non‐synonymous SNPs identified in the DNase 1L1 and 1L2 genes may exert no influence on the activity levels of DNases 1L1 and 1L2 in human populations.     

Multi objective SNP selection using pareto optimality

Biomarker discovery is a challenging task of bioinformatics especially when targeting high dimensional problems such as SNP (single nucleotide polymorphism) datasets. Various types of feature selection methods can be applied to accomplish this task. Typically, using features versus class labels of samples in the training dataset, these methods aim at selecting feature subsets with maximal classification accuracies. Although finding such class-discriminative features is crucial, selection of relevant SNPs for maximizing other properties that exist in the nature of population genetics such as the correlation between genetic diversity and geographical distance of ethnic groups can also be equally important. In this work, a methodology using a multi objective optimization technique called Pareto Optimal is utilized for selecting SNP subsets offering both high classification accuracy and correlation between genomic and geographical distances. In this method, discriminatory power of an SNP is determined using mutual information and its contribution to the genomic–geographical correlation is estimated using its loadings on principal components. Combining these objectives, the proposed method identifies SNP subsets that can better discriminate ethnic groups than those obtained with sole mutual information and yield higher correlation than those obtained with sole principal components on the Human Genome Diversity Project (HGDP) SNP dataset.     

A new set of DIP‐SNP markers for detection of unbalanced and degraded DNA mixtures

Unbalanced and degraded mixtures (UDM) are frequently encountered during forensic DNA analysis. For example, forensic DNA units regularly encounter DNA mixture signal where the DNA signal from the alleged offender is masked or swamped by high quantities of DNA from the victim. Our previous data presented a new kind of DNA markers that composed of a deletion/insertion polymorphism (DIP) and a SNP and we termed this new kind of microhaplotypes DIP‐SNP (combination of DIP and SNP). Since such markers could be designed short enough for degraded DNA amplification, we hypothesized that DIP‐SNP markers are applicable for typing of UDM. In this study, we developed a new set of DIP‐SNPs with short amplicons which were complement to our prior developed system. The multiplex PCR and SNaPshot assay were established for 20 DIP‐SNPs in a Chinese Han population. The DIP‐SNPs were capable of detecting the minor contributor's allele in home‐made DNA mixture with sensitivities from 1:100 to 1:1000 with a total of 1 –10 ng input DNA. Moreover, this system successfully typed the degraded DNA whether it came from the single source or mixture samples. In Chinese population, the system showed an average informative value of 0.293 and combined informative value of 0.998363862. Our results demonstrated that DIP‐SNPs may serve as a valuable tool in detection of UDM in forensic medicine.     

Global analysis of single nucleotide polymorphisms in the exons of human deoxyribonuclease I-like 1 and 2 genes

Several SNPs in the deoxyribonuclease I-like 1 (DNase 1L1) and DNase 1L2 were investigated. In the present study, the genotype distributions of three synonymous SNPs (V59V, rs1050095; P67P, rs1130929; A277A, rs17849495) in the DNase 1L1 gene and four non-synonymous SNPs, V122I (rs34952165), Q170H (rs6643670), and D227A (rs5987256) in the DNase 1L1 gene, as well as D197A (rs62621282) in the DNase 1L2 gene were investigated in 13 populations. In all the populations, no variation was found in four SNPs (V59V, Q170H, D227A, and A277A) in DNASE1L1 or in D197A in DNASE1L2. As for V122I, only the German population showed a low degree of polymorphism. The SNP V122I in DNASE1L1 was monoallelic for the G-allele in all of the Asian and African populations examined, with no polymorphism being evident. Since the A-allele in SNP V122I was distributed in only the Caucasian populations, not in the other ethnic groups, it was confirmed that the A-allele in SNP V122I was Caucasian-specific. On the other hand, only P67P in DNASE1L1 was polymorphic among three synonymous SNPs. The effect of nucleotide substitution corresponding to polymorphic SNP P67P on DNase 1L1 activity was examined: the corresponding nucleotide substitution in polymorphic SNP P67P has little effect on the DNase activity.     

A new approach to detect a set of SNP-SNP markers: Combining ARMS-PCR with SNaPshot technology

Microhaplotypes are a new promising type of forensic genetic marker. Without the interference of stutter and high mutation rates as for STRs, and with short amplification lengths and a higher degree of polymorphism than single SNP, microhaplotypes composed of two SNPs, SNP–SNP, have a strong application potential. Currently, the most common method to detect microhaplotypes is massive parallel sequencing. However, the cost and extensive use of instruments limit its wide application in forensic laboratories. In this study, we screened 23 new SNP–SNP loci and established a new detection method by combining a multiplex amplification refractory mutation system-based PCR (ARMS-PCR) and SNaPshot technology based on CE. First, we introduced an additional deliberate mismatch at the antepenultimate base from the 3′ end of primers when designing ARMS-PCR for SNP 1 (the first SNP of the SNP–SNP). Then, single base extension primers for SNaPshot assay were designed next to the position of SNP 2 (the second SNP). Finally, 15 loci were successfully built into four panels and these loci showed a relatively high level of polymorphism in the Southwest Chinese Han population. All the loci had an average probability of informative genotypes (I value) of 0.319 and a combined discrimination power of 0.999999999. Therefore, this new detection system will provide a valuable supplement to current methods.     

Confirmation that SNPs in the high mobility group‐A2 gene (HMGA2) are associated with adult height in the Japanese population; wide‐ranging population survey of height‐related SNPs in HMGA2

Adult height is a highly heritable trait in that multiple genes are involved. Recent genome‐wide association studies have identified a novel single‐nucleotide polymorphism (SNP) rs1042725 in the high mobility group‐A2 gene (HMGA2) and shown it to be associated with human height in Caucasian populations. We performed a replication study to examine the associations between SNPs in HMGA2 and adult height in the Japanese population based on autopsy cases. Although we could not confirm a significant association between rs1042725 in HMGA2 and adult height, another SNP, rs7968902, in the gene achieved significance for its association in the same populations, and the effect was the same as that documented previously. These findings permit us to conclude that the SNPs in HMGA2 are common variants influencing human height across different populations. Moreover, a worldwide population study of these SNPs using 14 different populations including Asians, Africans and Caucasians demonstrated that both haplotypes and genotypes for three height‐related SNPs (rs1042725, rs7968682 and rs7968902) in HMGA2 were distributed in an ethnicity‐dependent manner. This information will be useful for clarifying the genetic basis of human height.     

Intramolecular controls for the generation of clinical-quality SNP genotypes and the assessment of normal heterozygote imbalance

We present an inventive method that generates intramolecular controls for SNP analysis, termed Mirror SNPs. Using the ovine diseases of callipyge and scrapie as examples, we describe the PCR-driven production of balanced heterozygote copies of the various SNPs implicated in these diseases. In the absence of a callipyge-positive control DNA, we generated a balanced heterozygote Mirror SNP that represents both the wild-type and mutant forms of the causal polymorphism. Simultaneous analysis of this artificial Mirror SNP and the Real (target) SNP was used to prove the absence of the mutant form of the nucleotide at the Real SNP position in tested samples. Scrapie susceptibility was assessed using a PCR-driven system which generated four separate Mirror SNPs, and these enabled the confirmation of an apparent departure from 'balanced' heterozygote appearance at Real SNPs tested. Mirror SNP technology is generic and will enable the accurate assessment of rare and medically important SNP variants, more accurate frequency determinations, and the potential assessment of SNPs in 'mixed template' samples common in forensic analyses.     

Using propensity score adjustment method in genetic association studies

BackgroundThe statistical tests for single locus disease association are mostly under-powered. If a disease associated causal single nucleotide polymorphism (SNP) operates essentially through a complex mechanism that involves multiple SNPs or possible environmental factors, its effect might be missed if the causal SNP is studied in isolation without accounting for these unknown genetic influences. In this study, we attempt to address the issue of reduced power that is inherent in single point association studies by accounting for genetic influences that negatively impact the detection of causal variant in single point association analysis. In our method we use propensity score (PS) to adjust for the effect of SNPs that influence the marginal association of a candidate marker. These SNPs might be in linkage disequilibrium (LD) and/or epistatic with the target-SNP and have a joint interactive influence on the disease under study. We therefore propose a propensity score adjustment method (PSAM) as a tool for dimension reduction to improve the power for single locus studies through an estimated PS to adjust for influence from these SNPs while regressing disease status on the target-genetic locus. The degree of freedom of such a test is therefore always restricted to 1.ResultsWe assess PSAM under the null hypothesis of no disease association to affirm that it correctly controls for the type-I-error rate (<0.05). PSAM displays reasonable power (>70%) and shows an average of 15% improvement in power as compared with commonly-used logistic regression method and PLINK under most simulated scenarios. Using the open-access multifactor dimensionality reduction dataset, PSAM displays improved significance for all disease loci. Through a whole genome study, PSAM was able to identify 21 SNPs from the GAW16 NARAC dataset by reducing their original trend-test p-values from within 0.001 and 0.05 to p-values less than 0.0009, and among which 6 SNPs were further found to be associated with immunity and inflammation.ConclusionsPSAM improves the significance of single-locus association of causal SNPs which have had marginal single point association by adjusting for influence from other SNPs in a dataset. This would explain part of the missing heritability without increasing the complexity of the model due to huge multiple testing scenarios. The newly reported SNPs from GAW16 data would provide evidences for further research to elucidate the etiology of rheumatoid arthritis. PSAM is proposed as an exploratory tool that would be complementary to other existing methods. A downloadable user friendly program, PSAM, written in SAS, is available for public use.     

Developing a new nonbinary SNP fluorescent multiplex detection system for forensic application in China

Nonbinary single‐nucleotide polymorphisms (SNPs) are potential forensic genetic markers because their discrimination power is greater than that of normal binary SNPs, and that they can detect highly degraded samples. We previously developed a nonbinary SNP multiplex typing assay. In this study, we selected additional 20 nonbinary SNPs from the NCBI SNP database and verified them through pyrosequencing. These 20 nonbinary SNPs were analyzed using the fluorescent‐labeled SNaPshot multiplex SNP typing method. The allele frequencies and genetic parameters of these 20 nonbinary SNPs were determined among 314 unrelated individuals from Han populations from China. The total power of discrimination was 0.9999999999994, and the cumulative probability of exclusion was 0.9986. Moreover, the result of the combination of this 20 nonbinary SNP assay with the 20 nonbinary SNP assay we previously developed demonstrated that the cumulative probability of exclusion of the 40 nonbinary SNPs was 0.999991 and that no significant linkage disequilibrium was observed in all 40 nonbinary SNPs. Thus, we concluded that this new system consisting of new 20 nonbinary SNPs could provide highly informative polymorphic data which would be further used in forensic application and would serve as a potentially valuable supplement to forensic DNA analysis.     

Fluorescence polarization for single nucleotide polymorphism genotyping

Single nucleotide polymorphisms (SNPs) are the most abundant variations in the human genome and have become the primary markers for genetic studies for mapping and identifying susceptible genes for complex diseases. Methods that genotype SNPs quickly and economically are of high values for these studies because they require a large amount of genotyping. Fluorescence polarization (FP) is a robust technique that can detect products without separation and purification and it has been applied for SNP genotyping. In this article the applications of FP in SNP genotyping are reviewed and one of the methods, the FP-TDI assay, is discussed in details. It is hoped that readers could get useful information for the applications of FP in SNP genotyping and some insights of the FP-TDI assay.     

Generating SNP barcode to evaluate SNP-SNP interaction of disease by particle swarm optimization

Genome-wide association analysis involved many single-nucleotide polymorphisms (SNPs) data is challenging mathematically and computationally. Hence, we propose the odds ratio-based discrete binary particle swarm optimization (OR-DBPSO) method that uses the OR as a new quantitative measure of disease risk among many SNP combinations with genotypes called "SNP barcode". DBPSO are applied to generate SNP barcode, which computes the maximal difference of occurrence between the case and control groups, to predict disease susceptibility such as osteoporosis. Different SNP barcode patterns may occur several times in either low or high bone mineral density (BMD) groups. Our results showed that a DBPSO can effectively identify a specific SNP barcode with an optimized fitness value. SNP barcodes with a low fitness value will naturally be discarded from the population. A representative SNP barcode with a variable number of SNPs is processed to OR analysis to determine the maximum difference between the low and high BMD groups in statistics manner. Therefore, this paper introduces a powerful procedure to analyze disease-associated SNP-SNP interaction in genome-wide genes.     

Seven nonsynonymous SNPs in the gene encoding human deoxyribonuclease II may serve as a functional SNP potentially implicated in autoimmune dysfunction

Many nonsynonymous SNPs in the human DNase II gene (DNASE2), potentially relevant to autoimmunity in conditions such as rheumatoid arthritis, have been identified, but only limited population data are available and no studies have evaluated whether such SNPs are functional. Genotyping of all the 15 nonsynonymous human DNase II SNPs was performed in three ethnic groups including 16 different populations using the PCR‐restriction fragment length polymorphism technique. A series of constructs corresponding to each SNP was examined. Fifteen nonsynonymous SNPs in the gene, except for p.Val206Ile in a Korean population, exhibited a mono‐allelic distribution in all of the populations. On the basis of alterations in the activity levels resulting from the corresponding amino acid substitutions, four activity‐abolishing and five activity‐reducing SNPs were confirmed to be functional. The amino acid residues in activity‐abolishing SNPs were conserved in animal DNase II. All the nonsynonymous SNPs that affected the catalytic activity of human DNase II showed extremely low genetic heterogeneity. However, a minor allele of seven SNPs producing a loss‐of‐function or extremely low activity‐harboring variant could serve as a genetic risk factor for autoimmune dysfunction. These functional SNPs in DNASE2 may have clinical implications in relation to the prevalence of autoimmune diseases.     

Fractal structure of asphaltene aggregates

A photographic technique coupled with image analysis was used to measure the size and fractal dimension of asphaltene aggregates formed in toluene-heptane solvent mixtures. First, asphaltene aggregates were examined in a Couette device and the fractal-like aggregate structures were quantified using boundary fractal dimension. The evolution of the floc structure with time was monitored. The relative rates of shear-induced aggregation and fragmentation/restructuring determine the steady-state floc structure. The average floc structure became more compact or more organized as the floc size distribution attained steady state. Moreover, the higher the shear rate is, the more compact the floc structure is at steady state. Second, the fractal dimensions of asphaltene aggregates were also determined in a free-settling test. The experimentally determined terminal settling velocities and characteristic lengths of the aggregates were utilized to estimate the 2D and 3D fractal dimensions. The size-density fractal dimension (D(3)) of the asphaltene aggregates was estimated to be in the range from 1.06 to 1.41. This relatively low fractal dimension suggests that the asphaltene aggregates are highly porous and very tenuous. The aggregates have a structure with extremely low space-filling capacity.     

SSCP markers provide a useful alternative to microsatellites in genotyping and estimating genetic diversity in populations and germplasm collections of plant specialty crops

For well-studied plant species with whole genome sequence or extensive EST data, SNP markers are the logical choice for both genotyping and whole genome association studies. However, SNP markers may not address the needs of researchers working on specialty crops with limited available genomic information. Microsatellite markers have been frequently employed due to their robustness, but marker development can be difficult and may result in few polymorphic markers. SSCP markers, such as microsatellites, are PCR-based and scored by electrophoretic mobility but, because they are based on SNPs rather than length differences, occur more frequently and are easier to develop than microsatellites. We have examined how well correlated the estimation of genetic diversity and genetic distance are in a population or germplasm collection when measured by 13 highly polymorphic microsatellite markers or 20 SSCP markers. We observed a significant correlation in pairwise genetic distances of 82 individuals in an international cacao germplasm collection (Mantel test Rxy=0.59, p<0.0001 for 10 000 permutations). Both sets of markers could distinguish each individual in the population. These data provide strong support for the use of SSCP markers in the genotyping of plant species where development of microsatellites would be difficult or expensive.     

Simultaneous detection of single nucleotide polymorphisms and copy number variations in the CYP2D6 gene by multiplex polymerase chain reaction combined with capillary electrophoresis

CYP2D6 (cytochrome P450 2D6) is one of the most important enzymes involved in drug metabolism, and CYP2D6 gene variants may cause toxic effects of therapeutic drugs or treatment failure. In this research, a rapid and simple method for genotyping the most common mutant alleles in the Asian population (CYP2D6*1/*1, CYP2D6*1/*10, CYP2D6*10/*10, CYP2D6*1/*5, CYP2D6*5/*10, and CYP2D6*5/*5) was developed by allele-specific polymerase chain reaction (AS-PCR) combined with capillary electrophoresis (CE). We designed a second mismatch nucleotide next to the single nucleotide polymorphism (SNP) site in allele-specific primers to increase the difference in PCR amplification. Besides, we established simulation equations to predict the CYP2D6 genotypes by analyzing the DNA patterns in the CE chromatograms. The multiplex PCR combined with CE method was applied to test 50 patients, and all of the test results were compared with the DNA sequencing method, long-PCR method and real-time PCR method. The correlation of the analytical results between the proposed method and other methods were higher than 90%, and the proposed method is superior to other methods for being able to simultaneous detection of SNPs and copy number variations (CNV). Furthermore, we compared the plasma concentration of aripiprazole (a CYP2D6 substrate) and its major metabolites with the genotype of 25 patients. The results demonstrate the proposed genotyping method is effective for estimating the activity of the CYP2D6 enzyme and shows potential for application in personalized medicine. Similar approach can be applied to simultaneous detection of SNPs and CNVs of other genes.     

A standard protocol for single nucleotide primer extension in the human genome using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry

Analysis of single nucleotide polymorphisms (SNPs) has become an increasingly important area of research, with numerous applications in medical genetics, population genetics, forensic science, and agricultural biotechnology. Large-scale SNP analyses require the development of methodologies that are economical, flexible, accurate and capable of automation. Primer extension in conjunction with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS) is currently emerging as a potential method for high-throughput SNP genotyping. We have evaluated a number of published primer extension methods and refined a simple and robust protocol to analyze human autosomal disease-causing mutations and population genetic markers on the Y-chromosome. Twelve different variant sites were examined, and homozygotes, heterozygotes and hemizygotes were accurately typed. A 100% concordance was observed between SNP genotypes obtained using the MALDI-TOFMS technique and alternative genotyping methods, such as restriction fragment length polymorphism (RFLP) assays and denaturing high-performance liquid chromatography (DHPLC). Since multiple polymorphisms can be detected in single reactions, the method provides a cost-effective approach for SNP analysis. The protocol is also extremely flexible (able to accommodate new markers) and can be adapted to a number of platforms without the use of commercial kits.     

Relation between aggregation kinetics and the structure of kaolinite aggregates

Dynamic and static light scattering were applied to the determination of the stability ratio and fractal dimension of kaolinite (KGa-2) at different kaolinite or/and electrolyte concentrations at pH 9.5. Dynamic light scattering was used to measure the kinetics of early stage aggregation to determine the stability ratio, W, as well as the cluster sizes which determine the fractal regime. Static light scattering was used to measure the fractal dimension, D(f). Results show that the two classes of "universality" (Lin et al. Nature 1989, 339, 360) characterizing the diffusion- and reaction-limited regimes of cluster-cluster aggregation do apply to colloidal kaolinite as limit cases when W approximately 1 or W > 100, respectively. In the intermediate regime where 5 < W < 100, the growth of the aggregate radius showed a power-law behavior similar to diffusion-limited cluster aggregation. For the intermediate aggregation regime, a scaling relation between fractal dimension and stability ratio, reflecting a continuous increase in particle packing density in the aggregate as the sticking probability of particles was reduced, was demonstrated.     

Exploring of tri-allelic SNPs using pyrosequencing and the SNaPshot methods for forensic application

Tri-allelic single nucleotide polymorphisms (SNPs) are potential forensic markers for DNA analysis. Currently, only a limited number of tri-allelic SNP loci have been proved to be fit for forensic application. In this study, we aimed to develop an effective method to select and genotype tri-allelic SNPs based on both Pyrosequencing (PSQ) and the SNaPshot methods. 50 candidate SNPs were chosen from NCBI's dbSNP database and were analyzed by PSQ. The results revealed that 20 SNPs were tri-allelic and were located on 16 autosomal chromosomes. Then 20 SNP loci were combined in one multiplex polymerase chain reaction to develop a single base extension (SBE)-based SNP-typing assay. A total of 100 unrelated Chinese individuals were genotyped by this assay and allele frequencies were estimated. The total discrimination power was 0.999999999975 and the cumulative probability of exclusion was 0.9937. These data demonstrated that the strategy is a rapid and effective method for seeking and typing tri-allelic SNPs. In addition, the 20 tri-allelic SNP multiplex typing assay may be used to supplement paternity testing and human identification.