Novel Genetic Variants Associated with Lumbar Spondylosis in Koreans: A Genome-Wide Association Study

OBJECTIVE: The aim of this study was to identify the susceptibility genes responsible for lumbar spondylosis (LS) in Korean patients.METHODS: Data from 1427 subjects were made available for radiographic grading and genome wide association studies (GWAS) analysis. Lateral lumbar spine radiographs were obtained and the various degrees of degenerative change were semi-quantitatively scored. A pilot GWAS was performed using the AffymetrixGenome-Wide Human single-nucleotide polymorphisms (SNPs), 500K array. A total of 352228 SNPs were analyzed and the association between the SNPs and case-control status was analyzed by stepwise logistic regression analyses.RESULTS: The top 100 SNPs with a cutoff p-value of less than 3.7×10⁻⁴ were selected for joint space narrowing, while a cutoff p-value of 6.0×10⁻⁴ was applied to osteophytes and the Kellgren-Lawrence (K-L) osteoarthritis grade. The SNPs with the strongest effect on disc space narrowing, osteophytes, and K-L grade were serine incorporator 1 (rs155467, odds ratio [OR]=17.58, p=1.6×10⁻⁴), stromal interaction molecule 2 (STIM1, rs210781, OR=5.53, p=5×10⁻⁴), and transient receptor potential cation channel, subfamily C (rs11224760, OR=3.99, p=4.8×10⁻⁴), respectively. Leucine-rich repeat-containing G protein-coupled receptor 4 was significantly associated with both disc space narrowing and osteophytes (rs1979400, OR=2.01, p=1.1×10⁻⁴ for disc space narrowing, OR=1.79, p=3×10⁻⁴ for osteophytes), while zinc finger and BTB domain containing 7C was significantly and negatively associated with both osteophytes and a K-L grade >2 (rs12457004,OR=0.25, p=5.8×10⁻⁴ and OR=0.27, p=5.3×10⁻⁴, respectively).CONCLUSION: We identified SNPs that potentially contribute to the pathogenesis of LS. This is the first report of a GWAS in an Asian population.

Novel Genetic Variants Associated with Lumbar Spondylosis in Koreans: A Genome-Wide Association Study

OBJECTIVE: The aim of this study was to identify the susceptibility genes responsible for lumbar spondylosis (LS) in Korean patients. METHODS: Data from 1427 subjects were made available for radiographic grading and genome wide association studies (GWAS) analysis. Lateral lumbar spine radiographs were obtained and the various degrees of degenerative change were semi-quantitatively scored. A pilot GWAS was performed using the AffymetrixGenome-Wide Human single-nucleotide polymorphisms (SNPs), 500K array. A total of 352228 SNPs were analyzed and the association between the SNPs and case-control status was analyzed by stepwise logistic regression analyses. RESULTS: The top 100 SNPs with a cutoff p-value of less than 3.7×10⁻⁴ were selected for joint space narrowing, while a cutoff p-value of 6.0×10⁻⁴ was applied to osteophytes and the Kellgren-Lawrence (K-L) osteoarthritis grade. The SNPs with the strongest effect on disc space narrowing, osteophytes, and K-L grade were serine incorporator 1 (rs155467, odds ratio [OR]=17.58, p=1.6×10⁻⁴), stromal interaction molecule 2 (STIM1, rs210781, OR=5.53, p=5×10⁻⁴), and transient receptor potential cation channel, subfamily C (rs11224760, OR=3.99, p=4.8×10⁻⁴), respectively. Leucine-rich repeat-containing G protein-coupled receptor 4 was significantly associated with both disc space narrowing and osteophytes (rs1979400, OR=2.01, p=1.1×10⁻⁴ for disc space narrowing, OR=1.79, p=3×10⁻⁴ for osteophytes), while zinc finger and BTB domain containing 7C was significantly and negatively associated with both osteophytes and a K-L grade >2 (rs12457004,OR=0.25, p=5.8×10⁻⁴ and OR=0.27, p=5.3×10⁻⁴, respectively). CONCLUSION: We identified SNPs that potentially contribute to the pathogenesis of LS. This is the first report of a GWAS in an Asian population.

