Study Profile of the Tsuruoka Metabolomics Cohort Study (TMCS)

The Tsuruoka Metabolomics Cohort Study (TMCS) is an ongoing population-based cohort study being conducted in the rural area of Yamagata Prefecture, Japan. This study aimed to enhance the precision prevention of multi-factorial, complex diseases, including non-communicable and aging-associated diseases, by improving risk stratification and prediction measures. At baseline, 11,002 participants aged 35-74 years were recruited in Tsuruoka City, Yamagata Prefecture, Japan, between 2012 and 2015, with an ongoing follow-up survey. Participants underwent various measurements, examinations, tests, and questionnaires on their health, lifestyle, and social factors. This study used an integrative approach with deep molecular profiling to identify potential biomarkers linked to phenotypes that underpin disease pathophysiology and provide better mechanistic insights into social health determinants. The TMCS incorporates multi-omics data, including genetic and metabolomic analyses of 10,933 participants and comprehensive data collection ranging from physical, psychological, behavioral, and social to biological data. The metabolome is used as a phenotypic probe because it is sensitive to changes in physiological and external conditions. The TMCS focuses on collecting outcomes for cardiovascular disease, cancer incidence and mortality, disability, functional decline due to aging and disease sequelae, and the variation in health status within the body represented by omics analysis that lies between exposure and disease. It contains several sub-studies on aging, heated tobacco products, and women's health. This study is notable for its robust design, high participation rate (89%), and long-term repeated surveys. Moreover, it contributes to precision prevention in Japan and East Asia as a well-established multi-omics platform.

East Asian-specific and cross-ancestry genome-wide meta-analyses provide mechanistic insights into peptic ulcer disease

Peptic ulcer disease (PUD) refers to acid-induced injury of the digestive tract, occurring mainly in the stomach (gastric ulcer (GU)) or duodenum (duodenal ulcer (DU)). In the present study, we conducted a large-scale, cross-ancestry meta-analysis of PUD combining genome-wide association studies with Japanese and European studies (52,032 cases and 905,344 controls), and discovered 25 new loci highly concordant across ancestries. An examination of GU and DU genetic architecture demonstrated that GUs shared the same risk loci as DUs, although with smaller genetic effect sizes and higher polygenicity than DUs, indicating higher heterogeneity of GUs. Helicobacter pylori (HP)-stratified analysis found an HP-related host genetic locus. Integrative analyses using bulk and single-cell transcriptome profiles highlighted the genetic factors of PUD being enriched in the highly expressed genes in stomach tissues, especially in somatostatin-producing D cells. Our results provide genetic evidence that gastrointestinal cell differentiations and hormone regulations are critical in PUD etiology.

Integrated analysis of human DNA methylation, gene expression, and genomic variation in iMETHYL database using kernel tensor decomposition-based unsupervised feature extraction

Integrating gene expression, DNA methylation, and genomic variants simultaneously without location coincidence (i.e., irrespective of distance from each other) or pairwise coincidence (i.e., direct identification of triplets of gene expression, DNA methylation, and genomic variants, and not integration of pairwise coincidences) is difficult. In this study, we integrated gene expression, DNA methylation, and genome variants from the iMETHYL database using the recently proposed kernel tensor decomposition-based unsupervised feature extraction method with limited computational resources (i.e., short CPU time and small memory requirements). Our methods do not require prior knowledge of the subjects because they are fully unsupervised in that unsupervised tensor decomposition is used. The selected genes and genomic variants were significantly targeted by transcription factors that were biologically enriched in KEGG pathway terms as well as in the intra-related regulatory network. The proposed method is promising for integrated analyses of gene expression, methylation, and genomic variants with limited computational resources.

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.

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.

Evaluation of short-term epigenetic age fluctuation

