Generation of chimpanzee induced pluripotent stem cell lines for cross-species comparisons.

As humans' closest living relatives, chimpanzees offer valuable insights into human evolution. However, technical and ethical limitations hinder investigations into the molecular and cellular foundations that distinguish chimpanzee and human traits. Recently, induced pluripotent stem cells (iPSCs) have emerged as a novel model for functional comparative studies and provided a non-invasive alternative for studying embryonic phenomena. In this study, we generated five new chimpanzee iPSC lines from peripheral blood cells and skin fibroblasts with SeV vectors carrying four reprogramming factors (human OCT3/4, SOX2, KLF4, and L-MYC) and characterized their pluripotency and differentiation potential. We also examined the expression of a human-specific non-coding RNA, HSTR1, which is predicted to be involved in human brain development. Our results show that the chimpanzee iPSCs possess pluripotent characteristics and can differentiate into various cell lineages. Moreover, we found that HSTR1 is expressed in human iPSCs and their neural derivatives but not in chimpanzee counterparts, supporting its possible role in human-specific brain development. As iPSCs are inherently variable due to genetic and epigenetic differences in donor cells or reprogramming procedures, it is essential to expand the number of chimpanzee iPSC lines to comprehensively capture the molecular and cellular properties representative of chimpanzees. Hence, our cells provide a valuable resource for investigating the function and regulation of human-specific transcripts such as HSTR1 and for understanding human evolution more generally.

Metagenomic analyses of 7000 to 5500 years old coprolites excavated from the Torihama shell-mound site in the Japanese archipelago.

Coprolites contain various kinds of ancient DNAs derived from gut micro-organisms, viruses, and foods, which can help to determine the gut environment of ancient peoples. Their genomic information should be helpful in elucidating the interaction between hosts and microbes for thousands of years, as well as characterizing the dietary behaviors of ancient people. We performed shotgun metagenomic sequencing on four coprolites excavated from the Torihama shell-mound site in the Japanese archipelago. The coprolites were found in the layers of the Early Jomon period, corresponding stratigraphically to 7000 to 5500 years ago. After shotgun sequencing, we found that a significant number of reads showed homology with known gut microbe, viruses, and food genomes typically found in the feces of modern humans. We detected reads derived from several types of phages and their host bacteria simultaneously, suggesting the coexistence of viruses and their hosts. The food genomes provide biological evidence for the dietary behavior of the Jomon people, consistent with previous archaeological findings. These results indicate that ancient genomic analysis of coprolites is useful for understanding the gut environment and lifestyle of ancient peoples.

Infectious diseases may have arrested the southward advance of microblades in Upper Palaeolithic East Asia.

An unsolved archaeological puzzle of the East Asian Upper Palaeolithic is why the southward expansion of an innovative lithic technology represented by microblades stalled at the Qinling-Huaihe Line. It has been suggested that the southward migration of foragers with microblades stopped there, which is consistent with ancient DNA studies showing that populations to the north and south of this line had differentiated genetically by 19 000 years ago. Many infectious pathogens are believed to have been associated with hominins since the Palaeolithic, and zoonotic pathogens in particular are prevalent at lower latitudes, which may have produced a disease barrier. We propose a mathematical model to argue that mortality due to infectious diseases may have arrested the wave-of-advance of the technologically advantaged foragers from the north.

Endocrine, inflammatory and immune responses and individual differences in acute hypobaric hypoxia in lowlanders.

When lowlanders are exposed to environments inducing hypobaric hypoxia (HH) such as high mountains, hemodynamic changes occur to maintain oxygen levels in the body. However, changes to other physiological functions under such conditions have yet to be clarified. This study investigated changes in endocrine, inflammatory and immune parameters and individual differences during acute HH exposure using a climatic chamber (75 min of exposure to conditions mimicking 3500 m) in healthy lowlanders. Aldosterone and cortisol were significantly decreased and interleukin (IL)-6, IL-8 and white blood cell (WBC) counts were significantly increased after HH. Lower peripheral oxygen saturation (SpO2) was associated with higher IL-6 and WBC counts, and higher IL-8 was associated with higher cortisol. These findings suggest that endocrine, inflammatory and immune responses are evoked even with a short 75-min exposure to HH and individuals with lower SpO2 seemed to show more pronounced responses. Our results provide basic data for understanding the physiological responses and interactions of homeostatic systems during acute HH.

Genomes from Verteba cave suggest diversity within the Trypillians in Ukraine.

The transition to agriculture occurred relatively late in Eastern Europe, leading researchers to debate whether it was a gradual, interactive process or a colonisation event. In the forest and forest-steppe regions of Ukraine, farming appeared during the fifth millennium BCE, associated with the Cucuteni-Trypillia cultural complex (CTCC, ~ 5000-3000 BCE). Across Europe, the Neolithisation process was highly variable across space and over time. Here, we investigate the population dynamics of early agriculturalists from the eastern forest-steppe region based on the analyses of 20 ancient genomes from the site of Verteba Cave (3935-825 cal BCE). Results reveal that the CTCC individuals' ancestry is related to both western hunter-gatherers and Near Eastern farmers, has no local ancestry associated with Ukrainian Neolithic hunter-gatherers and has steppe ancestry. An Early Bronze Age individual has an ancestry profile related to the Yamnaya expansions but with 20% of ancestry related to the other Trypillian individuals, which suggests admixture between the Trypillians and the incoming populations carrying steppe-related ancestry. A Late Bronze Age individual dated to 980-825 cal BCE has a genetic profile indicating affinity to Beaker-related populations, detected close to 1000 years after the end of the Bell Beaker phenomenon during the third millennium BCE.

