# Immune-Related Genetic Overlap Between Regional Gray Matter Reductions and Psychiatric Symptoms in Adolescents, and Gene-Set Validation in a Translational Model

Lukas Penninck, El Chérif Ibrahim, Eric Artiges, Victor Gorgievski, Sylvane Desrivières, Severine Farley, Irina Filippi, Carlos E A de Macedo, Raoul Belzeaux, Tobias Banaschewski, et al.

► To cite this version:

Lukas Penninck, El Chérif Ibrahim, Eric Artiges, Victor Gorgievski, Sylvane Desrivières, et al. Immune-Related Genetic Overlap Between Regional Gray Matter Reductions and Psychiatric Symptoms in Adolescents, and Gene-Set Validation in a Translational Model. Frontiers in Systems Neuroscience, 2021, 15, 10.3389/fnsys.2021.725413. [https://doi.org/10.3389/fnsys.2021.725413](https://doi.org/10.3389/fnsys.2021.725413). [arXiv:2106.07254](https://arxiv.org/abs/2106.07254)

HAL Id: hal-03801037

https://hal.science/hal-03801037

Submitted on 6 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Immune-Related Genetic Overlap Between Regional Gray Matter Reductions and Psychiatric Symptoms in Adolescents, and Gene-Set Validation in a Translational Model

## OPEN ACCESS

**Edited by:**
Wen-Jun Gao,
Drexel University, United States

**Reviewed by:**
Wang Jiesi,
University of Chinese Academy
of Sciences, China
Adriana Caballero,
University of Illinois at Chicago,
United States

**\*Correspondence:**
Jean-Luc Martinot
jean-luc.martinot@inserm.fr

†These authors have contributed
equally to this work

**Received:** 15 June 2021
**Accepted:** 31 August 2021
**Published:** 30 September 2021

**Citation:**
Penninck L, Ibrahim EC, Artiges E, Gorgievski V, Desrivières S, Farley S, Filippi I, de Macedo CEA, Belzeaux R, Banaschewski T, Bokde ALW, Quinlan EB, Flor H, Grigis A, Garavan H, Gowland P, Heinz A, Brühl R, Nees F, Papadopoulos Orfanos D, Paus T, Poustka L, Fröhner JH, Smolka MN, Walter H, Whelan R, Grenier J, Schumann G, Paillère Martinot M-L, Tzavara ET and Martinot J-L (2021) *Immune-Related Genetic Overlap Between Regional Gray Matter Reductions and Psychiatric Symptoms in Adolescents, and Gene-Set Validation in a Translational Model*. *Front. Syst. Neurosci.* 15:725413.
doi: 10.3389/fnsys.2021.725413*

Lukas Penninck<fnref n="1" />†, El Chérif Ibrahim<fnref n="2" />†, Eric Artiges<fnref n="1" />,³, Victor Gorgievski<fnref n="4" />, Sylvane Desrivières⁵, Severine Farley<fnref n="4" />, Irina Filippi<fnref n="1" />, Carlos E. A. de Macedo<fnref n="4" />, Raoul Belzeaux<fnref n="2" />,<fnref n="6" />, Tobias Banaschewski⁷, Arun L. W. Bokde⁸, Erin Burke Quinlan⁵, Herta Flor⁹,¹⁰, Antoine Grigis¹¹, Hugh Garavan<fnref n="12" />, Penny Gowland¹³, Andreas Heinz¹⁴, Rüdiger Brühl¹⁵, Frauke Nees⁷,⁹,¹⁶, Dimitri Papadopoulos Orfanos¹¹, Tomáš Paus¹⁷,¹⁸, Luise Poustka¹⁹†, Juliane H. Fröhner²⁰, Michael N. Smolka²⁰, Henrik Walter¹⁴, Robert Whelan²¹, Julien Grenier²², Gunter Schumann⁵,²³,²⁴, Marie-Laure Paillère Martinot<fnref n="1" />,²⁵†, Eleni T. Tzavara<fnref n="4" />,<fnref n="6" />,²⁶† and Jean-Luc Martinot<fnref n="1" />*† for the IMAGEN Consortium

Adolescence is a period of vulnerability for the maturation of gray matter (GM) and also for the onset of psychiatric disorders such as major depressive disorder (MDD), bipolar disorder and schizophrenia. Chronic neuroinflammation is considered to play a role in the etiology of these illnesses. However, the involvement of neuroinflammation in the observed link between regional GM volume reductions and psychiatric symptoms is not

<footnote n="1">¹ Institut National de la Santé et de la Recherche Médicale, INSERM U1299 "Trajectoires Développementales en Psychiatrie", Université Paris-Saclay, Ecole Normale Supérieure Paris-Saclay, CNRS, Centre Borelli, Gif-sur-Yvette, France,</footnote>

<footnote n="2">² Aix Marseille Univ, CNRS, INT, Inst Neurosci Timone, Marseille, France, ³ EPS Barthelemy Durand, Etampes, France,</footnote>

<footnote n="4">⁴ University of Paris, CNRS, INCC, Paris, France, ⁵ Centre for Population Neuroscience and Precision Medicine (PONS),</footnote>

<footnote>Institute of Psychiatry, Psychology & Neuroscience, SGDP Centre, King's College London, London, United Kingdom,</footnote>

<footnote n="6">⁶ AP-HM, Hôpital Sainte Marguerite, Pôle de Psychiatrie Universitaire Solaris, Marseille, France, ⁷ Department of Child</footnote>

<footnote>and Adolescent Psychiatry and Psychotherapy, Central Institute of Mental Health, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany, ⁸ Discipline of Psychiatry, School of Medicine and Trinity College Institute of Neuroscience,</footnote>

<footnote>Trinity College Dublin, Dublin, Ireland, ⁹ Institute of Cognitive and Clinical Neuroscience, Central Institute of Mental Health,</footnote>

<footnote>Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany, ¹⁰ Department of Psychology, School of Social</footnote>

<footnote>Sciences, University of Mannheim, Mannheim, Germany, ¹¹ NeuroSpin, CEA, Université Paris-Saclay, Gif-sur-Yvette, France,</footnote>

<footnote n="12">¹² Department of Psychiatry and Psychology, University of Vermont, Burlington, VT, United States, ¹³ Sir Peter Mansfield</footnote>

<footnote>Imaging Centre, School of Physics and Astronomy, University of Nottingham, Nottingham, United Kingdom, ¹⁴ Department</footnote>

<footnote>of Psychiatry and Psychotherapy, CCM, Charité – Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin,</footnote>

<footnote>Humboldt-Universität zu Berlin, and Berlin Institute of Health, Berlin, Germany, ¹⁵ Physikalisch-Technische Bundesanstalt</footnote>

<footnote>(PTB), Braunschweig, Germany, ¹⁶ Institute of Medical Psychology and Medical Sociology, University Medical Center</footnote>

<footnote>Schleswig Holstein, Kiel University, Kiel, Germany, ¹⁷ Department of Psychiatry, Faculty of Medicine and Centre Hospitalier</footnote>

<footnote>Universitaire Sainte-Justine, University of Montreal, Montreal, QC, Canada, ¹⁸ Department of Psychology and Psychiatry,</footnote>

<footnote>University of Toronto, Toronto, ON, Canada, ¹⁹ Department of Child and Adolescent Psychiatry and Psychotherapy,</footnote>

<footnote>University Medical Centre Göttingen, Göttingen, Germany, ²⁰ Department of Psychiatry and Neuroimaging Center,</footnote>

<footnote>Technische Universität Dresden, Dresden, Germany, ²¹ School of Psychology and Global Brain Health Institute, Trinity</footnote>

<footnote>College Dublin, Dublin, Ireland, ²² Université de Paris, INSERM UMRS 1124, Paris, France, ²³ PONS Research Group,</footnote>

<footnote>Department of Psychiatry and Psychotherapy, Humboldt University, Berlin and Leibniz Institute for Neurobiology,</footnote>

<footnote>Magdeburg, Germany, ²⁴ Institute for Science and Technology of Brain-inspired Intelligence (ISTBI), Fudan University,</footnote>

<footnote>Shanghai, China, ²⁵ AP-HP-Sorbonne Université, Department of Child and Adolescent Psychiatry, Pitié-Salpêtrière Hospital,</footnote>

<footnote>Paris, France, ²⁶ Fondation Fondamental, Créteil, France</footnote>
established yet. Here, we investigated a possible common immune-related genetic link between these two phenomena in european adolescents recruited from the community. Hippocampal and medial prefrontal cortex (mPFC) were defined a priori as regions of interest (ROIs). Their GM volumes were extracted in 1,563 14-year-olds from the IMAGEN database. We found a set of 26 SNPs that correlated with the hippocampal volumes and 29 with the mPFC volumes at age 14. We formed two ROI-Related Immune-gene scores (RRI) with the inflammation SNPs that correlated to hippocampal GM volume and to mPFC GM volume. The predictive ability of both RRIs with regards to the presence of psychiatric symptoms at age 18 was investigated by correlating the RRIs with psychometric questionnaires obtained at age 18. The RRIs (but not control scores constructed with random SNPs) correlated with the presence of depressive symptoms, positive psychotic symptoms, and externalizing symptoms in later adolescence. In addition, the effect of childhood maltreatment, one of the major environmental risk factors for depression and other mental disorders, interacted with the RRI effect. We next sought to validate this finding by investigating our set of inflammatory genes in a translational animal model of early life adversity. Mice were subjected to a protocol of maternal separation at an early post-natal age. We evaluated depressive behaviors in separated and non-separated mice at adolescence and their correlations with the concomitant expression of our genes in whole blood samples. We show that in mice, early life adversity affected the expression of our set of genes in peripheral blood, and that levels of expression correlated with symptoms of negative affect in adolescence. Overall, our translational findings in adolescent mice and humans provide a novel validated gene-set of immune-related genes for further research in the early stages of mood disorders.

**Keywords: immunity genes, psychiatric symptoms, adolescence, MRI, childhood maltreatment**

INTRODUCTION

Some large-scale studies combining genetic and brain structural data hypothesize the existence of shared neurobiological mechanisms underlying prevalent psychiatric disorders, such as major depressive disorder (MDD), attention-deficit/hyperactivity disorder, and schizophrenia (Parker et al., 2020; Patel et al., 2021). The central findings supporting this hypothesis are the associations found between disorder-specific regional differences in brain structure (e.g., cortical thickness or regional volumes) and common clusters of genes involved in brain development or maturation. Although they do not establish causality, these observations point to the interplay of genetic and brain structural underpinnings in the pathophysiology of psychiatric illnesses. Here, we aim to further contribute to this endeavor by applying a targeted approach, i.e., by focusing on a limited number of genes and brain regions. The main advantage of such an approach is the improvement of statistical power to identify associations in smaller samples. Specifically, we will investigate the possible association between neuroinflammatory-related genes and regional gray matter (GM) volumes in the hippocampus and medial prefrontal cortex (mPFC) in the development of mood disorders (MDD and bipolar disorder) and schizophrenia.

Although a wide range of structural abnormalities has been associated with psychiatric disorders (e.g., Patel et al., 2021), herein we focused on volumetric GM measurements in the hippocampus and mPFC as regions of interest (ROI). Indeed, reduced hippocampal and PFC volume are among the most replicated findings in MRI studies of depression (MacQueen et al., 2003; Schmaal et al., 2020). Lower hippocampal volumes have been associated with adolescent onset MDD (Chen et al., 2010; Schmaal et al., 2016), and ventral medial PFC maturation has been related to negative affect in the developing brain (Ducharme et al., 2014). GM reductions in these two regions have been put forward as indicators of the severity and stage of MDD (Belleau et al., 2019). Still, in their review, Belleau et al. (2019) speculate that, although hippocampal and mPFC GM reductions have been associated with MDD, these reductions are neither necessary nor sufficient for inducing a depressive episode. Instead, these structural abnormalities should be regarded as intermediary effectors driving the progression and recurrence of depression. Consistently, our group has reported lower volumes in both regions as variables of interest for tackling irregular sleep habits paving the way to psychiatric symptoms in adolescents (Lapidaire et al., 2021).
Our second reduction in scope, i.e., focusing on a carefully
selected set of immune-related genes rather than a genome-
wide paradigm, is founded on research putting forward
chronic neuroinflammation as a neurobiological characteristic
of MDD driving GM loss (Kubera et al., 2011; Barnes et al.,
2017). This theory, originally called the “inflammatory and
neurodegenerative hypothesis of depression” by Maes et al.
(2009), is based on multiple pieces of evidence. The first paper
demonstrating a close connection between depression and the
immune system was published in 1990 and found that MDD
was often associated with a significantly higher number of
activated T-cells (Maes et al., 1990). Since then, there have been
many consistent findings of increased levels of proinflammatory
cytokines in the cerebrospinal fluid (CSF) of patients with
depression, the most prominent being interleukin-1 (IL-1), IL-
2, IL-6, IL-8, IL-12, interferon-γ and tumor necrosis factor-
α. In addition, elevations in peripheral blood concentrations
of chemokines, adhesion molecules, acute phase proteins and
inflammatory mediators such as prostaglandins have been
observed. Lastly, depression could be induced by administrating
cytokines (see the comprehensive reviews by Raison et al., 2006
and Miller and Raison, 2015).

To our knowledge, there is no report investigating the putative association between immune-related genetic variation and MDD-related hippocampal and mPFC GM reductions. In order to do so, we will combine the benefits of a candidate gene approach and a polygenic approach. More specifically, we will construct two “polygenic” scores using only inflammation-related common genetic polymorphisms: one in association with hippocampal GM volume and another in association with mPFC GM volume. Hence, these scores will be referred to as ROI-related immune-gene scores (RRI-scores). The IMAGEN database will be used to access genetic and T1 imaging data in 14-year-olds, an age of particular interest as confounding effects due to the use of medication can be considered minimal. A second notable asset of the IMAGEN database is the availability of follow-up data at the age of 18, including a wide range of psychometric data. We will use this follow-up to investigate both RRI-scores with regards to the participants’ psychiatric symptoms later in adolescence. Thus, we will examine whether the putative association between neuroinflammatory genes and regional GM reductions plays a role in the development of psychiatric symptoms.

Moreover, we will test the translational hypothesis that genetic predisposition influences the capacity of an environmental risk factor to induce a psychiatric disorder (Caspi and Moffitt, 2006; Bagot et al., 2014). First, we will assess whether there is an interaction between the RRIs and the degree of childhood maltreatment (CM) explaining negative affects at adolescence, in the IMAGEN database. Second, we will use a translational approach employing an animal model of early life adversity. We will assess whether (i) mice subjected to early life adversity display depressive-like behaviors at adolescence; (ii) the expression of the constructed hippocampal RRI gene-set is altered in peripheral blood in these mice; (iii) transcript levels correlate with the severity of depression-related behavioral scores in adolescent mice. The advantage of this combined approach is that we use transcriptional profiling, which measures the expression

of genes and is sensitive to both genotype and environment,
to gain insight toward the (patho)physiological link between
inflammatory pathways, childhood trauma, and depression
symptoms in adolescence.

MATERIALS AND METHODS

Participants

Participants' datasets were drawn from the IMAGEN project, a European multi-center collaboration combining genetic, neuro-imaging and neuropsychological data from 2223 adolescents at baseline (BL; 14 years old). Participants were followed up 2 years (follow-up 1; FU1) and 4 years later (FU2). An initial sample of 1563 14-year-old adolescents was defined, for which genetic information, T1-weighted MRI images passing the different quality control procedures and the multiple necessary variables were available. In order to perform correlational analyses with psychometric measurements, subgroups of the initial sample were constructed with participants for whom the necessary psychometric data were available. Recruitment procedures have been previously described (Schumann et al., 2010). Written informed consent was obtained from all participants and their legal guardians and verbal assent was obtained from the adolescents.

Neuro-Imaging Data

T1-Weighted MRI

High-resolution T1-weighted anatomical MR images were obtained by means of three Tesla scanners (Philips, Siemens, and General Electric), using a standardized 3D T1-weighted magnetization prepared rapid acquisition gradient echo (MPRAGE) sequence based on the ADNI protocol.<fnref n="1" /> The full details of the MRI acquisition protocols and quality checks have been described previously (Schumann et al., 2010). Image preprocessing was performed with Statistical Parametric Mapping 12 software (SPM12) and its toolbox extension Computational Anatomy Toolbox 12 software (CAT12). In summary, T1-weighted images were segmented and normalized using customized tissue probability maps. Then, the normalized, segmented, and modulated gray matter (GM) and white matter (WM) images were smoothed using a 8-mm full-width at half-maximum Gaussian kernel. Total GM, WM, and cerebrospinal fluid (CSF) volumes were computed for each participant. Total intracranial volume (TIV) was defined as the sum of GM, WM, and CSF volumes. Correct segmentation by CAT12 was verified through visual evaluation of the outliers determined by the automated quality control step “Check Sample Homogeneity” available in CAT12.

Extraction of Regional Gray Matter Volumes

The matlab-script "get_totalsm"<fnref n="2" /> was used to extract the
hippocampal and mPFC GM volumes from the baseline GM-
segmented T1 MRI images. In order to do so, two masks were

<footnote n="1">¹http://adni.loni.usc.edu/methods/mri-analysis/mri-acquisition/</footnote>

<footnote n="2">²http://www0.cs.ucl.ac.uk/staff/g.ridgway/vbm/get_totalsm</footnote>
designed (see **Supplementary Figure 1**) using WFU Pickatlas software (a SPM12 toolbox extension): a bilateral hippocampal mask, available in Pickatlas, and a mPFC mask, composed of the Brodmann areas (BA) 10, 11, 12, 14, 24, 25, 32, and 33. Here the medial prefrontal cortex was defined in its widest sense. For instance, BA11 was added because it pertains to the orbito-frontal cortex in its most medial part; BA10 pertains to the anterior prefrontal pole but it also includes a medial part.

## Genetic Data

### SNP Genotyping

The DNA purification and genotyping procedures implemented in the IMAGEN study have been previously described (Desrivières et al., 2015). Population homogeneity was verified with the Structure software using HapMap populations as reference groups (Pritchard et al., 2000). Further correction for population stratification through principal component analysis was deemed unnecessary. After the quality control measures, genotypic data for a total of 466 125 SNPs were considered. The software PLINK was used to extract the SNP genotypes.

### Construction of RRI-Scores

Based on an extensive literature screening, using the keywords “chronic,” “neuroinflammation,” and “review” in PubMed, genes encoding direct and indirect contributors to neuroinflammation were characterized. To this end, every gene (or protein) that was found to be related to neuroinflammation according to at least three reviews was added to the list. Albeit not systematic, we did consider this procedure to be appropriate for this exploratory study. Including more genes by loosening the constraints might be worth exploring in future research, but will not necessarily lead to RRIs with higher predictive power. The SNPs in and around (±5 kb) the listed immune-related genes were obtained by means of the UCSC Genome Browser<fnref n="3" /> and only those genotyped in the IMAGEN database were selected, a total of 674 SNPs. The methodology used in this study to construct the two RRIs (one associated with hippocampal GM volume, the other with mPFC GM volume) was based on the recent guide put forward by Choi et al. (2020). The effect of every SNP on either the hippocampal or the mPFC GM volume was assessed by performing linear regression analyses in R, using data from the initial sample and the standard *lm* function in R. The dependent variable was the BL GM volume of either the hippocampus or the mPFC, the independent variable was the major allele count for the SNP of interest (0 = homozygous for the minor allele, 1 = heterozygous, 2 = homozygous for the major allele). The regression was controlled for the covariates sex, Puberty Developmental Scale (PDS) score as a proxy of age, TIV and scanner type. Next, the SNPs that correlated (*p* < 0.1 without correction for multiple comparisons) with the hippocampal or mPFC GM volume were selected. In order to control for linkage disequilibrium (LD) and avoid redundancy in the SNPs included in the score, a manual procedure analogous to SNP pruning was carried out. More specifically, the online application LDmatrix developed by the

National Institute of Health<fnref n="4" /> was used to study the LD between the SNPs. Groups of SNPs that were found in high LD (r² > 0.5) were replaced by the most significantly correlated SNP of that group as representative, eliminating the other SNP(s) of that group. Lastly, for every participant, a hippocampal RRI (HRRI) and mPFC RRI (MRRI) were calculated based on the participant's genotype for the group of independent SNPs that were correlated with the hippocampal and mPFC GM volume, respectively. More precisely, the score was defined by the sum of minor alleles for the included SNPs, weighted by the effect size of those alleles individually. A normalization of the HRRI and MRRI values was performed in order to obtain two scores ranging from 0 to 10 that could easily be compared and combined.

### Construction of Control Scores

A collection of 674 SNPs available in the IMAGEN database was randomly selected. Using an identical methodology as described above, a control score explaining hippocampal GM volume and a control score explaining mPFC GM volume were designed using these 674 random SNPs.

## Questionnaire Data

### Questionnaires

Five questionnaire measurements were extracted from the IMAGEN database. First, the algorithmically calculated scores (ranging from 0 to 5) representing the probability of depression according to the DSM-IV (referred to as DepBand) were obtained for participants at FU2 through the Development and Well-Being Assessment (DAWBA), a self-administered diagnostic questionnaire.<fnref n="5" /> Second, the Community Assessment Psychic Experiences-42 (CAPE-42) questionnaire was used to obtain three scores evaluating the presence of depressive symptoms (referred to as the Depressive Dimension Score; DDS) as well as psychotic experiences, both positive (Positive Dimension Score; PDS) and negative (Negative Dimension Score; NDS) at FU2. Third, the self-reported Strengths and Difficulties Questionnaire (SDQ) at FU2 was used to construct two scores representing externalizing and internalizing behaviors; the Externalizing Score (ES) by summing the Conduct Problems Score and Hyperactivity Score, the Internalizing Score (IS) by summing the Emotional Symptoms Score and the Peer Problems Score. Fourth, a score representing childhood maltreatment (CM) was constructed based on information from the Childhood Trauma Questionnaire (CTQ). This retrospective recall-based questionnaire was administered to the participants at BL and, as described in the manual, produces a score (ranging from 0 to 4) representing the endured stress with regards to six categories. Participants were subsequently categorized in five groups based on the highest score in the six subscales. Ultimately, the three highest groups were merged, resulting in a CM score ranging from 0 to 2. Fifth, the Alcohol Use Disorders Identification Test (AUDIT) Score at FU2 was extracted from the database in order to be included as a covariate in linear regression analyses modeling FU2 GM volumes.

<footnote n="3">³https://genome.ucsc.edu</footnote>

<footnote n="4">⁴https://ldlink.nci.nih.gov/?tab=ldmatrix</footnote>

<footnote n="5">⁵http://www.dawba.info</footnote>
## Correlational Analyses

The above-mentioned psychometric measurements (DepBand, DDS, PDS, NDS, ES, and IS) were modeled separately in function of the HRRI, the MRRI and the sum of both scores (HMRRI), as well as the control scores. These regression analyses were performed in R. Since the dependent variables DepBand, ES and IS displayed probability distributions similar to a Poisson distribution, Poisson regression (the log-linear type of the generalized linear model) was opted for. The dependent variables DDS, PDS, and NDS were found normally distributed. However, a log-transformation of the dependent variable was implemented in order to correct for the positive skewness. All regression analyses were controlled for the covariates gender and CM. Also, the interaction between CM and the score was evaluated. P-values were corrected for multiple comparisons through the Benjamini-Hochberg false discovery rate (FDR) procedure, using the *p-adjust* function in R.

## Animals

All experiments on mice were carried out according to policies on the care and use of laboratory animals of European Community legislation 2010/63/EU. The local Ethics Committee (Comité d'éthique en expérimentation animale Charles Darwin N°5) approved the protocols used in this study (protocol number 01486).

The mice were kept under standard conditions at 22 ± 1°C, and a 12-h light-dark cycle with food and water available *ad libitum*.

## Maternal Separation/Maternal Stress Protocol

Pregnant dams (BALB/c Jico) were purchased from Centre d'Elevage Janvier, (Le Genest St Isle, France) to arrive in our facility 5 days before expected delivery. Dams and their respective litters were divided into two groups. The first group (MS, *Maternaly separated*; *n* = 3 dams) was subjected to maternal separation/maternal stress procedure; the second group (NS, *No Separation*; *n* = 2 dams) of dams and respective litters was kept in standard housing conditions as controls.

The maternal separation/maternal stress protocol was adapted from Franklin et al. (2010). The protocol combined (i) physical separation of the pups from the mother and among them; (ii) a short maternal stress at the end of the separation period; (iii) unpredictability regarding the timing of the separation and the maternal stressors.

For maternal separation the pups were placed in separate clean compartments inside a temperature- and humidity-controlled terrarium, to avoid any physical distress of the pups; the mother was placed in a clean novel cage with bedding, food, and water. Maternal separation lasted for 3 h and was applied once daily from post-natal day 1 (P1) to P14; the timing was unpredictable. Maternal stress was applied to the dam at the end of the 3 h separation period and consisted of one of the following: 20 min contention in a plastic perforated tube; 10 min forced swimming stress; 10 min tail-suspension stress. Stressors were alternated pseudorandomly. Both MS and NS groups were left undisturbed

between P14 and P21 (weaning). At weaning the sex of the pups was determined (for the present cohort: 6 male and 17 female) and they were subsequently assigned to social groups of 3-4 mice per cage, composed of animals of same sex, and subjected to the same protocol (MS or NS), but from more than one dams to avoid litter effect. The sex ratio per group was, for MS: 4M/9F; and for NS: 3M/7F. A Fisher's exact test applied to these sex ratios is not statistically significant (*p* > 0.999).

## Behavioral Characterization of the Pups at Late Adolescence

We evaluated behaviors associated with depression (anhedonia, anxiety) in the separated (*n* = 13) and non-separated (*n* = 10) pups at late adolescence (P52-59). Behavioral dimensions were assessed with the Sucrose preference (anhedonia; P52) and Dark-light tests (anxiety; P59).

### Sucrose Preference

For the sucrose preference test mice were first habituated to drink from two graduated pipettes one filled with water, and the other with sucrose solution for 3 days, the side of the sucrose pipette being alternated each day. On day 4 and after an overnight (15 h) deprivation of water, the two pipettes were presented again; one was filled with water and the other with 4% sucrose. The water and sucrose solution consumed over a 3 h-period, were measured. The sucrose preference index is defined as (sucrose consumed)/(sucrose consumed + water consumed) × 100 (percentage index).

### Dark-Light

The apparatus consisted of one box divided in two compartments, an illuminated one (30 × 20 × 20 cm), which is open, and a dark one (15 × 20 × 20 cm), which is covered with a lid. A small aperture (width of 5.5 cm and height of 7 cm) allows the mouse to freely move between the dark compartment and the illuminated one. At the beginning of the experiment, the mouse was placed in the illuminated box, facing the aperture. Time spent in the lit box was measured during a 10 min period.

Next, we performed a Z-normalization. For this an individual z-score was calculated for each test and for each animal as follows: $z-score = [(individual data for observed parameter) - (mean of control group)]/(standard deviation of control group)$. For both sucrose preference and dark/light z-scores were multiplied by -1, as decreased sucrose preference and decreased time in lit compartment measure depressive/anxiety-like behaviors. For the computation of means and standard deviation of control groups, control groups were defined as NS mice; note that for control groups the mean of z-score by behavioral dimension is equal to zero. Subsequently, a global "depression-index" for each animal was calculated by averaging the z-scores of the two individual tests as previously described (Apazoglou et al., 2018).

## Blood Collection

At the end of the behavioral evaluations (P60), 0.25 ml of blood was collected from the submandibular vein and stabilized with 1.3 ml RNAlater® solution (Life Technologies, Ambion,
Austin, TX). Mice were euthanized several months later by
pentobarbital injection.

RNA Isolation

Total RNAs were purified from the blood using the Mouse
RiboPure-Blood RNA isolation kit (Invitrogen), according to
manufacturer's recommendations. After washings, total RNAs
were eluted with 0.1 mM EDTA and were subsequently submitted
to DNase treatment (DNA-free™ kit, Life Technologies,
Ambion, Austin, TX). RNA concentration was determined using
a nanodrop ND-1000 spectrophotometer (Thermo Scientific,
Waltham, MA). Three samples (one female MS/two male NS)
did not provide enough quantities and qualities and were further
excluded from the gene expression analysis.

Candidate mRNA Expression
Quantification

850 ng of total RNA was reverse transcribed using the High-
Capacity cDNA Reverse Transcription kit (Life Technologies,
Applied Biosystems, Foster City, CA). 800 ng of the resulting
cDNA were combined with a TaqMan® Fast advanced Master
Mix (Thermo Fisher Scientific) and real-time PCR reactions were
simultaneously run in triplicate in a thermocycler under the
following conditions: 2 min at 50°C, 10 min at 92°C, 45 cycles of
1 s at 95°C, and 20 s at 60°C, in custom array microfluidic cards
(Applied Biosystems, Pleasanton, CA) using a QuantStudio™
7 system and data collected using QuantStudio™ Real-Time
PCR Software v1.1 (Applied Biosystems). For each of the 16
candidate genes tested (related to the SNPs included in the
HGPS), primer-probe sets were selected using the web portal
of the manufacturer (Applied Biosystems, see Supplementary
Table 5). In addition, Rab5a was universally used as a reference
gene (Hervé et al., 2017). Raw Ct values were obtained with
manual baseline settings on the ThermoFisher cloud RQ software
(Applied Biosystems), and then the relative expression level
of each mRNA was quantified by using the 2⁻²ΔΔCt method
(Livak and Schmittgen, 2001). In this method, each candidate
gene is quantified relative to the expression of Rab5a and each
amplification is also compared to a calibrator sample (the mean
of the samples from the NS mice).

Statistics

Behavioral Evaluation

The z-scores for anhedonia, dark light, and the global depression-
like index were analyzed separately. For each parameter,
the comparison of two independent groups (MS vs. NS)
involved Student's t-test that were performed using the
STATISTICA software. The results are expressed as mean ± SEM
(standard error).

Gene Expression

For gene expression comparisons, after observing non-
homogeneous variances for each candidate gene in each
subgroup (through the [R] function *levene_test*) and absence of
normality of residuals (through the [R] function *shapiro_test*)
of a parametric model (through the [R] function *lm*), a non-
parametric (by permutation) equivalent of a two-way factorial

ANOVA was performed through the [R] function *aovp* in the
*lmPerm* library. When ANOVA effects were significant, multiple
group comparisons for each gene were performed through the
[R] function *pairwise.perm.t.test* in the *RVAideMemoire* library
to provide FDR *p*-value adjustment.

Behavior and Gene Expression Correlations

Correlograms, allowing visualization of behavior and gene
expression data into correlation matrices were implemented
through the [R] function *corrplot* in the *corrplot* library and
variables were ordered according to first principal components.
Linear regressions with 95% confidence intervals were plotted
through the [R] functions *ggplot*, *geom_point*, and *geom.smooth*
(with "lm" method) in the *ggplot2* library.

RESULTS

Selection of Genes and SNPs

**Supplementary Table 1** lists the 90 immune-related genes selected for this study and classifies them in 6 categories: “Cytokines and Cytokine Receptors” (43 genes), “Oxidative Stress Effectors” (4), “Monocytosis and Granulopoiesis” (14), “Inflammatory Signaling Pathway” (21), “Kynurenine Pathway” (5), and “Phospholipases” (3). For these 90 genes, 674 related common SNPs were found genotyped in the IMAGEN database.

The Single Effects of Immune-Related
SNPs on Brain Structure

In our initial sample of 1563 14-year-old participants, the mean bilateral hippocampal GM volume was 1.20 ± 0.109 ml and the mean mPFC GM volume was 30.1 ± 3.42 ml. Other volumetric data from the 14-year-olds (BL) and 18-year-olds (FU2) can be found in **Supplementary Table 2**. Sex, PDS, TIV and scanner type were all found to be significantly reacted with GM volume (*p* < 0.001), regardless of the region. The individual effects of the selected 674 SNPs on BL hippocampal and mPFC GM volume were assessed through linear regression analyses, controlled for sex, PDS, TIV, and scanner type. No correlation surpassed the significance threshold of *p* < 0.05 after Bonferroni correction for multiple comparisons.

Construction of the RRIs

As described in the section “Materials and Methods,” we
constructed two scores, one explaining the hippocampal GM
volume at BL (HRRI), the other explaining the mPFC GM volume
at BL (MRRI), using only immune-related SNPs. We found 26
“independent” SNPs that were considerably correlated with the
hippocampal volume and thus incorporated in the HRRI; 29
SNPs were combined in the MRRI (Supplementary Tables 3, 4).

Correlation of the RRIs With
Psychometric Data

The first psychometric measurement of interest was the
DepBand, representing depression probability at FU2 (Table 1).
Controlling for the covariates sex and CM, DepBand, could
<table><caption>**TABLE 1** | Correlations between depression probability (DAWBA score) at age 18 and the genetic scores.</caption><thead><tr><th>Independent variable</th><th>χ<sup>2</sup> (df)</th><th>B</th><th>z</th><th>p</th><th>p (corr)</th></tr></thead><tbody><tr><td rowspan="5">DepBand at FU2</td><td>HRRI</td><td>1192 (857)</td><td>0.0207</td><td>0.770</td><td>0.441</td><td>0.721</td></tr><tr><td>MRRI</td><td>1185 (857)</td><td>0.0764</td><td>2.811</td><td>0.00494</td><td>0.0296*</td></tr><tr><td>HMRRI</td><td>1186 (857)</td><td>0.0477</td><td>2.518</td><td>0.0119</td><td>0.0536</td></tr><tr><td>CM:HRRI</td><td>1180 (856)</td><td>0.122</td><td>3.57</td><td>0.000364</td><td>0.00327*</td></tr><tr><td>CM:MRRI</td><td>1184 (856)</td><td>0.0271</td><td>0.735</td><td>0.463</td><td>0.556</td></tr><tr><td></td><td>CM:HMRRI</td><td>1179 (856)</td><td>0.0679</td><td>2.80</td><td>0.00516</td><td>0.0310*</td></tr></tbody></table>

The DepBand score (DAWBA questionnaire) obtained at FU2 was correlated with the three RRIs (ROI-Related-Immune-gene-score) as well as the three interactions between the scores and childhood maltreatment (CM) by means of a Poisson regression (N = 861).

χ² (df) = residual deviance (degrees of freedom), used for Chi-Squared goodness-of-fit test; B = unstandardized regression coefficient.

* p < 0.05 (FDR-corrected).

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Correlations between RRIs (determined from their link with GM regional volumes at age 14 as described in the text) and psychometric measurements at age 18. MPGS = MRRI, HPGS = HRRI, HMPGS = HMRRI. The colored lines are proposed regression lines for the different levels of childhood maltreatment (CM), after controlling for CM and sex. Red: CM = 0; green: CM = 1; blue: CM = 2. The blue and green line overlap in the bottom left graph.</figcaption></figure>

not be correlated with the HRRI [$p(\text{corr}) = 0.721$]. The MRRI, however, positively and significantly covaried with DepBand [$p(\text{corr}) = 0.0296$; **Figure 1**]. Since both scores were not found redundant ($p = 0.436$), we also created the HMRRI by adding

up the HRRI and MRRI. The HMRRI positively covaried with DepBand [$p(\text{corr}) = 0.0536$] but did not survive FDR-correction. Next, the interactions between CM and the scores in relation to DepBand were evaluated. The number of participants with
a CM score of 0 was 588, 283 had a score of 1, and 108 had a score of 2. A positive and significant interaction was observed between the HRRI and childhood maltreatment score CM [p corr) = 0.00327]. This suggests that, as the level of endured trauma increases, the association between MRRI and presence of depressive symptoms at age 18 increases as well. It is important to note that all goodness-of-fit Chi-Squared tests for the Poisson regressions were found significant, suggesting that the data do not fit the model perfectly well. Alternative models were considered but not found better.

The second psychometric tool of interest was the CAPE-42, which consists out of three subscores (Table 2). A positive and significant correlation was observed between the log-transformed Positive Dimension Score, representing the amount of positive psychotic symptoms at FU2, and the HMRRI [p corr) = 0.0296; Figure 1]. As no significant interactions were found between CM and the RRIs in explaining one of the three CAPE-42 subscores, these results are not shown in Table 2.

