Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network.

Genes disrupted in schizophrenia may be revealed by de novo mutations in affected persons from otherwise healthy families. Furthermore, during normal brain development, genes are expressed in patterns specific to developmental stage and neuroanatomical structure. We identified de novo mutations in persons with schizophrenia and then mapped the responsible genes onto transcriptome profiles of normal human brain tissues from age 13 weeks gestation to adulthood. In the dorsolateral and ventrolateral prefrontal cortex during fetal development, genes harboring damaging de novo mutations in schizophrenia formed a network significantly enriched for transcriptional coexpression and protein interaction. The 50 genes in the network function in neuronal migration, synaptic transmission, signaling, transcriptional regulation, and transport. These results suggest that disruptions of fetal prefrontal cortical neurogenesis are critical to the pathophysiology of schizophrenia. These results also support the feasibility of integrating genomic and transcriptome analyses to map critical neurodevelopmental processes in time and space in the brain.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth.

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Genomic loci influence patterns of structural covariance in the human brain.

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

Regression and alignment for functional data and network topology.

In the brain, functional connections form a network whose topological organization can be described by graph-theoretic network diagnostics. These include characterizations of the community structure, such as modularity and participation coefficient, which have been shown to change over the course of childhood and adolescence. To investigate if such changes in the functional network are associated with changes in cognitive performance during development, network studies often rely on an arbitrary choice of preprocessing parameters, in particular the proportional threshold of network edges. Because the choice of parameter can impact the value of the network diagnostic, and therefore downstream conclusions, we propose to circumvent that choice by conceptualizing the network diagnostic as a function of the parameter. As opposed to a single value, a network diagnostic curve describes the connectome topology at multiple scales-from the sparsest group of the strongest edges to the entire edge set. To relate these curves to executive function and other covariates, we use scalar-on-function regression, which is more flexible than previous functional data-based models used in network neuroscience. We then consider how systematic differences between networks can manifest in misalignment of diagnostic curves, and consequently propose a supervised curve alignment method that incorporates auxiliary information from other variables. Our algorithm performs both functional regression and alignment via an iterative, penalized, and nonlinear likelihood optimization. The illustrated method has the potential to improve the interpretability and generalizability of neuroscience studies where the goal is to study heterogeneity among a mixture of function- and scalar-valued measures.

ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion.

Arterial spin labeled (ASL) magnetic resonance imaging (MRI) is the primary method for noninvasively measuring regional brain perfusion in humans. We introduce ASLPrep, a suite of software pipelines that ensure the reproducible and generalizable processing of ASL MRI data.

Polygenic risk of social isolation behavior and its influence on psychopathology and personality.

Social isolation has been linked to a range of psychiatric issues, but the behavioral component that drives it is not well understood. Here, a genome-wide associations study (GWAS) was carried out to identify genetic variants that contribute specifically to social isolation behavior (SIB) in up to 449,609 participants from the UK Biobank. 17 loci were identified at genome-wide significance, contributing to a 4% SNP-based heritability estimate. Using the SIB GWAS, polygenic risk scores (PRS) were derived in ALSPAC, an independent, developmental cohort, and used to test for association with self-reported friendship scores, comprising items related to friendship quality and quantity, at age 12 and 18 to determine whether genetic predisposition manifests during childhood development. At age 18, friendship scores were associated with the SIB PRS, demonstrating that the genetic factors can predict related social traits in late adolescence. Linkage disequilibrium (LD) score correlation using the SIB GWAS demonstrated genetic correlations with autism spectrum disorder (ASD), schizophrenia, major depressive disorder (MDD), educational attainment, extraversion, and loneliness. However, no evidence of causality was found using a conservative Mendelian randomization approach between SIB and any of the traits in either direction. Genomic Structural Equation Modeling (SEM) revealed a common factor contributing to SIB, neuroticism, loneliness, MDD, and ASD, weakly correlated with a second common factor that contributes to psychiatric and psychotic traits. Our results show that SIB contributes a small heritable component, which is associated genetically with other social traits such as friendship as well as psychiatric disorders.

Psychosis superspectrum I: Nosology, etiology, and lifespan development.

