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.

Natural Selection Signatures in the Hondo and Ryukyu Japanese Subpopulations

Natural selection signatures across Japanese subpopulations are under-explored. Here we conducted genome-wide selection scans with 622,926 single nucleotide polymorphisms for 20,366 Japanese individuals, who were recruited from the main-islands of Japanese Archipelago (Hondo) and the Ryukyu Archipelago (Ryukyu), representing two major Japanese subpopulations. The integrated haplotype score (iHS) analysis identified several signals in one or both subpopulations. We found a novel candidate locus at IKZF2, especially in Ryukyu. Significant signals were observed in the major histocompatibility complex region in both subpopulations. The lead variants differed and demonstrated substantial allele frequency differences between Hondo and Ryukyu. The lead variant in Hondo tags HLA-A*33:03-C*14:03-B*44:03-DRB1*13:02-DQB1*06:04-DPB1*04:01, a haplotype specific to Japanese and Korean. While in Ryukyu, the lead variant tags DRB1*15:01-DQB1*06:02, which had been recognized as a genetic risk factor for narcolepsy. In contrast, it is reported to confer protective effects against type 1 diabetes and human T lymphotropic virus type 1-associated myelopathy/tropical spastic paraparesis. The FastSMC analysis identified 8 loci potentially affected by selection within the past 20-150 generations, including 2 novel candidate loci. The analysis also showed differences in selection patterns of ALDH2 between Hondo and Ryukyu, a gene recognized to be specifically targeted by selection in East Asian. In summary, our study provided insights into the selection signatures within the Japanese and nominated potential sources of selection pressure.

Genomic insights into a tripartite ancestry in the Southern Ryukyu Islands

A tripartite structure for the genetic origin of Japanese populations states that present-day populations are descended from three main ancestors: (1) the indigenous Jomon hunter-gatherers; (2) a Northeast Asian component that arrived during the agrarian Yayoi period; and (3) a major influx of East Asian ancestry in the imperial Kofun period. However, the genetic heterogeneity observed in different regions of the Japanese archipelago highlights the need to assess the applicability and suitability of this model. Here, we analyse historic genomes from the southern Ryukyu Islands, which have unique cultural and historical backgrounds compared with other parts of Japan. Our analysis supports the tripartite structure as the best fit in this region, with significantly higher estimated proportions of Jomon ancestry than mainland Japanese. Unlike the main islands, where each continental ancestor was directly brought by immigrants from the continent, those who already possessed the tripartite ancestor migrated to the southern Ryukyu Islands and admixed with the prehistoric people around the eleventh century AD, coinciding with the emergence of the Gusuku period. These results reaffirm the tripartite model in the southernmost extremes of the Japanese archipelago and show variability in how the structure emerged in diverse geographic regions.

Factors associated with bone thickness: Comparison of the cranium and humerus

Cortical bone thickness is important for the mechanical function of bone. Ontogeny, aging, sex, body size, hormone levels, diet, behavior, and genetics potentially cause variations in postcranial cortical robusticity. However, the factors associated with cranial cortical robusticity remain poorly understood. Few studies have examined cortical robusticity in both cranial and postcranial bones jointly. In the present study, we used computed tomography (CT) images to measure cortical bone thicknesses in the cranial vault and humeral diaphysis. This study clearly showed that females have a greater cranial vault thickness and greater age-related increase in cranial vault thickness than males. We found an age-related increase in the full thickness of the temporal cranial vault and the width of the humeral diaphysis, as well as an age-related decrease in the cortical thickness of the frontal cranial vault and the cortical thickness of the humeral diaphysis, suggesting that the mechanisms of bone modeling in cranial and long bones are similar. A positive correlation between cortical indices in the cranial vault and humeral diaphysis also suggested that common factors affect cortical robusticity. We also examined the association of polymorphisms in the WNT16 and TNFSF11 genes with bone thickness. However, no significant associations were observed. The present study provides fundamental knowledge about similarities and differences in the mechanisms of bone modeling between cranial and postcranial bones.

Cultural transmission of traditional songs in the Ryukyu Archipelago

