HaplotypeCN: Copy Number Haplotype Inference with Hidden Markov Model and Localized Haplotype Clustering

Copy number variation (CNV) has been reported to be associated with disease and various cancers. Hence, identifying the accurate position and the type of CNV is currently a critical issue. There are many tools targeting on detecting CNV regions, constructing haplotype phases on CNV regions, or estimating the numerical copy numbers. However, none of them can do all of the three tasks at the same time. This paper presents a method based on Hidden Markov Model to detect parent specific copy number change on both chromosomes with signals from SNP arrays. A haplotype tree is constructed with dynamic branch merging to model the transition of the copy number status of the two alleles assessed at each SNP locus. The emission models are constructed for the genotypes formed with the two haplotypes. The proposed method can provide the segmentation points of the CNV regions as well as the haplotype phasing for the allelic status on each chromosome. The estimated copy numbers are provided as fractional numbers, which can accommodate the somatic mutation in cancer specimens that usually consist of heterogeneous cell populations. The algorithm is evaluated on simulated data and the previously published regions of CNV of the 270 HapMap individuals. The results were compared with five popular methods: PennCNV, genoCN, COKGEN, QuantiSNP and cnvHap. The application on oral cancer samples demonstrates how the proposed method can facilitate clinical association studies. The proposed algorithm exhibits comparable sensitivity of the CNV regions to the best algorithm in our genome-wide study and demonstrates the highest detection rate in SNP dense regions. In addition, we provide better haplotype phasing accuracy than similar approaches. The clinical association carried out with our fractional estimate of copy numbers in the cancer samples provides better detection power than that with integer copy number states.     

Genetic Structure of Japanese Scallop Population in Hokkaido, Analyzed by Mitochondrial Haplotype Distribution

To examine the genetic structure of Japanese scallop populations (Mizuhopecten yessoensis) in Hokkaido prefecture, Japan, and compare it with those in the Aomori prefecture, we applied a method for lineage analysis based on sequence variation in a mitochondrial DNA segment (NcR2). After showing that there was a low probability of doubly uniparental inheritance of mitochondrial DNA in the scallop, we sequenced the NcR2 regions of 914 individuals from 15 populations (13 in Hokkaido and 2 in Aomori). In total, 103 different haplotypes were detected. Results of homogeneity tests for pairwise populations and the fixation indices indicated that significant heterogeneity (P < 0.0005) and structuring (pairwise fixation index F_ST = 0.1606–0.4444, P = 0.0000; fixation index among groups F_CT = 0.1549, P = 0.0078) could be inferred between the Hokkaido and Aomori groups, but not among populations within the groups. Moreover, heterogeneity of the haplotype distribution between populations of the 1980s and 1990s or 2000s at the 4 culturing areas was not observed (P > 0.05), and the haplotype diversity between them was not significant (P = 0.05), suggesting that the culture operations had not imparted a significant effect on the genetic structure during these periods.     

Models for Haplotype Evolution in a Nonstationary Population

Haplotype mapping has emerged in the past few years as a powerful tool for the fine mapping of disease genes. It is typically carried out on a sample of affected individuals from a population isolate. If the chromosome neighborhood of a disease gene is saturated with markers, then each new mutation in the population or existing mutation introduced by a population founder exhibits a unique haplotype signature at the time of its introduction. Partial disruptions of these signatures by recombination can be visualized in affecteds and provide important clues to the location of the disease gene. The current paper models haplotype evolution with the intention of clarifying the most favorable circumstances for haplotype mapping. Comparisons with linkage mapping are stressed. For dominant diseases, both deterministic and stochastic models are suggested. Numerical examples based on Finnish population parameters illustrate the general theory in the presence of the complications of selection, mutation, and slow, exponential growth of the isolate.     

Haplotype Inference for Population Data with Genotyping Errors

Inference of haplotypes is important in genetic epidemiology studies. However, all large genotype data sets have errors due to the use of inexpensive genotyping machines that are fallible and shortcomings in genotyping scoring softwares, which can have an enormous impact on haplotype inference. In this article, we propose two novel strategies to reduce the impact induced by genotyping errors in haplotype inference. The first method makes use of double sampling. For each individual, the “GenoSpectrum” that consists of all possible genotypes and their corresponding likelihoods are computed. The second method is a genotype clustering algorithm based on multi‐genotyping data, which also assigns a “GenoSpectrum” for each individual. We then describe two hybrid EM algorithms (called DS‐EM and MG‐EM) that perform haplotype inference based on “GenoSpectrum” of each individual obtained by double sampling and multi‐genotyping data. Both simulated data sets and a quasi real‐data set demonstrate that our proposed methods perform well in different situations and outperform the conventional EM algorithm and the HMM algorithm proposed by Sun, Greenwood, and Neal (2007, Genetic Epidemiology 31 , 937–948) when the genotype data sets have errors.     

