Ji Wan Park

A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits

To identify genetic factors influencing quantitative traits of biomedical importance, we conducted a genome-wide association study in 8,842 samples from population-based cohorts recruited in Korea. For height and body mass index, most variants detected overlapped those reported in European samples. For the other traits examined, replication of promising GWAS signals in 7,861 independent Korean samples identified six previously unknown loci. For pulse rate, signals reaching genome-wide significance mapped to chromosomes 1q32 (rs12731740, P = 2.9 × 10−9) and 6q22 (rs12110693, P = 1.6 × 10−9), with the latter ∼400 kb from the coding sequence of GJA1. For systolic blood pressure, the most compelling association involved chromosome 12q21 and variants near the ATP2B1 gene (rs17249754, P = 1.3 × 10−7). For waist-hip ratio, variants on chromosome 12q24 (rs2074356, P = 7.8 × 10−12) showed convincing associations, although no regional transcript has strong biological candidacy. Finally, we identified two loci influencing bone mineral density at multiple sites. On chromosome 7q31, rs7776725 (within the FAM3C gene) was associated with bone density at the radius (P = 1.0 × 10−11), tibia (P = 1.6 × 10−6) and heel (P = 1.9 × 10−10). On chromosome 7p14, rs1721400 (mapping close to SFRP4, a frizzled protein gene) showed consistent associations at the same three sites (P = 2.2 × 10−3, P = 1.4 × 10−7 and P = 6.0 × 10−4, respectively). This large-scale GWA analysis of well-characterized Korean population-based samples highlights previously unknown biological pathways.

Sample Size and Statistical Power Calculation in Genetic Association Studies

A sample size with sufficient statistical power is critical to the success of genetic association studies to detect causal genes of human complex diseases. Genome-wide association studies require much larger sample sizes to achieve an adequate statistical power. We estimated the statistical power with increasing numbers of markers analyzed and compared the sample sizes that were required in case-control studies and case-parent studies. We computed the effective sample size and statistical power using Genetic Power Calculator. An analysis using a larger number of markers requires a larger sample size. Testing a single-nucleotide polymorphism (SNP) marker requires 248 cases, while testing 500,000 SNPs and 1 million markers requires 1,206 cases and 1,255 cases, respectively, under the assumption of an odds ratio of 2, 5% disease prevalence, 5% minor allele frequency, complete linkage disequilibrium (LD), 1:1 case/control ratio, and a 5% error rate in an allelic test. Under a dominant model, a smaller sample size is required to achieve 80% power than other genetic models. We found that a much lower sample size was required with a strong effect size, common SNP, and increased LD. In addition, studying a common disease in a case-control study of a 1:4 case-control ratio is one way to achieve higher statistical power. We also found that case-parent studies require more samples than case-control studies. Although we have not covered all plausible cases in study design, the estimates of sample size and statistical power computed under various assumptions in this study may be useful to determine the sample size in designing a population-based genetic association study.

Adiponectin concentrations: a genome-wide association study.

Adiponectin is associated with obesity and insulin resistance. To date, there has been no genome-wide association study (GWAS) of adiponectin levels in Asians. Here we present a GWAS of a cohort of Korean volunteers. A total of 4,001 subjects were genotyped by using a genome-wide marker panel in a two-stage design (979 subjects initially and 3,022 in a second stage). Another 2,304 subjects were used for follow-up replication studies with selected markers. In the discovery phase, the top SNP associated with mean log adiponectin was rs3865188 in CDH13 on chromosome 16 (p = 1.69 × 10(-15) in the initial sample, p = 6.58 × 10(-39) in the second genome-wide sample, and p = 2.12 × 10(-32) in the replication sample). The meta-analysis p value for rs3865188 in all 6,305 individuals was 2.82 × 10(-83). The association of rs3865188 with high-molecular-weight adiponectin (p = 7.36 × 10(-58)) was even stronger in the third sample. A reporter assay that evaluated the effects of a CDH13 promoter SNP in complete linkage disequilibrium with rs3865188 revealed that the major allele increased expression 2.2-fold. This study clearly shows that genetic variants in CDH13 influence adiponectin levels in Korean adults.

Association between IRF6 and nonsyndromic cleft lip with or without cleft palate in four populations