This review describes the Hierarchical Taxonomy of Psychopathology (HiTOP) model of psychosis-related psychopathology, the psychosis superspectrum. The HiTOP psychosis superspectrum was developed to address shortcomings of traditional diagnoses for psychotic disorders and related conditions including low reliability, arbitrary boundaries between psychopathology and normality, high symptom co-occurrence, and heterogeneity within diagnostic categories. The psychosis superspectrum is a transdiagnostic dimensional model comprising two spectra-psychoticism and detachment-which are in turn broken down into fourteen narrow components, and two auxiliary domains-cognition and functional impairment. The structure of the spectra and their components are shown to parallel the genetic structure of psychosis and related traits. Psychoticism and detachment have distinct patterns of association with urbanicity, migrant and ethnic minority status, childhood adversity, and cannabis use. The superspectrum also provides a useful model for describing the emergence and course of psychosis, as components of the superspectrum are relatively stable over time. Changes in psychoticism predict the onset of psychosis-related psychopathology, whereas changes in detachment and cognition define later course. Implications of the superspectrum for genetic, socio-environmental, and longitudinal research are discussed. A companion review focuses on neurobiology, treatment response, and clinical utility of the superspectrum, and future research directions.

Neurocognitive profiles of 22q11.2 and 16p11.2 deletions and duplications.

Rare recurrent copy number variants (CNVs) at chromosomal loci 22q11.2 and 16p11.2 are genetic disorders with lifespan risk for neuropsychiatric disorders. Microdeletions and duplications are associated with neurocognitive deficits, yet few studies compared these groups using the same measures to address confounding measurement differences. We report a prospective international collaboration applying the same computerized neurocognitive assessment, the Penn Computerized Neurocognitive Battery (CNB), administered in a multi-site study on rare genomic disorders: 22q11.2 deletions (n = 492); 22q11.2 duplications (n = 106); 16p11.2 deletion (n = 117); and 16p11.2 duplications (n = 46). Domains examined include executive functions, episodic memory, complex cognition, social cognition, and psychomotor speed. Accuracy and speed for each domain were included as dependent measures in a mixed-model repeated measures analysis. Locus (22q11.2, 16p11.2) and Copy number (deletion/duplication) were grouping factors and Measure (accuracy, speed) and neurocognitive domain were repeated measures factors, with Sex and Site as covariates. We also examined correlation with IQ. We found a significant Locus × Copy number × Domain × Measure interaction (p = 0.0004). 22q11.2 deletions were associated with greater performance accuracy deficits than 22q11.2 duplications, while 16p11.2 duplications were associated with greater specific deficits than 16p11.2 deletions. Duplications at both loci were associated with reduced speed compared to deletions. Performance profiles differed among the groups with particularly poor memory performance of the 22q11.2 deletion group while the 16p11.2 duplication group had greatest deficits in complex cognition. Average accuracy on the CNB was moderately correlated with Full Scale IQ. Deletions and duplications of 22q11.2 and 16p11.2 have differential effects on accuracy and speed of neurocognition indicating locus specificity of performance profiles. These profile differences can help inform mechanistic substrates to heterogeneity in presentation and outcome, and can only be established in large-scale international consortia using the same neurocognitive assessment. Future studies could aim to link performance profiles to clinical features and brain function.

Sex differences in the functional topography of association networks in youth.

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

Immune status of the murine 22q11.2 deletion syndrome model.

Mice modeling the hemizygous deletion of chromosome 22q11.2 (22qMc) have been utilized to address various clinical phenotypes associated with the disease, including cardiac malformations, altered neural circuitry, and behavioral deficits. Yet, the status of the T cell compartment, an important clinical concern among 22q11.2 deletion syndrome (22qDS) patients, has not been addressed. While infancy and early childhood in 22qDS are associated with deficient T cell numbers and thymic hypoplasia, which can be severe in a small subset of patients, studies suggest normalization of the T cell counts by adulthood. We found that adult 22qMc do not exhibit thymic hypoplasia or altered thymic T cell development. Our findings that immune cell counts and inflammatory T cell activation are unaffected in 22qMc lend support to the hypothesis that human 22qDS immunodeficiencies are secondary to thymic hypoplasia, rather than intrinsic effects due to the deletion. Furthermore, the 22q11.2 deletion does not impact the differentiation capacity of T cells, nor their activity and response during inflammatory activation. Thus, 22qMc reflects the T cell compartment in adult 22qDS patients, and our findings suggest that 22qMc may serve as a novel model to address experimental and translational aspects of immunity in 22qDS.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Source-based morphometry reveals structural brain pattern abnormalities in 22q11.2 deletion syndrome.

