Genome-wide meta-analysis in Japanese populations identifies novel variants at the TMC6-TMC8 and SIX3-SIX2 loci associated with HbA1c

Comparison of the loci associated with HbA1c and blood glucose levels identified by a genome-wide association study in the Japanese population

Aims

Hemoglobin A1c (HbA1c) levels are widely employed to diagnose diabetes. However, estimates of the heritability of HbA1c and glucose levels are different. Therefore, we explored HbA1c- and blood glucose-associated loci in a non-diabetic Japanese population.

Methods

We conducted a two-stage genome-wide association study (GWAS) on variants associated with HbA1c and blood glucose levels in a Japanese population. In the initial stage, data of 4911 participants of the Japan Multi-Institutional Collaborative Cohort (J-MICC) were subjected to discovery analysis. In the second stage, two datasets from the Tohoku Medical Megabank project, with 8175 and 40,519 participants, were used for the replication study. Association of the imputed variants with HbA1c and blood glucose levels was determined via linear regression analyses adjusted for age, sex, body mass index (BMI), smoking, and genetic principal components (PC1-PC10). Moreover, we performed a BMI-stratified GWAS on HbA1c levels in the J-MICC. The discovery analysis and BMI-stratified GWAS results were validated with re-analyses of normalized HbA1c levels adjusted for site in addition to the above, and blood glucose adjusted for fasting time as an additional covariate.

Results

Genetic variants associated with HbA1c levels were identified in KCNQ1 and TMC6. None of the genetic variants associated with blood glucose levels in the discovery analysis were replicated. Association of rs2299620 in KCNQ1 with HbA1c levels showed heterogeneity between individuals with BMI ≥ 25 kg/m2 and BMI < 25 kg/m2.

Conclusions

The variant rs2299620 in KCNQ1 might affect HbA1c levels differentially based on BMI grouping in the Japanese population.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13340-023-00618-0.

Impact of diet and host genetics on the murine intestinal mycobiome

The mammalian gut is home to a diverse microbial ecosystem, whose composition affects various physiological traits of the host. Next-generation sequencing-based metagenomic approaches demonstrated how the interplay of host genetics, bacteria, and environmental factors shape complex traits and clinical outcomes. However, the role of fungi in these complex interactions remains understudied. Here, using 228 males and 363 females from an advanced-intercross mouse line, we provide evidence that fungi are regulated by host genetics. In addition, we map quantitative trait loci associated with various fungal species to single genes in mice using whole genome sequencing and genotyping. Moreover, we show that diet and its' interaction with host genetics alter the composition of fungi in outbred mice, and identify fungal indicator species associated with different dietary regimes. Collectively, in this work, we uncover an association of the intestinal fungal community with host genetics and a regulatory role of diet in this ecological niche.

Development and validation of genome-wide polygenic risk scores for predicting breast cancer incidence in Japanese females: a population-based case-cohort study

PURPOSE

This study aimed to develop an ancestry-specific polygenic risk scores (PRSs) for the prediction of breast cancer events in Japanese females and validate it in a longitudinal cohort study.

METHODS

Using publicly available summary statistics of female breast cancer genome-wide association study (GWAS) of Japanese and European ancestries, we, respectively, developed 31 candidate genome-wide PRSs using pruning and thresholding (P + T) and LDpred methods with varying parameters. Among the candidate PRS models, the best model was selected using a case-cohort dataset (63 breast cancer cases and 2213 sub-cohorts of Japanese females during a median follow-up of 11.9 years) according to the maximal predictive ability by Harrell's C-statistics. The best-performing PRS for each derivation GWAS was evaluated in another independent case-cohort dataset (260 breast cancer cases and 7845 sub-cohorts of Japanese females during a median follow-up of 16.9 years).

RESULTS

For the best PRS model involving 46,861 single nucleotide polymorphisms (SNPs; P + T method with P_T = 0.05 and R² = 0.2) derived from Japanese-ancestry GWAS, the Harrell's C-statistic was 0.598 ± 0.018 in the evaluation dataset. The age-adjusted hazard ratio for breast cancer in females with the highest PRS quintile compared with those in the lowest PRS quintile was 2.47 (95% confidence intervals, 1.64-3.70). The PRS constructed using Japanese-ancestry GWAS demonstrated better predictive performance for breast cancer in Japanese females than that using European-ancestry GWAS (Harrell's C-statistics 0.598 versus 0.586).

CONCLUSION

This study developed a breast cancer PRS for Japanese females and demonstrated the usefulness of the PRS for breast cancer risk stratification.

Concordance between fasting plasma glucose and HbA1c in the diagnosis of diabetes in black South African adults: a cross-sectional study

OBJECTIVES

