Single Nucleotide Polymorphism Coverage Research Articles

Assessment of genotype imputation performance using 1000 Genomes in African American studies.

Genotype imputation, used in genome-wide association studies to expand coverage of single nucleotide polymorphisms (SNPs), has performed poorly in African Americans compared to less admixed populations. Overall, imputation has typically relied on HapMap reference haplotype panels from Africans (YRI), European Americans (CEU), and Asians (CHB/JPT). The 1000 Genomes project offers a wider range of reference populations, such as African Americans (ASW), but their imputation performance has had limited evaluation. Using 595 African Americans genotyped on Illumina’s HumanHap550v3 BeadChip, we compared imputation results from four software programs (IMPUTE2, BEAGLE, MaCH, and MaCH-Admix) and three reference panels consisting of different combinations of 1000 Genomes populations (February 2012 release): (1) 3 specifically selected populations (YRI, CEU, and ASW); (2) 8 populations of diverse African (AFR) or European (AFR) descent; and (3) all 14 available populations (ALL). Based on chromosome 22, we calculated three performance metrics: (1) concordance (percentage of masked genotyped SNPs with imputed and true genotype agreement); (2) imputation quality score (IQS; concordance adjusted for chance agreement, which is particularly informative for low minor allele frequency [MAF] SNPs); and (3) average r2hat (estimated correlation between the imputed and true genotypes, for all imputed SNPs). Across the reference panels, IMPUTE2 and MaCH had the highest concordance (91%–93%), but IMPUTE2 had the highest IQS (81%–83%) and average r2hat (0.68 using YRI+ASW+CEU, 0.62 using AFR+EUR, and 0.55 using ALL). Imputation quality for most programs was reduced by the addition of more distantly related reference populations, due entirely to the introduction of low frequency SNPs (MAF≤2%) that are monomorphic in the more closely related panels. While imputation was optimized by using IMPUTE2 with reference to the ALL panel (average r2hat = 0.86 for SNPs with MAF>2%), use of the ALL panel for African American studies requires careful interpretation of the population specificity and imputation quality of low frequency SNPs.

A fine-mapping study of 7 top scoring genes from a GWAS for major depressive disorder.

Major depressive disorder (MDD) is a psychiatric disorder that is characterized -amongst others- by persistent depressed mood, loss of interest and pleasure and psychomotor retardation. Environmental circumstances have proven to influence the aetiology of the disease, but MDD also has an estimated 40% heritability, probably with a polygenic background. In 2009, a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. A non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became only nominally significant after post-hoc analysis with an Australian cohort which used similar ascertainment. The absence of genome-wide significance may be caused by low SNP coverage of genes. To increase SNP coverage to 100% for common variants (m.a.f.>0.1, r2>0.8), we selected seven genes from the GAIN-MDD GWAS: PCLO, GZMK, ANPEP, AFAP1L1, ST3GAL6, FGF14 and PTK2B. We genotyped 349 SNPs and obtained the lowest P-value for rs2715147 in PCLO at P = 6.8E−7. We imputed, filling in missing genotypes, after which rs2715147 and rs2715148 showed the lowest P-value at P = 1.2E−6. When we created a haplotype of these SNPs together with the non-synonymous coding SNP rs2522833, the P-value decreased to P = 9.9E−7 but was not genome wide significant. Although our study did not identify a more strongly associated variant, the results for PCLO suggest that the causal variant is in high LD with rs2715147, rs2715148 and rs2522833.

Hot topic: Performance of bovine high-density genotyping platforms in Holsteins and Jerseys

