Skip to main content

Mendelian randomization studies of brain MRI yield insights into the pathogenesis of neuropsychiatric disorders

Abstract

Background

Observational studies have identified various associations between neuroimaging alterations and neuropsychiatric disorders. However, whether such associations could truly reflect causal relations remains still unknown.

Results

Here, we leveraged genome-wide association studies (GWAS) summary statistics for (1) 11 psychiatric disorders (sample sizes varied from n = 9,725 to 1,331,010); (2) 110 diffusion tensor imaging (DTI) measurement (sample size n = 17,706); (3) 101 region-of-interest (ROI) volumes, and investigate the causal relationship between brain structures and neuropsychiatric disorders by two-sample Mendelian randomization. Among all DTI-Disorder combinations, we observed a significant causal association between the superior longitudinal fasciculus (SLF) and the risk of Anorexia nervosa (AN) (Odds Ratio [OR] = 0.62, 95 % confidence interval: 0.50 ~ 0.76, P = 6.4 × 10− 6). Similar significant associations were also observed between the body of the corpus callosum (fractional anisotropy) and Alzheimer’s disease (OR = 1.07, 95 % CI: 1.03 ~ 1.11, P = 4.1 × 10− 5). By combining all observations, we found that the overall p-value for DTI − Disorder associations was significantly elevated compared to the null distribution (Kolmogorov-Smirnov P = 0.009, inflation factor λ = 1.37), especially for DTI − Bipolar disorder (BP) (λ = 2.64) and DTI − AN (λ = 1.82). In contrast, for ROI-Disorder combinations, we only found a significant association between the brain region of pars triangularis and Schizophrenia (OR = 0.48, 95 % CI: 0.34 ~ 0.69, P = 5.9 × 10− 5) and no overall p-value elevation for ROI-Disorder analysis compared to the null expectation.

Conclusions

As a whole, we show that SLF degeneration may be a risk factor for AN, while DTI variations could be causally related to some neuropsychiatric disorders, such as BP and AN. Also, the white matter structure might have a larger impact on neuropsychiatric disorders than subregion volumes.

Background

Neuroimaging study is the most widely used procedure for studying brain disorders [1]. Many neuroimaging studies in the past quarter-century have revealed brain abnormalities in neuropsychiatric disorders [1, 2], which have served as the basis for biomarker discovery, clinical guidance, and investigations into the mechanisms of neuropsychiatric disorders [2, 3]. However, it is unclear whether such associations reflect disease causality [4]. One concern is that the spurious correlations [5] could emerge from indirect correlations with confounders such as medication, circadian and dietary changes, or false-positive events. Another issue is the direction of causality: the neurotoxicity hypothesis [6] suggests that psychiatric illnesses have toxic effects on the central nervous system, which leads to structural alterations following disease onset [7]. This theory gained support from several observations in which neuroimaging abnormalities exhibited dynamic progression during the neuropsychiatric disorder course [6, 7]. Under such theories, these associations should be utilized as clinical biomarkers rather than mechanism identifiers.

Several case-control studies with large sample sizes [8, 9] have found a significant correlation between neuroimaging alteration and neuropsychiatric disorder, yet still unable to distinguish the causality from the correlation. Longitudinal analyses have partially overcome the limitations of cross-sectional observational studies [10] for investigating the disease causality, such as detecting neuroimaging alterations in participants with the so-called high-risk status along with their progression [11, 12]. However, these studies are still limited by our current definition of high-risk cohorts. For example, the onset of neuropsychiatric disorders may be far earlier than a clinically recognizable high-risk status, such that the associated abnormalities may include changes that occurred after the primary pathology.

Despite the longitudinal analysis focusing on the high-risk status, an alternative method for addressing the challenge of causality is Mendelian Randomization (MR) [13]. MR has been used to derive the relationships between peripheral inflammatory markers and schizophrenia [14] and between physical activity and depression [15] with great successes. By selecting genetic locus that is strongly associated with exposure, so-called instruments, MR separates subjects into high and low lifetime exposure groups according to their genotypes on the instruments [16], then compares the prevalence of outcomes between the high and low exposure groups. Such grouping is considered unbiased as genotypes are randomly determined during meiosis. Furthermore, a two-sample Mendelian Randomization (MR) estimates instrument-exposure and instrument–outcome relationships in different cohorts to infer the exposure-outcome relationship, without the need for information on individuals [16]. The explosive growth of genome-wide association studies (GWAS) offers opportunities for applying MR to solve the debate of causality in different clinical medicine fields. Recently, two GWAS from UK-Biobank (UKBB) [17, 18] revealed the genetic basis of brain structural measurements, providing an opportunity to address clinical neuroimaging studies’ causality. In the present study, with the summary-statistics GWAS data of Magnetic resonance imaging (MRI) (N = 17,706 and 19,629) and twelve neuropsychiatric disorders (N = 9,725 to 1,331,010), we implemented a two-sample MR approach to detect the causal relationship between white matter (WM) structures (diffusion tensor imaging [DTI]), brain subregion volumes (region of interest [ROI]), and neuropsychiatric disorders. An overview of the study design is illustrated in Fig. 1.

