Causal effect of C-reactive protein and vitamin D on human cerebral anatomy observed among genetically correlated biomarkers in blood

Psychiatric disorders such as schizophrenia are commonly associated with structural brain alterations affecting the cortex, which frequently vary with clinically relevant factors including antipsychotic treatment, duration of illness and age of onset. While the underlying variables mediating these structural changes are poorly understood, recent genetic evidence suggests circulating metabolites and other biochemical traits play a causal role in a number of psychiatric disorders which could be mediated by changes in the cerebral cortex. In the current study, we leveraged publicly available genome-wide association study (GWAS) data to explore shared genetic architecture and evidence for causal relationships between a panel of 50 biochemical traits and measures of cortical thickness and surface area at both the global and regional levels. Linkage disequilibrium score regression identified a total of 20 significant and 156 suggestive genetically correlated biochemical-cortical trait pairings, of which six exhibited strong evidence for causality in a latent causal variable model. Interestingly, a negative causal relationship was identified between a unit increase in serum C-reactive protein levels and thickness of the lingual and lateral occipital regions that was also supported by Mendelian randomisation, while circulating vitamin D (25-hydroxyvitamin D) levels exhibited a positive causal effect on temporal pole thickness. Taken together, our findings suggest a subset of biochemical traits exhibit shared genetic architecture and potentially causal relationships with cortical thickness in functionally distinct regions, which may contribute to alteration of cortical structure in psychiatric disorders.

and evidence for causal relationships between a panel of 50 biochemical traits and measures of cortical thickness 23 and surface area at both the global and regional levels. Linkage disequilibrium score regression identified a total 24 of 20 significant and 156 suggestive genetically correlated biochemical-cortical trait pairings, of which six 25 exhibited strong evidence for causality in a latent causal variable model. Interestingly, a negative causal 26 relationship was identified between a unit increase in serum C-reactive protein levels and thickness of the lingual 27 and lateral occipital regions that was also supported by Mendelian randomisation, while circulating vitamin D 28 (25-hydroxyvitamin D) levels exhibited a positive causal effect on temporal pole thickness. Taken together, our 29 findings suggest a subset of biochemical traits exhibit shared genetic architecture and potentially causal 30 relationships with cortical thickness in functionally distinct regions, which may contribute to alteration of cortical 31 structure in psychiatric disorders.

32
. CC-BY 4.0 International license It is made available under a perpetuity.
is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 NOTE: This preprint reports new research that has not been certified by peer review and should not be used to guide clinical practice.

33
The pathogenesis of psychiatric disorders is underpinned by a complex interplay of genetic and environmental 34 risk factors. Large-scale genetic studies have identified a strong genetic component of these disorders, 35 characterised by a vast polygenic burden across the genome arising from variants ranging in frequency from 36 common to ultra-rare, as well as large effect size structural variants [1,2]. Cross-disorder analyses have 37 additionally revealed substantial proportions of shared genetic architecture between psychiatric disorders [3][4][5], 38 potentially contributing to similarities in clinical presentation and the frequency of observed psychiatric 39 comorbidities. In conjunction with shared genetic risk, many psychiatric disorders are also associated with 40 alteration of brain structure which may account for a range of psychiatric symptoms and offer utility in dissecting 41 underlying pathophysiological mechanisms. Widespread dysregulation of cortical structures in particular has 42 emerged as a consistent feature, with recent neuroanatomical meta-analyses observing grey matter reductions in 43 overlapping cortical regions for schizophrenia, major depressive disorder, bipolar disorder and obsessive-44 compulsive disorder, amongst others, while structural alterations unique to each disorder have also been identified 45 [6][7][8]. Although further investigation of these neuroanatomical changes may assist in the discovery of clinically 46 actionable pathways and biomarkers, the fundamental genetic and environmental factors associated with 47 alterations to brain structure remain poorly understood.

48
Dysregulation of circulating biochemical factors is one broad mechanism through which variations in cortical 49 structure may arise in psychiatric and neurodegenerative disorders. During CNS development, alteration of these 50 circulating factors could plausibly interfere with patterns of neuronal differentiation and migration important for 51 cortical development, while dysregulation in the adult brain may disrupt neuronal cytoarchitecture and integrity.

