Shared heritability of human face and brain shape

Loss of KAT6B causes premature ossification and promotes osteoblast differentiation during development

The MYST family histone acetyltransferase gene, KAT6B (MYST4, MORF, QKF) is mutated in two distinct human congenital disorders characterised by intellectual disability, facial dysmorphogenesis and skeletal abnormalities; the Say-Barber-Biesecker-Young-Simpson variant of Ohdo syndrome and Genitopatellar syndrome. Despite its requirement in normal skeletal development, the cellular and transcriptional effects of KAT6B in skeletogenesis have not been thoroughly studied. Here, we show that germline deletion of the Kat6b gene in mice causes premature ossification in vivo, resulting in shortened craniofacial elements and increased bone density, as well as shortened tibias with an expanded pre-hypertrophic layer, as compared to wild type controls. Mechanistically, we show that the loss of KAT6B in mesenchymal progenitor cells promotes transition towards an osteoblast-progenitor state with upregulation of gene targets of RUNX2, a master regulator of osteoblast development and concomitant downregulation of SOX9, a critical gene in chondrocyte development. Moreover, we find that compound heterozygosity at Kat6b and Runx2 loci partially rescues the reduction in ossification of Runx2 heterozygous, but not homozygous mice, suggesting that KAT6B may limit the action of RUNX2, possibly through a role in maintaining progenitors in an undifferentiated state. Moreover, our results show that KAT6B has essential roles in regulating the expression of a large number of genes involved in skeletogenesis and bone development.

Enhanced insights into the genetic architecture of 3D cranial vault shape using pleiotropy-informed GWAS

Large-scale GWAS studies have uncovered hundreds of genomic loci linked to facial and brain shape variation, but only tens associated with cranial vault shape, a largely overlooked aspect of the craniofacial complex. Surrounding the neocortex, the cranial vault plays a central role during craniofacial development and understanding its genetics are pivotal for understanding craniofacial conditions. Experimental biology and prior genetic studies have generated a wealth of knowledge that presents opportunities to aid further genetic discovery efforts. Here, we use the conditional FDR method to leverage GWAS data of facial shape, brain shape, and bone mineral density to enhance SNP discovery for cranial vault shape. This approach identified 120 independent genomic loci at 1% FDR, nearly tripling the number discovered through unconditioned analysis and implicating crucial craniofacial transcription factors and signaling pathways. These results significantly advance our genetic understanding of cranial vault shape and craniofacial development more broadly.

Optimized phenotyping of complex morphological traits: enhancing discovery of common and rare genetic variants

Genotype-phenotype (G-P) analyses for complex morphological traits typically utilize simple, predetermined anatomical measures or features derived via unsupervised dimension reduction techniques (e.g. principal component analysis (PCA) or eigen-shapes). Despite the popularity of these approaches, they do not necessarily reveal axes of phenotypic variation that are genetically relevant. Therefore, we introduce a framework to optimize phenotyping for G-P analyses, such as genome-wide association studies (GWAS) of common variants or rare variant association studies (RVAS) of rare variants. Our strategy is two-fold: (i) we construct a multidimensional feature space spanning a wide range of phenotypic variation, and (ii) within this feature space, we use an optimization algorithm to search for directions or feature combinations that are genetically enriched. To test our approach, we examine human facial shape in the context of GWAS and RVAS. In GWAS, we optimize for phenotypes exhibiting high heritability, estimated from either family data or genomic relatedness measured in unrelated individuals. In RVAS, we optimize for the skewness of phenotype distributions, aiming to detect commingled distributions that suggest single or few genomic loci with major effects. We compare our approach with eigen-shapes as baseline in GWAS involving 8246 individuals of European ancestry and in gene-based tests of rare variants with a subset of 1906 individuals. After applying linkage disequilibrium score regression to our GWAS results, heritability-enriched phenotypes yielded the highest SNP heritability, followed by eigen-shapes, while commingling-based traits displayed the lowest SNP heritability. Heritability-enriched phenotypes also exhibited higher discovery rates, identifying the same number of independent genomic loci as eigen-shapes with a smaller effective number of traits. For RVAS, commingling-based traits resulted in more genes passing the exome-wide significance threshold than eigen-shapes, while heritability-enriched phenotypes lead to only a few associations. Overall, our results demonstrate that optimized phenotyping allows for the extraction of genetically relevant traits that can specifically enhance discovery efforts of common and rare variants, as evidenced by their increased power in facial GWAS and RVAS.

Sex-stratified genome-wide meta-analysis identifies novel loci for cognitive decline in older adults

INTRODUCTION

Many complex traits and diseases show sex-specific biases in clinical presentation and prevalence.

METHODS

To understand sex-specific genetic architecture of cognitive decline across five cognitive domains (attention, memory, executive function, language, and visuospatial function) and global cognitive function, we performed sex-stratified genome-wide meta-analysis in 3021 older adults aged ≥ 65 years (female = 1545, male = 1476) from three prospective cohorts. Gene-based and gene-set enrichment analyses were conducted for each cognitive trait.

RESULTS

We identified a novel genome-wide significant (GWS) intergenic locus for decline of memory in males near RPS23P3 on chromosome 4 (rs6851574: minor allele frequency [MAF] = 0.39, Pmale = 4.10E-08, βmale = 0.19; Pinteraction = 3.76E-04). We also identified a subthreshold GWS locus for decline of executive function on chromosome 12 in females near the NDUFA12 gene, involved in oxidative phosphorylation (rs11107823: MAF = 0.12, Pfemale  = 9.35E-08, βfemale = 0.28; Pinteraction = 7.42E-06).