Fig. 1
figure1

Flowchart of the study. Integrative analysis: p value distribution of MR result

Results

Causal associations between DTI and neuropsychiatric Disorder

After data harmonization, 1253/1320 (110*12) DTI-Disorder pairs and 973/1212 (101*12) ROI-Disorder pairs had at least one strong instrument (SNP with P < 5 × 10− 8 in DTI GWAS [17], Additional file 1) and were analyzed by MR (Additional file 2 and 3). We started by evaluating the DTI-Disorder association. In the inverse variance–weighted (IVW) analysis (outer layer of Fig. 2a), which required at least two instruments, one DTI-Disorder pair achieved study-wide significance (SW, P < 0.05/1320): superior longitudinal fasciculus Axial Diffusivity (SLF.AxD) – the risk of Anorexia nervosa (AN) (Fig. 2b; Table 1, odds ratio [OR] = 0.62, 95 % confidence interval [CI], 0.50 to 0.76, P = 6.4 × 10− 6). Another DTI-Disorder pair, body of corpus callosum Fractional Anisotropy (BCC.FA) − Alzheimer disease (AD), was close to reaching SW significance (Fig. 1c and Table 1, OR = 1.07, 95 % CI, 1.03 to 1.11, P = 4.1 × 10− 5). Both associations had relatively consistent results across instruments (Fig. 2b, c), in accordance with the fact that no outlier was detected (Additional file 2). Five other DTI-Disorder pairs also reached single-disease significance (P < 0.05/110) (Table 1). No significant result was found for the Wald ratio (second layer of Fig. 1a) applied to DTI-Disorder pairs with only one instrument.

Fig. 2
figure2

MR results for DTI-Disease associations. a: Each radius represents one DTI measures (110 in total: 22 white matter tracts × 5 DTI parameters). From outer to inner layer: -log10(p) for MR-IVW; -log10(p) for MR-Wald Ratio; -log10(p) for intercepts (int.) of MR-Egger regression; heritability of each DTI measures. The dotted grey line indicated the nominal p threshold (0.05); solid grey line indicated study-wide significance p threshold (0.05/1253). CST: Corticospinal tract. b&c: Forest plots showing MR effect of Superior longitudinal fasciculus (axial diusivities) on Anorexia nervosa (SLF.AxD-AN) and body of corpus callosum (fractional anisotropy) on Alzheimer Disease (BDD.FA-AD). Each line showed the single SNP MR effects (95 % confidence interval) estimated by Wald Ratio, and the last line showed the meta-analysis results calculated by IVW. Vertical dashed line indicated the Egger estimation of MR effect

Table 1 All DTI-Disorder pairs reaching Single Disease significant threshold (p < 0.05 after correction)

A basic assumption of MR is that MR genetic instruments should only impact the outcome via the exposure and not any other pathway (horizontal pleiotropy) [16]. The existence of heterogeneity, which is introduced by outlier instruments, can also bias the MR estimation [19]. By applying various sensitivity tests, we confirmed that our results were not impacted by horizontal pleiotropy (Egger intercept P > 0.05; the third layer of Fig. 1a and Additional file 2) or heterogeneity (modified Cochran’s Q test P > 0.05; Table 1), or bidirectional effects (reverse MR P > 0.05; Additional file 2). They also showed consistent trends across different MR methods, which are robust against pleiotropy and measurement error (see Method for detail) (except for Fornix.FA − bipolar disorder [BP]; Table 1). They also showed little directional pleiotropy, as indicated by funnel plots (Additional file 4). Neither MR-PRESSO nor leave-one-out tests found the impact of outliers on these DTI-Disorder pairs (Additional files 2 and 4). Taken together, these results confirmed the casual relation between SLF.AxD-AN, and suggested potential relations of other five DTI-Disorder pairs.

Overall contribution of DTI on neuropsychiatric disorders

