Genome-wide association study identifies genetic risk loci for adiposity in a Taiwanese population.

Overweight and obese are risk factors for various diseases. In Taiwan, the combined prevalence of overweight and obesity has increased dramatically. Here, we conducted a genome-wide association study (GWAS) on four adiposity traits, including body-mass index (BMI), body fat percentage (BF%), waist circumference (WC), and waist-hip ratio (WHR), using the data for more than 21,000 subjects in Taiwan Biobank. Associations were evaluated between 6,546,460 single-nucleotide polymorphisms (SNPs) and adiposity traits, yielding 13 genome-wide significant (GWS) adiposity-associated trait-loci pairs. A known gene, FTO, as well as two BF%-associated loci (GNPDA2-GABRG1 [4p12] and RNU6-2-PIAS1 [15q23]) were identified as pleiotropic effects. Moreover, RALGAPA1 was found as a specific genetic predisposing factor to high BMI in a Taiwanese population. Compared to other populations, a slightly lower heritability of the four adiposity traits was found in our cohort. Surprisingly, we uncovered the importance of neural pathways that might influence BF%, WC and WHR in the Taiwanese (East Asian) population. Additionally, a moderate genetic correlation between the WHR and BMI (γg = 0.52; p = 2.37×10-9) was detected, suggesting different genetic determinants exist for abdominal adiposity and overall adiposity. In conclusion, the obesity-related genetic loci identified here provide new insights into the genetic underpinnings of adiposity in the Taiwanese population.

Detecting Genetic Ancestry and Adaptation in the Taiwanese Han People.

The Taiwanese people are composed of diverse indigenous populations and the Taiwanese Han. About 95% of the Taiwanese identify themselves as Taiwanese Han, but this may not be a homogeneous population because they migrated to the island from various regions of continental East Asia over a period of 400 years. Little is known about the underlying patterns of genetic ancestry, population admixture, and evolutionary adaptation in the Taiwanese Han people. Here, we analyzed the whole-genome single-nucleotide polymorphism genotyping data from 14,401 individuals of Taiwanese Han collected by the Taiwan Biobank and the whole-genome sequencing data for a subset of 772 people. We detected four major genetic ancestries with distinct geographic distributions (i.e., Northern, Southeastern, Japonic, and Island Southeast Asian ancestries) and signatures of population mixture contributing to the genomes of Taiwanese Han. We further scanned for signatures of positive natural selection that caused unusually long-range haplotypes and elevations of hitchhiked variants. As a result, we identified 16 candidate loci in which selection signals can be unambiguously localized at five single genes: CTNNA2, LRP1B, CSNK1G3, ASTN2, and NEO1. Statistical associations were examined in 16 metabolic-related traits to further elucidate the functional effects of each candidate gene. All five genes appear to have pleiotropic connections to various types of disease susceptibility and significant associations with at least one metabolic-related trait. Together, our results provide critical insights for understanding the evolutionary history and adaption of the Taiwanese Han population.

A functional variant near XCL1 gene improves breast cancer survival via promoting cancer immunity.

Most genome-wide association studies (GWASs) identify genetic variants for breast cancer occurrence. In contrast, few are for recurrence and mortality. We conducted a GWAS on breast cancer survival after diagnosis in estrogen receptor-positive patients, including 953 Taiwanese patients with 159 events. Through Cox proportional hazard models estimation, we identified 24 risk SNPs with p < 1 × 10-5 . Based on imputation and integrated analysis, one SNP, rs1024176 (located in 1q24.2, p = 2.43 × 10-5 ) was found to be a functional variant associated with breast cancer survival and XCL1 gene expression. A series of experimental approaches, including cell-based analyses and CRISPR/Cas9 genome-editing system, were then used and identified the transcription factor MYBL2 was able to discriminately bind to the A allele of rs1024176, the protective variant for breast cancer survival, which promoted XCL1 expression, but not to the G allele of rs1024176. The chemokine XCL1 attracts type 1 dendritic cells (DC1s) to the tumor microenvironment. In breast cancer tissues, we applied a two-step Mendelian randomization analysis, using expression quantitative trait loci as instrumental variables, to confirm higher XCL1 expression was correlated with higher DC1 signatures and favorable disease progression, through the causal effect of rs1024176-A allele. Our study supports the genetic effect on preventing breast cancer survival through XCL1-induced DC1 recruitment in tumor microenvironment.

A high-resolution HLA imputation system for the Taiwanese population: a study of the Taiwan Biobank.