DISCUSSION

Sex-aware genetic studies can help in the identification of novel genetic loci and enhance sex-specific understanding of cognitive decline.

HIGHLIGHTS

Our genome-wide meta-analysis of single variants identified two new genetic associations, one in males and one in females. The novel male association was observed with the decline of memory in the intergenic region near the RPS23P3 gene on chromosome 4. This intergenic region has previously been implicated in several brain and cognition related traits, including anatomical brain aging, brain shape, and educational attainment. The novel female-specific association was observed with decline in executive function on chromosome 12 near the NDUFA12 gene, which is involved in oxidative phosphorylation. Sex-stratified analyses offer sex-specific understanding of biological pathways involved in cognitive aging.

GenomicSEM Modelling of Diverse Executive Function GWAS Improves Gene Discovery

Previous research has supported the use of latent variables as the gold-standard in measuring executive function. However, for logistical reasons genome-wide association studies (GWAS) of executive function have largely eschewed latent variables in favour of singular task measures. As low correlations have traditionally been found between individual executive function (EF) tests, it is unclear whether these GWAS have truly been measuring the same construct. In this study, we addressed this question by performing a factor analysis on summary statistics from eleven GWAS of EF taken from five studies, using GenomicSEM. Models demonstrated a bifactor structure consistent with previous research, with factors capturing common EF and working memory- specific variance. Furthermore, the GWAS performed on this model identified 20 new genomic risk loci for common EF and 4 for working memory reaching genome-wide significance beyond what was found in the constituent GWAS, together resulting in 29 newly mapped EF genes. These results help to clarify the underlying genetic structure of EF and support the idea that EF GWAS are capable of measuring genetic variance related to latent EF constructs even when not using factor scores. Furthermore, they demonstrate that GenomicSEM can combine GWAS with divergent and non-ideal measures of the same phenotype to improve statistical power.

The shared genetic architecture and evolution of human language and musical rhythm

This study aimed to test theoretical predictions over biological underpinnings of previously documented phenotypic correlations between human language-related and musical rhythm traits. Here, after identifying significant genetic correlations between rhythm, dyslexia and various language-related traits, we adapted multivariate methods to capture genetic signals common to genome-wide association studies of rhythm (N = 606,825) and dyslexia (N = 1,138,870). The results revealed 16 pleiotropic loci (P < 5 × 10-8) jointly associated with rhythm impairment and dyslexia, and intricate shared genetic and neurobiological architectures. The joint genetic signal was enriched for foetal and adult brain cell-specific regulatory regions, highlighting complex cellular composition in their shared underpinnings. Local genetic correlation with a key white matter tract (the left superior longitudinal fasciculus-I) substantiated hypotheses about auditory-motor connectivity as a genetically influenced, evolutionarily relevant neural endophenotype common to rhythm and language processing. Overall, we provide empirical evidence of multiple aspects of shared biology linking language and musical rhythm, contributing novel insight into the evolutionary relationships between human musicality and linguistic communication traits.

The X chromosome's influences on the human brain

Genes on the X chromosome are extensively expressed in the human brain. However, little is known for the X chromosome's impact on the brain anatomy, microstructure, and functional networks. We examined 1045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X chromosome interactions while proposing an atlas outlining dosage compensation for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we found unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Genetic variants associated with age-related episodic memory decline implicate distinct memory pathologies

BACKGROUND

Approximately 40% of people aged ≥ 65 experience memory loss, particularly in episodic memory. Identifying the genetic basis of episodic memory decline is crucial for uncovering its underlying causes.

METHODS

We investigated common and rare genetic variants associated with episodic memory decline in 742 (632 for rare variants) Ashkenazi Jewish individuals (mean age 75) from the LonGenity study. All-atom molecular dynamics simulations were performed to uncover mechanistic insights underlying rare variants associated with episodic memory decline.

RESULTS

In addition to the common polygenic risk of Alzheimer's disease, we identified and replicated rare variant associations in ITSN1 and CRHR2. Structural analyses revealed distinct memory pathologies mediated by interfacial rare coding variants such as impaired receptor activation of corticotropin releasing hormone and dysregulated L-serine synthesis.

DISCUSSION

Our study uncovers novel risk loci for episodic memory decline. The identified underlying mechanisms point toward heterogenous memory pathologies mediated by rare coding variants.

HIGHLIGHTS

We demonstrated the contribution of the common polygenic risk of Alzheimer's disease to episodic memory decline. We discovered and replicated two risk genes associated with episodic memory decline implicated by rare variants, were discovered and replicated. We demonstrated molecular mechanisms and potential novel memory pathologies underlying interfacial rare coding variants. Molecular dynamics simulations were performed to understand the downstream effects of risk rare coding variants.

X-chromosome-wide association study for Alzheimer's disease

Due to methodological reasons, the X-chromosome has not been featured in the major genome-wide association studies on Alzheimer's Disease (AD). To address this and better characterize the genetic landscape of AD, we performed an in-depth X-Chromosome-Wide Association Study (XWAS) in 115,841 AD cases or AD proxy cases, including 52,214 clinically-diagnosed AD cases, and 613,671 controls. We considered three approaches to account for the different X-chromosome inactivation (XCI) states in females, i.e. random XCI, skewed XCI, and escape XCI. We did not detect any genome-wide significant signals (P ≤ 5 × 10-8) but identified seven X-chromosome-wide significant loci (P ≤ 1.6 × 10-6). The index variants were common for the Xp22.32, FRMPD4, DMD and Xq25 loci, and rare for the WNK3, PJA1, and DACH2 loci. Overall, this well-powered XWAS found no genetic risk factors for AD on the non-pseudoautosomal region of the X-chromosome, but it identified suggestive signals warranting further investigations.