Despite separate DTI-Disorder pairs that reached the significance threshold, we were also interested in whether DTI as a whole made a causal contribution to the diseases. To answer this question, we pooled the IVW P value for all DTI-Disorder pairs and compared them to the null uniform distribution (Fig. 3a and Additional file 2). The distribution of IVW P was significantly inflated (KS P = 0.009, λ = 1.37); this result persisted after removing heterogeneous DTI-Disorder pairs (those with Cochran’s P < 0.05) (KS P = 0.002, λ = 1.48) or outlier instruments (SNPs with MR-PRESSO P < 0.05) (KS P = 8.0 × 10− 5, λ = 1.49). When we analyzed each disease separately (“Original” in Fig. 3b), BP (Fig. 3c) and AN showed the most significant inflation (Table 2). These results were also relatively stable against heterogeneity and outlier removal (Fig. 3b and Additional file 4). The permutation test confirmed that the results for BP (permutation P value [pp] for KS test: 0.012; pp for λ: 0.002) and AN (pp for λ: 0.012) were not due to the bias inherent in the data or method. Additionally, heterogeneous DTI-Disorder pairs generally had non-significant MR results and did not contribute to inflation (Fig. 3d and Additional file 2). Thus, we concluded that DTI polymorphisms made an overall causal contribution to neuropsychiatric disorders, especially BP and AN.

Fig. 3
figure3

General contribution of DTI polymorphism to neuropsychiatric disorders. a: Quantile-Quantile (QQ) plot showing the distribution of all MR p values for the DTI-Disease (DD) association. b: Disease-specific inflation factor (λ). Solid points showed λ with no adjustment and were corresponded to λ in Table 2. Triangular points showed λ after removal of all DTI-Disorder pairs with significant heterogeneity. Cross points showed λ of p values which were calculated after removing all outlier SNPs (detected by MR-PRESSO) for each DTI-Disorder pair. Vertical error bars indicated a 95 % confidence interval. c: QQ plot for DTI-BP associations. d: Rank-Rank overlaps between MR effect and heterogeneity. The color of each grid corresponded to the proportion of each MR effect rank (each row summed up to 1). The number in each grid showed the exact number of DTI-Disorder pairs. NS: non-significant. Nominal: p < 0.05. SD: Single Disease level significance, p < 0.05/110. SW: Study-Wide significance, p < 0.05/1253

Table 2 Disease-specific p value distribution

Causal associations between brain volume and neuropsychiatric Disorders

Similar to the DTI-Disorder analysis, we also assessed the ROI-Disorder association. 973/1212 (101*12) RD pairs had at least one strong instrument (SNP with P < 5 × 10− 8) and were analyzed by MR. In the IVW analysis, we found no SW significant results (outer layer of Fig. 3a and Additional file 3). The Wald ratio revealed a marginal SW result for pars triangularis (PT)-SCZ (OR = 0.48, P = 5.93 × 10− 5) (second layer of Fig. 4a), which was driven by a single SNP (rs2279829). The only trait associated with rs2279829 in PhenoScanner [20] was daytime dozing or sleeping (P = 2.02 × 10− 6). The SMR-HEIDI test detected a significant MR effect in the same direction (OR = 0.48, P = 7.57 × 10− 4) with no evidence of colocalization pleiotropy (HEIDI P = 0.84). We estimated that per 1-SD increment in normalized PT volume, the risk of SCZ decreased by 52 % (OR = 0.48). In conclusion, these results suggested a potential causal relationship between PT and SCZ. However, the validation of this relationship requires more work on potential pleiotropy.

Fig. 4
figure4

MR results for ROI-Disease associations. a: Circos plot for ROI-Disease (RD) relations, similar to Fig. 1a. Arrow indicated result for pars triangularis-schizophrenia association, which was described in detail in the main text. b: QQ plot for all IVW p values of RD pairs, similar to Fig. 2a. c: Disease-specific inflation factor (λ), similar to Fig. 2b

The overall distribution of the IVW p-value did not significantly differ from the null distribution (Fig. 4b), even after removing heterogeneity (KS P = 0.12) and outliers (KS P = 0.07). However, IVW p-value for SCZ showed significant inflation (λ = 1.90, KS P = 0.001), and the result was significantly impacted by heterogeneity (Fig. 4c). Additionally, 10 % (6/58) of the nominally significant RD pairs showed heterogeneity, and 3 contained outlier SNPs (Additional files 3 and 4). In conclusion, there was no evidence that ROI polymorphisms had a universal contribution to neuropsychiatric disorders.

Discussion

The issue of causality has long beleaguered clinical neuroimaging studies [4]. Confirmation of a causal change can provide insights into disease mechanisms at the circuit and region levels and may reveal useful biomarkers for predicting prognosis. Many current neuroimaging studies cannot directly draw this conclusion [1, 4], which partly limits the translation of their findings to clinical practice. In this study, we conducted hypothesis-free, data-driven MR analyses to assess the causal relationship between neuroimaging polymorphisms and neuropsychiatric disorders in an unbiased manner.