Lastly, the Externalizing Score (ES) and Internalizing Score (IS) at FU2 were modeled in function of the scores, again using Poisson regression (Table 3). A significant inverse correlation was observed between the ES and the HRRI [p corr) = 0.0294; Figure 1]. We also found a significant positive interaction between HRRI and CM in explaining the ES [p corr) = 0.000857].

## Correlation of the Control Scores With Psychometric Data

No significant correlations were found between the control scores and the different psychiatric symptom measurements. In addition, we did not observe any significant interactions between childhood maltreatment and the control scores in function of the psychiatric symptoms.

## Behavioral Effects of Early Life Stress in Adolescent Animals

**Supplementary Figures 2, 3** as well as **Figure 2** show that early life stress profoundly affects anhedonia and anxiety at late adolescence in mice. At late adolescence, mice subjected to early life stress between P1 and P14 show increased anhedonia and

anxiety as compared to NS mice (raw data are presented in **Supplementary Figures 2, 3** and z-score in **Figure 2**). The global depression-index (**Figure 2**) is also significantly elevated in MS mice indicating that early life adversity (P1-P14) induces long-lasting negative affects still present in late adolescence (P52-P59).

## Effect of Early Life Stress on Candidate Gene Expression in Adolescent Animals

We next sought to investigate the effect of early life stress on the transcriptional expression of HRRI set of genes. Among the 17 candidate genes, 3 were not analyzed further (IKBKG, IL12B, IL13) because of low quality results. The remaining 14 candidate genes were well expressed in mouse blood and subsequently constituted the focus of our analysis as “mouse HRRI.” **Figure 3** shows that maternal separation has a significant effect on the global transcript level of mouse HRRI (p < 2.0E-16). Individually, post-hoc analysis demonstrated significant dysregulation for Ikbkb, Il10ra, Il10rb, Il18, Pla2g6, and Ptgs1. Among those Ptgs1 was increased while all the others were decreased in MS mice.