Syndrome-informed phenotyping identifies a polygenic background for achondroplasia-like facial variation in the general population

Human craniofacial shape is highly variable yet highly heritable with numerous genetic variants interacting through multiple layers of development. Here, we hypothesize that Mendelian phenotypes represent the extremes of a phenotypic spectrum and, using achondroplasia as an example, we introduce a syndrome-informed phenotyping approach to identify genomic loci associated with achondroplasia-like facial variation in the general population. We compare three-dimensional facial scans from 43 individuals with achondroplasia and 8246 controls to calculate achondroplasia-like facial scores. Multivariate GWAS of the control scores reveals a polygenic basis for facial variation along an achondroplasia-specific shape axis, identifying genes primarily involved in skeletal development. Jointly modeling these genes in two independent control samples, both human and mouse, shows craniofacial effects approximating the characteristic achondroplasia phenotype. These findings suggest that both complex and Mendelian genetic variation act on the same developmentally determined axes of facial variation, providing insights into the genetic intersection of complex traits and Mendelian disorders.

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

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

Gene-environment interactions in the influence of maternal education on adolescent neurodevelopment using ABCD study

Maternal education was strongly correlated with adolescent brain morphology, cognitive performances, and mental health. However, the molecular basis for the effects of maternal education on the structural neurodevelopment remains unknown. Here, we conducted gene-environment-wide interaction study using the Adolescent Brain Cognitive Development cohort. Seven genomic loci with significant gene-environment interactions (G×E) on regional gray matter volumes were identified, with enriched biological functions related to metabolic process, inflammatory process, and synaptic plasticity. Additionally, genetic overlapping results with behavioral and disease-related phenotypes indicated shared biological mechanism between maternal education modified neurodevelopment and related behavioral traits. Finally, by decomposing the multidimensional components of maternal education, we found that socioeconomic status, rather than family environment, played a more important role in modifying the genetic effects on neurodevelopment. In summary, our study provided analytical evidence for G×E effects regarding adolescent neurodevelopment and explored potential biological mechanisms as well as social mechanisms through which maternal education could modify the genetic effects on regional brain development.

A novel classification framework for genome-wide association study of whole brain MRI images using deep learning

Genome-wide association studies (GWASs) have been widely applied in the neuroimaging field to discover genetic variants associated with brain-related traits. So far, almost all GWASs conducted in neuroimaging genetics are performed on univariate quantitative features summarized from brain images. On the other hand, powerful deep learning technologies have dramatically improved our ability to classify images. In this study, we proposed and implemented a novel machine learning strategy for systematically identifying genetic variants that lead to detectable nuances on Magnetic Resonance Images (MRI). For a specific single nucleotide polymorphism (SNP), if MRI images labeled by genotypes of this SNP can be reliably distinguished using machine learning, we then hypothesized that this SNP is likely to be associated with brain anatomy or function which is manifested in MRI brain images. We applied this strategy to a catalog of MRI image and genotype data collected by the Alzheimer's Disease Neuroimaging Initiative (ADNI) consortium. From the results, we identified novel variants that show strong association to brain phenotypes.

Loss of normal facial asymmetry in schizophrenia and bipolar disorder: Implications for development of brain asymmetry in psychotic illness

Development of the craniofacies occurs in embryological intimacy with development of the brain and both show normal left-right asymmetries. While facial dysmorphology occurs to excess in psychotic illness, facial asymmetry has yet to be investigated as a putative index of brain asymmetry. Ninety-three subjects (49 controls, 22 schizophrenia, 22 bipolar disorder) received 3D laser surface imaging of the face. On geometric morphometric analysis with (x, y, z) visualisations of statistical models for facial asymmetries, in controls the upper face and periorbital region, which share embryological intimacy with the forebrain, showed marked asymmetries. Their geometry included: along the x-axis, rightward asymmetry in its dorsal-medial aspects and leftward asymmetry in its ventral-lateral aspects; along the z-axis, anterior protrusion in its right ventral-lateral aspect. In both schizophrenia and bipolar disorder these normal facial asymmetries were diminished, with residual retention of asymmetries in bipolar disorder. This geometry of normal facial asymmetries shows commonalities with that of normal frontal lobe asymmetries. These findings indicate a trans-diagnostic process that involves loss of facial asymmetries in both schizophrenia and bipolar disorder. Embryologically, they implicate loss of face-brain asymmetries across gestational weeks 7-14 in processes that involve genes previously associated with risk for schizophrenia.

Face-brain correlates as potential sex-specific biomarkers for schizophrenia and bipolar disorder

Given the shared ectodermal origin and integrated development of the face and the brain, facial biomarkers emerge as potential candidates to assess vulnerability for disorders in which neurodevelopment is compromised, such as schizophrenia (SZ) and bipolar disorder (BD). The sample comprised 188 individuals (67 SZ patients, 46 BD patients and 75 healthy controls (HC)). Using a landmark-based approach on 3D facial reconstructions, we quantified global and local facial shape differences between SZ/BD patients and HC using geometric morphometrics. We also assessed correlations between facial and brain cortical measures. All analyses were performed separately by sex. Diagnosis explained 4.1 % - 5.9 % of global facial shape variance in males and females with SZ, and 4.5 % - 4.1 % in BD. Regarding local facial shape, we detected 43.2 % of significantly different distances in males and 47.4 % in females with SZ as compared to HC, whereas in BD the percentages decreased to 35.8 % and 26.8 %, respectively. We detected that brain area and volume significantly explained 2.2 % and 2 % of facial shape variance in the male SZ - HC sample. Our results support facial shape as a neurodevelopmental marker for SZ and BD and reveal sex-specific pathophysiological mechanisms modulating the interplay between the brain and the face.

