Functional analysis of HapMap SNPs

Genome-wide association studies (GWAS) have successfully identified many genetic variants associated with complex diseases and traits. However, functional consequence of genetic variants studied in GWAS is not yet fully investigated, which would hinder the application of GWAS. We therefore performed a systematic functional analysis of HapMap SNPs, which have been most commonly used as the reference panel for GWAS. Our study highlights several characteristics of HapMap SNPs and identifies subsets of genetic variants with interesting functional implication. The results show that HapMap SNPs have good coverage within RefSeq genes, especially within known disease-related genes. On the other hand, only a small percentage of SNPs are non-synonymous SNPs while many SNPs are actually located at gene deserts. Moreover, many functionally important variants are not yet still interrogated. A redesigned SNP reference panel with additional functionally important variants would be useful to identify disease-causal variants in the future genome-wide studies.     

Effective selection of informative SNPs and classification on the HapMap genotype data

Background  

Since the single nucleotide polymorphisms (SNPs) are genetic variations which determine the difference between any two unrelated individuals, the SNPs can be used to identify the correct source population of an individual. For efficient population identification with the HapMap genotype data, as few informative SNPs as possible are required from the original 4 million SNPs. Recently, Park et al. (2006) adopted the nearest shrunken centroid method to classify the three populations, i.e., Utah residents with ancestry from Northern and Western Europe (CEU), Yoruba in Ibadan, Nigeria in West Africa (YRI), and Han Chinese in Beijing together with Japanese in Tokyo (CHB+JPT), from which 100,736 SNPs were obtained and the top 82 SNPs could completely classify the three populations.     

Genomic Scans of Zygotic Disequilibrium and Epistatic SNPs in HapMap Phase III Populations

Previous theory indicates that zygotic linkage disequilibrium (LD) is more informative than gametic or composite digenic LD in revealing natural population history. Further, the difference between the composite digenic and maximum zygotic LDs can be used to detect epistatic selection for fitness. Here we corroborate the theory by investigating genome-wide zygotic LDs in HapMap phase III human populations. Results show that non-Africa populations have much more significant zygotic LDs than do Africa populations. Africa populations (ASW, LWK, MKK, and YRI) possess more significant zygotic LDs for the double-homozygotes (D_AABB) than any other significant zygotic LDs (D_AABb, D_AaBB, and D_AaBb), while non-Africa populations generally have more significant D_AaBb’s than any other significant zygotic LDs (D_AABB, D_AABb, and D_AaBB). Average r-squares for any significant zygotic LDs increase generally in an order of populations YRI, MKK, CEU, CHB, LWK, JPT, CHD, TSI, GIH, ASW, and MEX. Average r-squares are greater for D_AABB and D_AaBb than for D_AaBB and D_AABb in each population. YRI and MKK can be separated from LWK and ASW in terms of the pattern of average r-squares. All population divergences in zygotic LDs can be interpreted with the model of Out of Africa for modern human origins. We have also detected 19735-95921 SNP pairs exhibiting strong signals of epistatic selection in different populations. Gene-gene interactions for some epistatic SNP pairs are evident from empirical findings, but many more epistatic SNP pairs await evidence. Common epistatic SNP pairs rarely exist among all populations, but exist in distinct regions (Africa, Europe, and East Asia), which helps to understand geographical genomic medicine.     

Comparing genetic variants detected in the 1000 genomes project with SNPs determined by the International HapMap Consortium

Single-nucleotide polymorphisms (SNPs) determined based on SNP arrays from the international HapMap consortium (HapMap) and the genetic variants detected in the 1000 genomes project (1KGP) can serve as two references for genomewide association studies (GWAS). We conducted comparative analyses to provide a means for assessing concerns regarding SNP array-based GWAS findings as well as for realistically bounding expectations for next generation sequencing (NGS)-based GWAS. We calculated and compared base composition, transitions to transversions ratio, minor allele frequency and heterozygous rate for SNPs from HapMap and 1KGP for the 622 common individuals. We analysed the genotype discordance between HapMap and 1KGP to assess consistency in the SNPs from the two references. In 1KGP, 90.58% of 36,817,799 SNPs detected were not measured in HapMap. More SNPs with minor allele frequencies less than 0.01 were found in 1KGP than HapMap. The two references have low discordance (generally smaller than 0.02) in genotypes of common SNPs, with most discordance from heterozygous SNPs. Our study demonstrated that SNP array-based GWAS findings were reliable and useful, although only a small portion of genetic variances were explained. NGS can detect not only common but also rare variants, supporting the expectation that NGS-based GWAS will be able to incorporate a much larger portion of genetic variance than SNP arrays-based GWAS.     