Our results showed that, in general, WM connectivity was more closely associated with the risk of neuropsychiatric disorders than gray matter volume (GMV) (Fig. 3a and 4b). Among neuropsychiatric disorders, BP showed the most significant association with genetically determined connectivity polymorphisms (Fig. 3b,c). These results support the dysconnectivity theory [21] of psychiatry, positing that major psychiatric disorders such as SCZ and BP have common WM abnormalities in their pathology [22]. According to this hypothesis, anatomic and neurodevelopmental changes that arise from neurotoxicity [6] are a consequence rather than a cause of the illness. Indirect evidence from a functional study of neuromodulation and myelination [21] and a case-control study of the high-risk state [12] supports the dysconnectivity theory. Research interest has now shifted from the region-of-origin to a connectome’s concept [22], in which connections between brain regions rather than the regions themselves cause BP and other mental illnesses. Our finding that DTI is more closely associated with the onset of neuropsychiatric disorders than ROI provides supportive evidence for this paradigm shift.

Our study’s top MR result was a novel risk factor for AN—namely, decreased SLF.AxD. A few studies have reported a decreased SLF integrity in AN patients [23, 24], but the results were inconsistent with the marginal effect size. As neuroimaging findings in AN patients are influenced by dietary and metabolic alterations [25], a causal change may not be manifested as a visible signal. Nonetheless, the SLF may contribute to body image distortion in the pathology of AN [24], probably through its connection to areas responsible for body image perception (prefrontal and parietal networks) and self-perception (inferior parietal lobe) [26]. Although abnormalities in both GMV [26] and WM connectivity [23, 24] have been observed, our results suggest that the latter is a primary cause. In contrast, the former is a consequence of neuromodulator mechanisms such as activity-dependent pruning [21].

The roles of BCC in AD and PT in SCZ—2 marginally significant results from our MR analysis—have received more attention in the literature than SLF in AN. However, both MR results should be interpreted with caution. BCC atrophy is widely observed in AD patients even at an early stage and reflects Wallerian degeneration and myelin breakdown [27]. However, our MR analysis revealed a reverse association: the FA of BCC was positively associated with AD risk (β = 0.07). One possible explanation for this discrepancy is that an enlarged BCC in early life is a risk factor for AD development at an older age, causing BCC atrophy after disease onset. Confirmation of such a complex theory requires more robust evidence from large-scale longitudinal studies. As for the PT-SCZ association, although it was validated by several additional analyses such as SMR-HEIDI and reverse MR, a single SNP-driven MR result is by nature suspect due to the unexplored pleiotropy [16, 28]. Because the volume reduction of PT has been demonstrated in high-risk psychosis and first-episode SCZ patients [29], we suggest that the inferred causality between PT and SCZ is plausible.

There were some limitations to this study. Firstly, classic MR methods largely depend on high heritability and strong exposure instruments [16, 28]. However, both heritability and number of instruments [17, 18] vary across the tested neuroimaging parameters, such that the power of MR is inconsistent across all DTI-Disorder and RD pairs. In fact, 67 DTI-Disorder and 239 RD pairs were discarded at the beginning of our analysis due to the absence of instruments. Even if there were causal links among them, they would not have been detected in our study. Thus, negative results for DD/RD pairs with limited instruments are not as convincing as positive results for those with adequate instruments. For the positive results, it should also be noticed that only the genetically regulated proportion of polymorphisms are associated with the disorders. Secondly, the original GWAS sample sizes may have impacted the MR results since estimation accuracy (i.e., standard error for effect size) for the instrument–outcome relationship is directly linked to the confidence interval of the MR effect estimate [16, 28]. Since GWAS OCD and TS recruited fewer than 10,000 cases, their MR results were underpowered. In fact, several DD/RD pairs for OCD and TS had a large MR effect, but their wide CI range resulted in non-significant P values. Future GWAS with a larger sample size, both for neuroimaging polymorphism and neuropsychiatric disorders, will provide a better chance to improve our understanding in this field.

Conclusions

In conclusion, our analysis results demonstrate that, in general, WM structures make a more significant contribution to the etiology of neuropsychiatric disorders—especially BP and AN—than brain subregion volumes. SLF.AxD was causally related to AN; marginally significant relationships were also found between BCC.FA and AD and between PT and SCZ.

Methods

We obtained publicly available GWAS summary statistics for DTI, ROI, and neuropsychiatric without collecting any individual information. Ethics approval was obtained in each of the original studies; therefore, no further ethics approval was needed for the current study.

Data collection and preprocessing