Genetic architecture of brain morphology and overlap with neuropsychiatric traits

Uncovering the genetic architectures of brain morphology offers valuable insights into brain development and disease. Genetic association studies of brain morphological phenotypes have discovered thousands of loci. However, interpretation of these loci presents a significant challenge. One potential solution is exploring the genetic overlap between brain morphology and disorders, which can improve our understanding of their complex relationships, ultimately aiding in clinical applications. In this review, we examine current evidence on the genetic associations between brain morphology and neuropsychiatric traits. We discuss the impact of these associations on the diagnosis, prediction, and treatment of neuropsychiatric diseases, along with suggestions for future research directions.

Discovering non-additive heritability using additive GWAS summary statistics

LD score regression (LDSC) is a method to estimate narrow-sense heritability from genome-wide association study (GWAS) summary statistics alone, making it a fast and popular approach. In this work, we present interaction-LD score (i-LDSC) regression: an extension of the original LDSC framework that accounts for interactions between genetic variants. By studying a wide range of generative models in simulations, and by re-analyzing 25 well-studied quantitative phenotypes from 349,468 individuals in the UK Biobank and up to 159,095 individuals in BioBank Japan, we show that the inclusion of a cis-interaction score (i.e. interactions between a focal variant and proximal variants) recovers genetic variance that is not captured by LDSC. For each of the 25 traits analyzed in the UK Biobank and BioBank Japan, i-LDSC detects additional variation contributed by genetic interactions. The i-LDSC software and its application to these biobanks represent a step towards resolving further genetic contributions of sources of non-additive genetic effects to complex trait variation.

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

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

Author Summary

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

Novel early-onset Alzheimer-associated genes influence risk through dysregulation of glutamate, immune activation, and intracell signaling pathways

Alzheimer Disease (AD) is a highly polygenic disease that presents with relatively earlier onset (≤70yo; EOAD) in about 5% of cases. Around 90% of these EOAD cases remain unexplained by pathogenic mutations. Using data from EOAD cases and controls, we performed a genome-wide association study (GWAS) and trans-ancestry meta-analysis on non-Hispanic Whites (NHW, NCase=6,282, NControl=13,386), African Americans (AA NCase=782, NControl=3,663) and East Asians (NCase=375, NControl=838 CO). We identified eight novel significant loci: six in the ancestry-specific analyses and two in the trans-ancestry analysis. By integrating gene-based analysis, eQTL, pQTL and functional annotations, we nominate four novel genes that are involved in microglia activation, glutamate production, and signaling pathways. These results indicate that EOAD, although sharing many genes with LOAD, harbors unique genes and pathways that could be used to create better prediction models or target identification for this type of AD.

Comparative analysis of craniofacial shape in two mouse models of Down syndrome: Ts65Dn and TcMAC21

Mouse models are central to studying and understanding the genotypic-to-phenotypic outcomes of Down syndrome (DS), a complex condition caused by an extra copy of the long arm of human chromosome 21. The recently developed TcMAC21-a transchromosomic mouse strain with comparable gene dosage to human chromosome 21 (Hsa21)-includes more Hsa21 genes than any other model of DS. Recent studies on TcMAC21 have provided valuable insight into the molecular, physiological, and neuroanatomical aspects of the model. However, relatively little is known about the craniofacial phenotype of TcMAC21 mice, particularly as it compares to the widely studied Ts65Dn model. Here we conducted a quantitative study of the cranial morphology of TcMAC21 and Ts65Dn mice and their respective unaffected littermates. Our comparative data comprise forty three-dimensional cranial measurements taken on micro-computed tomography scans of the heads of TcMAC21 and Ts65Dn mice. Our results show that TcMAC21 exhibit similar patterns of craniofacial change to Ts65Dn. However, the DS-specific morphology is more pronounced in Ts65Dn mice. Specifically, Ts65Dn present with more medio-lateral broadening and retraction of the snout compared to TcMAC21. Our findings reveal the complexity of potential gene interaction in the production of craniofacial phenotypes.

Rare genetic coding variants associated with age-related episodic memory decline implicate distinct memory pathologies in the hippocampus

Background

Approximately 40% of people aged 65 or older experience memory loss, particularly in episodic memory. Identifying the genetic basis of episodic memory decline is crucial for uncovering its underlying causes.

Methods

We investigated common and rare genetic variants associated with episodic memory decline in 742 (632 for rare variants) Ashkenazi Jewish individuals (mean age 75) from the LonGenity study. All-atom MD simulations were performed to uncover mechanistic insights underlying rare variants associated with episodic memory decline.

Results

In addition to the common polygenic risk of Alzheimer's Disease (AD), we identified and replicated rare variant association in ITSN1 and CRHR2 . Structural analyses revealed distinct memory pathologies mediated by interfacial rare coding variants such as impaired receptor activation of corticotropin releasing hormone and dysregulated L-serine synthesis.

Discussion

Our study uncovers novel risk loci for episodic memory decline. The identified underlying mechanisms point toward heterogeneous memory pathologies mediated by rare coding variants.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Regulatory elements in SEM1-DLX5-DLX6 (7q21.3) locus contribute to genetic control of coronal nonsyndromic craniosynostosis and bone density-related traits

Purpose

