DOI: 10.5281/zenodo.7948162
title: Changes in patterns of age-related network connectivity are associated with risk for schizophrenia

description: This is the online data repository accompanying the following manuscript: <strong>Changes in patterns of age-related network connectivity are associated with risk for schizophrenia</strong> Roberta Passiatore<sup>a,b,c</sup>, Linda A. Antonucci<sup>a</sup>, Thomas P. DeRamus<sup>b</sup>, Leonardo Fazio<sup>d</sup>, Giuseppe Stolfa<sup>a</sup>, Leonardo Sportelli<sup>a,e</sup>, Gianluca C. Kikidis<sup>a,e</sup>, Giuseppe Blasi<sup>a,f</sup>, Qiang Chen<sup>e</sup>, Juergen Dukart<sup>c,g</sup>, Aaron L. Goldman<sup>e</sup>, Venkata S. Mattay<sup>e,h</sup>, Teresa Popolizio<sup>i</sup>, Antonio Rampino<sup>a,f</sup>, Fabio Sambataro<sup>j</sup>, Pierluigi Selvaggi<sup>a,f</sup>, William Ulrich<sup>e</sup>, Apulian Network on Risk for Psychosis<sup>a,k,l,m,n,o</sup>, Daniel R. Weinberger<sup>e,h,p,q,r</sup>, Alessandro Bertolino<sup>a,f1</sup>, Vince D. Calhoun<sup>b1</sup>, Giulio Pergola<sup>a,e,p1</sup> <sup>a</sup> Department of Translational Biomedicine and Neuroscience, University of Bari Aldo Moro, 70124 Bari, Italy <sup>b</sup> Tri-institutional Center for Translational Research in Neuroimaging and Data Science, Georgia State University, Georgia Institute of Technology, and Emory University, 30303 Atlanta, GA <sup>c</sup> Institute of Neuroscience and Medicine, Brain and Behavior, Research Centre Jülich, 52428 Jülich, Germany <sup>d</sup> Department of Medicine and Surgery, Libera Università Mediterranea Giuseppe Degennaro, 70010 Casamassima, Italy <sup>e </sup>Lieber Institute for Brain Development, Johns Hopkins Medical Campus, 21205 Baltimore, MD <sup>f</sup> Psychiatric Unit, University Hospital, 70124 Bari, Italy <sup>g</sup> Institute of Systems Neuroscience, Medical Faculty, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany <sup>h</sup> Department of Neurology and Radiology, Johns Hopkins Medical Campus, 21287 Baltimore, MD <sup>i</sup> Neuroradiology Unit, Scientific Institute for Research, Hospitalization and Health Care, Casa Sollievo della Sofferenza, 71013 San Giovanni Rotondo, Foggia, Italy <sup>j</sup> Section of Psychiatry, Department of Neuroscience, University of Padova, 35121 Padova, Italy <sup>k</sup> Department of Mental Health, Azienda Sanitaria Locale Foggia, 71121 Foggia, Italy <sup>l </sup>Department of Clinical and Experimental Medicine, University of Foggia, 71122 Foggia, Italy <sup>m </sup>Department of Mental Health, Azienda Sanitaria Locale Barletta-Andria-Trani, 76123 Andria, Italy <sup>n </sup>Department of Mental Health, Azienda Sanitaria Locale Bari, 70132 Bari, Italy <sup>o </sup>Department of Mental Health, Azienda Sanitaria Locale Brindisi, 72100 Brindisi, Italy <sup>p</sup> Department of Psychiatry and Behavioral Sciences, Johns Hopkins University School of Medicine, 21205 Baltimore, MD <sup>q</sup> Department of Neuroscience, Johns Hopkins University School of Medicine, 21287 Baltimore, MD<br> <sup>r</sup> Department of Genetic Medicine, Johns Hopkins University School of Medicine, 21287 Baltimore, MD Members of the Apulian Network on Risk for Psychosis include Ileana Andriola, Lucia Mare, Nicolò Parente, Alessandra Raio, Veronica D. Toro, (Department of Translational Biomedicine and Neuroscience – University of Bari Aldo Moro, Bari, IT), Mario Altamura, Marimar Castrigno, Melania Difino (Department of Mental Health, ASL Foggia, Foggia, IT; Department of Clinical and Experimental Medicine, University of Foggia, Foggia, IT), Flora Brudaglio, Domenico S. Savino, Rossana Vista (Department of Mental Health, ASL Barletta-Andria-Trani, Andria, IT), Angela Carofiglio, Barbara Gelao, Marina Mancini (Department of Mental Health, ASL Bari, Bari, IT), Alessandro Saponaro, Anna Manzari, Domenico Suma (Department of Mental Health, ASL Brindisi, Brindisi, IT). <strong>Abstract</strong> Alterations in fMRI-based brain functional network connectivity (FNC) are associated with schizophrenia (SCZ) and the genetic risk or subthreshold clinical symptoms preceding the onset of SCZ, which often occurs in early adulthood. Thus, age-sensitive FNC changes may be relevant to SCZ risk-related FNC. We used independent component analysis to estimate FNC from childhood to adulthood in 9,236 individuals. To capture individual brain features more accurately than single- session fMRI, we studied an average of three fMRI scans per individual. To identify potential familial risk-related FNC changes, we compared age-related FNC in first-degree relatives of SCZ patients mostly including unaffected siblings (SIB) with neurotypical controls (NC) at the same age-stage. Then, we examined how polygenic risk scores for SCZ influenced risk-related FNC patterns. Finally, we investigated the same risk-related FNC patterns in adult SCZ patients (oSCZ) and young individuals with subclinical psychotic symptoms (PSY). Age-sensitive risk-related FNC patterns emerge during adolescence and early adulthood, but not before. Young SIB always followed older NC patterns, with decreased FNC in a cerebellar-occipitoparietal circuit and increased FNC in two prefrontal-sensorimotor circuits when compared to young NC. Two of these FNC alterations were also found in oSCZ, with one exhibiting reversed pattern. All were linked to polygenic risk for SCZ in unrelated individuals (R2 varied from 0.02 to 0.05). Young PSY showed FNC alterations in the same direction as SIB when compared to NC. These results suggest that age-related neurotypical FNC correlate with genetic risk for SCZ and are detectable with MRI in young participants. <em><strong>DOI: 10.1073/pnas.2221533120</strong></em> Files: - <strong>NeuroMark_1.0_template </strong> <strong><em>NeuroMark_1.0.nii</em>: </strong> 4D Nifti containing the 100 Independent Components (ICs) spatial maps that was used for individual IC extraction based on Du Y, Fu Z, Sui J, Gao S, Xing Y, Lin D, Salman M, Abrol A, Rahaman MA, Chen J, Hong LE, Kochunov P, Osuch EA, Calhoun VD; Alzheimer's Disease Neuroimaging Initiative. NeuroMark: An automated and adaptive ICA based pipeline to identify reproducible fMRI markers of brain disorders. Neuroimage Clin. 2020;28:102375. doi: 10.1016/j.nicl.2020.102375. Epub 2020 Aug 11. PMID: 32961402; PMCID: PMC7509081. <em><strong>NeuroMark_1.0_Label.doc</strong></em>: document containing labels for each IC extracted including brain regions, peak coordinates of ICs, corresponding functional network in which ICs are arranged, and the corresponding image number in the NeuroMark_1.0.nii file. NOTE - <strong>AveragedMatrix</strong> <em><strong>vectorized_FNC_LIBD.RData, </strong></em><strong><em>vectorized_FNC_UNIBA1.RData, </em></strong><strong><em>vectorized_FNC_UNIBA1.RData:</em></strong> <strong><em>FNC_all_averaged (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on all subjects averaged across groups, sessions and cohorts <strong><em>FNC_all_rest (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on all subjects during resting state fMRI averaged across groups and cohorts <strong><em>FNC_all_enc (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on all subjects during encoding (RISE/PEAR) fMRI averaged across groups and cohorts <strong><em>FNC_all_ret (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on all subjects during retrieval (RISE/PEAR) fMRI averaged across groups and cohorts <strong><em>FNC_all_nback (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on all subjects during working memory (2-Back) fMRI averaged across groups and cohorts <strong><em>FNC_all_emo (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix across subjects during emotion recognition (Faces/FMT) fMRI averaged across groups and cohorts <strong><em>FNC_yNC_averaged (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on yNC averaged across sessions <strong><em>FNC_yNC_rest (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on yNC during resting state fMRI <strong><em>FNC_yNC_enc (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on yNC during encoding (RISE/PEAR) fMRI <strong><em>FNC_yNC_ret (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on yNC during retrieval (RISE/PEAR) fMRI <strong><em>FNC_yNC_nback (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on yNC during working memory (2-Back) fMRI <strong><em>FNC_yNC_emo (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on yNC during emotion recognition (Faces/FMT) fMRI <strong><em>FN_oNC_averaged (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oNC averaged across sessions <strong><em>FNC_oNC_rest (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oNC during resting state fMRI <strong><em>FNC_oNC_enc (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oNC during encoding (RISE/PEAR) fMRI <strong><em>FNC_oNC_ret (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oNC during retrieval (RISE/PEAR) fMRI <strong><em>FNC_oNC_nback (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oNC during working memory (2-Back) fMRI <strong><em>FNC_oNC_emo (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oNC during emotion recognition (Faces/FMT) fMRI <strong><em>FNC_ySIB_averaged (UNIBA1, LIBD): </em></strong>mean FNC matrix on ySIB averaged across sessions <strong><em>FNC_ySIB_rest (LIBD): </em></strong>mean FNC matrix on ySIB during resting state fMRI <strong><em>FNC_ySIB_enc (UNIBA1, LIBD): </em></strong>mean FNC matrix on ySIB during encoding (PEAR) fMRI <strong><em>FNC_ySIB_ret (LIBD): </em></strong>mean FNC matrix on ySIB during retrieval (PEAR) fMRI <strong><em>FNC_ySIB_nback (UNIBA1, LIBD): </em></strong>mean FNC matrix on ySIB during working memory (2-Back) fMRI <strong><em>FNC_ySIB_emo (UNIBA1, LIBD): </em></strong>mean FNC matrix on ySIB during emotion recognition (Faces/FMT) fMRI <strong><em>FNC_oSIB_averaged (UNIBA1, LIBD): </em></strong>mean FNC matrix on oSIB averaged across sessions <strong><em>FNC_oSIB_rest (LIBD): </em></strong>mean FNC matrix on oSIB during resting state fMRI <strong><em>FNC_oSIB_enc (UNIBA1, LIBD): </em></strong>mean FNC matrix on oSIB during encoding (PEAR) fMRI <strong><em>FNC_oSIB_ret (LIBD): </em></strong>mean FNC matrix on oSIB during retrieval (PEAR) fMRI <strong><em>FNC_oSIB_nback (UNIBA1, LIBD): </em></strong>mean FNC matrix on oSIB during working memory (2-Back) fMRI <strong><em>FNC_oSIB_emo (LIBD): </em></strong>mean FNC matrix on oSIB during emotion recognition (FMT) fMRI <strong><em>FNC_yPSY_averaged (UNIBA2): </em></strong>mean FNC matrix on yPSY averaged across sessions <strong><em>FNC_yPSY_rest (UNIBA2): </em></strong>mean FNC matrix on yPSY during resting state fMRI <strong><em>FNC_yPSY_enc (UNIBA2): </em></strong>mean FNC matrix on yPSY during encoding (RISE) fMRI <strong><em>FNC_yPSY_ret (UNIBA2): </em></strong>mean FNC matrix on yPSY during retrieval (RISE) fMRI <strong><em>FNC_yPSY_nback (UNIBA2): </em></strong>mean FNC matrix on yPSY during working memory (2-Back) fMRI <strong><em>FNC_yPSY_emo (UNIBA2): </em></strong>mean FNC matrix on yPSY during emotion recognition (Faces) fMRI <strong><em>FNC_oSCZ_averaged (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oSCZ averaged across sessions <strong><em>FNC_oSCZ_rest (UNIBA2, LIBD): </em></strong>mean FNC matrix on oSCZ during resting state fMRI <strong><em>FNC_oSCZ_enc (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oSCZ during encoding (PEAR/RISE) fMRI <strong><em>FNC_oSCZ_ret (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oSCZ during retrieval (PEAR/RISE) fMRI <strong><em>FNC_oSCZ_nback (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oSCZ during working memory (2-Back) fMRI <strong><em>FNC_oSCZ_emo (UNIBA1, UNIBA2, LIBD): </em></strong>mean FNC matrix on oSCZ during emotion recognition (Faces/FMT) fMRI <strong>PlotMatrixFNC.m (Script): </strong>code used to generate FNC plots. Example of FNC plots are included in Figure 1. - <strong>Connectogram</strong> <em><strong>MeanMatrix.mat: </strong></em>matrix representing the difference (delta) in functional network connectivity (FNC) along the 54 ICs between oNC and yNC. This matrix was used to generate Figure 2. Matlab object. <em><strong>MeanMatrix.txt: </strong></em>matrix representing the difference (delta) in functional network connectivity (FNC) along the 54 ICs between oNC and yNC. This matrix was used to generate Figure 2. Text file. <em><strong>ConnectogramNeuroMark_1.0.m (Script)</strong></em>: code used to generate Figure 2. Graphical and statistical options are provided. Please refer to https://trendscenter.wpengine.com/trends/software/gift/docs/v4.0b_gica_manual.pdf for further details. - <strong>ResFNCAgeDifferences</strong> <em><strong>data_example.RData (Data): </strong></em>example data frame for <em>AgeDiff </em>code. <strong><em>AgeDiff_real.R (Script)</em></strong>: code to test age differences through Linear mixed effect model on multiple sessions across age-groups. <em><strong>AgeDiff_permutation.R (Script): </strong></em>code to permute age differences through Linear mixed effect model on multiple sessions across age-groups. <strong><em>ResFNCAgeDifference.RData: </em></strong>summary statistics obtained on UNIBA1, UNIBA2, LIBD cohorts included in the manuscript. <em><strong>boxplot_permutation_p_LIBD.pdf, boxplot_permutation_p_UNIBA1.pdf, boxplot_permutation_p_UNIBA2.pdf: </strong></em>error bars showing results obtained after permutation (10,000 iteration) across IC pairs. <strong>- dataFNC</strong> <strong><em>data_final.RData: </em></strong>de-identified individuals' aggregate FNC data for each cohort (UNIBA1, UNIBA2, LIBD, PNC, ABCD, UKB) <em><strong>plotData.R (Script): </strong></em>code used to plot FNC differences across groups (Figure 3). Code used to assess group differences - Wilcoxon Rank sum test - is included. <strong>- ResFNCRiskDifferences</strong> <em><strong>res_cerebellar_occipitoparietal_ICs.csv: </strong></em>summary statistics of differences across groups assessed through the Wilcoxon Rank-sum test (p.adjusted<sub>FDR</sub>&lt;0.05) on the cerebellar-occipitoparietal FNC <strong><em>res_dorsolateralPFC_sensorimotor_ICs.csv: </em></strong>summary statistics of differences across groups assessed through the Wilcoxon Rank-sum test (p.adjusted<sub>FDR</sub>&lt;0.05) on medial PFC-sensorimotor FNC <strong><em>res_medialPFC_sensorimotor_ICs.csv: </em></strong>summary statistics of differences across groups assessed through the Wilcoxon Rank-sum test (p.adjusted<sub>FDR</sub>&lt;0.05) on dorsolateral PFC-sensorimotor FNC - <strong>ResFNCxPRS</strong> <strong><em>ResMetanalysisPRS_FNC_R.csv: </em></strong>summary statistics of metanalysis on the averaged FNC x PRS association across PNC, UNIBA1, UNIBA2, LIBD, UKB cohorts on younger and older groups <strong><em>ResMetanalysisPRS_FNC_children_R.csv: </em></strong>summary statistics of metanalysis on the averaged FNC x PRS association across PNC and ABCD cohorts on children. <strong><em>t2r.R (Script):</em></strong> code for t to R<sup>2</sup> transform <em><strong>MetanalysisAveragedFNC_PRS.R (Script): </strong></em>code used perform metanalysis on the averaged FNC x PRS association across PNC, UNIBA1, UNIBA2, LIBD, UKB cohorts on younger and older groups. <strong><em>MetanalysisAveragedFNC_PRS_children.R (Script): </em></strong>code used perform metanalysis on the averaged FNC x PRS association across PNC and ABCD cohorts on children. <strong><em>MetanalysisPRS_FNC_sessions_R.R (Script): </em></strong>code used perform metanalysis on the session-specific FNC x PRS association across PNC, UNIBA1, UNIBA2, LIBD, UKB cohorts on younger and older groups. The present work includes data from UK Biobank (ID #41655), Adolescence Brain Cognitive Development (ABCD; ID #12036) Study, and Philadelphia Neurodevelopmental Cohort (PNC; ID #20998). Data are publicly accessible under approved request at the following links: - UK Biobank: https://www.ukbiobank.ac.uk/ - ABCD: https://abcdstudy.org/ - PNC: https://www.med.upenn.edu/bbl/philadelphianeurodevelopmentalcohort.html The ABCD data used in this report came from NIMH Data Archive DOI: https://doi.org/10.15154/1528793. Additional code for MRI data processing and Independent Component Analysis are publicly available at the following link: https://trendscenter.org/software/ Please, check https://doi.org/10.5281/zenodo.7853030 for future updates. For any data inquiries please contact: <strong>Roberta Passiatore: </strong>Roberta.Passiatore@uniba.it <strong>Giulio Pergola: </strong>Giulio.Pergola@uniba.it - Giulio.Pergola@libd.org <strong>Vince D Calhoun:</strong> VCalhoun@gsu.edu <strong>Alessandro Bertolino:</strong> Alessandro.Bertolino@uniba.it <strong>Funding</strong> This work received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 798181 awarded to G.P. (PI) and A.B., D.R.W.; from the European Union funding within the MUR PNRR Extended Partnership initiative on Neuroscience and Neuropharmacology (Project no. PE00000006 CUP H93C22000660006 “MNESYS, A multiscale integrated approach to the study of the nervous system in health and disease”) to G.B., A.R., A.B., and G.P.; from the funding initiative Horizon Europe Seeds 2021 (Next Generation 8 EU–MUR D.M. 737/2021) for the project S68 CUP: H99J21017550006 to G.P. (PI) and L.A.A., G.B., A.R., A.B.; from the Apulian regional government for the project: “Early Identification of Psychosis Risk” to A.B.; from the Research Projects of National Relevance (PRIN) 2017 Prot. 2017K2NEF4 awarded to G.P.; from a Collaboration Grant from Ex