Inferring human neutral genetic variation from craniodental phenotypes.

There is a growing consensus that global patterns of modern human cranial and dental variation are shaped largely by neutral evolutionary processes, suggesting that craniodental features can be used as reliable proxies for inferring population structure and history in bioarchaeological, forensic, and paleoanthropological contexts. However, there is disagreement on whether certain types of data preserve a neutral signature to a greater degree than others. Here, we address this unresolved question and systematically test the relative neutrality of four standard metric and nonmetric craniodental data types employing an extensive computational genotype-phenotype comparison across modern populations from around the world. Our computation draws on the largest existing data sets currently available, while accounting for geographically structured environmental variation, population sampling uncertainty, disparate numbers of phenotypic variables, and stochastic variation inherent to a neutral model of evolution. Our results reveal that the four data types differentially capture neutral genomic variation, with highest signals preserved in dental nonmetric and cranial metric data, followed by cranial nonmetric and dental metric data. Importantly, we demonstrate that combining all four data types together maximizes the neutral genetic signal compared with using them separately, even with a limited number of phenotypic variables. We hypothesize that this reflects a lower level of genetic integration through pleiotropy between, compared to within, the four data types, effectively forming four different modules associated with relatively independent sets of loci. Therefore, we recommend that future craniodental investigations adopt holistic combined data approaches, allowing for more robust inferences about underlying neutral genetic variation.

Dental size variation in admixed Latin Americans: Effects of age, sex and genomic ancestry.

Dental size variation in modern humans has been assessed from regional to worldwide scales, especially under microevolutionary and forensic contexts. Despite this, populations of mixed continental ancestry such as contemporary Latin Americans remain unexplored. In the present study we investigated a large Latin American sample from Colombia (N = 804) and obtained buccolingual and mesiodistal diameters and three indices for maxillary and mandibular teeth (except third molars). We evaluated the correlation between 28 dental measurements (and three indices) with age, sex and genomic ancestry (estimated using genome-wide SNP data). In addition, we explored correlation patterns between dental measurements and the biological affinities, based on these measurements, between two Latin American samples (Colombians and Mexicans) and three putative parental populations: Central and South Native Americans, western Europeans and western Africans through PCA and DFA. Our results indicate that Latin Americans have high dental size diversity, overlapping the variation exhibited by the parental populations. Several dental dimensions and indices have significant correlations with sex and age. Western Europeans presented closer biological affinities with Colombians, and the European genomic ancestry exhibited the highest correlations with tooth size. Correlations between tooth measurements reveal distinct dental modules, as well as a higher integration of postcanine dentition. The effects on dental size of age, sex and genomic ancestry is of relevance for forensic, biohistorical and microevolutionary studies in Latin Americans.

Global patterns of the cranial form of modern human populations described by analysis of a 3D surface homologous model.

This study assessed the regional diversity of the human cranial form by using geometric homologous models based on scanned data from 148 ethnic groups worldwide. This method adopted a template-fitting technique for a nonrigid transformation via the iterative closest point algorithm to generate the homologous meshes. Through the application of principal component analysis to 342 sampled homologous models, the largest variation was detected in overall size, and small South Asian crania were clearly verified. The next greatest diversity was found in the length/breadth proportion of the neurocranium, which showed the contrast between the elongated crania of Africans and the globular crania of Northeast Asians. Notably, this component was slightly correlated with the facial profile. Well-known facial features, such as the forward projection of the cheek among Northeast Asians and compaction of the European maxilla, were reconfirmed. These facial variations were highly correlated with the calvarial outline, particularly the degree of frontal and occipital inclines. An allometric pattern was detected in facial proportions in relation to overall cranial size; in larger crania, the facial profiles tend to be longer and narrower, as demonstrated among many American natives and Northeast Asians. Although our study did not include data on environmental variables that are likely to affect cranial morphology, such as climate or dietary conditions, the large datasets of homologous cranial models will be usefully available for seeking various attributions to phenotypic skeletal characteristics.

Female craniometrics support the 'two-layer model' of human dispersal in Eastern Eurasia.