The etiopathogenesis of coronal nonsyndromic craniosynostosis (cNCS), a congenital condition defined by premature fusion of 1 or both coronal sutures, remains largely unknown.

Methods

We conducted the largest genome-wide association study of cNCS followed by replication, fine mapping, and functional validation of the most significant region using zebrafish animal model.

Results

Genome-wide association study identified 6 independent genome-wide-significant risk alleles, 4 on chromosome 7q21.3 SEM1-DLX5-DLX6 locus, and their combination conferred over 7-fold increased risk of cNCS. The top variants were replicated in an independent cohort and showed pleiotropic effects on brain and facial morphology and bone mineral density. Fine mapping of 7q21.3 identified a craniofacial transcriptional enhancer (eDlx36) within the linkage region of the top variant (rs4727341; odds ratio [95% confidence interval], 0.48[0.39-0.59]; P = 1.2E-12) that was located in SEM1 intron and enriched in 4 rare risk variants. In zebrafish, the activity of the transfected human eDlx36 enhancer was observed in the frontonasal prominence and calvaria during skull development and was reduced when the 4 rare risk variants were introduced into the sequence.

Conclusion

Our findings support a polygenic nature of cNCS risk and functional role of craniofacial enhancers in cNCS susceptibility with potential broader implications for bone health.

A common cis-regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis

Genome-wide association studies (GWAS) identified thousands of genetic variants linked to phenotypic traits and disease risk. However, mechanistic understanding of how GWAS variants influence complex morphological traits and can, in certain cases, simultaneously confer normal-range phenotypic variation and disease predisposition, is still largely lacking. Here, we focus on rs6740960, a single nucleotide polymorphism (SNP) at the 2p21 locus, which in GWAS studies has been associated both with normal-range variation in jaw shape and with an increased risk of non-syndromic orofacial clefting. Using in vitro derived embryonic cell types relevant for human facial morphogenesis, we show that this SNP resides in an enhancer that regulates chondrocytic expression of PKDCC - a gene encoding a tyrosine kinase involved in chondrogenesis and skeletal development. In agreement, we demonstrate that the rs6740960 SNP is sufficient to confer chondrocyte-specific differences in PKDCC expression. By deploying dense landmark morphometric analysis of skull elements in mice, we show that changes in Pkdcc dosage are associated with quantitative changes in the maxilla, mandible, and palatine bone shape that are concordant with the facial phenotypes and disease predisposition seen in humans. We further demonstrate that the frequency of the rs6740960 variant strongly deviated among different human populations, and that the activity of its cognate enhancer diverged in hominids. Our study provides a mechanistic explanation of how a common SNP can mediate normal-range and disease-associated morphological variation, with implications for the evolution of human facial features.

Identification and validation of supervariants reveal novel loci associated with human white matter microstructure

As an essential part of the central nervous system, white matter coordinates communications between different brain regions and is related to a wide range of neurodegenerative and neuropsychiatric disorders. Previous genome-wide association studies (GWASs) have uncovered loci associated with white matter microstructure. However, GWASs suffer from limited reproducibility and difficulties in detecting multi-single-nucleotide polymorphism (multi-SNP) and epistatic effects. In this study, we adopt the concept of supervariants, a combination of alleles in multiple loci, to account for potential multi-SNP effects. We perform supervariant identification and validation to identify loci associated with 22 white matter fractional anisotropy phenotypes derived from diffusion tensor imaging. To increase reproducibility, we use United Kingdom (UK) Biobank White British (n = 30,842) data for discovery and internal validation, and UK Biobank White but non-British (n = 1927) data, Europeans from the Adolescent Brain Cognitive Development study (n = 4399) data, and Europeans from the Human Connectome Project (n = 319) data for external validation. We identify 23 novel loci on the discovery set that have not been reported in the previous GWASs on white matter microstructure. Among them, three supervariants on genomic regions 5q35.1, 8p21.2, and 19q13.32 have P-values lower than 0.05 in the meta-analysis of the three independent validation data sets. These supervariants contain genetic variants located in genes that have been related to brain structures, cognitive functions, and neuropsychiatric diseases. Our findings provide a better understanding of the genetic architecture underlying white matter microstructure.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

Syndrome-informed phenotyping identifies a polygenic background for achondroplasia-like facial variation in the general population

Human craniofacial shape is highly variable yet highly heritable with genetic variants interacting through multiple layers of development. Here, we hypothesize that Mendelian phenotypes represent the extremes of a phenotypic spectrum and, using achondroplasia as an example, we introduce a syndrome-informed phenotyping approach to identify genomic loci associated with achondroplasia-like facial variation in the normal population. We compared three-dimensional facial scans from 43 individuals with achondroplasia and 8246 controls to calculate achondroplasia-like facial scores. Multivariate GWAS of the control scores revealed a polygenic basis for normal facial variation along an achondroplasia-specific shape axis, identifying genes primarily involved in skeletal development. Jointly modeling these genes in two independent control samples showed craniofacial effects approximating the characteristic achondroplasia phenotype. These findings suggest that both complex and Mendelian genetic variation act on the same developmentally determined axes of facial variation, providing new insights into the genetic intersection of complex traits and Mendelian disorders.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

X chromosome-wide association study of quantitative biomarkers from the Alzheimer's Disease Neuroimaging Initiative study

Introduction

Alzheimer's disease (AD) is a complex neurodegenerative disease with high heritability. Compared to autosomes, a higher proportion of disorder-associated genes on X chromosome are expressed in the brain. However, only a few studies focused on the identification of the susceptibility loci for AD on X chromosome.

Methods