Two high-density single nucleotide polymorphism (SNP) genotyping arrays have recently become available for bovine genomic analyses, the Illumina High-Density Bovine BeadChip Array (777,962 SNP) and the Affymetrix Axiom Genome-Wide BOS 1 Array (648,874 SNP). These products each have unique design and chemistry attributes, and the extent of marker overlap and their potential utility for quantitative trait loci fine mapping, detection of copy number variation, and multibreed genomic selection are of significant interest to the cattle community. This is the first study to compare the performance of these 2 arrays. Deoxyribonucleic acid samples from 16 dairy cattle (10 Holstein, 6 Jersey) were used for the comparison. An independent set of DNA samples taken from 46 Jersey cattle and 18 Holstein cattle were used to ascertain the amount of SNP variation accounted by the 16 experimental samples. Data were analyzed with SVS7 software (Golden Helix Inc., Bozeman, MT) to remove SNP having a call rate less than 90%, and linkage disequilibrium pruning was used to remove linked SNP (r2 ≥ 0.9). Maximum, average, and median gaps were calculated for each analysis based on genomic position of SNP on the bovine UMD3.1 genome assembly. All samples were successfully genotyped (≥98% SNP genotyped) with both platforms. The average number of genotyped SNP in the Illumina platform was 775,681 and 637,249 for the Affymetrix platform. Based on genomic position, a total of 107,896 SNP were shared between the 2 platforms; however, based on genotype concordance, only 96,031 SNP had complete concordance at these loci. Both Affymetrix BOS 1 and Illumina BovineHD genotyping platforms are well designed and provide high-quality genotypes and similar coverage of informative SNP. Despite fewer total SNP on BOS 1, 19% more SNP remained after linkage disequilibrium pruning, resulting in a smaller gap size (5.2 vs. 6.9kb) in Holstein and Jersey samples relative to BovineHD. However, only 224,115 Illumina and 241,038 Affymetrix SNP remained following removal of SNP with a minor allele frequency of zero in Holstein and Jersey samples, resulting in an average gap size of 11,887 bp and 11,018 bp, respectively. Combining the 354,348 informative (r2 ≥ 0.9), polymorphic (minor allele frequency ≥ 0), unique SNP data from both platforms decreased the average gap size to 7,560 bp. Genome-wide copy number variant analyses were performed using intensity files from both platforms. The BovineHD platform provided an advantage to the copy number variant data compared with the BOS 1 because of the larger number of SNP, higher intensity signals, and lower background effects. The combined use of both platforms significantly improved coverage over either platform alone and decreased the gap size between SNP, providing a valuable tool for fine mapping quantitative trait loci and multibreed animal evaluation.

Meta analysis of candidate gene variants outside the LPA locus with Lp(a) plasma levels in 14,500 participants of six White European cohorts

Background Both genome-wide association studies and candidate gene studies have reported that the major determinant of plasma levels of the Lipoprotein (a) [Lp(a)] reside within the LPA locus on chromosome 6. We have used data from the HumanCVD BeadChip to explore the contribution of other candidate genes determining Lp(a) levels. Methods 48,032 single nucleotide polymorphisms (SNPs) from the Illumina HumanCVD BeadChip were genotyped in 5059 participants of the Whitehall II study (WHII) of randomly ascertained healthy men and women. SNPs showing association with Lp(a) levels of p < 10 −4 outside the LPA locus were selected for replication in a total of an additional 9463 participants of five European based studies (EAS, EPIC-Norfolk, NPHSII, PROCARDIS, and SAPHIR). Results In Whitehall II, apart from the LPA locus (where p values for several SNPs were <10 −30) there was significant association at four loci GALNT2, FABP1, PPARGC1A and TNFRSFF11A. However, a meta-analysis of the six studies did not confirm any of these findings. Conclusion Results from this meta analysis of 14,522 participants revealed no candidate genes from the HumanCVD BeadChip outside the LPA locus to have an effect on Lp(a) levels. Further studies with genome-wide and denser SNP coverage are required to confirm or refute this finding.

Abstract 3811: A germline polymorphism in the ANRIL (CDKN2BAS) locus is associated with susceptibility to Philadelphia-positive acute lymphoblastic leukemia (ALL)

