A common variation in EDAR is a genetic determinant of shovel-shaped incisors

Genomic imputation of ancient Asian populations contrasts local adaptation in pre- and post-agricultural Japan

Early modern humans lived as hunter-gatherers for millennia before agriculture, yet the genetic adaptations of these populations remain a mystery. Here, we investigate selection in the ancient hunter-gatherer-fisher Jomon and contrast pre- and post-agricultural adaptation in the Japanese archipelago. Building on the successful validation of imputation with ancient Asian genomes, we identify selection signatures in the Jomon, particularly robust signals from KITLG variants, which may have influenced dark pigmentation evolution. The Jomon lacks well-known adaptive variants (EDAR, ADH1B, and ALDH2), marking their emergence after the advent of farming in the archipelago. Notably, the EDAR and ADH1B variants were prevalent in the archipelago 1,300 years ago, whereas the ALDH2 variant could have emerged later due to its absence in other ancient genomes. Overall, our study underpins local adaptation unique to the Jomon population, which in turn sheds light on post-farming selection that continues to shape contemporary Asian populations.

Malaria resistance-related biological adaptation and complex evolutionary footprints inferred from one integrative Tai-Kadai-related genomic resource

Pathogen‒host adaptative interactions and complex population demographical processes, including admixture, drift, and Darwen selection, have considerably shaped the Neolithic-to-Modern Western Eurasian population structure and genetic susceptibility to modern human diseases. However, the genetic footprints of evolutionary events in East Asia remain unknown due to the underrepresentation of genomic diversity and the design of large-scale population studies. We reported one aggregated database of genome-wide SNP variations from 796 Tai-Kadai (TK) genomes, including that of Bouyei first reported here, to explore the genetic history, population structure, and biological adaptative features of TK people from southern China and Southeast Asia. We found geography-related population substructure among TK people using the state-of-the-art population genetic structure reconstruction techniques based on the allele frequency spectrum and haplotype-resolved phased fragments. We found that the northern TK people from Guizhou harbored one TK-dominant ancestry maximized in the Bouyei people, and the southern TK people from Thailand were more influenced by Southeast Asians and indigenous people. We reconstructed fitted admixture models and demographic graphs, which showed that TK people received gene flow from ancient southern rice farmer-related lineages related to the Hmong-Mien and Austroasiatic people and from northern millet farmers associated with the Sino-Tibetan people. Biological adaptation focused on our identified unique TK lineages related to Bouyei, which showed many adaptive signatures conferring Malaria resistance and low-rate lipid metabolism. Further gene enrichment, the allele frequency distribution of derived alleles, and their correlation with the incidence of Malaria further confirmed that CR1 played an essential role in the resistance of Malaria in the ancient "Baiyue" tribes.

Finite element analysis of Neanderthal and early Homo sapiens maxillary central incisor

Neanderthal anterior teeth are very large and have a distinctive morphology characterized by robust 'shovel-shaped' crowns. These features are frequently seen as adaptive responses in dissipating heavy mechanical loads resulting from masticatory and non-masticatory activities. Although the long-standing debate surrounding this hypothesis has played a central role in paleoanthropology, is still unclear if Neanderthal anterior teeth can resist high mechanical loads or not. A novel way to answer this question is to use a multidisciplinary approach that considers together tooth architecture, dental wear and jaw movements. The aim of this study is to functionally reposition the teeth of Le Moustier 1 (a Neanderthal adolescent) and Qafzeh 9 (an early Homo sapiens adolescent) derived from wear facet mapping, occlusal fingerprint analysis and physical dental restoration methods. The restored dental arches are then used to perform finite element analysis on the left central maxillary incisor during edge-to-edge occlusion. The results show stress distribution differences between Le Moustier 1 and Qafzeh 9, with the former displaying higher tensile stress in enamel around the lingual fossa but lower concentration of stress in the lingual aspect of the root surface. These results seem to suggest that the presence of labial convexity, lingual tubercle and of a large root surface in Le Moustier 1 incisor helps in dissipating mechanical stress. The absence of these dental features in Qafzeh 9 is compensated by the presence of a thicker enamel, which helps in reducing the stress in the tooth crown.

Differentiated adaptative genetic architecture and language-related demographical history in South China inferred from 619 genomes from 56 populations

BACKGROUND

The underrepresentation of human genomic resources from Southern Chinese populations limited their health equality in the precision medicine era and complete understanding of their genetic formation, admixture, and adaptive features. Besides, linguistical and genetic evidence supported the controversial hypothesis of their origin processes. One hotspot case was from the Chinese Guangxi Pinghua Han people (GPH), whose language was significantly similar to Southern Chinese dialects but whose uniparental gene pool was phylogenetically associated with the indigenous Tai-Kadai (TK) people. Here, we analyzed genome-wide SNP data in 619 people from four language families and 56 geographically different populations, in which 261 people from 21 geographically distinct populations were first reported here.

RESULTS

We identified significant population stratification among ethnolinguistically diverse Guangxi populations, suggesting their differentiated genetic origin and admixture processes. GPH shared more alleles related to Zhuang than Southern Han Chinese but received more northern ancestry relative to Zhuang. Admixture models and estimates of genetic distances showed that GPH had a close genetic relationship with geographically close TK compared to Northern Han Chinese, supporting their admixture origin hypothesis. Further admixture time and demographic history reconstruction supported GPH was formed via admixture between Northern Han Chinese and Southern TK people. We identified robust signatures associated with lipid metabolisms, such as fatty acid desaturases (FADS) and medically relevant loci associated with Mendelian disorder (GJB2) and complex diseases. We also explored the shared and unique selection signatures of ethnically different but linguistically related Guangxi lineages and found some shared signals related to immune and malaria resistance.

CONCLUSIONS

Our genetic analysis illuminated the language-related fine-scale genetic structure and provided robust genetic evidence to support the admixture hypothesis that can explain the pattern of observed genetic diversity and formation of GPH. This work presented one comprehensive analysis focused on the population history and demographical adaptative process, which provided genetic evidence for personal health management and disease risk prediction models from Guangxi people. Further large-scale whole-genome sequencing projects would provide the entire landscape of southern Chinese genomic diversity and their contributions to human health and disease traits.

Peopling of the Americas: A new approach to assessing dental morphological variation in Asian and Native American populations

