A genome-wide association study identifies multiple loci for variation in human ear morphology

A multi-ancestry GWAS meta-analysis of facial features and its application in predicting archaic human features

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.

Human cytomegalovirus infection coopts chromatin organization to diminish TEAD1 transcription factor activity

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

Unveiling facial kinship: The BioKinVis dataset for facial kinship verification and genetic association studies

Facial image-based kinship verification represents a burgeoning frontier within the realms of computer vision and biomedicine. Recent genome-wide association studies have underscored the heritability of human facial morphology, revealing its predictability based on genetic information. These revelations form a robust foundation for advancing facial image-based kinship verification. Despite strides in computer vision, there remains a discernible gap between the biomedical and computer vision domains. Notably, the absence of family photo datasets established through biological paternity testing methods poses a significant challenge. This study addresses this gap by introducing the biological kinship visualization dataset, encompassing 5773 individuals from 2412 families with biologically confirmed kinship. Our analysis delves into the distribution and influencing factors of facial similarity among parent-child pairs, probing the potential association between forensic short tandem repeat polymorphisms and facial similarity. Additionally, we have developed a machine learning model for facial image-based kinship verification, achieving an accuracy of 0.80 in the dataset. To facilitate further exploration, we have established an online tool and database, accessible at http://120.55.161.230:88/.

Genome-wide analysis of spina bifida risk variants in a case-control study from Bangladesh

BACKGROUND

Human studies of genetic risk factors for neural tube defects, severe birth defects associated with long-term health consequences in surviving children, have predominantly been restricted to a subset of candidate genes in specific biological pathways including folate metabolism.

METHODS

In this study, we investigated the association of genetic variants spanning the genome with risk of spina bifida (i.e., myelomeningocele and meningocele) in a subset of families enrolled from December 2016 through December 2022 in a case-control study in Bangladesh, a population often underrepresented in genetic studies. Saliva DNA samples were analyzed using the Illumina Global Screening Array. We performed genetic association analyses to compare allele frequencies between 112 case and 121 control children, 272 mothers, and 128 trios.

RESULTS

In the transmission disequilibrium test analyses with trios only, we identified three novel exonic spina bifida risk loci, including rs140199800 (SULT1C2, p = 1.9 × 10-7), rs45580033 (ASB2, p = 4.2 × 10-10), and rs75426652 (LHPP, p = 7.2 × 10-14), after adjusting for multiple hypothesis testing. Association analyses comparing cases and controls, as well as models that included their mothers, did not identify genome-wide significant variants.

CONCLUSIONS

This study identified three novel single nucleotide polymorphisms involved in biological pathways not previously associated with neural tube defects. The study warrants replication in larger groups to validate findings and to inform targeted prevention strategies.

Selection-recombination-mutation dynamics: Gradient, limit cycle, and closed invariant curve

In this paper, the replicator dynamics of the two-locus two-allele system under weak mutation and weak selection is investigated in a generation-wise nonoverlapping unstructured population of individuals mating at random. Our main finding is that the dynamics is gradient-like when the point mutations at the two loci are independent. This is in stark contrast to the case of one-locus-multi-allele where the existence gradient behavior is contingent on a specific relationship between the mutation rates. When the mutations are not independent in the two-locus-two-allele system, there is the possibility of nonconvergent outcomes, like asymptotically stable oscillations, through either the Hopf bifurcation or the Neimark-Sacker bifurcation depending on the strength of the weak selection. The results can be straightforwardly extended for multilocus-two-allele systems.

Population diversity and equity in the genomic era: going global to return to the local

Advances in precision medicine depend on the quantity and quality of available genomic information. Various articles alert about the current disparities between the world's regions regarding the amount of genomic information available and the negative impact this will have on global health. The objective of this paper is to review these articles to describe what aspects they emphasize and highlight some issues that remain to be analyzed from the perspective of a "peripheral" country. Most of these articles come from central countries, where the need for more diversity in genomics is already detected. Several authors analyze lack of human diversity with focus on national, while others analyze the problem from a global perspective. Depending on the country of origin of the research, the claim for greater diversity has different meanings. Broadly, high-income countries advocate for better coverage looking within the boundaries of their own countries. In other regions of the world, where this field of research has not yet been massively developed, the same need for greater inclusiveness of origins in population genomics studies is not detected. An under-analyzed aspect is the unequal starting point between regions regarding the economic resources available for the development of this field of medicine, and for science and health in general. Although this macroeconomic and social aspect is usually absent in scientific analyses, without it solved, it will be impossible to guarantee that all world populations are equally represented in the panels or genomic databases that serve as input for precision medicine development.

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.