22q11.2 deletion syndrome (22q11DS) is the most frequently occurring microdeletion in humans. It is associated with a significant impact on brain structure, including prominent reductions in gray matter volume (GMV), and neuropsychiatric manifestations, including cognitive impairment and psychosis. It is unclear whether GMV alterations in 22q11DS occur according to distinct structural patterns. Then, 783 participants (470 with 22q11DS: 51% females, mean age [SD] 18.2 [9.2]; and 313 typically developing [TD] controls: 46% females, mean age 18.0 [8.6]) from 13 datasets were included in the present study. We segmented structural T1-weighted brain MRI scans and extracted GMV images, which were then utilized in a novel source-based morphometry (SBM) pipeline (SS-Detect) to generate structural brain patterns (SBPs) that capture co-varying GMV. We investigated the impact of the 22q11.2 deletion, deletion size, intelligence quotient, and psychosis on the SBPs. Seventeen GMV-SBPs were derived, which provided spatial patterns of GMV covariance associated with a quantitative metric (i.e., loading score) for analysis. Patterns of topographically widespread differences in GMV covariance, including the cerebellum, discriminated individuals with 22q11DS from healthy controls. The spatial extents of the SBPs that revealed disparities between individuals with 22q11DS and controls were consistent with the findings of the univariate voxel-based morphometry analysis. Larger deletion size was associated with significantly lower GMV in frontal and occipital SBPs; however, history of psychosis did not show a strong relationship with these covariance patterns. 22q11DS is associated with distinct structural abnormalities captured by topographical GMV covariance patterns that include the cerebellum. Findings indicate that structural anomalies in 22q11DS manifest in a nonrandom manner and in distinct covarying anatomical patterns, rather than a diffuse global process. These SBP abnormalities converge with previously reported cortical surface area abnormalities, suggesting disturbances of early neurodevelopment as the most likely underlying mechanism.

Brain-age prediction: Systematic evaluation of site effects, and sample age range and size.