OBJECTIVES

Through biodistance analyses, anthropologists have used dental morphology to elucidate how people moved into and throughout the Americas. Here, we apply a method that focuses on individuals rather than sample frequencies through the application rASUDAS2, based on a naïve Bayes' algorithm.

MATERIALS AND METHODS

Using the database of C.G. Turner II, we calculated the probability that an individual could be assigned to one of seven biogeographic groups (American Arctic, North & South America, East Asia, Southeast Asia & Polynesia, Australo-Melanesia, Western Eurasia, & Sub-Saharan Africa) through rASUDAS2. The frequency of classifications for each biogeographic group was determined for 1418 individuals from six regions across Asia and the Americas.

RESULTS

Southeast Asians show mixed assignments but rarely to American Arctic or "American Indian." East Asians are assigned to East Asia half the time while 30% are assigned as Native American. People from the American Arctic and North & South America are assigned to Arctic America or non-Arctic America 75%-80% of the time, with 10%-15% classified as East Asian.

DISCUSSION

All Native American groups have a similar degree of morphological affinity to East Asia, as 10%-15% are classified as East Asian. East Asians are classified as Native American in 30% of cases. Individuals in the Western Hemisphere are decreasingly classified as Arctic the farther south they are located. Equivalent levels of classification as East Asian across all Native American groups suggests one divergence between East Asians and the population ancestral to all Native Americans. Non-arctic Native American groups are derived from the Arctic population, which represents the Native American founder group.

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.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Genetic/Protein Association of Atopic Dermatitis and Tooth Agenesis

Atopic dermatitis and abnormalities in tooth development (including hypomineralization, hypodontia and microdontia) have been observed to co-occur in some patients. A common pathogenesis pathway that involves genes and protein interactions has been hypothesized. This review aims to first provide a description of the key gene mutations and signaling pathways associated with atopic dermatitis and tooth agenesis (i.e., the absence of teeth due to developmental failure) and identify the possible association between the two diseases. Second, utilizing a list of genes most commonly associated with the two diseases, we conducted a protein–protein network interaction analysis using the STRING database and identified a novel association between the Wnt/β-catenin signaling pathway (major pathway responsible for TA) and desmosomal proteins (component of skin barrier that affect the pathogenesis of AD). Further investigation into the mechanisms that may drive their co-occurrence and underlie the development of the two diseases is warranted.

Additional predictions for forensic DNA phenotyping of externally visible characteristics using the ForenSeq and Imagen kits

Multiplex DNA typing methods using massively parallel sequencing can be used to predict externally visible characteristics (EVCs) in forensic DNA phenotyping through the analysis of single-nucleotide polymorphisms. The focus of EVC determination has focused on hair color, eye color, and skin tone as well as visible biogeographical ancestry features. In this study, we researched off-label applications beyond what is currently marketed by the manufacturer of the Verogen ForenSeq kit primer set B and Imagen primer set E SNP loci. We investigated additional EVC predictions by examining published genome wide sequencing studies and reported allele-specific gene expression and predictive values. We have identified 15 SNPs included in the ForenSeq kit panel and Imagen kits that have additional EVC prediction capabilities beyond what is published in the Verogen manuals. The additional EVCs that can be predicted include hair graying, ephelides hyperpigmented spots, dermatoheliosis, facial pigmented spots, standing height, pattern balding, helix-rolling ear morphology, hair shape, hair thickness, facial morphology, eyebrow thickness, sarcoidosis, obesity, vitiligo, and tanning propensity. The loci can be used to augment and refine phenotype predictions with software such as MetaHuman for missing persons, cold case, and historic case investigations.

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.

KDF1 Novel Variant Causes Unique Dental and Oral Epithelial Defects

Keratinocyte differentiation factor 1 (KDF1) is a recently identified and rare candidate gene for human tooth agenesis; however, KDF1-related morphological characteristics and pathological changes in dental tissue and the oral epithelium remain largely unknown. Here, we employed whole-exome sequencing (WES) and Sanger sequencing to screen for the suspected variants in a cohort of 151 tooth agenesis patients, and we segregated a novel KDF1 heterozygous missense variation, c.920G>C (p.R307P), in a non-syndromic tooth agenesis family. Essential bioinformatics analyses and tertiary structural predictions were performed to analyze the structural changes and functional impacts of the novel KDF1 variant. The subsequent functional assessment using a TOP-flash/FOP-flash luciferase reporter system demonstrated that KDF1 variants suppressed the activation of canonical Wnt signaling in 293T cells. To comprehensively investigate the KDF1-related oral morphological anomalies, we performed scanning electron microscopy and ground section of the lower right lateral deciduous incisor extracted from #285 proband, and histopathological assessment of the gingiva. The phenotypic analyses revealed a series of tooth morphological anomalies related to the KDF1 variant R307P, including a shovel-shaped lingual surface of incisors and cornicione-shaped marginal ridges with anomalous morphological occlusal grooves of premolars and molars. Notably, keratinized gingival epithelium abnormalities were revealed in the proband and characterized by epithelial dyskeratosis with residual nuclei, indistinct stratum granulosum, epithelial hyperproliferation, and impaired epithelial differentiation. Our findings revealed new developmental anomalies in the tooth and gingival epithelium of a non-syndromic tooth agenesis individual with a novel pathogenic KDF1 variant, broadening the phenotypic spectrum of KDF1-related disorders and providing new evidence for the crucial role of KDF1 in regulating human dental and oral epithelial development.

Evolutionary adaptation highlights the interconnection of fatty acids, sunlight, inflammation and epithelial adhesion

Gene variants that influence human biology today reflect thousands of years of evolution. Genetic effects on infant health are a major point of selective pressure, given that childhood survival is essential to evolutionary success. Knowledge of this evolutionary history can have implications for paediatric research. CONCLUSION: An episode of human adaptation to the extremely low ultraviolet radiation environment of the Arctic 20,000 years ago implicates the Ectodysplasin A Receptor (EDAR) and the Fatty Acid Desaturases (FADS) in human lactation and epithelial inflammation.

Evidence of selection in the ectodysplasin pathway among endangered aquatic mammals

Synopsis

