A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation

PITX2 expression and Neanderthal introgression in HS3ST3A1 contribute to variation in tooth dimensions in modern humans

Dental morphology varies greatly throughout evolution, including in the human lineage, but little is known about the biology of this variation. Here, we use multiomics analyses to examine the genetics of variation in tooth crown dimensions. In a human cohort with mixed continental ancestry, we detected genome-wide significant associations at 18 genome regions. One region includes EDAR, a gene known to impact dental features in East Asians. Furthermore, we find that EDAR variants increase the mesiodistal diameter of all teeth, following an anterior-posterior gradient of decreasing strength. Among the 17 novel-associated regions, we replicate 7/13 in an independent human cohort and find that 4/12 orthologous regions affect molar size in mice. Two association signals point to compelling candidate genes. One is ∼61 kb from PITX2, a major determinant of tooth development. Another overlaps HS3ST3A1, a paralogous neighbor of HS3ST3B1, a tooth enamel knot factor. We document the expression of Pitx2 and Hs3st3a1 in enamel knot and dental epithelial cells of developing mouse incisors. Furthermore, associated SNPs in PITX2 and HS3ST3A1 overlap enhancers active in these cells, suggesting a role for these SNPs in gene regulation during dental development. In addition, we document that Pitx2 and Hs3st3a1/Hs3st3b1 knockout mice show alterations in dental morphology. Finally, we find that associated SNPs in HS3ST3A1 are in a DNA tract introgressed from Neanderthals, consistent with an involvement of HS3ST3A1 in tooth size variation during human evolution.

Mapping genes for human face shape: Exploration of univariate phenotyping strategies

Human facial shape, while strongly heritable, involves both genetic and structural complexity, necessitating precise phenotyping for accurate assessment. Common phenotyping strategies include simplifying 3D facial features into univariate traits such as anthropometric measurements (e.g., inter-landmark distances), unsupervised dimensionality reductions (e.g., principal component analysis (PCA) and auto-encoder (AE) approaches), and assessing resemblance to particular facial gestalts (e.g., syndromic facial archetypes). This study provides a comparative assessment of these strategies in genome-wide association studies (GWASs) of 3D facial shape. Specifically, we investigated inter-landmark distances, PCA and AE-derived latent dimensions, and facial resemblance to random, extreme, and syndromic gestalts within a GWAS of 8,426 individuals of recent European ancestry. Inter-landmark distances exhibit the highest SNP-based heritability as estimated via LD score regression, followed by AE dimensions. Conversely, resemblance scores to extreme and syndromic facial gestalts display the lowest heritability, in line with expectations. Notably, the aggregation of multiple GWASs on facial resemblance to random gestalts reveals the highest number of independent genetic loci. This novel, easy-to-implement phenotyping approach holds significant promise for capturing genetically relevant morphological traits derived from complex biomedical imaging datasets, and its applications extend beyond faces. Nevertheless, these different phenotyping strategies capture different genetic influences on craniofacial shape. Thus, it remains valuable to explore these strategies individually and in combination to gain a more comprehensive understanding of the genetic factors underlying craniofacial shape and related traits.

Archaic hominin admixture and its consequences for modern humans

As anatomically modern humans dispersed out of Africa, they encountered and mated with now extinct hominins, including Neanderthals and Denisovans. It is now well established that all non-African individuals derive approximately 2% of their genome from Neanderthal ancestors and individuals of Melanesian and Australian aboriginal ancestry inherited an additional 2%-5% of their genomes from Denisovan ancestors. Attention has started to shift from documenting amounts of archaic admixture and identifying introgressed segments to understanding their molecular, phenotypic, and evolutionary consequences and refining models of human history. Here, we review recent insights into admixture between modern and archaic humans, emphasizing methodological innovations and the functional and phenotypic effects Neanderthal and Denisovan sequences have in contemporary individuals.

Genetics of Latin American Diversity Project: Insights into population genetics and association studies in admixed groups in the Americas

Latin Americans are underrepresented in genetic studies, increasing disparities in personalized genomic medicine. Despite available genetic data from thousands of Latin Americans, accessing and navigating the bureaucratic hurdles for consent or access remains challenging. To address this, we introduce the Genetics of Latin American Diversity (GLAD) Project, compiling genome-wide information from 53,738 Latin Americans across 39 studies representing 46 geographical regions. Through GLAD, we identified heterogeneous ancestry composition and recent gene flow across the Americas. Additionally, we developed GLAD-match, a simulated annealing-based algorithm, to match the genetic background of external samples to our database, sharing summary statistics (i.e., allele and haplotype frequencies) without transferring individual-level genotypes. Finally, we demonstrate the potential of GLAD as a critical resource for evaluating statistical genetic software in the presence of admixture. By providing this resource, we promote genomic research in Latin Americans and contribute to the promises of personalized medicine to more people.

High heritability of human facial traits reveals associations with CNTLN, BRCA1, and TMPRSS6 loci in Korean families