Triangulation supports agricultural spread of the Transeurasian languages.

The origin and early dispersal of speakers of Transeurasian languages-that is, Japanese, Korean, Tungusic, Mongolic and Turkic-is among the most disputed issues of Eurasian population history1-3. A key problem is the relationship between linguistic dispersals, agricultural expansions and population movements4,5. Here we address this question by 'triangulating' genetics, archaeology and linguistics in a unified perspective. We report wide-ranging datasets from these disciplines, including a comprehensive Transeurasian agropastoral and basic vocabulary; an archaeological database of 255 Neolithic-Bronze Age sites from Northeast Asia; and a collection of ancient genomes from Korea, the Ryukyu islands and early cereal farmers in Japan, complementing previously published genomes from East Asia. Challenging the traditional 'pastoralist hypothesis'6-8, we show that the common ancestry and primary dispersals of Transeurasian languages can be traced back to the first farmers moving across Northeast Asia from the Early Neolithic onwards, but that this shared heritage has been masked by extensive cultural interaction since the Bronze Age. As well as marking considerable progress in the three individual disciplines, by combining their converging evidence we show that the early spread of Transeurasian speakers was driven by agriculture.

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.

Exploring correlations in genetic and cultural variation across language families in northeast Asia.

Culture evolves in ways that are analogous to, but distinct from, genomes. Previous studies examined similarities between cultural variation and genetic variation (population history) at small scales within language families, but few studies have empirically investigated these parallels across language families using diverse cultural data. We report an analysis comparing culture and genomes from in and around northeast Asia spanning 11 language families. We extract and summarize the variation in language (grammar, phonology, lexicon), music (song structure, performance style), and genomes (genome-wide SNPs) and test for correlations. We find that grammatical structure correlates with population history (genetic history). Recent contact and shared descent fail to explain the signal, suggesting relationships that arose before the formation of current families. Our results suggest that grammar might be a cultural indicator of population history while also demonstrating differences among cultural and genetic relationships that highlight the complex nature of human history.

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.

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.

Parallel Evolution of Two dmrt1-Derived Genes, dmy and dm-W, for Vertebrate Sex Determination.

Animal sex-determining genes, which bifurcate for female and male development, are diversified even among closely related species. Most of these genes emerged independently from various sex-related genes during species diversity as neofunctionalization-type genes. However, the common mechanisms of this divergent evolution remain poorly understood. Here, we compared the molecular evolution of two sex-determining genes, the medaka dmy and the clawed frog dm-W, which independently evolved from the duplication of the transcription factor-encoding masculinization gene dmrt1. Interestingly, we detected parallel amino acid substitutions, from serine (S) to threonine (T), on the DNA-binding domains of both ancestral DMY and DM-W, resulting from positive selection. Two types of DNA-protein binding experiments and a luciferase reporter assay demonstrated that these S-T substitutions could strengthen the DNA-binding abilities and enhance the transcriptional regulation function. These findings suggest that the parallel S-T substitutions may have contributed to the establishment of dmy and dm-W as sex-determining genes.

Medaka Population Genome Structure and Demographic History Described via Genotyping-by-Sequencing.

Medaka is a model organism in medicine, genetics, developmental biology and population genetics. Lab stocks composed of more than 100 local wild populations are available for research in these fields. Thus, medaka represents a potentially excellent bioresource for screening disease-risk- and adaptation-related genes in genome-wide association studies. Although the genetic population structure should be known before performing such an analysis, a comprehensive study on the genome-wide diversity of wild medaka populations has not been performed. Here, we performed genotyping-by-sequencing (GBS) for 81 and 12 medakas captured from a bioresource and the wild, respectively. Based on the GBS data, we evaluated the genetic population structure and estimated the demographic parameters using an approximate Bayesian computation (ABC) framework. The genome-wide data confirmed that there were substantial differences between local populations and supported our previously proposed hypothesis on medaka dispersal based on mitochondrial genome (mtDNA) data. A new finding was that a local group that was thought to be a hybrid between the northern and the southern Japanese groups was actually an origin of the northern Japanese group. Thus, this paper presents the first population-genomic study of medaka and reveals its population structure and history based on chromosomal genetic diversity.

A new targeted capture method using bacterial artificial chromosome (BAC) libraries as baits for sequencing relatively large genes.