Abstract Background: The CDKN2A/B locus is inactivated in several haeamatologic malignancies mainly due to deletion, aberrant repression or epigenetic silencing, but little is known about the potential role of its single nucleotide polymorphisms (SNPs) in leukemia susceptibility. Patients and methods: To investigate whether polymorphisms within this locus can correlate with increased susceptibility to haeamatologic malignancies, an association study was performed by genotyping 23 SNPs spanning the MTAP, CDKN2A/B and CDKN2BAS loci, as well as relative intergenic regions, in a case-control cohort made up of 332 samples: 149 leukemia patients, including Philadelphia positive (Ph+) ALL (n=92) and acute myeloid leukemia (AML) samples (n=57), and 183 unrelated healthy controls. The median age was 52 years (18-78 years) and 53 years (21-71 years) in Ph+ ALL and AML patients, respectively. AML cases included FAB M0-M5, miscellaneous cytogenetic abnormalities and normal karyotype subtypes. 6 SNPs were selected on the basis of their previous association with several diseases, such as coronary artery disease, type 2 diabetes mellitus, frailty. The remaining 17 SNPs were selected among those included in the Affymetrix Genome-Wide Human SNP Array 6.0 to deepen the SNPs coverage for the examined region. Genotyping was performed using iPLEX Gold technology and MassARRAY high-throughput DNA analysis with Matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (Sequenom, Inc., San Diego, CA). Results: 17 SNPs, spanning the 9p genomic interval that encompasses the MTAP, CDKN2A/B and CDKN2BAS loci, were successfully genotyped and used for investigating their potential associations with the leukemia phenotypes. Potential population stratification affecting the control sample was ruled out as its genotypes distribution satisfies the Hardy-Weinberg equilibrium criterion. The rs564398, mapping to the CDKN2BAS locus that encodes for ANRIL antisense non-coding RNA, showed a statistically significant correlation with the ALL phenotype, with a risk pattern that was compatible with an overdominant model of disease susceptibility and an Odd Ratio (OR) of 2 (95% CI, 1.20 to 3.33; p = 7.1 × 10-3). Since a co-ordinated regulation of ANRIL and p14/ARF, p16/CDKN2A, p15/CDKN2B transcription has been already observed in both physiologic and pathologic conditions, we hypothesized that rs564398 association reflects a condition of high linkage disequilibrium between such polymorphism and a causative variant that is able to alter CDKN2A/B expression profiles by changing ANRIL dosage, thus leading to abnormal proliferative boosts and consequent increased ALL susceptibility. Supported by European LeukemiaNet, AIL, AIRC, FIRB 2006, Ateneo RFO grants, Project of integrated program, Programma di Ricerca Regione – Università 2007 – 2009. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3811. doi:10.1158/1538-7445.AM2011-3811

Rationale for an international consortium to study inherited genetic susceptibility to childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia is the major pediatric cancer in developed countries. To date most association studies of acute lymphoblastic leukemia have been based on the candidate gene approach and have evaluated a restricted number of polymorphisms. Such studies have served to highlight difficulties in conducting statistically and methodologically rigorous investigations into acute lymphoblastic leukemia risk. Recent genome-wide association studies of childhood acute lymphoblastic leukemia have provided robust evidence that common variation at four genetic loci confers a modest increase in risk. The accumulated experience to date and relative lack of success of initial efforts to identify novel acute lymphoblastic leukemia predisposition loci emphasize the need for alternative study designs and methods. The International Childhood Acute Lymphoblastic Leukaemia Genetics Consortium includes 12 research groups in Europe, Asia, the Middle East and the Americas engaged in studying the genetics of acute lymphoblastic leukemia. The initial goal of this consortium is to identify and characterize low-penetrance susceptibility variants for acute lymphoblastic leukemia through association-based analyses. Efforts to develop genome-wide association studies of acute lymphoblastic leukemia, in terms of both sample size and single nucleotide polymorphism coverage, and to increase the number of single nucleotide polymorphisms taken forward to large-scale replication should lead to the identification of additional novel risk variants for acute lymphoblastic leukemia. Ethnic differences in the risk of acute lymphoblastic leukemia are well recognized and thus in assessing the interplay between inherited and non-genetic risk factors, analyses using different population cohorts with different incidence rates are likely to be highly informative. Given that the frequency of many acute lymphoblastic leukemia subgroups is small, identifying differential effects will realistically only be possible through multi-center pooled analyses. Here, we review the rationale for identifying genetic risk variants for acute lymphoblastic leukemia and our proposed strategy for establishing the International Childhood Acute Lymphoblastic Leukaemia Genetics Consortium.