Facial features are determined by interactions between genetic and environmental factors. However, genes underlying facial similarities in individuals from the same family remain less explored. To identify genetic variants associated with heritable facial features, we investigated familial (parent-offspring) associations and estimated familial correlation and heritability using 39 facial measurements in 408 individuals from 117 Korean families. Facial trait heritability ranged from 0.124 to 0.669. Longitudinal facial growth-related traits were highly heritable, including distances from the nasion to right alare (h 2  = 0.668898), pogonion to midendocanthion (h 2  = 0.661557), subnasale to midendocanthion (h 2  = 0.656882), and morphological facial height (h 2  = 0. 654376). We identified the top three significant genome-wide associated variants in the eye, nose, and lip-jaw regions. CNTLN (rs10511632: beta = -0.02696, p = 1.146 × 10-9) and BRCA1 (rs397509305: beta = 0.02741, p = 7.17 × 10-9) loci were associated with distance from the nasion to the right alare. The TMPRSS6 (rs228913: beta = 0.05101, p = 3.68 × 10-9) locus was associated with the distance from the labiale superius to the pogonion and lower facial height. These associations were maintained in an independent unrelated population. In conclusion, we identified new gene variants associated with longitudinal facial morphology that may affect individual facial differences, which has important implications for clinical and forensic applications.

Genetic associations with disease in populations with Indigenous American ancestries

The genetic architecture of complex diseases affecting populations with Indigenous American ancestries is poorly understood due to their underrepresentation in genomics studies. While most of the genetic diversity associated with disease trait variation is shared among worldwide populations, a fraction of this component is expected to be unique to each continental group, including Indigenous Americans. Here, I describe the current state of knowledge from genome-wide association studies on Indigenous populations, as well as non-Indigenous populations with partial Indigenous ancestries from the American continent, focusing on disease susceptibility and anthropometric traits. While some studies identified risk alleles unique to Indigenous populations, their effects on trait variation are mostly small. I suggest that the associations rendered by many inter-population studies are probably inflated due to the absence of socio-cultural-economic covariates in the association models. I encourage the inclusion of admixed individuals in future GWAS studies to control for inter-ancestry differences in environmental factors. I suggest that some complex diseases might have arisen as trade-off costs of adaptations to past evolutionary selective pressures. Finally, I discuss how expanding panels with Indigenous ancestries in GWAS studies is key to accurately assess genetic risk in populations from the American continent, thus decreasing global health disparities.

Association between thermogenic brown fat and genes under positive natural selection in circumpolar populations

BACKGROUND

Adaptation to cold was essential for human migration across Eurasia. Non-shivering thermogenesis through brown adipose tissue (BAT) participates in cold adaptation because some genes involved in the differentiation and function of BAT exhibit signatures of positive natural selection in populations at high latitudes. Whether these genes are associated with the inter-individual variability in BAT thermogenesis remains unclear. In this study, we evaluated the potential associations between BAT activity and single nucleotide polymorphisms (SNPs) in candidate gene regions in East Asian populations.

METHODS

BAT activity induced by mild cold exposure was measured in 399 healthy Japanese men and women using fluorodeoxyglucose-positron emission tomography and computed tomography (FDG-PET/CT). The capacity for cold-induced thermogenesis and fat oxidation was measured in 56 men. Association analyses with physiological traits were performed for 11 SNPs at six loci (LEPR, ANGPTL8, PLA2G2A, PLIN1, TBX15-WARS2, and FADS1) reported to be under positive natural selection. Associations found in the FDG-PET/CT population were further validated in 84 healthy East Asian men and women, in whom BAT activity was measured using infrared thermography. Associations between the SNP genotypes and BAT activity or other related traits were tested using multiple logistic and linear regression models.

RESULTS

Of the 11 putative adaptive alleles of the six genes, two intronic SNPs in LEPR (rs1022981 and rs12405556) tended to be associated with higher BAT activity. However, these did not survive multiple test comparisons. Associations with lower body fat percentage, plasma triglyceride, insulin, and HOMA-IR levels were observed in the FDG-PET/CT population (P < 0.05). Other loci, including TBX15-WARS2, which is speculated to mediate cold adaptation in Greenland Inuits, did not show significant differences in BAT thermogenesis.

CONCLUSIONS

Our results suggest a marginal but significant association between LEPR SNPs. However, robust supporting evidence was not established for the involvement of other loci under positive natural selection in cold adaptation through BAT thermogenesis in East Asian adults. Given the pleiotropic function of these genes, factors other than cold adaptation through BAT thermogenesis, such as diet adaptation, may contribute to positive natural selection at these loci.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Selection scan in Native Americans of Mexico identifies FADS2 rs174616: Evidence of gene-diet interactions affecting lipid levels and Delta-6-desaturase activity

