The allele frequency of ALDH2*Glu504Lys and ADH1B*Arg47His for the Ryukyu islanders and their history of expansion among East Asians

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Exploring the genetic diversity of the Japanese population: Insights from a large-scale whole genome sequencing analysis

The Japanese archipelago is a terminal location for human migration, and the contemporary Japanese people represent a unique population whose genomic diversity has been shaped by multiple migrations from Eurasia. We analyzed the genomic characteristics that define the genetic makeup of the modern Japanese population from a population genetics perspective from the genomic data of 9,287 samples obtained by high-coverage whole-genome sequencing (WGS) by the National Center Biobank Network. The dataset comprised populations from the Ryukyu Islands and other parts of the Japanese archipelago (Hondo). The Hondo population underwent two episodes of population decline during the Jomon period, corresponding to the Late Neolithic, and the Edo period, corresponding to the Early Modern era, while the Ryukyu population experienced a population decline during the shell midden period of the Late Neolithic in this region. Haplotype analysis suggested increased allele frequencies for genes related to alcohol and fatty acid metabolism, which were reported as loci that had experienced positive natural selection. Two genes related to alcohol metabolism were found to be 12,500 years out of phase with the time when they began to increase in the allele frequency; this finding indicates that the genomic diversity of Japanese people has been shaped by events closely related to agriculture and food production.

Demographic history of Ryukyu islanders at the southern part of the Japanese Archipelago inferred from whole-genome resequencing data

The Ryukyu Islands are located in the southernmost part of the Japanese Archipelago and consist of several island groups. Each island group has its own history and culture, which differ from those of mainland Japan. People of the Ryukyu Islands are genetically subdivided; however, their detailed demographic history remains unclear. We report the results of a whole-genome sequencing analysis of a total of 50 Ryukyu islanders, focusing on genetic differentiation between Miyako and Okinawa islanders. We confirmed that Miyako and Okinawa islanders cluster differently in principal component analysis and ADMIXTURE analysis and that there is a population structure among Miyako islanders. The present study supports the hypothesis that population differentiation is primarily caused by genetic drift rather than by differences in the rate of migration from surrounding regions, such as the Japanese main islands or Taiwan. In addition, the genetic cline observed among Miyako and Okinawa islanders can be explained by recurrent migration beyond the bounds of these islands. Our analysis also suggested that the presence of multiple subpopulations during the Neolithic Ryukyu Jomon period is not crucial to explain the modern Ryukyu populations. However, the assumption of multiple subpopulations during the time of admixture with mainland Japanese is necessary to explain the modern Ryukyu populations. Our findings add insights that could help clarify the complex history of populations in the Ryukyu Islands.

Association of drinking behaviors with TXNIP DNA methylation levels in leukocytes among the general Japanese population

Background: Thioredoxin-interacting protein (TXNIP) controls the cellular redox balance by binding to and inhibiting the expression and function of thioredoxin. DNA methylation of the TXNIP gene is involved in the regulation of TXNIP mRNA expression. Changes in TXNIP DNA methylation levels are associated with the development of various diseases such as type 2 diabetes mellitus (T2DM). However, few studies have focused on the influence of lifestyle factors such as alcohol intake on TXNIP DNA methylation.Objectives: This research examines the association of drinking behaviors with TXNIP DNA methylation levels in the general Japanese population.Methods: We conducted a cross-sectional study of 404 subjects (176 males and 228 females) who were divided into non-, moderate and heavy drinkers based on self-reported drinking behaviors. TXNIP DNA methylation levels in leukocytes were determined using a pyrosequencing assay.Results: The mean TXNIP DNA methylation level in heavy drinkers (74.2%) was significantly lower than that in non- and moderate drinkers (non: 77.7%, p < .001; moderate: 76.6%, p = .011). Multivariable linear regression analysis showed that log-transformed values of daily (b = -1.34; p < .001) and cumulative (b = -1.06; p = .001) alcohol consumption were associated with decreased TXNIP DNA methylation levels.Conclusion: TXNIP DNA methylation levels in heavy drinkers was lower than in non- and- moderate drinkers. Decreased TXNIP DNA methylation level increases the expression of TXNIP and elevates the risk of developing of diseases such as T2DM. Therefore, decreasing alcohol use in heavy drinkers may lessen the likelihood of some alcohol-related illnesses moderated through TXNIP DNA methylation.