Susceptibility to Philadelphia-Positive acute Lymphoblastic Leukemia (ALL) Is Associated with a Germline Polymorphism In the ANRIL (CDKN2BAS) Locus.

AbstractAbstract 1670▪Among genes that are supposed to be involved in altered molecular pathways leading to leukemogenesis, the CDKN2A/B locus is particularly noteworthy since it encodes for the INK4-class cyclin dependent kinase inhibitors p15 INK4B, p16 INK4A and p14 ARF, which act as tumor suppressors. The CDKN2A/B locus has been found to be inactivated in several haeamatologic malignancies mainly due to deletion, aberrant repression or epigenetic silencing, whereas little is known about the potential role of its single nucleotide polymorphisms (SNPs) in leukemia susceptibility. To investigate whether polymorphisms within this broad genomic region can correlate with increased susceptibility to haeamatologic malignancies, an association study was performed by genotyping 23 SNPs spanning the MTAP, CDKN2A/B and CDKN2BAS loci, as well as relative intergenic regions, in a case-control cohort made up of 332 samples: 149 leukemia patients, including Philadelphia positive (Ph+) ALL (n=92) and acute myeloid leukemia (AML) samples (n=57), and 183 unrelated healthy controls. Ph+ ALL patients showed a median age of 52 years (18-78 years); 24 cases (26%) expressed the p210 oncoprotein, 62 (67%) the p190 and 6 (7%) both the proteins. AML patients showed a median age of 53 years (21-71 years) and they included FAB-M0, M1, M2, M3, M4, M5, miscellaneous cytogenetic abnormalities and normal karyotype subtypes. 6 SNPs were selected on the basis of their previous association with several diseases, such as coronary artery disease (rs2891168, rs518394, rs564398, rs10757278), type 2 diabetes mellitus (rs564398), frailty (rs2811712). The remaining 17 SNPs were selected among those included in the Affymetrix Genome-Wide Human SNP Array 6.0 to deepen the SNPs coverage for the examined region. Genotyping was performed using iPLEX Gold technology and MassARRAY high-throughput DNA analysis with Matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (Sequenom, Inc., San Diego, CA). Furthermore, all leukemia samples were also previously genotyped using the GeneChip® Human Mapping 250K NspI and Genome-Wide Human SNP 6.0 (Affymetrix). A total of 17 SNPs, spanning the 9p genomic interval that encompasses the MTAP, CDKN2A/B and CDKN2BAS loci, were successfully genotyped and used for investigating their potential associations with the leukemia phenotypes. Five SNPs (rs1012713, rs10965179, rs34011899, rs3731232, rs3218010) with Minor Allele Frequency (MAF) <0.05 in cases and controls, as well as one SNP (rs3931609) showing >20% missing call rates, were instead excluded from the association analysis. Potential population stratification affecting the control sample was ruled out as its genotypes distribution satisfies the Hardy-Weinberg equilibrium criterion. Among the 17 SNPs, rs564398, mapping to the CDKN2BAS locus that encodes for ANRIL antisense non-coding RNA, showed a statistically significant correlation with the ALL phenotype, with a risk pattern that was compatible with an overdominant model of disease susceptibility and an Odd Ratio (OR) of 2 (95% CI, 1.20 to 3.33; p= 7.1 × 10-3). Since a co-ordinated regulation of ANRIL and p14/ARF, p16/CDKN2A, p15/CDKN2B transcription has been already observed in both physiologic and pathologic conditions, we hypothesized that rs564398 association reflects a condition of high linkage disequilibrium between such polymorphism and a causative variant that is able to alter CDKN2A/B expression profiles by changing ANRIL dosage, thus leading to abnormal proliferative boosts and consequent increased ALL susceptibility. Supported by European LeukemiaNet, AIL, AIRC, FIRB 2006, PRIN 2008, Ateneo RFO grants, Project of integrated program (PIO), Programma di Ricerca Regione – Università 2007 – 2009. Disclosures:Baccarani:Novartis: Consultancy, Honoraria; BMS: Consultancy, Honoraria. Martinelli:Novartis: Consultancy, Honoraria; BMS: Honoraria; Pfizer: Consultancy.