The genetic instruments for DTI measurements have been previously described [17]. Briefly, the ENIGMA-DTI pipeline [30] was used to analyze UKBB diffusion MRI data for 17,706 European participants and generate 110 DTI parameters—namely, fractional anisotropy (FA), axial diffusivity (AxD), mean diffusivity, mode of anisotropy, and radial diffusivity, of 21 WM tracts as well as their mean values. The genetic instrument for ROI volumes was also obtained from UKBB GWAS [18], which included 19,629 European participants and used the standard OASIS-30 Atropos template for registration and Mindboggle-101 atlas for labeling [31].

We collected GWAS summary statistics from the following neuropsychiatric studies on European cohorts: (1) Alzheimer disease (AD) [32]; (2) Attention-deficit/ hyperactivity disorder (ADHD) [33]; (3) Anorexia Nervosa (AN) [34]; (4) anxiety disorders (Categorical phenotype) [35]; (5) autism spectrum disorder (ASD) [36]; 6)bipolar disorder (BP) [37]; 7) insomnia [38] 8) major depression disorders (MD) [39]; 9) obsessive-compulsory disorder (OCD) [40]; 10) posttraumatic stress disorder (PTSD) [41]; 11) schizophrenia (SCZ) [42]; 12) Tourette disorder (TD) [43] (Fig. 1). For GWAS from Psychiatric Genetic Consortium (PGC), there were few samples from UKB. For AD GWAS, the UKBB participants were not included in the case-control analysis. For other GWAS, we could not quantify the extent of sample overlap due to the lack of individual information. Since all GWAS used in the current study were conducted in European ancestry, we did not further adjust for the impact of population stratification.

For each DTI and ROI measurement, we retained single nucleotide polymorphisms (SNPs) with P < 5 × 10− 8 as strong instruments for MR; measurements without a strong instrument were discarded. We removed SNPs with linkage disequilibrium (LD) r2 ≥ 0.001 for each measurement using reference LD data from the 1000 Genomes Project [44]. Data harmonization was applied independently for each DTI-disease (DD) and ROI-disease (RD) pair with the TwoSampleMR R package [45]. Since many of the GWAS summary statistics analyzed in this study did not provide allele frequency information, we did not exclude SNPs based on ambiguous strand error. For all binary phenotypes, we log-transformed the odds ratio to generate the β value.

Power calculation

We calculated the variance in phenotype explained by each instrument by

$$ < mathdollar>{R}^2=\frac{2\ast EAF\ast \left(1- EAF\right)\ast {\beta}^2}{2\ast EAF\ast \left(1- EAF\right)\ast {\beta}^2+2\ast EAF\ast \left(1- EAF\right)\ast N\ast se{\left(\beta \right)}^2} $$

Where EAF was the effect allele frequency, β was the effect size, N was the sample size, and se(β) was the standard error of effect size. The F statistic was then denoted as

$$F=\frac{{R}^{2}*(N-2)}{1-{R}^{2}}$$

R2 and F were used to evaluate power for each instrument. For each RD and DTI-Disorder pair, we calculated the overall MR power using mRnd tool [46], assuming OR = 1.3 and type I error = 0.05. The assumption of OR was based on the actual MR effects passing the significance threshold. We took this assumption because there is limited observational estimation of OR that is currently available. MR power and the number of valid instruments for each pair are recorded in Additional files 2 and 3. Details for all instruments are shown in Additional file 1.

Calculation of MR effects

For DD/RD pairs with at least two instruments, we performed a meta-analysis of each instrument’s MR effect using the inverse variance–weighted (IVW) method. The results were considered preliminary results and were used for downstream analyses. For the top IVW findings, we additionally applied weighted mode [47], weighted median [13], and MR-Egger regression [48] approaches—which are relatively robust against horizontal pleiotropy [15]—to further confirm the validity of the MR effect. For DD/RD pairs with only one instrument, estimates based on the Wald ratio were considered preliminary results. For the top Wald ratio finding, we used summary data-based MR (SMR) [49] to confirm the existence of the MR effect and heterogeneity in dependent instruments (HEIDI) [49] to rule out the probability that the MR effect was driven by colocalization of the instrument with the effective locus.

Since both IVW and Wald ratio results were taken into account, the p-value was adjusted by the Bonferroni method by numbers of all DTI-Disorder pairs (1320) or ROI-Disorder pairs (1212).

Sensitivity analysis

The intercept of the Egger regression was used as an indicator of potential horizontal pleiotropy, while modified Cochran’s Q for IVW and Rucker’s Q for Egger regression [50] were used as an indicator of heterogeneity. We used MR pleiotropy residual sum and outlier (MR-PRESSO) [19] with the number of permutations = 2500 to detect potential outlier SNPs for each DD/RD pair and generate an overall p-value for heterogeneity; those with outlier(s) were reanalyzed by IVW after removing the outlier(s). We also applied leave-one-out tests for all top findings to further evaluate the effects of unknown outliers.