## Correlations Between Gene Expression and Behavior in Adolescent Mice

**Figure 4** shows that among the 14 mouse HRRI genes profiled, 6 were correlated with the global depression-index (a composite score of anhedonia and anxiety subscores). We noted that all the genes that correlated with behavior, also highly correlated with each other. **Figure 5** illustrates the tight link between one of the most significant altered genes, Il10rb, and depressive-like behavior. Regression analyses for behavior and Ikbkb, Il10rb, Il18, Pla2g6, and Ptgs1 are shown in **Supplementary Figure 4**.

## DISCUSSION

In this translational study, we explored the link between inflammation-related genes and brain structure, along with early life adversity and emergence of psychiatric symptoms. We investigated imaging genetics in a large database of community adolescents. We formed scores aggregating inflammation genes related with hippocampal or mPFC volumes at age 14. We found

<table><caption>**TABLE 2** | Correlations between the clinical dimensions (CAPE-42 subscores) and the genetic scores.</caption><thead><tr><th></th><th>Independent variable</th><th>R<sup>2</sup></th><th>B</th><th>t</th><th>p</th><th>p (corr)</th></tr></thead><tbody><tr><td rowspan="3">Positive Dimension (log)</td><td>HRRI</td><td>0.0449</td><td>0.0115</td><td>2.01</td><td>0.0444</td><td>0.114</td></tr><tr><td>MRRI</td><td>0.0466</td><td>0.0134</td><td>2.38</td><td>0.0175</td><td>0.0630</td></tr><tr><td>HMRRI</td><td>0.0506</td><td>0.0124</td><td>3.11</td><td>0.00192</td><td>0.0296*</td></tr><tr><td rowspan="3">Negative Dimension (log)</td><td>HRRI</td><td>0.0431</td><td>-0.00887</td><td>-1.06</td><td>0.289</td><td>0.520</td></tr><tr><td>MRRI</td><td>0.0420</td><td>0.0087</td><td>0.226</td><td>0.821</td><td>0.921</td></tr><tr><td>HMRRI</td><td>0.0423</td><td>-0.00342</td><td>-0.583</td><td>0.560</td><td>0.840</td></tr><tr><td rowspan="3">Depressive Dimension (log)</td><td>HRRI</td><td>0.091</td><td>0.000741</td><td>0.091</td><td>0.928</td><td>0.928</td></tr><tr><td>MRRI</td><td>0.0933</td><td>0.0124</td><td>1.54</td><td>0.123</td><td>0.277</td></tr><tr><td>HMRRI</td><td>0.0952</td><td>0.00667</td><td>1.161</td><td>0.246</td><td>0.492</td></tr></tbody></table>

