Behaviour-correlated profiles of cerebellar-cerebral functional connectivity observed in independent neurodevelopmental disorder cohorts

The cerebellum, through its connectivity with the cerebral cortex, plays an integral role in regulating cognitive and affective processes, and its dysregulation can result in neurodevelopmental disorder (NDD)-related behavioural deficits. Identifying cerebellar-cerebral functional connectivity (FC) profiles in children with NDDs can provide insight into common connectivity profiles and their correlation to NDD-related behaviours. 479 participants from the Province of Ontario Neurodevelopmental Disorders (POND) network (typically developing = 93, Autism Spectrum Disorder = 172, Attention Deficit/Hyperactivity Disorder = 161, Obsessive-Compulsive Disorder = 53, mean age = 12.2) underwent resting-state functional magnetic resonance imaging and behaviour testing (Social Communication Questionnaire, Toronto Obsessive-Compulsive Scale, and Child Behaviour Checklist - Attentional Problems Subscale). FC components maximally correlated to behaviour were identified using canonical correlation analysis. Results were then validated by repeating the investigation in 556 participants from an independent NDD cohort provided from a separate consortium (Healthy Brain Network (HBN)). Replication of canonical components was quantified by correlating the feature vectors between the two cohorts. The two cerebellar-cerebral FC components that replicated to the greatest extent were correlated to, respectively, obsessive-compulsive behaviour (behaviour feature vectors, rPOND-HBN = -0.97; FC feature vectors, rPOND-HBN = -0.68) and social communication deficit contrasted against attention deficit behaviour (behaviour feature vectors, rPOND-HBN = -0.99; FC feature vectors, rPOND-HBN = -0.78). The statistically stable (|z| > 1.96) features of the FC feature vectors, measured via bootstrap re-sampling, predominantly comprised of correlations between cerebellar attentional and control network regions and cerebral attentional, default mode, and control network regions. In both cohorts, spectral clustering on FC loading values resulted in subject clusters mixed across diagnostic categories, but no cluster was significantly enriched for any given diagnosis as measured via chi-squared test (p > 0.05). Overall, two behaviour-correlated components of cerebellar-cerebral functional connectivity were observed in two independent cohorts. This suggests the existence of generalizable cerebellar network differences that span across NDD diagnostic boundaries.

Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Buffering of transcription rate by mRNA half-life is a conserved feature of Rett syndrome models

Transcriptional changes in Rett syndrome (RTT) are assumed to directly correlate with steady-state mRNA levels, but limited evidence in mice suggests that changes in transcription can be compensated by post-transcriptional regulation. We measure transcription rate and mRNA half-life changes in RTT patient neurons using RATEseq, and re-interpret nuclear and whole-cell RNAseq from Mecp2 mice. Genes are dysregulated by changing transcription rate or half-life and are buffered when both change. We utilized classifier models to predict the direction of transcription rate changes and find that combined frequencies of three dinucleotides are better predictors than CA and CG. MicroRNA and RNA-binding Protein (RBP) motifs are enriched in 3'UTRs of genes with half-life changes. Nuclear RBP motifs are enriched on buffered genes with increased transcription rate. We identify post-transcriptional mechanisms in humans and mice that alter half-life or buffer transcription rate changes when a transcriptional modulator gene is mutated in a neurodevelopmental disorder.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Contribution of rare variant associations to neurodegenerative disease presentation

Genetic factors contribute to neurodegenerative diseases, with high heritability estimates across diagnoses; however, a large portion of the genetic influence remains poorly understood. Many previous studies have attempted to fill the gaps by performing linkage analyses and association studies in individual disease cohorts, but have failed to consider the clinical and pathological overlap observed across neurodegenerative diseases and the potential for genetic overlap between the phenotypes. Here, we leveraged rare variant association analyses (RVAAs) to elucidate the genetic overlap among multiple neurodegenerative diagnoses, including Alzheimer's disease, amyotrophic lateral sclerosis, frontotemporal dementia (FTD), mild cognitive impairment, and Parkinson's disease (PD), as well as cerebrovascular disease, using the data generated with a custom-designed neurodegenerative disease gene panel in the Ontario Neurodegenerative Disease Research Initiative (ONDRI). As expected, only ~3% of ONDRI participants harboured a monogenic variant likely driving their disease presentation. Yet, when genes were binned based on previous disease associations, we observed an enrichment of putative loss of function variants in PD genes across all ONDRI cohorts. Further, individual gene-based RVAA identified significant enrichment of rare, nonsynonymous variants in PARK2 in the FTD cohort, and in NOTCH3 in the PD cohort. The results indicate that there may be greater heterogeneity in the genetic factors contributing to neurodegeneration than previously appreciated. Although the mechanisms by which these genes contribute to disease presentation must be further explored, we hypothesize they may be a result of rare variants of moderate phenotypic effect contributing to overlapping pathology and clinical features observed across neurodegenerative diagnoses.

Emotional face processing across neurodevelopmental disorders: a dynamic faces study in children with autism spectrum disorder, attention deficit hyperactivity disorder and obsessive-compulsive disorder

Autism spectrum disorder (ASD) is classically associated with poor face processing skills, yet evidence suggests that those with obsessive-compulsive disorder (OCD) and attention deficit hyperactivity disorder (ADHD) also have difficulties understanding emotions. We determined the neural underpinnings of dynamic emotional face processing across these three clinical paediatric groups, including developmental trajectories, compared with typically developing (TD) controls. We studied 279 children, 5-19 years of age but 57 were excluded due to excessive motion in fMRI, leaving 222: 87 ASD, 44 ADHD, 42 OCD and 49 TD. Groups were sex- and age-matched. Dynamic faces (happy, angry) and dynamic flowers were presented in 18 pseudo-randomized blocks while fMRI data were collected with a 3T MRI. Group-by-age interactions and group difference contrasts were analysed for the faces vs. flowers and between happy and angry faces. TD children demonstrated different activity patterns across the four contrasts; these patterns were more limited and distinct for the NDDs. Processing happy and angry faces compared to flowers yielded similar activation in occipital regions in the NDDs compared to TDs. Processing happy compared to angry faces showed an age by group interaction in the superior frontal gyrus, increasing with age for ASD and OCD, decreasing for TDs. Children with ASD, ADHD and OCD differentiated less between dynamic faces and dynamic flowers, with most of the effects seen in the occipital and temporal regions, suggesting that emotional difficulties shared in NDDs may be partly attributed to shared atypical visual information processing.

Structural neuroimaging biomarkers for obsessive-compulsive disorder in the ENIGMA-OCD consortium: medication matters

No diagnostic biomarkers are available for obsessive-compulsive disorder (OCD). Here, we aimed to identify magnetic resonance imaging (MRI) biomarkers for OCD, using 46 data sets with 2304 OCD patients and 2068 healthy controls from the ENIGMA consortium. We performed machine learning analysis of regional measures of cortical thickness, surface area and subcortical volume and tested classification performance using cross-validation. Classification performance for OCD vs. controls using the complete sample with different classifiers and cross-validation strategies was poor. When models were validated on data from other sites, model performance did not exceed chance-level. In contrast, fair classification performance was achieved when patients were grouped according to their medication status. These results indicate that medication use is associated with substantial differences in brain anatomy that are widely distributed, and indicate that clinical heterogeneity contributes to the poor performance of structural MRI as a disease marker.