Practical Calling Approach for Exome Array-Based Genome-Wide Association Studies in Korean Population

A genome wide association study for lung function in the Korean population using an exome array

BACKGROUND/AIMS

Lung function is an objective indicator of diagnosis and prognosis of respiratory diseases. Many common genetic variants have been associated with lung function in multiple ethnic populations. We looked for coding variants associated with forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC) in the Korean general population.

METHODS

We carried out exome array analysis and lung function measurements of the FEV1 and FEV1/FVC in 7,524 individuals of the Korean population. We evaluated single variants with minor allele frequency greater than 0.5%. We performed look-ups for candidate coding variants associations in the UK Biobank, SpiroMeta, and CHARGE consortia.

RESULTS

We identified coding variants in the SMIM29 (C6orf1) (p = 1.2 × 10-5) and HMGA1 locus on chromosome 6p21, the GIT2 (p = 6.5 × 10-5) locus on chromosome 12q24, and the ARHGEF40 (p = 9.9 × 10-5) locus on chromosome 14q11 as having a significant association with lung function (FEV1). We also confirmed a previously reported association with lung function and chronic obstructive pulmonary disease in the FAM13A (p = 4.54 × 10-6) locus on chromosome 4q22, in TNXB (p = 1.30 × 10-6) and in AGER (p = 1.09 × 10-8) locus on chromosome 6p21.

CONCLUSION

Our exome array analysis identified that several protein coding variants were associated with lung function in the Korean population. Common coding variants in SMIM29 (C6orf1), HMGA1, GIT2, FAM13A, TNXB, AGER and low-frequency variant in ARHGEF40 potentially affect lung function, which warrant further study.

An exome-wide rare variant analysis of Korean men identifies three novel genes predisposing to prostate cancer

Since prostate cancer is highly heritable, common variants associated with prostate cancer have been studied in various populations, including those in Korea. However, rare and low-frequency variants have a significant influence on the heritability of the disease. The contributions of rare variants to prostate cancer susceptibility have not yet been systematically evaluated in a Korean population. In this work, we present a large-scale exome-wide rare variant analysis of 7,258 individuals (985 cases with prostate cancer and 6,273 controls). In total, 19 rare variant loci spanning 7 genes contributed to an association with prostate cancer susceptibility. In addition to replicating previously known susceptibility genes (e.g., CDYL2, MST1R, GPER1, and PARD3B), 3 novel genes were identified (FDR q < 0.05), including the non-coding RNAs ENTPD3-AS1, LOC102724438, and protein-coding gene SPATA3. Additionally, 6 pathways were identified based on identified variants and genes, including estrogen signaling pathway, signaling by MST1, IL-15 production, MSP-RON signaling pathway, and IL-12 signaling and production in macrophages, which are known to be associated with prostate cancer. In summary, we report novel genes and rare variants that potentially play a role in prostate cancer susceptibility in the Korean population. These observations demonstrated a path towards one of the fundamental goals of precision medicine, which is to identify biomarkers for a subset of the population with a greater risk of disease than others.

Exome chip-driven association study of lipidemia in >14,000 Koreans and evaluation of genetic effect on identified variants between different ethnic groups

Lipid levels in blood are widely used to diagnose and monitor chronic diseases. It is essential to identify the genetic traits involved in lipid metabolism for understanding chronic diseases. However, the influence of genetic traits varies depending on race, sex, age, and ethnicity. Therefore, research focusing on populations of individual countries is required, and the results can be used as a basis for comparison of results of other studies at the cross-racial and cross-country levels. In the present study, we selected lipid-related variants and evaluated their effects on lipid-related diseases in more than 14,000 subjects of three cohorts using the Illumina Human Exome Beadchip. A genome-wide association study was conducted using EPACTs after adjusting for age, sex, and recruitment area. A genome-wide significance cutoff was defined as p < 5E-08 in all the three cohorts. Sixteen variants represented the lipid traits and were classified as vulnerable to borderline hypertriglyceridemia, hyper-LDL-cholesterolemia, or hypo-HDL-cholesterolemia. Moreover, we compared the genetic effects of the 16 variants between ethnic groups and identified the missense variants in apolipoprotein A-V, cholesterol ester transfer protein, and apolipoprotein E as Asian-specific. Our study provides candidate genes as markers for chronic diseases through the evaluation of genetic effects.