The ectodysplasin pathway has been a target of evolution repeatedly. Genetic variation in the key genes of this pathway (EDA, EDAR, and EDARADD) results in a rich source of pleiotropic effects across ectodermally-derived structures, including teeth, hair, sweat glands, and mammary glands. In addition, a non-canonical Wnt pathway has a very similar functional role, making variation in the WNT10A gene also of evolutionary significance. The adaptation of mammals to aquatic environments has occurred independently in at least 4 orders, whose species occupy a wide geographic range (from equatorial to polar regions) and exhibit great phenotypic variation in ectodermally-derived structures, including the presence or absence of fur and extreme lactational strategies. The role of the ectodysplasin pathway in the adaptation to aquatic environments has been never explored in mammalian species. In the present study, we analyze the genetic variation in orthologous coding sequences from EDA, EDAR, EDARADD, and WNT10A genes together with ectodermally-derived phenotypic variation from 34 aquatic and non-aquatic mammalian species to assess signals of positive selection, gene-trait coevolution, and genetic convergence. Our study reveals strong evidence of positive selection in a proportion of coding sites in EDA and EDAR genes in 3 endangered aquatic mammals (the Hawaiian monk seal, the Yangtze finless porpoise, and the sea otter). We hypothesize functional implications potentially related to the adaptation to the low-latitude aquatic environment in the Hawaiian monk seal and the freshwater in the Yangtze finless porpoise. The signal in the sea otter is likely the result of an increased genetic drift after an intense bottleneck and reduction of genetic diversity. Besides positive selection, we have not detected robust signals of gene-trait coevolution or convergent amino acid shifts in the ectodysplasin pathway associated with shared phenotypic traits among aquatic mammals. This study provides new evidence of the evolutionary role of the ectodysplasin pathway and encourages further investigation, including functional studies, to fully resolve its relationship with mammalian aquatic adaptation.

Spanish

La vía de la ectodisplasina ha sido objeto de la evolución repetidamente. La variación genética en los principales genes de esta vía (EDA, EDAR y EDARADD) da como resultado una gran diversidad de efectos pleiotrópicos en las estructuras derivadas del ectodermo, incluidos los dientes, el cabello, las glándulas sudoríparas y las glándulas mamarias. Además, una vía wnt no canónica tiene un papel funcional muy similar, por lo que la variación en el gen WNT10A también tiene importancia evolutiva. La adaptación de los mamíferos a los entornes acuáticos se ha producido de forma independiente en al menos cuatro órdenes, cuyas especies ocupan un amplio rango geográfico (desde regiones ecuatoriales a polares) y presentan una gran variación fenotípica en las estructuras derivadas del ectodermo, incluyendo la presencia o ausencia de pelaje y estrategias de lactancia muy diferentes. El papel de la vía de la ectodisplasina en la adaptación a entornos acuáticos no se ha explorado nunca en especies de mamíferos. En este estudio, analizamos la variación genética en las secuencias codificantes ortólogas de los genes EDA, EDAR, EDARADD y WNT10A junto con la variación fenotípica derivada del ectodermo de 34 especies de mamíferos acuáticos y no acuáticos para evaluar señales de selección positiva, coevolución gen-rasgo y convergencia genética. Nuestro estudio revela señales de selección positiva en regiones de las secuencias codificantes de los genes EDA y EDAR en tres mamíferos acuáticos en peligro de extinción (la foca monje de Hawái, la marsopa lisa y la nutria marina). Estas señales podrían tener implicaciones funcionales potencialmente relacionadas con la adaptación al entorno acuático de baja latitud en la foca monje de Hawái y el agua dulce en la marsopa lisa. La señal en la nutria marina es probablemente el resultado de una mayor deriva genética tras un intenso un cuello de botella y una reducción de la diversidad genética. A parte de selección positiva, no hemos detectado señales sólidas de coevolución gen-rasgo o cambios convergentes de aminoácidos en la vía de la ectodisplasina asociados a rasgos fenotípicos compartidos entre mamíferos acuáticos. Este estudio proporciona nuevas evidencias del papel evolutivo de la vía de la ectodisplasina y quiere promover futuras investigaciones con estudios funcionales para acabar de resolver la relación de esta vía con la adaptación acuática de los mamíferos.

Similarity-Based Analysis of Allele Frequency Distribution among Multiple Populations Identifies Adaptive Genomic Structural Variants

Structural variants have a considerable impact on human genomic diversity. However, their evolutionary history remains mostly unexplored. Here, we developed a new method to identify potentially adaptive structural variants based on a similarity-based analysis that incorporates genotype frequency data from 26 populations simultaneously. Using this method, we analyzed 57,629 structural variants and identified 576 structural variants that show unusual population differentiation. Of these putatively adaptive structural variants, we further showed that 24 variants are multiallelic and overlap with coding sequences, and 20 variants are significantly associated with GWAS traits. Closer inspection of the haplotypic variation associated with these putatively adaptive and functional structural variants reveals deviations from neutral expectations due to: 1) population differentiation of rapidly evolving multiallelic variants, 2) incomplete sweeps, and 3) recent population-specific negative selection. Overall, our study provides new methodological insights, documents hundreds of putatively adaptive variants, and introduces evolutionary models that may better explain the complex evolution of structural variants.

Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.

Biometry of width between labial transitional line angles in anterior teeth: an observational study

PURPOSE

The maximum width between the mesial and distal labial transitional line angles, described as “esthetic width” herein, could significantly influence the visual perception of the teeth and smile. This study aimed to conduct biometric research on esthetic width and to explore whether regular distribution exists in the esthetic width of human teeth.

MATERIALS AND METHODS

A total of 4,264 maxillary and mandibular anterior teeth were measured using the Geomagic studio software program. The proportions of maxillary to mandibular homonymous teeth and proportions between the adjacent teeth were calculated. Bilateral symmetry and the correlation between the esthetic and mesiodistal widths were both accounted for during the measurement procedures.

RESULTS

The mean esthetic widths were 6.773 ± 0.518 mm and 4.329 ± 0.331 mm for maxillary and mandibular central incisors, respectively, 5.451 ± 0.487 mm and 5.008 ± 0.351 mm for maxillary and mandibular lateral incisors, respectively, and 3.340 ± 0.353 mm and 5.958 ± 0.415 mm for maxillary and mandibular canines, respectively. Except for the mandibular canines, no significant difference in esthetic width was found among homonymous teeth from the same jaw. A high linear correlation was found between the esthetic and mesiodistal widths of the same tooth, except for the maxillary canines. Esthetic width proportions among different tooth categories showed some regular patterns, which were similar to those of the mesiodistal width.