GLIDERS - A web-based search engine for genome-wide linkage disequilibrium between HapMap SNPs

Background  

A number of tools for the examination of linkage disequilibrium (LD) patterns between nearby alleles exist, but none are available for quickly and easily investigating LD at longer ranges (>500 kb). We have developed a web-based query tool (GLIDERS: Genome-wide LInkage DisEquilibrium Repository and Search engine) that enables the retrieval of pairwise associations with r² ≥ 0.3 across the human genome for any SNP genotyped within HapMap phase 2 and 3, regardless of distance between the markers.     

FstSNP-HapMap3: a database of SNPs with high population differentiation for HapMap3

The International HapMap Project has recently made available genotypes and frequency data for phase 3 (NCBI build 36, dbSNPb129) of the HapMap providing an enriched genotype dataset for approximately 1.6 million single nucleotide polymorphisms (SNPs) from 1,115 individuals with ancestry from parts of Africa, Asia, Europe, North America and Mexico. In the present study, we aim to facilitate pharmacogenetics studies by providing a database of SNPs with high population differentiation through a genomewide test on allele frequency variation among 11 HapMap3 samples. Common SNPs with minor allele frequency greater than 5¢ from each of 11 HapMap3 samples were included in the present analysis. The population differentiation is measured in terms of fixation index (Fst), and the SNPs with Fst values over 0.5 were defined as highly differentiated SNPs. Our tests were carried out between all pairs of the 11 HapMap3 samples or among subgroups with the same continental ancestries. Altogether we carried out 64 genomewide Fst tests and identified 28,215 highly differentiated SNPs for 49 different combinations of HapMap3 samples in the current database.     

SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap

SUMMARY: The interpretation of genome-wide association results is confounded by linkage disequilibrium between nearby alleles. We have developed a flexible bioinformatics query tool for single-nucleotide polymorphisms (SNPs) to identify and to annotate nearby SNPs in linkage disequilibrium (proxies) based on HapMap. By offering functionality to generate graphical plots for these data, the SNAP server will facilitate interpretation and comparison of genome-wide association study results, and the design of fine-mapping experiments (by delineating genomic regions harboring associated variants and their proxies). AVAILABILITY: SNAP server is available at http://www.broad.mit.edu/mpg/snap/.     

An evaluation of the performance of tag SNPs derived from HapMap in a Caucasian population