Using the data from the Alzheimer's Disease Neuroimaging Initiative Study, we conducted an X chromosome-wide association study between 16 AD quantitative biomarkers and 19,692 single nucleotide polymorphisms (SNPs) based on both the cross-sectional and longitudinal studies.

Results

We identified 15 SNPs statistically significantly associated with different quantitative biomarkers of the AD. For the cross-sectional study, six SNPs (rs5927116, rs4596772, rs5929538, rs2213488, rs5920524, and rs5945306) are located in or near to six genes DMD, TBX22, LOC101928437, TENM1, SPANXN1, and ZFP92, which have been reported to be associated with schizophrenia or neuropsychiatric diseases in literature. For the longitudinal study, four SNPs (rs4829868, rs5931111, rs6540385, and rs763320) are included in or near to two genes RAC1P4 and AFF2, which have been demonstrated to be associated with brain development or intellectual disability in literature, while the functional annotations of other five novel SNPs (rs12157031, rs428303, rs5953487, rs10284107, and rs5955016) have not been found.

Discussion

15 SNPs were found statistically significantly associated with the quantitative biomarkers of the AD. Follow-up study in molecular genetics is needed to verify whether they are indeed related to AD. The findings in this article expand our understanding of the role of the X chromosome in exploring disease susceptibility, introduce new insights into the molecular genetics behind the AD, and may provide a mechanistic clue to further AD-related studies.

The shared genetic architecture and evolution of human language and musical rhythm

Rhythm and language-related traits are phenotypically correlated, but their genetic overlap is largely unknown. Here, we leveraged two large-scale genome-wide association studies performed to shed light on the shared genetics of rhythm (N=606,825) and dyslexia (N=1,138,870). Our results reveal an intricate shared genetic and neurobiological architecture, and lay groundwork for resolving longstanding debates about the potential co-evolution of human language and musical traits.

Genome-wide association study of hippocampal blood-oxygen-level-dependent-cerebral blood flow correlation in Chinese Han population

Correlation between blood-oxygen-level-dependent (BOLD) and cerebral blood flow (CBF) has been used as an index of neurovascular coupling. Hippocampal BOLD-CBF correlation is associated with neurocognition, and the reduced correlation is associated with neuropsychiatric disorders. We conducted the first genome-wide association study of the hippocampal BOLD-CBF correlation in 4,832 Chinese Han subjects. The hippocampal BOLD-CBF correlation had an estimated heritability of 16.2-23.9% and showed reliable genome-wide significant association with a locus at 3q28, in which many variants have been linked to neuroimaging and cerebrospinal fluid markers of Alzheimer's disease. Gene-based association analyses showed four significant genes (GMNC, CRTC2, DENND4B, and GATAD2B) and revealed enrichment for mast cell calcium mobilization, microglial cell proliferation, and ubiquitin-related proteolysis pathways that regulate different cellular components of the neurovascular unit. This is the first unbiased identification of the association of hippocampal BOLD-CBF correlation, providing fresh insights into the genetic architecture of hippocampal neurovascular coupling.

Analyzing Facial Asymmetry in Alzheimer's Dementia Using Image-Based Technology

Several studies have demonstrated accelerated brain aging in Alzheimer's dementia (AD). Previous studies have also reported that facial asymmetry increases with age. Because obtaining facial images is much easier than obtaining brain images, the aim of this work was to investigate whether AD exhibits accelerated aging patterns in facial asymmetry. We developed new facial asymmetry measures to compare Alzheimer's patients with healthy controls. A three-dimensional camera was used to capture facial images, and 68 facial landmarks were identified using an open-source machine-learning algorithm called OpenFace. A standard image registration method was used to align the three-dimensional original and mirrored facial images. This study used the registration error, representing landmark superimposition asymmetry distances, to examine 29 pairs of landmarks to characterize facial asymmetry. After comparing the facial images of 150 patients with AD with those of 150 age- and sex-matched non-demented controls, we found that the asymmetry of 20 landmarks was significantly different in AD than in the controls (p < 0.05). The AD-linked asymmetry was concentrated in the face edge, eyebrows, eyes, nostrils, and mouth. Facial asymmetry evaluation may thus serve as a tool for the detection of AD.

CEGA: a method for inferring natural selection by comparative population genomic analysis across species

We developed maximum likelihood method for detecting positive selection or balancing selection using multilocus or genomic polymorphism and divergence data from two species. The method is especially useful for investigating natural selection in noncoding regions. Simulations demonstrate that the method outperforms existing methods in detecting both positive and balancing selection. We apply the method to population genomic data from human and chimpanzee. The list of genes identified under selection in the noncoding regions is prominently enriched in pathways related to the brain and nervous system. Therefore, our method will serve as a useful tool for comparative population genomic analysis.

Postnatal growth and spatial conformity of the cranium, brain, eyeballs and masseter muscles in the macaque (Macaca mulatta)

Spatial growth constraints in the head region can lead to coordinated patterns of morphological variation that pleiotropically modify genetically defined phenotypes as the tissues compete for space. Here we test for such architectural modifications during rhesus macaque (Macaca mulatta) postnatal ontogeny. We captured cranium and brain shape from 153 MRI datasets spanning 13 to 1090 postnatal days and tested for patterns of covariation with measurements of relative brain, eyeball, and masseter muscle size as well as callosal tract length. We find that the shape of the infant (<365 days) macaque cranium was most closely aligned to masseter muscle and brain size measured relative to face size. Infant brain and juvenile (365-1090 days) cranium shape were more closely linked with brain size relative to basicranium and face size. Meanwhile, the juvenile macaque brain shape was dominated by the size of the brain relative to that of the basicranium. Associations with relative eyeball size and commissural tract lengths were weaker. Our results are consistent with a spatial-packing regime operating during postnatal macaque ontogeny, in which relative growth of the masseter, face and basicranium have a greater influence than brain growth on the overall shape of the cranium and brain.