Whole-Genome Sequencing of a 900-Year-Old Human Skeleton Supports Two Past Migration Events from the Russian Far East to Northern Japan

Recent studies on paleogenomics have reported some Paleolithic and Neolithic genomes that have provided new insights into the human population history in East and Northeast Asia. However, there remain some cases where more recent migration events need to be examined to elucidate the detailed formation process of local populations. Although the area around northern Japan is one of the regions archaeologically suggested to have been affected by migration waves after the Neolithic period, the genetic source of these migrations are still unclear. Thus, genomic data from such past migrant populations would be highly informative to clarify the detailed formation process of local populations in this region. Here, we report the genome sequence of a 900-year-old adult female (NAT002) belonging to the prehistoric Okhotsk people, who have been considered to be the past migrants to northern Japan after the Neolithic period. We found a close relationship between NAT002 and modern Lower Amur populations and past admixture events between the Amur, Jomon, and Kamchatka ancestries. The admixture dating suggested migration of Amur-related ancestry at approximately 1,600 BP, which is compatible with the archaeological evidence regarding the settlement of the Okhotsk people. Our results also imply migration of Kamchatka-related ancestry at approximately 2,000 BP. In addition, human leukocyte antigen (HLA) typing detected the HLA-B*40 allele, which is reported to increase the risk of arthritis, suggesting the genetic vulnerability of NAT002 to hyperostosis, which was observed around her chest clavicle.

Genome-Wide Natural Selection Signatures Are Linked to Genetic Risk of Modern Phenotypes in the Japanese Population

Elucidation of natural selection signatures and relationships with phenotype spectra is important to understand adaptive evolution of modern humans. Here, we conducted a genome-wide scan of selection signatures of the Japanese population by estimating locus-specific time to the most recent common ancestor using the ascertained sequentially Markovian coalescent (ASMC), from the biobank-based large-scale genome-wide association study data of 170,882 subjects. We identified 29 genetic loci with selection signatures satisfying the genome-wide significance. The signatures were most evident at the alcohol dehydrogenase (ADH) gene cluster locus at 4q23 (P_ASMC = 2.2 × 10⁻³⁶), followed by relatively strong selection at the FAM96A (15q22), MYOF (10q23), 13q21, GRIA2 (4q32), and ASAP2 (2p25) loci (P_ASMC < 1.0 × 10⁻¹⁰). The additional analysis interrogating extended haplotypes (integrated haplotype score) showed robust concordance of the detected signatures, contributing to fine-mapping of the genes, and provided allelic directional insights into selection pressure (e.g., positive selection for ADH1B-Arg48His and HLA-DPB1*04:01). The phenome-wide selection enrichment analysis with the trait-associated variants identified a variety of the modern human phenotypes involved in the adaptation of Japanese. We observed population-specific evidence of enrichment with the alcohol-related phenotypes, anthropometric and biochemical clinical measurements, and immune-related diseases, differently from the findings in Europeans using the UK Biobank resource. Our study demonstrated population-specific features of the selection signatures in Japanese, highlighting a value of the natural selection study using the nation-wide biobank-scale genome and phenotype data.

A genome-wide association study on confection consumption in a Japanese population: the Japan Multi-Institutional Collaborative Cohort Study