The three subscores of the CAPE-42 questionnaire obtained at FU2 were correlated with the three RRIs after log-transformation of the dependent variables (N = 930).

R² = coefficient of determination; B = unstandardized regression coefficient.

* p < 0.05 (FDR-corrected).
<table><caption>**TABLE 3** | Correlations between Externalizing and Internalizing Score (SDQ) at age 18 and the genetic scores.</caption><thead><tr><th></th><th>Independent variable</th><th>χ<sup>2</sup> (df)</th><th>B</th><th>z</th><th>p</th><th>p (corr)</th></tr></thead><tbody><tr><td rowspan="5">Externalizing at FU2</td><td>HRRI</td><td>1563 (951)</td><td>-0.0275</td><td>-2.90</td><td>0.00368</td><td>0.0294*</td></tr><tr><td>MRRI</td><td>1572 (951)</td><td>-0.00262</td><td>-0.282</td><td>0.778</td><td>0.921</td></tr><tr><td>HMRRI</td><td>1567 (951)</td><td>-0.0146</td><td>-2.22</td><td>0.0263</td><td>0.079</td></tr><tr><td>CM:HRRI</td><td>1547 (950)</td><td>0.0538</td><td>4.07</td><td>4.76 × 10<sup>-5</sup></td><td>0.000857*</td></tr><tr><td>CM:MRRI</td><td>1572 (950)</td><td>-0.00832</td><td>-0.610</td><td>0.542</td><td>0.588</td></tr><tr><td rowspan="5">Internalizing at FU2</td><td>CM:HMRRI</td><td>1561 (950)</td><td>0.0228</td><td>2.44</td><td>0.01473</td><td>0.0529</td></tr><tr><td>HRRI</td><td>1824 (951)</td><td>-0.00465</td><td>-0.471</td><td>0.637</td><td>0.844</td></tr><tr><td>MRRI</td><td>1825 (951)</td><td>-0.00158</td><td>-0.164</td><td>0.870</td><td>0.921</td></tr><tr><td>HMRRI</td><td>1824 (951)</td><td>-0.00305</td><td>-0.445</td><td>0.657</td><td>0.845</td></tr><tr><td>CM:HRRI</td><td>1819 (950)</td><td>0.0298</td><td>2.21</td><td>0.0269</td><td>0.0807</td></tr><tr><td rowspan="2"></td><td>CM:MRRI</td><td>1824 (950)</td><td>-0.00118</td><td>-0.0850</td><td>0.932</td><td>0.932</td></tr><tr><td>CM:HMRRI</td><td>1822 (950)</td><td>0.0143</td><td>1.51</td><td>0.132</td><td>0.237</td></tr></tbody></table>

