Rare variants discovery by extensive whole-genome sequencing of the Han Chinese population in Taiwan: Applications to cardiovascular medicine.

Introduction: A population-specific genomic reference is important for research and clinical practice, yet it remains unavailable for Han Chinese (HC) in Taiwan. Objectives: We report the first whole genome sequencing (WGS) database of HC (1000 Taiwanese genome (1KTW-WGS)) and demonstrate several applications to cardiovascular medicine. Methods: Whole genomes of 997 HC were sequenced to at least 30X depth. A total of 20,117 relatively healthy HC individuals were genotyped using a customized Axiom GWAS array. We performed a genome-wide genotype imputation technique using IMPUTE2. Results: We identified 26.7 million single-nucleotide variants (SNVs) and 4.2 million insertions-deletions. Of the SNVs, 16.1% were novel relative to dbSNP (build 152), and 34.2% were novel relative to gnomAD. A total of 18,450 healthy HC individuals were genotyped using a customized Genome-Wide Association Study (GWAS) array. We identified hypertension-associated variants and developed a hypertension prediction model based on the correlation between the WGS data and GWAS data (combined clinical and genetic models, AUC 0.887), and also identified 3 novel hyperlipidemia-associated variants. Each individual carried an average of 16.42 (SD = 3.72) disease-causing variants. Additionally, we established an online SCN5A (an important cardiac gene) database that can be used to explore racial differences. Finally, pharmacogenetics studies identified HC population-specific SNVs in genes (CYP2C9 and VKORC1) involved in drug metabolism and blood clotting. Conclusion: This research demonstrates the benefits of constructing a population-specific genomic reference database for precision medicine.

Genetic profiles of 103,106 individuals in the Taiwan Biobank provide insights into the health and history of Han Chinese.

Personalized medical care focuses on prediction of disease risk and response to medications. To build the risk models, access to both large-scale genomic resources and human genetic studies is required. The Taiwan Biobank (TWB) has generated high-coverage, whole-genome sequencing data from 1492 individuals and genome-wide SNP data from 103,106 individuals of Han Chinese ancestry using custom SNP arrays. Principal components analysis of the genotyping data showed that the full range of Han Chinese genetic variation was found in the cohort. The arrays also include thousands of known functional variants, allowing for simultaneous ascertainment of Mendelian disease-causing mutations and variants that affect drug metabolism. We found that 21.2% of the population are mutation carriers of autosomal recessive diseases, 3.1% have mutations in cancer-predisposing genes, and 87.3% carry variants that affect drug response. We highlight how TWB data provide insight into both population history and disease burden, while showing how widespread genetic testing can be used to improve clinical care.

Integrative Genome-Wide Association Studies of eQTL and GWAS Data for Gout Disease Susceptibility.

There is a paucity of genome-wide association study on Han Chinese gout patients. We performed a genome-wide association meta-analysis on two Taiwanese cohorts consisting of 758 gout cases and 14166 controls of Han Chinese ancestry. All the participants were recruited from the Taiwan Biobank. For pathway analysis, we applied ICSNPathway (Identify candidate Causal SNPs and Pathways) analysis, and to investigate whether expression-associated genetic variants contribute to gout susceptibility, we systematically integrated lymphoblastoid expression quantitative trait loci (eQTL) and genome-wide association data of gout using Sherlock, a Bayesian statistical frame-work. In the meta-analysis, we found 4 SNPs that reached genome-wide statistical significance (P < 5.0 × 10-8). These SNPs are in or close to ABCG2, PKD2 and NUDT9 gene on chromosome 4. ICSNPathway analysis identified rs2231142 as the candidate causal SNP, and ABCG2 as the candidate gene. Sherlcok analysis identified three genes, which were significantly associated with the risk of gout (PKD2, NUTD9, and NAP1L5). To conclude, we reported novel susceptible loci for gout that has not been previously addressed in the literature.

Population structure of Han Chinese in the modern Taiwanese population based on 10,000 participants in the Taiwan Biobank project.