An imputation algorithm for human leukocyte antigen (HLA) is helpful for exploring novel disease associations. However, population-specific HLA imputation references are essential for achieving high imputation accuracy. In this study, a subset of 1012 individuals from the Taiwan Biobank (TWB) who underwent both whole-genome SNP array and NGS-based HLA typing were used to establish Taiwanese HLA imputation references. The HIBAG package was used to generate the imputation references for eight HLA loci at a two- and three-field resolution. Internal validation was carried out to evaluate the call threshold and accuracy for each HLA gene. HLA class II genes found to be associated with rheumatoid arthritis (RA) were validated in this study by the imputed HLA alleles. Our Taiwanese population-specific references achieved average HLA imputation accuracies of 98.11% for two-field and 98.08% for three-field resolution. The frequency distribution of imputed HLA alleles among 23,972 TWB subjects were comparable with PCR-based HLA alleles in general Taiwanese reported in the allele frequency net database. We replicated four common HLA alleles (HLA-DRB1*03:01, DRB1*04:05, DQA1*03:03, and DQB1*04:01) significantly associated with RA. The population-specific references provide an informative tool to investigate the associations of HLA variants and human diseases in large-scale population-based studies.

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.

FGFR2 regulates Mre11 expression and double-strand break repair via the MEK-ERK-POU1F1 pathway in breast tumorigenesis.

The association between breast cancer risk and genetic variants of fibroblast growth factor receptor 2 (FGFR2) has been identified and repeatedly confirmed; however, the mechanism underlying FGFR2 in breast tumorigenesis remains obscure. Given that breast tumorigenesis is particularly related to DNA double-strand-break-repair (DSBR), we examined the hypothesis that FGFR2 is involved in DSBR. Our results show that expression of Mre11, a vital exonuclease in DSBR, is downregulated by FGFR2, which is further linked to decreased DSBR. Analysis of the Mre11 promoter revealed that POU1F1 mediates FGFR2-induced Mre11 downregulation. Furthermore, ERK, downstream of FGFR2, directly interacts with and phosphorylates POU1F1, increasing POU1F1 binding capacity to the Mre11 promoter and repressing Mre11 expression, which consequently affects DSBR and sensitizes breast cancer cells to chemotherapeutic treatments. The importance of the FGFR2-Mre11-DSBR link in cancer progression is suggested by the finding that genotypes of FGFR2 and Mre11 are associated with survival of breast cancer patients and that FGFR2 expression correlates with cancer prognosis specifically in patients receiving chemotherapy. This study yields important insight into the role of FGFR2 in breast tumorigenesis and may facilitate development of a useful therapeutic approach for breast cancer.

A novel estrogen receptor-microRNA 190a-PAR-1-pathway regulates breast cancer progression, a finding initially suggested by genome-wide analysis of loci associated with lymph-node metastasis.

To identify microRNAs that are important in regulating breast cancer progression, the present study used data for the 199 961 single-nucleotide polymorphisms (SNPs) in 837 breast cancer patients genotyped in a recent genome-wide association study to identify loci associated with lymph node metastasis (LNM). SNPs tagging the 15q22.2 locus showed a significant association with LNM and miR-190a was found to be the only microRNA in this region. The role of miR-190a in LNM was supported by the findings that increased miR-190a expression inhibited cell migration and invasiveness and that the target of miR-190a was protease-activated-receptor 1 (PAR-1), which is a metastasis promoting protein in several cancers. In addition, the promoter region of miR-190a was defined and found to contain half of an estrogen response element, suggesting that miR-190a is regulated by estrogen receptor (ER) signaling. This was confirmed by the findings that miR-190a expression was activated by 17ß-estradiol and that ERα bound directly to this promoter. The importance of this ERα-miR190a-PAR-1 link in breast tumorigenesis is suggested by the findings of (i) an association between genetic polymorphism of the miR-190a-containing region and LNM that is modified by SNPs of PAR-1 and is particularly significant in ERα-positive patients and (ii) a combined effect of ERα and miR-190a expression on tumor grade/cancer stage. More importantly, the level of miR-190a expression in primary breast carcinomas correlated with overall survival. These findings suggest a novel pathway in which ERα signaling regulates miR-190a expression, causing inhibition of PAR-1 expression, correlated with inhibition of cancer metastasis.

Functional variants at the 21q22.3 locus involved in breast cancer progression identified by screening of genome-wide estrogen response elements.