山东荣成人群线粒体DNA多态性研究

人类线粒体DNA（mtDNA)COⅡ/tRNA^lys区有两个9－bp(CCCCCTCTA)的串联重复序列,此重复序列中一个重复单位的缺失,在亚态地区人群中很普遍。对210名山东荣成人的mtDNA COⅡ/tRNA^lys区的9-bp 缺失情况进行了检测,并从中随机选取95个样本,利用PCR－RFLP法对另外6个区进行了多态性分析,以确定其单位型。结果表明,荣成人9－bp缺失频率为12.4%,相对于已检测的中国其他群体,此缺失频率处于中等水平。同时多态性分析也表明在95个被检测对象中存在27种不同的单倍型。此外还发现了两个未报道过的新酶切位点,序列分析表明是由点突变造成的。     

GABRB2 Haplotype Association with Heroin Dependence in Chinese Population

Substance dependence is a frequently observed comorbid disorder in schizophrenia, but little is known about genetic factors possibly shared between the two psychotic disorders. GABRB2, a schizophrenia candidate gene coding for GABA_A receptor β₂ subunit, is examined for possible association with heroin dependence in Han Chinese population. Four single nucleotide polymorphisms (SNPs) in GABRB2, namely rs6556547 (S1), rs1816071 (S3), rs18016072 (S5), and rs187269 (S29), previously associated with schizophrenia, were examined for their association with heroin dependence. Two additional SNPs, rs10051667 (S31) and rs967771 (S32), previously associated with alcohol dependence and bipolar disorder respectively, were also analyzed. The six SNPs were genotyped by direct sequencing of PCR amplicons of target regions for 564 heroin dependent individuals and 498 controls of Han Chinese origin. Interestingly, it was found that recombination between the haplotypes of all-derived-allele (H1; OR = 1.00) and all-ancestral-allele (H2; OR = 0.74) at S5-S29 junction generated two recombinants H3 (OR = 8.51) and H4 (OR = 5.58), both conferring high susceptibility to heroin dependence. Additional recombination between H2 and H3 haplotypes at S1-S3 junction resulted in a risk-conferring haplotype H5 (OR = 1.94x10⁹). In contrast, recombination between H1 and H2 haplotypes at S3-S5 junction rescued the risk-conferring effect of recombination at S5-S29 junction, giving rise to the protective haplotype H6 (OR = 0.68). Risk-conferring effects of S1-S3 and S5-S29 crossovers and protective effects of S3-S5 crossover were seen in both pure heroin dependent and multiple substance dependence subgroups. In conclusion, significant association was found with haplotypes of the S1-S29 segment in GABRB2 for heroin dependence in Han Chinese population. Local recombination was an important determining factor for switching haplotypes between risk-conferring and protective statuses. The present study provide evidence for the schizophrenia candidate gene GABRB2 to play a role in heroin dependence, but replication of these findings is required.     

HLA-G UTR Haplotype Conservation in the Malian Population: Association with Soluble HLA-G

The HLA-G molecule plays an important role in immunomodulation. In a previous study carried out on a southern French population our team showed that HLA-G haplotypes, defined by SNPs in the coding region and specific SNPs located in 5′URR and 3′UTR regulatory regions, are associated with differential soluble HLA-G expression (sHLA-G). Furthermore, the structure of these HLA-G haplotypes appears to be conserved in geographically distant populations.The aim of our study is to confirm these expectations in a sub-Saharan African population and to explore additional factors, such as HLA-A alleles, that might influence sHLA-G expression.DNA and plasma samples were collected from 229 Malians; HLA-G and HLA-A genotyping were respectively performed by the Snap Shot® method and by Luminex™ technology. sHLA-G dosage was performed using an ELISA kit. HLA-G and HLA-A allelic and haplotypic frequencies were estimated using an EM algorithm from the Gene[Rate] program. Associations between genetic and non genetic parameters with sHLA-G were performed using a non-parametric test with GRAPH PAD Prism 5.Our results reveal a good conservation of the HLA-G UTR haplotype structure in populations with different origins and demographic histories. These UTR haplotypes appear to be involved in different sHLA-G expression patterns. Specifically, the UTR-2 haplotype was associated with low sHLA-G levels, displaying a dominant negative effect. Furthermore, an allelic effect of both HLA-G and HLA-A, as well as non genetic parameters, such as age and gender possibly linked to osteogenesis and sexual hormones, also seem to be involved in the modulation of sHLA-G.These data suggest that further investigation in larger cohorts and in populations from various ethnical backgrounds is necessary not only to detect new functional polymorphism in HLA-G regulatory regions, but also to reveal the extent of biological phenomena that influence sHLA-G secretion and this might therefore have an impact on transplantation practice.     