52
Systemic effects have been widely studied in psychiatry as many of these biochemical variables can be modulated 53 through existing drugs or lifestyle intervention, potentially informing novel treatment interventions. For example, 54 observational evidence suggests elevated C-reactive protein (CRP) in the serum of individuals with schizophrenia 55 is associated with cortical thinning in the frontal, insula and temporal regions [9]. An array of studies has identified 56 many other blood-based biochemical traits with potential diagnostic or prognostic utility [10,11], however a major 57 drawback of these observational studies is their inability to discriminate between correlations and causal 58 relationships, thus any biological effects resulting in alteration of brain structure are difficult to detect. To this 59 end, genome-wide association studies (GWAS) are proving increasingly valuable for examining and 60 distinguishing genetic correlations and genetically informed causal relationships amongst traits of interest. In 61 particular, germline genetic variants are largely fixed at birth which helps to mitigate effects due to reverse 62 causation in observational studies [12]. Indeed, recent studies utilising GWAS-guided methods of causal inference 63 such as Mendelian randomization have uncovered putative causal relationships between blood-based biochemical 64 traits and psychiatric disorders, using genome-wide significant SNPs as genetic proxies for biochemical exposures 84 variables correcting for the use of multiple scanners.

85
We additionally obtained summary statistics for a panel of 50 blood-based biochemical traits produced from a 86 cohort of > 300,000 individuals in the UK biobank (http://www.nealelab.is/uk-biobank). These traits included 87 blood cell counts, metabolites, enzymes, lipids and other biomarkers, all of which have high or medium confidence 88 SNP heritability estimates significantly different from zero. See Table S1 for further details.

90
Genetic correlation 91 Genetic correlations amongst cortical and biochemical traits were examined via linkage disequilibrium score 92 regression (LDSR), as described in detail previously [17]. Briefly, LDSR estimates genetic covariance between 93 two traits by regressing SNP-level c 2 values -the product of SNP marginal effect sizes from both traits (Z1Z2) -94 with respect to LD scores, which estimate the total LD associated with the SNP of interest. Genetic correlation 95 (rg) is then determined by normalising the genetic covariance by their respective trait heritabilities, with standard 96 errors estimated by jackknifing over 200 blocks of adjacent SNPs. All LDSR analyses were conducted using the 97 ldsc python package [18]. We initially analysed biochemical traits with respect to total cortical SA and mean TH, 98 after which regional genetic correlations were explored. Notably, LDSR is robust to sample overlap, thus enabling 99 accurate rg estimates in the current study despite the presence of UKBB samples in both the biochemical and 100 cortical GWAS. Prior to LDSR analysis, all summary statistics were converted to a standardised "munged" format,

113
To identify genetically correlated trait pairings with evidence for a causal relationship, we employed a latent 114 causal variable (LCV) model to estimate genetic causality as described previously [20]. Briefly, the LCV model 115 assumes a latent variable, L, mediates the genetic correlation between two traits such that if trait 1 is strongly 116 correlated with L, it is assumed to be partially genetically causal for trait 2. To test for partial genetic causality, 117 the mixed fourth moments (cokurtosis) of SNP marginal effect size distributions for each trait are compared to 118 determine whether SNPs affecting trait 1 have proportional effects on trait 2, but not vice versa. For interpretation, 119 the LCV model reports a posterior mean genetic causality proportion ( $ ) of trait 1 on trait 2, wherein $ = 0 120 suggest no causal relationship, $ > 0 implies trait 1 is partially genetically causal for trait 2 and $ < 0 121 suggests trait 2 is partially genetically causal for trait 1. We consider significant | $ | estimates > 0.6 as strong 122 evidence for partial genetic causality as shown previously [20], noting that $ summarises the strength of 123 evidence for a genetically causal relationship, rather than the magnitude of a causal effect. We note that, like 124 LDSR, the LCV model is resistant to statistical inflation arising from sample overlap due to the incorporation of 125 the LDSR intercept. All summary statistics subjected to LCV modelling were "munged" prior to analysis as 126 recommended.

128
Mendelian randomisation 129 Two-sample Mendelian randomization (MR) was employed to further explore and quantify the magnitude of 130 causal relationships amongst trait pairings with evidence for partial genetic causality via LCV modelling. MR 131 leverages genetic instrumental variables (IVs) rigorously associated with an exposure -most commonly 132 independent genome-wide significant SNPs -to estimate the causal effect of the exposure on an outcome. We 133 specifically tested the effect of C-reactive protein (CRP) and vitamin D (measured as 25-hydroxyvitamin D) 134 exposures on regional cortical TH measurements utilizing LD-clumped, non-palindromic IVs sourced from well-135 powered non-UK Biobank GWAS [21,22], as most two-sample MR methods are sensitive to sample overlap. For 136 both GWAS, the IV F-statistic was > 10, indicating IVs for both exposures are sufficiently powered [23]. All 137 exposure-outcome trait pairings were analysed via five MR models with differing underlying assumptions 138 regarding the validity of using SNPs as IVs. Specifically, the inverse variance weighted model with multiplicative 139 random effects (IVWm) was used as our principal model, which is generally considered the most-well powered 140 model but has a zero percent breakdown point, and thus, assumes all IVs are valid [24]. We also incorporated an 141 IVW model with fixed effects (IVWf), which is a less conservative iteration of the IVW estimator that accounts 142 for inter-instrument heterogeneity in a more-general manner, and thus is more statistically valid at the risk of 143 additional bias in the presence of heterogeneity. To account for potential violations of the "all-valid" assumption,

144
we additionally utilised a weighted median model, which through the plurality valid assumption is said to be an 145 unbiased estimator if at least 50% of the weight in the model is derived from valid IVs [25]. A weighted mode 146 model was additionally employed, which also tests the plurality valid assumption and is more-robust to outlying 147 IVs at the expense of power [26]. Finally, the MR Egger approach was also utilised which does not constrain the 148 intercept as zero, with a significantly non-zero intercept considered to be evidence of potential unbalanced is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 of global pleiotropy that is also related to heterogeneity [29]. Finally, the Steiger directionality test was 153 implemented to evaluate evidence that the assumed causal direction (biochemical trait → cortical property) was

184
We then explored site-specific genetic correlations throughout the cortex by analysing measures of SA and TH 185 for 34 distinct regions. For regional SA, 18 traits pairings survived multiple testing correction while an additional 186 120 suggestive pairings were identified, involving a total of 18 biochemical traits and 32 cortical regions (Fig. 1b, 187 Table 2, Table S3). Notably, 95.8% of these correlations were negative, with particularly strong representation 188 from white blood cell (WBC) count (26), neutrophil count (23) and monocyte count (12), while the superior frontal 189 (10), middle temporal (9), medial orbitofrontal (8) and inferior parietal (8) regions were most consistently affected.

190
The extent of genetic correlation for regional TH was comparatively modest, with two significant and 33 191 suggestive trait pairings identified involving 11 biochemical traits and 15 cortical phenotypes (Fig. 1b is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint 6 Table S3). Although 72.7% of correlations were positive, many of these involved haematocrit percentage (8), 193 haemoglobin (6) and red blood cell count (5), which were also suggestively correlated with global mean cortical 194 TH (Fig. 1a). C-reactive protein, however, exhibited no response at the global level but was negatively correlated  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint

215
(bottom) of 34 cortical regions as defined by the Desikan-Killiany atlas.  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint

217
The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint 9 reticulocyte traits) tended to cluster in both the SA and TH analyses (Fig. S1). In contrast to the cortical 230 phenotypes, these traits exhibited a wider range of both positive and negative correlations.