INTRODUCTION: Estrogen forms a complex with the estrogen receptor (ER) that binds to estrogen response elements (EREs) in the regulatory region of estrogen-responsive genes and regulates their transcription. Sequence variants in the regulatory regions have the potential to affect the transcription factor-regulatory sequence interaction, resulting in altered expression of target genes. This study explored the association between single-nucleotide polymorphisms (SNPs) within the ERE-associated sequences and breast cancer progression. METHODS: The ERE-associated sequences throughout the whole genome that have been demonstrated to bind ERα in vivo were blasted against online information from SNP data sets and 54 SNPs located adjacent to estrogen-responsive genes were selected for genotyping in two independent cohorts of breast cancer patients: 779 patients in the initial screening stage and another 888 in the validation stage. Deaths due to breast cancer or recurrence of breast cancer were defined as the respective events of interest, and the hazard ratios of individual SNPs were estimated based on the Cox proportional hazards model. Furthermore, functional assays were performed, and information from publicly available genomic data and bioinformatics platforms were used to provide additional evidence for the associations identified in the association analyses. RESULTS: The SNPs at 21q22.3 ERE were significantly associated with overall survival and disease-free survival of patients. Furthermore, these 21q22.3 SNPs (rs2839494 and rs1078272) could affect the binding of this ERE-associated sequence to ERα or Rad21 (an ERα coactivator), respectively, which resulted in a difference in ERα-activated expression of the reporter gene. CONCLUSION: These findings support the idea that functional variants in the ERα-regulating sequence at 21q22.3 are important in determining breast cancer progression.

Phenomic landscape and pharmacogenomic implications for HLA region in a Taiwan Han Chinese population.

BACKGROUND: The human leukocyte antigen (HLA) genes, exhibiting significant genetic diversity, are associated with susceptibility to various clinical diseases and diverse in drug responses. High costs of HLA sequencing and the population-specific architecture of this genetic region necessitate the establishment of a population-specific HLA imputation reference panel. Moreover, there is a lack of understanding about the genetic and phenotypic landscape of HLA variations within the Taiwanese population. METHODS: We created models for a Taiwanese-specific HLA imputation reference panel. These models were trained with the array genotype data and HLA sequencing data from 845 Taiwanese subjects. HLA imputation was applied for 59,448 Taiwanese subjects to characterize the HLA allele and haplotype frequencies. Additionally, a phenome-wide association study (PheWAS) was conducted to identify the phenotypes associated with HLA variations. The association of the biallelic HLA variants with the binary and quantitative traits were evaluated with additive logistic and linear regression models, respectively. Furthermore, an omnibus test with likelihood-ratio test was applied for each HLA amino acid position in the multiallelic HLA amino acid polymorphisms to compare the difference between a fitted model and a null model following a χ2 distribution of n-1 degree of freedom at a position with n residues. Finally, we estimated the prevalence of adverse drug reactions (ADR)-related HLA alleles in the Taiwanese population. RESULTS: In this study, the reference panel models displayed remarkable accuracy, with averages of 99.3%, 98.9%, and 99.1% for 2-, 4-, 6-digit alleles of the eight classical HLA genes, respectively. For PheWAS, a total of 18,136 significant associations with HLA variants across 26 phenotypes are identified (p < 5×10-8), highlighting the pleiotropy feature of the HLA region. Among the independent signals, 15 are novel, including the association of HLA-B pos 138 variation with ankylosing spondylitis (AS), and rs9266290 and rs9266292 with allergy. Through an analysis spanning the entire HLA region, we identified clusters of phenotype correlations. Finally, the carriers of pharmacogenomic related HLA alleles, including HLA-C*01:02 (35.86%), HLA-B*58:01 (20.9%), and HLA-B*15:02 (8.38%), were characterized in the Taiwanese general population. CONCLUSIONS: We successfully delivered the HLA imputation for 59,448 Taiwanese subjects and characterized the genetic and phenotypic landscapes of the HLA variations. In addition, we quantified the estimated prevalence of the ADR-related HLA alleles in the Taiwanese population. The developed HLA imputation reference panel could be used for estimation of population HLA allele frequencies, which can facilitate further studies in the role of HLA variants in a wider range of phenotypes in the population.

Cancer carrier screening in the general population using whole-genome sequencing.