Differences in individual eating habits may be influenced by genetic factors, in addition to cultural, social or environmental factors. Previous studies suggested that genetic variants within sweet taste receptor genes family were associated with sweet taste perception and the intake of sweet foods. The aim of this study was to conduct a genome-wide association study (GWAS) to find genetic variations that affect confection consumption in a Japanese population. We analysed GWAS data on confection consumption using 14 073 participants from the Japan Multi-Institutional Collaborative Cohort study. We used a semi-quantitative FFQ to estimate food intake that was validated previously. Association of the imputed variants with confection consumption was performed by linear regression analysis with adjustments for age, sex, total energy intake and principal component analysis components 1-3. Furthermore, the analysis was repeated adjusting for alcohol intake (g/d) in addition to the above-described variables. We found 418 SNP located in 12q24 that were associated with confection consumption. SNP with the ten lowest P-values were located on nine genes including at the BRAP, ACAD10 and aldehyde dehydrogenase 2 regions on 12q24.12-13. After adjustment for alcohol intake, no variant was associated with confections intake with genome-wide significance. In conclusion, we found a significant number of SNP located on 12q24 genes that were associated with confections intake before adjustment for alcohol intake. However, all of them lost statistical significance after adjustment for alcohol intake.

Evolutionary history of disease-susceptibility loci identified in longitudinal exome-wide association studies

Background

Our longitudinal exome‐wide association studies previously detected various genetic determinants of complex disorders using ~26,000 single‐nucleotide polymorphisms (SNPs) that passed quality control and longitudinal medical examination data (mean follow‐up period, 5 years) in 4884–6022 Japanese subjects. We found that allele frequencies of several identified SNPs were remarkably different among four ethnic groups. Elucidating the evolutionary history of disease‐susceptibility loci may help us uncover the pathogenesis of the related complex disorders.

Methods

In the present study, we conducted evolutionary analyses such as extended haplotype homozygosity, focusing on genomic regions containing disease‐susceptibility loci and based on genotyping data of our previous studies and datasets from the 1000 Genomes Project.

Results

Our evolutionary analyses suggest that derived alleles of rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs closely located at 12q24.1 associated with type 2 diabetes mellitus, obesity, dyslipidemia, and three complex disorders (hypertension, hyperuricemia, and dyslipidemia), respectively, rapidly expanded after the human dispersion from Africa (Out‐of‐Africa). Allele frequencies of GGA3 and six SNPs at 12q24.1 appeared to have remarkably changed in East Asians, whereas the derived alleles of rs34902660 of SLC17A3 and rs7656604 at 4q13.3 might have spread across Japanese and non‐Africans, respectively, although we cannot completely exclude the possibility that allele frequencies of disease‐associated loci may be affected by demographic events.

Conclusion

Our findings indicate that derived allele frequencies of nine disease‐associated SNPs (rs78338345 of GGA3, rs7656604 at 4q13.3, rs34902660 of SLC17A3, and six SNPs at 12q24.1) identified in the longitudinal exome‐wide association studies largely increased in non‐Africans after Out‐of‐Africa.

Biomedical and genetic characteristics of the Ryukyuans: demographic history, diseases and physical and physiological traits

Context: The Ryukyu Islands stretch across a southwestern area of the Japanese Archipelago. Because of their unique geographical and historical backgrounds, Ryukyuans have their own genetic and phenotypic characteristics, which have been disclosed in previous anthropological and biomedical studies. Objective: The history, peopling and biomedical and genetic characteristics of Ryukyuans are reviewed and future research directions are discussed. Conclusion: Morphological and genetic studies have suggested the complex demographic history of Ryukyuans and their relationships with other Asian populations. Knowledge of population formation processes is important to understand the distribution of pathogens. In viral infectious diseases, some strains that may be associated with disease symptoms are specific to Ryukyuans. Dramatic changes in diet have played an important role among Ryukyuans in terms of increases in lifestyle-related diseases and mortality risks. To achieve a better understanding of pathogenic disease factors, further integration of findings regarding the genetic and biomedical characteristics of the Ryukyuans is needed.

Deep whole-genome sequencing reveals recent selection signatures linked to evolution and disease risk of Japanese