Structural neuroimaging data have been used to compute an estimate of the biological age of the brain (brain-age) which has been associated with other biologically and behaviorally meaningful measures of brain development and aging. The ongoing research interest in brain-age has highlighted the need for robust and publicly available brain-age models pre-trained on data from large samples of healthy individuals. To address this need we have previously released a developmental brain-age model. Here we expand this work to develop, empirically validate, and disseminate a pre-trained brain-age model to cover most of the human lifespan. To achieve this, we selected the best-performing model after systematically examining the impact of seven site harmonization strategies, age range, and sample size on brain-age prediction in a discovery sample of brain morphometric measures from 35,683 healthy individuals (age range: 5-90 years; 53.59% female). The pre-trained models were tested for cross-dataset generalizability in an independent sample comprising 2101 healthy individuals (age range: 8-80 years; 55.35% female) and for longitudinal consistency in a further sample comprising 377 healthy individuals (age range: 9-25 years; 49.87% female). This empirical examination yielded the following findings: (1) the accuracy of age prediction from morphometry data was higher when no site harmonization was applied; (2) dividing the discovery sample into two age-bins (5-40 and 40-90 years) provided a better balance between model accuracy and explained age variance than other alternatives; (3) model accuracy for brain-age prediction plateaued at a sample size exceeding 1600 participants. These findings have been incorporated into CentileBrain (https://centilebrain.org/#/brainAGE2), an open-science, web-based platform for individualized neuroimaging metrics.

Similarity-based multimodal regression.

To better understand complex human phenotypes, large-scale studies have increasingly collected multiple data modalities across domains such as imaging, mobile health, and physical activity. The properties of each data type often differ substantially and require either separate analyses or extensive processing to obtain comparable features for a combined analysis. Multimodal data fusion enables certain analyses on matrix-valued and vector-valued data, but it generally cannot integrate modalities of different dimensions and data structures. For a single data modality, multivariate distance matrix regression provides a distance-based framework for regression accommodating a wide range of data types. However, no distance-based method exists to handle multiple complementary types of data. We propose a novel distance-based regression model, which we refer to as Similarity-based Multimodal Regression (SiMMR), that enables simultaneous regression of multiple modalities through their distance profiles. We demonstrate through simulation, imaging studies, and longitudinal mobile health analyses that our proposed method can detect associations between clinical variables and multimodal data of differing properties and dimensionalities, even with modest sample sizes. We perform experiments to evaluate several different test statistics and provide recommendations for applying our method across a broad range of scenarios.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Domain-specific associations between psychopathology and neurocognitive functioning.

BACKGROUND: Neurocognitive dysfunction is a transdiagnostic finding in psychopathology, but relationships among cognitive domains and general and specific psychopathology dimensions remain unclear. This study aimed to examine associations between cognition and psychopathology dimensions in a large youth cohort. METHOD: The sample (N = 9350; age 8-21 years) was drawn from the Philadelphia Neurodevelopmental Cohort. Data from structured clinical interviews were modeled using bifactor confirmatory factor analysis (CFA), resulting in an overall psychopathology ('p') factor score and six orthogonal psychopathology dimensions: dysphoria/distress, obsessive-compulsive, behavioral/externalizing, attention-deficit/hyperactivity, phobias, and psychosis. Neurocognitive data were aggregated using correlated-traits CFA into five factors: executive functioning, memory, complex cognition, social cognition, and sensorimotor speed. We examined relationships among specific and general psychopathology dimensions and neurocognitive factors. RESULTS: The final model showed both overall and specific associations between cognitive functioning and psychopathology, with acceptable fit (CFI = 0.91; TLI = 0.90; RMSEA = 0.024; SRMR = 0.054). Overall psychopathology and most psychopathology dimensions were negatively associated with neurocognitive functioning (phobias [p < 0.0005], behavioral/externalizing [p < 0.0005], attention-deficit/hyperactivity [p < 0.0005], psychosis [p < 0.0005 to p < 0.05]), except for dysphoria/distress and obsessive-compulsive symptoms, which were positively associated with complex cognition (p < 0.05 and p < 0.01, respectively). CONCLUSION: By modeling a broad range of cognitive and psychopathology domains in a large, diverse sample of youth, we found aspects of neurocognitive functioning shared across clinical phenotypes, as well as domain-specific patterns. Findings support transdiagnostic examination of cognitive performance to parse variability in the link between neurocognitive functioning and clinical phenotypes.

Immune and oxidative stress biomarkers in pediatric psychosis and psychosis-risk: Meta-analyses and systematic review.

OBJECTIVE: While genetic and cohort studies suggest immune and reduction/oxidation (redox) alterations occur in psychosis, less is known about potential alterations in children and adolescents. METHODS: We conducted a systematic review to identify immune and redox biomarker studies in children and adolescents (mean age ≤ 18 years old) across the psychosis spectrum: from psychotic like experiences, which are common in children, to threshold psychotic disorders like schizophrenia. We conducted meta-analyses when at least three studies measured the same biomarker. RESULTS: The systematic review includes 38 pediatric psychosis studies. The meta-analyses found that youth with threshold psychotic disorders had higher neutrophil/lymphocyte ratio (Hedge's g = 0.40, 95 % CI 0.17 - 0.64), tumor necrosis factor (Hedge's g = 0.38, 95 % CI 0.06 - 0.69), C-reactive protein (Hedge's g = 0.38, 95 % CI 0.05 - 0.70), interleukin-6 (Hedge's g = 0.35; 95 % CI 0.11 - 0.64), and total white blood cell count (Hedge's g = 0.29, 95 % CI 0.12 - 0.46) compared to youth without psychosis. Other immune and oxidative stress meta-analytic findings were very heterogeneous. CONCLUSION: Results from several studies are consistent with the hypothesis that signals often classified as "proinflammatory" are elevated in threshold pediatric psychotic disorders. Data are less clear for immune markers in subthreshold psychosis and redox markers across the subthreshold and threshold psychosis spectrum. Immune and redox biomarker intervention studies are lacking, and research investigating interventions targeting the immune system in threshold pediatric psychosis is especially warranted.

Traumatic stress load and stressor reactivity score associated with accelerated gray matter maturation in youths indexed by normative models.

Understanding how traumatic stress affects typical brain development during adolescence is critical to elucidate underlying mechanisms related to both maladaptive functioning and resilience after traumatic exposures. The current study aimed to map deviations from normative ranges of brain gray matter for youths with traumatic exposures. For each cortical and subcortical gray matter region, normative percentiles of variations were established using structural MRI from typically developing youths without any traumatic exposure (n = 245; age range = 8-23) from the Philadelphia Neurodevelopmental Cohort (PNC). The remaining PNC participants with neuroimaging data (n = 1129) were classified as either within the normative range (5-95%), delayed (>95%) or accelerated (<5%) maturational ranges for each region using the normative model. An averaged quantile regression index was calculated across all regions. Mediation models revealed that high traumatic stress load was positively associated with poorer cognitive functioning and greater psychopathology, and these associations were mediated by accelerated gray matter maturation. Furthermore, higher stressor reactivity scores, which represent a less resilient response under traumatic stress, were positively correlated with greater acceleration of gray matter maturation (r = 0.224, 95% CI = [0.17, 0.28], p < 0.001), suggesting that more accelerated maturation was linked to greater stressor response regardless of traumatic stress load. We conclude that traumatic stress is a source of deviation from normative brain development associated with poorer cognitive functioning and more psychopathology in the long run.

Alprazolam modulates persistence energy during emotion processing in first-degree relatives of individuals with schizophrenia: a network control study.

Schizophrenia is marked by deficits in facial affect processing associated with abnormalities in GABAergic circuitry, deficits also found in first-degree relatives. Facial affect processing involves a distributed network of brain regions including limbic regions like amygdala and visual processing areas like fusiform cortex. Pharmacological modulation of GABAergic circuitry using benzodiazepines like alprazolam can be useful for studying this facial affect processing network and associated GABAergic abnormalities in schizophrenia. Here, we use pharmacological modulation and computational modeling to study the contribution of GABAergic abnormalities toward emotion processing deficits in schizophrenia. Specifically, we apply principles from network control theory to model persistence energy - the control energy required to maintain brain activation states - during emotion identification and recall tasks, with and without administration of alprazolam, in a sample of first-degree relatives and healthy controls. Here, persistence energy quantifies the magnitude of theoretical external inputs during the task. We find that alprazolam increases persistence energy in relatives but not in controls during threatening face processing, suggesting a compensatory mechanism given the relative absence of behavioral abnormalities in this sample of unaffected relatives. Further, we demonstrate that regions in the fusiform and occipital cortices are important for facilitating state transitions during facial affect processing. Finally, we uncover spatial relationships (i) between regional variation in differential control energy (alprazolam versus placebo) and (ii) both serotonin and dopamine neurotransmitter systems, indicating that alprazolam may exert its effects by altering neuromodulatory systems. Together, these findings provide a new perspective on the distributed emotion processing network and the effect of GABAergic modulation on this network, in addition to identifying an association between schizophrenia risk and abnormal GABAergic effects on persistence energy during threat processing.

Rare coding variants as risk modifiers of the 22q11.2 deletion implicate postnatal cortical development in syndromic schizophrenia.

22q11.2 deletion is one of the strongest known genetic risk factors for schizophrenia. Recent whole-genome sequencing of schizophrenia cases and controls with this deletion provided an unprecedented opportunity to identify risk modifying genetic variants and investigate their contribution to the pathogenesis of schizophrenia in 22q11.2 deletion syndrome. Here, we apply a novel analytic framework that integrates gene network and phenotype data to investigate the aggregate effects of rare coding variants and identified modifier genes in this etiologically homogenous cohort (223 schizophrenia cases and 233 controls of European descent). Our analyses revealed significant additive genetic components of rare nonsynonymous variants in 110 modifier genes (adjusted P = 9.4E-04) that overall accounted for 4.6% of the variance in schizophrenia status in this cohort, of which 4.0% was independent of the common polygenic risk for schizophrenia. The modifier genes affected by rare coding variants were enriched with genes involved in synaptic function and developmental disorders. Spatiotemporal transcriptomic analyses identified an enrichment of coexpression between modifier and 22q11.2 genes in cortical brain regions from late infancy to young adulthood. Corresponding gene coexpression modules are enriched with brain-specific protein-protein interactions of SLC25A1, COMT, and PI4KA in the 22q11.2 deletion region. Overall, our study highlights the contribution of rare coding variants to the SCZ risk. They not only complement common variants in disease genetics but also pinpoint brain regions and developmental stages critical to the etiology of syndromic schizophrenia.