Geographic patterns of cultural variations are affected by how cultural traits are transmitted within and between populations. It has been argued that cultural traits are transmitted in different manners depending on their characteristics; for example, words for basic concepts are less liable to horizontal transmission between populations (i.e., borrowing) than other words. Here we examine the geographic variation of traditional songs in the Ryukyu Archipelago, southwestern islands of Japan, to explore cultural evolution of music with a focus on different social contexts in which songs are sung. Published scores of 1,342 traditional songs are coded using the CantoCore song classification scheme and distances between the songs are calculated from the codings. Neighbor-Net graphs of regions/islands are generated on the basis of the musical distances, and delta scores are obtained to examine the treelikeness of the networks. We also perform analysis of molecular variance (AMOVA) to evaluate the extent of musical diversification among regions/islands. Our results suggest that horizontal transmission between populations has played a greater role in the formation of musical diversity than that of linguistic diversity in the Ryukyu Archipelago and that the social context in which songs are sung has an effect on how they are transmitted within and between populations. In addition, we compare the observed patterns of song diversity among regions/islands with those of lexical and mitochondrial-DNA (mtDNA) diversity, showing that the variation of songs sung in the "work" context are associated with the linguistic variation, whereas no association is found between the musical and genetic variation.

Genetic Alterations in Adult T-Cell Leukemia/Lymphoma: Novel Discoveries with Clinical and Biological Significance

Adult T-cell leukemia/lymphoma (ATLL) is a refractory T-cell neoplasm that develops in human T-cell leukemia virus type-I (HTLV-1) carriers. Large-scale comprehensive genomic analyses have uncovered the landscape of genomic alterations of ATLL and have identified several altered genes related to prognosis. The genetic alterations in ATLL are extremely enriched in the T-cell receptor/nuclear factor-κB pathway, suggesting a pivotal role of deregulation in this pathway in the transformation of HTLV-1-infected cells. Recent studies have revealed the process of transformation of HTLV-1-infected cells by analyzing longitudinal samples from HTLV-1 carriers and patients with overt ATLL, an endeavor that might enable earlier ATLL diagnosis. The latest whole-genome sequencing study discovered 11 novel alterations, including CIC long isoform, which had been overlooked in previous studies employing exome sequencing. Our study group performed the targeted sequencing of ATLL in Okinawa, the southernmost island in Japan and an endemic area of HTLV-1, where the comprehensive genetic alterations had never been analyzed. We found associations of genetic alterations with HTLV-1 strains phylogenetically classified based on the tax gene, an etiological virus factor in ATLL. This review summarizes the genetic alterations in ATLL, with a focus on their clinical significance, geographical heterogeneity, and association with HTLV-1 strains.

Epidemiology of developmental dysplasia of the hip: analysis of Japanese national database

BackgroundDevelopmental dysplasia of the hip (DDH) is a cluster of hip development disorders that affects infants. The incidence of DDH-related dislocation (DDH-dislocation) is reportedly 0.1-0.3%; however, the nationwide incidence of DDH-dislocation in Japan has not been previously reported. The primary aim of this study was to report the nationwide incidence of DDH-dislocation in Japan using the National Database of Health Insurance Claims and Specific Health Checkups of Japan (NDB), and to examine its regional variation across Japan.MethodsThis was a retrospective birth cohort study using the NDB. Data on patients born between 2011 and 2013 and assigned DDH-dislocation-related disease codes during 2011-2018 were extracted. Among these, patients who underwent treatment for DDH-dislocation between 2011 and 2018 were defined as patients with DDH-dislocation.ResultsAcross the 2011, 2012, and 2013 birth cohorts, 2367 patients were diagnosed with DDH-dislocation, yielding the nationwide incidence of 0.076%. Region-specific incidence rates were almost similar across Japan. Secondary analyses revealed that 273 (11.5%) patients were diagnosed at the age of ≥1 year. The effect of birth during the cold months on the incidence of DDH-dislocation was significant (relative risk [RR]=1.89, 95% confidence interval [CI]: 1.75-2.06). The risk of DDH-dislocation among girls was approximately seven times higher than that among boys.ConclusionsThis is the first study to report the nationwide incidence of DDH-dislocation in Japan, which was estimated at 0.076%. The regional variation was trivial and unlikely to be clinically significant. Thus, the incidence rates were approximately equal across all regions in Japan.

Fine-Scale Genetic Structure and Demographic History in the Miyako Islands of the Ryukyu Archipelago