We investigated concordance between haemoglobin A1c (HbA₁c)-defined diabetes and fasting plasma glucose (FPG)-defined diabetes in a black South African population with a high prevalence of obesity.

DESIGN

Cross-sectional study.

SETTING

Rural South African population-based cohort.

PARTICIPANTS

765 black individuals aged 40-70 years and with no history of diabetes.

PRIMARY AND SECONDARY OUTCOME MEASURES

The primary outcome measure was concordance between HbA_1c-defined diabetes and FPG-defined diabetes. Secondary outcome measures were differences in anthropometric characteristics, fat distribution and insulin resistance (measured using Homoeostatic Model Assessment of Insulin Resistance (HOMA-IR)) between those with concordant and discordant HbA_1c/FPG classifications and predictors of HbA_1c variance.

RESULTS

The prevalence of HbA_1c-defined diabetes was four times the prevalence of FPG-defined diabetes (17.5% vs 4.2%). Classification was discordant in 15.7% of participants, with 111 individuals (14.5%) having HbA_1c-only diabetes (kappa 0.23; 95% CI 0.14 to 0.31). Median body mass index, waist and hip circumference, waist-to-hip ratio, subcutaneous adipose tissue and HOMA-IR in participants with HbA_1c-only diabetes were similar to those in participants who were normoglycaemic by both biomarkers and significantly lower than in participants with diabetes by both biomarkers (p<0.05). HOMA-IR and fat distribution explained additional HbA_1c variance beyond glucose and age only in women.

CONCLUSIONS

Concordance was poor between HbA_1c and FPG in diagnosis of diabetes in black South Africans, and participants with HbA_1c-only diabetes phenotypically resembled normoglycaemic participants. Further work is necessary to determine which of these parameters better predicts diabetes-related morbidities in this population and whether a population-specific HbA_1c threshold is necessary.

SIX2 Regulates Human β Cell Differentiation from Stem Cells and Functional Maturation In Vitro

Generation of insulin-secreting β cells in vitro is a promising approach for diabetes cell therapy. Human embryonic stem cells (hESCs) and human induced pluripotent stem cells (hiPSCs) are differentiated to β cells (SC-β cells) and mature to undergo glucose-stimulated insulin secretion, but molecular regulation of this defining β cell phenotype is unknown. Here, we show that maturation of SC-β cells is regulated by the transcription factor SIX2. Knockdown (KD) or knockout (KO) of SIX2 in SC-β cells drastically limits glucose-stimulated insulin secretion in both static and dynamic assays, along with the upstream processes of cytoplasmic calcium flux and mitochondrial respiration. Furthermore, SIX2 regulates the expression of genes associated with these key β cell processes, and its expression is restricted to endocrine cells. Our results demonstrate that expression of SIX2 influences the generation of human SC-β cells in vitro.

Velazco-Cruz et al. characterize the role of SIX2 in stem cell differentiation to β cells. SIX2 expression is restricted to late-stage endocrine cells. Generation of β cells does not require SIX2, but lack of SIX2 impairs maturation, as assessed by glucose-stimulated insulin secretion, calcium flux, mitochondrial respiration, and gene expression.

CRISPR-based genome editing in primary human pancreatic islet cells

Gene targeting studies in primary human islets could advance our understanding of mechanisms driving diabetes pathogenesis. Here, we demonstrate successful genome editing in primary human islets using clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9). CRISPR-based targeting efficiently mutated protein-coding exons, resulting in acute loss of islet β-cell regulators, like the transcription factor PDX1 and the KATP channel subunit KIR6.2, accompanied by impaired β-cell regulation and function. CRISPR targeting of non-coding DNA harboring type 2 diabetes (T2D) risk variants revealed changes in ABCC8, SIX2 and SIX3 expression, and impaired β-cell function, thereby linking regulatory elements in these target genes to T2D genetic susceptibility. Advances here establish a paradigm for genetic studies in human islet cells, and reveal regulatory and genetic mechanisms linking non-coding variants to human diabetes risk.

Modeling different types of diabetes using human pluripotent stem cells

Diabetes mellitus (DM) is a metabolic disease characterized by chronic hyperglycemia as a result of progressive loss of pancreatic β cells, which could lead to several debilitating complications. Different paths, triggered by several genetic and environmental factors, lead to the loss of pancreatic β cells and/or function. Understanding these many paths to β cell damage or dysfunction could help in identifying therapeutic approaches specific for each path. Most of our knowledge about diabetes pathophysiology has been obtained from studies on animal models, which do not fully recapitulate human diabetes phenotypes. Currently, human pluripotent stem cell (hPSC) technology is a powerful tool for generating in vitro human models, which could provide key information about the disease pathogenesis and provide cells for personalized therapies. The recent progress in generating functional hPSC-derived β cells in combination with the rapid development in genomic and genome-editing technologies offer multiple options to understand the cellular and molecular mechanisms underlying the development of different types of diabetes. Recently, several in vitro hPSC-based strategies have been used for studying monogenic and polygenic forms of diabetes. This review summarizes the current knowledge about different hPSC-based diabetes models and how these models improved our current understanding of the pathophysiology of distinct forms of diabetes. Also, it highlights the progress in generating functional β cells in vitro, and discusses the current challenges and future perspectives related to the use of the in vitro hPSC-based strategies.