Purpose: The interferon regulatory factor 6 (IRF6), the gene that causes van der Woude syndrome has been shown to be associated with nonsyndromic cleft lip with or without palate in several populations. This study aimed to confirm the contribution of IRF6 to cleft lip with or without palate risk in additional Asian populations.Methods: A set of 13 single nucleotide polymorphisms was tested for association with cleft lip with or without palate in 77 European American, 146 Taiwanese, 34 Singaporean, and 40 Korean case-parent trios using both the transmission disequilibrium test and conditional logistic regression models.Results: Evidence of linkage and association was observed among all four populations; and two specific haplotypes [GC composed of rs2235373-rs2235371 (p.V274I) and AAG of rs599021-rs2235373-rs595918] showed the most significant over- and undertransmission among Taiwanese cases (P = 9 × 10−6 and P = 5 × 10−6, respectively). The AGC/CGC diplotype composed of rs599021-rs2235373-rs2013162 showed almost a 7-fold increase in risk among the Taiwanese sample (P < 10−3). These results confirmed the contribution of this gene to susceptibility of oral clefts across different populations; however, the specific single nucleotide polymorphisms showing statistical significance differed among ethnic groups.Conclusion: The high-risk genotypes and diplotypes identified here may provide a better understanding of the etiological role of this gene in oral clefts and potential options for genetic counseling.

High throughput SNP and expression analyses of candidate genes for non-syndromic oral clefts

Background: Recent work suggests that multiple genes and several environmental risk factors influence risk for non-syndromic oral clefts, one of the most common birth defects in humans. Advances in high-throughput genotyping technology now make it possible to test multiple markers in many candidate genes simultaneously. Methods: We present findings from family based association tests of single nucleotide polymorphism (SNP) markers in 64 candidate genes genotyped using the BeadArray approach in 58 case-parent trios from Maryland (USA) to illustrate how multiple markers in multiple genes can be analysed. To assess whether these genes were expressed in human craniofacial structures relevant to palate and lip development, we also analysed data from the Craniofacial and Oral Gene Expression Network (COGENE) consortium, and searched public databases for expression profiles of these genes. Results: Thirteen candidate genes showed significant evidence of linkage in the presence of disequilibrium, and ten of these were found to be expressed in relevant embryonic tissues: SP100, MLPH, HDAC4, LEF1, C6orf105, CD44, ALX4, ZNF202, CRHR1, and MAPT. Three other genes showing statistical evidence (ADH1C, SCN3B, and IMP5) were not expressed in the embryonic tissues examined here. Conclusions: This approach demonstrates how statistical evidence on large numbers of SNP markers typed in case-parent trios can be combined with expression data to identify candidate genes for complex disorders. Many of the genes reported here have not been previously studied as candidates for oral clefts and warrant further investigation.

Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations

Isolated oral clefts, including cleft lip with/without cleft palate (CL/P) and cleft palate (CP), have a complex and heterogeneous etiology. Case-parent trios from three populations were used to study genes spanning chromosome 2, where single nucleotide polymorphic (SNP) markers were analyzed individually and as haplotypes. Case-parent trios from three populations (74 from Maryland, 64 from Singapore and 95 from Taiwan) were genotyped for 962 SNPs in 104 genes on chromosome 2, including two well-recognized candidate genes: TGFA and SATB2. Individual SNPs and haplotypes (in sliding windows of 2–5 SNPs) were used to test for linkage and disequilibrium separately in CL/P and CP trios. A novel candidate gene (ZNF533) showed consistent evidence of linkage and disequilibrium in all three populations for both CL/P and CP. SNPs in key regions of ZNF533 showed considerable variability in estimated genotypic odds ratios and their significance, suggesting allelic heterogeneity. Haplotype frequencies for regions of ZNF533 were estimated and used to partition genetic variance into among-and within-population components. Wright’s fixation index, a measure of genetic diversity, showed little difference between Singapore and Taiwan compared with Maryland. The tensin-1 gene (TNS1) also showed evidence of linkage and disequilibrium among both CL/P and CP trios in all three populations, albeit at a lower level of significance. Additional genes (VAX2, GLI2, ZHFX1B on 2p; WNT6–WNT10A and COL4A3–COL4A4 on 2q) showed consistent evidence of linkage and disequilibrium only among CL/P trios in all three populations, and TGFA showed significant evidence in two of three populations.

Tooth loss, hypertension and risk for stroke in a Korean population

Tooth loss has been suggested as a potential risk factor for stroke. We conducted a prospective cohort study of stroke in Korea on hypertension, diabetes, smoking, and tooth loss to characterize their independent effects and interactions. The overall risk of stroke and the risk of different subtypes of stroke were evaluated in relation to tooth loss using Cox proportional hazards models among 867,256 Korean men and women, aged 30-95 years, who received health insurance from the National Health Insurance Corporation and were medically evaluated between 1992 and 1995, with tooth loss measured. The overall prevalence of having at least one tooth removed among the people in the study was 29.8% (31.9% for men and 22.3% for women). During a 14-year follow-up, 28,258 strokes with 5105 fatal strokes occurred. For men and women, tooth loss was associated with total stroke and stroke subtypes. In a multivariable model adjusting for selected covariates, a graded association between higher tooth loss and higher risk of total stroke was observed in men [> or =7 lost teeth versus 0 (hazard ratio (HR)=1.3; 95% confidence interval (CI), 1.2-1.4)] and in women (HR=1.2; 95% CI, 1.0-1.3). The HRs for ischemic and hemorrhagic stroke were also similar in men and women. There was evidence of interaction of hemorrhagic stroke risk with hypertension and tooth loss. Tooth loss is independently associated with increased risk of stroke and hypertension does interact antagonistically, particularly for hemorrhagic stroke.