The Ryukyu Archipelago is located in the southwest of the Japanese islands and is composed of dozens of islands, grouped into the Miyako Islands, Yaeyama Islands, and Okinawa Islands. Based on the results of principal component analysis on genome-wide single-nucleotide polymorphisms, genetic differentiation was observed among the island groups of the Ryukyu Archipelago. However, a detailed population structure analysis of the Ryukyu Archipelago has not yet been completed. We obtained genomic DNA samples from 1,240 individuals living in the Miyako Islands, and we genotyped 665,326 single-nucleotide polymorphisms to infer population history within the Miyako Islands, including Miyakojima, Irabu, and Ikema islands. The haplotype-based analysis showed that populations in the Miyako Islands were divided into three subpopulations located on Miyakojima northeast, Miyakojima southwest, and Irabu/Ikema. The results of haplotype sharing and the D statistics analyses showed that the Irabu/Ikema subpopulation received gene flows different from those of the Miyakojima subpopulations, which may be related with the historically attested immigration during the Gusuku period (900 − 500 BP). A coalescent-based demographic inference suggests that the Irabu/Ikema population firstly split away from the ancestral Ryukyu population about 41 generations ago, followed by a split of the Miyako southwest population from the ancestral Ryukyu population (about 16 generations ago), and the differentiation of the ancestral Ryukyu population into two populations (Miyako northeast and Okinawajima populations) about seven generations ago. Such genetic information is useful for explaining the population history of modern Miyako people and must be taken into account when performing disease association studies.

An analysis of the demographic history of the risk allele R4810K in RNF213 of moyamoya disease

BACKGROUND

Ring finger protein 213 (RNF213) is a susceptibility gene of moyamoya disease (MMD). A previous case-control study and a family analysis demonstrated a strong association of the East Asian-specific variant, R4810K (rs112735431), with MMD. Our aim is to uncover evolutionary history of R4810K in East Asian populations.

METHODS

The RNF213 locus of 24 MMD patients in Japan were sequenced using targeted-capture sequencing. Based on the sequence data, we conducted population genetic analysis and estimated the age of R4810K using coalescent simulation.

RESULTS

The diversity of the RNF213 gene was higher in Africans than non-Africans, which can be explained by bottleneck effect of the out-of-Africa migration. Coalescent simulation showed that the risk variant was born in East Asia 14,500-5100 years ago and came to the Japanese archipelago afterward, probably in the period when the known migration based on archaeological evidences occurred.

CONCLUSIONS

Although clinical data show that the symptoms varies, all sequences harboring the risk allele are almost identical with a small number of exceptions, suggesting the MMD phenotypes are unaffected by the variants of this gene and rather would be more affected by environmental factors.

Genetic profile of adult T-cell leukemia/lymphoma in Okinawa: Association with prognosis, ethnicity, and HTLV-1 strains

Genetic alterations in adult T-cell leukemia/lymphoma (ATLL), a T-cell malignancy associated with HTLV-1, and their clinical impacts, especially from the perspective of viral strains, are not fully elucidated. We employed targeted next-generation sequencing and single nucleotide polymorphism array for 89 patients with ATLL in Okinawa, the southernmost islands in Japan, where the frequency of HTLV-1 tax subgroup-A (HTLV-1-taxA) is notably higher than that in mainland Japan, where most ATLL cases have HTLV-1-taxB, and compared the results with previously reported genomic landscapes of ATLL in mainland Japan and the USA. Okinawan patients exhibited similar mutation profiles to mainland Japanese patients, with frequent alterations in TCR/NF-ĸB (eg, PRKCB, PLCG1, and CARD11) and T-cell trafficking pathways (CCR4 and CCR7), in contrast with North American patients who exhibited a predominance of epigenome-associated gene mutations. Some mutations, especially GATA3 and RHOA, were detected more frequently in Okinawan patients than in mainland Japanese patients. Compared to HTLV-1-taxB, HTLV-1-taxA was significantly dominant in Okinawan patients with these mutations (GATA3, 34.1% vs 14.6%, P = .044; RHOA, 24.4% vs 6.3%, P = .032), suggesting the contribution of viral strains to these mutation frequencies. From a clinical viewpoint, we identified a significant negative impact of biallelic inactivation of PRDM1 (P = .027) in addition to the previously reported PRKCB mutations, indicating the importance of integrated genetic analysis. This study suggests that heterogeneous genetic abnormalities in ATLL depend on the viral strain as well as on the ethnic background. This warrants the need to develop therapeutic interventions considering regional characteristics.