Considerable effort has been spent on lowering and maintaining the epigenetic age. However, the extent to which epigenetic age fluctuates under normal conditions is poorly understood. Therefore, we analyzed methylation data from monocytes and peripheral blood mononuclear cells collected from two Japanese men. The ranges of the Pan-tissue, Skin and blood, and DNAm PhenoAge epigenetic age during 3 months were ≥ 5.62, ≥ 3.04, and ≥ 8.23 years, and the maximum daily changes were 5.21, 3.20, and 6.53 years, respectively. These fluctuations were not suppressed by correcting for cell-type composition. Although the underlying biological mechanism remains unclear, there was a nonnegligible degree of age fluctuation which should inform personalized clinical applications.

Low MICA Gene Expression Confers an Increased Risk of Graves' Disease: A Mendelian Randomization Study

Background: Expression of natural killer group 2 member D (NKG2D) ligand (NKG2DL) plays a major role as a "danger signal" on stressed cells to promote removal of the latter by NKG2D-expressing cytotoxic lymphocytes. NKG2DL expression has been found in peripheral immune cells as well, such as in macrophages; however, the effect of this expression is yet to be determined. Methods: We determined instrumental variables (IVs; R² <0.01 in linkage disequilibrium), explaining the major variance in major histocompatibility complex class I chain-related protein A (MICA) and B (MICB) gene expression levels from the expression-quantitative trait locus (eQTL) of NKG2DLs based on the RNA-seq analysis of peripheral blood mononuclear cells (PBMCs) from 381 Japanese. Simultaneously, the target outcomes were filtered by PheWAS from 58 health risks, using a community-based cohort study composed of 44,739 Japanese residents. Finally, we estimated the causal effect of gene expression levels on the outcomes using the Mendelian randomization approach. Results: We determined nine and four IVs, explaining 87.6% and 33.0% of MICA and MICB gene expression levels, respectively. In the association test, we identified 10 or 13 significant outcomes associated with the MICA or MICB eQTLs, respectively, as well as the causal effect of MICA expression on Graves' disease (GD) (p = 4.2 × 10^-3; odds ratio per 1 S.D. difference in the expression: 0.983 [confidence interval: 0.971-0.995]), using the weighted median estimator, without significant pleiotropy (p > 0.05), and the results were consistent across the sensitivity analyses. Conclusions: Our study provide novel evidence associating NKG2DL expression with GD, an autoimmune thyroiditis; direction of the effect indicated the immunoregulatory role of MICA expression in PBMCs, suggesting the importance of further functional assays in inflammatory diseases.

Longitudinal DNA methylation dynamics as a practical indicator in clinical epigenetics

BACKGROUND

One of the fundamental assumptions of DNA methylation in clinical epigenetics is that DNA methylation status can change over time with or without interplay with environmental and clinical conditions. However, little is known about how DNA methylation status changes over time under ordinary environmental and clinical conditions. In this study, we revisited the high frequency longitudinal DNA methylation data of two Japanese males (24 time-points within three months) and characterized the longitudinal dynamics.

RESULTS

The results showed that the majority of CpGs on Illumina HumanMethylation450 BeadChip probe set were longitudinally stable over the time period of three months. Focusing on dynamic and stable CpGs extracted from datasets, dynamic CpGs were more likely to be reported as epigenome-wide association study (EWAS) markers of various traits, especially those of immune- and inflammatory-related traits; meanwhile, the stable CpGs were enriched in metabolism-related genes and were less likely to be EWAS markers, indicating that the stable CpGs are stable both in the short-term within individuals and under various environmental and clinical conditions.

CONCLUSIONS

This study indicates that CpGs with different stabilities are involved in different functions and traits, and thus, they are potential indicators that can be applied for clinical epigenetic studies to outline underlying mechanisms.

Characterization of a novel nicotinamide adenine dinucleotide-cytochrome b5 reductase mutation associated with canine hereditary methemoglobinemia