SIX2 and SIX3 coordinately regulate functional maturity and fate of human pancreatic β cells

The physiological functions of many vital tissues and organs continue to mature after birth, but the genetic mechanisms governing this postnatal maturation remain an unsolved mystery. Human pancreatic β cells produce and secrete insulin in response to physiological cues like glucose, and these hallmark functions improve in the years after birth. This coincides with expression of the transcription factors SIX2 and SIX3, whose functions in native human β cells remain unknown. Here, we show that shRNA-mediated SIX2 or SIX3 suppression in human pancreatic adult islets impairs insulin secretion. However, transcriptome studies revealed that SIX2 and SIX3 regulate distinct targets. Loss of SIX2 markedly impaired expression of genes governing β-cell insulin processing and output, glucose sensing, and electrophysiology, while SIX3 loss led to inappropriate expression of genes normally expressed in fetal β cells, adult α cells, and other non-β cells. Chromatin accessibility studies identified genes directly regulated by SIX2. Moreover, β cells from diabetic humans with impaired insulin secretion also had reduced SIX2 transcript levels. Revealing how SIX2 and SIX3 govern functional maturation and maintain developmental fate in native human β cells should advance β-cell replacement and other therapeutic strategies for diabetes.

The KLF14 Variant is Associated with Type 2 Diabetes and HbA1C Level

The purpose of this study was to scan variants in coding region of Krȕppel like factor14 (KLF14) locus and assess association related to type 2 diabetes (T2D) in Iranian population. We sequenced the coding region of KLF14 to scan variants in case-sibling study (92 individuals with T2D and 92 healthy older siblings). To confirm, we analyzed rs76603546 association with T2D in a larger unrelated case-control study by PCR-RFLP (475 cases and 512 controls). We analyzed the association of rs76603546 with HbA_1C, BMI, fat mass, waist circumference, fasting glucose, cholesterol and HOMA-IR (Homeostatic Model Assessment for Insulin Resistance) using one-way ANOVA analysis. Also, association of genotypes with T2D adjusted for confounding variables was analyzed using logistic regression. HaploReg v 4.1 was used to predict rs76603546 possible function. Sequencing results analysis revealed the association of C allele of rs76603546, synonymous variant C>T, [OR 2.10 (1.38-3.20), P value < 0.001] and CC genotype of rs76603546 [OR 4.3 (1.79-10.23), P value = 0.001] with susceptibility to T2D. PCR-Restriction Fragment Length Polymorphism (RFLP) results analysis confirmed the association of rs76603546 with T2D [C allele, OR 1.91 (1.59-2.29), P value = 0.002, CC genotype, OR 3.27 (2.26-4.73), P value = 0.002 and TC genotype, OR 1.74 (1.31-2.31), P value = 0.001]. The CC genotype of rs76603546 is associated with HbA_1C level (P value < 0.001) and BMI (P value = 0.02). After adjustment with confounding variables, we observed association of CC genotype with T2D [OR 2.542 (1.25-3.77), P value = 0.03]. Among over 220 SNPs, rs76603546 was associated with T2D, HbA_1C and BMI in our study.

A genome-wide association study for highly sensitive cardiac troponin T levels identified a novel genetic variation near a RBAK-ZNF890P locus in the Japanese general population

BACKGROUND

Cardiovascular disease (CVD) is a major cause of mortality worldwide. High-sensitivity cardiac troponin T (hs-cTnT) is released into the bloodstream due to cardiomyocyte damage and is associated with a high CVD risk. This study aimed to investigate hs-cTnT-related genetic variation and to examine whether this is an associated risk factor for CVD in the Japanese general population.

METHODS

This was a genome-wide association study (GWAS) based on a cohort from the 2013 Tohoku Medical Megabank Project community study. The GWAS was performed using a HumanOmniExpressExome BeadChip array with 914,035 autosomal single-nucleotide polymorphisms. The Framingham Risk Score and the Suita score were used to evaluate the future risk of CVD.

RESULTS