Population genomics of East Asian ethnic groups

East Asia constitutes one-fifth of the global population and exhibits substantial genetic diversity. However, genetic investigations on populations in this region have been largely under-represented compared with European populations. Nonetheless, the last decade has seen considerable efforts and progress in genome-wide genotyping and whole-genome sequencing of the East-Asian ethnic groups. Here, we review the recent studies in terms of ancestral origin, population relationship, genetic differentiation, and admixture of major East- Asian groups, such as the Chinese, Korean, and Japanese populations. We mainly focus on insights from the whole-genome sequence data and also include the recent progress based on mitochondrial DNA (mtDNA) and Y chromosome data. We further discuss the evolutionary forces driving genetic diversity in East-Asian populations, and provide our perspectives for future directions on population genetics studies, particularly on underrepresented indigenous groups in East Asia.

Ancient Jomon genome sequence analysis sheds light on migration patterns of early East Asian populations

Anatomically modern humans reached East Asia more than 40,000 years ago. However, key questions still remain unanswered with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of the region may elucidate a more detailed picture of the peopling of East Eurasia. Here, we analyze the whole-genome sequence of a 2,500-year-old individual (IK002) from the main-island of Japan that is  characterized with a typical Jomon culture. The phylogenetic analyses support multiple waves of migration, with IK002 forming a basal lineage to the East and Northeast Asian genomes examined, likely representing some of the earliest-wave migrants who went north from Southeast Asia to East Asia. Furthermore, IK002 shows strong genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration. This study highlights the power of ancient genomics to provide new insights into the complex history of human migration into East Eurasia.

Takashi Gakuhari, Shigeki Nakagome et al. report the genomic analysis on a 2.5 kya individual from the ancient Jomon culture in present-day Japan. Phylogenetic analysis with comparison to other Eurasian sequences suggests early migration patterns in Asia and provides insight into the genetic affinities between peoples of the region.

Evolutionary History of the Risk of SNPs for Diffuse-Type Gastric Cancer in the Japanese Population

A genome wide association study reported that the T allele of rs2294008 in a cancer-related gene, PSCA, is a risk allele for diffuse-type gastric cancer. This allele has the highest frequency (0.63) in Japanese in Tokyo (JPT) among 26 populations in the 1000 Genomes Project database. F_ST ≈ 0.26 at this single nucleotide polymorphism is one of the highest between JPT and the genetically close Han Chinese in Beijing (CHB). To understand the evolutionary history of the alleles in PSCA, we addressed: (i) whether the C non-risk allele at rs2294008 is under positive selection, and (ii) why the mainland Japanese population has a higher T allele frequency than other populations. We found that haplotypes harboring the C allele are composed of two subhaplotypes. We detected that positive selection on both subhaplotypes has occurred in the East Asian lineage. However, the selection on one of the subhaplotypes in JPT seems to have been relaxed or ceased after divergence from the continental population; this may have caused the elevation of T allele frequency. Based on simulations under the dual structure model (a specific demography for the Japanese) and phylogenetic analysis with ancient DNA, the T allele at rs2294008 might have had high frequency in the Jomon people (one of the ancestral populations of the modern Japanese); this may explain the high T allele frequency in the extant Japanese.

The evolving Japanese: the dual structure hypothesis at 30

The population history of Japan has been one of the most intensively studied anthropological questions anywhere in the world, with a huge literature dating back to the nineteenth century and before. A growing consensus over the 1980s that the modern Japanese comprise an admixture of a Neolithic population with Bronze Age migrants from the Korean peninsula was crystallised in Kazurō Hanihara's influential 'dual structure hypothesis' published in 1991. Here, we use recent research in biological anthropology, historical linguistics and archaeology to evaluate this hypothesis after three decades. Although the major assumptions of Hanihara's model have been supported by recent work, we discuss areas where new findings have led to a re-evaluation of aspects of the hypothesis and emphasise the need for further research in key areas including ancient DNA and archaeology.

Analysis of whole Y-chromosome sequences reveals the Japanese population history in the Jomon period

