dbCNS: A New Database for Conserved Noncoding Sequences.

We developed dbCNS (http://yamasati.nig.ac.jp/dbcns), a new database for conserved noncoding sequences (CNSs). CNSs exist in many eukaryotes and are assumed to be involved in protein expression control. Version 1 of dbCNS, introduced here, includes a powerful and precise CNS identification pipeline for multiple vertebrate genomes. Mutations in CNSs may induce morphological changes and cause genetic diseases. For this reason, many vertebrate CNSs have been identified, with special reference to primate genomes. We integrated ∼6.9 million CNSs from many vertebrate genomes into dbCNS, which allows users to extract CNSs near genes of interest using keyword searches. In addition to CNSs, dbCNS contains published genome sequences of 161 species. With purposeful taxonomic sampling of genomes, users can employ CNSs as queries to reconstruct CNS alignments and phylogenetic trees, to evaluate CNS modifications, acquisitions, and losses, and to roughly identify species with CNSs having accelerated substitution rates. dbCNS also produces links to dbSNP for searching pathogenic single-nucleotide polymorphisms in human CNSs. Thus, dbCNS connects morphological changes with genetic diseases. A test analysis using 38 gnathostome genomes was accomplished within 30 s. dbCNS results can evaluate CNSs identified by other stand-alone programs using genome-scale data.

Inferring intelligence of ancient people based on modern genomic studies.

Quantification of ancient human intelligence has become possible with recent advances in polygenic prediction. Intelligence is a complex trait that has both environmental and genetic components and high heritability. Large-scale genome-wide association studies based on ~270,000 individuals have demonstrated highly significant single-nucleotide polymorphisms (SNPs) associated with intelligence in present-day humans. We utilized those previously reported 12,037 SNPs to estimate a genetic component of intelligence in ancient Funadomari Jomon individual from 3700 years BP as well as four individuals of Afanasievo nuclear family from about 4100 years BP and who are considered anatomically modern humans. We have demonstrated that ancient individuals could have been not inferior in intelligence compared to present-day humans through assessment of the genetic component of intelligence. We have also confirmed that alleles associated with intelligence tend to spread equally between ancestral and derived origin suggesting that intelligence may be a neutral trait in human evolution.

Genome-wide SNP data of Izumo and Makurazaki populations support inner-dual structure model for origin of Yamato people.

The "Dual Structure" model on the formation of the modern Japanese population assumes that the indigenous hunter-gathering population (symbolized as Jomon people) admixed with rice-farming population (symbolized as Yayoi people) who migrated from the Asian continent after the Yayoi period started. The Jomon component remained high both in Ainu and Okinawa people who mainly reside in northern and southern Japan, respectively, while the Yayoi component is higher in the mainland Japanese (Yamato people). The model has been well supported by genetic data, but the Yamato population was mostly represented by people from Tokyo area. We generated new genome-wide SNP data using Japonica Array for 45 individuals in Izumo City of Shimane Prefecture and for 72 individuals in Makurazaki City of Kagoshima Prefecture in Southern Kyushu, and compared these data with those of other human populations in East Asia, including BioBank Japan data. Using principal component analysis, phylogenetic network, and f4 tests, we found that Izumo, Makurazaki, and Tohoku populations are slightly differentiated from Kanto (including Tokyo), Tokai, and Kinki regions. These results suggest the substructure within Mainland Japanese maybe caused by multiple migration events from the Asian continent following the Jomon period, and we propose a modified version of "Dual Structure" model called the "Inner-Dual Structure" model.

House mouse Mus musculus dispersal in East Eurasia inferred from 98 newly determined complete mitochondrial genome sequences.

The Eurasian house mouse Mus musculus is useful for tracing prehistorical human movement related to the spread of farming. We determined whole mitochondrial DNA (mtDNA) sequences (ca. 16,000 bp) of 98 wild-derived individuals of two subspecies, M. m. musculus (MUS) and M. m. castaneus (CAS). We revealed directional dispersals reaching as far as the Japanese Archipelago from their homelands. Our phylogenetic analysis indicated that the eastward movement of MUS was characterised by five step-wise regional extension events: (1) broad spatial expansion into eastern Europe and the western part of western China, (2) dispersal to the eastern part of western China, (3) dispersal to northern China, (4) dispersal to the Korean Peninsula and (5) colonisation and expansion in the Japanese Archipelago. These events were estimated to have occurred during the last 2000-18,000 years. The dispersal of CAS was characterised by three events: initial divergences (ca. 7000-9000 years ago) of haplogroups in northernmost China and the eastern coast of India, followed by two population expansion events that likely originated from the Yangtze River basin to broad areas of South and Southeast Asia, including Sri Lanka, Bangladesh and Indonesia (ca. 4000-6000 years ago) and to Yunnan, southern China and the Japanese Archipelago (ca. 2000-3500). This study provides a solid framework for the spatiotemporal movement of the human-associated organisms in Holocene Eastern Eurasia using whole mtDNA sequences, reliable evolutionary rates and accurate branching patterns. The information obtained here contributes to the analysis of a variety of animals and plants associated with prehistoric human migration.

Chance, Finiteness, and History.

Importance of chance, finiteness, and history in evolution is pointed out with special reference to the neutral theory.

Genomic Locations of Conserved Noncoding Sequences and Their Proximal Protein-Coding Genes in Mammalian Expression Dynamics.

Experimental studies have found the involvement of certain conserved noncoding sequences (CNSs) in the regulation of the proximal protein-coding genes in mammals. However, reported cases of long range enhancer activities and inter-chromosomal regulation suggest that proximity of CNSs to protein-coding genes might not be important for regulation. To test the importance of the CNS genomic location, we extracted the CNSs conserved between chicken and four mammalian species (human, mouse, dog, and cattle). These CNSs were confirmed to be under purifying selection. The intergenic CNSs are often found in clusters in gene deserts, where protein-coding genes are in paucity. The distribution pattern, ChIP-Seq, and RNA-Seq data suggested that the CNSs are more likely to be regulatory elements and not corresponding to long intergenic noncoding RNAs. Physical distances between CNS and their nearest protein coding genes were well conserved between human and mouse genomes, and CNS-flanking genes were often found in evolutionarily conserved genomic neighborhoods. ChIP-Seq signal and gene expression patterns also suggested that CNSs regulate nearby genes. Interestingly, genes with more CNSs have more evolutionarily conserved expression than those with fewer CNSs. These computationally obtained results suggest that the genomic locations of CNSs are important for their regulatory functions. In fact, various kinds of evolutionary constraints may be acting to maintain the genomic locations of CNSs and protein-coding genes in mammals to ensure proper regulation.

A partial nuclear genome of the Jomons who lived 3000 years ago in Fukushima, Japan.

The Jomon period of the Japanese Archipelago, characterized by cord-marked 'jomon' potteries, has yielded abundant human skeletal remains. However, the genetic origins of the Jomon people and their relationships with modern populations have not been clarified. We determined a total of 115 million base pair nuclear genome sequences from two Jomon individuals (male and female each) from the Sanganji Shell Mound (dated 3000 years before present) with the Jomon-characteristic mitochondrial DNA haplogroup N9b, and compared these nuclear genome sequences with those of worldwide populations. We found that the Jomon population lineage is best considered to have diverged before diversification of present-day East Eurasian populations, with no evidence of gene flow events between the Jomon and other continental populations. This suggests that the Sanganji Jomon people descended from an early phase of population dispersals in East Asia. We also estimated that the modern mainland Japanese inherited <20% of Jomon peoples' genomes. Our findings, based on the first analysis of Jomon nuclear genome sequence data, firmly demonstrate that the modern mainland Japanese resulted from genetic admixture of the indigenous Jomon people and later migrants.

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.

Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway.

The mammalian Y chromosome has unique characteristics compared with the autosomes or X chromosomes. Here we report the finished sequence of the chimpanzee Y chromosome (PTRY), including 271 kb of the Y-specific pseudoautosomal region 1 and 12.7 Mb of the male-specific region of the Y chromosome. Greater sequence divergence between the human Y chromosome (HSAY) and PTRY (1.78%) than between their respective whole genomes (1.23%) confirmed the accelerated evolutionary rate of the Y chromosome. Each of the 19 PTRY protein-coding genes analyzed had at least one nonsynonymous substitution, and 11 genes had higher nonsynonymous substitution rates than synonymous ones, suggesting relaxation of selective constraint, positive selection or both. We also identified lineage-specific changes, including deletion of a 200-kb fragment from the pericentromeric region of HSAY, expansion of young Alu families in HSAY and accumulation of young L1 elements and long terminal repeat retrotransposons in PTRY. Reconstruction of the common ancestral Y chromosome reflects the dynamic changes in our genomes in the 5-6 million years since speciation.

A SNP in the ABCC11 gene is the determinant of human earwax type.

Human earwax consists of wet and dry types. Dry earwax is frequent in East Asians, whereas wet earwax is common in other populations. Here we show that a SNP, 538G --> A (rs17822931), in the ABCC11 gene is responsible for determination of earwax type. The AA genotype corresponds to dry earwax, and GA and GG to wet type. A 27-bp deletion in ABCC11 exon 29 was also found in a few individuals of Asian ancestry. A functional assay demonstrated that cells with allele A show a lower excretory activity for cGMP than those with allele G. The allele A frequency shows a north-south and east-west downward geographical gradient; worldwide, it is highest in Chinese and Koreans, and a common dry-type haplotype is retained among various ethnic populations. These suggest that the allele A arose in northeast Asia and thereafter spread through the world. The 538G --> A SNP is the first example of DNA polymorphism determining a visible genetic trait.

Denisova admixture and the first modern human dispersals into Southeast Asia and Oceania.

It has recently been shown that ancestors of New Guineans and Bougainville Islanders have inherited a proportion of their ancestry from Denisovans, an archaic hominin group from Siberia. However, only a sparse sampling of populations from Southeast Asia and Oceania were analyzed. Here, we quantify Denisova admixture in 33 additional populations from Asia and Oceania. Aboriginal Australians, Near Oceanians, Polynesians, Fijians, east Indonesians, and Mamanwa (a "Negrito" group from the Philippines) have all inherited genetic material from Denisovans, but mainland East Asians, western Indonesians, Jehai (a Negrito group from Malaysia), and Onge (a Negrito group from the Andaman Islands) have not. These results indicate that Denisova gene flow occurred into the common ancestors of New Guineans, Australians, and Mamanwa but not into the ancestors of the Jehai and Onge and suggest that relatives of present-day East Asians were not in Southeast Asia when the Denisova gene flow occurred. Our finding that descendants of the earliest inhabitants of Southeast Asia do not all harbor Denisova admixture is inconsistent with a history in which the Denisova interbreeding occurred in mainland Asia and then spread over Southeast Asia, leading to all its earliest modern human inhabitants. Instead, the data can be most parsimoniously explained if the Denisova gene flow occurred in Southeast Asia itself. Thus, archaic Denisovans must have lived over an extraordinarily broad geographic and ecological range, from Siberia to tropical Asia.

The functional A allele was resurrected via recombination in the human ABO blood group gene.

Functional A and B alleles are distinguished at two critical sites in exon 7 of the human ABO blood group gene. The most frequent nonfunctional O alleles have one-base deletion in exon 6 producing a frameshift, and it has the A type signature in two critical sites in exon 7. Previous studies indicated that B and O alleles were derived from A allele in human lineage. In this study, we conducted a phylogenetic network analysis using six representative haplotypes: A101, A201, B101, O01, O02, and O09. The result indicated that the A allele, possibly once extinct in the human lineage a long time ago, was resurrected by a recombination between B and O alleles less than 300,000 years ago.

Evolutionary history of continental southeast Asians: "early train" hypothesis based on genetic analysis of mitochondrial and autosomal DNA data.

The population history of the indigenous populations in island Southeast Asia is generally accepted to have been shaped by two major migrations: the ancient "Out of Africa" migration ∼50,000 years before present (YBP) and the relatively recent "Out of Taiwan" expansion of Austronesian agriculturalists approximately 5,000 YBP. The Negritos are believed to have originated from the ancient migration, whereas the majority of island Southeast Asians are associated with the Austronesian expansion. We determined 86 mitochondrial DNA (mtDNA) complete genome sequences in four indigenous Malaysian populations, together with a reanalysis of published autosomal single-nucleotide polymorphism (SNP) data of Southeast Asians to test the plausibility and impact of those migration models. The three Austronesian groups (Bidayuh, Selatar, and Temuan) showed high frequencies of mtDNA haplogroups, which originated from the Asian mainland ∼30,000-10,000 YBP, but low frequencies of "Out of Taiwan" markers. Principal component analysis and phylogenetic analysis using autosomal SNP data indicate a dichotomy between continental and island Austronesian groups. We argue that both the mtDNA and autosomal data suggest an "Early Train" migration originating from Indochina or South China around the late-Pleistocene to early-Holocene period, which predates, but may not necessarily exclude, the Austronesian expansion.

Regional DNA methylation differences between humans and chimpanzees are associated with genetic changes, transcriptional divergence and disease genes.

Changes in gene expression have been proposed to have an important role in the evolutionary changes in phenotypes. Interspecific changes in gene expression can result not only from genetic changes in regulatory regions but also from epigenetic changes in such regions. Here we report the identification of genomic regions showing differences in DNA methylation between humans and chimpanzees (termed S-DMRs for species-specific differentially methylated regions) on chromosomes 21 and 22. These regional methylation differences are frequently associated with genes, including those relevant to a disease, such as Alzheimer's disease, diabetes mellitus or cancer. Methylation differences are often correlated with changes in promoter activity or alternative splicing. Comparative studies including other great ape species provide evidence for the contribution of genetic changes to some of these S-DMRs. Genetic changes responsible for the S-DMRs include gain or loss of CTCF-binding site and changes in CpG density in microsatellite repeats. Our results suggest that DNA methylation changes, often caused by small sequence changes, contribute to transcriptional and phenotypic diversification in hominid evolution.

The PNarec method for detection of ancient recombinations through phylogenetic network analysis.

Recombinations are known to disrupt bifurcating tree structure of gene genealogies. Although recently occurred recombinations are easily detectable by using conventional methods, recombinations may have occurred at any time. We devised a new method for detecting ancient recombinations through phylogenetic network analysis, and detected five ancient recombinations in gibbon ABO blood group genes [Kitano et al., 2009. Mol. Phylogenet. Evol., 51, 465-471]. We present applications of this method, now named as "PNarec", to various virus sequences as well as HLA genes.

Admixture patterns and genetic differentiation in negrito groups from West Malaysia estimated from genome-wide SNP data.

Southeast Asia houses various culturally and linguistically diverse ethnic groups. In Malaysia, where the Malay, Chinese, and Indian ethnic groups form the majority, there exist minority groups such as the "negritos" who are believed to be descendants of the earliest settlers of Southeast Asia. Here we report patterns of genetic substructure and admixture in two Malaysian negrito populations (Jehai and Kensiu), using ~50,000 genome-wide single-nucleotide polymorphism (SNP) data. We found traces of recent admixture in both the negrito populations, particularly in the Jehai, with the Malay through principal component analysis and STRUCTURE analysis software, which suggested that the admixture was as recent as one generation ago. We also identified significantly differentiated nonsynonymous SNPs and haplotype blocks related to intracellular transport, metabolic processes, and detection of stimulus. These results highlight the different levels of admixture experienced by the two Malaysian negritos. Delineating admixture and differentiated genomic regions should be of importance in designing and interpretation of molecular anthropology and disease association studies.

A Study of the Genetic Structure of Hybrid Camels in Kazakhstan.

Camel farming is gaining scientific interest due to its unique agricultural characteristics. Camels are versatile for milk and meat production, wool, racing, transport, and tourism. To use their full potential, it is essential to improve our understanding of the genetic structure of these animals. One-humped and two-humped camels have received detailed genetic descriptions, while there is no such information for their hybrids, which outperform their parent species in several agricultural characteristics. Thus, in this study, for the first time, the whole genome sequencing data (WGS) of five hybrid camels bred in the Almaty region of Kazakhstan are presented in comparison with the WGS data of one-humped, two-humped, and wild camels. A total of 43,552,164 single-nucleotide polymorphisms were found across the studied groups. Further comparison of these SNPs showed the following number of private SNPs among the populations: hybrid camels (3,271,083), wild camels (2,515,591), Bactrians (1,244,694), and dromedaries (531,224). The genetic structure of the studied animals was described, and a phylogenetic tree was built to assess their genetic distance. It was found that the studied hybrids are genetically closer to dromedaries since they were on the close branch of the phylogenetic tree.

Loss-of-function mutation in a repressor module of human-specifically activated enhancer HACNS1.

The cis-regulatory element contributed to gaining humanness is of great interest in human evolutionary studies. A human-accelerated region exceeding neutral evolutionary rates, termed HACNS1, was recently reported as a positively selected sequence acquiring novel TF-binding sites responsible for human-specific gain of limb enhancer function. However, another possibility is loss of function in repressor element in HACNS1. Signature of the human substitutions in the 81-bp region infers that a GC-biased gene conversion (BGC) might create these seemingly excessive substitutions. To evaluate the 81-bp function, we performed transgenic mouse assay of the HACNS1 construct lacking the 81-bp region. The deleted construct showed similar enhancer activity to the intact human HACNS1, suggesting that the function of the human 81-bp region is not an activating enhancer but rather a disrupted repressor. This result infers that loss of function in the HACNS1 81-bp region, possibly via a BGC, played an important role in human-specific evolution.

The history of human populations in the Japanese Archipelago inferred from genome-wide SNP data with a special reference to the Ainu and the Ryukyuan populations.

The Japanese Archipelago stretches over 4000 km from north to south, and is the homeland of the three human populations; the Ainu, the Mainland Japanese and the Ryukyuan. The archeological evidence of human residence on this Archipelago goes back to >30 000 years, and various migration routes and root populations have been proposed. Here, we determined close to one million single-nucleotide polymorphisms (SNPs) for the Ainu and the Ryukyuan, and compared these with existing data sets. This is the first report of these genome-wide SNP data. Major findings are: (1) Recent admixture with the Mainland Japanese was observed for more than one third of the Ainu individuals from principal component analysis and frappe analyses; (2) The Ainu population seems to have experienced admixture with another population, and a combination of two types of admixtures is the unique characteristics of this population; (3) The Ainu and the Ryukyuan are tightly clustered with 100% bootstrap probability followed by the Mainland Japanese in the phylogenetic trees of East Eurasian populations. These results clearly support the dual structure model on the Japanese Archipelago populations, though the origins of the Jomon and the Yayoi people still remain to be solved.

Insights into Mus musculus Population Structure across Eurasia Revealed by Whole-Genome Analysis.

For more than 100 years, house mice (Mus musculus) have been used as a key animal model in biomedical research. House mice are genetically diverse, yet their genetic background at the global level has not been fully understood. Previous studies have suggested that they originated in South Asia and diverged into three major subspecies, almost simultaneously, approximately 110,000-500,000 years ago; however, they have spread across the world with the migration of modern humans in prehistoric and historic times (∼10,000 years ago to the present day) and have undergone secondary contact, which has complicated the genetic landscape of wild house mice. In this study, we sequenced the whole-genome sequences of 98 wild house mice collected from Eurasia, particularly East Asia, Southeast Asia, and South Asia. Although wild house mice were found to consist of three major genetic groups corresponding to the three major subspecies, individuals representing admixtures between subspecies were more prevalent in East Asia than has been previously recognized. Furthermore, several samples exhibited an incongruent pattern of genealogies between mitochondrial and autosomal genomes. Using samples that likely retained the original genetic components of subspecies with the least admixture, we estimated the pattern and timing of divergence among the subspecies. The estimated divergence time of the three subspecies was 187,000-226,000 years ago. These results will help us to understand the genetic diversity of wild mice on a global scale, and the findings will be particularly useful in future biomedical and evolutionary studies involving laboratory mice established from such wild mice.