The Externalizing and Internalizing Score (SDQ) obtained at FU2 were correlated with the three RRIs as well as the three interactions between the RRIs and childhood maltreatment (CM) by means of a Poisson regression (N = 955).

χ²(df) = residual deviance (degrees of freedom), used for Chi-Squared goodness-of-fit test; B = unstandardized regression coefficient.

*p* < 0.05 (FDR-corrected).

that these RRIs related with psychiatric symptoms at age 18. Also, we found a set of inflammation genes that were related to gray matter volume of hippocampal regions, and to childhood maltreatment score in these adolescents. Expression levels of inflammation genes associated with psychiatric symptoms (genes from HRRI) were subsequently examined in a murine model of early life adversity. Interestingly, expression of these genes was significantly altered after maternal separation in mice.

Excessive activation of the immune system as well as abnormalities in brain structure have been associated with

depression. Furthermore, depression as well as its associated structural abnormalities are relatively heritable. Meta-analyses have found a heritability of 37% for depression (Sullivan et al., 2000) and a heritability between 40 and 70% for hippocampal volume (Gu and Kanai, 2014). Yet, characterization of the contributing individual genetic factors has proven to be very difficult. In a large genome wide association (GWA) study that investigated how common genetic variants affect the structure of subcortical regions, only a couple of genetic loci could be significantly correlated (Hibar et al., 2015). Similarly, despite considerable success within other illnesses such as diabetes and rheumatoid arthritis, GWA analyses of MDD have overall failed to produce results at the SNP level (Bogdan et al., 2017). Explanations for this could be found in the phenotypic heterogeneity, the lack of very large sample sizes and the complex functional architecture of the genetic polymorphisms. This last issue has been addressed by exploring gene-environment interactions and applying polygenic approaches. For example, in a recent study, structural abnormalities in schizophrenia were explored by creating a polygenic risk score (PRS) based on the weighted effects of SNPs found associated with schizophrenia in a prior GWA study (Alnæs et al., 2019). Also, PRSs for bipolar disorder and PRSs for schizophrenia both were found to have certain predictive power with regards to depression, corroborating prior evidence that these disorders share some common genetic overlap (Musliner et al., 2019).