The Jomon and the Yayoi are considered to be the two major ancestral populations of the modern mainland Japanese. The Jomon people, who inhabited mainland Japan, admixed with Yayoi immigrants from the Asian continent. To investigate the population history in the Jomon period (14,500–2,300 years before present [YBP]), we analyzed whole Y-chromosome sequences of 345 Japanese males living in mainland Japan. A phylogenetic analysis of East Asian Y chromosomes identified a major clade (35.4% of mainland Japanese) consisting of only Japanese Y chromosomes, which seem to have originated from indigenous Jomon people. A Monte Carlo simulation indicated that ~70% of Jomon males had Y chromosomes in this clade. The Bayesian skyline plots of 122 Japanese Y chromosomes in the clade detected a marked decrease followed by a subsequent increase in the male population size from around the end of the Jomon period to the beginning of the Yayoi period (2,300 YBP). The colder climate in the Late to Final Jomon period may have resulted in critical shortages of food for the Jomon people, who were hunter-gatherers, and the rice farming introduced by Yayoi immigrants may have helped the population size of the Jomon people to recover.

The geographical configuration of a language area influences linguistic diversity

Like the transfer of genetic variation through gene flow, language changes constantly as a result of its use in human interaction. Contact between speakers is most likely to happen when they are close in space, time, and social setting. Here, we investigated the role of geographical configuration in this process by studying linguistic diversity in Japan, which comprises a large connected mainland (less isolation, more potential contact) and smaller island clusters of the Ryukyuan archipelago (more isolation, less potential contact). We quantified linguistic diversity using dialectometric methods, and performed regression analyses to assess the extent to which distance in space and time predict contemporary linguistic diversity. We found that language diversity in general increases as geographic distance increases and as time passes-as with biodiversity. Moreover, we found that (I) for mainland languages, linguistic diversity is most strongly related to geographic distance-a so-called isolation-by-distance pattern, and that (II) for island languages, linguistic diversity reflects the time since varieties separated and diverged-an isolation-by-colonisation pattern. Together, these results confirm previous findings that (linguistic) diversity is shaped by distance, but also goes beyond this by demonstrating the critical role of geographic configuration.

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.

Fine-tuning of Approximate Bayesian Computation for human population genomics

Approximate Bayesian Computation (ABC) is a flexible statistical tool widely applied to addressing a variety of questions regarding the origin and evolution of humans. The significant growth of genomic scale data from diverse geographic populations has facilitated the use of ABC in modelling the complex processes that underlie human demography and local adaptation. However, a fundamental issue still remains in how to efficiently capture patterns of genetic variation with a set of summary statistics in order to achieve better approximation of Bayesian inference. Here, we review recent advances in ABC methodology and its applications for human population genomics, with a particular focus on optimal tuning of ABC approaches for different types of genetic data and different sets of evolutionary parameters.

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.

Japaneseplex: A forensic SNP assay for identification of Japanese people using Japanese-specific alleles

It is sometimes necessary to determine whether a forensic biological sample came from a Japanese person. In this study, we developed a 60-locus SNP assay designed for the differentiation of Japanese people from other East Asians using entirely and nearly Japanese-specific alleles. This multiplex assay consisted of 6 independent PCR reactions followed by single nucleotide extension. The average number and standard deviation of Japanese-specific alleles possessed by an individual were 0.81 ± 0.93 in 108 Koreans from Seoul, 8.87 ± 2.89 in 103 Japanese from Tottori, 17.20 ± 3.80 in 88 Japanese from Okinawa, and 0 in 220 Han Chinese from Wuxi and Changsha. The Koreans had 0-4 Japanese-specific alleles per individual, whereas the Japanese had 4-26 Japanese-specific alleles. Almost all Japanese were distinguished from the Koreans and other people by the factorial correspondence and principal component analyses. The Snipper program was also useful to estimate the degree of Japaneseness. The method described here was successfully applied to the differentiation of Japanese from non-Japanese people in forensic cases. This Japanese-specific SNP assay was named Japaneseplex.

High Prevalence of Distinct Human Herpesvirus 8 Contributes to the High Incidence of Non-acquired Immune Deficiency Syndrome-Associated Kaposi's Sarcoma in Isolated Japanese Islands

Background

Non-acquired immune deficiency syndrome (AIDS) Kaposi's sarcoma (KS) is extremely rare in Japan but highly endemic in Okinawa, especially in Miyako Islands. We aimed to elucidate the exact incidence and cause of this high prevalence.

Methods