231
We subsequently employed finite Gaussian mixture modelling (GMM) to identify latent genetic relationships 232 between cortical measures and biochemical traits. Two distinct components (clusters) of cortical regions with 233 respect to their biochemical correlations were identified in the TH analysis, with cluster 1 (17 regions) largely 234 localised to the frontal and temporal lobes, while cluster 2 (17 regions) was predominantly parietal and occipital 235 ( Fig. 2b & c, Table S4). Although two clusters were also observed in the SA analysis, the spatial distribution was 236 less clear, as cluster 1 (19 regions) contained a series of contiguous regions spanning a number of cortical lobes, 237 whereas cluster 2 (15 regions) consisted of regions with relatively diffuse spatial localisation ( Fig. 2d & e, Table   238 S5). These results therefore suggest that genetic correlations between regional cortical structure and biochemical 239 traits exhibit a degree of spatial organisation.

240
GMM applied to cortical correlation profiles of the biochemical traits failed to identify discrete clusters of 241 biomarkers for either global SA or TH, as the optimum BIC value yielded only a single component. Therefore,

242
we conducted GMM for the biochemical traits after combining regional SA and TH correlation profiles to capture 243 groups of biomarkers with similar effects on both cortical properties. Interestingly, seven clusters of biochemical 244 traits were observed (Fig. 2f, Table S6). Clusters one and six were largely composed of traits with negative LDSR 245 Z-scores for regional SA and positive scores for TH, including cholesterol, testosterone and RBC-related traits, is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint

287
Strong evidence for causal relationships between blood-based biomarkers and regional cortical thickness