Searching for positive selection signals across genomes has identified functional genetic variants responding to environmental change. In Native Americans of Mexico, we used the fixation index (Fst) and population branch statistic (PBS) to identify SNPs suggesting positive selection. The 103 most differentiated SNPs were tested for associations with metabolic traits, the most significant association was FADS2/rs174616 with body mass index (BMI). This variant lies within a linkage disequilibrium (LD) block independent of previously reported FADS selection signals and has not been clearly associated with metabolic phenotypes. We tested this variant in two independent cohorts with cardiometabolic data. In the Genetics of Atherosclerotic Disease (GEA) cohort, the derived allele (T) was associated with increased BMI, lower LDL-C levels and a decreased risk of subclinical atherosclerosis in women. Significant gene-diet interactions affected lipid, apolipoprotein and adiponectin levels with differences according to sex, involving mainly total and complex dietary carbohydrate%. In the Genotype-related Effects of PUFA trial, the derived allele was associated with lower Δ-6 desaturase activity and erythrocyte membrane dihomo-gamma-linolenic acid (DGLA) levels, and with increased Δ-5 desaturase activity and eicosapentaenoic acid levels. This variant interacted with dietary carbohydrate% affecting Δ-6 desaturase activity. Notably, the relationship of DGLA and other erythrocyte membrane LC-PUFA indices with HOMA-IR differed according to rs174616 genotype, which has implications regarding how these indices should be interpreted. In conclusion, this observational study identified rs174616 as a signal suggesting selection in an independent linkage disequilibrium block, was associated with cardiometabolic and erythrocyte measurements of LC-PUFA in two independent Mexican cohorts and showed significant gene-diet interactions.

ZSCAN10 deficiency causes a neurodevelopmental disorder with characteristic oto-facial malformations

Neurodevelopmental disorders are major indications for genetic referral and have been linked to more than 1500 loci including genes encoding transcriptional regulators. The dysfunction of transcription factors often results in characteristic syndromic presentations; however, at least half of these patients lack a genetic diagnosis. The implementation of machine learning approaches has the potential to aid in the identification of new disease genes and delineate associated phenotypes. Next generation sequencing was performed in seven affected individuals with neurodevelopmental delay and dysmorphic features. Clinical characterization included reanalysis of available neuroimaging datasets and 2D portrait image analysis with GestaltMatcher. The functional consequences of ZSCAN10 loss were modelled in mouse embryonic stem cells (mESCs), including a knockout and a representative ZSCAN10 protein truncating variant. These models were characterized by gene expression and western blot analyses, chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR) and immunofluorescence staining. Zscan10 knockout mouse embryos were generated and phenotyped. We prioritized bi-allelic ZSCAN10 loss-of-function variants in seven affected individuals from five unrelated families as the underlying molecular cause. RNA-sequencing analyses in Zscan10-/- mESCs indicated dysregulation of genes related to stem cell pluripotency. In addition, we established in mESCs the loss-of-function mechanism for a representative human ZSCAN10 protein truncating variant by showing alteration of its expression levels and subcellular localization, interfering with its binding to DNA enhancer targets. Deep phenotyping revealed global developmental delay, facial asymmetry and malformations of the outer ear as consistent clinical features. Cerebral MRI showed dysplasia of the semicircular canals as an anatomical correlate of sensorineural hearing loss. Facial asymmetry was confirmed as a clinical feature by GestaltMatcher and was recapitulated in the Zscan10 mouse model along with inner and outer ear malformations. Our findings provide evidence of a novel syndromic neurodevelopmental disorder caused by bi-allelic loss-of-function variants in ZSCAN10.

Exploring the pathological mechanisms underlying Cohen syndrome

Cohen Syndrome (CS) is a rare autosomal recessive disorder caused by biallelic mutations in the VPS13B gene. It is characterized by multiple clinical features, including acquired microcephaly, developmental delay, intellectual disability, neutropenia, and retinal degeneration. VPS13B is part of the bridge-like lipid transport (BLTP) protein family, which in mammals also includes VPS13A, -C, and -D. The proteins of this family are peripheral membrane proteins with different sub-cellular localization, but all share similar structural features and have been proposed to act as lipid transport proteins at organellar membrane contact sites. VPS13B is localized at the Golgi apparatus and is essential for the maintenance of organelle architecture. Here we present a review of the experimental data on the function of the protein at the cellular level, discussing the potential link with disease phenotype and review the studies on animal models recapitulating features of the human disease.

Mapping genes for human face shape: exploration of univariate phenotyping strategies

Human facial shape, while strongly heritable, involves both genetic and structural complexity, necessitating precise phenotyping for accurate assessment. Common phenotyping strategies include simplifying 3D facial features into univariate traits such as anthropometric measurements (e.g., inter-landmark distances), unsupervised dimensionality reductions (e.g., principal component analysis (PCA) and auto-encoder (AE) approaches), and assessing resemblance to particular facial gestalts (e.g., syndromic facial archetypes). This study provides a comparative assessment of these strategies in genome-wide association studies (GWASs) of 3D facial shape. Specifically, we investigated inter-landmark distances, PCA and AE-derived latent dimensions, and facial resemblance to random, extreme, and syndromic gestalts within a GWAS of 8,426 individuals of recent European ancestry. Inter-landmark distances exhibit the highest SNP-based heritability as estimated via LD score regression, followed by AE dimensions. Conversely, resemblance scores to extreme and syndromic facial gestalts display the lowest heritability, in line with expectations. Notably, the aggregation of multiple GWASs on facial resemblance to random gestalts reveals the highest number of independent genetic loci. This novel, easy-to-implement phenotyping approach holds significant promise for capturing genetically relevant morphological traits derived from complex biomedical imaging datasets, and its applications extend beyond faces. Nevertheless, these different phenotyping strategies capture different genetic influences on craniofacial shape. Thus, it remains valuable to explore these strategies individually and in combination to gain a more comprehensive understanding of the genetic factors underlying craniofacial shape and related traits.

