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.
Melatonin exerts a wide range of effects among various tissues and organs. However, there is currently no study to investigate the genetic determinants of melatonin secretion. Here, we conducted a ...genome-wide association study (GWAS) for melatonin secretion using morning urine 6-hydroxymelatonin sulfate-to-creatinine ratio (UMCR). We initially enrolled 5000 participants from Taiwan Biobank in this study. After excluding individuals that did not have their urine collected in the morning, those who had history of neurological or psychiatric disorder, and those who failed to pass quality control, association of single nucleotide polymorphisms with log-transformed UMCR adjusted for age, sex and principal components of ancestry were analyzed. A second model additionally adjusted for estimated glomerular filtration rate (eGFR). A total of 2373 participants underwent the genome-wide analysis. Five candidate loci associated with log UMCR (P value ranging from 6.83 × 10
to 3.44 × 10
) encompassing ZFHX3, GALNT15, GALNT13, LDLRAD3 and intergenic between SEPP1 and FLJ32255 were identified. Similar results were yielded with further adjustment for eGFR. Interestingly, the identified genes are associated with circadian behavior, neuronal differentiation, motor disorders, anxiety, and neurodegenerative diseases. We conducted the first GWAS for melatonin secretion and identified five candidate genetic loci associated with melatonin level. Replication and functional studies are needed in the future.
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.
What's new?
Although many prognostic predictors have been identified, the factors that determine patient differences in breast cancer progression and survival are not fully understood. This study addresses the important role of cancer immune control in promoting human breast cancer survival. Using genetic, cell‐based, and multi‐omic analyses, the authors find that one of the most significant SNP, rs1024176 at 1q24.2 could be a functional variant regulating XCL1 gene expression and thus recruiting type 1 dendritic cells to the breast tumor microenvironment. These findings reveal a pathway to prevent breast cancer progression through XCL1‐induced type 1 dendritic cells recruitment in the tumor microenvironment.
In a three-stage genome-wide association study among East Asian women including 22,780 cases and 24,181 controls, we identified 3 genetic loci newly associated with breast cancer risk, including ...rs4951011 at 1q32.1 (in intron 2 of the ZC3H11A gene; P=8.82×10(-9)), rs10474352 at 5q14.3 (near the ARRDC3 gene; P=1.67×10(-9)) and rs2290203 at 15q26.1 (in intron 14 of the PRC1 gene; P=4.25×10(-8)). We replicated these associations in 16,003 cases and 41,335 controls of European ancestry (P=0.030, 0.004 and 0.010, respectively). Data from the ENCODE Project suggest that variants rs4951011 and rs10474352 might be located in an enhancer region and transcription factor binding sites, respectively. This study provides additional insights into the genetics and biology of breast cancer.
Abstract
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.
Fibroblast growth factors (FGFs) 21 and 23 are recently identified hormones regulating metabolism of glucose, lipid, phosphate and vitamin D. Here we conducted a genome-wide association study (GWAS) ...for circulating FGF21 and FGF23 concentrations to identify their genetic determinants. We enrolled 5,000 participants from Taiwan Biobank for this GWAS. After excluding participants with diabetes mellitus and quality control, association of single nucleotide polymorphisms (SNPs) with log-transformed FGF21 and FGF23 serum concentrations adjusted for age, sex and principal components of ancestry were analyzed. A second model additionally adjusted for body mass index (BMI) and a third model additionally adjusted for BMI and estimated glomerular filtration rate (eGFR) were used. A total of 4,201 participants underwent GWAS analysis. rs67327215, located within RGS6 (a gene involved in fatty acid synthesis), and two other SNPs (rs12565114 and rs9520257, located between PHC2-ZSCAN20 and ARGLU1-FAM155A respectively) showed suggestive associations with serum FGF21 level (P = 6.66 × 10
, 6.00 × 10
and 6.11 × 10
respectively). The SNPs rs17111495 and rs17843626 were significantly associated with FGF23 level, with the former near PCSK9 gene and the latter near HLA-DQA1 gene (P = 1.04 × 10
and 1.80 × 10
respectively). SNP rs2798631, located within the TGFB2 gene, was suggestively associated with serum FGF23 level (P = 4.97 × 10
). Additional adjustment for BMI yielded similar results. For FGF23, further adjustment for eGFR had similar results. We conducted the first GWAS of circulating FGF21 levels to date. Novel candidate genetic loci associated with circulating FGF21 or FGF23 levels were found. Further replication and functional studies are needed to support our findings.
The key signature of cancer genomes is the accumulation of DNA mutations, the most abundant of which is the cytosine-to-thymine (C-to-T) transition that results from cytosine deamination. Analysis of ...The Cancer Genome Atlas (TCGA) database has demonstrated that this transition is caused mainly by upregulation of the cytosine deaminase APOBEC3B (A3B), but the mechanism has not been completely characterized. We found that B-Myb (encoded by MYBL2) binds the A3B promoter, causing transactivation, and this is responsible for the C-to-T transitions and DNA hypermutation in breast cancer cells. Analysis of TCGA database yielded similar results, supporting that MYBL2 and A3B are upregulated and putatively promote C-to-T transitions in multiple cancer types. Moreover, blockade of EGF receptor with afatinib attenuated B-Myb-A3B signaling, suggesting a clinically relevant means of suppressing mutagenesis. Our results suggest that B-Myb-A3B contributes to DNA damage and could be targeted by inhibiting EGF receptor.
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.
MicroRNAs play critical roles in regulating various physiological processes, including growth and development. Previous studies have shown that microRNA-124 (miR-124) participates not only in ...regulation of early neurogenesis but also in suppression of tumorigenesis. In the present study, we found that overexpression of miR-124 was associated with reduced DNA repair capacity in cultured cancer cells and increased sensitivity of cells to DNA-damaging anti-tumor drugs, specifically those that cause the formation of DNA strand-breaks (SBs). We then examined which DNA repair-related genes, particularly the genes of SB repair, were regulated by miR-124. Two SB repair-related genes, encoding ATM interactor (ATMIN) and poly (ADP-ribose) polymerase 1 (PARP1), were strongly affected by miR-124 overexpression, by binding of miR-124 to the 3¢-untranslated region of their mRNAs. As a result, the capacity of cells to repair DNA SBs, such as those resulting from homologous recombination, was significantly reduced upon miR-124 overexpression. A particularly important therapeutic implication of this finding is that overexpression of miR-124 enhanced cell sensitivity to multiple DNA-damaging agents via ATMIN- and PARP1-mediated mechanisms. The translational relevance of this role of miR-124 in anti-tumor drug sensitivity is suggested by the finding that increased miR-124 expression correlates with better breast cancer prognosis, specifically in patients receiving chemotherapy. These findings suggest that miR-124 could potentially be used as a therapeutic agent to improve the efficacy of chemotherapy with DNA-damaging agents.
To evaluate the association between serum urate and risk of incident chronic kidney disease (CKD) and to assess whether serum urate plays a causal role in CKD.
We conducted a prospective cohort study ...and Mendelian randomization analysis that analyzed longitudinal data from the Taiwan Biobank between January 1, 2012, and December 31, 2021.
A total of 34,831 individuals met the inclusion criteria, of which 4697 (13.5%) had hyperuricemia. After a median (interquartile range) follow-up of 4.1 (3.1-4.9) years, 429 participants developed CKD. After adjustment for age, sex, and comorbid conditions, each mg/dL increase in serum urate was associated with a 15% higher risk of incident CKD (HR, 1.15; 95% CI, 1.08 to 1.24; P<.001). The genetic risk score and seven Mendelian randomization methods revealed no significant association between serum urate levels and the risk of incident CKD (HR, 1.03; 95% CI, 0.72 to 1.46; P=0.89; all P>.05 for 7 Mendelian randomization methods).
This prospective, population-based cohort study showed that elevated serum urate is a significant risk factor for incident CKD; however, Mendelian randomization analyses failed to provide evidence that serum urate had a causal effect on CKD in the East Asian population.