Non-AIDS KS cases in Okinawa Prefecture over the past 31 years were reviewed, and human herpesvirus 8 (HHV8) seroprevalence in Miyako Islands was determined. We examined whole-genome sequences of 3 HHV8 strains and performed whole-exome sequencing of 4 male patients from Miyako Islands.

Results

Approximately half of the non-AIDS KS cases in Okinawa Prefecture were from Miyako Islands. The age-adjusted incidence rate was 0.87/105 per year for Miyako Islands and 0.056/105 per year for the rest of Okinawa. Human herpesvirus 8 seroprevalence was 15.4% in Miyako Islands. The 3 HHV8 genomes isolated from Miyako islanders formed a phylogenetically branch distinct from those of previously sequenced HHV8 strains and shared specific mutations in 9 proteins. These mutations were verified in Okinawan patients other than those from Miyako Islands. Whole-exome sequencing of the 4 male Miyako Islanders did not reveal shared pathogenic mutations.

Conclusions

Miyako Islands are an endemic area of non-AIDS KS. The high rate of a distinct HHV8 may contribute to the high incidence of KS in the region.

An assessment of postcranial indices, ratios, and body mass versus eco-geographical variables of prehistoric Jomon, Yayoi agriculturalists, and Kumejima Islanders of Japan

OBJECTIVES

We explore variations in body and limb proportions of the Jomon hunter-gatherers (14,000-2500 BP), the Yayoi agriculturalists (2500-1700 BP) of Japan, and the Kumejima Islanders of the Ryukyus (1600-1800 AD) with 11 geographically diverse skeletal postcranial samples from Africa, Europe, Asia, Australia, and North America using brachial-crural indices, femur head-breadth-to-femur length ratio, femur head-breadth-to-lower-limb-length ratio, and body mass as indicators of phenotypic climatic adaptation. Specifically, we test the hypothesis that variation in limb proportions seen in Jomon, Yayoi, and Kumejima is a complex interaction of genetic adaptation; development and allometric constraints; selection, gene flow and genetic drift with changing cultural factors (i.e., nutrition) and climate.

METHODS

The skeletal data (1127 individuals) were subjected to principle components analysis, Manly's permutation multiple regression tests, and Relethford-Blangero analysis.

RESULTS

The results of Manly's tests indicate that body proportions and body mass are significantly correlated with latitude, and minimum and maximum temperatures while limb proportions were not significantly correlated with these climatic variables. Principal components plots separated "climatic zones:" tropical, temperate, and arctic populations. The indigenous Jomon showed cold-adapted body proportions and warm-adapted limb proportions. Kumejima showed cold-adapted body proportions and limbs. The Yayoi adhered to the Allen-Bergmann expectation of cold-adapted body and limb proportions. Relethford-Blangero analysis showed that Kumejima experienced gene flow indicated by high observed variances while Jomon experienced genetic drift indicated by low observed variances.

CONCLUSIONS

The complex interaction of evolutionary forces and development/nutritional constraints are implicated in the mismatch of limb and body proportions.

Experimental evidence reveals the UCP1 genotype changes the oxygen consumption attributed to non-shivering thermogenesis in humans

Humans have spread out all over the world adapting to many different cold environments. Recent worldwide genome analyses and animal experiments have reported dozens of genes associated with cold adaptation. The uncoupling protein 1 (UCP1) gene enhances thermogenesis reaction in a physiological process by blocking ATP (adenosine triphosphate) synthesis on a mitochondrial membrane in brown adipose tissues. To our knowledge, no previous studies have shown an association between variants of the UCP1 gene and physiological phenotypes concerning non-shivering thermogenesis (NST) under the condition of low temperature in humans. We showed that the degree of NST for healthy subjects in an artificial climate chamber is significantly different among UCP1 genotypes. Defining the haplotypes covering the UCP1 region (39.4 kb), we found that the frequency of the haplotype with the highest NST was significantly correlated with latitudes and ambient temperature. Thus, the data in this study provide the first evidence that the UCP1 genotype alters the efficiency of NST in humans, and likely supports the hypothesis that the UCP1 gene has been related to cold adaptation in human evolutionary history.

Confounding effects of microbiome on the susceptibility of TNFSF15 to Crohn's disease in the Ryukyu Islands