Poor replication of candidate genes for major depressive disorder using genome-wide association data

Data from the Genetic Association Information Network (GAIN) genome-wide association study (GWAS) in major depressive disorder (MDD) were used to explore previously reported candidate gene and single-nucleotide polymorphism (SNP) associations in MDD. A systematic literature search of candidate genes associated with MDD in case-control studies was performed before the results of the GAIN MDD study became available. Measured and imputed candidate SNPs and genes were tested in the GAIN MDD study encompassing 1738 cases and 1802 controls. Imputation was used to increase the number of SNPs from the GWAS and to improve coverage of SNPs in the candidate genes selected. Tests were carried out for individual SNPs and the entire gene using different statistical approaches, with permutation analysis as the final arbiter. In all, 78 papers reporting on 57 genes were identified, from which 92 SNPs could be mapped. In the GAIN MDD study, two SNPs were associated with MDD: C5orf20 (rs12520799; P=0.038; odds ratio (OR) AT=1.10, 95% CI 0.95-1.29; OR TT=1.21, 95% confidence interval (CI) 1.01-1.47) and NPY (rs16139; P=0.034; OR C allele=0.73, 95% CI 0.55-0.97), constituting a direct replication of previously identified SNPs. At the gene level, TNF (rs76917; OR T=1.35, 95% CI 1.13-1.63; P=0.0034) was identified as the only gene for which the association with MDD remained significant after correction for multiple testing. For SLC6A2 (norepinephrine transporter (NET)) significantly more SNPs (19 out of 100; P=0.039) than expected were associated while accounting for the linkage disequilibrium (LD) structure. Thus, we found support for involvement in MDD for only four genes. However, given the number of candidate SNPs and genes that were tested, even these significant may well be false positives. The poor replication may point to publication bias and false-positive findings in previous candidate gene studies, and may also be related to heterogeneity of the MDD phenotype as well as contextual genetic or environmental factors.

Robust QTL fine mapping by applying a quantitative transmission disequilibrium test to the Mendelian sampling term

In many farm animal populations, high-density single nucleotide polymorphism (SNP) genotypes are becoming available on a large scale, and routine estimation of breeding values is implemented for a multiplicity of traits. We propose to apply the basic principle of the quantitative transmission disequilibrium test (QTDT) to estimated Mendelian sampling terms. A two-step procedure is suggested, where in the first step additive breeding values are estimated with a mixed linear model and the Mendelian sampling terms are calculated from the estimated breeding values. In the second step, the QTDT is applied to these estimated Mendelian sampling terms. The resulting test is expected to yield significant results if the SNP is in sufficient linkage disequilibrium and linkage with quantitative trait loci (QTL). This principle is illustrated with a simulated data set comprising 4665 individuals genotyped for 6000 SNP and 15 true QTL. Thirteen of the fifteen QTL were significant on a genome-wide 0.1% error level. Results for the empirical power are derived from repeated samples of 1000 and 3000 genotyped individuals, respectively. General properties and potential extensions of the methodology are indicated. Owing to its computational simplicity and speed, the suggested procedure is well suited to scan whole genomes with high-density SNP coverage in samples of substantial size and for a multiplicity of different traits.

