One of the features that distinguishes modern humans from our extinct relatives and ancestors is a globular shape of the braincase [1-4]. As the endocranium closely mirrors the outer shape of the brain, these differences might reflect altered neural architecture [4, 5]. However, in the absence of fossil brain tissue, the underlying neuroanatomical changes as well as their genetic bases remain elusive. To better understand the biological foundations of modern human endocranial shape, we turn to our closest extinct relatives: the Neandertals. Interbreeding between modern humans and Neandertals has resulted in introgressed fragments of Neandertal DNA in the genomes of present-day non-Africans [6, 7]. Based on shape analyses of fossil skull endocasts, we derive a measure of endocranial globularity from structural MRI scans of thousands of modern humans and study the effects of introgressed fragments of Neandertal DNA on this phenotype. We find that Neandertal alleles on chromosomes 1 and 18 are associated with reduced endocranial globularity. These alleles influence expression of two nearby genes, UBR4 and PHLPP1, which are involved in neurogenesis and myelination, respectively. Our findings show how integration of fossil skull data with archaic genomics and neuroimaging can suggest developmental mechanisms that may contribute to the unique modern human endocranial shape.

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.