Statistical and functional convergence of common and rare genetic influences on autism at chromosome 16p

Cortical gradient perturbation in attention deficit hyperactivity disorder correlates with neurotransmitter-, cell type-specific and chromosome- transcriptomic signatures

AIMS

This study aimed to illuminate the neuropathological landscape of attention deficit hyperactivity disorder (ADHD) by a multiscale macro-micro-molecular perspective from in vivo neuroimaging data.

METHODS

The "ADHD-200 initiative" repository provided multi-site high-quality resting-state functional connectivity (rsfc-) neuroimaging for ADHD children and matched typically developing (TD) cohort. Diffusion mapping embedding model to derive the functional connectome gradient detecting biologically plausible neural pattern was built, and the multivariate partial least square method to uncover the enrichment of neurotransmitomic, cellular and chromosomal gradient-transcriptional signatures of AHBA enrichment and meta-analytic decoding.

RESULTS

Compared to TD, ADHD children presented connectopic cortical gradient perturbations in almost all the cognition-involved brain macroscale networks (all pBH <0.001), but not in the brain global topology. As an intermediate phenotypic variant, such gradient perturbation was spatially enriched into distributions of GABAA/BZ and 5-HT2A receptors (all pBH <0.01) and co-varied with genetic transcriptional expressions (e.g. DYDC2, ATOH7, all pBH <0.01), associated with phenotypic variants in episodic memory and emotional regulations. Enrichment models demonstrated such gradient-transcriptional variants indicated the risk of both cell-specific and chromosome- dysfunctions, especially in enriched expression of oligodendrocyte precursors and endothelial cells (all pperm <0.05) as well enrichment into chromosome 18, 19 and X (pperm <0.05).

CONCLUSIONS

Our findings bridged brain macroscale neuropathological patterns to microscale/cellular biological architectures for ADHD children, demonstrating the neurobiologically pathological mechanism of ADHD into the genetic and molecular variants in GABA and 5-HT systems as well brain-derived enrichment of specific cellular/chromosomal expressions.

Polygenic scores for autism are associated with neurite density in adults and children from the general population

Genetic variants linked to autism are thought to change cognition and behaviour by altering the structure and function of the brain. Although a substantial body of literature has identified structural brain differences in autism, it is unknown whether autism-associated common genetic variants are linked to changes in cortical macro- and micro-structure. We investigated this using neuroimaging and genetic data from adults (UK Biobank, N = 31,748) and children (ABCD, N = 4,928). Using polygenic scores and genetic correlations we observe a robust negative association between common variants for autism and a magnetic resonance imaging derived phenotype for neurite density (intracellular volume fraction) in the general population. This result is consistent across both children and adults, in both the cortex and in white matter tracts, and confirmed using polygenic scores and genetic correlations. There were no sex differences in this association. Mendelian randomisation analyses provide no evidence for a causal relationship between autism and intracellular volume fraction, although this should be revisited using better powered instruments. Overall, this study provides evidence for shared common variant genetics between autism and cortical neurite density.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

Strategies for dissecting the complexity of neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) are associated with a wide range of clinical features, affecting multiple pathways involved in brain development and function. Recent advances in high-throughput sequencing have unveiled numerous genetic variants associated with NDDs, which further contribute to disease complexity and make it challenging to infer disease causation and underlying mechanisms. Herein, we review current strategies for dissecting the complexity of NDDs using model organisms, induced pluripotent stem cells, single-cell sequencing technologies, and massively parallel reporter assays. We further highlight single-cell CRISPR-based screening techniques that allow genomic investigation of cellular transcriptomes with high efficiency, accuracy, and throughput. Overall, we provide an integrated review of experimental approaches that can be applicable for investigating a broad range of complex disorders.

Dissecting 16p11.2 hemi-deletion to study sex-specific striatal phenotypes of neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) are polygenic in nature and copy number variants (CNVs) are ideal candidates to study the nature of this polygenic risk. The disruption of striatal circuits is considered a central mechanism in NDDs. The 16p11.2 hemi-deletion (16p11.2 del/+) is one of the most common CNVs associated with NDD, and 16p11.2 del/+ mice show sex-specific striatum-related behavioral phenotypes. However, the critical genes among the 27 genes in the 16p11.2 region that underlie these phenotypes remain unknown. Previously, we applied a novel strategy to identify candidate genes associated with the sex-specific phenotypes of 16p11.2 del/+ mice and highlighted three genes within the deleted region: thousand and one amino acid protein kinase 2 (Taok2), seizure-related 6 homolog-like 2 (Sez6l2), and major vault protein (Mvp). Using CRISPR/Cas9, we generated mice carrying null mutations in Taok2, Sez6l2, and Mvp (3 gene hemi-deletion (3g del/+)). Hemi-deletion of these 3 genes recapitulates sex-specific behavioral alterations in striatum-dependent behavioral tasks observed in 16p11.2 del/+ mice, specifically male-specific hyperactivity and impaired motivation for reward seeking. Moreover, RNAseq analysis revealed that 3g del/+ mice exhibit gene expression changes in the striatum similar to 16p11.2 del/+ mice exclusively in males. Subsequent analysis identified translation dysregulation and/or extracellular signal-regulated kinase signaling as plausible molecular mechanisms underlying male-specific, striatum-dependent behavioral alterations. Interestingly, ribosomal profiling supported the notion of translation dysregulation in both 3g del/+ and 16p11.2 del/+ male mice. However, mice carrying a 4-gene deletion (with an additional deletion of Mapk3) exhibited fewer phenotypic similarities with 16p11.2 del/+ mice. Together, the mutation of 3 genes within the 16p11.2 region phenocopies striatal sex-specific phenotypes of 16p11.2 del/+ mice. These results support the importance of a polygenic approach to study NDDs and underscore that the effects of the large genetic deletions result from complex interactions between multiple candidate genes.