The GWAS identified 10 loci reaching suggestive significance in the discovery cohort. A replication analysis confirmed that one of the 10 loci, rs7798496, is associated with elevated hs-cTnT levels. The combined P value in the discovery and replication cohorts for the association between the rs7798496 and hs-cTnT levels was 3.4 × 10^-8, which indicates that the novel variant reached genome-wide significance. The rs7798496 loci was located at an intergenic region between the retinoblastoma gene product (RB)-associated Krüppell-associated box (KRAB) zinc finger, zinc finger protein 890, and pseudogene (ZNF890P). Logistic regression analysis revealed that the presence of the rs7798496 T allele was strongly associated with a high risk for CVD.

CONCLUSIONS

This study provides insights into a link between a novel genetic variant, T allele of rs7798269, and elevated hs-cTnT levels as a future risk for CVD in the general Japanese population.

A Genome-Wide Association Study Identifies the Association between the 12q24 Locus and Black Tea Consumption in Japanese Populations

Several genome-wide association studies (GWASs) have reported the association between genetic variants and the habitual consumption of foods and drinks; however, no association data are available regarding the consumption of black tea. The present study aimed to identify genetic variants associated with black tea consumption in 12,258 Japanese participants. Data on black tea consumption were collected by a self-administered questionnaire, and genotype data were obtained from a single nucleotide polymorphism array. In the discovery GWAS, two loci met suggestive significance (p < 1.0 × 10^-6). Three genetic variants (rs2074356, rs144504271, and rs12231737) at 12q24 locus were also significantly associated with black tea consumption in the replication stage (p < 0.05) and during the meta-analysis (p < 5.0 × 10^-8). The association of rs2074356 with black tea consumption was slightly attenuated by the additional adjustment for alcohol drinking frequency. In conclusion, genetic variants at the 12q24 locus were associated with black tea consumption in Japanese populations, and the association is at least partly mediated by alcohol drinking frequency.

ALDH2 genotype modulates the association between alcohol consumption and AST/ALT ratio among middle-aged Japanese men: a genome-wide G × E interaction analysis

Liver tests (LT), especially to measure AST, ALT and GGT levels, are widely used to evaluate the risk of alcohol-related liver disease (ALD). In this study, we investigated the potential genetic factors that modulate the association between LTs and alcohol consumption. We conducted a genome-wide interaction meta-analysis in 7856 Japanese subjects from Tohoku Medical Megabank Community-Based Cohort (TMM CommCohort) study recruited in 2013, and identified 2 loci (12q24 and 2p16) with genome-wide significance (P > 5 × 10-8). The significant variants in the 12q24 included rs671, a variant associated with alcohol intolerance and located at a coding exon of ALDH2. We found that the amount of alcohol consumption was associated with increased level AST/ALT ratio among the subjects with the rs671 GA genotype. The elevated AST/ALT ratio among subjects with moderate-to-high levels of drinking behavior and the rs671 GA genotype was due to decreased levels of ALT, which was not accompanied with significant differences in AST levels. Although the interaction effect was significant in both men and women, the effect was much larger in men. Our results suggest that the impact of alcohol consumption on LT varies according to the ALDH2 genotype, providing an insight for the accurate screening of ALD in drinkers with the rs671 GA genotype.

Machine learning for effectively avoiding overfitting is a crucial strategy for the genetic prediction of polygenic psychiatric phenotypes

The accuracy of previous genetic studies in predicting polygenic psychiatric phenotypes has been limited mainly due to the limited power in distinguishing truly susceptible variants from null variants and the resulting overfitting. A novel prediction algorithm, Smooth-Threshold Multivariate Genetic Prediction (STMGP), was applied to improve the genome-based prediction of psychiatric phenotypes by decreasing overfitting through selecting variants and building a penalized regression model. Prediction models were trained using a cohort of 3685 subjects in Miyagi prefecture and validated with an independently recruited cohort of 3048 subjects in Iwate prefecture in Japan. Genotyping was performed using HumanOmniExpressExome BeadChip Arrays. We used the target phenotype of depressive symptoms and simulated phenotypes with varying complexity and various effect-size distributions of risk alleles. The prediction accuracy and the degree of overfitting of STMGP were compared with those of state-of-the-art models (polygenic risk scores, genomic best linear-unbiased prediction, summary-data-based best linear-unbiased prediction, BayesR, and ridge regression). In the prediction of depressive symptoms, compared with the other models, STMGP showed the highest prediction accuracy with the lowest degree of overfitting, although there was no significant difference in prediction accuracy. Simulation studies suggested that STMGP has a better prediction accuracy for moderately polygenic phenotypes. Our investigations suggest the potential usefulness of STMGP for predicting polygenic psychiatric conditions while avoiding overfitting.

Genome-wide association study identifies new loci for albuminuria in the Japanese population

BACKGROUND