BACKGROUND: Cancer is a major cause of death, and its early identification and intervention have potential for clinical actionability and benefits for human health. The studies using whole-genome sequencing (WGS) and large samples analysis of cancer-related genes have been rarely done. METHODS: We performed WGS to explore germline mutations in coding and non-coding areas of cancer-related genes and non-coding driver genes and regulatory areas. Structural variants (SVs) was also analyzed. We used several tools and a subgrouping method to analyze the variants in 1491 healthy participants. Moreover, 275 cancer-related genes sequencing was carried out in 125 cancer patients. RESULTS: The incidence of familial cancer in the Taiwanese general population is 8.79% (131/1491). Cancer carrier rate of cancer-related genes is about 7.04% (105/1491) for pathogenic/likely pathogenic variants (P/LP) on ClinVar database only, and 28.24% (421/1491) for P/LP and loss of function variants. The carrier frequencies of cancer-related genes P/LP on ClinVar database were as follows: 8.40% (11/131), 7.11% (28/394), and 6.83% (66/966) in FC, 1MC, and nMC, respectively. The SVs and non-coding driver gene variants are uncommon. There are 1.54% (23/1491) of actionable cancer genes in American College of Medical Genetics and Genomics (ACMG), and the germline mutation rate of 275 cancer-related genes is 7.2% (9/125) in cancer patients including 4.0% (5/125) of actionable cancer genes in ACMG. After analyzing the frequencies of P/LP variants on GJB2 and SLC25A13 genes, we suggest that these two genes may not be cancer-related genes and need be re-evaluated. CONCLUSIONS: WGS analysis can completely detect germline mutations in cancer carriers. This study use subgrouping approach for samples provides a strategy to study whether a gene or variant is a cancer-related gene or variant in the future studies.

MicroRNA-30a inhibits cell migration and invasion by downregulating vimentin expression and is a potential prognostic marker in breast cancer.

Tumor recurrence and metastasis result in an unfavorable prognosis for cancer patients. Recent studies have suggested that specific microRNAs (miRNAs) may play important roles in the development of cancer cells. However, prognostic markers and the outcome prediction of the miRNA signature in breast cancer patients have not been comprehensively assessed. The aim of this study was to identify miRNA biomarkers relating to clinicopathological features and outcome of breast cancer. A miRNA microarray analysis was performed on breast tumors of different lymph node metastasis status and with different progression signatures, indicated by overexpression of cyclin D1 and ß-catenin genes, to identify miRNAs showing a significant difference in expression. The functional interaction between the candidate miRNA, miR-30a, and the target gene, Vim, which codes for vimentin, a protein involved in epithelial-mesenchymal transition, was examined using the luciferase reporter assay, western blotting, and migration and invasion assays. The association between the decreased miR-30a levels and breast cancer progression was examined in a survival analysis. miR-30a negatively regulated vimentin expression by binding to the 3'-untranslated region of Vim. Overexpression of miR-30a suppressed the migration and invasiveness phenotypes of breast cancer cell lines. Moreover, reduced tumor expression of miR-30a in breast cancer patients was associated with an unfavorable outcome, including late tumor stage, lymph node metastasis, and worse progression (mortality and recurrence) (p < 0.05). In conclusion, these findings suggest a role for miR-30a in inhibiting breast tumor invasiveness and metastasis. The finding that miR-30a downmodulates vimentin expression might provide a therapeutic target for the treatment of breast cancer.

Taiwan Biobank: A rich biomedical research database of the Taiwanese population.

The Taiwan Biobank (TWB) is an ongoing prospective study of >150,000 individuals aged 20-70 in Taiwan. A comprehensive list of phenotypes was collected for each consented participant at recruitment and follow-up visits through structured interviews and physical measurements. Biomarkers and genetic data were generated from blood and urine samples. We present here an overview of TWB's genetic data quality, population structure, and familial relationship, which consists of predominantly Han Chinese ancestry, and highlight its important attributes and genetic findings thus far. A linkage to Taiwan's National Health Insurance database of >25 years and other registries is underway to enrich the phenotypic spectrum and enable deep and longitudinal genetic investigations. TWB provides one of the largest biobank resources for biomedical and public health research in East Asia that will contribute to our understanding of the genetic basis of health and disease in global populations through collaborative studies with other biobanks.

Genomic and Metabolomic Landscape of Right-Sided and Left-Sided Colorectal Cancer: Potential Preventive Biomarkers.