Since psychiatric disorders as well as their associated structural abnormalities seem to involve a genetic contribution, we hypothesized the existence of an immune-related genetic overlap between GM structural reductions and psychiatric symptoms. This was explored by constructing two RRIs based solely on SNPs related to inflammatory genes: one predicting hippocampal GM volume in 14-year-olds (HRRI), the other predicting mPFC GM volume (MRRI). We found that both scores were correlated with the presence of different symptoms later in adolescence.

We observed that the RRIs describing the genetic variation in less than 30 inflammatory SNPs had small yet significant

**FIGURE 2** | Early life stress induces depressive-like behaviors in adolescent mice. Newborn mice were either subjected to a maternal separation paradigm between P1 (post-natal day 1) and P14 (MS, black dots) or were left undisturbed (controls, NS, white dots). NS and MS mice were evaluated for anhedonia and anxiety (measured in the sucrose preference and the dark-light tests, respectively) in late adolescence, between P52 and P59. Anhedonia and anxiety scores were z-transformed and a composite depression-index (global z-score) was averaged. Two-tailed Student's t-test shows increased levels of anhedonia (df = 21; t = 4.155) and anxiety (df = 21; t = 2.202) as well as increased depression index (df = 21, t = 4.797) in MS mice as compared to NS. *p* < 0.05; **p* < 0.0001.
<figure><img src="image_4.png" /><figcaption>**FIGURE 3** | Transcriptional expression of HRRI candidate genes in adolescent mice. Total RNA from blood of NS and MS mice at P60 was profiled by RT-qPCR for 14 HRRI candidate gene transcripts. The expression of each transcript was quantified relative to the expression of a reference gene, *Rab5a*, whereas the mean of NS mice was used as a calibrator. Statistical analysis was realized using a permutation-based non-parametric factorial ANOVA. Yellow and blue dashed lines indicate the mean values of all NS and MS, respectively, and highlight a significant “gene” effect between NS and MS animals (*p* = 0.00295). Significant multigroup comparisons for each gene were performed by pairwise permutation *t*-tests and are as follows, *Ptgs1* (*p* = 0.022), *Ikbkb* (*p* = 0.028), *Pla2g6* (*p* = 0.022), *Il10ra* (*p* = 0.008), *Il10rb* (*p* = 0.002), *Il18* (*p* = 0.032). \* *p* < 0.05; \*\* *p* < 0.01.</figcaption></figure>

predictive power regarding certain psychometric measurements obtained at age 18. As expected, effect sizes were consistently relatively small. The MRRI was found correlated with the presence of depressive symptoms. Secondly, both HRRI and MRRI were correlated with positive psychotic symptoms. Schizophrenia, the main disorder linked with these symptoms, has been consistently associated with neuroanatomical abnormalities such as reductions in GM volume (Weinberger, 1987; Keshavan et al., 2005; Bakhshi and Chance, 2015). Thirdly, we found a negative relationship between the HRRI and externalizing symptoms at age 18. This seems unexpected at first sight. A recent review addressing neuro-imaging findings in two of the major externalizing disorders, conduct disorder and oppositional defiant disorder, did not describe any studies reporting increases of hippocampal GM volume (Noordermeer et al., 2016). However, functional deficiencies in the amygdala, common in externalizing disorders, could be explained by abnormalities in the neighboring hippocampal complex (Yang and Wang, 2017). Lastly, we observed positive interactions between the scores and childhood maltreatment. This means that the probability of developing certain psychiatric symptoms due to a history of childhood maltreatment will be larger in the context of a specific genetic predisposition, in this case a high RRI. Gene-environment interactions in psychiatric diseases have been described repeatedly. For example, a polymorphism in the promoter region of the serotonin transporter gene was reported to moderate the influence of stressful life events on depression (Caspi et al., 2003).

The ability of the RRIs to predict to a limited extent the presence of psychiatric symptoms suggests the existence of the proposed genetic overlap. Indeed, the same variation in immune-related genes was found to explain both GM volume reductions in the hippocampus or mPFC and the degree of certain psychiatric symptoms. However, it could be argued that the RRIs predictive ability is solely due to the fact that the RRIs are constructed in such a way that they represent a portion of the GM volume variance. As the link between GM volumes and psychiatric illnesses is already established, the ability of an alternative score representing those structural abnormalities to predict psychiatric symptoms would not be surprising. In order to investigate this, we performed a control study by constructing two scores on the basis of random SNPs. These control scores did not significantly correlate with any psychometric measurement, nor did they display any interactions with childhood maltreatment. This higher predictive power of the non-random RRIs points at the involvement of the immune system. We thus not only corroborate prior evidence for the link between structural GM reductions and psychiatric illnesses, but also provide pioneering evidence strongly suggesting an immune-related genetic overlap between regional GM volumes and psychiatric symptoms, and define a novel combination of genes involved in this link. In order to further investigate this suggested causality, it would be interesting to perform a longitudinal study in which brain structural changes during adolescence are associated with the development of psychiatric symptoms.
<figure><img src="image_5.png" /></figure>

**FIGURE 4** | Correlation between transcriptional expression of HRRI candidate genes and depressive-like behavior in adolescent mice. Correlation matrix graph, correlogram, highlighting correlation between qualitative (Z_Behavior, a depression-index reflecting a composite score of anhedonia and anxiety subscores) and quantitative (transcriptional relative level of HRRI candidate genes obtained by RT-qPCR from blood samples) variables in adolescent mice. The variables are ordered according to first principal components. Positive correlations are displayed in red and negative correlations in blue color. The intensity of the color and the size of the squares are proportional to the Pearson correlation coefficients. Only significant correlations are indicated by a colored square.

To further investigate the functional importance of the novel gene-set that we defined we employed a translational approach linking genotype and gene expression analyses. Indeed, the phenotype during adolescence is likely to be modulated by both genotype and environment, so that genotype analyses alone probably cannot account for their interaction (Kent et al., 2012). In contrast, transcriptional profiling that measures the expression of genes is sensitive to both genotype and environment and therefore may offer insights in pathophysiology. We focused on blood transcriptomics, since blood signature demonstrated that it could represent a surrogate for brain gene expression and may predict stress-induced behaviors (Sullivan et al., 2006; van Heerden et al., 2009; Rollins et al., 2010; Tylee et al., 2013; Witt et al., 2013; Luykx et al., 2016; Hervé et al., 2017).

We implemented an animal model of early life adversity and measured depression-like behaviors in adolescence. Based on the above reported effects in human subjects we hypothesized that early life adversity would affect not only behavior but also the expression of our set of genes, and that expression would correlate with symptoms of negative affects at adolescence.

Mice were subjected to a protocol of early life adversity (maternal separation) at an early post-natal age (P1-P15). We evaluated behaviors associated with depression (anhedonia: sucrose preference; anxiety: dark-light box) in MS and NS mice at adolescence (P52-P59). We measured the expression of our genes in whole blood samples collected at the same time-point (P60). We specifically focused on the hippocampal

gene-set since the hippocampus is a region consistently implicated in depression and depression-like phenotypes in humans and mice (Vythilingam et al., 2002; Turecki et al., 2012; Apazoglou et al., 2018).

Our results showed that mice subjected to early life adversity displayed negative affects at adolescence (Figure 2). The expression of our mouse HRRI gene-set was altered in mice subjected to early life adversity (Figure 3), and transcript levels inside this gene-set correlated with depression-related behavioral score at adolescence (Figure 4).

Within the examined gene-set, **Figure 3** shows a significant decrease for *Ikbkb*, *Il10ra*, *Il10rb*, *Il18*, *Pla2g6*, and an increase for *Ptgs1* transcripts.

**Figure 4** shows that among these genes, three (*Ikbkb*, *Il10rb*, *Pla2g6*) were significantly correlated with the global depression-index in MS mice. Notably the three also highly correlated with each other.

Interestingly, *Ikbkb*, *Il10rb*, and *Pla2g6*, which were decreased in MS mice are implicated in inflammatory homeostasis. *Ikbkb* is a regulator of the canonical NF-Kappa-B pathway a key-pathway in immunity/inflammation (Schmid and Birbach, 2008); *Il10rb* encodes for the anti-inflammatory cytokine IL10 (Shouval et al., 2014), and mice invalidated for IL10 show increased depressive-like behaviors (Mesquita et al., 2008); *Pla2g6* encodes for the iPLA2β protein, which regulates an overall anti-inflammatory response and whose dysregulation is associated with neurogenerative disorders (Guo et al., 2018).
<figure><img src="image_6.png" /><figcaption>**FIGURE 5** | Link between *Il10rb* relative level of transcriptional expression and depressive-like behavior. Linear regressions with 95% confidence intervals (in gray) are plotted between the depression-index and the *Il10rb* transcriptional relative level obtained by RT-qPCR from blood in MS (red circles) and NS (green triangles) adolescent mice. The Pearson correlation coefficient and the associated *p*-value are indicated.</figcaption></figure>