Transmit Radiant Individuality to Offspring (TRIO) study: investigating intergenerational transmission effects on brain development

Intergenerational transmission is a crucial aspect of human development. Although prior studies have demonstrated the continuity of psychopathology and maladaptive upbringing environments between parents and offspring, the underlying neurobiological mechanisms remain unclear. We have begun a novel neuroimaging research project, the Transmit Radiant Individuality to Offspring (TRIO) study, which focuses on biological parent-offspring trios. The participants of the TRIO study were Japanese parent-offspring trios consisting of offspring aged 10-40 and their biological mother and father. Structural and functional brain images of all participants were acquired using magnetic resonance imaging (MRI). Saliva samples were collected for DNA analysis. We obtained psychosocial information, such as intelligence, mental health problems, personality traits, and experiences during the developmental period from each parent and offspring in the same manner as much as possible. By April 2023, we completed data acquisition from 174 trios consisting of fathers, mothers, and offspring. The target sample size was 310 trios. However, we plan to conduct genetic and epigenetic analyses, and the sample size is expected to be expanded further while developing this project into a multi-site collaborative study in the future. The TRIO study can challenge the elucidation of the mechanism of intergenerational transmission effects on human development by collecting diverse information from parents and offspring at the molecular, neural, and behavioral levels. Our study provides interdisciplinary insights into how individuals' lives are involved in the construction of the lives of their descendants in the subsequent generation.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

Genetic insights into human cortical organization and development through genome-wide analyses of 2,347 neuroimaging phenotypes

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

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.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Genotype Data and Derived Genetic Instruments of Adolescent Brain Cognitive Development Study® for Better Understanding of Human Brain Development

The data release of Adolescent Brain Cognitive Development® (ABCD) Study represents an extensive resource for investigating factors relating to child development and mental wellbeing. The genotype data of ABCD has been used extensively in the context of genetic analysis, including genome-wide association studies and polygenic score predictions. However, there are unique opportunities provided by ABCD genetic data that have not yet been fully tapped. The diverse genomic variability, the enriched relatedness among ABCD subsets, and the longitudinal design of the ABCD challenge researchers to perform novel analyses to gain deeper insight into human brain development. Genetic instruments derived from the ABCD genetic data, such as genetic principal components, can help to better control confounds beyond the context of genetic analyses. To facilitate the use genomic information in the ABCD for inference, we here detail the processing procedures, quality controls, general characteristics, and the corresponding resources in the ABCD genotype data of release 4.0.

Microvascular and macrovascular complications of type 2 diabetes mellitus: Exome wide association analyses

Background

Type 2 diabetes mellitus (T2DM) is a chronic, metabolic disorder in which concomitant insulin resistance and β-cell impairment lead to hyperglycemia, influenced by genetic and environmental factors. T2DM is associated with long-term complications that have contributed to the burden of morbidity and mortality worldwide. The objective of this manuscript is to conduct an Exome-Wide Association Study (EWAS) on T2DM Emirati individuals to improve our understanding on diabetes-related complications to improve early diagnostic methods and treatment strategies.

Methods

This cross-sectional study recruited 310 Emirati participants that were stratified according to their medically diagnosed diabetes-related complications: diabetic retinopathy, diabetic neuropathy, diabetic nephropathy, and cardiovascular complications. The Illumina's Infinium Exome-24 array was used and 39,840 SNPs remained for analysis after quality control.

Findings

The analysis revealed the associations of various genes with each complication category: 1) diabetic retinopathy was associated to SHANK3 gene in locus 22q13.33 (SNP rs9616915; p=5.18 x10-4), ZSCAN5A gene in locus 19q13.43 (SNP rs7252603; p=7.55 x10-4), and DCP1B gene in locus 12p13.33 (SNPs rs715146, rs1044950, rs113147414, rs34730825; p=7.62 x10-4); 2) diabetic neuropathy was associated to ADH4 gene in locus 4q23 (SNP rs4148883; p=1.23 x10-4), SLC11A1 gene in locus 2q35 (SNP rs17235409; p=1.85 x10-4), and MATN4 gene in locus 20q13.12 (SNP rs2072788; p=2.68 x10-4); 3) diabetic nephropathy was associated to PPP1R3A gene in locus 7q31.1 (SNP rs1799999; p=1.91 x10-4), ZNF136 gene in locus 19p13.2 (SNP rs140861589; p=2.80 x10-4), and HSPA12B gene in locus 20p13 (SNP rs6076550; p=2.86 x10-4); and 4) cardiovascular complications was associated to PCNT gene in locus 21q22.3 (SNPs rs7279204, rs6518289, rs2839227, rs2839223; p=2.18 x10-4,3.04 x10-4,4.51 x10-4,5.22 x10-4 respectively), SEPT14 gene in locus 7p11.2 (SNP rs146350220; p=2.77 x10-4), and WDR73 gene in locus 15q25.2 (SNP rs72750868; p=4.47 x10-4).

Interpretation

We have identified susceptibility loci associated with each category of T2DM-related complications in the Emirati population. Given that only 16% of the markers from the Illumina's Infinium Exome chip passed quality control assessment, this demonstrates that multiple variants were, either, monomorphic in the Arab population or were not genotyped due to the use of a Euro-centric EWAS array that limits the possibility of including targeted ethnic-specific SNPs. Our results suggest the alarming possibility that lack of representation in reference panels could inhibit discovery of functionally important loci associated to T2DM complications. Further effort must be conducted to improve the representation of diverse populations in genotyping and sequencing studies.