Patient Brain Organoids Identify a Link between the 16p11.2 Copy Number Variant and the RBFOX1 Gene

Copy number variants (CNVs) that delete or duplicate 30 genes within the 16p11.2 genomic region give rise to a range of neurodevelopmental phenotypes with high penetrance in humans. Despite the identification of this small region, the mechanisms by which 16p11.2 CNVs lead to disease are unclear. Relevant models, such as human cortical organoids (hCOs), are needed to understand the human-specific mechanisms of neurodevelopmental disease. We generated hCOs from 17 patients and controls, profiling 167,958 cells with single-cell RNA-sequencing analysis, which revealed neuronal-specific differential expression of genes outside the 16p11.2 region that are related to cell-cell adhesion, neuronal projection growth, and neurodevelopmental disorders. Furthermore, 16p11.2 deletion syndrome organoids exhibited reduced mRNA and protein levels of RBFOX1, a gene that can also harbor CNVs linked to neurodevelopmental phenotypes. We found that the genes previously shown to be regulated by RBFOX1 are also perturbed in organoids from patients with the 16p11.2 deletion syndrome and thus identified a novel link between independent CNVs associated with neuronal development and autism. Overall, this work suggests convergent signaling, which indicates the possibility of a common therapeutic mechanism across multiple rare neuronal diseases.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.

Dissecting 16p11.2 hemi-deletion to study sex-specific striatal phenotypes of neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) are polygenic in nature and copy number variants (CNVs) are ideal candidates to study the nature of this polygenic risk. The disruption of striatal circuits is considered a central mechanism in NDDs. The 16p11.2 hemi-deletion (16p11.2 del) is one of the most common CNVs associated with NDD, and 16p11.2 del/+ mice show sex-specific striatum-related behavioral phenotypes. However, the critical genes among the 27 genes in the 16p11.2 region that underlie these phenotypes remain unknown. Previously, we applied a novel strategy to identify candidate genes associated with the sex-specific phenotypes of 16p11.2 del/+ mice and identified 3 genes of particular importance within the deleted region: thousand and one amino acid protein kinase 2 (Taok2), seizure-related 6 homolog-like 2 (Sez6l2), and major vault protein (Mvp). Using the CRISPR/Cas9 technique, we generated 3 gene hemi-deletion (3g del/+) mice carrying null mutations in Taok2, Sez6l2, and Mvp. We assessed striatum-dependent phenotypes of these 3g del/+ mice in behavioral, molecular, and imaging studies. Hemi-deletion of Taok2, Sez6l2, and Mvp induces sex-specific behavioral alterations in striatum-dependent behavioral tasks, specifically male-specific hyperactivity and impaired motivation for reward seeking, resembling behavioral phenotypes of 16p11.2 del/+ mice. Moreover, RNAseq analysis revealed that 3g del/+ mice exhibit gene expression changes in the striatum similar to 16p11.2 del/+ mice, but only in males. Pathway analysis identified ribosomal dysfunction and translation dysregulation as molecular mechanisms underlying male-specific, striatum-dependent behavioral alterations. Together, the mutation of 3 genes within the 16p11.2 region phenocopies striatal sex-specific phenotypes of 16p11.2 del/+ mice, unlike single gene mutation studies. These results support the importance of a polygenic approach to study NDDs and our novel strategy to identify genes of interest using gene expression patterns in brain regions, such as the striatum, which are impacted in these disorders.

Dissecting 16p11.2 hemi-deletion to study sex-specific striatal phenotypes of neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) are polygenic in nature and copy number variants (CNVs) are ideal candidates to study the nature of this polygenic risk. The disruption of striatal circuits is considered a central mechanism in NDDs. The 16p11.2 hemi-deletion (16p11.2 del) is one of the most common CNVs associated with NDD, and 16p11.2 del/+ mice show sex-specific striatum-related behavioral phenotypes. However, the critical genes among the 27 genes in the 16p11.2 region that underlie these phenotypes remain unknown. Previously, we applied a novel strategy to identify candidate genes associated with the sex-specific phenotypes of 16p11.2 del/+ mice and identified 3 genes of particular importance within the deleted region: thousand and one amino acid protein kinase 2 ( Taok2 ), seizure-related 6 homolog-like 2 ( Sez6l2 ), and major vault protein ( Mvp ). Using the CRISPR/Cas9 technique, we generated 3 gene hemi-deletion (3g del/+) mice carrying null mutations in Taok2, Sez6l2 , and Mvp . We assessed striatum-dependent phenotypes of these 3g del/+ mice in behavioral, molecular, and imaging studies. Hemi-deletion of Taok2, Sez6l2 , and Mvp induces sex-specific behavioral alterations in striatum-dependent behavioral tasks, specifically male-specific hyperactivity and impaired motivation for reward seeking, resembling behavioral phenotypes of 16p11.2 del/+ mice. Moreover, RNAseq analysis revealed that 3g del/+ mice exhibit gene expression changes in the striatum similar to 16p11.2 del/+ mice, but only in males. Pathway analysis identified ribosomal dysfunction and translation dysregulation as molecular mechanisms underlying male-specific, striatum-dependent behavioral alterations. Together, the mutation of 3 genes within the 16p11.2 region phenocopies striatal sex-specific phenotypes of 16p11.2 del/+ mice, unlike single gene mutation studies. These results support the importance of a polygenic approach to study NDDs and our novel strategy to identify genes of interest using gene expression patterns in brain regions, such as the striatum, which are impacted in these disorders.