CONCLUSION

Esthetic width is regularly distributed among the teeth in the Chinese population. This could provide an important reference for anterior dental restorations and dimension recovery in esthetic reconstruction of anterior teeth.

Prioritising positively selected variants in whole-genome sequencing data using FineMAV

Background

In population genomics, polymorphisms that are highly differentiated between geographically separated populations are often suggestive of Darwinian positive selection. Genomic scans have highlighted several such regions in African and non-African populations, but only a handful of these have functional data that clearly associates candidate variations driving the selection process. Fine-Mapping of Adaptive Variation (FineMAV) was developed to address this in a high-throughput manner using population based whole-genome sequences generated by the 1000 Genomes Project. It pinpoints positively selected genetic variants in sequencing data by prioritizing high frequency, population-specific and functional derived alleles.

Results

We developed a stand-alone software that implements the FineMAV statistic. To graphically visualise the FineMAV scores, it outputs the statistics as bigWig files, which is a common file format supported by many genome browsers. It is available as a command-line and graphical user interface. The software was tested by replicating the FineMAV scores obtained using 1000 Genomes Project African, European, East and South Asian populations and subsequently applied to whole-genome sequencing datasets from Singapore and China to highlight population specific variants that can be subsequently modelled. The software tool is publicly available at https://github.com/fadilla-wahyudi/finemav.

Conclusions

The software tool described here determines genome-wide FineMAV scores, using low or high-coverage whole-genome sequencing datasets, that can be used to prioritize a list of population specific, highly differentiated candidate variants for in vitro or in vivo functional screens. The tool displays these scores on the human genome browsers for easy visualisation, annotation and comparison between different genomic regions in worldwide human populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04506-9.

Variation in Middle Stone Age mandibular molar enamel-dentine junction topography at Klasies River Main Site assessed by diffeomorphic surface matching

The morphology and variability of the Middle Stone Age (MSA) hominin fossils from Klasies River Main Site have been the focus of investigation for more than four decades. The mandibular remains have figured prominently in discussions relating to robusticity, size dimorphism, and symphyseal morphology. Variation in corpus size between the robust SAM-AP 6223 and the diminutive SAM-AP 6225 mandibles is particularly impressive, and the difference between the buccolingual diameters of their M₂s significantly exceeds recent human sample variation. SAM-AP 6223 and SAM-AP 6225 are the only Klasies specimens with homologous teeth (M₂ and M₃) that permit comparisons of crown morphology. While the differences in dental trait expression at the outer enamel surfaces of these molars are slight, diffeomorphic surface analyses of their underlying enamel-dentine junction (EDJ) topographies reveal differences that are well beyond the means of pairwise differences among comparative samples of Later Stone Age (LSA) Khoesan and recent African homologues. The EDJs of both SAM-AP 6225 molars and the SAM-AP 6223 M₃ fall outside the envelopes that define the morphospace of these two samples. Although the radiocarbon dated LSA individuals examined here differ by a maximum of some 7000 years, and the two Klasies jaws may differ by perhaps as much as 18,000 years, it is difficult to ascribe their differences to time alone. With reference to the morphoscopic traits by which the SAM-AP 6223 and SAM-AP 6225 EDJs differ, the most striking is the expression of the protoconid cingulum. This is very weakly developed on the SAM-AP 6223 molars and distinct in SAM-AP 6225. As such, this diminutive fossil exhibits a more pronounced manifestation of what is likely a plesiomorphic feature, thus adding to the morphological mosaicism that is evident in the Klasies hominin assemblage. Several possible explanations for the variation and mosaicism in this MSA sample are discussed.

Association of EDARV370A with breast density and metabolic syndrome in Latinos

The ectodysplasin receptor (EDAR) is a tumor necrosis factor receptor (TNF) superfamily member. A substitution in an exon of EDAR at position 370 (EDARV370A) creates a gain of function mutant present at high frequencies in Asian and Indigenous American populations but absent in others. Its frequency is intermediate in populations of Mexican ancestry. EDAR regulates the development of ectodermal tissues, including mammary ducts. Obesity and type 2 diabetes mellitus are prevalent in people with Indigenous and Latino ancestry. Latino patients also have altered prevalence and presentation of breast cancer. It is unknown whether EDARV370A might connect these phenomena. The goals of this study were to determine 1) whether EDARV370A is associated with metabolic phenotypes and 2) if there is altered breast anatomy in women carrying EDARV370A. Participants were from two Latino cohorts, the Arizona Insulin Resistance (AIR) registry and Sangre por Salud (SPS) biobank. The frequency of EDARV370A was 47% in the Latino cohorts. In the AIR registry, carriers of EDARV370A (GG homozygous) had significantly (p < 0.05) higher plasma triglycerides, VLDL, ALT, 2-hour post-challenge glucose, and a higher prevalence of prediabetes/diabetes. In a subset of the AIR registry, serum levels of ectodysplasin A2 (EDA-A2) also were associated with HbA1c and prediabetes (p < 0.05). For the SPS biobank, participants that were carriers of EDARV370A had lower breast density and higher HbA1c (both p < 0.05). The significant associations with measures of glycemia remained when the cohorts were combined. We conclude that EDARV370A is associated with characteristics of the metabolic syndrome and breast density in Latinos.

The need to incorporate human variation and evolutionary theory in forensic anthropology: A call for reform