Author Summary

Advancements linking variation in the human genome to phenotypes have rapidly evolved in recent decades and have revealed that most human traits are influenced by genetic variants to at least some degree. While many traits, such as stature, are straightforward to acquire and investigate, the multivariate and multipartite nature of facial shape makes quantification more challenging. In this study, we compared the impact of different facial phenotyping approaches on gene mapping outcomes. Our findings suggest that the choice of facial phenotyping method has an impact on apparent trait heritability and the ability to detect genetic association signals. These results offer valuable insights into the importance of phenotyping in genetic investigations, especially when dealing with highly complex morphological traits.

Facial modeling and measurement based upon homologous topographical features

Measurement of human faces is fundamental to many applications from recognition to genetic phenotyping. While anthropometric landmarks provide a conventional set of homologous measurement points, digital scans are increasingly used for facial measurement, despite the difficulties in establishing their homology. We introduce an alternative basis for facial measurement, which 1) provides a richer information density than discrete point measurements, 2) derives its homology from shared facial topography (ridges, folds, etc.), and 3) quantifies local morphological variation following the conventions and practices of anatomical description. A parametric model that permits matching a broad range of facial variation by the adjustment of 71 parameters is demonstrated by modeling a sample of 80 adult human faces. The surface of the parametric model can be adjusted to match each photogrammetric surface mesh generally to within 1 mm, demonstrating a novel and efficient means for facial shape encoding. We examine how well this scheme quantifies facial shape and variation with respect to geographic ancestry and sex. We compare this analysis with a more conventional, landmark-based geometric morphometric (GMM) study with 43 landmarks placed on the same set of scans. Our multivariate statistical analysis using the 71 attribute values separates geographic ancestry groups and sexes with a high degree of reliability, and these results are broadly similar to those from GMM, but with some key differences that we discuss. This approach is compared with conventional, non-parametric methods for the quantification of facial shape, including generality, information density, and the separation of size and shape. Potential uses for phenotypic and dysmorphology studies are also discussed.

The genetic changes that shaped Neandertals, Denisovans, and modern humans

Modern human ancestors diverged from the ancestors of Neandertals and Denisovans about 600,000 years ago. Until about 40,000 years ago, these three groups existed in parallel, occasionally met, and exchanged genes. A critical question is why modern humans, and not the other two groups, survived, became numerous, and developed complex cultures. Here, we discuss genetic differences among the groups and some of their functional consequences. As more present-day genome sequences become available from diverse groups, we predict that very few, if any, differences will distinguish all modern humans from all Neandertals and Denisovans. We propose that the genetic basis of what constitutes a modern human is best thought of as a combination of genetic features, where perhaps none of them is present in each and every present-day individual.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

A genome-wide association scan reveals novel loci for facial traits of Koreans

DNA-based prediction of externally visible characteristics (EVC) with SNPs is one of the research areas of interest in the forensic field. Based on a previous study performing GWAS on facial traits in a Korean population, herein, we present results stemming from GWA analysis with KoreanChip and novel genetic loci satisfying genome-wide significant level. We discovered a total of 20 signals and 12 loci were found to have novel associations with facial traits, including six loci located in intergenic regions and six loci located at UBE2O, HECTD2, CCDC108, TPK1, FCN2, and FRMPD1. Additionally, we performed a polygenic score analysis for 33 distance-related traits in facial phenotyping and determined genetic relationships between facial traits and SNPs using the GCTA program. The results of the current study offer an understanding of how facial morphology is influenced by complex genetic structures and provide insights into forensic investigation and population genetics.

Neanderthal introgression in SCN9A impacts mechanical pain sensitivity

The Nav1.7 voltage-gated sodium channel plays a key role in nociception. Three functional variants in the SCN9A gene (encoding M932L, V991L, and D1908G in Nav1.7), have recently been identified as stemming from Neanderthal introgression and to associate with pain symptomatology in UK BioBank data. In 1000 genomes data, these variants are absent in Europeans but common in Latin Americans. Analysing high-density genotype data from 7594 Latin Americans, we characterized Neanderthal introgression in SCN9A. We find that tracts of introgression occur on a Native American genomic background, have an average length of ~123 kb and overlap the M932L, V991L, and D1908G coding positions. Furthermore, we measured experimentally six pain thresholds in 1623 healthy Colombians. We found that Neanderthal ancestry in SCN9A is significantly associated with a lower mechanical pain threshold after sensitization with mustard oil and evidence of additivity of effects across Nav1.7 variants. Our findings support the reported association of Neanderthal Nav1.7 variants with clinical pain, define a specific sensory modality affected by archaic introgression in SCN9A and are consistent with independent effects of the Neanderthal variants on Nav1.7 function.