To analyze a specific genome region using next-generation sequencing technologies, the enrichment of DNA libraries with targeted capture methods has been standardized. For enrichment of mitochondrial genome, a previous study developed an original targeted capture method that use baits constructed from long-range polymerase chain reaction (PCR) amplicons, common laboratory reagents, and equipment. In this study, a new targeted capture method is presented, that of bacterial artificial chromosome (BAC) double capture (BDC), modifying the previous method, but using BAC libraries as baits for sequencing a relatively large gene. We applied the BDC approach for the 214 kb autosomal region, ring finger protein 213, which is the susceptibility gene of moyamoya disease (MMD). To evaluate the reliability of BDC, cost and data quality were compared with those of a commercial kit. While the ratio of duplicate reads was higher, the cost was less than that of the commercial kit. The data quality was sufficiently the same as that of the kit. Thus, BDC can be an easy, low-cost, and useful method for analyzing individual genome regions with substantial length.

Abnormal nuclear morphology is independent of longevity in a zmpste24-deficient fish model of Hutchinson-Gilford progeria syndrome (HGPS).

Lamin is an intermediate protein underlying the nuclear envelope and it plays a key role in maintaining the integrity of the nucleus. A defect in the processing of its precursor by a metalloprotease, ZMPSTE24, results in the accumulation of farnesylated prelamin in the nucleus and causes various diseases, including Hutchinson-Gilford progeria syndrome (HGPS). However, the role of lamin processing is unclear in fish species. Here, we generated zmpste24-deficient medaka and evaluated their phenotype. Unlike humans and mice, homozygous mutants did not show growth defects or lifespan shortening, despite lamin precursor accumulation. Gonadosomatic indices, blood glucose levels, and regenerative capacity of fins were similar in 1-year-old mutants and their wild-type (WT) siblings. Histological examination showed that the muscles, subcutaneous fat tissues, and gonads were normal in the mutants at the age of 1 year. However, the mutants showed hypersensitivity to X-ray irradiation, although p53target genes, p21 and mdm2, were induced 6 h after irradiation. Immunostaining of primary cultured cells from caudal fins and visualization of nuclei using H2B-GFP fusion proteins revealed an abnormal nuclear shape in the mutants both in vitro and in vivo. The telomere lengths were significantly shorter in the mutants compared to WT. Taken together, these results suggest that zmpste24-deficient medaka phenocopied HGPS only partially and that abnormal nuclear morphology and lifespan shortening are two independent events in vertebrates.

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.

An ancestral haplotype of the human PERIOD2 gene associates with reduced sensitivity to light-induced melatonin suppression.

Humans show various responses to the environmental stimulus in individual levels as "physiological variations." However, it has been unclear if these are caused by genetic variations. In this study, we examined the association between the physiological variation of response to light-stimulus and genetic polymorphisms. We collected physiological data from 43 subjects, including light-induced melatonin suppression, and performed haplotype analyses on the clock genes, PER2 and PER3, exhibiting geographical differentiation of allele frequencies. Among the haplotypes of PER3, no significant difference in light sensitivity was found. However, three common haplotypes of PER2 accounted for more than 96% of the chromosomes in subjects, and 1 of those 3 had a significantly low-sensitive response to light-stimulus (P < 0.05). The homozygote of the low-sensitive PER2 haplotype showed significantly lower percentages of melatonin suppression (P < 0.05), and the heterozygotes of the haplotypes varied their ratios, indicating that the physiological variation for light-sensitivity is evidently related to the PER2 polymorphism. Compared with global haplotype frequencies, the haplotype with a low-sensitive response was more frequent in Africans than in non-Africans, and came to the root in the phylogenetic tree, suggesting that the low light-sensitive haplotype is the ancestral type, whereas the other haplotypes with high sensitivity to light are the derived types. Hence, we speculate that the high light-sensitive haplotypes have spread throughout the world after the Out-of-Africa migration of modern humans.

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.

An Approach to Elucidate NBS1 Function in DNA Repair Using Frequent Nonsynonymous Polymorphism in Wild Medaka (Oryzias latipes) Populations.

Nbs1 is one of the genes responsible for Nijmegen breakage syndrome, which is marked with high radiosensitivity. In human NBS1 (hNBS1), Q185E polymorphism is known as the factor to cancer risks, although its DSB repair defect has not been addressed. Here we investigated the genetic variations in medaka (Oryzias latipes) wild populations, and found 40 nonsynonymous single nucleotide polymorphisms (SNPs) in medaka nbs1 (olnbs1) gene within 5 inbred strains. A mutation to histidine in Q170 residue in olNbs1, which corresponds to Q185 residue of hNBS1, was widely distributed in the closed colonies derived from the eastern Korean population of medaka. Overexpression of H170 type olNbs1 in medaka cultured cell lines resulted in the increased accumulation of olNbs1 at laser-induced DSB sites. Autophosphorylation of DNA-dependent protein kinase at T2609 was suppressed after the γ-ray irradiation, which was followed by prolonged formation of γ-H2AX foci and delayed DSB repair. These findings suggested that the nonsynonymous SNP (Q170H) in olnbs1, which induced DSB repair defects, is specifically distributed in the eastern Korean population of medaka. Furthermore, examination using the variation within wild populations might provide a novel method to characterize a driving force to spread the disease risk alleles.

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.