Association between MSX1 SNPs and Nonsyndromic Cleft Lip with or without Cleft Palate in the Korean Population

The purpose of this study was to investigate the contribution of MSX1 gene to the risk of nonsyndromic cleft lip with or without cleft palate (NS-CL +/- P) in the Korean population. The samples consisted of 142 NS-CL +/- P families (9 with cleft lip, 26 with cleft lip and alveolus, and 107 with cleft lip and palate; 76 trios and 66 dyads). Three single nucleotide polymorphisms (SNPs: rs3821949, rs12532, and rs4464513) were tested for association with NS-CL +/- P case-parent trios using transmission disequilibrium test (TDT) and conditional logistic regression models (CLRMs). Minor allele frequency, heterozygosity, chi2 test for Hardy-Weinberg equilibrium, and pairwise linkage disequilibrium (LD) at each SNP were computed. The family- and haplotype-based association test programs were used to perform allelic and genotypic TDTs for individual SNPs and to fabricate sliding windows of haplotypes. Genotypic odds ratios (GORs) were obtained from CLRMs using R software. Although the family-based TDT indicated a meaningful association for rs3821949 (P = 0.028), the haplotype analysis did not reveal any significant association with rs3821949, rs12532, or rs4464513. The A allele at rs3821949 had a significant increased risk of NS-CL +/- P (GOR, 1.64; 95% confidence interval,1.03-2.63; P = 0.038, additive model). A positive association is suggested between MSX1 rs3821949 and NS-CL +/- P in the Korean population.

Genetic Risk Prediction for Normal-Karyotype Acute Myeloid Leukemia Using Whole-Exome Sequencing

Normal-karyotype acute myeloid leukemia (NK-AML) is a highly malignant and cytogenetically heterogeneous hematologic cancer. We searched for somatic mutations from 10 pairs of tumor and normal cells by using a highly efficient and reliable analysis workflow for whole-exome sequencing data and performed association tests between the NK-AML and somatic mutations. We identified 21 nonsynonymous single nucleotide variants (SNVs) located in a coding region of 18 genes. Among them, the SNVs of three leukemia-related genes (MUC4, CNTNAP2, and GNAS) reported in previous studies were replicated in this study. We conducted stepwise genetic risk score (GRS) models composed of the NK-AML susceptible variants and evaluated the prediction accuracy of each GRS model by computing the area under the receiver operating characteristic curve (AUC). The GRS model that was composed of five SNVs (rs75156964, rs56213454, rs6604516, rs10888338, and rs2443878) showed 100% prediction accuracy, and the combined effect of the three reported genes was validated in the current study (AUC, 0.98; 95% confidence interval, 0.92 to 1.00). Further study with large sample sizes is warranted to validate the combined effect of these somatic point mutations, and the discovery of novel markers may provide an opportunity to develop novel diagnostic and therapeutic targets for NK-AML.

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.

ADIPOQ Gene Variants Associated with Susceptibility to Obesity and Low Serum Adiponectin Levels in Healthy Koreans / 한국역학회지

OBJECTIVES: This study aimed to measure the association between the adiponectin, C1Q and collagen domain-containing (ADIPOQ) gene variants and obesity in Koreans. METHODS: Three single nucleotide polymorphisms located in the ADIPOQ gene were genotyped in a population-based cross-sectional study of 986 healthy Koreans. Three different case-control groups (i.e. G1, G2, and G3) were defined according to body mass index (BMI) and serum adiponectin levels. Allelic and genotypic associations of this gene with obesity were measured using multivariate logistic regression analyses in each group. RESULTS: The G allele of -11377C>G, a polymorphism located in the promoter region of the ADIPOQ gene (odds ratio (OR), 1.48; 95% confidence interval, 1.13-1.94) and most haplotypes including this allele significantly increased the risk for obesity. However, the OR decreased from 3.98 (G1 group) to 2.90 (G2 group) and 2.30 (G3 group) when a less strict definition of obesity was used. Most haplotypes, including this allele, significantly increased the risk of obesity. The statistical evidence from the GG genotype of -11377C>G (OR, 3.98) and the GT/GT diplotype composed of -11377G>C and +45T>G (OR, 5.20) confirmed the contribution of the G allele toward a predisposition for obesity. CONCLUSION: These results suggest the contribution of the ADIPOQ gene toward susceptibility to obesity in healthy Koreans. The high-risk genotypes and haplotypes identified here may provide more information for identifying individuals who are at risk of obesity.