Characterization of Vps13b-mutant mice reveals neuroanatomical and behavioral phenotypes with females less affected

The vacuolar protein sorting-associated protein 13B (VPS13B) is a large and highly conserved protein. Disruption of VPS13B causes the autosomal recessive Cohen syndrome, a rare disorder characterized by microcephaly and intellectual disability among other features, including developmental delay, hypotonia, and friendly-personality. However, the underlying mechanisms by which VPS13B disruption leads to brain dysfunction still remain unexplained. To gain insights into the neuropathogenesis of Cohen syndrome, we systematically characterized brain changes in Vps13b-mutant mice and compared murine findings to 235 previously published and 17 new patients diagnosed with VPS13B-related Cohen syndrome. We showed that Vps13b is differentially expressed across brain regions with the highest expression in the cerebellum, the hippocampus and the cortex with postnatal peak. Half of the Vps13b-/- mice die during the first week of life. The remaining mice have a normal lifespan and display the core phenotypes of the human disease, including microcephaly, growth delay, hypotonia, altered memory, and enhanced sociability. Systematic 2D and 3D brain histo-morphological analyses reveal specific structural changes in the brain starting after birth. The dentate gyrus is the brain region with the most prominent reduction in size, while the motor cortex is specifically thinner in layer VI. The fornix, the fasciculus retroflexus, and the cingulate cortex remain unaffected. Interestingly, these neuroanatomical changes implicate an increase of neuronal death during infantile stages with no progression in adulthood suggesting that VPS13B promotes neuronal survival early in life. Importantly, whilst both sexes were affected, some neuroanatomical and behavioral phenotypes were less pronounced or even absent in females. We evaluate sex differences in Cohen patients and conclude that females are less affected both in mice and patients. Our findings provide new insights about the neurobiology of VPS13B and highlight previously unreported brain phenotypes while defining Cohen syndrome as a likely new entity of non-progressive infantile neurodegeneration.

Data-driven trait heritability-based extraction of human facial phenotypes

A genome-wide association study (GWAS) of a complex, multi-dimensional morphological trait, such as the human face, typically relies on predefined and simplified phenotypic measurements, such as inter-landmark distances and angles. These measures are predominantly designed by human experts based on perceived biological or clinical knowledge. To avoid use handcrafted phenotypes (i.e., a priori expert-identified phenotypes), alternative automatically extracted phenotypic descriptors, such as features derived from dimension reduction techniques (e.g., principal component analysis), are employed. While the features generated by such computational algorithms capture the geometric variations of the biological shape, they are not necessarily genetically relevant. Therefore, genetically informed data-driven phenotyping is desirable. Here, we propose an approach where phenotyping is done through a data-driven optimization of trait heritability, defined as the degree of variation in a phenotypic trait in a population that is due to genetic variation. The resulting phenotyping process consists of two steps: 1) constructing a feature space that models shape variations using dimension reduction techniques, and 2) searching for directions in the feature space exhibiting high trait heritability using a genetic search algorithm (i.e., heuristic inspired by natural selection). We show that the phenotypes resulting from the proposed trait heritability-optimized training differ from those of principal components in the following aspects: 1) higher trait heritability, 2) higher SNP heritability, and 3) identification of the same number of independent genetic loci with a smaller number of effective traits. Our results demonstrate that data-driven trait heritability-based optimization enables the automatic extraction of genetically relevant phenotypes, as shown by their increased power in genome-wide association scans.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Jaw size variation is associated with a novel craniofacial function for galanin receptor 2 in an adaptive radiation of pupfishes

Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology that also provides an opportunity to uncover new genes and regulatory networks with potential clinical relevance. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island in the Bahamas. We confirmed the loss of a putative Sry transcription factor binding site in the upstream region of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage and premaxilla using in situ hybridization chain reaction (HCR). We then experimentally demonstrated a novel function for Galr2 in craniofacial development and jaw elongation by exposing embryos to drugs that inhibit Galr2 activity. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a postulated Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

Recent advances in Forensic DNA Phenotyping of appearance, ancestry and age

Forensic DNA Phenotyping (FDP) comprises the prediction of a person's externally visible characteristics regarding appearance, biogeographic ancestry and age from DNA of crime scene samples, to provide investigative leads to help find unknown perpetrators that cannot be identified with forensic STR-profiling. In recent years, FDP has advanced considerably in all of its three components, which we summarize in this review article. Appearance prediction from DNA has broadened beyond eye, hair and skin color to additionally comprise other traits such as eyebrow color, freckles, hair structure, hair loss in men, and tall stature. Biogeographic ancestry inference from DNA has progressed from continental ancestry to sub-continental ancestry detection and the resolving of co-ancestry patterns in genetically admixed individuals. Age estimation from DNA has widened beyond blood to more somatic tissues such as saliva and bones as well as new markers and tools for semen. Technological progress has allowed forensically suitable DNA technology with largely increased multiplex capacity for the simultaneous analysis of hundreds of DNA predictors with targeted massively parallel sequencing (MPS). Forensically validated MPS-based FDP tools for predicting from crime scene DNA i) several appearance traits, ii) multi-regional ancestry, iii) several appearance traits together with multi-regional ancestry, and iv) age from different tissue types, are already available. Despite recent advances that will likely increase the impact of FDP in criminal casework in the near future, moving reliable appearance, ancestry and age prediction from crime scene DNA to the level of detail and accuracy police investigators may desire, requires further intensified scientific research together with technical developments and forensic validations as well as the necessary funding.