Design and evaluation of a panel of single-nucleotide polymorphisms in microRNA genomic regions for association studies in human disease

MicroRNAs (miRNA) are recognized posttranscriptional gene repressors involved in the control of almost every biological process. Allelic variants in these regions may be an important source of phenotypic diversity and contribute to disease susceptibility. We analyzed the genomic organization of 325 human miRNAs (release 7.1, miRBase) to construct a panel of 768 single-nucleotide polymorphisms (SNPs) covering approximately 1 Mb of genomic DNA, including 131 isolated miRNAs (40%) and 194 miRNAs arranged in 48 miRNA clusters, as well as their 5-kb flanking regions. Of these miRNAs, 37% were inside known protein-coding genes, which were significantly associated with biological functions regarding neurological, psychological or nutritional disorders. SNP coverage analysis revealed a lower SNP density in miRNAs compared with the average of the genome, with only 24 SNPs located in the 325 miRNAs studied. Further genotyping of 340 unrelated Spanish individuals showed that more than half of the SNPs in miRNAs were either rare or monomorphic, in agreement with the reported selective constraint on human miRNAs. A comparison of the minor allele frequencies between Spanish and HapMap population samples confirmed the applicability of this SNP panel to the study of complex disorders among the Spanish population, and revealed two miRNA regions, hsa-mir-26a-2 in the CTDSP2 gene and hsa-mir-128-1 in the R3HDM1 gene, showing geographical allelic frequency variation among the four HapMap populations, probably because of differences in natural selection. The designed miRNA SNP panel could help to identify still hidden links between miRNAs and human disease.

Association and Interaction Analyses of GABBR1 and GABBR2 with Nicotine Dependence in European- and African-American Populations

Previous studies have demonstrated that the γ-aminobutyric acid type B (GABAB) receptor plays an essential role in modulating neurotransmitter release and regulating the activity of ion channels and adenyl cyclase. However, whether the naturally occurring polymorphisms in the two GABAB receptor subunit genes interact with each other to alter susceptibility to nicotine dependence (ND) remains largely unknown. In this study, we genotyped 5 and 33 single nucleotide polymorphisms (SNPs) for GABAB receptor subunit 1 and 2 genes (GABBR1, GABBR2), respectively, in a sample of 2037 individuals from 602 nuclear families of African- American (AA) or European-American (EA) origin. We conducted association analyses to determine (1) the association of each subunit gene with ND at both the individual SNP and haplotype levels and (2) the collective effect(s) of SNPs in both GABAB subunits on the development of ND. Several individual SNPs and haplotypes in GABBR2 were significantly associated with ND in both ethnic samples. Two haplotypes in AAs and one haplotype in EAs showed a protective effect against ND, whilst two other haplotypes in AAs and three haplotypes in EAs showed a risk effect for developing ND. Interestingly, these significant haplotypes were confined to two regions of GABBR2 in the AA and EA samples. Additionally, we found two minor haplotypes in GABBR1 to be positively associated with Heaviness of Smoking Index (HSI) in the EA sample. Finally, we demonstrated the presence of epistasis between GABBR1 and GABBR2 for developing ND. The variants of GABBR1 and GABBR2 are significantly associated with ND, and the involvement of GABBR1 is most likely through its interaction with GABBR2, whereas GABBR2 polymorphisms directly alter susceptibility to ND. Future studies are needed with more dense SNP coverage of GABBR1 and GABBR2 to verify the epistatic effects of the two subunit genes.

Association analysis between polymorphisms in the conserved dopamine neurotrophic factor (CDNF) gene and cocaine dependence