288
We next examined whether any correlated trait pairings exhibited evidence for partial genetic causality using a  Table 3, Table S7). To ensure these results were 301 not affected by global measures, we repeated this analysis using GWAS summary statistics covaried for mean 302 global cortical TH. Strong evidence for partial genetic causality remained after this adjustment for CRP on lingual 303 TH ( $ = 0.69, SE = 0.19, P = 1.41x10 -73 ) and lateral occipital TH ( $ = 0.6, SE = 0.21, P = 5.93x10 -6 ) and 304 vitamin D on temporal pole TH ( $ = 0.72, SE = 0.18, P = 2.63x10 -7 ; Table 3, Table S7). Considering the sign 305 of genetic correlations between these traits, it is therefore likely that increased CRP levels are associated with is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint Table 3, Table S7). Modest $ estimates were also obtained for these biochemical markers with respect to global 311 cortical TH: CRP ( $ = 0.34, SE = 0.28, P = 0.33), RBC count ( $ = 0.51, SE = 0.16, P = 1.09x10 -6 ), and 312 haematocrit percentage ( $ = 0.39, SE = 0.26, P = 0.16; Table S8). We thus suspect the strong $ estimates 313 obtained for these biochemical traits with respect to uncorrected regional TH GWAS quantified the combined 314 effect on regional and global TH.  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint

330
The total effect of CRP levels (natural log transformed mg/L) on lingual and lateral occipital TH (in mm) was 331 next quantified via two-sample Mendelian randomisation (MR). Due to the inclusion of UKBB samples in the 332 ENIGMA metanalysis, we obtained instrumental variables (IVs) for CRP utilising summary statistics sourced 333 from a non-UKBB cohort [21]. Furthermore, we employed an inverse variance weighted estimator with 334 multiplicative random effects (IVWm) as our primary test, followed by four other methods with different 335 underlying assumptions regarding IV validity. Across all analyses, a total of 48 CRP-associated IVs were 336 identified after harmonisation and removal of palindromic SNPs. Interestingly, we identified evidence suggesting 337 each natural log transformed mg/L increase in CRP was associated with a significant reduction in lingual (bIVWm=      Table 4, Table S9). Across 347 all analyses, we identified no substantial evidence for outlier SNPs using MR-PRESSO and leave-one-out 348 analyses, although the IV exposure-outcome effects exhibited evidence for SNP heterogeneity in both instances 349 (Cochran's Q ≥ 69.86, P ≤ 0.01; Tables S10 & S11). This heterogeneity, however, is not unexpected given the 350 underlying biological complexity of these traits. Moreover, the Egger intercepts were not significantly different 351 than zero (P ≥ 0.26), therefore providing statistical evidence of no unbalanced pleiotropy (Table S12).  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 14 370 371 is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 Scatter plots comparing IV effect sizes in exposure and outcome GWAS. Each trendline corresponds to one of 375 five MR methods utilized for sensitivity analysis. (d) As in (c), except using cortical GWAS corrected for global 376 TH. Table 4. Causal relationships between biochemical traits and regional cortical TH estimated via two-sample 379 Mendelian randomisation.

384
Limited transcriptome-wide overlap amongst causally linked trait pairings

385
Transcriptome-wide association studies were performed for CRP, vitamin D and their respective causally 386 associated cortical phenotypes to compare and contrast predicted gene expression profiles. Briefly, this approach 387 utilises tissue-specific cis-eQTL expression weights to impute gene expression profiles associated with the 388 phenotype of interest. We specifically employed pre-computed SNP-expression weights for cortical tissue and 389 whole blood, while liver expression weights were additionally included as CRP is hepatically synthesised (see 390   Table S13 for full results). At a PBonferroni < 0.05 cut-off, the average number of associated genes per trait across 391 all tissues was as follows: CRP = 157, vitamin D = 37, lateral occipital TH = 1, lingual TH = 4 and temporal pole 392 TH = 2 (see Fig. 5a for representative Miami plots, see Table S14 for full results). Given the sparse identification 393 of associated genes for the cortical traits, no overlapping genes were identified (Table 5, Table S15). However,

394
after correcting for global TH, RP11-148021.6, was found to overlap between CRP (ZTWAS = 6.03) and lateral 395 occipital TH (ZTWAS = -4.46) using cortical expression weights ( is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint Z-scores in the CRP TWAS and positive Z-scores in the lingual TH TWAS. To more-broadly explore 401 transcriptomic overlap amongst trait pairings with evidence for causality, we utilised RHOGE to estimate 402 transcriptome-wide correlation (rGE). For all pairings and tissues, no significant correlations were identified after 403 multiple testing correction (Fig. 5b, Table S16), however, at a nominal P < 0.05 cut-off vitamin D and temporal 404 pole TH were positively correlated after correction for global measures using cortical weights (rGE = 0.45, P = 405 0.023, SE = 0.187). Collectively, these findings suggest there is minimal predicted transcriptomic overlap between 406 causal biochemical-cortical trait pairings. However, we also identified only a limited number of transcriptome-407 wide associated signals for the cortical properties suggesting that greater sample sizes are required, as well as 408 more tissue and cell-type specific expression weights. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint

410
The copyright holder for this this version posted September 15, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint   is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021 (mCRP) at the site of insult, which increase the abundance of M1 macrophages and Th1 T-cells, thereby activating 462 a robust immune response with potentially deleterious effects on host tissues [46]. Persistent elevation of CRP is 463 therefore strongly associated with chronic, low-level inflammation thought to negatively impact the structural 464 integrity of brain regions. For instance, elevated CRP is positively correlated with cortical thinning [9, 47] and 465 decreased grey matter volume [48,49], while deposition of mCRP has been identified at sites of 466 neurodegeneration [50]. Although CRP-mediated thinning of the lingual and lateral occipital regions has not been 467 explicitly reported in previous studies, recent work has shown elevated CRP is associated with lower regional

484
CRP and lingual TH. We additionally note that reverse causality is unlikely for these putative relationships given 485 the strong directional evidence identified via LCV.

486
Mounting observational evidence suggests psychiatric disorders including schizophrenia, major depressive 487 disorder and bipolar disorder are associated with elevated CRP levels [9,52,53] is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint 20 a protective effect against underlying chronic inflammation. This is an issue for observational inferences in 502 particular, which employ tests such as the mean difference in CRP between cases and controls, as the magnitude 503 of CRP elevation can often be relatively small and more indicative of this moderate state. Regardless of the 504 underlying mechanism, these findings collectively advocate further investigation of the direct and indirect effects 505 of CRP on neuronal integrity to further reconcile the impact of CRP levels on cortical structure and psychiatric 506 disorders.

507
We additionally identified a positive causal relationship for vitamin D on temporal pole TH after correcting for 508 global TH. A nominally significant transcriptome-wide correlation between these traits was also observed after 509 correcting for global TH. Together, these findings suggest any effect of vitamin D was highly specific to the 510 temporal pole, broadly consistent with a recent cross-sectional study wherein total vitamin D intake and vitamin 511 D supplementation were specifically associated with enhanced temporal lobe TH in cognitively normal, older

517
[62], and exhibits neuroprotective effects with respect to excessive calcium, reactive oxygen species and 518 corticosterone [63], indicating a particularly important role throughout the brain. While only one significant 519 relationship was identified for vitamin D in the present study, we consequently suspect protective effects of 520 vitamin D on cortical regions other than the temporal pole may be mediated through complex multivariate 521 interactions not readily detectable in our analyses.

522
The proposed relationship between vitamin D and temporal pole TH is nonetheless significant for general is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ; https://doi.org/10.1101/2021.09.11.21263146 doi: medRxiv preprint on brain development and schizophrenia risk [76] will likely prove critical in mapping the impact of vitamin D in 541 development of key cortical regions such as the temporal pole.

542
In summary, our findings suggest subsets of biochemical exposures and cortical structural properties share genetic 543 architecture and, in some cases, exhibit evidence for causal relationships. We acknowledge a number of limitations 544 and caveats with respect to interpretation of the presented data. Firstly, all analyses in the current study are 545 inherently subject to limitations and biases potentially associated with the utilised summary statistics, such as 546 population stratification [78] and selection bias [79]. The UKBB cohort in particular is predominantly composed 547 individuals over the age of 40, thus age-related modulation of variants may affect genetic correlations and causal 548 estimates. Secondly, while many of the reported genetic correlations likely reflect shared genetic architecture, 549 horizontal pleiotropy may mediate associations between some trait pairings, however we note these are still useful 550 in guiding identification of genes that impact both biochemical traits and cortical structure. Finally, we caution 551 that all reported causal relationships require validation in randomised controlled trials to confirm the putative 552 causal effects. Despite the limitations of these analyses, genetically informed causal inference represents an 553 exciting opportunity to screen and prioritise biochemical traits en masse to guide future investigation of these 554 exposures in the context of neuronal function, brain cytoarchitecture and psychiatric illness.  is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted September 15, 2021. ;https://doi.org/10.1101https://doi.org/10. /2021