ABSTRACT
Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.
Subject(s)
Phosphoric Monoester Hydrolases , Tungsten Compounds , Catalysis , Mutation , Kinetics , Binding SitesABSTRACT
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.
Subject(s)
Genetic Loci , Genome-Wide Association Study , Humans , Male , Animals , Mice , Genome-Wide Association Study/methods , Phenotype , Asia , Polymorphism, Single Nucleotide/geneticsABSTRACT
BACKGROUND: The profound influence of orthodontic treatments on facial aesthetics has been a topic of increasing interest. This study delves into the intricate interplay between orthodontic treatments, facial feature alterations, and aesthetic perceptions. METHODS: A total of 73 patients who had undergone orthodontic treatment were included in this study. Facial photographs were taken before and after treatment. Ten orthodontists provided facial aesthetic ratings (FAR) for each patient's frontal, profile, and overall views. 48 facial landmarks were manually placed by the orthodontists and normalized using Generalized Procrustes analysis (GPA). Two types of phenotypes were derived from facial landmarks. Global facial phenotypes were then extracted using principal component analysis (PCA). Additionally, 37 clinical features related to aesthetics and orthodontics were extracted. The association between facial features and changes in FAR after orthodontic treatment was determined using these two types of phenotypes. RESULTS: The FAR exhibited a high correlation among orthodontic experts, particularly in the profile view. The FAR increased after orthodontic treatment, especially in profile views. Extraction of premolars and orthognathic surgery were found to result in higher FAR change. For global facial phenotypes, the most noticeable changes in the frontal and profile views associated with FAR occurred in the lip area, characterized by inward retraction of the lips and slight chin protrusion in the profile view, as well as a decrease in lip height in the frontal view. The changes observed in the profile view were statistically more significant than those in the frontal view. These facial changes were consistent with the changes from orthodontic treatment. For clinical features, two profile features, namely pg.sm.hori and pg.n.ls, were found to be associated with FAR following orthodontic treatment. The highest FAR scores were achieved when pg.sm.hori was at 80° and pg.n.ls was at 8°. On the other hand, frontal clinical features had a subtle effect on FAR during orthodontic treatment. CONCLUSIONS: This study demonstrated that orthodontic treatment improves facial aesthetics, particularly at lip aera in the profile view. Profile clinical features, such as pg.sm.hori and pg.n.ls, are essential in orthodontic treatment which could increase facial aesthetics.
Subject(s)
Esthetics, Dental , Face , Humans , Retrospective Studies , Lip , ChinABSTRACT
Two series of novel 4-phenoxypyridine derivatives containing imidazole-4-carboxamide and 4-methyl-5-oxo-4,5-dihydro-1,2,4-triazole-3-carboxamide moieties were synthesized and evaluated for their in vitro inhibitory activities against c-Met kinase and antiproliferative activities against MKN-45, A549 and H460 cancer cell lines. The results indicated that most of the compounds showed moderate to good antitumor activities. The most promising compound T14 (with c-Met IC50 value of 0.012 µM) showed remarkable antiproliferative activities against MKN-45, A549 and H460 cell lines with IC50 values of 0.64 µM, 1.92 µM and 2.68 µM, respectively. Their preliminary structure-activity relationships (SARs) studies indicate that imidazole-4-carboxamide was more preferred as linker part, and electron-withdrawing groups (especially halogen groups) on the terminal phenyl rings were beneficial for improving the antitumor activities.
Subject(s)
Antineoplastic Agents , Quinolines , Aminoimidazole Carboxamide/analogs & derivatives , Cell Line, Tumor , Cell Proliferation , Drug Design , Drug Screening Assays, Antitumor , Protein Kinase Inhibitors , Proto-Oncogene Proteins c-met , Quinolines/pharmacology , Structure-Activity Relationship , TriazolesABSTRACT
Pancreatic cancer is highly lethal with limited effective treatments. This study explores the therapeutic effects of eupalinolide B (EB) from Eupatorium lindleyanum DC on pancreatic cancer cells. Through cellular functional assays, we observed that EB effectively inhibits cell viability, proliferation, migration, and invasion. In a xenograft mouse model, EB treatment resulted in reduced pancreatic cancer tumor growth and decreased expression of Ki-67. Mechanistically, EB induces apoptosis, elevates reactive oxygen species (ROS) levels, and disrupts copper homeostasis. RNA sequencing (RNA-seq) and gene set enrichment analysis (GSEA) identified copper ion binding pathways and potential involvement in cuproptosis. Furthermore, EB enhances the cytotoxic effects of elesclomol (ES), increasing ROS levels in a copper-dependent manner and exhibiting synergistic cytotoxicity. These findings suggest that EB, either alone or in combination with ES, represents a promising strategy for targeting metal ion dysregulation and inducing potential cuproptosis in pancreatic cancer, offering a potential improvement in therapeutic outcomes.
ABSTRACT
Glioblastoma (GBM) remains the most lethal primary brain tumor, characterized by dismal survival rates. Novel molecular targets are urgently required to enhance therapeutic outcomes. A combination of bioinformatics analysis and experimental validation was employed to investigate the role of EGFLAM in GBM. The Chinese Glioma Genome Atlas provided a platform for gene expression profiling, while siRNA-mediated knockdown and overexpression assays in GBM cell lines, alongside in vivo tumorigenesis models, facilitated functional validation. EGFLAM was found to be significantly overexpressed in GBM tissues, correlating with adverse prognostic factors and higher tumor grades, particularly in patients over the age of 41. Functional assays indicated that EGFLAM is vital for maintaining GBM cell proliferation, viability, and invasiveness. Knockdown of EGFLAM expression led to a marked decrease in tumorigenic capabilities. Proteomic interactions involving EGFLAM, such as with NUP205, were implicated in cell cycle regulation, providing insight into its oncogenic mechanism. In vivo studies further demonstrated that silencing EGFLAM expression could inhibit tumor growth, underscoring its therapeutic potential. The study identifies EGFLAM as a pivotal oncogenic factor in GBM, serving as both a prognostic biomarker and a viable therapeutic target. These findings lay the groundwork for future research into EGFLAM-targeted therapies, aiming to improve clinical outcomes for GBM patients.
Subject(s)
Biomarkers, Tumor , Brain Neoplasms , Cell Proliferation , Glioblastoma , Humans , Glioblastoma/genetics , Glioblastoma/therapy , Glioblastoma/metabolism , Glioblastoma/pathology , Animals , Cell Line, Tumor , Biomarkers, Tumor/genetics , Biomarkers, Tumor/metabolism , Brain Neoplasms/genetics , Brain Neoplasms/therapy , Brain Neoplasms/metabolism , Prognosis , Female , Male , Gene Expression Regulation, Neoplastic , Mice , Mice, Nude , Middle Aged , Adult , Xenograft Model Antitumor Assays , Molecular Targeted TherapyABSTRACT
RATIONALE AND OBJECTIVES: Gliomas are aggressive brain tumors with a poor prognosis. Assessing treatment response is challenging because magnetic resonance imaging (MRI) may not distinguish true progression (TP) from pseudoprogression (PsP). This review aims to discuss imaging techniques and liquid biopsies used to distinguish TP from PsP. MATERIALS AND METHODS: This review synthesizes existing literature to examine advances in imaging techniques, such as magnetic resonance diffusion imaging (MRDI), perfusion-weighted imaging (PWI) MRI, and liquid biopsies, for identifying TP or PsP through tumor markers and tissue characteristics. RESULTS: Advanced imaging techniques, including MRDI and PWI MRI, have proven effective in delineating tumor tissue properties, offering valuable insights into glioma behavior. Similarly, liquid biopsy has emerged as a potent tool for identifying tumor-derived markers in biofluids, offering a non-invasive glimpse into tumor evolution. Despite their promise, these methodologies grapple with significant challenges. Their sensitivity remains inconsistent, complicating the accurate differentiation between TP and PSP. Furthermore, the absence of standardized protocols across platforms impedes the reliability of comparisons, while inherent biological variability adds complexity to data interpretation. CONCLUSION: Their potential applications have been highlighted, but gaps remain before routine clinical use. Further research is needed to develop and validate these promising methods for distinguishing TP from PsP in gliomas.
Subject(s)
Brain Neoplasms , Disease Progression , Glioma , Humans , Glioma/diagnostic imaging , Glioma/pathology , Liquid Biopsy/methods , Brain Neoplasms/diagnostic imaging , Brain Neoplasms/pathology , Diagnosis, Differential , Biomarkers, Tumor , Magnetic Resonance Imaging/methods , Diffusion Magnetic Resonance Imaging/methodsABSTRACT
Facial morphology, a complex trait influenced by genetics, holds great significance in evolutionary research. However, due to limited fossil evidence, the facial characteristics of Neanderthals and Denisovans have remained largely unknown. In this study, we conducted a large-scale multi-ethnic meta-analysis of the genome-wide association study (GWAS), including 9674 East Asians and 10,115 Europeans, quantitatively assessing 78 facial traits using 3D facial images. We identified 71 genomic loci associated with facial features, including 21 novel loci. We developed a facial polygenic score (FPS) that enables the prediction of facial features based on genetic information. Interestingly, the distribution of FPSs among populations from diverse continental groups exhibited relevant correlations with observed facial features. Furthermore, we applied the FPS to predict the facial traits of seven Neanderthals and one Denisovan using ancient DNA and aligned predictions with the fossil records. Our results suggested that Neanderthals and Denisovans likely shared similar facial features, such as a wider but shorter nose and a wider endocanthion distance. The decreased mouth width was characterized specifically in Denisovans. The integration of genomic data and facial trait analysis provides valuable insights into the evolutionary history and adaptive changes in human facial morphology.
ABSTRACT
Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.
Subject(s)
DNA Methylation , East Asian People , Quantitative Trait Loci , Female , Humans , Male , East Asian People/genetics , Genetic Predisposition to Disease , Genome-Wide Association Study/methods , Multifactorial Inheritance , Polymorphism, Single NucleotideABSTRACT
Background: Colon cancer (CC) is a highly heterogeneous malignancy associated with high morbidity and mortality. Pyroptosis is a type of programmed cell death characterized by an inflammatory response that can affect the tumor immune microenvironment and has potential prognostic and therapeutic value. The aim of this study was to evaluate the association between pyroptosis-related gene (PRG) expression and CC. Methods: Based on the expression profiles of PRGs, we classified CC samples from The Cancer Gene Atlas and Gene Expression Omnibus databases into different clusters by unsupervised clustering analysis. The best prognostic signature was screened and established using least absolute shrinkage and selection operator (LASSO) and multivariate COX regression analyses. Subsequently, a nomogram was established based on multivariate COX regression analysis. Next, gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were performed to explore the potential molecular mechanisms between the high- and low-risk groups and to explore the differences in clinicopathological characteristics, gene mutation characteristics, abundance of infiltrating immune cells, and immune microenvironment between the two groups. We also evaluated the association between common immune checkpoints and drug sensitivity using risk scores. The immunohistochemistry staining was utilized to confirm the expression of the selected genes in the prognostic model in CC. Results: The 1163 CC samples were divided into two clusters (clusters A and B) based on the expression profiles of the 33 PRGs. Genes with prognostic value were screened from the DEGs between the two clusters, and an eight PRGs prognostic model was constructed. GSEA and GSVA of the high- and low-risk groups revealed that they were mainly enriched in inflammatory response-related pathways. Compared to those in the low-risk group, patients in the high-risk group had worse overall survival, an immunosuppressive microenvironment, and worse sensitivity to immunotherapy and drug treatment. Conclusion: Our findings provide a foundation for future research targeting pyroptosis and new insights into prognosis and immunotherapy from the perspective of pyroptosis in CC.
ABSTRACT
With the irrational use of antibiotics and the increasing abuse of oral antibiotics, the drug resistance of gastrointestinal pathogens has become a prominent problem in clinical practice. Gut microbiota plays an important role in maintaining human health, and the change of microbiota also affects the activity of pathogenic bacteria. Interfering with antibiotic resistant bacteria by affecting gut microbiota has also become an important regulatory signal. In clinical application, due to the unique advantages of traditional Chinese medicine in sterilization and drug resistance, it is possible for traditional Chinese medicine to improve the gut microbial microenvironment. This review discusses the strategies of traditional Chinese medicine for the treatment of drug-resistant bacterial infections by changing the gut microenvironment, unlocking the interaction between microbiota and drug resistance of pathogenic bacteria.
Subject(s)
Gastrointestinal Microbiome , Microbiota , Humans , Anti-Bacterial Agents/pharmacology , Anti-Bacterial Agents/therapeutic use , Medicine, Chinese Traditional , BacteriaABSTRACT
Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.
ABSTRACT
Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363 K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.
Subject(s)
X-Ray Diffraction , Temperature , Crystallography , Binding Sites , Protein DomainsABSTRACT
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.
Subject(s)
Neanderthals , Humans , Animals , Mice , Neanderthals/genetics , Genome-Wide Association Study , Nose , Cell DifferentiationABSTRACT
Cryo-cooling has been nearly universally adopted to mitigate X-ray damage and facilitate crystal handling in protein X-ray crystallography. However, cryo X-ray crystallographic data provide an incomplete window into the ensemble of conformations that is at the heart of protein function and energetics. Room-temperature (RT) X-ray crystallography provides accurate ensemble information, and recent developments allow conformational heterogeneity (the experimental manifestation of ensembles) to be extracted from single-crystal data. Nevertheless, high sensitivity to X-ray damage at RT raises concerns about data reliability. To systematically address this critical issue, increasingly X-ray-damaged high-resolution data sets (1.02-1.52â Å resolution) were obtained from single proteinase K, thaumatin and lysozyme crystals at RT (277â K). In each case a modest increase in conformational heterogeneity with X-ray damage was observed. Merging data with different extents of damage (as is typically carried out) had negligible effects on conformational heterogeneity until the overall diffraction intensity decayed to â¼70% of its initial value. These effects were compared with X-ray damage effects in cryo-cooled crystals by carrying out an analogous analysis of increasingly damaged proteinase K cryo data sets (0.9-1.16â Å resolution). X-ray damage-associated heterogeneity changes were found that were not observed at RT. This property renders it difficult to distinguish real from artefactual conformations and to determine the conformational response to changes in temperature. The ability to acquire reliable heterogeneity information from single crystals at RT, together with recent advances in RT data collection at accessible synchrotron beamlines, provides a strong motivation for the widespread adoption of RT X-ray crystallography to obtain conformational ensemble information.
Subject(s)
Endopeptidase K/chemistry , Proteins , Crystallography, X-Ray , Proteins/chemistry , Reproducibility of Results , Temperature , X-RaysABSTRACT
Composite solid-state electrolytes (CSEs) have been developed rapidly in recent years owing to their high electrochemical stability, low cost, and easy processing characteristics. Most CSEs, however, require high temperatures or flammable liquid solvents to exhibit their acceptable electrochemical performance. Room-temperature all-solid-state batteries without liquid electrolytes are still unsatisfactory and under development. Herein, we have prepared a composite solid electrolyte with excellent performance using a polymer electrolyte poly(vinylidene fluoride-hexafluoropropylene) and an inorganic electrolyte Li6.4La3Zr1.4Ta0.6O12. With the assistance of lithium salts and plasticizers, the prepared CSE achieves a high ionic conductivity of 4.05 × 10-4 S·cm-1 at room temperature. The Li/CSE/Li symmetric cell can be stably cycled for more than 1000 h at 0.1 mA/cm2 without short circuits. The all-solid-state lithium metal battery using a LiFePO4 cathode displays a high discharge capacity of 148.1 mAh·g-1 and a capacity retention of 90.21% after 100 cycles. Moreover, the high electrochemical window up to 4.7 V of the CSE makes it suitable for high-voltage service environments. The all-solid-state battery using a lithium nickel-manganate cathode shows a high discharge specific capacity of 197.85 mAh·g-1 with good cycle performance. This work might guide the improvement of future CSEs and the exploration of flexible all-solid-state lithium metal batteries.
ABSTRACT
Vaccinated or not? This is an attitude survey for 'approach-avoidance conflict' under uncertainty. Therefore, measuring people's attitude toward COVID-19 vaccination is relatively distinctive from an attitude over a general conflict. An online survey of 3123 respondents from 30 provincial-level regions - out of 31 - on the Chinese mainland was conducted from January 22 to 27, 2021 to measure their willingness to be vaccinated. We found that over half of the respondents chose the options 'not to be vaccinated now' and 'wait and see before making a vaccination decision,' thereby indicating that people's willingness to be vaccinated is not as optimistic as anticipated in the early stage of vaccination in China. Hence, investigators should carefully select the measuring method to assess the 'true' levels of willingness to accept COVID-19 vaccines. Lastly, the relevant departments should fully predict obstacles to achieve immunity coverage and prepare for the 'vaccine hesitancy' of people in need.
Subject(s)
COVID-19 Vaccines , COVID-19 , COVID-19/prevention & control , Humans , SARS-CoV-2 , Surveys and Questionnaires , UncertaintyABSTRACT
Craniosynostosis (CRS) is a disease with prematurely fused cranial sutures. In the last decade, the whole-exome sequencing (WES) was widely used in Caucasian populations. The WES largely contributed in genetic diagnosis and exploration on new genetic mechanisms of CRS. In this study, we enrolled 264 CRS patients in China. After a 17-gene-panel sequencing designed in the previous study, 139 patients were identified with pathogenic/likely pathogenic (P/LP) variants according to the ACMG guideline as positive genetic diagnosis. WES was then performed on 102 patients with negative genetic diagnosis by panel. Ten P/LP variants were additionally identified in ten patients, increasing the genetic diagnostic yield by 3.8% (10/264). The novel variants in ANKH, H1-4, EIF5A, SOX6, and ARID1B expanded the mutation spectra of CRS. Then we designed a compatible research pipeline (RP) for further exploration. The RP could detect all seven P/LP SNVs and InDels identified above, in addition to 15 candidate variants found in 13 patients with worthy of further study. In sum, the 17-gene panel and WES identified positive genetic diagnosis for 56.4% patients (149/264) in 16 genes. At last, in our estimation, the genetic testing strategy of "Panel-first" saves 24.3% of the cost compared with "WES only", suggesting the "Panel-first" is an economical strategy.
ABSTRACT
Facial and cranial variation represent a multidimensional set of highly correlated and heritable phenotypes. Little is known about the genetic basis explaining this correlation. We develop a software package ALoSFL for simultaneous localization of facial and cranial landmarks from head computed tomography (CT) images, apply it in the analysis of head CT images of 777 Han Chinese women, and obtain a set of phenotypes representing variation in face, skull and facial soft tissue thickness (FSTT). Association analysis of 301 single nucleotide polymorphisms (SNPs) from 191 distinct genomic loci previously associated with facial variation reveals an unexpected larger number of loci showing significant associations (P < 1e-3) with cranial phenotypes than expected under the null (O/E = 3.39), suggesting facial and cranial phenotypes share a substantial proportion of genetic components. Adding FSTT to a SNP-only model shows a large impact in explaining facial variance. A gene ontology analysis reveals that bone morphogenesis and osteoblast differentiation likely underlie our cranial-significant findings. Overall, this study simultaneously investigates the genetic effects on both facial and cranial variation of the same sample, supporting that facial variation is a composite phenotype of cranial variation and FSTT.
Subject(s)
Face , Forensic Anthropology , Female , Animals , Face/anatomy & histology , Anatomic Landmarks , Skull/diagnostic imaging , Skull/anatomy & histology , PhenotypeABSTRACT
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.