This study reports a cranio-morphometric analysis of female human remains from seven archaeological sites in China, Vietnam and Taiwan that date between 16,000 and 5300 BP. The aim of the analysis is to test the "two-layer" model of human dispersal in eastern Eurasia, using previously unanalysed female remains to balance the large sample of previously-analysed males. The resulting craniometric data indicate that the examined specimens all belong to the "first layer" of dispersal, and share a common ancestor with recent Australian and Papuan populations, and the ancient Jomon people of Japan. The analysed specimens pre-date the expansion of agricultural populations of East/Northeast Asian origin-that is, the "second layer" of human dispersal proposed by the model. As a result of this study, the two-layer model, which has hitherto rested on evidence only from male skeletons, is now strongly supported by female-derived data. Further comparisons reveal that the people of the first layer were closer in terms of their facial morphology to modern Africans and Sri Lankan Veddah than to modern Asians and Europeans, suggesting that the Late Pleistocene through Middle Holocene hunter-gatherers examined in this study were direct descendants of the anatomically modern humans who first migrated out of Africa through southern Eurasia.

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.

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.

The prehistoric peopling of Southeast Asia.

The human occupation history of Southeast Asia (SEA) remains heavily debated. Current evidence suggests that SEA was occupied by Hòabìnhian hunter-gatherers until ~4000 years ago, when farming economies developed and expanded, restricting foraging groups to remote habitats. Some argue that agricultural development was indigenous; others favor the "two-layer" hypothesis that posits a southward expansion of farmers giving rise to present-day Southeast Asian genetic diversity. By sequencing 26 ancient human genomes (25 from SEA, 1 Japanese Jomon), we show that neither interpretation fits the complexity of Southeast Asian history: Both Hòabìnhian hunter-gatherers and East Asian farmers contributed to current Southeast Asian diversity, with further migrations affecting island SEA and Vietnam. Our results help resolve one of the long-standing controversies in Southeast Asian prehistory.

Reconstructing human population history from dental phenotypes.

Dental phenotypic data are often used to reconstruct biological relatedness among past human groups. Teeth are an important data source because they are generally well preserved in the archaeological and fossil record, even when associated skeletal and DNA preservation is poor. Furthermore, tooth form is considered to be highly heritable and selectively neutral; thus, teeth are assumed to be an excellent proxy for neutral genetic data when none are available. However, to our knowledge, no study to date has systematically tested the assumption of genetic neutrality of dental morphological features on a global scale. Therefore, for the first time, this study quantifies the correlation of biological affinities between worldwide modern human populations, derived independently from dental phenotypes and neutral genetic markers. We show that population relationship measures based on dental morphology are significantly correlated with those based on neutral genetic data (on average r = 0.574, p < 0.001). This relatively strong correlation validates tooth form as a proxy for neutral genomic markers. Nonetheless, we suggest caution in reconstructions of population affinities based on dental data alone because only part of the dental morphological variation among populations can be explained in terms of neutral genetic differences.

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.

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.

Testing modern human out-of-Africa dispersal models and implications for modern human origins.

The modern human expansion process out of Africa has important implications for understanding the genetic and phenotypic structure of extant populations. While intensely debated, the primary hypotheses focus on either a single dispersal or multiple dispersals out of the continent. Here, we use the human fossil record from Africa and the Levant, as well as an exceptionally large dataset of Holocene human crania sampled from Asia, to model ancestor-descendant relationships along hypothetical dispersal routes. We test the spatial and temporal predictions of competing out-of-Africa models by assessing the correlation of geographical distances between populations and measures of population differentiation derived from quantitative cranial phenotype data. Our results support a model in which extant Australo-Melanesians are descendants of an initial dispersal out of Africa by early anatomically modern humans, while all other populations are descendants of a later migration wave. Our results have implications for understanding the complexity of modern human origins and diversity.

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.

Natural allelic variations of xenobiotic-metabolizing enzymes affect sexual dimorphism in Oryzias latipes.