Detecting Local Haplotype Sharing and Haplotype Association

A novel haplotype association method is presented, and its power is demonstrated. Relying on a statistical model for linkage disequilibrium (LD), the method first infers ancestral haplotypes and their loadings at each marker for each individual. The loadings are then used to quantify local haplotype sharing between individuals at each marker. A statistical model was developed to link the local haplotype sharing and phenotypes to test for association. We devised a novel method to fit the LD model, reducing the complexity from putatively quadratic to linear (in the number of ancestral haplotypes). Therefore, the LD model can be fitted to all study samples simultaneously, and, consequently, our method is applicable to big data sets. Compared to existing haplotype association methods, our method integrated out phase uncertainty, avoided arbitrariness in specifying haplotypes, and had the same number of tests as the single-SNP analysis. We applied our method to data from the Wellcome Trust Case Control Consortium and discovered eight novel associations between seven gene regions and five disease phenotypes. Among these, GRIK4, which encodes a protein that belongs to the glutamate-gated ionic channel family, is strongly associated with both coronary artery disease and rheumatoid arthritis. A software package implementing methods described in this article is freely available at http://www.haplotype.org.     

Haplotype Parsing

While the shared consensus genetic sequence of our species contains a great deal of information about our common biology, there is also much to be learned from the subtle genetic variations across our species. These variations are believed to be generally of little or no direct functional significance and predominantly reflect the chance accumulation of small genetic changes since our emergence as a species. Therefore, they carry little useful information when observed in a single individual. When tallied across a whole population though, these chance mutations can teach us a great deal about our evolutionary history and the patterns of inheritance in particular individuals. In particular, frequently observed patterns of single nucleotide polymorphisms (SNPs) in a population can identify segments of chromosome that have been passed down largely intact through long stretches of our evolution. Finding these frequently conserved chromosomal segments, or haplotypes, and developing methods to identify haplotype patterns in particular individuals, will in turn help us to identify those particular segments that carry genetic factors influencing risk for many common human diseases. To make the best use of this data, we will need to develop new models for the encoding of information in genome variations--the "language of genetic variation"--and new algorithms for fitting datasets to those models. This article surveys past work by the author and colleagues on this problem, utilising computational methods for locating frequent patterns in haploid sequence data, and "parsing" sequences so as to optimally explain them given the knowledge of the general population structure. The author's recent work in this area has been compiled into a set of computational tools available at http://www-2.cs.cmu.edu/~russells/software/hapmotif.html.     

广州地区106例汉族人群MICA和MICB微卫星基因座 单体型研究