Cocaine-induced neuroplasticity changes in the mesocorticolimbic dopamine systems are thought to be involved in the pathophysiology of cocaine dependence. Since neurotrophic factors have been observed to prevent/reverse and mimic cocaine-induced neurobiological changes in the brain, related genes are plausible candidates for susceptibility to cocaine dependence. The novel conserved dopamine neurotrophic factor protein (CDNF) promotes the survival, growth, and function of dopamine-specific neurons and is expressed in brain regions that undergo cocaine-induced neuroplasticity. In this study, we hypothesize that polymorphisms in the CDNF gene (CDNF/ARMETL1) contribute to increased risk for cocaine dependence. Cocaine dependent individuals (n=351) and unaffected controls (n=257) of African descent were genotyped for four single nucleotide polymorphisms (SNPs) in the CDNF gene (rs11259365, rs7094179, rs7900873, rs2278871). We observed no significant differences in allele, genotype, or haplotype frequencies between cases and controls for any of the tested SNPs. Our study suggests that there is no association between variants in the CDNF gene and cocaine dependence. However, additional studies using larger sample sizes, comprehensive SNP coverage, and clinically homogenous populations are necessary before confidently excluding CDNF as a significant genetic risk factor for cocaine dependence.

Maternal and offspring genetic variants of AKR1C3 and the risk of childhood leukemia

The aldo-keto reductase 1C3 (AKR1C3) gene located on chromosome 10p15-p14, a regulator of myeloid cell proliferation and differentiation, represents an important candidate gene for studying human carcinogenesis. In a prospectively enrolled population-based case-control study of Han Chinese conducted in Kaohsiung in southern Taiwan, a total of 114 leukemia cases and 221 controls <20 years old were recruited between November 1997 and December 2005. The present study set out to evaluate the association between childhood leukemia and both maternal and offspring's genotypes. To do so, we conducted a systematic assessment of common single-nucleotide polymorphisms (SNPs) at the 5' flanking 10 kb to 3' UTR of AKR1C3 gene. Gln5His and three tagSNPs (rs2245191, rs10508293 and rs3209896) and one multimarker (rs2245191, rs10508293 and rs3209896) were selected with average 90% coverage of untagged SNPs by using the HapMap II data set. Odds ratios and 95% confidence intervals were adjusted for age and gender. After correcting for multiple comparisons, we observed that risk of developing childhood leukemia is significantly associated with rs10508293 polymorphism on intron 4 of the AKR1C3 gene in both offspring alone and in the combined maternal and offspring genotypes (nominal P < 0.0001, permutation P < 0.005). The maternal methylenetetrahydrofolate reductase A1298C polymorphism was found to be an effect modifier of the maternal intron 4 polymorphism of the AKR1C3 gene (rs10508293) and the childhood leukemia risk. In conclusion, this study suggests that AKR1C3 polymorphisms may be important predictive markers for childhood leukemia susceptibility.

Polymorphisms in Interleukin-15 Gene on Chromosome 4q31.2 Are Associated with Psoriasis Vulgaris in Chinese Population

Through a series of linkage analyses in a large Chinese family cohort of psoriasis, we previously identified and confirmed a non-HLA psoriasis linkage locus PSORS9 within a small region at 4q31.2-32.1. Within the critical region of the PSORS9 locus, IL-15 has been long recognized as a strong candidate gene for psoriasis. In this study, we investigated the association between IL-15 genetic polymorphisms and psoriasis in a large Chinese sample. Highly significant evidence for association was identified at a single-nucleotide polymorphism (SNP) (g.96516A --> T) within the 3'-untranslated region (UTR) of the IL-15 gene (P=0.00006, after correction for multiple testing). Haplotype analysis using the SNPs within the 3'UTR region also provided strong supporting evidence for association (P=0.00005), where we identified a haplotype of the 3'UTR region of IL-15 associated with increased risk to psoriasis (odds ratio=1.65). This association was also supported by the results of our expression activity analyses, where we demonstrated that the identified risk haplotype is associated with an increased activity of IL-15. Therefore, we provided early evidence for the important role of IL-15 genetic variants in the pathogenesis of psoriasis, probably by increasing interleukin production and inflammation in the lesions of psoriasis.