GRHL2 and AP2a coordinate early surface ectoderm lineage commitment during development

Ectodermal dysplasias including skin abnormalities and cleft lip/palate result from improper surface ectoderm (SE) patterning. However, the connection between SE gene regulatory networks and disease remains poorly understood. Here, we dissect human SE differentiation with multiomics and establish GRHL2 as a key mediator of early SE commitment, which acts by skewing cell fate away from the neural lineage. GRHL2 and master SE regulator AP2a balance early cell fate output, with GRHL2 facilitating AP2a binding to SE loci. In turn, AP2a restricts GRHL2 DNA binding away from de novo chromatin contacts. Integration of these regulatory sites with ectodermal dysplasia-associated genomic variants annotated within the Biomedical Data Commons identifies 55 loci previously implicated in craniofacial disorders. These include ABCA4/ARHGAP29 and NOG regulatory regions where disease-linked variants directly affect GRHL2/AP2a binding and gene transcription. These studies elucidate the logic underlying SE commitment and deepen our understanding of human oligogenic disease pathogenesis.

Perspectives on the future of dysmorphology

The field of clinical genetics and genomics continues to evolve. In the past few decades, milestones like the initial sequencing of the human genome, dramatic changes in sequencing technologies, and the introduction of artificial intelligence, have upended the field and offered fascinating new insights. Though difficult to predict the precise paths the field will follow, rapid change may continue to be inevitable. Within genetics, the practice of dysmorphology, as defined by pioneering geneticist David W. Smith in the 1960s as "the study of, or general subject of abnormal development of tissue form" has also been affected by technological advances as well as more general trends in biomedicine. To address possibilities, potential, and perils regarding the future of dysmorphology, a group of clinical geneticists, representing different career stages, areas of focus, and geographic regions, have contributed to this piece by providing insights about how the practice of dysmorphology will develop over the next several decades.

MicroRNA-eQTLs in the developing human neocortex link miR-4707-3p expression to brain size

Expression quantitative trait loci (eQTL) data have proven important for linking non-coding loci to protein-coding genes. But eQTL studies rarely measure microRNAs (miRNAs), small non-coding RNAs known to play a role in human brain development and neurogenesis. Here, we performed small-RNA sequencing across 212 mid-gestation human neocortical tissue samples, measured 907 expressed miRNAs, discovering 111 of which were novel, and identified 85 local-miRNA-eQTLs. Colocalization of miRNA-eQTLs with GWAS summary statistics yielded one robust colocalization of miR-4707-3p expression with educational attainment and brain size phenotypes, where the miRNA expression increasing allele was associated with decreased brain size. Exogenous expression of miR-4707-3p in primary human neural progenitor cells decreased expression of predicted targets and increased cell proliferation, indicating miR-4707-3p modulates progenitor gene regulation and cell fate decisions. Integrating miRNA-eQTLs with existing GWAS yielded evidence of a miRNA that may influence human brain size and function via modulation of neocortical brain development.

Increasing diversity in connectomics with the Chinese Human Connectome Project

Cultural differences and biological diversity play important roles in shaping human brain structure and function. To date, most large-scale multimodal neuroimaging datasets have been obtained primarily from people living in Western countries, omitting the crucial contrast with populations living in other regions. The Chinese Human Connectome Project (CHCP) aims to address these resource and knowledge gaps by acquiring imaging, genetic and behavioral data from a large sample of participants living in an Eastern culture. The CHCP collected multimodal neuroimaging data from healthy Chinese adults using a protocol comparable to that of the Human Connectome Project. Comparisons between the CHCP and Human Connectome Project revealed both commonalities and distinctions in brain structure, function and connectivity. The corresponding large-scale brain parcellations were highly reproducible across the two datasets, with the language processing task showing the largest differences. The CHCP dataset is publicly available in an effort to facilitate transcultural and cross-ethnic brain-mind studies.

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

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

Heritability enrichment in context-specific regulatory networks improves phenotype-relevant tissue identification

Systems genetics holds the promise to decipher complex traits by interpreting their associated SNPs through gene regulatory networks derived from comprehensive multi-omics data of cell types, tissues, and organs. Here, we propose SpecVar to integrate paired chromatin accessibility and gene expression data into context-specific regulatory network atlas and regulatory categories, conduct heritability enrichment analysis with genome-wide association studies (GWAS) summary statistics, identify relevant tissues, and estimate relevance correlation to depict common genetic factors acting in the shared regulatory networks between traits. Our method improves power upon existing approaches by associating SNPs with context-specific regulatory elements to assess heritability enrichments and by explicitly prioritizing gene regulations underlying relevant tissues. Ablation studies, independent data validation, and comparison experiments with existing methods on GWAS of six phenotypes show that SpecVar can improve heritability enrichment, accurately detect relevant tissues, and reveal causal regulations. Furthermore, SpecVar correlates the relevance patterns for pairs of phenotypes and better reveals shared SNP-associated regulations of phenotypes than existing methods. Studying GWAS of 206 phenotypes in UK Biobank demonstrates that SpecVar leverages the context-specific regulatory network atlas to prioritize phenotypes' relevant tissues and shared heritability for biological and therapeutic insights. SpecVar provides a powerful way to interpret SNPs via context-specific regulatory networks and is available at https://github.com/AMSSwanglab/SpecVar, copy archived at swh:1:rev:cf27438d3f8245c34c357ec5f077528e6befe829.