The present findings suggest an inflammatory network of genes that most likely is involved in “depression-associated” neuroinflammatory adaptations in the periphery and CNS. We propose that early stressors like adversity can trigger an imbalance between anti-inflammatory and pro-inflammatory transcripts that may be at the origin of psychiatric symptoms in adolescence. These transcripts might provide both clinical biomarkers and novel targets in understanding and preventing individual developmental trajectories of psychiatric vulnerability.

## DATA AVAILABILITY STATEMENT

The data analyzed in this study is subject to the following licenses/restrictions: The datasets used for analysis can be acceded through requests to the last two authors. Requests to access these datasets should be directed to J-LM and ET, Jean-luc.martinot@inserm.fr and eleni.tzavara@inserm.fr.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by the Ethics committee CPP IDF 7, le Kremlin Bicêtre, France. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin. The animal study was reviewed and approved by the Comité d'éthique en expérimentation animale Charles Darwin N°5, Institut du Cerveau et de la Moelle Epinière, Paris.

## AUTHOR CONTRIBUTIONS

LPe, J-LM, ET, EI, and SD: study conception. J-LM, ET, RBe, EQ, and AH: study administration. IF, VG, TB, AB, HG, RBr, PG, FN, TP, MS, HW, RW, and JG: human study data acquisition. DP, AG, LPo, and JF: human study data QC. VG, EI, SF, CdM, and JG: mouse data acquisition. LPe, EI, EA, and ET: data analysis. LPe, EI, M-LP, ET, and J-LM: manuscript writing. M-LP, ET, EI, and J-LM: manuscript reviewing and editing. All authors contributed to the article and approved the submitted version.

## FUNDING

This work received support from the following sources: the European Union-funded FP6 Integrated Project IMAGEN (Reinforcement-related behaviour in normal brain function and psychopathology) (LSHM-CT-2007-037286), the Horizon 2020 funded ERC Advanced Grant “STRATIFY” (Brain network based stratification of reinforcement-related disorders) (695313), Human Brain Project (HBP SGA 2, 785907, and HBP SGA 3, 945539), the Medical Research Council Grant “c-VEDA” (Consortium on Vulnerability to Externalizing Disorders and Addictions) (MR/N000390/1), the National Institute of Health (NIH) (R01DA049238, A decentralized macro and micro gene-by-environment interaction analysis of substance use behavior and its brain biomarkers), the National Institute for Health Research (NIHR) Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King's College London, the Bundesministerium für Bildung und Forschung (BMBF grants 01GS08152; 01EV0711; Forschungsnetz AERIAL 01EE1406A, 01EE1406B; Forschungsnetz IMAC-Mind 01GL1745B), the Deutsche Forschungsgemeinschaft (DFG grants SM 80/7-2, SFB 940, TRR 265, NE 1383/14-1), the Medical Research Foundation and Medical Research Council (grants MR/R00465X/1 and MR/S020306/1), the National Institutes of Health (NIH) funded ENIGMA (grants 5U54EB020403-05 and 1R56AG058854-01). Further support was provided by grants from: – the ANR (ANR-12-SAMA-0004, AAPG2019 – GeBra, ANR-19-NEUR-ANTaRES), the Eranet Neuron (AF12-NEUR0008-01 – WM2NA; and ANR-18-NEUR00002-01 – ADORe), the Fondation de France (00081242), the Fondation pour la Recherche Médicale (DPA20140629802), the Mission Interministérielle de Lutte-contre-les-Drogues-et-les-Conduites-Addictives (MILDECA), the Assistance-Publique-Hôpitaux-de-Paris and INSERM (interface grant), Paris Sud University IDEX 2012, the Fondation de l'Avenir (grant AP-RM-17-013), the Fédération pour la Recherche sur le Cerveau; the National Institutes of Health, Science Foundation Ireland (16/ERCD/3797), United States (Axon, Testosterone and Mental Health during Adolescence; RO1 MH085772-01A1), and by NIH Consortium grant U54 EB020403, supported by a cross-NIH alliance that funds Big Data to Knowledge Centres of Excellence.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnsys.2021.725413/full#supplementary-material
REFERENCES

Alnæs, D., Kaufmann, T., van der Meer, D., Coirdova-Palomera, A., Rokicki, J., Moberget, T., et al. (2019). Brain heterogeneity in schizophrenia and its association with polygenic risk. *JAMA Psychiatry* 76, 739–748.

Apazoglou, K., Farley, S., Gorgievski, V., Belzeaux, R., Lopez, J. P., Grenier, J., et al. (2018). Antidepressive effects of targeting ELK-1 signal transduction. *Nat. Med.* 24, 591–597.

Bagot, R. C., Labontè, B., Penpa, C. J., and Nestler, E. J. (2014). Epigenetic signaling in psychiatric disorders: stress and depression. *Dialogues Clin. Neurosci.* 16, 281–295. doi: 10.31887/dcns.2014.16.3/rbagot

Bakhshi, K., and Chance, S. A. (2015). The neuropathology of schizophrenia: a selective review of past studies and emerging themes in brain structure and cytoarchitecture. *Neuroscience* 303, 82–102. doi: 10.1016/j.neuroscience.2015.06.028

Barnes, J., Mondelli, V., and Pariante, C. M. (2017). Genetic contributions of inflammation to depression. *Neuropsychopharmacology* 42, 81–98. doi: 10.1038/npp.2016.169

Belleau, E. L., Treadway, M. T., and Pizzagalli, D. A. (2019). The impact of stress and major depressive disorder on hippocampal and medial prefrontal cortex morphology. *Biol. Psychiatry* 85, 443–453. doi: 10.1016/j.biopsych.2018.09.031

Bogdan, R., Salmeron, B. J., Carey, C. E., Agrawal, A., Calhoun, V. D., Garavan, H., et al. (2017). Imaging genetics and genomics in psychiatry: a critical review of progress and potential. *Biol. Psychiatry* 82, 165–175. doi: 10.1016/j.biopsych.2016.12.030

Caspi, A., and Moffitt, T. E. (2006). Gene-environment interactions in psychiatry: joining forces with neuroscience. *Nat. Rev. Neurosci.* 7, 583–590. doi: 10.1038/nrn1925

Caspi, A., Sugden, K., Moffitt, T. E., Taylor, A., Craig, I. W., Harrington, H., et al. (2003). Influence of life stress on depression: moderation by a polymorphism in the 5-HTT gene. *Science* 301, 386–389. doi: 10.1126/science.1083968

Chen, M. C., Hamilton, J. P., and Gotlib, I. H. (2010). Decreased hippocampal volume in healthy girls at risk of depression. *Arch. Gen. Psychiatry* 67, 270–276. doi: 10.1001/archgenpsychiatry.2009.202

Choi, S. W., Mak, T. S. H., and O'Reilly, P. F. (2020). Tutorial: a guide to performing polygenic risk score analyses. *Nat. Protoc.* 15, 2759–2772. doi: 10.1038/s41596-020-0353-1

Desrivières, S., Lourdusamy, A., Tao, C., Toro, R., Jia, T., Loth, E., et al. (2015). Single nucleotide polymorphism in the neuroplastin locus associates with cortical thickness and intellectual ability in adolescents. *Mol. Psychiatry* 20, 263–274. doi: 10.1038/mp.2013.197

Ducharme, S., Albaugh, M. D., Hudziak, J. J., Botteron, K. N., Nguyen, T. V., Truong, C., et al. (2014). Anxious/depressed symptoms are linked to right ventromedial prefrontal cortical thickness maturation in healthy children and young adults. *Cereb. Cortex* 11, 2941–2950.

Franklin, T. B., Russig, H., Weiss, I. C., Gräff, J., Linder, N., Michalon, A., et al. (2010). Epigenetic transmission of the impact of early stress across generations. *Biol. Psychiatry* 68, 408–415. doi: 10.1016/j.biopsych.2010.05.036

Gu, J., and Kanai, R. (2014). What contributes to individual differences in brain structure? *Front. Hum. Neurosci.* 8:262. doi: 10.3389/fnhum.2014.00262

Guo, Y. P., Tang, B. S., and Guo, J. F. (2018). PLA2G6-associated neurodegeneration (PLAN): review of clinical phenotypes and genotypes. *Front. Neurol.* 9:1100. doi: 10.3389/fneur.2018.01100

Hervé, M., Bergon, A., Le Guisquet, A. M., Leman, S., Consoloni, J. L., Fernandez-Nunez, N., et al. (2017). Translational identification of transcriptional signatures of major depression and antidepressant response. *Front. Mol. Neurosci.* 10:248. doi: 10.3389/fnmol.2017.00248

Hibar, D. P., Stein, J. L., Renteria, M. E., Arias-Vasquez, A., Desrivières, S., Jahanshad, N., et al. (2015). Common genetic variants influence human subcortical brain structures. *Nature* 520, 224–229.

Kent, J. W. Jr., Göring, H. H., Charlesworth, J. C., Drigalenko, E., Diego, V. P., Curran, J. E., et al. (2012). Genotype × age interaction in human transcriptional ageing. *Mech. Ageing Dev.* 133, 581–590.