Long divergent haplotypes introgressed from wild sheep are associated with distinct morphological and adaptive characteristics in domestic sheep

The worldwide sheep population comprises more than 1000 breeds. Together, these exhibit a considerable morphological diversity, which has not been extensively investigated at the molecular level. Here, we analyze whole-genome sequencing individuals of 1,098 domestic sheep from 154 breeds, and 69 wild sheep from seven Ovis species. On average, we detected 6.8%, 1.0% and 0.2% introgressed sequence in domestic sheep originating from Iranian mouflon, urial and argali, respectively, with rare introgressions from other wild species. Interestingly, several introgressed haplotypes contributed to the morphological differentiations across sheep breeds, such as a RXFP2 haplotype from Iranian mouflon conferring the spiral horn trait, a MSRB3 haplotype from argali strongly associated with ear morphology, and a VPS13B haplotype probably originating from urial and mouflon possibly associated with facial traits. Our results reveal that introgression events from wild Ovis species contributed to the high rate of morphological differentiation in sheep breeds, but also to individual variation within breeds. We propose that long divergent haplotypes are a ubiquitous source of phenotypic variation that allows adaptation to a variable environment, and that these remain intact in the receiving population probably due to reduced recombination.

Aesthetic considerations when treating the Latin American patient: Thriving in diversity international roundtable series

BACKGROUND

The Hispanic/Latin American population is the fastest growing non-Caucasian group in the United States. Within this group, demand for aesthetic procedures is on the rise. High ethnic variability among these patients influences treatment approaches and patient priorities. Understanding these ethnic differences is central to providing optimal care.

AIMS

To discuss similarities and differences in anatomy and treatment preferences of Hispanic/Latin American patients both within the United States and internationally and explore how these differences may influence or inform aesthetic practices.

PATIENTS/METHODS

In support of clinicians who wish to serve a diverse patient population, a 6-part, international roundtable series focused on diversity in aesthetics was conducted from August 24, 2021 to May 16, 2022. In this roundtable, held in Medellin, Columbia, expert clinicians from across Latin America and the United States were invited to contribute and share best practices.

RESULTS

The results of the second roundtable in the series, the Latin American Patient, are described here. A special emphasis is placed on procedures that address the most commonly encountered concerns in these patients.

CONCLUSIONS

Hispanic and Latino patients represent a broad demographic with unique anatomical features, aesthetic preferences, and treatment priorities. Clinicians should consider these differences when treating this patient population.

Advancement in Human Face Prediction Using DNA

The rapid improvements in identifying the genetic factors contributing to facial morphology have enabled the early identification of craniofacial syndromes. Similarly, this technology can be vital in forensic cases involving human identification from biological traces or human remains, especially when reference samples are not available in the deoxyribose nucleic acid (DNA) database. This review summarizes the currently used methods for predicting human phenotypes such as age, ancestry, pigmentation, and facial features based on genetic variations. To identify the facial features affected by DNA, various two-dimensional (2D)- and three-dimensional (3D)-scanning techniques and analysis tools are reviewed. A comparison between the scanning technologies is also presented in this review. Face-landmarking techniques and face-phenotyping algorithms are discussed in chronological order. Then, the latest approaches in genetic to 3D face shape analysis are emphasized. A systematic review of the current markers that passed the threshold of a genome-wide association (GWAS) of single nucleotide polymorphism (SNP)-face traits from the GWAS Catalog is also provided using the preferred reporting items for systematic reviews and meta-analyses (PRISMA), approach. Finally, the current challenges in forensic DNA phenotyping are analyzed and discussed.

Combining genome-wide association studies highlight novel loci involved in human facial variation

Standard genome-wide association studies (GWASs) rely on analyzing a single trait at a time. However, many human phenotypes are complex and composed by multiple correlated traits. Here we introduce C-GWAS, a method for combining GWAS summary statistics of multiple potentially correlated traits. Extensive computer simulations demonstrated increased statistical power of C-GWAS compared to the minimal p-values of multiple single-trait GWASs (MinGWAS) and the current state-of-the-art method for combining single-trait GWASs (MTAG). Applying C-GWAS to a meta-analysis dataset of 78 single trait facial GWASs from 10,115 Europeans identified 56 study-wide suggestively significant loci with multi-trait effects on facial morphology of which 17 are novel loci. Using data from additional 13,622 European and Asian samples, 46 (82%) loci, including 9 (53%) novel loci, were replicated at nominal significance with consistent allele effects. Functional analyses further strengthen the reliability of our C-GWAS findings. Our study introduces the C-GWAS method and makes it available as computationally efficient open-source R package for widespread future use. Our work also provides insights into the genetic architecture of human facial appearance.