The Taiwan Biobank (TWB) aims to build a nationwide research database that integrates genomic/epigenomic profiles, lifestyle patterns, dietary habits, environmental exposure history and long-term health outcomes of 300,000 residents of Taiwan. We describe here an investigation of the population structure of Han Chinese on this Pacific island using genotype data of 591,048 SNPs in an initial freeze of 10,801 unrelated TWB participants. In addition to the North-South cline reported in other Han Chinese populations, we find the Taiwanese Han Chinese clustered into three cline groups: 5% were of northern Han Chinese ancestry, 79.9% were of southern Han Chinese ancestry, and 14.5% belonged to a third (T) group. We also find that this T group is genetically distinct from neighbouring Southeast Asians and Austronesian tribes but similar to other southern Han Chinese. Interestingly, high degree of LD between HLA haplotype A*33:03-B*58:01, an MHC allele being of pathological relevance, and SNPs across the MHC region was observed in subjects with T origin, but not in other Han Chinese. This suggested the T group individuals may have experienced evolutionary events independent from the other southern Han Chinese. Based on the newly-discovered population structure, we detect different loci susceptible to type II diabetes in individuals with southern and northern Han Chinese ancestries. Finally, as one of the largest dataset currently available for the Chinese population, genome-wide statistics for the 10,810 subjects are made publicly accessible through Taiwan View (https://taiwanview.twbiobank.org.tw/index; date last accessed October 14, 2016) to encourage future genetic research and collaborations with the island Taiwan.

Identification of susceptibility gene associated with female primary Sjögren's syndrome in Han Chinese by genome-wide association study.

Primary Sjögren's syndrome (PSS) is an autoimmune disease targeting exocrine glands. It ten times more dominantly affects women than men with an onset peak at menopause. The genetic factor predisposing women to PSS remains unclear. Therefore, we aimed to identify susceptibility loci for PSS in women. We performed genome-wide association study (GWAS) using 242 female PSS patients and 1444 female control in Han Chinese population residing in Taiwan. Replication was conducted in an independent cohort of 178 female PSS and 14,432 control subjects. We identified rs117026326 on GTF2I with GWAS significance (P = 1.10 × 10-15) and rs13079920 on RBMS3 with suggestive significance (P = 2.90 × 10-5) associating with PSS in women. The association of RBMS3 was further evidenced by imputation in which rs13072846 (P = 4.89 × 10-5) was identified and confirmed as female PSS associating SNP within the same LD with rs13079920. PSS pathogenesis involves both immune (effector) and exocrine (target) system. We suggested that while GTF2I is a previously reported associating gene which may function in immune system, RBMS3 is a novel susceptibility gene that predisposes women to PSS potentially through modulating acinar apoptosis and TGF-ß signaling in target exocrine system.

Identification of Susceptible Loci and Enriched Pathways for Bipolar II Disorder Using Genome-Wide Association Studies.

BACKGROUND: This study aimed to identify susceptible loci and enriched pathways for bipolar disorder subtype II. METHODS: We conducted a genome-wide association scan in discovery samples with 189 bipolar disorder subtype II patients and 1773 controls, and replication samples with 283 bipolar disorder subtype II patients and 500 controls in a Taiwanese Han population using Affymetrix Axiom Genome-Wide CHB1 Array. We performed single-marker and gene-based association analyses, as well as calculated polygeneic risk scores for bipolar disorder subtype II. Pathway enrichment analyses were employed to reveal significant biological pathways. RESULTS: Seven markers were found to be associated with bipolar disorder subtype II in meta-analysis combining both discovery and replication samples (P<5.0×10-6), including markers in or close to MYO16, HSP90AB3P, noncoding gene LOC100507632, and markers in chromosomes 4 and 10. A novel locus, ETF1, was associated with bipolar disorder subtype II (P<6.0×10-3) in gene-based association tests. Results of risk evaluation demonstrated that higher genetic risk scores were able to distinguish bipolar disorder subtype II patients from healthy controls in both discovery (P=3.9×10-4~1.0×10-3) and replication samples (2.8×10-4~1.7×10-3). Genetic variance explained by chip markers for bipolar disorder subtype II was substantial in the discovery (55.1%) and replication (60.5%) samples. Moreover, pathways related to neurodevelopmental function, signal transduction, neuronal system, and cell adhesion molecules were significantly associated with bipolar disorder subtype II. CONCLUSION: We reported novel susceptible loci for pure bipolar subtype II disorder that is less addressed in the literature. Future studies are needed to confirm the roles of these loci for bipolar disorder subtype II.

Novel susceptibility gene for nonfamilial hypokalemic periodic paralysis.

OBJECTIVE: To identify susceptibility genes to nonfamilial hypokalemic periodic paralysis (hypoKPP) consisting of thyrotoxic periodic paralysis (TPP) and sporadic periodic paralysis (SPP) and explore the potential pathogenic mechanisms. METHODS: We enrolled patients with nonfamilial hypoKPP not carrying mutations in CACNA1S, SCN4A, KCNJ18, or KCNJ2 and conducted genome-wide association analyses comparing 77 patients with TPP and 32 patients with SPP with 1,730 controls in a Han Chinese population in Taiwan. Replication was performed using an independent Han Chinese cohort of 50 patients with TPP, 22 patients with SPP, and 376 controls. RESULTS: We identified 4 single nucleotide polymorphisms (rs312692, rs312736, rs992072, rs393743) located about 100 Kb downstream of KCNJ2 on chromosome 17q24.3 associated with both TPP and SPP reaching genome-wide significance (p < 9 × 10(-8)). rs312736 was mapped to CTD-2378E21.1, a lincRNA, and direct sequencing revealed an exon variant rs312732 (risk allele A) highly associated with both TPP (p = 1.81 × 10(-12); odds ratio [OR] 3.22 [95% confidence interval (CI) 2.36-4.40]) and SPP (p = 8.6 × 10(-12); OR 5.4 [95% CI 3.17-9.18]). Overexpression of C (normal allele) CTD-2378E21.1 in C2C12 skeletal muscle cell, but not A (risk allele) CTD-2378E21.1, showed significantly decreased Kcnj2 expression, indicating A-type CTD-2378E21.1 has lost the ability to regulate Kcnj2. CONCLUSIONS: Our study reveals a shared genetic predisposition between TPP and SPP. CTD-2378E21.1 is a novel disease-associated gene for both TPP and SPP and may negatively regulate KCNJ2 expression. These findings provide new insights into the pathogenesis of nonfamilial hypoKPP.