• Genomics & Informatics
• /
• v.17 no.4
• /
• pp.48.1-48.6
• /
• 2019

Over the last decade, genome-wide association studies (GWASs) have provided an unprecedented amount of genetic variations that are associated with various phenotypes. However, previous GWAS were mostly conducted in European populations, and these biased results for non-Europeans may result in a significant reduction in risk prediction for non-Europeans. An issue with the early GWAS was the winner's curse problem, which led to misleading results when constructing the polygenic risk scores (PRS). Therefore, more non-European population-based studies are needed to validate reported variants and improve genetic risk assessment across diverse populations. In this study, we validated 422 variants independently associated with glycemic indexes, liver enzymes, and type 2 diabetes in 125,872 samples from a Korean population, and further validated the results by assessing publicly available summary statistics from European GWAS (n = 898,130). Among the 422 independently associated variants, 284, 320, and 361 variants were replicated in Koreans, Europeans, and either one of the two populations. In addition, the effect sizes for Koreans and Europeans were moderately correlated (r = 0.33-0.68). However, 61 variants were not replicated in both Koreans and Europeans. Our findings provide valuable information on effect sizes and statistical significance, which is essential to improve the assessment of disease risk using PRS analysis.

• Asian-Australasian Journal of Animal Sciences
• /
• v.32 no.5
• /
• pp.607-613
• /
• 2019

Objective: Intramuscular fat (IMF) content plays an important role in meat quality. Identification of single nucleotide polymorphisms (SNPs) and genes related to pig IMF, especially using pig populations with high IMF content variation, can help to establish novel molecular breeding tools for optimizing IMF in pork and unveil the mechanisms that underlie fat metabolism. Methods: We collected muscle samples of 453 Chinese Lulai black pigs, measured IMF content by Soxhlet petroleum-ether extraction method, and genotyped genome-wide SNPs using GeneSeek Genomic Profiler Porcine HD BeadChip. Then a genome-wide association study was performed using a linear mixed model implemented in the GEMMA software. Results: A total of 43 SNPs were identified to be significantly associated with IMF content by the cutoff p<0.001. Among these significant SNPs, the greatest number of SNPs (n = 19) were detected on Chr.9, and two linkage disequilibrium blocks were formed among them. Additionally, 17 significant SNPs are mapped to previously reported quantitative trait loci (QTLs) of IMF and confirmed previous QTLs studies. Forty-two annotated genes centering these significant SNPs were obtained from Ensembl database. Overrepresentation test of pathways and gene ontology (GO) terms revealed some enriched reactome pathways and GO terms, which mainly involved regulation of basic material transport, energy metabolic process and signaling pathway. Conclusion: These findings improve our understanding of the genetic architecture of IMF content in pork and facilitate the follow-up study of fine-mapping genes that influence fat deposition in muscle.

• Genomics & Informatics
• /
• v.15 no.4
• /
• pp.178-182
• /
• 2017

Next-generation sequencing (NGS) technology has become a trend in the genomics research area. There are many software programs and automated pipelines to analyze NGS data, which can ease the pain for traditional scientists who are not familiar with computer programming. However, downstream analyses, such as finding differentially expressed genes or visualizing linkage disequilibrium maps and genome-wide association study (GWAS) data, still remain a challenge. Here, we introduce a dockerized web application written in R using the Shiny platform to visualize pre-analyzed RNA sequencing and GWAS data. In addition, we have integrated a genome browser based on the JBrowse platform and an automated intermediate parsing process required for custom track construction, so that users can easily build and navigate their personal genome tracks with in-house datasets. This application will help scientists perform series of downstream analyses and obtain a more integrative understanding about various types of genomic data by interactively visualizing them with customizable options.

• Journal of Genetic Medicine
• /
• v.14 no.1
• /
• pp.8-17
• /
• 2017

Purpose: Pulse wave velocity (PWV) is an indicator of arterial stiffness, and is considered a marker of vascular damage. However, a genome-wide association study analyzing single nucleotide polymorphisms (SNPs) associated with brachial-ankle PWV (baPWV) has not been conducted in healthy populations. We performed this study to identify SNPs associated with baPWV in healthy populations in Korea. Materials and Methods: Genomic SNPs data for 2,407 individuals from three sites were analyzed as part of the Korean Genomic Epidemiologic Study. Without replication samples, we performed multivariable analysis as a post hoc analysis to verify the findings in site adjusted analysis. Healthy subjects aged between 40 and 70 years without self-reported history or diagnosis of hypertension, diabetes, hyperlipidemia, heart disease, cerebrovascular disease and cancer were included. We excluded subjects with a creatinine level >1.4 mg/dL (men) and 1.2 mg/dL (women). Results: In the site-adjusted association analysis, significant associations (P<$5{\times}10^{-8}$) with baPWV were detected for only 5 SNPs with low minor allele frequency. In multivariable analysis adjusted by age, sex, height, body mass index, mean arterial pressure, site, smoking, alcohol, and exercise, 11 SNPs were found to be associated (P<$5{\times}10^{-8}$) with baPWV. The 5 SNPs (P<$5{\times}10^{-8}$) linked to three genes (OPCML, PRR35 and RAB40C) were common between site-adjusted analysis and multivariable analysis. However, meta-analysis of the result from three sites for the 11 SNPs showed no significant associations. Conclusion: Using the recent standard for genome-wide association study, we did not find any evidence of significant association signals with baPWV.