利用聚合酶链反应和荧光（6－FAM）自动化检测技术对广东地区汉族106例无亲缘关系样本进行MICA基因外显子5和MICB基因内含子1微卫星基因座多态性及其单体型分布调查。根据群体资料估算两者间的单体型频率、连锁不平衡参数、相对连锁不平衡参数。结果显示,广州地区汉族人群MICA和MICB微卫星基因座基因型分布符合Hardy-Weinberg平衡法则,共检出MICA微卫星基因座 5个等位基因, MICB微卫星基因座14个等位基因。其中MICA A5基因频率最高（0.2877）,A4基因频率最低（0.1321）。MICB CA14等位基因频率最高（0.3255）,CA19、CA28等位基因频率最低（0.0047）,未检出CA27。21种MICA－MICB单体型频率大于1%(连锁不平衡参数>0), 其中单体型A5－CA14 (16.73%), A5.1－CA18 (8.75%), A4－CA26(3.76%),A9－CA15(3.66%)和A6－CA21(2.61%)为强连锁常见单体型（χ2>3.84, P<0.05）。广州地区汉族人群MICA和MICB微卫星基因座多态性和单体型分布有其自身特点,MICA和MICB微卫星基因座适合做为遗传标志,用于人类学、遗传疾病基因连锁分析、法医学亲子鉴定和个体识别等研究领域。Abstract: This study is to investigate genetic polymorphisms and haplotypes of microsatellite locus in the exon 5 of the MICA gene and intron 1 of the MICB gene based on 106 samples of Guangzhou Han Population by polymerase chain reaction and fluorescent technique (6-FAM). The corresponding haplotype frequencies, linkage disequilibria values and relative linkage disequilibria values were estimated based on population data. The results show that the genotype distributions of MICA and MICB microsatellite meet Hardy-Weinberg equilibrium in Guangdong Han population. In total, 5 alleles of MICA microsatellite locus and 14 alleles of MICB microsatellite locus were observed. MICA A5 was the most common allele (0.2877), whereas A4 was the least popular one (0.1321). MICB CA14 was the most common allele (0.3255), and CA19 and CA28 were the least popular ones (0.0047). CA27 was not observed. Twenty-one kinds of MICA-MICB haplotypes occurred at frequencies of more than 1% (linkage disequilibria value>0). The common MICA-MICB haplotypes were A5-CA14（16.73％）, A5.1- CA18 (8.75%), A4- CA26(3.76%),A9-CA15(3.66%) and A6-CA21(2.61%)（χ2>3.84, P<0.05）, and they were strong linkage disequilibria. The polymorphisms and haplotypes distributions of MICA and MICB microsatellite locus in Guangzhou Han population have their own genetic characteristics. The microsatellite locus of the exon5 of the MICA gene and intron 1 of the MICB gene could be used as the genetic markers in the studies of anthropology, linkage analysis of genetic disease genes, individual identification and paternity test in forensic medicine.     

Single-Particle Spectroscopy of Supported Silver Clusters on Silicon: Substrate Effects on Localized Surface Plasmons

The presence of a surrounding medium strongly affects the spectral properties of localized surface plasmons at metallic nanoparticles. Vice versa, plasmonic resonances have large impact on the electric polarization in a surrounding or supporting material. For applications, e.g., in light-converting devices, the coupling of localized surface plasmons with polarizations in semiconducting substrates is of particular importance. Using photoemission electron microscopy with tunable laser excitation, we perform single-particle spectroscopy of silver nanoclusters directly grown on Si(100). Two distinct localized surface plasmon modes are observed as resonances in the two-photon photoemission signals from individual silver clusters. The strengths of these resonances strongly depend on the polarization of the exciting electric field, which allows us to assign them to plasmon modes with polarizations parallel and perpendicular, respectively, to the supporting silicon substrate. Our mode assignment is supported by simulations which provide insight into the mutual interaction of charge oscillations at the particle surface with electric polarizations at the silver/silicon interface.

     

Laplacian Eigenfunctions Learn Population Structure

Principal components analysis has been used for decades to summarize genetic variation across geographic regions and to infer population migration history. More recently, with the advent of genome-wide association studies of complex traits, it has become a commonly-used tool for detection and correction of confounding due to population structure. However, principal components are generally sensitive to outliers. Recently there has also been concern about its interpretation. Motivated from geometric learning, we describe a method based on spectral graph theory. Regarding each study subject as a node with suitably defined weights for its edges to close neighbors, one can form a weighted graph. We suggest using the spectrum of the associated graph Laplacian operator, namely, Laplacian eigenfunctions, to infer population structure. In simulations and real data on a ring species of birds, Laplacian eigenfunctions reveal more meaningful and less noisy structure of the underlying population, compared with principal components. The proposed approach is simple and computationally fast. It is expected to become a promising and basic method for population genetics and disease association studies.     

Nucleotide Diversity and Haplotype Structure of the Human Angiotensinogen Gene in Two Populations