Evaluation of loci to predict ear morphology using two SNaPshot assays

Human ear morphology prediction with SNP-based genotypes is growing in forensic DNA phenotyping and is scarcely explored in Pakistan as a part of EVCs (externally visible characteristics). The ear morphology prediction assays with 21 SNPs were assessed for their potential utility in forensic identification of population. The SNaPshot™ multiplex chemistries, capillary electrophoresis methods and GeneMapper™ software were used for obtaining genotypic data. A total of 33 ear phenotypes were categorized with digital photographs of 300 volunteers. SHEsis software was applied to make LD plot. Ordinal and multinomial logistic regression was implemented for association testing. Multinomial logistic regression was executed to construct the prediction model in 90% training and 10% testing subjects. Several influential SNPs for ear phenotypic variation were found in association testing. The model based on genetic markers predicted ear phenotypes with moderate to good predictive accuracies demonstrated with the area under curve (AUC), sensitivity and specificity of predicted phenotypes. As an additional EVC, the estimated ear phenotypic profiles have the possibility of determining the human ear morphology differences in unknown biological samples found in crimes that do not result in a criminal database hit. Furthermore, this can help in facial reconstruction and act as an investigational lead.

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 risk factors for postoperative atrial fibrillation-a nationwide genome-wide association study (GWAS)

Background

Atrial fibrillation (AF) is a major cause of morbidity with a high prevalence among the elderly and has an established genetic disposition. Surgery is a well-known risk factor for AF; however, it is currently not recognized how much common genetic variants influence the postoperative risk. The purpose of this study was to identify Single Nucleotide Polymorphisms associated with postoperative AF.

Methods

The UK Biobank was utilized to conduct a Genome-Wide Association Study (GWAS) to identify variants associated with AF after surgery. An initial discovery GWAS was performed in patients that had undergone surgery with subsequent replication in a unique non-surgical cohort. In the surgical cohort, cases were defined as newly diagnosed AF within 30 days after surgery. The threshold for significance was set at 5 × 10^-8.

Results

After quality control, 144,196 surgical patients with 254,068 SNPs were left for analysis. Two variants (rs17042171 (p = 4.86 × 10^-15) and rs17042081 (p = 7.12 × 10^-15)) near the PITX2-gene reached statistical significance. These variants were replicated in the non-surgical cohort (1.39 × 10^-101 and 1.27 × 10^-93, respectively). Several other loci were significantly associated with AF in the non-surgical cohort.

Conclusion

In this GWAS-analysis of a large national biobank, we identified 2 variants that were significantly associated with postoperative AF. These variants were subsequently replicated in a unique non-surgical cohort. These findings bring new insight in the genetics of postoperative AF and may help identify at-risk patients and guide management.

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.

Functional Pathway and Process Enrichment Analysis of Genes Associated With Morphological Abnormalities of the Outer Ear

Congenital anomalies of the outer ear are common birth defects, including a variety of congenital deformities or malformations ranging from mild structural anomalies to total absence of the ear. Despite its high incidence and detrimental impact on patients, the etiology of outer ear abnormalities remains poorly understood. The goal of this study was to summarize the related genes and improve our understanding of the genetic etiology of morphological abnormalities of the outer ear. Human Phenotype Ontology (HPO) database, Mouse Genome Informatics (MGI) database, and PubMed search engine were used to acquire the genes associated with abnormal human or mouse outer ear. Metascape was employed on the genes above to conduct functional annotation, pathway and process enrichment analysis, protein-protein interaction network analysis, and MCODE component analysis. After a comprehensive review of the databases and literature, we identified 394 human genes and 148 mouse genes that have been associated with abnormal phenotypes of the outer ear, and we identified several biological pathways for human and mouse respectively. Especially, the analysis of common genes shared by human and mouse emphasized the importance of certain genes ( PAX6 , PBX1 , HOXA1 , HOXA2 , TBX1 , TBX15 , PRRX1 , and HMX1 ) in the embryonic development of the external ear. Through our analysis of genes associated with morphological abnormalities of the outer ear, the authors have shown that embryonic development pathways take important roles in the morphogenesis of abnormal external ear and highlighted some potential genetic drivers.

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.

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.

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.