Colorectal cancer (CRC) is the third most common cancer worldwide. The incidence and mortality rates of CRC are significantly higher in Taiwan than in other developed countries. Genes involved in CRC tumorigenesis differ depending on whether the tumor occurs on the left or right side of the colon, and genomic analysis is a keystone in the study and treatment of CRC subtypes. However, few studies have focused on the genetic landscape of Taiwanese patients with CRC. This study comprehensively analyzed the genomes of 141 Taiwanese patients with CRC through whole-exome sequencing. Significant genomic differences related to the site of CRC development were observed. Blood metabolomic profiling and polygenic risk score analysis were performed to identify potential biomarkers for the early identification and prevention of CRC in the Taiwanese population. Our findings provide vital clues for establishing population-specific treatments and health policies for CRC prevention in Taiwan.

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.

The Causal Relationship of Circulating Triglyceride and Glycated Hemoglobin: A Mendelian Randomization Study.

CONTEXT: The association between circulating triglyceride (TG) and glycated hemoglobin A1c (HbA1c), a biomarker for type 2 diabetes, has been widely addressed, but the causal direction of the relationship is still ambiguous. OBJECTIVE: To confirm the causal relationship between TG and HbA1c by using bidirectional and 2-step Mendelian randomization (MR) approaches. METHODS: We carried out a bidirectional MR approach using the summarized results from the public database to examine any potential causal effects between serum TG and HbA1c in 16 000 individuals of the Taiwan Biobank cohort. We used the MR estimate and the MR inverse variance-weighted method to reveal that relationship between TG and HbA1c. To further determine whether the DNA methylation at specific sequences mediate the causal pathway between TG and HbA1c, using the 2-step MR approach. RESULTS: We identified that a single-unit increase in TG measured via log transformation of mg/dL data was associated with a significant increase of 10 units of HbA1c (95% CI = 1.05-18.95, P = 0.029). In contrast, the genetic determinants of HbA1c do not contribute to the amount of circulating TG (beta = 1.75, 95% CI = -11.50 to 14.90). Sensitivity analyses, included the weighted-median approach and MR-Egger regression, were performed to confirm no pleiotropic effect among these instrumental variables. Furthermore, we identified the genetic variant, rs1823200, is associated with both methylation of the CpG site adjacent to CADPS gene and HbA1c level. CONCLUSION: Our study suggests that higher circulating TG can have an affect on genomic methylation status, ultimately causing elevated level of circulating HbA1c.

Blood multiomics reveal insights into population clusters with low prevalence of diabetes, dyslipidemia and hypertension.

Diabetes, dyslipidemia and hypertension are important metabolic diseases that impose a great burden on many populations worldwide. However, certain population strata have reduced prevalence for all three diseases, but the underlying mechanisms are poorly understood. We sought to identify the phenotypic, genomic and metabolomic characteristics of the low-prevalence population to gain insights into possible innate non-susceptibility against metabolic diseases. We performed k-means cluster analysis of 16,792 subjects using anthropometric and clinical biochemistry data collected by the Taiwan Biobank. Nuclear magnetic resonance spectra-based metabolome analysis was carried out for 217 subjects with normal body mass index, good exercise habits and healthy lifestyles. We found that the gene APOA5 was significantly associated with reduced prevalence of disease, and lesser associations included the genes HIF1A, LIMA1, LPL, MLXIPL, and TRPC4. Blood plasma of subjects belonging to the low disease prevalence cluster exhibited lowered levels of the GlycA inflammation marker, very low-density lipoprotein and low-density lipoprotein cholesterol, triglycerides, valine and leucine compared to controls. Literature mining revealed that these genes and metabolites are biochemically linked, with the linkage between lipoprotein metabolism and inflammation being particularly prominent. The combination of phenomic, genomic and metabolomic analysis may also be applied towards the study of metabolic disease prevalence in other populations.

Large-scale genome-wide association study identifies HLA class II variants associated with chronic HBV infection: a study from Taiwan Biobank.