Variation in the angiotensinogen gene, AGT, has been associated with variation in plasma angiotensinogen levels. In addition, the T235M polymorphism in the AGT product is associated with an increased risk of essential hypertension in multiple populations, making AGT a good example of a quantitative-trait locus underlying susceptibility to a common disease. To better understand genetic variation in AGT, we sequenced a 14.4-kb genomic region spanning the entire AGT and identified 44 single-nucleotide polymorphisms (SNPs). Forty-two SNPs were observed both in 88 white and in 77 Japanese unselected subjects. Six major haplotypes accounted for most of the variation in this region, indicating less allelic complexity than in many other genomic regions. Although the two populations were found to share all of the major AGT haplotypes, there were substantial differences in haplotype frequencies. Pairwise linkage disequilibrium (LD), measured by the D', r(2), and d(2) statistics, demonstrated a general pattern of decline with increasing distance, but, as expected in a small genomic region, individual LD values were highly variable. LD between T235M and each of the other 39 SNPs was assessed in order to model the usefulness of LD to detect a disease-associated mutation. Among the Japanese subjects, 13 (33%) of the SNPs had r(2) values >0.1, whereas this statistic was substantially higher for the white subjects (occurring in 35/39 [90%]). LD between a hypertension-associated promoter mutation, A-6G, and 39 SNPs was also measured. Similar results were obtained, with 33% of the SNPs showing r(2)>0.1 in the Japanese subjects and 92% of the SNPs showing r(2)>0.1 in the white subjects. This difference, which occurs despite an overall similarity in LD patterns in the two populations, reflects a much higher frequency of the M235-associated haplotype in the white sample. These results have important implications for the usefulness of LD approaches in the mapping of genes underlying susceptibility to complex diseases.     

中国汉族人群MASP2基因单核苷酸多态性的连锁不平衡和单体型分析

目的：探讨中国汉族人群甘露糖结合凝集素相关的丝氨酸蛋白酶（MASP2）基因单核苷酸多态性（SNP）的连锁不平衡和单体型,为相关研究提供更为详实的数据。方法：选取30例广州汉族无关个体,用重新测序的方法进行MASP2基因SNP的发掘,构建连锁不平衡（LD）模式,与HapMap公布的其他4个人群数据相比,并选择4个标签SNP（tagSNP）在232例北京汉族个体和291例广州汉族个体中分析MASP2的多态性和单体型。结果和结论：对MASP2的重测序共检测出16个SNP,LD分析显示来自中国不同地区人群的LD模式基本相同,与来自美国和日本的人群存在差异;北京和广州地区汉族人群tagSNP的等位基因和基因型频率分布不存在显著差异。     

ADRB2 Haplotype Is Associated With Glucose Tolerance and Insulin Sensitivity in Obese Postmenopausal Women

The β₂‐adrenergic receptor (ADRB2) mediates obesity, cardiorespiratory fitness, and insulin resistance. We examined the hypothesis that ADRB2 Arg16Gly‐Gln27Glu haplotype is associated with body composition, glucose tolerance, and insulin sensitivity in obese, postmenopausal women. Obese (>35% body fat), postmenopausal (age 45–75 years) women (n = 123) underwent genotyping, dual‐energy X‐ray absorptiometry, and computed tomography scans, exercise testing (VO_2max), 2‐h oral glucose tolerance tests (OGTTs), and hyperinsulinemic‐euglycemic clamps (80 mU/m²/min). Analysis of covariance (ANCOVA) tested for differences among haplotypes, with race, % body fat, and VO_2max as covariates. We found that ADRB2 haplotype was independently associated with % body fat, abdominal fat distribution, VO_2max, insulin sensitivity (M/ΔInsulin), and glucose tolerance (ANOVA, P < 0.05 for all). Women homozygous for Gly16–Gln27 haplotype had the highest % body fat (52.7 ± 1.9%), high abdominal fat, low M/ΔInsulin (0.49 ± 0.08 mg/kg/min/pmol/l/10²), and impaired glucose tolerance (IGT) during an OGTT (G₁₂₀ = 10.2 ± 0.9 mmol/l). Women homozygous for Gly16–Glu27 haplotype also had low M/ΔInsulin (0.51 ± 0.05 mg/kg/min/pmol/l/10²) and IGT (G₁₂₀ = 8.2 ± 0.7 mmol/l). Subjects with Arg16–Gln27/Gly16–Gln27 haplotype combination had the highest VO_2max (1.84 ± 0.07 l/min) and M/ΔInsulin (0.7 ± 0.04 mg/kg/min/pmol/l/10²), and normal glucose tolerance (G₁₂₀ = 6.4 ± 0.4 mmol/l), despite being obese. These data show associations of the ADRB2 Arg16Gly‐Gln27Glu haplotype with VO_2max and body composition, and an independent association with glucose metabolism, which persists after controlling for body composition and fitness. This suggests that ADRB2 haplotypes may mediate insulin action, glucose tolerance, and potentially risk for type 2 diabetes mellitus (T2DM) in obese, postmenopausal women.