Crohn's disease (CD) involves chronic inflammation in the gastrointestinal tract due to dysregulation of the host immune response to the gut microbiome. Even though the host-microbiome interactions are likely contributors to the development of CD, a few studies have detected genetic variants that change bacterial compositions and increase CD risk. We focus on one of the well-replicated susceptible genes, tumor necrosis factor superfamily member 15 (TNFSF15), and apply statistical analyses for personal profiles of genotypes and salivary microbiota collected from CD cases and controls in the Ryukyu Islands, southernmost islands of the Japanese archipelago. Our association test confirmed the susceptibility of TNFSF15 in the Ryukyu Islands. We found that the recessive model was supported to fit the observed genotype frequency of risk alleles slightly better than the additive model, defining the genetic effect on CD if a pair of the chromosomes in an individual consists of all risk alleles. The combined analysis of haplotypes and salivary microbiome from a small set of samples showed a significant association of the genetic effect with the increase of Prevotella, which led to a significant increase of CD risk. However, the genetic effect on CD disappeared if the abundance of Prevotella was low, suggesting the genetic contribution to CD is conditionally independent given a fixed amount of Prevotella. Although our statistical power is limited due to the small sample size, these results support an idea that the genetic susceptibility of TNFSF15 to CD may be confounded, in part, by the increase of Prevotella.

Characterising private and shared signatures of positive selection in 37 Asian populations

The Asian Diversity Project (ADP) assembled 37 cosmopolitan and ethnic minority populations in Asia that have been densely genotyped across over half a million markers to study patterns of genetic diversity and positive natural selection. We performed population structure analyses of the ADP populations and divided these populations into four major groups based on their genographic information. By applying a highly sensitive algorithm haploPS to locate genomic signatures of positive selection, 140 distinct genomic regions exhibiting evidence of positive selection in at least one population were identified. We examined the extent of signal sharing for regions that were selected in multiple populations and observed that populations clustered in a similar fashion to that of how the ancestry clades were phylogenetically defined. In particular, populations predominantly located in South Asia underwent considerably different adaptation as compared with populations from the other geographical regions. Signatures of positive selection present in multiple geographical regions were predicted to be older and have emerged prior to the separation of the populations in the different regions. In contrast, selection signals present in a single population group tended to be of lower frequencies and thus can be attributed to recent evolutionary events.

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

OBJECTIVES

A cline of frequencies of the derived allele of the ALDH2 gene, which causes a deficiency of an enzyme and "facial flushing" in humans who drink alcohol, has been known among the people of the Japanese archipelago. This cline is conventionally explained by admixture with immigrants from the Asian continent occurring during the Yayoi period. Previous studies lack sufficient data from the peripheral regions of the indigenous Jomon people, and those data the ADH1B gene that is involved in the Class I ADH gene cluster and contains another variant leading to a functional change.

METHODS

We focused on the southwestern-most people from the Ryukyu Islands (n = 218) and those from northern Kyushu (n = 21) where the Yayoi immigrants likely arrived. We investigated both the Class I ADH and ALDH2 loci, as well as neutral genetic markers.

RESULTS

In the Ryukyu Islands, the frequencies of the ancestral alleles in both loci were always higher than those in mainland Japan, while the frequencies of ADH1B were less than those of the derived allele. A haplotype block was not observed in ALDH2 but was in Class I ADH.

DISCUSSION

Our data suggest that the derived allele of ALDH2 came with the Yayoi immigrants from the Asian continent to the Japanese archipelago. However, the derived allele of ADH1B is unlikely to be related to the Yayoi migration. Therefore, we postulate that the expansion of the derived allele of ADHIB in East Asia could be traced back to the last glacial period.

Genotyping of the c.1423C>T (p.P475S) polymorphism in the ADAMTS13 gene by APLP and HRM assays: Northeastern Asian origin of the mutant

ADAMTS13 is a von Willebrand factor-cleaving protease. The mutant types of p.P475S (c.1423C>T) polymorphism in ADAMTS13 have a reduced activity in comparison with the wild type. In the present study, we investigated the frequency of the C-to-T substitution in 2584 genomic DNA samples from 25 Asian, European, and African populations using APLP (amplified product length polymorphism) and/or HRM (high-resolution melting) assays. Allele T (ADAMTS13(∗)T) was detected only in Asian populations and its frequency was observed to decrease gradually from north to south in 24 East Asian populations. Almost all ADAMTS13(∗)T were associated with ABO(∗)O. These results suggested that ADAMTS13(∗)T had occurred on a chromosome with ABO(∗)O in a northern part of East Asia. This SNP is useful as an ancestry-informative marker, and the present genotyping techniques are applicable to the investigation of an association between this SNP and aortic dissection (Kobayashi et al., 2012).