Keshavan, M. S., Berger, G., Zipursky, R. B., Wood, S. J., and Pantelis, C. (2005). Neurobiology of early psychosis. *Br. J. Psychiatry (Supp.)* 48, s8–s18.

Kubera, M., Obuchowicz, E., Goehler, L., Brzeszcz, J., and Maes, M. (2011). In animal models, psychosocial stress-induced (neuro)inflammation, apoptosis and reduced neurogenesis are associated to the onset of depression. *Prog. Neuropsychopharmacol. Biol. Psychiatry.* 35, 744–759. doi: 10.1016/j.pnpbp.2010.08.026

Lapidaire, W., Urrila, A. S., Artiges, E., Miranda, R., Vulser, H., Bézivin-Frere, P., et al. (2021). IMAGEN consortium. Irregular sleep habits, regional grey matter volumes, and psychological functioning in adolescents. *PLoS One* 16:e0243720. doi: 10.1371/journal.pone.0243720

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. *Methods* 25, 402–408. doi: 10.1006/meth.2001.1262

Luykx, J. J., Olde Loohuis, L. M., Neeleman, M., Strengman, E., Bakker, S. C., Lentjes, E., et al. (2016). Peripheral blood gene expression profiles linked to monoamine metabolite levels in cerebrospinal fluid. *Transl. Psychiatry* 6:e983. doi: 10.1038/tp.2016.245

MacQueen, G. M., Campbell, S., McEwen, B. S., Macdonald, K., Amano, S., Joffe, R. T., et al. (2003). Course of illness, hippocampal function, and hippocampal volume in major depression. *Proc. Natl. Acad. Sci. U.S.A.* 100, 1387–1392. doi: 10.1073/pnas.0337481100

Maes, M., Bosmans, E., Suy, E., Vandervorst, C., De Jonckheere, C., and Raus, J. (1990). Immune disturbances during major depression: upregulated expression of interleukin-2 receptors. *Neuropsychobiology* 24, 115–120. doi: 10.1159/000119472

Maes, M., Yirmiya, R., Noraberg, J., Brene, S., Hibbeln, J., Perini, G., et al. (2009). The inflammatory & neurodegenerative (I&ND) hypothesis of depression: leads for future research and new drug developments in depression. *Metab. Brain Dis.* 24, 27–53.

Mesquita, A. R., Correia-Neves, M., Roque, S., Castro, A. G., Vieira, P., Pedrosa, J., et al. (2008). IL-10 modulates depressive-like behavior. *J. Psychiatr. Res.* 43, 89–97. doi: 10.1016/j.jpsychires.2008.02.004

Miller, A., and Raison, C. (2015). The role of inflammation in depression: from evolutionary imperative to modern treatment target. *Nat. Rev. Immunol.* 16, 22–34. doi: 10.1038/nri.2015.5

Musliner, K. L., Mortensen, P. B., McGrath, J. J., Suppli, N. P., Hougaard, D. M., Bybjerg-Grauholm, J., et al. (2019). Association of polygenic liabilities for major depression, bipolar disorder, and schizophrenia with risk for depression in the Danish population. *JAMA Psychiatry* 76, 516–525. doi: 10.1001/jamapsychiatry.2018.4166

Noordermeer, S. D., Luman, M., and Oosterlaan, J. (2016). A systematic review and meta-analysis of neuroimaging in oppositional defiant disorder (ODD) and conduct disorder (CD) taking attention-deficit hyperactivity disorder (ADHD) Into Account. *Neuropsychol. Rev.* 26, 44–72. doi: 10.1007/s11065-015-9315-8

Parker, N., Patel, Y., Jackowski, A. P., Pan, P. M., Salum, G. A., Pausova, Z., et al. (2020). Assessment of neurobiological mechanisms of cortical thinning during childhood and adolescence and their implications for psychiatric disorders. *JAMA Psychiatry* 77, 1127–1136. doi: 10.1001/jamapsychiatry.2020.1495

Patel, Y., Parker, N., Shin, J., Howard, D., French, L., Thomopoulos, S. I., et al. (2021). Virtual histology of cortical thickness and shared neurobiology in 6 psychiatric disorders. *JAMA Psychiatry* 78, 47–63.

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. *Genetics* 155, 945–959. doi: 10.1093/genes/155.2.945

Raison, C. L., Capuron, L., and Miller, A. H. (2006). Cytokines sing the blues: inflammation and the pathogenesis of depression. *Trends Immunol.* 27, 24–31. doi: 10.1016/j.it.2005.11.006

Rollins, B., Martin, M. V., Morgan, L., and Vawter, M. P. (2010). Analysis of whole genome biomarker expression in blood and brain. *Am. J. Med. Genet. B Neuropsychiatr. Genet.* 153b, 919–936.

Schmaal, L., Pozzi, E., CHo, T., van Velzen, L. S., Veer, I. M., Opel, N., et al. (2020). ENIGMA MDD: seven years of global neuroimaging studies of major depression through worldwide data sharing. *Transl. Psychiatry* 10:172.

Schmaal, L., Veltman, D. J., van Erp, T. G., Samann, P. G., Frodl, T., Jahanshad, N., et al. (2016). Subcortical brain alterations in major depressive disorder: findings from the ENIGMA major depressive disorder working group. *Mol. Psychiatry.* 21, 806–812.

Schmid, J. A., and Birbach, A. (2008). IkappaB kinase beta (IKKbeta/IKK2/IKBKB)-a key molecule in signaling to the transcription
factor NF-kappaB. *Cytokine Growth Factor Rev.* 19, 157-165. doi: 10.1016/j.cytogfr.2008.01.006

Schumann, G., Loth, E., Banaschewski, T., Barbot, A., Barker, G., Buüchel, C., et al. (2010). The IMAGEN study: reinforcement-related behaviour in normal brain function and psychopathology. *Mol. Psychiatry* 15, 1128–1139. doi: 10.1038/ mp.2010.4

Shouval, D. S., Biswas, A., Goettel, J. A., Mccann, K., Conaway, E., Redhu, N. S., et al. (2014). Interleukin-10 receptor signaling in innate immune cells regulates mucosal immune tolerance and anti-inflammatory macrophage function. *Immunity* 40, 706–719. doi: 10.1016/j.immuni.2014.03.011

Sullivan, P. F., Fan, C., and Perou, C. M. (2006). Evaluating the comparability of gene expression in blood and brain. *Am. J. Med. Genet. B Neuropsychiatr. Genet.* 141b, 261–268. doi: 10.1002/ajmg.b.30272

Sullivan, P. F., Neale, M. C., and Kendler, K. S. (2000). Genetic epidemiology of major depression: review and meta-analysis. *Am. J. Psychiatry* 157, 1552–1562. doi: 10.1176/appi.ajp.157.10.1552

Turecki, G., Ernst, C., Jollant, F., Labonté, B., and Mechawar, N. (2012). The neurodevelopmental origins of suicidal behavior. *Trends Neurosci.* 35, 14–23. doi: 10.1016/j.tins.2011.11.008

Tylee, D. S., Kawaguchi, D. M., and Glatt, S. J. (2013). On the outside, looking in: a review and evaluation of the comparability of blood and brain “-omes”. *Am. J. Med. Genet. B Neuropsychiatr. Genet.* 162b, 595–603. doi: 10.1002/ajmg.b.32150

van Heerden, J. H., Conesa, A., Stein, D. J., Montaner, D., Russell, V., and Illing, N. (2009). Parallel changes in gene expression in peripheral blood mononuclear cells and the brain after maternal separation in the mouse. *BMC Res. Notes* 2:195. doi: 10.1186/1756-0500-2-195

Vythilingam, M., Heim, C., Newport, J., Miller, A. H., Anderson, E., Bronen, R., et al. (2002). Childhood trauma associated with smaller hippocampal volume in women with major depression. *Am. J. Psychiatry* 159, 2072–2080. doi: 10.1176/appi.ajp.159.12.2072

Weinberger, D. R. (1987). Implications of normal brain development for the pathogenesis of schizophrenia. *Arch. Gen. Psychiatry* 44, 660–669. doi: 10.1001/ archpsyc.1987.01800190080012

Witt, S. H., Sommer, W. H., Hansson, A. C., Sticht, C., Rietschel, M., and Witt, C. C. (2013). Comparison of gene expression profiles in the blood, hippocampus and prefrontal cortex of rats. *Silico Pharmacol.* 1:15.

Yang, Y., and Wang, J. Z. (2017). From structure to behavior in basolateral amygdala-hippocampus circuits. *Front. Neural Circuits* 11:86. doi: 10.3389/ fncir.2017.00086

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Penninck, Ibrahim, Artiges, Gorgievski, Desrivières, Farley, Filippi, de Macedo, Belzeaux, Banaschewski, Bokde, Quinlan, Flor, Grigis, Garavan, Gowland, Heinz, Brühl, Nees, Papadopoulos Orfanos, Paus, Poustka, Fröhner, Smolka, Walter, Whelan, Grenier, Schumann, Paillère Martinot, Tzavara and Martinot. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.