Urinary albumin excretion (UAE) is a risk factor for cardiovascular diseases, metabolic syndrome, chronic kidney disease, etc. Only a few genome-wide association studies (GWAS) for UAE have been conducted in the European population, but not in the Asian population. Here we conducted GWAS and identified several candidate genes harboring single nucleotide polymorphisms (SNPs) responsible for UAE in the Japanese population.

METHODS

We conducted GWAS for UAE in 7805 individuals of Asian ancestry from health-survey data collected by Tohoku Medical Megabank Organization (ToMMo) and Iwate Tohoku Medical Megabank Organization (IMM). The SNP genotype data were obtained with a SNP microarray. After imputation using a haplotype panel consisting of 2000 genome sequencing, 4,962,728 SNP markers were used for the GWAS.

RESULTS

Eighteen SNPs at 14 loci (GRM7, EXOC1/NMU, LPA, STEAP1B/RAPGEF5, SEMA3D, PRKAG2, TRIQK, SERTM1, TPT1-AS1, OR5AU1, TSHR, FMN1/RYR3, COPRS, and BRD1) were associated with UAE in the Japanese individuals. A locus with particularly strong associations was observed on TSHR, chromosome 14 [rs116622332 (p = 3.99 × 10^-10)].

CONCLUSION

In this study, we successfully identified UAE-associated variant loci in the Japanese population. Further study is required to confirm this association.

Strong association between the 12q24 locus and sweet taste preference in the Japanese population revealed by genome-wide meta-analysis

The sweet taste preference of humans is an important adaptation to ensure the acquisition of carbohydrate nutrition; however, overconsumption of sweet foods can potentially lead to diseases such as obesity and diabetes. Although previous studies have suggested that interindividual variation of human sweet taste preference is heritable, genetic loci associated with the trait have yet to be fully elucidated. Here, we genotyped 12,312 Japanese participants using the HumanCore-12+ Custom BeadChip or the HumanCore-24 Custom BeadChip microarrays. The sweet taste preference of the participants was surveyed via an internet-based questionnaire, resulting in a five-point scale of sweet taste preference. The genome-wide meta-analysis of the Japanese participants revealed a strong association between the 12q24 locus and sweet taste preference scale (P = 2.8 × 10^-70). The lead variant rs671 is monoallelic in non-East Asian populations and is located in the aldehyde dehydrogenase (ALDH2) gene, encoding an enzyme involved in alcohol metabolism. The association between the minor allele of rs671 and sweet taste preference was attenuated by adjusting for alcohol drinking. The subgroup analysis showed that the effect of rs671 on sweet taste preference was greater in males than in females. In conclusion, we found an association between the 12q24 locus and sweet taste preference in the Japanese population, and showed that the adjustment for drinking habits attenuated the association. This novel genetic association may provide new clues to elucidate mechanisms determining sweet taste preferences.

Genome-Wide Polygenic Score and the Risk of Ischemic Stroke in a Prospective Cohort

Background and Purpose—

Environmental and genetic factors contribute to the development of ischemic stroke (IS). We recently developed a genome-wide polygenic risk score (PRS) for IS using case-control datasets from 4 large-scale observational studies conducted in Japan. Our objective in the present study was to confirm the association between the PRS and the risk of IS with data from an independent prospective cohort recruited from the general Japanese population.

Methods—

A total of 3038 subjects aged ≥40 years were followed up for 10 years (2002–2012). The genome-wide PRS was calculated using genotype data from >350 000 single-nucleotide polymorphisms. The PRS levels were divided into quintiles. High and low genetic risk groups were defined as top 60% and bottom 40% of PRS, respectively. The hazard ratio (HR) for the development of IS was estimated using a Cox proportional hazards model.

Results—

During the follow-up period, 91 cases developed first-ever IS. The age- and sex-adjusted HR for IS increased with higher PRS levels ( P for trend, 0.03). Subjects with the highest quintile level of PRS had a 2.44-fold (95% CI, 1.16–5.12) greater risk for IS than those with the lowest quintile level after adjusting for age and sex. A similar association was observed after adjusting for environmental risk factors ( P for trend, 0.03). As compared with low genetic risk group, the age- and sex-adjusted HR in high genetic risk group was 1.63 (95% CI, 1.04–2.55), which was comparable to the HR of hypertension (HR, 1.41), diabetes mellitus (HR, 1.72), and smoking (HR, 1.54). The age- and sex-adjusted HR increased with the number of environmental risk factors in both high and low genetic risk groups without significant interaction.

Conclusions—

A high genome-wide PRS was a significant risk factor for IS independent of environmental risk factors in a general Japanese population. This finding suggests that PRS may be useful to identify individuals at a high risk of IS.

Study Profile of the Tohoku Medical Megabank Community-Based Cohort Study

Background

We established a community-based cohort study to assess the long-term impact of the Great East Japan Earthquake on disaster victims and gene-environment interactions on the incidence of major diseases, such as cancer and cardiovascular diseases.