The Korea Biobank Array: Design and Identification of Coding Variants Associated with Blood Biochemical Traits

We introduce the design and implementation of a new array, the Korea Biobank Array (referred to as KoreanChip), optimized for the Korean population and demonstrate findings from GWAS of blood biochemical traits. KoreanChip comprised >833,000 markers including >247,000 rare-frequency or functional variants estimated from >2,500 sequencing data in Koreans. Of the 833 K markers, 208 K functional markers were directly genotyped. Particularly, >89 K markers were presented in East Asians. KoreanChip achieved higher imputation performance owing to the excellent genomic coverage of 95.38% for common and 73.65% for low-frequency variants. From GWAS (Genome-wide association study) using 6,949 individuals, 28 associations were successfully recapitulated. Moreover, 9 missense variants were newly identified, of which we identified new associations between a common population-specific missense variant, rs671 (p.Glu457Lys) of ALDH2, and two traits including aspartate aminotransferase (P = 5.20 × 10⁻¹³) and alanine aminotransferase (P = 4.98 × 10⁻⁸). Furthermore, two novel missense variants of GPT with rare frequency in East Asians but extreme rarity in other populations were associated with alanine aminotransferase (rs200088103; p.Arg133Trp, P = 2.02 × 10⁻⁹ and rs748547625; p.Arg143Cys, P = 1.41 × 10⁻⁶). These variants were successfully replicated in 6,000 individuals (P = 5.30 × 10⁻⁸ and P = 1.24 × 10⁻⁶). GWAS results suggest the promising utility of KoreanChip with a substantial number of damaging variants to identify new population-specific disease-associated rare/functional variants.

The SNP-set based association study identifies ITGA1 as a susceptibility gene of attention-deficit/hyperactivity disorder in Han Chinese

Genome-wide association studies, which detect the association between single-nucleotide polymorphisms (SNPs) and disease susceptibility, have been extensively applied to study attention-deficit/hyperactivity disorder (ADHD), but genome-wide significant associations have not been found yet. Genetic heterogeneity and insufficient genomic coverage may account for the missing heritability. We performed a two-stage association study for ADHD in the Han Chinese population. In the discovery stage, 1033 ADHD patients and 950 healthy controls were genotyped using both the Affymetrix Genome-Wide Human SNP Array 6.0 and the Illumina Infinium HumanExome BeadChip. The genotyped SNPs were combined to generate a powerful SNP set with better genomic coverage especially for the nonsynonymous variants. In addition to the association of single SNPs, we collected adjacent SNPs as SNP sets, which were determined by either genes or successive sliding windows, to evaluate their synergetic effect. The candidate susceptibility SNPs were further replicated in an independent cohort of 1441 ADHD patients and 1447 healthy controls. No genome-wide significant SNPs or gene-based SNP sets were found to be associated with ADHD. However, two continuous sliding windows located in ITGA1 (P-value=8.33E-7 and P-value=8.43E-7) were genome-wide significant. The quantitative trait analyses also demonstrated their association with ADHD core symptoms and executive functions. The association was further validated by follow-up replications for four selected SNPs: rs1979398 (P-value=2.64E-6), rs16880453 (P-value=3.58E-4), rs1531545 (P-value=7.62E-4) and rs4074793 (P-value=2.03E-4). Our results suggest that genetic variants in ITGA1 may be involved in the etiology of ADHD and the SNP-set based analysis is a promising strategy for the detection of underlying genetic risk factors.