Sexual dimorphisms, which are phenotypic differences between males and females, are driven by sexual selection. Interestingly, sexually selected traits show geographical variations within species despite strong directional selective pressures. This paradox has eluded many evolutionary biologists for some time, and several models have been proposed (e.g. 'indicator model' and 'trade-off model'). However, disentangling which of these theories explains empirical patterns remains difficult, because genetic polymorphisms that cause variation in sexual differences are still unknown. In this study, we show that polymorphisms in cytochrome P450 (CYP) 1B1, which encodes a xenobiotic-metabolizing enzyme, are associated with geographical differences in sexual dimorphism in the anal fin morphology of medaka fish (Oryzias latipes). Biochemical assays and genetic cross experiments show that high- and low-activity CYP1B1 alleles enhanced and declined sex differences in anal fin shapes, respectively. Behavioural and phylogenetic analyses suggest maintenance of the high-activity allele by sexual selection, whereas the low-activity allele possibly has experienced positive selection due to by-product effects of CYP1B1 in inferred ancestral populations. The present data can elucidate evolutionary mechanisms behind genetic variations in sexual dimorphism and indicate trade-off interactions between two distinct mechanisms acting on the two alleles with pleiotropic effects of xenobiotic-metabolizing enzymes.

Genome-wide SNP analysis reveals population structure and demographic history of the ryukyu islanders in the southern part of the Japanese archipelago.

The Ryukyu Islands are located to the southwest of the Japanese archipelago. Archaeological evidence has revealed the existence of prehistoric cultural differentiation between the northern Ryukyu islands of Amami and Okinawa, and the southern Ryukyu islands of Miyako and Yaeyama. To examine a genetic subdivision in the Ryukyu Islands, we conducted genome-wide single nucleotide polymorphism typing of inhabitants from the Okinawa Islands, the Miyako Islands, and the Yaeyama Islands. Principal component and cluster analyses revealed genetic differentiation among the island groups, especially between Okinawa and Miyako. No genetic affinity was observed between aboriginal Taiwanese and any of the Ryukyu populations. The genetic differentiation observed between the inhabitants of the Okinawa Islands and the Miyako Islands is likely to have arisen due to genetic drift rather than admixture with people from neighboring regions. Based on the observed genetic differences, the divergence time between the inhabitants of Okinawa and Miyako islands was dated to the Holocene. These findings suggest that the Pleistocene inhabitants, whose bones have been found on the southern Ryukyu Islands, did not make a major genetic contribution, if any, to the present-day inhabitants of the southern Ryukyu Islands.

The use of tooth crown measurements in the assessment of ancestry.

As the hardest tissue in the body, teeth have the potential to offer a wealth of biological information to the forensic anthropologist, which can include the assessment of ancestry. Using a large data set of dental measurements, the efficacy of mesiodistal and buccolingual tooth dimensions to discriminate between broad, geographically based groups is explored. A general pattern is identified: African populations have the largest teeth, Asians possess teeth of intermediate size, and Europeans have the smallest teeth. In a discriminant function analysis using crown measurements of all teeth (mandibular and maxillary and excluding the third molar), individuals were correctly classified in 71.3% of cases. When the sex of the individual is known, classification is improved up to 88.1% in females and 71.9% of males (cross-validated). Based on these results, we argue that dental metrics can be regularly employed as part of the development of the biological profile.

Comparability of multiple data types from the Bering Strait region: cranial and dental metrics and nonmetrics, mtDNA, and Y-chromosome DNA.

Different data types have previously been shown to have the same microevolutionary patterns in worldwide data sets. However, peopling of the New World studies have shown a difference in migration paths and timings using multiple types of data, spurring research to understand why this is the case. This study was designed to test the degree of similarity in evolutionary patterns by using cranial and dental metric and nonmetric data, along with Y-chromosome DNA and mtDNA. The populations used included Inuits from Alaska, Canada, Siberia, Greenland, and the Aleutian Islands. For comparability, the populations used for the cranial and molecular data were from similar geographic regions or had a shared population history. Distance, R and kinship matrices were generated for use in running Mantel tests, PROTEST analyses, and Procrustes analyses. A clear patterning was seen, with the craniometric data being most highly correlated to the mtDNA data and the cranial nonmetric data being most highly correlated with the Y-chromosome data, while the phenotypic data were also linked. This patterning is suggestive of a possible male or female inheritance, or the correlated data types are affected by the same or similar evolutionary forces. The results of this study indicate cranial traits have some degree of heritability. Moreover, combining data types leads to a richer knowledge of biological affinity. This understanding is important for bioarchaeological contexts, in particular, peopling of the New World studies where focusing on reconciling the results from comparing multiple data types is necessary to move forward.