Methods

We asked participants to join our cohort in the health check-up settings and assessment center based settings. Inclusion criteria were aged 20 years or over and living in Miyagi or Iwate Prefecture. We obtained information on lifestyle, effect of disaster, blood, and urine information (Type 1 survey), and some detailed measurements (Type 2 survey), such as carotid echography and calcaneal ultrasound bone mineral density. All participants agreed to measure genome information and to distribute their information widely.

Results

As a result, 87,865 gave their informed consent to join our study. Participation rate at health check-up site was about 70%. The participants in the Type 1 survey were more likely to have psychological distress than those in the Type 2 survey, and women were more likely to have psychological distress than men. Additionally, coastal residents were more likely to have higher degrees of psychological distress than inland residents, regardless of sex.

Conclusion

This cohort comprised a large sample size and it contains information on the natural disaster, genome information, and metabolome information. This cohort also had several detailed measurements. Using this cohort enabled us to clarify the long-term effect of the disaster and also to establish personalized prevention based on genome, metabolome, and other omics information.

The interaction between ABCA1 polymorphism and physical activity on the HDL-cholesterol levels in a Japanese population

Few studies have investigated the interactions between HDL-C-related SNPs identified by genome-wide association (GWA) study and physical activity (PA) on HDL-C. First, we conducted a sex-stratified GWA study in a discovery sample (2,231 men and 2,431 women) and replication sample (2,599 men and 3,109 women) to identify SNPs influencing log-transformed HDL-C in Japanese participants in the baseline survey of the Japan Multi-Institutional Collaborative Cohort Study. We also replicated previously reported HDL-C-related SNPs in a combined (discovery plus replication) sample (4,830 men and 5,540 women). We then analyzed the interactions of the HDL-C-related SNPs with PA on HDL-C. The sex-stratified GWA analyses identified 11 and 10 HDL-C-related SNPs in men and women as targets for an interaction analysis. Among these, only one interaction of ABCA1 rs1883025 with PA was statistically significant in men, after Bonferroni correction [P-interaction = 0.001 (α = 0.05/21 = 0.002)]. The per-major-allele (C allele) increase in log-transformed HDL-C was lost in men with low PA (β = 0.008) compared with those with medium (β = 0.032) or high PA (β = 0.034). These findings suggest that the benefit of carrying a C allele of ABCA1 rs1883025 on enhancing HDL-C may be attenuated in inactive men.

A Genome-Wide Association Study Identifies Blood Disorder-Related Variants Influencing Hemoglobin A1c With Implications for Glycemic Status in U.S. Hispanics/Latinos

OBJECTIVE

We aimed to identify hemoglobin A1c (HbA1c)-associated genetic variants and examine their implications for glycemic status evaluated by HbA1c in U.S. Hispanics/Latinos with diverse genetic ancestries.

RESEARCH DESIGN AND METHODS

We conducted a genome-wide association study (GWAS) of HbA1c in 9,636 U.S. Hispanics/Latinos without diabetes from the Hispanic Community Health Study/Study of Latinos, followed by a replication among 4,729 U.S. Hispanics/Latinos from three independent studies.

RESULTS

Our GWAS and replication analyses showed 10 previously known and novel loci associated with HbA1c at genome-wide significance levels (P < 5.0 × 10-8). In particular, two African ancestry-specific variants, HBB-rs334 and G6PD-rs1050828, which are causal mutations for sickle cell disease and G6PD deficiency, respectively, had ∼10 times larger effect sizes on HbA1c levels (β = -0.31% [-3.4 mmol/mol]) and -0.35% [-3.8 mmol/mol] per minor allele, respectively) compared with other HbA1c-associated variants (0.03-0.04% [0.3-0.4 mmol/mol] per allele). A novel Amerindian ancestry-specific variant, HBM-rs145546625, was associated with HbA1c and hematologic traits but not with fasting glucose. The prevalence of hyperglycemia (prediabetes and diabetes) defined using fasting glucose or oral glucose tolerance test 2-h glucose was similar between carriers of HBB-rs334 or G6PD-rs1050828 HbA1c-lowering alleles and noncarriers, whereas the prevalence of hyperglycemia defined using HbA1c was significantly lower in carriers than in noncarriers (12.2% vs. 28.4%, P < 0.001). After recalibration of the HbA1c level taking HBB-rs334 and G6PD-rs1050828 into account, the prevalence of hyperglycemia in carriers was similar to noncarriers (31.3% vs. 28.4%, P = 0.28).

CONCLUSIONS

This study in U.S. Hispanics/Latinos found several ancestry-specific alleles associated with HbA1c through erythrocyte-related rather than glycemic-related pathways. The potential influences of these nonglycemic-related variants need to be considered when the HbA1c test is performed.

