•Competitive null models may yield false positives from co-expression.•Self-contained null models may yield false positives from bimodal correlations.•Test statistics interact differently with two ...types of null models.•Supplementary analyses with various configurations support the findings.
Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.
Background and Aims
Recently, we demonstrated that a distinct pattern of structural covariance networks (SCN) from magnetic resonance imaging (MRI)‐derived measurements of brain cortical thickness ...characterized young adults with alcohol use disorder (AUD) and predicted current and future problematic drinking in adolescents relative to controls. Here, we establish the robustness and value of SCN for identifying heavy alcohol users in three additional independent studies.
Design and Setting
Cross‐sectional and longitudinal studies using data from the Pediatric Imaging, Neurocognition and Genetics (PING) study (n = 400, age range = 14–22 years), the National Consortium on Alcohol and Neurodevelopment in Adolescence (NCANDA) (n = 272, age range = 17–22 years) and the Human Connectome Project (HCP) (n = 375, age range = 22–37 years).
Cases
Cases were defined based on heavy alcohol use patterns or former alcohol use disorder (AUD) diagnoses: 50, 68 and 61 cases were identified. Controls had none or low alcohol use or absence of AUD: 350, 204 and 314 controls were selected.
Measurements
Graph theory metrics of segregation and integration were used to summarize SCN.
Findings
Mirroring our prior findings, and across the three data sets, cases had a lower clustering coefficient area under the curve (AUC) = −0.029, P = 0.002, lower modularity (AUC = −0.14, P = 0.004), lower average shortest path length (AUC = −0.078, P = 0.017) and higher global efficiency (AUC = 0.007, P = 0.010). Local efficiency differences were marginal (AUC = −0.017, P = 0.052). That is, cases exhibited lower network segregation and higher integration, suggesting that adjacent nodes (i.e. brain regions) were less similar in thickness whereas spatially distant nodes were more similar.
Conclusion
Structural covariance network (SCN) differences in the brain appear to constitute an early marker of heavy alcohol use in three new data sets and, more generally, demonstrate the utility of SCN‐derived metrics to detect brain‐related psychopathology.
Nicotine and illicit stimulants are very addictive substances. Although associations between grey matter and dependence on stimulants have been frequently reported, white matter correlates have ...received less attention.
Eleven international sites ascribed to the ENIGMA-Addiction consortium contributed data from individuals with dependence on cocaine (n = 147), methamphetamine (n = 132) and nicotine (n = 189), as well as non-dependent controls (n = 333). We compared the fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD) and mean diffusivity (MD) of 20 bilateral tracts. Also, we compared the performance of various machine learning algorithms in deriving brain-based classifications on stimulant dependence.
The cocaine and methamphetamine groups had lower regional FA and higher RD in several association, commissural, and projection white matter tracts. The methamphetamine dependent group additionally showed lower regional AD. The nicotine group had lower FA and higher RD limited to the anterior limb of the internal capsule. The best performing machine learning algorithm was the support vector machine (SVM). The SVM successfully classified individuals with dependence on cocaine (AUC = 0.70, p < 0.001) and methamphetamine (AUC = 0.71, p < 0.001) relative to non-dependent controls. Classifications related to nicotine dependence proved modest (AUC = 0.62, p = 0.014).
Stimulant dependence was related to FA disturbances within tracts consistent with a role in addiction. The multivariate pattern of white matter differences proved sufficient to identify individuals with stimulant dependence, particularly for cocaine and methamphetamine.
•Evidence of white matter differences in stimulant dependence is inconsistent.•Most studies are often underpowered and limited to a few a priori selected tracts.•We provide robust evidence of white matter differences in stimulant dependence.•Machine learning methods can classify stimulant dependence using DTI-derived data.
While there is substantial evidence that cannabis use is associated with differences in human brain development, most of this evidence is correlational in nature. Bayesian causal network (BCN) ...modeling attempts to identify probable causal relationships in correlational data using conditional probabilities to estimate directional associations between a set of interrelated variables. In this study, we employed BCN modeling in 637 adolescents from the IMAGEN study who were cannabis naïve at age 14 to provide evidence that the accelerated prefrontal cortical thinning found previously in adolescent cannabis users by Albaugh et al. 1 is a result of cannabis use causally affecting neurodevelopment. BCNs incorporated data on cannabis use, prefrontal cortical thickness, and other factors related to both brain development and cannabis use, including demographics, psychopathology, childhood adversity, and other substance use. All BCN algorithms strongly suggested a directional relationship from adolescent cannabis use to accelerated cortical thinning. While BCN modeling alone does not prove a causal relationship, these results are consistent with a body of animal and human research suggesting that adolescent cannabis use adversely affects brain development.
The neurobiological bases of the association between development and psychopathology remain poorly understood. Here, we identify a shared spatial pattern of cortical thickness (CT) in normative ...development and several psychiatric and neurological disorders. Principal component analysis (PCA) was applied to CT of 68 regions in the Desikan-Killiany atlas derived from three large-scale datasets comprising a total of 41,075 neurotypical participants. PCA produced a spatially broad first principal component (PC1) that was reproducible across datasets. Then PC1 derived from healthy adult participants was compared to the pattern of CT differences associated with psychiatric and neurological disorders comprising a total of 14,886 cases and 20,962 controls from seven ENIGMA disease-related working groups, normative maturation and aging comprising a total of 17,697 scans from the ABCD Study® and the IMAGEN developmental study, and 17,075 participants from the ENIGMA Lifespan working group, as well as gene expression maps from the Allen Human Brain Atlas. Results revealed substantial spatial correspondences between PC1 and widespread lower CT observed in numerous psychiatric disorders. Moreover, the PC1 pattern was also correlated with the spatial pattern of normative maturation and aging. The transcriptional analysis identified a set of genes including KCNA2, KCNS1 and KCNS2 with expression patterns closely related to the spatial pattern of PC1. The gene category enrichment analysis indicated that the transcriptional correlations of PC1 were enriched to multiple gene ontology categories and were specifically over-represented starting at late childhood, coinciding with the onset of significant cortical maturation and emergence of psychopathology during the prepubertal-to-pubertal transition. Collectively, the present study reports a reproducible latent pattern of CT that captures interregional profiles of cortical changes in both normative brain maturation and a spectrum of psychiatric disorders. The pubertal timing of the expression of PC1-related genes implicates disrupted neurodevelopment in the pathogenesis of the spectrum of psychiatric diseases emerging during adolescence.
Brain asymmetry reflects left‐right hemispheric differentiation, which is a quantitative brain phenotype that develops with age and can vary with psychiatric diagnoses. Previous studies have shown ...that substance dependence is associated with altered brain structure and function. However, it is unknown whether structural brain asymmetries are different in individuals with substance dependence compared with nondependent participants. Here, a mega‐analysis was performed using a collection of 22 structural brain MRI datasets from the ENIGMA Addiction Working Group. Structural asymmetries of cortical and subcortical regions were compared between individuals who were dependent on alcohol, nicotine, cocaine, methamphetamine, or cannabis (n = 1,796) and nondependent participants (n = 996). Substance‐general and substance‐specific effects on structural asymmetry were examined using separate models. We found that substance dependence was significantly associated with differences in volume asymmetry of the nucleus accumbens (NAcc; less rightward; Cohen's d = 0.15). This effect was driven by differences from controls in individuals with alcohol dependence (less rightward; Cohen's d = 0.10) and nicotine dependence (less rightward; Cohen's d = 0.11). These findings suggest that disrupted structural asymmetry in the NAcc may be a characteristic of substance dependence.
A mega‐analysis with 22 datasets from the ENIGMA Addiction Working Group was performed. Structural asymmetries of cortical and subcortical regions were compared between individuals who were dependent on alcohol, nicotine, cocaine, methamphetamine, or cannabis (n = 1,796) and nondependent participants (n = 996). Less rightward asymmetry of the nucleus accumbens was observed in participants with substance dependence as compared with nondependent participants, suggesting that disrupted structural asymmetry in the nucleus accumbens may be a characteristic of substance dependence.
Synaptotagmin-1 is an essential regulator of synaptic vesicle exocytosis, and its encoding gene (SYT1) is a genome and transcriptome-wide association hit in cognitive performance, personality and ...cocaine use disorder (CUD) studies. Additionally, in candidate gene studies the specific variant rs2251214 has been associated with attention-deficit/hyperactivity disorder (ADHD), antisocial personality disorder and other externalizing phenotypes in adults with ADHD, as well as with response to methylphenidate (MPH) treatment. In this context, we sought to evaluate, in an independent sample, the association of this variant with CUD, a phenotype that shares common biological underpinnings with the previously associated traits. We tested the association between SYT1-rs2251214 and CUD susceptibility and severity (addiction severity index) in a sample composed by 315 patients addicted to smoked cocaine and 769 non-addicted volunteers. SYT1-rs2251214 was significantly associated with susceptibility to CUD, where the G allele presented increased risk for the disorder in the genetic models tested (P = 0.0021, OR = 1.44, allelic; P = 0.0012, OR = 1.48, additive; P = 0.0127, OR = 1.41, dominant). This is the same allele that was associated with increased risk for ADHD and other externalizing behaviors, as well as poor response to MPH treatment in previous studies. These findings suggest that the neurotransmitter exocytosis pathway might play a critical role in the liability for psychiatric disorders, especially externalizing behaviors and CUD.
•SYT1-rs2251214 was associated with cocaine use disorder.•The same allele associated with CUD was previously linked to externalizing behaviors.•The same allele associated with CUD was previously linked to stimulant response.
ADHD is associated with smaller subcortical brain volumes and cortical surface area, with greater effects observed in children than adults. It is also associated with dysregulation of the HPA axis. ...Considering the effects of the glucocorticoid receptor (NR3C1) in neurophysiology, we hypothesize that the blurred relationships between brain structures and ADHD in adults could be partly explained by
NR3C1
gene variation. Structural T1-weighted images were acquired on a 3 T scanner (
N
= 166). Large-scale genotyping was performed, and it was followed by quality control and pruning procedures, which resulted in 48 independent
NR3C1
gene variants analyzed. After a stringent Bonferroni correction, two SNPs (rs2398631 and rs72801070) moderated the association between ADHD and accumbens and amygdala volumes in adults. The significant SNPs that interacted with ADHD appear to have a role in gene expression regulation, and they are in linkage disequilibrium with
NR3C1
variants that present well-characterized physiological functions. The literature-reported associations of ADHD with accumbens and amygdala were only observed for specific
NR3C1
genotypes. Our findings reinforce the influence of the
NR3C1
gene on subcortical volumes and ADHD. They suggest a genetic modulation of the effects of a pivotal HPA axis component in the neuroanatomical features of ADHD.
There is evidence that dopamine receptors D2 (
DRD2
) and D4 (
DRD4
) polymorphisms may influence substance use disorders (SUD) susceptibility both individually and through their influence in the ...formation of DRD2–DRD4 heteromers. The dopaminergic role on the vulnerability to addiction appears to be influenced by sex. A cross-sectional study with 307 crack cocaine addicts and 770 controls was conducted. The influence of
DRD2
rs2283265 and
DRD4
48 bp VNTR in exon 3 variants, as well as their interaction on crack cocaine addiction susceptibility and severity were evaluated in women and men separately. An association between the
DRD2
T allele and crack cocaine addiction was found in women. In this same group, interaction analysis demonstrated that the presence of
DRD2
-T allele and concomitant absence of
DRD4
-7R allele were associated with risk for crack cocaine addiction. No influence of
DRD2
and
DRD4
variants was observed in men regarding addiction severity. This study reinforces the role of dopaminergic genes in externalizing behaviors, especially the influence of DRD2–DRD4 interaction on SUD. This is the fourth sample that independently associated the DRD2–DRD4 interaction with SUD itself or related disorders. In addition, our findings point out to a potential difference of dopaminergic neurotransmission across sex influencing addiction susceptibility.
Mixed findings exist in studies comparing brain responses to reward in adolescents and adults. Here we examined the trajectories of brain response, functional connectivity and task-modulated network ...properties during reward processing with a large-sample longitudinal design. Participants from the IMAGEN study performed a Monetary Incentive Delay task during fMRI at timepoint 1 (T1; n = 1304, mean age=14.44 years old) and timepoint 2 (T2; n = 1241, mean age=19.09 years). The Alcohol Use Disorders Identification Test (AUDIT) was administrated at both T1 and T2 to assess a participant’s alcohol use during the past year. Voxel-wise linear mixed effect models were used to compare whole brain response as well as functional connectivity of the ventral striatum (VS) during reward anticipation (large reward vs no-reward cue) between T1 and T2. In addition, task-modulated networks were constructed using generalized psychophysiological interaction analysis and summarized with graph theory metrics. To explore alcohol use in relation to development, participants with no/low alcohol use at T1 but increased alcohol use to hazardous use level at T2 (i.e., participants with AUDIT≤2 at T1 and ≥8 at T2) were compared against those with consistently low scores (i.e., participants with AUDIT≤2 at T1 and ≤7 at T2). Across the whole sample, lower brain response during reward anticipation was observed at T2 compared with T1 in bilateral caudate nucleus, VS, thalamus, midbrain, dorsal anterior cingulate as well as left precentral and postcentral gyrus. Conversely, greater response was observed bilaterally in the inferior and middle frontal gyrus and right precentral and postcentral gyrus at T2 (vs. T1). Increased functional connectivity with VS was found in frontal, temporal, parietal and occipital regions at T2. Graph theory metrics of the task-modulated network showed higher inter-regional connectivity and topological efficiency at T2. Interactive effects between time (T1 vs. T2) and alcohol use group (low vs. high) on the functional connectivity were observed between left middle temporal gyrus and right VS and the characteristic shortest path length of the task-modulated networks. Collectively, these results demonstrate the utility of the MID task as a probe of typical brain response and network properties during development and of differences in these features related to adolescent drinking, a reward-related behaviour associated with heightened risk for future negative health outcomes.
•Imaging data during reward anticipation at T1 (age 14) and T2 (age 19) was compared.•Brain response decreased in subcortical areas and increased in cortical areas at T2.•Functional connectivity (FC) with the ventral striatum increased at T2.•Topological efficiency of task-modulated network increased at T2.•The developmental pattern was altered in those who increased drinking most at T2.