For all MR results, we tested the reverse MR effect (i.e., neuropsychiatric disorders as exposure and DTI/ROI as outcome) by selecting SNPs with P < 5 × 10− 8 for each neuropsychiatric disorder as instruments. The absence of a reverse MR effect (IVW P > 0.05) was considered as evidence for the validity of directionality.

Analysis of the general causal contribution

To assess the general contribution of DTI (ROI) polymorphisms to neuropsychiatric disorder, we pooled the IVW P values for all DTI-Disorder (RD) pairs and compared them to the null uniform distribution using quantile-quantile plots. A positive bias (inflation) from uniform distribution was considered as evidence for a general contribution. The significance of inflation was evaluated with the Kolmogorov–Smirnov (KS) test, while the extent of inflation was assessed with the inflation factor λ, which was calculated by chi-square regression using the GenABEL R package [51]. Since the inflation factor might be overestimated due to the small number of P values, we shuffled SNP labels for BP and AN 1,000 times to carry out a permutation test. Permutation P < 0.05 was considered evidence for significant inflation. These tests were also separately applied to each disease and repeated after removing heterogeneous results (those with Cochran’s P < 0.05) or outlier SNPs (those with MR-PRESSO P < 0.05).

Availability of data and materials

Disease GWAS summary data were downloaded from https://www.med.unc.edu/pgc/download-results/. GWAS summary of neuroimaging data were downloaded from https://github.com/BIG-S2/GWAS.

Abbreviations

2SMR:

2-sample Mendelian randomization.

AD:

Alzheimer Disease

AN:

Anorexia Nervosa

BP:

Bipolar disorder

DTI:

Diffusion tensor imaging

ROI:

Region-of-interest