GWAS of habitual coffee consumption reveals a sex difference in the genetic effect of the 12q24 locus in the Japanese population

Background

Studies on genetic effects of coffee consumption are scarce for Asian populations. We conducted a genome-wide association study (GWAS) of habitual coffee consumption in Japan using a self-reporting online survey.

Results

Candidate genetic loci associated with habitual coffee consumption were searched within a discovery cohort (N = 6,264) and confirmed in a replication cohort (N = 5,975). Two loci achieved genome-wide significance (P < 5 × 10^− 8) in a meta-analysis of the discovery and replication cohorts: an Asian population-specific 12q24 (rs79105258; P = 9.5 × 10^− 15), which harbors CUX2, and 7p21 (rs10252701; P = 1.0 × 10^− 14), in the upstream region of the aryl hydrocarbon receptor (AHR) gene, involved in caffeine metabolism. Subgroup analysis revealed a stronger genetic effect of the 12q24 locus in males (P for interaction = 8.2 × 10^− 5). Further, rs79105258 at the 12q24 locus exerted pleiotropic effects on body mass index (P = 3.5 × 10^− 4) and serum triglyceride levels (P = 8.7 × 10^− 3).

Conclusions

Our results consolidate the association of habitual coffee consumption with the 12q24 and 7p21 loci. The different effects of the 12q24 locus between males and females are a novel finding that improves our understanding of genetic influences on habitual coffee consumption.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0763-7) contains supplementary material, which is available to authorized users.

Identification of the 12q24 locus associated with fish intake frequency by genome-wide meta-analysis in Japanese populations

Background

Japan is traditionally a country with one of the highest levels of fish consumption worldwide, although the westernization of the Japanese diet has resulted in the reduction of fish consumption. A recent meta-analysis of genome-wide association studies (GWASs) on Western populations has identified a single nucleotide polymorphism (SNP) associated with fish intake frequency. Here, we examined the genetic basis for fish intake frequency among Japanese individuals.

Results

We conducted a meta-analysis of a GWAS including 12,603 Japanese individuals and identified a susceptibility locus for fish intake frequency at 12q24 (lead variant was rs11066015, P = 5.4 × 10^-11). rs11066015 was in a strong linkage disequilibrium with rs671, a well-known SNP related to alcohol metabolism. When adjusted for alcohol drinking, the association between rs11066015 and fish intake frequency was substantially attenuated. Subgroup analysis revealed that the effect of the 12q24 variant on fish intake frequency was stronger in males than in females (P for interaction = 0.007) and stronger in the older subgroup than in the younger subgroup (P for interaction = 0.006).

Conclusions

Our findings suggest that the 12q24 locus is associated with fish intake frequency via alcohol drinking. This study can help contribute to personalized nutrition information, suggesting that fish intake should be promoted to consumers who have the rs11066015 minor allele, which is genetically linked to low fish intake frequency, especially in male and older individuals.

Genomewide Association Study of Leisure-Time Exercise Behavior in Japanese Adults

Supplemental digital content is available in the text.

Purpose

Although several genetic factors may play a role in leisure-time exercise behavior, there is currently no evidence of a significant genomewide association, and candidate gene replication studies have produced inconsistent results.

Methods

We conducted a two-stage genomewide association study and candidate single-nucleotide polymorphisms (SNP) association study on leisure-time exercise behavior using 13,980 discovery samples from the Japan Multi-Institutional Collaborative Cohort (J-MICC) study, and 2036 replication samples from the Hospital-based Epidemiologic Research Program at Aichi Cancer Center-2 study. Leisure-time physical activity was measured using a self-administered questionnaire that inquired about the type, frequency and duration of exercise. Participants with ≥4 MET·h·wk⁻¹ of leisure-time physical activity were defined as exhibiting leisure-time exercise behavior. Association testing using mixed linear regression models was performed on the discovery and replication samples, after which the results were combined in a meta-analysis. In addition, we tested six candidate genetic variants derived from previous genomewide association study.

Results

We found that one novel SNP (rs10252228) located in the intergenic region between NPSR1 and DPY19L1 was significantly associated with leisure-time exercise behavior in discovery samples. This association was also significant in replication samples (combined P value by meta-analysis = 2.2 × 10⁻⁹). Several SNP linked with rs10252228 were significantly associated with gene expression of DPY19L1 and DP19L2P1 in skeletal muscle, heart, whole blood, and the nervous system. Among the candidate SNP, rs12612420 in DNAPTP6 demonstrated nominal significance in discovery samples but not in replication samples.

Conclusions

We identified a novel genetic variant associated with regular leisure-time exercise behavior. Further functional studies are required to validate the role of these variants in exercise behavior.