Mapping the genetic diversity of HLA haplotypes in the Japanese populations

Japan has often been viewed as an Asian country that possesses a genetically homogenous community. The basis for partitioning the country into prefectures has largely been geographical, although cultural and linguistic differences still exist between some of the districts/prefectures, especially between Okinawa and the mainland prefectures. The Major Histocompatibility Complex (MHC) region has consistently emerged as the most polymorphic region in the human genome, harbouring numerous biologically important variants; nevertheless the presence of population-specific long haplotypes hinders the imputation of SNPs and classical HLA alleles. Here, we examined the extent of genetic variation at the MHC between eight Japanese populations sampled from Okinawa, and six other prefectures located in or close to the mainland of Japan, specifically focusing at the haplotypes observed within each population, and what the impact of any variation has on imputation. Our results indicated that Okinawa was genetically farther to the mainland Japanese than were Gujarati Indians from Tamil Indians, while the mainland Japanese from six prefectures were more homogeneous than between northern and southern Han Chinese. The distribution of haplotypes across Japan was similar, although imputation was most accurate for Okinawa and several mainland prefectures when population-specific panels were used as reference.

Unique characteristics of the Ainu population in Northern Japan

Various genetic data (classic markers, mitochondrial DNAs, Y chromosomes and genome-wide single-nucleotide polymorphisms (SNPs)) have confirmed the coexistence of three major human populations on the Japanese Archipelago: Ainu in Hokkaido, Ryukyuans in the Southern Islands and Mainland Japanese. We compared genome-wide SNP data of the Ainu, Ryukyuans and Mainland Japanese, and found the following results: (1) the Ainu are genetically different from Mainland Japanese living in Tohoku, the northern part of Honshu Island; (2) using Ainu as descendants of the Jomon people and continental Asians (Han Chinese, Koreans) as descendants of Yayoi people, the proportion of Jomon genetic component in Mainland Japanese was ~18% and ~28% in Ryukyuans; (3) the time since admixture for Mainland Japanese ranged from 55 to 58 generations ago, and 43 to 44 generations ago for the Ryukyuans, depending on the number of Ainu individuals with varying rates of recent admixture with Mainland Japanese; (4) estimated haplotypes of some Ainu individuals suggested relatively long-term admixture with Mainland Japanese; and (5) highly differentiated genomic regions between Ainu and Mainland Japanese included EDAR and COL7A1 gene regions, which were shown to influence macroscopic phenotypes. These results clearly demonstrate the unique status of the Ainu and Ryukyuan people within East Asia.

Model-based verification of hypotheses on the origin of modern Japanese revisited by Bayesian inference based on genome-wide SNP data

Various hypotheses for the peopling of the Japanese archipelago have been proposed, which can be classified into three models: transformation, replacement, and hybridization. In recent years, one of the hybridization models ("dual-structure model") has been widely accepted. According to this model, Neolithic hunter-gatherers known as Jomon, who are assumed to have originated in southeast Asia and lived in the Japanese archipelago greater than 10,000 years ago, admixed with an agricultural people known as Yayoi, whom were migrants from the East Asian continent 2,000-3,000 years ago. Meanwhile, some anthropologists propose that rather, morphological differences between the Jomon and Yayoi people can be explained by microevolution following the lifestyle change. To resolve this controversy, we compared three demographic models by approximate Bayesian computation using genome-wide single nucleotide polymorphism (gwSNP) data from the Ainu people who are thought to be direct descendants of indigenous Jomon. If we assume Chinese people sampled in Beijing from HapMap have the same ancestry as Yayoi, then the hybridization model is predicted to be between 29 and 63 times more likely than the replacement and transformation models, respectively. Furthermore, our data provide strong support for a model in which the Jomon lineages had population structure diversified in local areas before the admixture event. Initial divergence between the Jomon and Yayoi ancestries was dated to late Pleistocene, followed by the divergence of Jomon lineages at early Holocene. These results suggest gwSNP data provides a detailed picture of the complex hybridization model for Japanese population history.