Copy Number Variants and Common Disorders: Filling the Gaps and Exploring Complexity in Genome-Wide Association Studies

Genome-wide association scans (GWASs) using single nucleotide polymorphisms (SNPs) have been completed successfully for several common disorders and have detected over 30 new associations. Considering the large sample sizes and genome-wide SNP coverage of the scans, one might have expected many of the common variants underpinning the genetic component of various disorders to have been identified by now. However, these studies have not evaluated the contribution of other forms of genetic variation, such as structural variation, mainly in the form of copy number variants (CNVs). Known CNVs account for over 15% of the assembled human genome sequence. Since CNVs are not easily tagged by SNPs, might have a wide range of copy number variability, and often fall in genomic regions not well covered by whole-genome arrays or not genotyped by the HapMap project, current GWASs have largely missed the contribution of CNVs to complex disorders. In fact, some CNVs have already been reported to show association with several complex disorders using candidate gene/region approaches, underpinning the importance of regions not investigated in current GWASs. This reveals the need for new generation arrays (some already in the market) and the use of tailored approaches to explore the full dimension of genome variability beyond the single nucleotide scale.

Single-Nucleotide Polymorphisms and Haplotypes of Bone Morphogenetic Protein Genes and Peripheral Bone Mineral Density in Young Korean Men and Women

Bone morphogenetic proteins (BMPs) play critical roles in osteoblast differentiation. To investigate the association between common single-nucleotide polymorphisms (SNPs) of BMPs and bone mineral density (BMD), a cross-sectional study was conducted in healthy Korean men (n = 237) and women (n = 276) aged 20-39 years. Calcaneus and distal radius BMD were measured by dual energy X-ray absorptiometry. SNPs of BMP2 (-1103C > A, c.584G > A, IVS1-2744A > G, c.893T > A), BMP4 (c.712T > C, IVS1-160C > T), and BMP6 (c.1283C > G, IVS4-6838A > G, IVS5 + 24C > T) were determined using the 5'-nuclease assay. Significant associations were observed between BMP2 c.584G > A, c.893T > A genotypes and male calcaneus as well as female distal radius BMD. Men with the BMP2 c.893 AA genotype had a 16% higher BMD at the calcaneus (P for trend = 0.014), whereas women with this genotype had a 7% lower BMD at the distal radius than the other genotypes (P for trend = 0.010). A significant association was also observed between BMP4 IVS1-160C > T and male calcaneus BMD (P for trend = 0.024). When the association between haplotypes and BMD was investigated, the AAGA haplotype of BMP2 was significantly associated with low bone mass in female distal radius (P for trend = 0.013). These results suggested that one or more SNPs of BMP2 and BMP4 are associated with peripheral BMD in Korean men and women. However, this association is dependent on anatomical sites and gender. Thus, larger studies with complete coverage of SNPs are needed in the future.

The usefulness of different density SNP maps for disease association studies of common variants.

Large-scale discovery and validation of single-nucleotide polymorphisms (SNPs) facilitates indirect association mapping. It has recently been estimated that, in Europeans, 77% of all SNPs with frequency of 10% or more could be ascertained through linkage disequilibrium (LD) by genotyping variants in the database dbSNP. Using a sampling approach from 73 genes with near complete SNP maps, we show here the usefulness of SNP maps at different densities and the large variability of SNP coverage in different genomic regions. While even sparse SNP maps are of some value to genetic mapping, in order to undertake disease association studies providing at least 80% of SNPs in 90% of genes, much denser maps need to be constructed, at more than one SNP per kb in some regions.