References

  1. 1.

    Etkin A. A reckoning and research agenda for neuroimaging in psychiatry. Am J Psychiatry. 2019;176:507–11.

    PubMed  Article  Google Scholar 

  2. 2.

    Lui S, Zhou XJ, Sweeney JA, Gong Q. Psychoradiology: The frontier of neuroimaging in psychiatry. Radiology. 2016;281:357–72.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Aydin O, Unal Aydin P, Arslan A. Development of Neuroimaging-Based Biomarkers in Psychiatry. Adv Exp Med Biol. 2019;1192:159–95.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Etkin A. Addressing the causality gap in human psychiatric neuroscience. JAMA Psychiatry. 2018;75:3–4.

    PubMed  Article  Google Scholar 

  5. 5.

    Haig BD. What Is a Spurious Correlation? Underst Stat. 2003;2:125–32.

    Article  Google Scholar 

  6. 6.

    Weinberger DR, McClure RK. Neurotoxicity, neuroplasticity, and magnetic resonance imaging morphometry: What is happening in the schizophrenic brain? Arch Gen Psychiatry. 2002;59:553–8.

    PubMed  Article  Google Scholar 

  7. 7.

    DeLisi LE. Defining the course of brain structural change and plasticity in schizophrenia. Psychiatry Res. 1999;92:1–9.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Hoogman M, Bralten J, Hibar DP, Mennes M, Zwiers MP, Schweren LSJ, et al. Subcortical brain volume differences in participants with attention deficit hyperactivity disorder in children and adults: a cross-sectional mega-analysis. The Lancet Psychiatry. 2017;4:310–9.

    PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    van Erp TGM, Walton E, Hibar DP, Schmaal L, Jiang W, Glahn DC, et al. Cortical Brain Abnormalities in 4474 Individuals With Schizophrenia and 5098 Control Subjects via the Enhancing Neuro Imaging Genetics Through Meta Analysis (ENIGMA) Consortium. Biol Psychiatry. 2018;84:644–54.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Pantelis C, Velakoulis D, McGorry PD, Wood SJ, Suckling J, Phillips LJ, et al. Neuroanatomical abnormalities before and after onset of psychosis: a cross-sectional and longitudinal MRI comparison. Lancet. 2003;361:281–8.

    PubMed  Article  PubMed Central  Google Scholar 

  11. 11.

    Tang Y, Pasternak O, Kubicki M, Rathi Y, Zhang T, Wang J, et al. Altered cellular white matter but not extracellular free water on diffusion MRI in individuals at clinical high risk for psychosis. Am J Psychiatry. 2019;176:820–8.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Roberts G, Perry A, Lord A, Frankland A, Leung V, Holmes-Preston E, et al. Structural dysconnectivity of key cognitive and emotional hubs in young people at high genetic risk for bipolar disorder. Mol Psychiatry. 2018;23:413–21.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  13. 13.

    Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016;40:304–14.

    PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Hartwig FP, Borges MC, Horta BL, Bowden J, Davey Smith G. Inflammatory Biomarkers and Risk of Schizophrenia. JAMA Psychiatry. 2017;74:1226.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Choi KW, Chen CY, Stein MB, Klimentidis YC, Wang MJ, Koenen KC, et al. Assessment of Bidirectional Relationships between Physical Activity and Depression among Adults: A 2-Sample Mendelian Randomization Study. JAMA Psychiatry. 2019;76:399–408.

    PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum Mol Genet. 2018;27:R195–208.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Zhao B, Zhang J, Ibrahim JG, Luo T, Santelli RC, Li Y, et al. Large-scale GWAS reveals genetic architecture of brain white matter microstructure and genetic overlap with cognitive and mental health traits (n = 17,706). Mol Psychiatry. 2019. https://doi.org/10.1038/s41380-019-0569-z.

  18. 18.

    Zhao B, Luo T, Li T, Li Y, Zhang J, Shan Y, et al. Genome-wide association analysis of 19,629 individuals identifies variants influencing regional brain volumes and refines their genetic co-architecture with cognitive and mental health traits. Nat Genet. 2019;51:1637–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50:693–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Staley JR, Blackshaw J, Kamat MA, Ellis S, Surendran P, Sun BB, et al. PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics. 2016;32:3207–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Friston K, Brown HR, Siemerkus J, Stephan KE. The dysconnection hypothesis (2016). Schizophr Res. 2016;176:83–94.

    PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Perry A, Roberts G, Mitchell PB, Breakspear M. Connectomics of bipolar disorder: a critical review, and evidence for dynamic instabilities within interoceptive networks. Mol Psychiatry. 2019;24:1296–318.

    PubMed  Article  Google Scholar 

  23. 23.

    Phillipou A, Carruthers SP, Di Biase MA, Zalesky A, Abel LA, Castle DJ, et al. White matter microstructure in anorexia nervosa. Hum Brain Mapp. 2018;39:4385–92.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Barona M, Brown M, Clark C, Frangou S, White T, Micali N. White matter alterations in anorexia nervosa: Evidence from a voxel-based meta-analysis. Neurosci Biobehav Rev. 2019;100:285–95.

    PubMed  Article  Google Scholar 

  25. 25.

    Treasure J, Zipfel S, Micali N, Wade T, Stice E, Claudino A, et al. Anorexia nervosa. Nat Rev Dis Prim. 2015;1:1–21.

    Google Scholar 

  26. 26.

    Gaudio S, Quattrocchi CC. Neural basis of a multidimensional model of body image distortion in anorexia nervosa. Neurosci Biobehav Rev. 2012;36:1839–47.

    PubMed  Article  Google Scholar 

  27. 27.

    Di Paola M, Spalletta G, Caltagirone C. In vivo structural neuroanatomy of corpus callosum in Alzheimer’s disease and mild cognitive impairment using different MRI techniques: A review. J Alzheimer’s Dis. 2010;20:67–95.

    Article  Google Scholar 

  28. 28.

    Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28:30–42.

    PubMed  Article  Google Scholar 

  29. 29.

    Iwashiro N, Suga M, Takano Y, Inoue H, Natsubori T, Satomura Y, et al. Localized gray matter volume reductions in the pars triangularis of the inferior frontal gyrus in individuals at clinical high-risk for psychosis and first episode for schizophrenia. Schizophr Res. 2012;137:124–31.

    PubMed  Article  Google Scholar 

  30. 30.

    Jahanshad N, Kochunov PV, Sprooten E, Mandl RC, Nichols TE, Almasy L, et al. Multi-site genetic analysis of diffusion images and voxelwise heritability analysis: A pilot project of the ENIGMA-DTI working group. Neuroimage. 2013;81:455–69.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Klein A, Ghosh SS, Bao FS, Giard J, Häme Y, Stavsky E, et al. Mindboggling morphometry of human brains. PLoS Comput Biol. 2017;13:e1005350.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  32. 32.

    Jansen IE, Savage JE, Watanabe K, Bryois J, Williams DM, Steinberg S, et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet. 2019;51:404–13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Demontis D, Walters RK, Martin J, Mattheisen M, Als TD, Agerbo E, et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. 2019;51:63–75.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Watson HJ, Yilmaz Z, Thornton LM, Hübel C, Coleman JRI, Gaspar HA, et al. Genome-wide association study identifies eight risk loci and implicates metabo-psychiatric origins for anorexia nervosa. Nat Genet. 2019;51:1207–14.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Otowa T, Hek K, Lee M, Byrne EM, Mirza SS, Nivard MG, et al. Meta-analysis of genome-wide association studies of anxiety disorders. Mol Psychiatry. 2016;21:1391–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Grove J, Ripke S, Als TD, Mattheisen M, Walters RK, Won H, et al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet. 2019;51:431–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Stahl EA, Breen G, Forstner AJ, McQuillin A, Ripke S, Trubetskoy V, et al. Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat Genet. 2019;51:793–803.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Jansen PR, Watanabe K, Stringer S, Skene N, Bryois J, Hammerschlag AR, et al. Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways. Nat Genet. 2019;51:394–403.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Howard DM, Adams MJ, Clarke TK, Hafferty JD, Gibson J, Shirali M, et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat Neurosci. 2019;22:343–52.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Arnold PD, Askland KD, Barlassina C, Bellodi L, Bienvenu OJ, Black D, et al. Revealing the complex genetic architecture of obsessive-compulsive disorder using meta-analysis. Mol Psychiatry. 2018;23:1181–8.

    CAS  Article  Google Scholar 

  41. 41.

    Nievergelt CM, Maihofer AX, Klengel T, Atkinson EG, Chen C-Y, Choi KW, et al. International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci. Nat Commun. 2019;10:4558.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  42. 42.

    Ripke S, Neale BM, Corvin A, Walters JTR, Farh KH, Holmans PA, et al. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421–7.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  43. 43.

    Yu D, Sul JH, Tsetsos F, Nawaz MS, Huang AY, Zelaya I, et al. Interrogating the Genetic Determinants of Tourette’s Syndrome and Other Tic Disorders Through Genome-Wide Association Studies. Am J Psychiatry. 2019;176:217–27.

    PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    1000 Genomes Project Consortium. A global reference for human genetic variation. Nature. 2015;526:68–74.

    Article  CAS  Google Scholar 

  45. 45.

    Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408.

  46. 46.

    Brion M-JA, Shakhbazov K, Visscher PM. Calculating statistical power in Mendelian randomization studies. Int J Epidemiol. 2013;42(5):1497-501.

  47. 47.

    Hartwig FP, Smith GD, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46:1985–98.

    PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.

    PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Wu Y, Zeng J, Zhang F, Zhu Z, Qi T, Zheng Z, et al. Integrative analysis of omics summary data reveals putative mechanisms underlying complex traits. Nat Commun. 2018;9:918.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  50. 50.

    Ucker GR¨, Schwarzer G, Carpenter JR, Binder H, Schumacher M. Treatment-effect estimates adjusted for small-study effects via a limit meta-analysis. Biostatistics. 2011;12:122–42.

    Article  Google Scholar 

  51. 51.

    Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: An R library for genome-wide association analysis. Bioinformatics. 2007;23:1294–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