Integration of multimodal data in the developing tooth reveals candidate regulatory loci driving human odontogenic phenotypes

Human odontogenic aberrations such as abnormal tooth number and delayed tooth eruption can occur as a symptom of rare syndromes or, more commonly, as nonsyndromic phenotypes. These phenotypes can require extensive and expensive dental treatment, posing a significant burden. While many dental phenotypes are heritable, most nonsyndromic cases have not been linked to causal genes. We demonstrate the novel finding that common sequence variants associated with human odontogenic phenotypes are enriched in developmental craniofacial enhancers conserved between human and mouse. However, the bulk nature of these samples obscures if this finding is due to the tooth itself or the surrounding tissues. We therefore sought to identify enhancers specifically active in the tooth anlagen and quantify their contribution to the observed genetic enrichments. We systematically identified 22,001 conserved enhancers active in E13.5 mouse incisors using ChIP-seq and machine learning pipelines and demonstrated biologically relevant enrichments in putative target genes, transcription factor binding motifs, and in vivo activity. Multi-tissue comparisons of human and mouse enhancers revealed that these putative tooth enhancers had the strongest enrichment of odontogenic phenotype-associated variants, suggesting a role for dysregulation of tooth developmental enhancers in human dental phenotypes. The large number of these regions genome-wide necessitated prioritization of enhancer loci for future investigations. As enhancers modulate gene expression, we prioritized regions based on enhancers' putative target genes. We predicted these target genes and prioritized loci by integrating chromatin state, bulk gene expression and coexpression, GWAS variants, and cell type resolved gene expression to generate a prioritized list of putative odontogenic phenotype-driving loci active in the developing tooth. These genomic regions are of particular interest for downstream experiments determining the role of specific dental enhancer:gene pairs in odontogenesis.

Decoding the Human Face: Progress and Challenges in Understanding the Genetics of Craniofacial Morphology

Variations in the form of the human face, which plays a role in our individual identities and societal interactions, have fascinated scientists and artists alike. Here, we review our current understanding of the genetics underlying variation in craniofacial morphology and disease-associated dysmorphology, synthesizing decades of progress on Mendelian syndromes in addition to more recent results from genome-wide association studies of human facial shape and disease risk. We also discuss the various approaches used to phenotype and quantify facial shape, which are of particular importance due to the complex, multipartite nature of the craniofacial form. We close by discussing how experimental studies have contributed and will further contribute to our understanding of human genetic variation and then proposing future directions and applications for the field.

Classification of early age facial growth pattern and identification of the genetic basis in two Korean populations

Childhood to adolescence is an accelerated growth period, and genetic features can influence differences of individual growth patterns. In this study, we examined the genetic basis of early age facial growth (EAFG) patterns. Facial shape phenotypes were defined using facial landmark distances, identifying five growth patterns: continued-decrease, decrease-to-increase, constant, increase-to-decrease, and continued-increase. We conducted genome-wide association studies (GWAS) for 10 horizontal and 11 vertical phenotypes. The most significant association for horizontal phenotypes was rs610831 (TRIM29; β = 0.92, p-value = 1.9 × 10⁻⁹) and for vertical phenotypes was rs6898746 (ZSWIM6; β = 0.1103, p-value = 2.5 × 10⁻⁸). It is highly correlated with genes already reported for facial growth. This study is the first to classify and characterize facial growth patterns and related genetic polymorphisms.

Genetic variants underlying differences in facial morphology in East Asian and European populations

Facial morphology-a conspicuous feature of human appearance-is highly heritable. Previous studies on the genetic basis of facial morphology were performed mainly in European-ancestry cohorts (EUR). Applying a data-driven phenotyping and multivariate genome-wide scanning protocol to a large collection of three-dimensional facial images of individuals with East Asian ancestry (EAS), we identified 244 variants in 166 loci (62 new) associated with typical-range facial variation. A newly proposed polygenic shape analysis indicates that the effects of the variants on facial shape in EAS can be generalized to EUR. Based on this, we further identified 13 variants related to differences between facial shape in EUR and EAS populations. Evolutionary analyses suggest that the difference in nose shape between EUR and EAS populations is caused by a directional selection, due mainly to a local adaptation in Europeans. Our results illustrate the underlying genetic basis for facial differences across populations.

The effects of Tbx15 and Pax1 on facial and other physical morphology in mice