Mendelian Randomization Analysis Support Causal Associations of HbA1c with Circulating Triglyceride, Total and Low-density Lipoprotein Cholesterol in a Chinese Population

Previous observational studies supported a positive association of glycated hemoglobin A1c (HbA1c) level with serum triglyceride (TG), total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C). However, the causal relationship between HbA1c and either one of them was unclear in the East Asians. We performed a Mendelian Randomization (MR) analysis in a community-based study sample in Shanghai, China (n = 11,935). To clarify the cause-and-effect relationships of HbA1c with the four interested lipids, an Expanded HbA1c genetic risk score (GRS) with 17 HbA1c-related common variants and a Conservative score by excluding 11 variants were built and adopted as the Instrumental Variables (IVs), respectively. The Expanded HbA1c-GRS was associated with 0.19 unit increment in log-TG (P = 0.009), 0.42 mmol/L TC (P = 0.01), and 0.33 mmol/L LDL-C (P = 0.01); while the Conservative HbA1c-GRS was associated with 0.22 unit in log-TG (P = 0.03), 0.60 mmol/L TC (P = 0.01), and 0.51 mmol/L LDL-C (P = 0.007). No causal relationship was detected for HDL-C. Sensitivity analysis supported the above findings. In conclusions, MR analysis supports a causal role of increased HbA1c level in increment of circulating TG, TC, and LDL-C in a Chinese population.

Genome-wide analysis of polymorphism × sodium interaction effect on blood pressure identifies a novel 3'-BCL11B gene desert locus

Excessive sodium intake is a global risk factor for hypertension. Sodium effects on blood pressure vary from person to person; hence, high-risk group targeting based on personal genetic information can play a complementary role to ongoing population preventive approaches to reduce sodium consumption. To identify genetic factors that modulate sodium effects on blood pressure, we conducted a population-based genome-wide interaction analysis in 8,768 Japanese subjects, which was >3 times larger than a similar previous study. We tested 7,135,436 polymorphisms in the discovery cohort, and loci that met suggestive significance were further examined in an independent replication cohort. We found that an interaction between a novel 3'-BCL11B gene desert locus and daily sodium consumption was significantly associated with systolic blood pressure in both discovery and replication cohorts under the recessive model. Further statistical analysis of rs8022678, the sentinel variant of the 3'-BCL11B gene desert locus, showed that differences in mean systolic blood pressure between high and low sodium consumption subgroups were 5.9 mm Hg (P = 8.8 × 10-12) in rs8022678 A carriers and -0.3 mm Hg (P = 0.27) in rs8022678 A non-carriers, suggesting that the rs8022678 genotype can classify persons into sodium-sensitive (A carriers) and sodium-insensitive (A non-carriers) subgroups. Our results implied that rs8022678 A carriers may receive a greater benefit from sodium-lowering interventions than non-carriers.

Identification of Novel Candidate Markers of Type 2 Diabetes and Obesity in Russia by Exome Sequencing with a Limited Sample Size

Type 2 diabetes (T2D) and obesity are common chronic disorders with multifactorial etiology. In our study, we performed an exome sequencing analysis of 110 patients of Russian ethnicity together with a multi-perspective approach based on biologically meaningful filtering criteria to detect novel candidate variants and loci for T2D and obesity. We have identified several known single nucleotide polymorphisms (SNPs) as markers for obesity (rs11960429), T2D (rs9379084, rs1126930), and body mass index (BMI) (rs11553746, rs1956549 and rs7195386) (p < 0.05). We show that a method based on scoring of case-specific variants together with selection of protein-altering variants can allow for the interrogation of novel and known candidate markers of T2D and obesity in small samples. Using this method, we identified rs328 in LPL (p = 0.023), rs11863726 in HBQ1 (p = 8 × 10⁻⁵), rs112984085 in VAV3 (p = 4.8 × 10⁻⁴) for T2D and obesity, rs6271 in DBH (p = 0.043), rs62618693 in QSER1 (p = 0.021), rs61758785 in RAD51B (p = 1.7 × 10⁻⁴), rs34042554 in PCDHA1 (p = 1 × 10⁻⁴), and rs144183813 in PLEKHA5 (p = 1.7 × 10⁻⁴) for obesity; and rs9379084 in RREB1 (p = 0.042), rs2233984 in C6orf15 (p = 0.030), rs61737764 in ITGB6 (p = 0.035), rs17801742 in COL2A1 (p = 8.5 × 10⁻⁵), and rs685523 in ADAMTS13 (p = 1 × 10⁻⁶) for T2D as important susceptibility loci in Russian population. Our results demonstrate the effectiveness of whole exome sequencing (WES) technologies for searching for novel markers of multifactorial diseases in cohorts of limited size in poorly studied populations.