In 1992, Norm Sauer called for a language shift in which practitioners would move away from the socially loaded term "race" and replace it with the less provocative term "ancestry." While many heeded the call and moved towards ancestry in their research and reports, the actual approach to research and analysis did not change. In response to this change, there was a large growth in ancestry estimation method development in the early decade of the 2000s. However, the practice of ancestry estimation did not adequately incorporate evolutionary theory in interpretation or trait selection and continued with little critical reflection. In the past decade, there has been an increase in ancestry validation methods with little critique of the "race" concept or discussion of modern human variation or reference samples. To advance, forensic anthropologists need to reckon with the practice of ancestry estimation as it is currently practiced. We are calling for another reform in the axiom focusing on evolutionary theory, population history, trait selection, and population-level reference samples. The practice needs to abandon the terms ancestry and race completely and recalibrate to an analysis of population affinity. Population affinity is a statistical approach based on the underlying population structure that would allow the understanding of how microevolutionary forces act in concert with historical events (e.g., colonization, the Transatlantic Slave Trade, etc.) to shape modern human variation. This is not to be confused with geographic ancestry that all too often can be perceived as interchangeable with social race and as an affirmation of the biological concept of race. It is time to critically evaluate the social and scientific implications of the current practice of ancestry estimation, and re-frame our approach to studying and analyzing modern human variation through a population structure approach.

The genetic architecture of anterior tooth morphology in a longitudinal sample of Australian twins and families

OBJECTIVE

This study presents a quantitative genetic analysis of human anterior dental morphology in a longitudinal sample of known genealogy. The primary aim of this work is to generate a suite of genetic correlations within and between deciduous and permanent characters to access patterns of integration across the diphyodont dental complex.

DESIGN

Data were recorded from casted tooth crowns representing participants of a long-term Australian twin and family study (deciduous n = 290, permanent n = 339). Morphological trait expression was observed and scored following Arizona State University Dental Anthropology System standards. Bivariate genetic correlations were estimated using maximum likelihood variance decomposition models in SOLAR v.8.1.1.

RESULTS

Genetic correlation estimates indicate high levels of integration between antimeres but low to moderate levels among traits within a tooth row. Only 9% of deciduous model comparisons were significant, while pleiotropy was indicated for one third of permanent trait pairs. Canine characters stood out as strongly integrated, especially in the deciduous dentition. For homologous characters across dentitions (e.g., deciduous i¹ shoveling and permanent I¹ shoveling), ∼70% of model comparisons yielded significant genetic correlations.

CONCLUSIONS

Patterns of genetic correlation suggest a morphological canine module that spans the primary and secondary dentition. Results also point to the existence of a genetic mechanism conserving morphology across the diphyodont dental complex, such that paired deciduous and permanent traits are more strongly integrated than characters within individual tooth rows/teeth.

Methods in Genetic Analysis for Evaluation Mandibular Shape and Size Variations in Human Mandible

ABSTRACT

The human mandible has been investigated from both clinical and evolutionary perspectives. Recent advances in genome science have identified the genetic regulation of human mandibular shape and size. Identification of genes that regulate mandibular shape and size would not only enhance our understanding of the mechanisms of mandibular growth and development but also help define a strategy to prevent mandibular dysplasia. This review provides a comprehensive summary of why and how the mandible was evaluated in the human mandible genome study. The variation in human mandibular shape and size has been progressively clarified, not only by focusing on the mandible alone but also by using extremely diverse approaches. The methods of data acquisition for evaluating human mandibular shape and size variation are well established. Furthermore, this review explains how to proceed with future research.

The human EDAR 370V/A polymorphism affects tooth root morphology potentially through the modification of a reaction-diffusion system

Morphological variations in human teeth have long been recognized and, in particular, the spatial and temporal distribution of two patterns of dental features in Asia, i.e., Sinodonty and Sundadonty, have contributed to our understanding of the human migration history. However, the molecular mechanisms underlying such dental variations have not yet been completely elucidated. Recent studies have clarified that a nonsynonymous variant in the ectodysplasin A receptor gene (EDAR 370V/A; rs3827760) contributes to crown traits related to Sinodonty. In this study, we examined the association between the EDAR polymorphism and tooth root traits by using computed tomography images and identified that the effects of the EDAR variant on the number and shape of roots differed depending on the tooth type. In addition, to better understand tooth root morphogenesis, a computational analysis for patterns of tooth roots was performed, assuming a reaction-diffusion system. The computational study suggested that the complicated effects of the EDAR polymorphism could be explained when it is considered that EDAR modifies the syntheses of multiple related molecules working in the reaction-diffusion dynamics. In this study, we shed light on the molecular mechanisms of tooth root morphogenesis, which are less understood in comparison to those of tooth crown morphogenesis.

Scanning the human genome for "signatures" of positive selection: Transformative opportunities and ethical obligations

The relationship history of evolutionary anthropology and genetics is complex. At best, genetics is a beautifully integrative part of the discipline. Yet this integration has also been fraught, with punctuated, disruptive challenges to dogma, periodic reluctance by some members of the field to embrace results from analyses of genetic data, and occasional over-assertions of genetic definitiveness by geneticists. At worst, evolutionary genetics has been a tool for reinforcing racism and colonialism. While a number of genetics/genomics papers have disproportionately impacted evolutionary anthropology, here we highlight the 2002 presentation of an elegantly powerful approach for identifying "signatures" of past positive selection from haplotype-based patterns of genetic variation. Together with technological advances in genotyping methods, this article transformed our field by facilitating genome-wide "scans" for signatures of past positive selection in human populations. This approach helped researchers test longstanding evolutionary anthropology hypotheses while simultaneously providing opportunities to develop entirely new ones. Genome-wide scans for signatures of positive selection have since been conducted in diverse worldwide populations, with striking findings of local adaptation and convergent evolution. Yet there are ethical considerations with respect to the ubiquity of these studies and the cross-application of the genome-wide scan approach to existing datasets, which we also discuss.

Characterisation of a second gain of function EDAR variant, encoding EDAR380R, in East Asia

Ectodysplasin A1 receptor (EDAR) is a TNF receptor family member with roles in the development and growth of hair, teeth and glands. A derived allele of EDAR, single-nucleotide variant rs3827760, encodes EDAR:p.(Val370Ala), a receptor with more potent signalling effects than the ancestral EDAR370Val. This allele of rs3827760 is at very high frequency in modern East Asian and Native American populations as a result of ancient positive selection and has been associated with straighter, thicker hair fibres, alteration of tooth and ear shape, reduced chin protrusion and increased fingertip sweat gland density. Here we report the characterisation of another SNV in EDAR, rs146567337, encoding EDAR:p.(Ser380Arg). The derived allele of this SNV is at its highest global frequency, of up to 5%, in populations of southern China, Vietnam, the Philippines, Malaysia and Indonesia. Using haplotype analyses, we find that the rs3827760 and rs146567337 SNVs arose on distinct haplotypes and that rs146567337 does not show the same signs of positive selection as rs3827760. From functional studies in cultured cells, we find that EDAR:p.(Ser380Arg) displays increased EDAR signalling output, at a similar level to that of EDAR:p.(Val370Ala). The existence of a second SNV with partly overlapping geographic distribution, the same in vitro functional effect and similar evolutionary age as the derived allele of rs3827760, but of independent origin and not exhibiting the same signs of strong selection, suggests a northern focus of positive selection on EDAR function in East Asia.