The Haplotype Map (HapMap) project recently generated genotype data for more than 1 million single-nucleotide polymorphisms (SNPs) in four population samples. The main application of the data is in the selection of tag single-nucleotide polymorphisms (tSNPs) to use in association studies. The usefulness of this selection process needs to be verified in populations outside those used for the HapMap project. In addition, it is not known how well the data represent the general population, as only 90–120 chromosomes were used for each population and since the genotyped SNPs were selected so as to have high frequencies. In this study, we analyzed more than 1,000 individuals from Estonia. The population of this northern European country has been influenced by many different waves of migrations from Europe and Russia. We genotyped 1,536 randomly selected SNPs from two 500-kbp ENCODE regions on Chromosome 2. We observed that the tSNPs selected from the CEPH (Centre d'Etude du Polymorphisme Humain) from Utah (CEU) HapMap samples (derived from US residents with northern and western European ancestry) captured most of the variation in the Estonia sample. (Between 90% and 95% of the SNPs with a minor allele frequency of more than 5% have an r² of at least 0.8 with one of the CEU tSNPs.) Using the reverse approach, tags selected from the Estonia sample could almost equally well describe the CEU sample. Finally, we observed that the sample size, the allelic frequency, and the SNP density in the dataset used to select the tags each have important effects on the tagging performance. Overall, our study supports the use of HapMap data in other Caucasian populations, but the SNP density and the bias towards high-frequency SNPs have to be taken into account when designing association studies.     

Constructing an initial map of transmission distortion based on high density HapMap SNPs across the human autosomes

Transmission distortion （TD） is a significant departure from Mendelian predictions of genes or chromosomes to offspring. While many biological processes have been implicated, there is still much to be understood about TD in humans. Here we present our findings from a genome-wide scan for evidence of TD using haplotype data of 60 trio families from the International HapMap Project. Fisher＇s exact test was applied to assess the extent of TD in 629,958 SNPs across the autosomes. Based on the empirical distribution of PFisher and further permutation tests, we identified 1,205 outlier loci and 224 candidate genes with TD. Using the PANTHER gene ontology database, we found 19 categories of biological processes with an enrichment of candidate genes. In particular, the “protein phosphorylation” category contained the largest number of candidates in both HapMap samples. Further analysis uncovered an intriguing non-synonymous change in PPPIR12B, a gene related to protein phosphorylation, which appears to influence the allele transmission from male parents in the YRI （Yoruba from Ibadan, Nigeria） population. Our findings also indicate an ethnicity-related property of TD signatures in HapMap samples and provide new clues for our understanding of TD in humans.     

猪CAPN1基因部分外显子及3'UTR区的SNPs检测

以野猪.民猪和大白猪为研究对象,根据网上公布的序列设计了7对引物,采用测序,PCR-SSCP和PCR-RFLP方法对CAPN1基因的部分外显子和3'UTR区进行了单核苷酸多态性检测和基因型分析,探讨CAPN1基因多态性与瘦肉率和嫩度的关系.研究发现11个SNPs,其中5个位于外显子,4个位于内含子,2个位于3'UTR区,外显子中的突变有一处是错义突变,导致了蛋白质多肽链第260位氨基酸发生了M/V的替代.群体遗传学分析表明,在所检测的各多态位点上,野猪、民猪、大白猪3个品种间不同基因型的分布都存在着极显著的差异(P<0.01),而野猪和民猪之间各基因型的分布差异不显著(P>0.05),民猪和大白猪之间各基因型的分布存在着极显著的差异(P<0.01).结合品种特性分析表明,P4、P6引物和3'UTR区Hinf1位点所检测的不同基因型和瘦肉率具有一定的相关性.     

猪CAPN1基因部分外显子及3′UTR区的SNPs检测

以野猪、民猪和大白猪为研究对象, 根据网上公布的序列设计了7对引物, 采用测序、PCR-SSCP和PCR-RFLP方法对CAPN1基因的部分外显子和3′UTR区进行了单核苷酸多态性检测和基因型分析, 探讨CAPN1基因多态性与瘦肉率和嫩度的关系。研究发现11个SNPs, 其中5个位于外显子, 4个位于内含子, 2个位于3′UTR区, 外显子中的突变有一处是错义突变, 导致了蛋白质多肽链第260位氨基酸发生了M/V的替代。群体遗传学分析表明, 在所检测的各多态位点上, 野猪、民猪、大白猪3个品种间不同基因型的分布都存在着极显著的差异(P<0.01), 而野猪和民猪之间各基因型的分布差异不显著(P>0.05), 民猪和大白猪之间各基因型的分布存在着极显著的差异(P<0.01)。结合品种特性分析表明, P4、P6引物和3′ UTR区HinfⅠ位点所检测的不同基因型和瘦肉率具有一定的相关性。     

Genome‐wide association of drought‐related and biomass traits with HapMap SNPs in Medicago truncatula

Improving drought tolerance of crop plants is a major goal of plant breeders. In this study, we characterized biomass and drought‐related traits of 220 Medicago truncatula HapMap accessions. Characterized traits included shoot biomass, maximum leaf size, specific leaf weight, stomatal density, trichome density and shoot carbon‐13 isotope discrimination (δ¹³C) of well‐watered M. truncatula plants, and leaf performance in vitro under dehydration stress. Genome‐wide association analyses were carried out using the general linear model (GLM), the standard mixed linear model (MLM) and compressed MLM (CMLM) in TASSEL, which revealed significant overestimation of P‐values by CMLM. For each trait, candidate genes and chromosome regions containing SNP markers were found that are in significant association with the trait. For plant biomass, a 0.5 Mbp region on chromosome 2 harbouring a plasma membrane intrinsic protein, PIP2, was discovered that could potentially be targeted to increase dry matter yield. A protein disulfide isomerase‐like protein was found to be tightly associated with both shoot biomass and leaf size. A glutamate‐cysteine ligase and an aldehyde dehydrogenase family protein with Arabidopsis homologs strongly expressed in the guard cells were two of the top genes identified by stomata density genome‐wide association studies analysis.     

Evaluating HapMap SNP data transferability in a large-scale genotyping project involving 175 cancer-associated genes

One of the many potential uses of the HapMap project is its application to the investigation of complex disease aetiology among a wide range of populations. This study aims to assess the transferability of HapMap SNP data to the Spanish population in the context of cancer research. We have carried out a genotyping study in Spanish subjects involving 175 candidate cancer genes using an indirect gene-based approach and compared results with those for HapMap CEU subjects. Allele frequencies were very consistent between the two samples, with a high positive correlation (R) of 0.91 (P<<1×10⁻⁶). Linkage disequilibrium patterns and block structures across each gene were also very similar, with disequilibrium coefficient (r ²) highly correlated (R=0.95, P<<1×10⁻⁶). We found that of the 21 genes that contained at least one block larger than 60 kb, nine (ATM, ATR, BRCA1, ERCC6, FANCC, RAD17, RAD50, RAD54B and XRCC4) belonged to the GO category “DNA repair”. Haplotype frequencies per gene were also highly correlated (mean R=0.93), as was haplotype diversity (R=0.91, P<<1×10⁻⁶). “Yin yang” haplotypes were observed for 43% of the genes analysed and 18% of those were identical to the ancestral haplotype (identified in Chimpazee). Finally, the portability of tagSNPs identified in the HapMap CEU data using pairwise r ² thresholds of 0.8 and 0.5 was assessed by applying these to the Spanish and current HapMap data for 66 genes. In general, the HapMap tagSNPs performed very well. Our results show generally high concordance with HapMap data in allele frequencies and haplotype distributions and confirm the applicability of HapMap SNP data to the study of complex diseases among the Spanish population. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Evaluating the transferability of Hapmap SNPs to a Singapore Chinese population

Background

Sexual reproduction relies on two key events: formation of cells with a haploid genome (the gametes) and restoration of diploidy after fertilization. Therefore the underlying mechanisms must have been evolutionary linked and there is a need for evidence that could support such a model.

Results

We describe the identification and the characterization of yem ¹ , the first yem-alpha mutant allele (V478E), which to some extent affects diploidy reduction and its restoration. Yem-alpha is a member of the Ubinuclein/HPC2 family of proteins that have recently been implicated in playing roles in chromatin remodeling in concert with HIRA histone chaperone. The yem ¹ mutant females exhibited disrupted chromosome behavior in the first meiotic division and produced very low numbers of viable progeny. Unexpectedly these progeny did not display paternal chromosome markers, suggesting that they developed from diploid gametes that underwent gynogenesis, a form of parthenogenesis that requires fertilization.

Conclusions

We focus here on the analysis of the meiotic defects exhibited by yem ¹ oocytes that could account for the formation of diploid gametes. Our results suggest that yem ¹ affects chromosome segregation presumably by affecting kinetochores function in the first meiotic division. This work paves the way to further investigations on the evolution of the mechanisms that support sexual reproduction.     

Evaluating GWAS-Identified SNPs for Age at Natural Menopause among Chinese Women

Background

Age at natural menopause (ANM) is a complex trait with high heritability and is associated with several major hormonal-related diseases. Recently, several genome-wide association studies (GWAS), conducted exclusively among women of European ancestry, have discovered dozens of genetic loci influencing ANM. No study has been conducted to evaluate whether these findings can be generalized to Chinese women.

Methodology/Principal Findings

We evaluated the index single nucleotide polymorphisms (SNPs) in 19 GWAS-identified genetic susceptibility loci for ANM among 3,533 Chinese women who had natural menopause. We also investigated 3 additional SNPs which were in LD with the index SNP in European-ancestry but not in Asian-ancestry populations. Two genetic risk scores (GRS) were calculated to summarize SNPs across multiple loci one for all SNPs tested (GRSall), and one for SNPs which showed association in our study (GRSsel). All 22 SNPs showed the same association direction as previously reported. Eight SNPs were nominally statistically significant with P≤0.05: rs4246511 (RHBDL2), rs12461110 (NLRP11), rs2307449 (POLG), rs12611091 (BRSK1), rs1172822 (BRSK1), rs365132 (UIMC1), rs2720044 (ASH2L), and rs7246479 (TMEM150B). Especially, SNPs rs4246511, rs365132, rs1172822, and rs7246479 remained significant even after Bonferroni correction. Significant associations were observed for GRS. Women in the highest quartile began menopause 0.7 years (P = 3.24×10⁻⁹) and 0.9 years (P = 4.61×10⁻¹¹) later than those in the lowest quartile for GRSsel and GRSall, respectively.

Conclusions

Among the 22 investigated SNPs, eight showed associations with ANM (P<0.05) in our Chinese population. Results from this study extend some recent GWAS findings to the Asian-ancestry population and may guide future efforts to identify genetic determination of menopause.     

On the challenges of the HapMap resource

The International HapMap Project provides a key resource of genotypic data on human lymphoblastoid cell lines derived from four major world populations of European, African, Chinese and Japanese ancestry for researchers to associate with various phenotypic data to find genes affecting health, disease and response to drugs. Recently, the HapMap resource has significantly benefited research areas such as gene expression variation studies. Besides some intrinsic limitations, there are a few challenges that should be considered in the next wave of research using this tremendous resource. We suggest that overcoming these challenges or considering the confounding variables in the interpretation of results can provide more insights into the current views of the human genome as well as complex traits such as drug response variation and susceptibility to common diseases.     

Rank and Order: Evaluating the Performance of SNPs for Individual Assignment in a Non-Model Organism

Single nucleotide polymorphisms (SNPs) are valuable tools for ecological and evolutionary studies. In non-model species, the use of SNPs has been limited by the number of markers available. However, new technologies and decreasing technology costs have facilitated the discovery of a constantly increasing number of SNPs. With hundreds or thousands of SNPs potentially available, there is interest in comparing and developing methods for evaluating SNPs to create panels of high-throughput assays that are customized for performance, research questions, and resources. Here we use five different methods to rank 43 new SNPs and 71 previously published SNPs for sockeye salmon: F_ST, informativeness (I_n), average contribution to principal components (LC), and the locus-ranking programs BELS and WHICHLOCI. We then tested the performance of these different ranking methods by creating 48- and 96-SNP panels of the top-ranked loci for each method and used empirical and simulated data to obtain the probability of assigning individuals to the correct population using each panel. All 96-SNP panels performed similarly and better than the 48-SNP panels except for the 96-SNP BELS panel. Among the 48-SNP panels, panels created from F_ST, I_n, and LC ranks performed better than panels formed using the top-ranked loci from the programs BELS and WHICHLOCI. The application of ranking methods to optimize panel performance will become more important as more high-throughput assays become available.     

HGDP and HapMap Analysis by Ancestry Mapper Reveals Local and Global Population Relationships

Knowledge of human origins, migrations, and expansions is greatly enhanced by the availability of large datasets of genetic information from different populations and by the development of bioinformatic tools used to analyze the data. We present Ancestry Mapper, which we believe improves on existing methods, for the assignment of genetic ancestry to an individual and to study the relationships between local and global populations. The principle function of the method, named Ancestry Mapper, is to give each individual analyzed a genetic identifier, made up of just 51 genetic coordinates, that corresponds to its relationship to the HGDP reference population. As a consequence, the Ancestry Mapper Id (AMid) has intrinsic biological meaning and provides a tool to measure similarity between world populations. We applied Ancestry Mapper to a dataset comprised of the HGDP and HapMap data. The results show distinctions at the continental level, while simultaneously giving details at the population level. We clustered AMids of HGDP/HapMap and observe a recapitulation of human migrations: for a small number of clusters, individuals are grouped according to continental origins; for a larger number of clusters, regional and population distinctions are evident. Calculating distances between AMids allows us to infer ancestry. The number of coordinates is expandable, increasing the power of Ancestry Mapper. An R package called Ancestry Mapper is available to apply this method to any high density genomic data set.