Generation of monoclonal antibodies against the Galß1-4Gal epitope: a key tool in studies of species-specific glycans expressed in fish, amphibians and birds.

Whereas the Galß1-4Gal epitope is rarely found in mammalian glycans, it has been found in glycans of various species of non-mammalian vertebrates, such as fish, amphibians and birds. Although glycans containing Galß1-4Gal in these vertebrates were detected by precise structural analysis of the glycans using mass spectrometry and/or NMR spectrometry, there are no convenient methods to detect Galß1-4Gal from various samples. To analyze systematically the distribution of Galß1-4Gal in nature, we generated mouse monoclonal antibodies (mAbs) specific for Galß1-4Gal using extracts of medaka eggs as an immunogen. Four mAbs (two immunoglobulin (Ig)Ms and two IgG1s) were obtained by enzyme-linked immunosorbent assay-based screening. The specificities of these mAbs were evaluated by frontal affinity chromatography using 142 kinds of 2-aminopyridine (PA)-derivatized oligosaccharides. While all mAbs interacted with (Galß1-4Gal)-containing oligosaccharides at their non-reducing termini with dissociation constants (K(d)) ranging from 1.0 x 10⁻5 to 2.8 x 10⁻4 M, no apparent interaction was observed with any other glycans. The number of branches containing Galß1-4Gal on N-glycans did not significantly affect K(d) of mAbs of IgG1 subclasses, but those of IgM mAbs were decreased by ∼1 order of magnitude, in increments of the number of branches present. Using the mAbs, we established that Galß1-4Gal is also expressed on glycoproteins in various tissues from the African clawed frog. Immunohistochemical staining of medaka sections revealed that Galß1-4Gal epitopes were expressed in the endothelium, epithelium and epidermis, which directly contact the external environment or invading organisms. Thus, these mAbs are useful for systematically investigating the species-specific expression of glycans, which may act as a barrier against infection.

Effects of an Asian-specific nonsynonymous EDAR variant on multiple dental traits.

Dental morphology is highly diverse among individuals and between human populations. Although it is thought that genetic factors mainly determine common dental variations, only a few such genetic factors have been identified. One study demonstrated that a nonsynonymous single-nucleotide polymorphism (370V/A, rs3827760) in the ectodysplasin A receptor gene (EDAR) is associated with shoveling and double-shoveling grades of upper first incisors and tooth crown size. Here, we examined the association of EDAR 370V/A with several dental characters in Korean and Japanese subjects. A meta-analysis that combined analyses of Korean and Japanese subjects revealed that the Asian-specific 370A allele is associated with an increase in the grades of shoveling and double shoveling, as previously found. We also showed a highly significant association between EDAR 370V/A genotype and crown size, especially mesiodistal diameters of anterior teeth. Moreover, we found that the 370A allele was associated with the presence of hypoconulids of lower second molars. These results indicated that the EDAR polymorphism is responsible, in part, for the Sinodonty and Sundadonty dichotomy in Asian populations, and clearly demonstrated that the EDAR polymorphism has pleiotropic effects on tooth morphology. As the 370A allele is known to be a most likely target of positive selection in Asian populations, some phenotypes associated with the variant may be 'hitchhiking phenotypes', while others may be actual targets of selection.

An Aboriginal Australian genome reveals separate human dispersals into Asia.

We present an Aboriginal Australian genomic sequence obtained from a 100-year-old lock of hair donated by an Aboriginal man from southern Western Australia in the early 20th century. We detect no evidence of European admixture and estimate contamination levels to be below 0.5%. We show that Aboriginal Australians are descendants of an early human dispersal into eastern Asia, possibly 62,000 to 75,000 years ago. This dispersal is separate from the one that gave rise to modern Asians 25,000 to 38,000 years ago. We also find evidence of gene flow between populations of the two dispersal waves prior to the divergence of Native Americans from modern Asian ancestors. Our findings support the hypothesis that present-day Aboriginal Australians descend from the earliest humans to occupy Australia, likely representing one of the oldest continuous populations outside Africa.