Ancient genomes from northern China suggest links between subsistence changes and human migration

Northern China harbored the world’s earliest complex societies based on millet farming, in two major centers in the Yellow (YR) and West Liao (WLR) River basins. Until now, their genetic histories have remained largely unknown. Here we present 55 ancient genomes dating to 7500-1700 BP from the YR, WLR, and Amur River (AR) regions. Contrary to the genetic stability in the AR, the YR and WLR genetic profiles substantially changed over time. The YR populations show a monotonic increase over time in their genetic affinity with present-day southern Chinese and Southeast Asians. In the WLR, intensification of farming in the Late Neolithic is correlated with increased YR affinity while the inclusion of a pastoral economy in the Bronze Age was correlated with increased AR affinity. Our results suggest a link between changes in subsistence strategy and human migration, and fuel the debate about archaeolinguistic signatures of past human migration.

Northern China contains some of the world’s earliest farming societies. Here, authors use 55 ancient genomes to trace the genetic history of human migrations across northern China for the last 7500 years, and document genetic changes mirroring shifts in subsistence strategy.

Testing the utility of dental morphological trait combinations for inferring human neutral genetic variation

Researchers commonly rely on human dental morphological features in order to reconstruct genetic affinities among past individuals and populations, particularly since teeth are often the best preserved part of a human skeleton. Tooth form is considered to be highly heritable and selectively neutral and, therefore, to be an excellent proxy for DNA when none is available. However, until today, it remains poorly understood whether certain dental traits or trait combinations preserve neutral genomic signatures to a greater degree than others. Here, we address this long-standing research gap by systematically testing the utility of 27 common dental traits and >134 million possible trait combinations in reflecting neutral genomic variation in a worldwide sample of modern human populations. Our analyses reveal that not all traits are equally well-suited for reconstructing population affinities. Whereas some traits largely reflect neutral variation and therefore evolved primarily as a result of genetic drift, others can be linked to nonstochastic processes such as natural selection or hominin admixture. We also demonstrate that reconstructions of population affinity based on many traits are not necessarily more reliable than those based on only a few traits. Importantly, we find a set of highly diagnostic trait combinations that preserve neutral genetic signals best (up to [Formula: see text] _r = 0.580; 95% r range = 0.293 to 0.758; P = 0.001). We propose that these trait combinations should be prioritized in future research, as they allow for more accurate inferences about past human population dynamics when using dental morphology as a proxy for DNA.

Prevalence of malocclusion in Chinese schoolchildren from 1991 to 2018: A systematic review and meta-analysis

BACKGROUND

Malocclusion is a common oral health problem in schoolchildren. Literature describing the prevalence of malocclusion varies substantially across China.

AIM

This study identified the epidemiological characteristics of malocclusion among Chinese schoolchildren from 1991 to 2018.

DESIGN

Six English and Chinese electronic databases were searched through November 2018. The search was supplemented by hand searching to identify relevant surveys. The overall prevalence of malocclusion was estimated by a random-effects meta-analysis model, and variations in different groups were assessed by subgroup meta-analysis.

RESULTS

Thirty-seven eligible articles describing 117 682 samples were investigated. The pooled national prevalence for malocclusion was 47.92% (95% CI: 58.6%-71.9%). For the Angle classification, the overall prevalence rates were 30.07% (95% CI: 25.37%-35.48%), 9.91% (95% CI: 7.41%-13.79%), and 4.76% (95% CI: 3.85%-6.54%) for Class I, Class II, and Class III malocclusion, respectively. A deep overbite (16.67%, 95% CI: 11.50%-23.08%) was shown to be the most common trait of malocclusion. When stratified by sex, males had a slightly higher prevalence than females (RR = 1.04, 95% CI: 1.01-1.06). More importantly, an ascending trend and substantial variations across the country were observed.

CONCLUSIONS

Our results confirmed that malocclusion has become a serious oral health problem in Chinese schoolchildren, highlighting the need for proactive interventions at an early age. Moreover, high-quality epidemiological studies on malocclusion are still required.

Novel genetic loci affecting facial shape variation in humans

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

Tooth dimensions and body size in a Pygmy population

Background: The relationship between tooth size and stature has been analysed extensively at the interspecies level but has received less attention at the intraspecies level. The relationship between these two parameters does not seem to be the same among modern human populations.Aim: The aim of this study is to analyse the relationship between tooth dimensions and body measurements in the Baka Pygmies.Subjects and methods: Height, weight, and tooth dimensions were obtained for 45 adult Baka females and 17 males from Le Bosquet (Cameroon). Correlations were obtained between the variables and compared to results for other human populations.Results: The Baka population is distinctive in the small number of significant correlations. Only two buccolingual diameters among Baka females show any significant correlation with height. The lack of significant correlations between tooth dimensions and body dimensions among the Baka means that changes in body size are accompanied by random variations in tooth dimensions.Conclusion: The absence of correlations may be accounted for by the impact of environmental effects on the somatic growth of the Baka producing a Pygmy phenotype adapted to live in the forest. It is worth noting that many correlations become significant when sexes are pooled.

An approximate full-likelihood method for inferring selection and allele frequency trajectories from DNA sequence data