BACKGROUND: Chronic hepatitis B virus (HBV) infection is a great health burden with geographical variations. AIMS: To explore genetic variants associated with chronic HBV infection. METHODS: The study included 15 352 participants seropositive for HBV core antibodies in Taiwan Biobank. Among them, 2591 (16.9%) seropositive for HBV surface antigen (HBsAg) were defined as having chronic HBV infection. All participants were examined for whole-genome genotyping by Axiom-Taiwan Biobank Array. The human leucocyte antigen (HLA) imputation was performed after identification of the variants within the region. Logistic regressions were used to estimate odds ratios (ORs) with 95% confidence intervals. Correlations of different HLA allele frequencies with HBsAg seroprevalence were evaluated across worldwide populations by Pearson correlation coefficients. Epitope prediction was performed for HLA alleles using NetMHCIIpan method. RESULTS: Located within a cluster of 450 single nucleotide polymorphisms in HLA class II, rs7770370 (P = 2.73 × 10-35 ) was significantly associated with HBV chronicity (Pcorrected  < 8.6 × 10-8 ). Imputation analyses showed that HLA-DPA1*02:02 and HLA-DPB1*05:01 were associated with chronic HBV, with adjusted ORs of 1.43 (1.09-1.89) and 1.61 (1.29-2.01). These allele frequencies were positively correlated with global HBsAg seroprevalence, with R of 0.75 and 0.62 respectively (P < 0.05). HLA-DRB1*13:02, HLA-DQA1* 01:02 and HLA-DQB1*06:09 associated with HBV chronicity negatively, with adjusted ORs of 0.31 (0.17-0.58), 0.70 (0.56-0.87) and 0.33 (0.18-0.63). These HLA alleles had various binding affinities to the predicted epitopes derived from HBV nucleocapsid protein. CONCLUSIONS: HLA class II variants are relevant for chronicity after HBV acquisition.

The rs1256328 (ALPL) and rs12654812 (RGS14) Polymorphisms are Associated with Susceptibility to Calcium Nephrolithiasis in a Taiwanese population.

Nephrolithiasis is a common disease affecting almost all populations, with an increasing prevalence over the past decades. Previous studies revealed several functional polymorphisms associated with the pathogenesis of nephrolithiasis. However, data on Asian populations are limited. In this study, three candidate polymorphisms were selected from previous studies to investigate the correlations with nephrolithiasis in a Taiwanese population. In total, 454 nephrolithiasis patients were recruited from Kaohsiung Medical University Hospital, with SNP frequency for 1513 subjects of general population from the Taiwan Biobank (TWB) as a genotypic reference. Results revealed that subjects with minor TT genotype at rs1256328 (alkaline phosphatase, liver/bone/kidney (ALPL)) have higher susceptibility to nephrolithiasis (odds ratio (OR) = 2.03, p = 0.0013). In addition, subjects carrying the minor AA genotype at rs12654812 (regulator of G protein signaling 14 (RGS14)) have higher susceptibility to nephrolithiasis (OR = 1.91, p = 0.0017). Among nephrolithiasis patients, subjects with GG at rs7627468 (calcium-sensing receptor (CASR)) have lower pH level in urine (p = 0.0088). Importantly, rs7627468 is associated with the expressions of IQCB1 and EAF2. rs12654812 could influence the expression of RGS14 itself, MXD3, and FGFR4. In summary, this study successfully validated the genetic roles of rs1256328 and rs12654812 in human nephrolithiasis.

Meta-Analysis of Genome-Wide Association Studies Identifies Three Loci Associated With Stiffness Index of the Calcaneus.

The stiffness index (SI) from quantitative ultrasound measurements is a good indicator of BMD and may be used to predict the risk of osteoporotic fracture. We conducted a genomewide association study (GWAS) for SI using 7742 individuals from the Taiwan Biobank, followed by a replication study in a Korean population (n = 2955). Approximately 6.1 million SNPs were subjected to association analysis, and SI-associated variants were identified. We further conducted a meta-analysis of Taiwan Biobank significant SNPs with a Korean population-based cohort. Candidate genes were prioritized according to epigenetic annotations, gene ontology, protein-protein interaction, GWAS catalog, and expression quantitative trait loci analyses. Our results revealed seven significant single-nucleotide polymorphisms (SNPs) within three loci: 7q31.31, 17p13.3, and 11q14.2. Conditional analysis showed that three SNPs, rs2536195 (CPED1/WNT16), rs1231207 (SMG6), and rs4944661 (LOC10050636/TMEM135), were the most important signals within these regions. The associations for the three SNPs were confirmed in a UK Biobank estimated BMD GWAS; these three cytobands were replicated successfully after a meta-analysis with a Korean population cohort as well. However, two SNPs were not replicated. After prioritization, we identified two novel genes, RAB15 and FNTB, as strong candidates for association with SI. Our study identified three SI-associated SNPs and two novel SI-related genes. Overall, these results provide further insight into the genetic architecture of osteoporosis. Further studies in larger East Asian populations are needed. © 2019 American Society for Bone and Mineral Research.