Meta-Analysis of Homogeneous Subgroups Reveals Association between PDE4D Gene Variants and Ischemic Stroke

Background: An Icelandic study showed a significant positive association between phosphodiesterase 4D (PDE4D) gene variants and stroke. However, subsequent studies reported conflicting results, possibly due to small sample sizes and the heterogeneity of the studies. Method: We performed a meta-analysis on 6 SNPs of the PDE4D gene to investigate the association between this gene and ischemic stroke by integrating the results of previous studies, comprising 11,834 cases and 15,233 controls. A pooled genotypic odds ratio (OR) for each SNP was determined under 3 genetic models (i.e. dominant, recessive, and codominant) using both fixed- and random-effects models with consideration for heterogeneity and publication bias across studies. Results: Among the SNPs included in this study, SNP56 (rs702553) showed the most significant association with ischemic stroke in a meta-analysis comprised of 7 homogenous studies. The overall OR of the TT genotype compared to the AA genotype was 1.29 (95% CI 1.03–1.61; p = 0.022). For SNP83 (rs966221), a protective effect of the ancestral allele T was observed only in Asian populations (ORTT 0.79, 95% CI 0.69–0.90; p = 0.0005). This meta-analysis revealed a significant association of PDE4D gene variants with the risk of ischemic stroke, and further investigations are warranted to evaluate possible ethnic-specific effects.

Association between genes on chromosome 4p16 and non-syndromic oral clefts in four populations

Isolated cleft lip with or without cleft palate and cleft palate are among the most common human birth defects. Several candidate gene studies on MSX1 have shown significant association between markers in MSX1 and risk of oral clefts, and re-sequencing studies have identified multiple mutations in MSX1 in a small minority of cases, which may account for 1–2% of all isolated oral clefts cases. We explored the 2-Mb region around MSX1, using a marker map of 393 single nucleotide polymorphisms (SNPs) in 297 cleft lip, with or without cleft palate, case–parent trios and 84 cleft palate trios from Maryland, Taiwan, Singapore, and Korea. Both individual markers and haplotypes of two to five SNPs showed several regions yielding statistical evidence for linkage and disequilibrium. Two genes (STK32B and EVC) yielded consistent evidence from cleft lip, with or without cleft palate, trios in all four populations. These two genes plus EVC2 also yielded suggestive evidence for linkage and disequilibrium among cleft palate trios. This analysis suggests that several genes, not just MSX1, in this region may influence risk of oral clefts.

Genetic risk assessment for cardiovascular disease with seven genes associated with plasma C-reactive protein concentrations in Asian populations.

Plasma C-reactive protein (CRP) level is a predictor of cardiovascular risk. We performed a meta-analysis on the effect of 12 single-nucleotide polymorphisms (SNPs) within 8 candidate loci in 36 752 Asians. In addition, we created weighted genetic risk scores (wGRSs) to evaluate the combined effects of genetic variants, which were suggested in the meta-analysis, for predicting the risks of elevated CRP levels as well as increased risks of hypertension and cardiovascular disease (CVD) in 748 Koreans. Nine SNPs located in seven genes, CRP, IL6R, GCKR, IL6, CYP17A1, HNF1A and APOE, were significantly associated with circulating CRP levels in this meta-analysis. Two SNPs, rs7310409 (HNF1A, P=3.4 × 10−23) and rs7553007 (CRP, P=3.4 × 10−17), had the most significant effects on CRP levels; and two SNPs, rs2097677 (IL6) and rs1004467 (CYP17A1) have never been found in the previous European meta-analysis. In Koreans, the subjects in the highest wGRS group had an ∼2.5-fold higher mean CRP level compared with those in the lowest wGRS group (P=2.1 × 10−5). We observed significant increases in the risks of hypertension (odds ratio=2.18, P=0.006) and CVD (odds ratio=9.59, P=3.2 × 10−6) among the subjects in the highest wGRS group. The wGRS models specific to Koreans may warrant further validation to be used as a proxy for the risk of CVD in Asians.