Most current methods for detecting natural selection from DNA sequence data are limited in that they are either based on summary statistics or a composite likelihood, and as a consequence, do not make full use of the information available in DNA sequence data. We here present a new importance sampling approach for approximating the full likelihood function for the selection coefficient. Our method CLUES treats the ancestral recombination graph (ARG) as a latent variable that is integrated out using previously published Markov Chain Monte Carlo (MCMC) methods. The method can be used for detecting selection, estimating selection coefficients, testing models of changes in the strength of selection, estimating the time of the start of a selective sweep, and for inferring the allele frequency trajectory of a selected or neutral allele. We perform extensive simulations to evaluate the method and show that it uniformly improves power to detect selection compared to current popular methods such as nSL and SDS, and can provide reliable inferences of allele frequency trajectories under many conditions. We also explore the potential of our method to detect extremely recent changes in the strength of selection. We use the method to infer the past allele frequency trajectory for a lactase persistence SNP (MCM6) in Europeans. We also infer the trajectory of a SNP (EDAR) in Han Chinese, finding evidence that this allele's age is much older than previously claimed. We also study a set of 11 pigmentation-associated variants. Several genes show evidence of strong selection particularly within the last 5,000 years, including ASIP, KITLG, and TYR. However, selection on OCA2/HERC2 seems to be much older and, in contrast to previous claims, we find no evidence of selection on TYRP1.

Biomedical and genetic characteristics of the Ryukyuans: demographic history, diseases and physical and physiological traits

Context: The Ryukyu Islands stretch across a southwestern area of the Japanese Archipelago. Because of their unique geographical and historical backgrounds, Ryukyuans have their own genetic and phenotypic characteristics, which have been disclosed in previous anthropological and biomedical studies. Objective: The history, peopling and biomedical and genetic characteristics of Ryukyuans are reviewed and future research directions are discussed. Conclusion: Morphological and genetic studies have suggested the complex demographic history of Ryukyuans and their relationships with other Asian populations. Knowledge of population formation processes is important to understand the distribution of pathogens. In viral infectious diseases, some strains that may be associated with disease symptoms are specific to Ryukyuans. Dramatic changes in diet have played an important role among Ryukyuans in terms of increases in lifestyle-related diseases and mortality risks. To achieve a better understanding of pathogenic disease factors, further integration of findings regarding the genetic and biomedical characteristics of the Ryukyuans is needed.

EDAR, LYPLAL1, PRDM16, PAX3, DKK1, TNFSF12, CACNA2D3, and SUPT3H gene variants influence facial morphology in a Eurasian population

In human society, the facial surface is visible and recognizable based on the facial shape variation which represents a set of highly polygenic and correlated complex traits. Understanding the genetic basis underlying facial shape traits has important implications in population genetics, developmental biology, and forensic science. A number of single nucleotide polymorphisms (SNPs) are associated with human facial shape variation, mostly in European populations. To bridge the gap between European and Asian populations in term of the genetic basis of facial shape variation, we examined the effect of these SNPs in a European-Asian admixed Eurasian population which included a total of 612 individuals. The coordinates of 17 facial landmarks were derived from high resolution 3dMD facial images, and 136 Euclidean distances between all pairs of landmarks were quantitatively derived. DNA samples were genotyped using the Illumina Infinium Global Screening Array and imputed using the 1000 Genomes reference panel. Genetic association between 125 previously reported facial shape-associated SNPs and 136 facial shape phenotypes was tested using linear regression. As a result, a total of eight SNPs from different loci demonstrated significant association with one or more facial shape traits after adjusting for multiple testing (significance threshold p < 1.28 × 10-3), together explaining up to 6.47% of sex-, age-, and BMI-adjusted facial phenotype variance. These included EDAR rs3827760, LYPLAL1 rs5781117, PRDM16 rs4648379, PAX3 rs7559271, DKK1 rs1194708, TNFSF12 rs80067372, CACNA2D3 rs56063440, and SUPT3H rs227833. Notably, the EDAR rs3827760 and LYPLAL1 rs5781117 SNPs displayed significant association with eight and seven facial phenotypes, respectively (2.39 × 10-5 < p < 1.28 × 10-3). The majority of these SNPs showed a distinct allele frequency between European and East Asian reference panels from the 1000 Genomes Project. These results showed the details of above eight genes influence facial shape variation in a Eurasian population.

Developmental dental anomalies and their potential role in establishing identity in post-mortem cases: a review

Teeth may provide useful forensic evidence owing to features like uniqueness, stability and comparability. Moreover, the human dentition is heterodont, i.e. all the teeth have different morphology - incisors, canines, premolars and molars. There are sometimes deviations from normal morphology, such as the presence of extra teeth, variation in their shape and size eg the presence of an extra cusp, fractured crown/root, Carabelli's cusp, peg laterals, transpositions, fusion, etc. These differences can help forensic personnel identify bodies, especially where other methods of identification like facial features, fingerprints or DNA typing cannot yield satisfactory results as in cases of badly decomposed bodies, burnt remains, mass disasters, etc. Identification from dentition is based on the direct comparison of post-mortem dental profiles with ante-mortem dental records of the deceased. This article aims to review these developmental and morphological dental traits and their role in post-mortem identification.

Assessment of Correlation of Growth Hormone Receptor Gene with Tooth Dimensions: A CBCT and Genotyping Study

Aim:

Assessment of root morphology, size, and amount of bone around tooth is essential before starting the orthodontic treatment. The present study aimed to assess the relationship between tooth dimensions with two gene variants of growth hormone (GH), namely rs6184 and rs6180.

Materials and Methods:

This study was conducted on 218 subjects (males: 104, females: 114) requiring orthodontic treatment. All underwent cone beam computed tomography (CBCT) scan for orthodontic treatment planning with Kodak CBCT machine. In all teeth, crown height (CH), root length (RL), and crown–root ratio were evaluated. Two growth hormone receptor (GHR) variants (rs6184 and rs6180) were genotyped using the TaqMan genotyping assay.

Results:

The mean CH and RL of all teeth, that is, maxillary and mandibular central incisors, lateral incisor, canine, first premolar, second premolar, first molar, and second molar, were measured. There was no significant difference in males and females (P > 0.05). Allele frequencies of GHR variants for rs6180 and rs6184 were 48.1% and 8.92%, respectively. Multiple regression analysis showed GHR rs6184 association with maxillary central incisor CH, maxillary canine RL, mandibular canine CH, and mandibular first premolar RL (P < 0.05).

Conclusion:

There was correlation of CH of maxillary and mandibular canine and RL of maxillary canine and mandibular first premolar with GHR rs6184.

Identification of a Dental Anomaly in a Tiriyo Indigenous Park Patient by Computed Tomography: A Case Report