About this supplement

This article has been published as part of BMC Bioinformatics Volume 22 Supplement 3, 2021: 19th International Conference on Bioinformatics 2020 (InCoB2020): genomics. The full contents of the supplement are available online at https://bmcgenomics.biomedcentral.com/articles/supplements/volume-22-supplement-3.

Funding

The study was supported by National Natural Science Foundation of China (No: 81971292 and 81671328); Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (Grant No. 1610000043); Innovation Research Plan supported by Shanghai Municipal Education Commission (Grant No. ZXWF082101). The funding bodies had no role in the study and collection, analysis, and interpretation of data and in writing the manuscript. Publication was funded by National Natural Science Foundation of China (No: 81971292).

Author information

Affiliations

Authors

Contributions

G.N.L designed and supervised the study. S.Y, W.W collected the data. W.Q and W.S preprocessed and analyzed the data. W.S, W.Q, W.W, and S.Y interpreted the data. W.S wrote the manuscript. All authors read and approved the manuscript.

Corresponding author

Correspondence to Guan Ning Lin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Summary statistics for all instrument variables after data harmonization.

Additional file 2.

Full results for DTI-Disease association analysis.

Additional file 3.

Full results for ROI-Disease association analysis.

Additional file 4.

Funnel plots for top MR findings, leave-one-out results for top DTI-Disorder pairs, general contribution of DTI to neuropsychiatric disease after sensitivity test adjustment, Rank-Rank overlaps between MR effect and heterogeneity for ROI-Disease analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Song, W., Qian, W., Wang, W. et al. Mendelian randomization studies of brain MRI yield insights into the pathogenesis of neuropsychiatric disorders. BMC Genomics 22, 342 (2021). https://doi.org/10.1186/s12864-021-07661-8

Download citation

Keywords

  • Neuroimaging
  • Neuropsychiatric disorders
  • Dysconnectivity
  • Anorexia nervosa
  • Superior longitudinal fasciculus