• Genomics & Informatics
• /
• v.20 no.1
• /
• pp.14.1-14.4
• /
• 2022

Permutation testing is a robust and popular approach for significance testing in genomic research that has the advantage of reducing inflated type 1 error rates; however, its computational cost is notorious in genome-wide association studies (GWAS). Here, we developed a supercomputing-aided approach to accelerate the permutation testing for GWAS, based on the message-passing interface (MPI) on parallel computing architecture. Our application, called MPI-GWAS, conducts MPI-based permutation testing using a parallel computing approach with our supercomputing system, Nurion (8,305 compute nodes, and 563,740 central processing units [CPUs]). For 10⁷ permutations of one locus in MPI-GWAS, it was calculated in 600 s using 2,720 CPU cores. For 10⁷ permutations of ~30,000-50,000 loci in over 7,000 subjects, the total elapsed time was ~4 days in the Nurion supercomputer. Thus, MPI-GWAS enables us to feasibly compute the permutation-based GWAS within a reason-able time by harnessing the power of parallel computing resources.

• Genomics & Informatics
• /
• v.14 no.4
• /
• pp.149-159
• /
• 2016

With the success of the genome-wide association studies (GWASs), many candidate loci for complex human diseases have been reported in the GWAS catalog. Recently, many disease prediction models based on penalized regression or statistical learning methods were proposed using candidate causal variants from significant single-nucleotide polymorphisms of GWASs. However, there have been only a few systematic studies comparing existing methods. In this study, we first constructed risk prediction models, such as stepwise linear regression (SLR), least absolute shrinkage and selection operator (LASSO), and Elastic-Net (EN), using a GWAS chip and GWAS catalog. We then compared the prediction accuracy by calculating the mean square error (MSE) value on data from the Korea Association Resource (KARE) with body mass index. Our results show that SLR provides a smaller MSE value than the other methods, while the numbers of selected variables in each model were similar.

• Journal of Animal Science and Technology
• /
• v.62 no.3
• /
• pp.306-312
• /
• 2020

Korean army dogs are raised for special purposes and have contributed much to society. However, several diseases occur in dogs. Canine hip dysplasia (CHD) is a musculoskeletal disorder that occurs frequently in Korean army dogs and interferes with their activities. If we could control CHD, this would have a positive effect on their performance. This study performed a genome-wide association study (GWAS) in 69 Korean army dogs to find significant loci for CHD using 170K single nucleotide polymorphisms (SNPs). CHD was classified according to the Norberg angle criterion. The control group comprised 62 dogs classified as relatively normal, and 7 dogs with severe CHD formed the case group. From the GWAS analysis, we concluded that SNPs present on chromosome 4 might have a significant impact on the overall expression of canine hip dysplasia.

• Proceedings of the Korean Society of Crop Science Conference
• /
• 2017.06a
• /
• pp.58-58
• /
• 2017

The advent of next generation sequencing technology has elicited plenty of sequencing data available in agriculturally relevant plant species. For most crop species, it is too expensive to obtain the whole genome sequence data with sufficient coverage. Thus, many approaches have been developed to bring down the cost of NGS. Genotyping-by-sequencing (GBS) is a cost-effective genotyping method for complex genetic populations. GBS can be used for the analysis of genomic selection (GS), genome-wide association study (GWAS) and constructing haplotype and genetic linkage maps in a variety of plant species. For efficiently dealing with plant GBS data, the TASSEL-GBS pipeline is one of the most popular choices for many researchers. TASSEL-GBS is JAVA based a software package to obtain genotyping data from raw GBS sequences. Here, we describe application of GBS and bioinformatics pipeline of TASSEL-GBS for analyzing plant genetics data.

• Asian-Australasian Journal of Animal Sciences
• /
• v.26 no.11
• /
• pp.1523-1528
• /
• 2013

To increase plumage color uniformity and understand the genetic background of Korean chickens, we performed a genome-wide association study of different plumage color in Korean native chickens. We analyzed 60K SNP chips on 279 chickens with GEMMA methods for GWAS and estimated the genetic heritability for plumage color. The estimated heritability suggests that plumage coloration is a polygenic trait. We found new loci associated with feather pigmentation at the genome-wide level and from the results infer that there are additional genetic effect for plumage color. The results will be used for selecting and breeding chicken for plumage color uniformity.

• Journal of Korean Neurosurgical Society
• /
• v.61 no.1
• /
• pp.66-74
• /
• 2018

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{\times}10^{-4}$ were selected for joint space narrowing, while a cutoff p-value of $6.0{\times}10^{-4}$ 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{\times}10^{-4}$), stromal interaction molecule 2 (STIM1, rs210781, OR=5.53, $p=5{\times}10^{-4}$), and transient receptor potential cation channel, subfamily C (rs11224760, OR=3.99, $p=4.8{\times}10^{-4}$), 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{\times}10^{-4}$ for disc space narrowing, OR=1.79, $p=3{\times}10^{-4}$ 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{\times}10^{-4}$ and OR=0.27, $p=5.3{\times}10^{-4}$, 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.