This report describes anatomical variations in an indigenous patient from the Brazilian Amazon. A 13-year-old indigenous girl attended the dental clinic for a routine examination. Clinically, a change in the coronary morphology of all upper incisors was observed; characterized by a shovel-shaped lingual surface-a feature considered a polygenic hereditary trait commonly found in native American people. The x-ray examination revealed the presence of a root anomaly in the left upper central incisor. A cone-beam computed tomography (CBCT) scan was performed, revealing the presence of a supernumerary root located on the lingual surface. A single wide canal, which bifurcated in the middle-third level into two canals with different foramina, was observed in the cervical portion. It is essential for dental surgeons to be aware of possible anatomical differences, especially considering the origin of the patient, to avoid interference in treatment success.

Ancient Fennoscandian genomes reveal origin and spread of Siberian ancestry in Europe

European population history has been shaped by migrations of people, and their subsequent admixture. Recently, ancient DNA has brought new insights into European migration events linked to the advent of agriculture, and possibly to the spread of Indo-European languages. However, little is known about the ancient population history of north-eastern Europe, in particular about populations speaking Uralic languages, such as Finns and Saami. Here we analyse ancient genomic data from 11 individuals from Finland and north-western Russia. We show that the genetic makeup of northern Europe was shaped by migrations from Siberia that began at least 3500 years ago. This Siberian ancestry was subsequently admixed into many modern populations in the region, particularly into populations speaking Uralic languages today. Additionally, we show that ancestors of modern Saami inhabited a larger territory during the Iron Age, which adds to the historical and linguistic information about the population history of Finland.

A geometric morphometric approach to the study of variation of shovel-shaped incisors

OBJECTIVES

The scoring and analysis of dental nonmetric traits are predominantly accomplished by using the Arizona State University Dental Anthropology System (ASUDAS), a standard protocol based on strict definitions and three-dimensional dental plaques. However, visual scoring, even when controlled by strict definitions of features, visual reference, and the experience of the observer, includes an unavoidable part of subjectivity. In this methodological contribution, we propose a new quantitative geometric morphometric approach to quickly and efficiently assess the variation of shoveling in modern human maxillary central incisors (UI1).

MATERIALS AND METHODS

We analyzed 87 modern human UI1s by means of virtual imaging and the ASU-UI1 dental plaque grades using geometric morphometrics by placing semilandmarks on the labial crown aspect. The modern human sample was composed of individuals from Europe, Africa, and Asia and included representatives of all seven grades defined by the ASUDAS method.

RESULTS

Our results highlighted some limitations in the use of the current UI1 ASUDAS plaque, indicating that it did not necessarily represent an objective gradient of expression of a nonmetric tooth feature. Rating of shoveling tended to be more prone to intra- and interobserver bias for the highest grades. In addition, our analyses suggest that the observers were strongly influenced by the depth of the lingual crown aspect when assessing the shoveling.

DISCUSSION

In this context, our results provide a reliable and reproducible framework reinforced by statistical results supporting the fact that open scale numerical measurements can complement the ASUDAS method.

Towards broadening Forensic DNA Phenotyping beyond pigmentation: Improving the prediction of head hair shape from DNA

Human head hair shape, commonly classified as straight, wavy, curly or frizzy, is an attractive target for Forensic DNA Phenotyping and other applications of human appearance prediction from DNA such as in paleogenetics. The genetic knowledge underlying head hair shape variation was recently improved by the outcome of a series of genome-wide association and replication studies in a total of 26,964 subjects, highlighting 12 loci of which 8 were novel and introducing a prediction model for Europeans based on 14 SNPs. In the present study, we evaluated the capacity of DNA-based head hair shape prediction by investigating an extended set of candidate SNP predictors and by using an independent set of samples for model validation. Prediction model building was carried out in 9674 subjects (6068 from Europe, 2899 from Asia and 707 of admixed European and Asian ancestries), used previously, by considering a novel list of 90 candidate SNPs. For model validation, genotype and phenotype data were newly collected in 2415 independent subjects (2138 Europeans and 277 non-Europeans) by applying two targeted massively parallel sequencing platforms, Ion Torrent PGM and MiSeq, or the MassARRAY platform. A binomial model was developed to predict straight vs. non-straight hair based on 32 SNPs from 26 genetic loci we identified as significantly contributing to the model. This model achieved prediction accuracies, expressed as AUC, of 0.664 in Europeans and 0.789 in non-Europeans; the statistically significant difference was explained mostly by the effect of one EDAR SNP in non-Europeans. Considering sex and age, in addition to the SNPs, slightly and insignificantly increased the prediction accuracies (AUC of 0.680 and 0.800, respectively). Based on the sample size and candidate DNA markers investigated, this study provides the most robust, validated, and accurate statistical prediction models and SNP predictor marker sets currently available for predicting head hair shape from DNA, providing the next step towards broadening Forensic DNA Phenotyping beyond pigmentation traits.

The transcription factor NKX2-3 mediates p21 expression and ectodysplasin-A signaling in the enamel knot for cusp formation in tooth development

Tooth morphogenesis is initiated by reciprocal interactions between the ectoderm and neural crest-derived mesenchyme. During tooth development, tooth cusps are regulated by precise control of proliferation of cell clusters, termed enamel knots, that are present among dental epithelial cells. The interaction of ectodysplasin-A (EDA) with its receptor, EDAR, plays a critical role in cusp formation by these enamel knots, and mutations of these genes is a cause of ectodermal dysplasia. It has also been reported that deficiency in Nkx2-3, encoding a member of the NK2 homeobox family of transcription factors, leads to cusp absence in affected teeth. However, the molecular role of NKX2-3 in tooth morphogenesis is not clearly understood. Using gene microarray analysis in mouse embryos, we found that Nkx2-3 is highly expressed during tooth development and increased during the tooth morphogenesis, especially during cusp formation. We also demonstrate that NKX2-3 is a target molecule of EDA and critical for expression of the cell cycle regulator p21 in the enamel knot. Moreover, NKX2-3 activated the bone morphogenetic protein (BMP) signaling pathway by up-regulating expression levels of Bmp2 and Bmpr2 in dental epithelium and decreased the expression of the dental epithelial stem cell marker SRY box 2 (SOX2). Together, our results indicate that EDA/NKX2-3 signaling is essential for enamel knot formation during tooth morphogenesis in mice.