Novel genetic associations with five aesthetic facial traits: A genome-wide association study in the Chinese population

Background: The aesthetic facial traits are closely related to life quality and strongly influenced by genetic factors, but the genetic predispositions in the Chinese population remain poorly understood.

Methods: A genome-wide association studies (GWAS) and subsequent validations were performed in 26,806 Chinese on five facial traits: widow’s peak, unibrow, double eyelid, earlobe attachment, and freckles. Functional annotation was performed based on the expression quantitative trait loci (eQTL) variants, genome-wide polygenic scores (GPSs) were developed to represent the combined polygenic effects, and single nucleotide polymorphism (SNP) heritability was presented to evaluate the contributions of the variants.

Results: In total, 21 genetic associations were identified, of which ten were novel: GMDS-AS1 (rs4959669, p = 1.29 × 10⁻⁴⁹) and SPRED2 (rs13423753, p = 2.99 × 10⁻¹⁴) for widow’s peak, a previously unreported trait; FARSB (rs36015125, p = 1.96 × 10⁻²¹) for unibrow; KIF26B (rs7549180, p = 2.41 × 10⁻¹⁵), CASC2 (rs79852633, p = 4.78 × 10⁻¹¹⁾, RPGRIP1L (rs6499632, p = 9.15 × 10⁻¹¹), and PAX1 (rs147581439, p = 3.07 × 10⁻⁸) for double eyelid; ZFHX3 (rs74030209, p = 9.77 × 10⁻¹⁴) and LINC01107 (rs10211400, p = 6.25 × 10⁻¹⁰) for earlobe attachment; and SPATA33 (rs35415928, p = 1.08 × 10⁻⁸) for freckles. Functionally, seven identified SNPs tag the missense variants and six may function as eQTLs. The combined polygenic effect of the associations was represented by GPSs and contributions of the variants were evaluated using SNP heritability.

Conclusion: These identifications may facilitate a better understanding of the genetic basis of features in the Chinese population and hopefully inspire further genetic research on facial development.

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.

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.

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.

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.

The impact of socioeconomic and phenotypic traits on self-perception of ethnicity in Latin America

Self-perception of ethnicity is a complex social trait shaped by both, biological and non-biological factors. We developed a comprehensive analysis of ethnic self-perception (ESP) on a large sample of Latin American mestizos from five countries, differing in age, socio-economic and education context, external phenotypic attributes and genetic background. We measured the correlation of ESP against genomic ancestry, and the influence of physical appearance, socio-economic context, and education on the distortion observed between both. Here we show that genomic ancestry is correlated to aspects of physical appearance, which in turn affect the individual ethnic self-perceived ancestry. Also, we observe that, besides the significant correlation among genomic ancestry and ESP, specific physical or socio-economic attributes have a strong impact on self-perception. In addition, the distortion among ESP and genomic ancestry differs across age ranks/countries, probably suggesting the underlying effect of past public policies regarding identity. Our results indicate that individuals' own ideas about its origins should be taken with caution, especially in aspects of modern life, including access to work, social policies, and public health key decisions such as drug administration, therapy design, and clinical trials, among others.

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.

Increased facial asymmetry in focal epilepsies associated with unilateral lesions

The epilepsies are now conceptualized as network disruptions: focal epilepsies are considered to have network alterations in the hemisphere of seizure onset, whilst generalized epilepsies are considered to have bi-hemispheric network changes. Increasingly, many epilepsies are also considered to be neurodevelopmental disorders, with early changes in the brain underpinning seizure biology. The development of the structure of the face is influenced by complex molecular interactions between surface ectoderm and underlying developing forebrain and neural crest cells. This influence is likely to continue postnatally, given the evidence of facial growth changes over time in humans until at least 18 years of age. In this case-control study, we hypothesized that people with lateralized focal epilepsies (i.e. unilateral network changes) have an increased degree of facial asymmetry, compared with people with generalized epilepsies or controls without epilepsy. We applied three-dimensional stereophotogrammetry and dense surface models to evaluate facial asymmetry in people with epilepsy, aiming to generate new tools to explore pathophysiological mechanisms in epilepsy. We analysed neuroimaging data to explore the correlation between face and brain asymmetry. We consecutively recruited 859 people with epilepsy attending the epilepsy clinics at a tertiary referral centre. We used dense surface modelling of the full face and signature analyses of three-dimensional facial photographs to analyse facial differences between 378 cases and 205 healthy controls. Neuroimaging around the time of the facial photograph was available for 234 cases. We computed the brain asymmetry index between contralateral regions. Cases with focal symptomatic epilepsy associated with unilateral lesions showed greater facial asymmetry compared to controls (P = 0.0001, two-sample t-test). This finding was confirmed by linear regression analysis after controlling for age and gender. We also found a significant correlation between duration of illness and the brain asymmetry index of total average cortical thickness (r = -0.19, P = 0.0075) but not for total average surface area (r = 0.06, P = 0.3968). There was no significant correlation between facial asymmetry and asymmetry of regional cortical thickness or surface area. We propose that the greater facial asymmetry in cases with focal epilepsy caused by unilateral abnormality might be explained by early unilateral network disruption, and that this is independent of underlying brain asymmetry. Three-dimensional stereophotogrammetry and dense surface modelling are a novel powerful phenotyping tool in epilepsy that may permit greater understanding of pathophysiology in epilepsy, and generate further insights into the development of cerebral networks underlying epilepsy, and the genetics of facial and neural development.

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.

Prediction of eye, hair and skin colour in Latin Americans

Here we evaluate the accuracy of prediction for eye, hair and skin pigmentation in a dataset of > 6500 individuals from Mexico, Colombia, Peru, Chile and Brazil (including genome-wide SNP data and quantitative/categorical pigmentation phenotypes - the CANDELA dataset CAN). We evaluated accuracy in relation to different analytical methods and various phenotypic predictors. As expected from statistical principles, we observe that quantitative traits are more sensitive to changes in the prediction models than categorical traits. We find that Random Forest or Linear Regression are generally the best performing methods. We also compare the prediction accuracy of SNP sets defined in the CAN dataset (including 56, 101 and 120 SNPs for eye, hair and skin colour prediction, respectively) to the well-established HIrisPlex-S SNP set (including 6, 22 and 36 SNPs for eye, hair and skin colour prediction respectively). When training prediction models on the CAN data, we observe remarkably similar performances for HIrisPlex-S and the larger CAN SNP sets for the prediction of hair (categorical) and eye (both categorical and quantitative), while the CAN sets outperform HIrisPlex-S for quantitative, but not for categorical skin pigmentation prediction. The performance of HIrisPlex-S, when models are trained in a world-wide sample (although consisting of 80% Europeans, https://hirisplex.erasmusmc.nl), is lower relative to training in the CAN data (particularly for hair and skin colour). Altogether, our observations are consistent with common variation of eye and hair colour having a relatively simple genetic architecture, which is well captured by HIrisPlex-S, even in admixed Latin Americans (with partial European ancestry). By contrast, since skin pigmentation is a more polygenic trait, accuracy is more sensitive to prediction SNP set size, although here this effect was only apparent for a quantitative measure of skin pigmentation. Our results support the use of HIrisPlex-S in the prediction of categorical pigmentation traits for forensic purposes in Latin America, while illustrating the impact of training datasets on its accuracy.

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.

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

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Minor Physical Anomalies in Bipolar Disorder

OBJECTIVE

High-arched palate is more frequent in schizophrenia and bipolar disorder (BD). Upto 40% of patients develop schizophrenia in 22q11.2 Deletion Syndrome manifested with cleft lip and palate, which originate from the first pharyngeal arch in embryo. The auricle also originates from the dorsal ends of the first and second pharyngeal arches; hence, we aimed to determine the associations between auricular anomalies and BD.

METHODS

We screened for 36 minor physical anomalies of the auricle in 146 patients with BD.

RESULTS

7 out of the of 36 assessed anomalies highly differed between healthy subjects and BD patients. A regression model including the differing anomalies predicted healthy subjects and BD-patients by 78.8% and 68.5%, respectively.

CONCLUSIONS

Assessing minor anomalies in psychiatric disorders may help to discover novel pathogenesis pathways and even new endophenotypes.

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues

Background

Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues.

Results

Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions.

Conclusions

The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

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.

Facial soft tissue thickness in Yangtze River delta Han population: Accurate assessment and comparative analysis utilizing Cone-Beam CT

Facial approximation plays a vital role in personal identification when other methods cannot be applied, and the study on facial soft tissue thickness (FSTT) is important to allow more accurate facial approximation. The Cone-Beam CT (CBCT) is considered to be an accurate approach in assessing FSTT, however only few studies have assessed the FSTT norms by CBCT in Chinese populations. The aim of this study was to (1) establish forensic measurement approach on a 3D model reconstructed by CBCT image; (2) obtain FSTT norms from the adults of Yangtze River delta Han population; (3) determine the impacts of sex, age, BMI as well as malocclusion type on FSTT and (4) form a FSTT database through CBCT. A total of 60 facial landmarks were primarily established from 3D reconstructed models from CBCT images of 424 healthy Han adults (130 males and 294 females) aged from 21 to 50 years in Yangtze River delta region and measurements were carried out. The multivariate analysis of variance test and multiple linear regressions test showed that the FSTT values on most landmarks were significantly different between sexes and BMI, and males showed greater values, expect for the landmarks distributed in malarzygomatic areas. Simultaneously, BMI mostly affected the landmarks located in zygomatic region. Besides, only a few landmarks were influenced by age. Interestingly, several periodontal landmarks were influenced by malocclusion types, thus the impact of malocclusion types on FSTT requires attention. Finally, the validation analysis suggested that measurement approach is an important factor for FSTT.

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.

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.

Genetic components of human pain sensitivity: a protocol for a genome-wide association study of experimental pain in healthy volunteers

INTRODUCTION

Pain constitutes a major component of the global burden of diseases. Recent studies suggest a strong genetic contribution to pain susceptibility and severity. Whereas most of the available evidence relies on candidate gene association or linkage studies, research on the genetic basis of pain sensitivity using genome-wide association studies (GWAS) is still in its infancy. This protocol describes a proposed GWAS on genetic contributions to baseline pain sensitivity and nociceptive sensitisation in a sample of unrelated healthy individuals of mixed Latin American ancestry.

METHODS AND ANALYSIS

A GWAS on genetic contributions to pain sensitivity in the naïve state and following nociceptive sensitisation will be conducted in unrelated healthy individuals of mixed ancestry. Mechanical and thermal pain sensitivity will be evaluated with a battery of quantitative sensory tests evaluating pain thresholds. In addition, variation in mechanical and thermal sensitisation following topical application of mustard oil to the skin will be evaluated.

ETHICS AND DISSEMINATION

This study received ethical approval from the University College London research ethics committee (3352/001) and from the bioethics committee of the Odontology Faculty at the University of Antioquia (CONCEPTO 01-2013). Findings will be disseminated to commissioners, clinicians and service users via papers and presentations at international conferences.

Comparison of methods for multivariate gene-based association tests for complex diseases using common variants

Complex diseases are usually associated with multiple correlated phenotypes, and the analysis of composite scores or disease status may not fully capture the complexity (or multidimensionality). Joint analysis of multiple disease-related phenotypes in genetic tests could potentially increase power to detect association of a disease with common SNPs (or genes). Gene-based tests are designed to identify genes containing multiple risk variants that individually are weakly associated with a univariate trait. We combined three multivariate association tests (O'Brien method, TATES, and MultiPhen) with two gene-based association tests (GATES and VEGAS) and compared performance (type I error and power) of six multivariate gene-based methods using simulated data. Data (n = 2000) for genetic sequence and correlated phenotypes were simulated by varying causal variant proportions and phenotype correlations for various scenarios. These simulations showed that two multivariate association tests (TATES and MultiPhen, but not O'Brien) paired with VEGAS have inflated type I error in all scenarios, while the three multivariate association tests paired with GATES have correct type I error. MultiPhen paired with GATES has higher power than competing methods if the correlations among phenotypes are low (r < 0.57). We applied these gene-based association methods to a GWAS dataset from the Alzheimer's Disease Genetics Consortium containing three neuropathological traits related to Alzheimer disease (neuritic plaque, neurofibrillary tangles, and cerebral amyloid angiopathy) measured in 3500 autopsied brains. Gene-level significant evidence (P < 2.7 × 10^-6) was identified in a region containing three contiguous genes (TRAPPC12, TRAPPC12-AS1, ADI1) using O'Brien and VEGAS. Gene-wide significant associations were not observed in univariate gene-based tests.

Latin Americans show wide-spread Converso ancestry and imprint of local Native ancestry on physical appearance

Historical records and genetic analyses indicate that Latin Americans trace their ancestry mainly to the intermixing (admixture) of Native Americans, Europeans and Sub-Saharan Africans. Using novel haplotype-based methods, here we infer sub-continental ancestry in over 6,500 Latin Americans and evaluate the impact of regional ancestry variation on physical appearance. We find that Native American ancestry components in Latin Americans correspond geographically to the present-day genetic structure of Native groups, and that sources of non-Native ancestry, and admixture timings, match documented migratory flows. We also detect South/East Mediterranean ancestry across Latin America, probably stemming mostly from the clandestine colonial migration of Christian converts of non-European origin (Conversos). Furthermore, we find that ancestry related to highland (Central Andean) versus lowland (Mapuche) Natives is associated with variation in facial features, particularly nose morphology, and detect significant differences in allele frequencies between these groups at loci previously associated with nose morphology in this sample.

Reevaluation of the earlobe types in Koreans

The shape of the earlobes has a variety of genetic significance. This study analyzed the frequencies of the earlobe shapes in the Korean population. Data were collected on randomly selected 500 males and 500 females in Daegu Metropolitan City, with all participant ages being in their twenties. Obtuse angled earlobes accounted for 41.2% of the earlobes observed, while acute angled earlobes prevalence was calculated at 38.8% and right angled earlobe was 20.0% of the total (sexes combined). In men, the acute angled earlobe was the most frequent type (43.0%), while the obtuse angled earlobe was the most frequent type in females (45.2%). These differences were statistically significant (p = 0.015). Overall, attached type earlobe (61.2%) was more frequent than free type earlobe. The attached type earlobe was more common in both sex groups (57.0% in male and 65.4% in female), and the proportion was significantly higher for females (p = 0.006). In conclusion, the findings in this study suggest that the attached earlobe type is the most common among Koreans, and the proportion of earlobe types among males and females is significantly different. Further studies are needed to understand the genetic background of earlobe types among Koreans.

Limited Evidence for Parallel Molecular Adaptations Associated with the Subterranean Niche in Mammals: A Comparative Study of Three Superorders

Among mammals, several lineages have independently adapted to a subterranean niche and possess similar phenotypic traits for burrowing (e.g., cylindrical bodies, short limbs, and absent pinnae). Previous research on mole-rats has revealed molecular adaptations for coping with reduced oxygen, elevated carbon dioxide, and the absence of light. In contrast, almost nothing is known regarding molecular adaptations in other subterranean lineages (e.g., true moles and golden moles). Therefore, the extent to which the recurrent phenotypic adaptations of divergent subterranean taxa have arisen via parallel routes of molecular evolution remains untested. To address these issues, we analyzed ∼8,000 loci in 15 representative subterranean taxa of four independent transitions to an underground niche for signatures of positive selection and convergent amino acid substitutions. Complementary analyses were performed in nonsubterranean "control" taxa to assess the biological significance of results. We found comparable numbers of positively selected genes in each of the four subterranean groups; however, correspondence in terms of gene identity between gene sets was low. Furthermore, we did not detect evidence of more convergent amino acids among subterranean species pairs compared with levels found between nonsubterranean controls. Comparisons with nonsubterranean taxa also revealed loci either under positive selection or with convergent substitutions, with similar functional enrichment (e.g., cell adhesion, immune response, and coagulation). Given the limited indication that positive selection and convergence occurred in the same loci, we conclude that selection may have acted on different loci across subterranean mammal lineages to produce similar phenotypes.

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.

Genome-wide association study in Japanese females identifies fifteen novel skin-related trait associations

Skin trait variation impacts quality-of-life, especially for females from the viewpoint of beauty. To investigate genetic variation related to these traits, we conducted a GWAS of various skin phenotypes in 11,311 Japanese women and identified associations for age-spots, freckles, double eyelids, straight/curly hair, eyebrow thickness, hairiness, and sweating. In silico annotation with RoadMap Epigenomics epigenetic state maps and colocalization analysis of GWAS and GTEx Project eQTL signals provided information about tissue specificity, candidate causal variants, and functional target genes. Novel signals for skin-spot traits neighboured AKAP1/MSI2 (rs17833789; P = 2.2 × 10-9), BNC2 (rs10810635; P = 2.1 × 10-22), HSPA12A (rs12259842; P = 7.1 × 10-11), PPARGC1B (rs251468; P = 1.3 × 10-21), and RAB11FIP2 (rs10444039; P = 5.6 × 10-21). HSPA12A SNPs were the only protein-coding gene eQTLs identified across skin-spot loci. Double edged eyelid analysis identified that a signal around EMX2 (rs12570134; P = 8.2 × 10-15) was also associated with expression of EMX2 and the antisense-RNA gene EMX2OS in brain putamen basal ganglia tissue. A known hair morphology signal in EDAR was associated with both eyebrow thickness (rs3827760; P = 1.7 × 10-9) and straight/curly hair (rs260643; P = 1.6 × 10-103). Excessive hairiness signals' top SNPs were also eQTLs for TBX15 (rs984225; P = 1.6 × 10-8), BCL2 (rs7226979; P = 7.3 × 10-11), and GCC2 and LIMS1 (rs6542772; P = 2.2 × 10-9). For excessive sweating, top variants in two signals in chr2:28.82-29.05 Mb (rs56089836; P = 1.7 × 10-11) were eQTLs for either PPP1CB or PLB1, while a top chr16:48.26-48.45 Mb locus SNP was a known ABCC11 missense variant (rs6500380; P = 6.8 × 10-10). In total, we identified twelve loci containing sixteen association signals, of which fifteen were novel. These findings will help dermatologic researchers better understand the genetic underpinnings of skin-related phenotypic variation in human populations.

Developmental pathways inferred from modularity, morphological integration and fluctuating asymmetry patterns in the human face

Facial asymmetries are usually measured and interpreted as proxies to developmental noise. However, analyses focused on its developmental and genetic architecture are scarce. To advance on this topic, studies based on a comprehensive and simultaneous analysis of modularity, morphological integration and facial asymmetries including both phenotypic and genomic information are needed. Here we explore several modularity hypotheses on a sample of Latin American mestizos, in order to test if modularity and integration patterns differ across several genomic ancestry backgrounds. To do so, 4104 individuals were analyzed using 3D photogrammetry reconstructions and a set of 34 facial landmarks placed on each individual. We found a pattern of modularity and integration that is conserved across sub-samples differing in their genomic ancestry background. Specifically, a signal of modularity based on functional demands and organization of the face is regularly observed across the whole sample. Our results shed more light on previous evidence obtained from Genome Wide Association Studies performed on the same samples, indicating the action of different genomic regions contributing to the expression of the nose and mouth facial phenotypes. Our results also indicate that large samples including phenotypic and genomic metadata enable a better understanding of the developmental and genetic architecture of craniofacial phenotypes.

Multiethnic GWAS Reveals Polygenic Architecture of Earlobe Attachment

The genetic basis of earlobe attachment has been a matter of debate since the early 20th century, such that geneticists argue both for and against polygenic inheritance. Recent genetic studies have identified a few loci associated with the trait, but large-scale analyses are still lacking. Here, we performed a genome-wide association study of lobe attachment in a multiethnic sample of 74,660 individuals from four cohorts (three with the trait scored by an expert rater and one with the trait self-reported). Meta-analysis of the three expert-rater-scored cohorts revealed six associated loci harboring numerous candidate genes, including EDAR, SP5, MRPS22, ADGRG6 (GPR126), KIAA1217, and PAX9. The large self-reported 23andMe cohort recapitulated each of these six loci. Moreover, meta-analysis across all four cohorts revealed a total of 49 significant (p < 5 × 10-8) loci. Annotation and enrichment analyses of these 49 loci showed strong evidence of genes involved in ear development and syndromes with auricular phenotypes. RNA sequencing data from both human fetal ear and mouse second branchial arch tissue confirmed that genes located among associated loci showed evidence of expression. These results provide strong evidence for the polygenic nature of earlobe attachment and offer insights into the biological basis of normal and abnormal ear development.