Hereditary methemoglobinemia associated with nicotinamide adenine dinucleotide-cytochrome b5 reductase (b5R) deficiency is a rare autosomal recessive disorder in animals. Recently, nonsynonymous b5R gene (CYB5R3) variants have been reported to be associated with canine and feline hereditary methemoglobinemia. However, the underlying molecular mechanisms of canine and feline methemoglobinemia caused by these nonsynonymous variants have not yet been reported. Previously, we reported a Pomeranian dog family with hereditary methemoglobinemia, carrying CYB5R3 mutation of an A>C transition at codon 194 in exon 7, replacing an isoleucine residue with leucine (p.Ile194Leu). In this study, we investigated the enzymatic and structural properties of the soluble form of wild-type and Ile194Leu canine b5Rs to characterize the effects of this missense mutation. Our results showed that the kinetic properties of the mutant enzyme were not affected by this amino acid substitution. The secondary structure of the wild-type and Ile194Leu b5Rs detected by circular dichroism showed a similar pattern. However, the mutant enzyme exhibited decreased heat stability and increased susceptibility to trypsin hydrolysis. Moreover, the thermostability and unfolding measurements indicated that the mutant enzyme was more sensitive to temperature-dependent denaturation than the wild-type b5R. We concluded from these results that unstable mutant enzyme properties with normal enzymatic activity would be associated with hereditary methemoglobinemia in the Pomeranian dog family.

Retrospective study on intercurrent pancreatitis with Babesia gibsoni infection in dogs

Babesia rossi infection has been reported to be associated with the high prevalence of pancreatitis in dogs. In this study, we retrospectively investigated whether pancreatitis occurs in B. gibsoni-infected dogs. The clinical manifestations, and hematological and serum biochemical examination results, including canine pancreatic-specific lipase (cPL), in 20 B. gibsoni-infected dogs were analyzed. The cPL concentration exceeded 400 µg/l in only 2 dogs, and they were suspected of having pancreatitis. Although the cPL concentration did not correlate with the degree of anemia or the level of parasitemia, it correlated with the band neutrophil count, platelet count, and blood urea nitrogen (BUN) level. Our study suggested that the prevalence of pancreatitis is lower among B. gibsoni-infected dogs than B. rossi-infected dogs.

Familial Congenital Methemoglobinemia in Pomeranian Dogs Caused by a Missense Variant in the NADH-Cytochrome B5 Reductase Gene

Background

In veterinary medicine, congenital methemoglobinemia associated with nicotinamide adenine dinucleotide (NADH)‐cytochrome b5 reductase (b5R) deficiency is rare. It has been reported in several breeds of dogs, but little information is available about its etiology.

Objectives

To analyze the NADH‐cytochrome b5 reductase gene, CYB5R3, in a Pomeranian dog family with methemoglobinemia suspected to be caused by congenital b5R deficiency.

Animals

Three Pomeranian dogs from a family with methemoglobinemia were analyzed. Five healthy beagles and 5 nonrelated Pomeranian dogs without methemoglobinemia were used as controls.

Methods

Methemoglobin concentration, b5R activity, and reduced glutathione (GSH) concentration were measured, and a turbidity index was used to evaluate Heinz body formation. The CYB5R3 genes of the affected dog and healthy dogs were analyzed by direct sequencing.

Results

Methemoglobin concentrations in erythrocytes of the affected dogs were remarkably higher than those of the control dogs. The b5R activity of the affected dogs was notably lower than that of the control dogs. DNA sequencing indicated that this Pomeranian family carried a CYB5R3 gene missense variant (ATC→CTC at codon 194) that resulted in the replacement of isoleucine (Ile) by leucine (Leu).

Conclusions and Clinical Importance

This dog family had familial congenital methemoglobinemia caused by b5R deficiency, which resulted from a nonsynonymous variant in the CYB5R3 gene. This variation (c.580A>C) led to an amino acid substitution (p.Ile194Leu), and Ile194 was located in the proximal region of the NADH‐binding motif. Our data suggested that this variant in the canine CYB5R3 gene would affect function of the b5R in erythrocytes.