DNA variants in or close to the human TBX15 and PAX1 genes have been repeatedly associated with facial morphology in independent genome-wide association studies, while their functional roles in determining facial morphology remain to be understood. We generated Tbx15 knockout (Tbx15 -/-) and Pax1 knockout (Pax1 -/-) mice by applying the one-step CRISPR/Cas9 method. A total of 75 adult mice were used for subsequent phenotype analysis, including 38 Tbx15 mice (10 homozygous Tbx15 -/-, 18 heterozygous Tbx15 +/-, 10 wild-type Tbx15 +/+ WT littermates) and 37 Pax1 mice (12 homozygous Pax1 -/-, 15 heterozygous Pax1 +/-, 10 Pax1 +/+ WT littermates). Facial and other physical morphological phenotypes were obtained from three-dimensional (3D) images acquired with the HandySCAN BLACK scanner. Compared to WT littermates, the Tbx15 -/- mutant mice had significantly shorter faces (p = 1.08E-8, R2 = 0.61) and their ears were in a significantly lower position (p = 3.54E-8, R2 = 0.62) manifesting a "droopy ear" characteristic. Besides these face alternations, Tbx15 -/- mutant mice displayed significantly lower weight as well as shorter body and limb length. Pax1 -/- mutant mice showed significantly longer noses (p = 1.14E-5, R2 = 0.46) relative to WT littermates, but otherwise displayed less obvious morphological alterations than Tbx15 -/- mutant mice did. We provide the first direct functional evidence that two well-known and replicated human face genes, Tbx15 and Pax1, impact facial and other body morphology in mice. The general agreement between our findings in knock-out mice with those from previous GWASs suggests that the functional evidence we established here in mice may also be relevant in humans.

Deep ancestral introgression shapes evolutionary history of dragonflies and damselflies

Introgression is an important biological process affecting at least 10% of the extant species in the animal kingdom. Introgression significantly impacts inference of phylogenetic species relationships where a strictly binary tree model cannot adequately explain reticulate net-like species relationships. Here we use phylogenomic approaches to understand patterns of introgression along the evolutionary history of a unique, non-model insect system: dragonflies and damselflies (Odonata). We demonstrate that introgression is a pervasive evolutionary force across various taxonomic levels within Odonata. In particular, we show that the morphologically "intermediate" species of Anisozygoptera (one of the three primary suborders within Odonata besides Zygoptera and Anisoptera), which retain phenotypic characteristics of the other two suborders, experienced high levels of introgression likely coming from zygopteran genomes. Additionally, we find evidence for multiple cases of deep inter-superfamilial ancestral introgression.

Haplotype-based membership inference from summary genomic data

Motivation

The availability of human genomic data, together with the enhanced capacity to process them, is leading to transformative technological advances in biomedical science and engineering. However, the public dissemination of such data has been difficult due to privacy concerns. Specifically, it has been shown that the presence of a human subject in a case group can be inferred from the shared summary statistics of the group, e.g. the allele frequencies, or even the presence/absence of genetic variants (e.g. shared by the Beacon project) in the group. These methods rely on the availability of the target’s genome, i.e. the DNA profile of a target human subject, and thus are often referred to as the membership inference method.

Results

In this article, we demonstrate the haplotypes, i.e. the sequence of single nucleotide variations (SNVs) showing strong genetic linkages in human genome databases, may be inferred from the summary of genomic data without using a target’s genome. Furthermore, novel haplotypes that did not appear in the database may be reconstructed solely from the allele frequencies from genomic datasets. These reconstructed haplotypes can be used for a haplotype-based membership inference algorithm to identify target subjects in a case group with greater power than existing methods based on SNVs.

Availability and implementation

The implementation of the membership inference algorithms is available at https://github.com/diybu/Haplotype-based-membership-inferences.

A genome-wide association study identifies novel gene associations with facial skin wrinkling and mole count in Latin Americans

BACKGROUND

Genome-wide association studies (GWASs) have identified genes influencing skin ageing and mole count in Europeans, but little is known about the relevance of these (or other genes) in non-Europeans.

OBJECTIVES

To conduct a GWAS for facial skin ageing and mole count in adults < 40 years old, of mixed European, Native American and African ancestry, recruited in Latin America.

METHODS

Skin ageing and mole count scores were obtained from facial photographs of over 6000 individuals. After quality control checks, three wrinkling traits and mole count were retained for genetic analyses. DNA samples were genotyped with Illumina's HumanOmniExpress chip. Association testing was performed on around 8 703 729 single-nucleotide polymorphisms (SNPs) across the autosomal genome.

RESULTS

Genome-wide significant association was observed at four genome regions: two were associated with wrinkling (in 1p13·3 and 21q21·2), one with mole count (in 1q32·3) and one with both wrinkling and mole count (in 5p13·2). Associated SNPs in 5p13·2 and in 1p13·3 are intronic within SLC45A2 and VAV3, respectively, while SNPs in 1q32·3 are near the SLC30A1 gene, and those in 21q21·2 occur in a gene desert. Analyses of SNPs in IRF4 and MC1R are consistent with a role of these genes in skin ageing.

CONCLUSIONS

We replicate the association of wrinkling with variants in SLC45A2, IRF4 and MC1R reported in Europeans. We identify VAV3 and SLC30A1 as two novel candidate genes impacting on wrinkling and mole count, respectively. We provide the first evidence that SLC45A2 influences mole count, in addition to variants in this gene affecting melanoma risk in Europeans.