Understanding natural selection is crucial to unveiling evolution of modern humans. Here, we report natural selection signatures in the Japanese population using 2234 high-depth whole-genome sequence (WGS) data (25.9×). Using rare singletons, we identify signals of very recent selection for the past 2000–3000 years in multiple loci (ADH cluster, MHC region, BRAP-ALDH2, SERHL2). In large-scale genome-wide association study (GWAS) dataset (n = 171,176), variants with selection signatures show enrichment in heterogeneity of derived allele frequency spectra among the geographic regions of Japan, highlighted by two major regional clusters (Hondo and Ryukyu). While the selection signatures do not show enrichment in archaic hominin-derived genome sequences, they overlap with the SNPs associated with the modern human traits. The strongest overlaps are observed for the alcohol or nutrition metabolism-related traits. Our study illustrates the value of high-depth WGS to understand evolution and their relationship with disease risk.

Recent natural selection left signals in human genomes. Here, Okada et al. generate high-depth whole-genome sequence (WGS) data (25.9×) from 2,234 Japanese people of the BioBank Japan Project (BBJ), and identify signals of recent natural selection which overlap variants associated with human traits.

[Importance of an Aldehyde Dehydrogenase 2 Polymorphism in Preventive Medicine]

Unlike genetic alterations in other aldehyde dehydrogenase (ALDH) isozymes, a defective ALDH2 polymorphism (rs671), which is carried by almost half of East Asians, does not show a clear phenotype such as a shortened life span. However, impacts of a defective ALDH2 allele, ALDH2*2, on various disease risks have been reported. As ALDH2 is responsible for the detoxification of endogenous aldehydes, a negative effect of this polymorphism is predicted, but bidirectional effects have been actually observed and the mechanisms underlying such influences are often complex. One reason for this complexity may be the existence of compensatory aldehyde detoxification systems and the secondary effects of these systems. There are many issues to be addressed with regard to the ALDH2 polymorphism in the field of preventive medicine, including the following concerns. First, ALDH2 in the fetal stage plays a role in aldehyde detoxification; therefore, prenatal health effects of environmental aldehyde exposure are of concern for ALDH2*2-carrying fetuses. Second, ALDH2*2 carriers are at high risk of drinking-related cancers. However, their drinking habits result in less worsening of physiological findings, such as energy metabolism index and liver functions, compared with non-ALDH2*2 carriers, and therefore opportunities to detect excessive drinking can be lost. Third, personalized medicine such as personalized prescriptions for ALDH2*2 carriers will be required in the clinical setting, and accumulation of evidence is awaited. Lastly, since the ALDH2 polymorphism is not considered in workers' limits of exposure to aldehydes and their precursors, efforts to lower exposure levels beyond legal standards are required.

Analysis of mainland Japanese and Okinawan Japanese populations using the precision ID Ancestry Panel

We typed 165 AIMs in 49 mainland Japanese and 47 Okinawa Japanese using the Precision ID Ancestry Panel (Thermo Fisher Scientific). None of the 165 SNPs showed significant deviation from Hardy-Weinberg equilibrium in the mainland Japanese. One SNP (rs3943253) showed significant deviation from Hardy-Weinberg equilibrium in Okinawa Japanese. Fisher's exact tests showed that the genotype frequencies of 14 loci were significantly different (p<0.05) between the two populations before correction for multiple testing. After Bonferroni correction, only rs671 remained statistically significant (p<0.0003). This SNP is located in the ALDH2 gene. The mutant A allele is associated with increased side effects after alcohol intake. The frequency of the GG genotype (wild type) was higher in the Okinawa Japanese (78.7%) than in mainland Japanese (34.7%; Bonferroni corrected P<0.001). For 31 (63.3%) of the mainland Japanese and 42 (89.4%) of Okinawa Japanese, the highest population likelihood was obtained with the Japanese reference population. However, only in a few individuals, the likelihoods were significantly different from those calculated using reference data from neighboring populations. The likelihoods for mainland Japanese and Okinawa Japanese were not significantly different from each other for any of the investigated individuals. STRUCTURE and PCA analyses showed that mainland Japanese, Okinawa Japanese, and East Asians could not be differentiated with the Precision ID Ancestry Panel.