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

A variant in HMMR/HMMR-AS1 is associated with serum alanine aminotransferase levels in the Ryukyu population

The Ryukyu archipelago is located southwest of the Japanese islands, and people originally from this region, the Ryukyu population, have a unique genetic background distinct from that of other populations, including people from mainland Japan. However, few genetic studies have focused on the Ryukyu population. In this study, we performed genome-wide association studies (GWAS) on the serum levels of alanine aminotransferase (ALT, n = 15,224), aspartate aminotransferase (AST, n = 15,203), and gamma-glutamyl transferase (GGT, n = 14,496) in the Ryukyu population. We found 13 loci with a genome-wide significant association (P < 5 × 10-8), three for ALT, four for AST, and six for GGT, including one novel locus associated with ALT: rs117595134-A in HMMR/HMMR-AS1, ß =  - 0.131, standard error = 0.024, P = 4.90 × 10-8. Rs117595134-A is common in the Japanese population but is not observed in other ethnic populations in the 1000 genomes database. Additionally, 77 of 80 loci derived from Korean GWAS and 541 of 716 loci from European GWAS showed the same directions of effect (P = 1.41 × 10-19, P = 2.50 × 10-44, binomial test), indicating that most of susceptibility loci are shared between the Ryukyu population and other ethnic populations.

Genome-wide association studies for pelvic organ prolapse in the Japanese population

Pelvic organ prolapse (POP) affects approximately 40% of elderly women, characterized by the descent of the pelvic organs into the vaginal cavity. Here we present the results of a genome-wide association study (GWAS) for susceptibility to POP comprising 771 cases and 76,625 controls in the Japanese population. We identified a significant association of WT1 locus with POP in the Japanese population; rs10742277; odds ratio (OR) = 1.48, 95% confidence interval (CI), 1.29-1.68, P = 6.72 × 10-9. Subsequent cross-ancestry GWAS meta-analysis combining the Japanese data and previously reported European data, including 28,857 cases and 622,916 controls, identified FGFR2 locus as a novel susceptibility locus to POP (rs7072877; OR = 1.06, 95% CI, 1.04-1.08, P = 4.11 × 10-8). We also observed consistent directions of the effects for 21 out of 24 European GWAS derived loci (binomial test P = 2.8 × 10-4), indicating that most of susceptibility loci for POP are shared across the Japanese and European populations.

Region-Specific Genetic Diversity of Black Rats (Rattus rattus Complex) in Southeast and East Asia Shaped by Rapid Population Expansion Events

Among the six mitochondrial DNA lineages of the black rat (Rattus rattus Complex; RrC), lineages II and IV are widespread in Southeast and East Asia. This study explored their demographic history using 17 new sequences from the Miyako Islands in the Ryukyu archipelago, together with 178 publicly available cytochrome b sequences. We defined six and two haplotype groups showing rapid population expansion signals in Lineages II and IV, respectively. The six haplotype groups of Lineage II were represented by haplotypes from 1) Myanmar/Bangladesh/Northeast India, 2) Laos, 3) Thailand, 4) Indonesia/Philippines, 5) Vietnam/southern China, and 6) the Ryukyu archipelago. These expansion times were estimated using time-dependent evolutionary rates to be 115,300 years ago (ya), 128,500 ya, 9600 ya, 10,600 ya, 7200 ya, and 1400 ya, respectively, although all had large confidence intervals. The two groups of Lineage IV were recovered from the mainland and islands of Southeast Asia with predicted expansion times of 197,000 ya and 5800 ya, respectively. These results suggest that climatic fluctuations during the last 200,000 years of the Quaternary, affected the population dynamics in subtropical areas at different times. Furthermore, the results of the younger rapid expansion events of RrC suggest the possibility of agricultural advancement and dispersal of Neolithic farmers to different areas within the mainland and islands of Southeast Asia during the Holocene. A subset of rats from the Miyako Islands were found to have the same lineage IV haplotypes as those in Southeast Asia, suggesting a recent introduction of these new lineages.

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.

The Newfoundland and Labrador mosaic founder population descends from an Irish and British diaspora from 300 years ago

The founder population of Newfoundland and Labrador (NL) is a unique genetic resource, in part due to its geographic and cultural isolation, where historical records describe a migration of European settlers, primarily from Ireland and England, to NL in the 18th and 19th centuries. Whilst its historical isolation, and increased prevalence of certain monogenic disorders are well appreciated, details of the fine-scale genetic structure and ancestry of the population are lacking. Understanding the genetic origins and background of functional, disease causing, genetic variants would aid genetic mapping efforts in the Province. Here, we leverage dense genome-wide SNP data on 1,807 NL individuals to reveal fine-scale genetic structure in NL that is clustered around coastal communities and correlated with Christian denomination. We show that the majority of NL European ancestry can be traced back to the south-east and south-west of Ireland and England, respectively. We date a substantial population size bottleneck approximately 10-15 generations ago in NL, associated with increased haplotype sharing and autozygosity. Our results reveal insights into the population history of NL and demonstrate evidence of a population conducive to further genetic studies and biomarker discovery.

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.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

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.