Genome-Wide Association Analyses on Blood Pressure Using Three Different Phenotype Definitions

Hypertension is the most prevalent disease worldwide and is itself a risk factor for cerebral, cardiac, and renal diseases. The inconsistency of candidate genes suggested by previous genomewide association studies (GWASs) may be due to not only differences in study design and genetic or environmental background but also the difference in the power of analysis between continuous traits and discrete traits. We analyzed 352,228 single nucleotide polymorphisms (SNPs) in 8842 unrelated Koreans obtained from Ansan and Ansung cohorts. We performed a series of GWA analyses using three different phenotype models; young hypertensive cases (278 subjects) versus elderly normotensive controls (680 subjects); the upper 25% (2211 hypertensive cases) versus the lower 25% of the SBP distribution (2211 hypotensive controls); and finally SBP and DBP as continuous traits (8842 subjects). The numbers of young hypertensive cases and elderly normotensive controls were not large enough to achieve genomewide significance. The model comparing the upper 25% subjects to the lower 25% of subjects showed a power that was approximate to that of QTL analysis. Two neighboring SNPs of the ATP2B1 gene, rs17249754 (SBP, p=2.53-10; DBP, p=1.28x10-8) and rs7136259 (SBP, p=1.30x10-9; DBP, p=6.41x10-8), were associated with both SBP and DBP. Interestingly, a SNP of the RPL6 gene, rs11066280, revealed a significant genomewide association with SBP in men only (p=3.85x10-8), and four SNPs located near the MAN2A1 gene showed a strong association with DBP only in elderly men aged 60-70 years (e.g., rs6421827, p=4.86x10-8). However, we did not observe any gene variant attaining genome-wide significance consistently in the three phenotype models except for the ATP2B1 gene variants. In general, the association signal with blood pressure was stronger in women than in men. Genes identified in GWASs are expected to open the way for prevention, early diagnosis, and personalized treatment of hypertension.

Candidate gene polymorphisms for diabetes mellitus, cardiovascular disease and cancer are associated with longevity in Koreans

Long-lived people may have a unique genetic makeup that makes them more resistant than the general population to prevalent age-related diseases; however, not much is known about genes involved in the longevity. To identify susceptibility variants controlling longevity, we performed a high-throughput candidate gene study using 137 Koreans over 90 yr old and 213 young healthy Koreans. We evaluated 463 informative markers located in 176 candidate genes mostly for diabetes mellitus, cardiovascular disease and cancer under five genetic models. We estimated the odds ratios for each allele, genotype, haplotype, and gene-gene interaction using logistic regression analysis. Associations between 13 genes and longevity were detected at a P-value less than 0.01. Particularly, the rs671 (A) allele of the aldehyde dehydrogenase 2 family (mitochondrial) (ALDH2) gene was associated with longevity only in men (OR 2.11, P = 0.008). Four genes, proprotein convertase subtilisin/kexin type 1 (PCSK1, P = 0.008), epidermal growth factor receptor (EGFR, P = 0.003), paired box 4 (PAX4, P = 0.008), and V-yes-1 Yamaguchi sarcoma viral related oncogene homolog (LYN, P = 0.002) consistently yielded statistical evidence for association with longevity. The findings of the current study may provide a starting point for future studies to unravel genetic factors controlling longevity in Koreans.