Striatal dopamine gene network moderates the effect of early adversity on the risk for adult psychiatric and cardiometabolic comorbidity

Barbara Barth, Danusa Mar Arcego, Euclides José de Mendonça Filho, Randriely Merscher Sobreira de Lima, Carine Parent, Carla Dalmaz, André Krumel Portella, Irina Pokhvisneva, Michael J Meaney, Patricia Pelufo Silveira
doi: https://doi.org/10.1101/2022.04.23.22274209

Abstract

Cardiometabolic and psychiatric disorders often co-exist and share common early life risk factors, such as low birth weight. However, the biological pathways linking early adversity to adult cardiometabolic/psychiatric comorbidity remain unknown. Dopamine (DA) neurotransmission in the striatum is sensitive to early adversity and influences the development of both cardiometabolic and psychiatric diseases. Here we show that a co-expression based polygenic score (ePGS) reflecting individual variations in the expression of the striatal dopamine transporter gene (SLC6A3) network significantly interacts with birth weight to predict psychiatric and cardiometabolic comorbidities in both adults (UK Biobank, N= 225,972) and adolescents (ALSPAC, N= 1188). Decreased birth weight is associated with an increased risk for psychiatric and cardiometabolic comorbidities, but the effect is dependent on a striatal SLC6A3 ePGS, that reflects individual variation in gene expression of genes coexpressed with the SLC6A3 gene in the striatum. Neuroanatomical analyses revealed that SNPs from the striatum SLC6A3 ePGS were significantly associated with prefrontal cortex gray matter density, suggesting a neuroanatomical basis for the link between early adversity and psychiatric and cardiometabolic comorbidity. Our study reveals that psychiatric and cardiometabolic diseases share common developmental pathways and underlying neurobiological mechanisms that includes dopamine signaling in the prefrontal cortex.

Introduction

The co-occurrence of more than one chronic disease1 has high prevalence in primary care settings2, inflating health care utilization and functional disability3. Psychiatric and cardiometabolic disorders, which are highly comorbid4,5, rank amongst the leading global causes of disability-adjusted life years worldwide6,7. Prospective studies show a bi-directional relationship between psychiatric and cardiometabolic conditions8. Meta-analytic evidence from longitudinal studies indicates that diabetes increases the risk for depression by approximately 25% and that depression increases the risk for type 2 diabetes by 40-60%7,9. The odds for depression also increase with one or more non-psychiatric coexisting chronic conditions, especially coronary artery disease, chronic arthritis, and stroke10. Anxiety is also associated with 41% increased risk of developing cardiovascular disease11. Among adult patients with schizophrenia, the prevalence of diabetes averages of 15%, which is higher than the 10% prevalence of diabetes in the general population. This association persists even after controlling for obesity and use of antipsychotic drugs12.

The underlying mechanism for these comorbidities remains unknown, but an emerging explanation is that psychiatric and cardiometabolic disorders share common developmental pathways. For example, low birth weight broadly reflects an unhealthy fetal environment and is considered a prevalent form of early life adversity13 associated with increased morbimortality throughout the life course14 and is specifically linked to cardiometabolic14–18 and psychiatric disorders19–24. An obvious question concerns the biological mechanisms that underlie such a developmental trajectory involved in the development of cardiometabolic and psychiatric comorbidities.

The brain dopamine system is highly sensitive to early adversity25 and proposed as a mechanism underlying developmental pathways to multiple psychiatric and metabolic comorbidities26,27. Early life adversity such as fetal growth restriction that leads to low birth weight alters dopaminergic signaling 28–30. Dysfunction of dopamine neurotransmission in both the ventral and dorsal striatum associates with depression31, as well as dysregulated food intake and altered energy homeostasis29,30,32. Striatal dopamine signaling also appears to regulate systemic glucose metabolism in humans.33 The striatum harbors dopaminergic neurons34 and the striatal dopamine transporter (DAT) is a critical regulator of striatal dopamine release and reuptake.35 Dopamine signaling is influenced by core metabolic hormones such as leptin and insulin, through their actions on the expression and function of DAT 32,36, which is encoded by the SLC6A3 (solute carrier family 6 member 3) gene.

Based on the large evidence supporting the relation between metabolism, mental health and striatal dopaminergic neurotransmission, as well as the effects of early adversity on striatal dopamine function, we hypothesized that the striatal SLC6A3 gene network underlies the association between early life adversity and the comorbidity between psychiatric and cardiometabolic disorders in humans. We therefore aimed to test if individual differences in the function of a striatal SLC6A3 gene network might moderate the effects of early life adversity on psychiatric and cardiometabolic comorbidities in adults and adolescents. To achieve this, we created a SLC6A3 striatal co-expression-based polygenic score (striatum SLC6A3 ePGS) reflecting the genetic capacity for expression of the striatal DAT1 gene network (possibly influencing dopamine signaling) and analyzed the effect of its interaction with birth weight on the comorbidity of psychiatric and cardiometabolic conditions in adults (UK Biobank) and adolescents (Avon Longitudinal Study of Parents and Children, ALSPAC).

Methods

Participants

We used genomic and phenotypic data from two cohorts, one from adults (Uk Biobank), and one from adolescents (Avon Longitudinal Study on Parents and Children, ALSPAC). Adult cohort: The UK Biobank is a large population-based study from the United Kingdom37. Participants, aged 37-73, were recruited between 2006 and 2010 resulting in 502,543 subjects. Detailed description of the inclusion/exclusion criteria for the current analysis and the corresponding sample size at each step can be found in Supplementary information, Supplementary Figure 1. After all exclusions and inclusion criteria, the number of subjects that remained for the analysis was 225,972 (mean age = 55.22, SD = 8.08) (Table 1). We used all the data available for the brain imaging analysis considering the inclusion/exclusion criteria (Supplementary information, Supplementary Figure 1, N=11,167, mean age = 53.86, SD = 7.39).

View this table:
Table 1. Description of the baseline characteristics in UK Biobank sample for high and low striatum SLC6A3 ePGS groups, which were defined by median split of the ePGS score.

Adolescent cohort

To explore our findings in an earlier developmental time point we used data from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort38–40. This is a transgenerational prospective observational cohort that recruited 14, 541 pregnant women residents in Avon County, UK. Additional recruitment (N=913) was performed later during Phases II, III and IV respectively, bringing the total sample size of prospective mother-child dyads to 15,658. For more information on ALSPAC variables, please see http://www.bristol.ac.uk/alspac/researchers/our-data/. Data from the adolescent offspring aged between 15.5 and 17.5 were used in this study. Only subjects with available phenotypic data of interest, early life adversity measure, in this case birth weight and genotyping data were considered for the analyses (N=1,188) (Table 2).

View this table:
Table 2. Description of the baseline characteristics of the ALSPAC sample used according to high and low striatum SLC6A3 ePGS. Groups were divided by the median of the ePGS score.

Detailed description of the inclusion/exclusion criteria and the corresponding sample size at each step can be found in Supplementary information, Supplementary Figure 2. See Supplementary Information Supplementary Methods for detailed description of the genotyping procedure for each cohort.

Ethics approval and consent

UK Biobank

Informed consent was obtained from each participant, and the project has been approved by the North-West Multicentre Research 580 Ethics Committee (REC reference 11/NW/0382), the National Information Governance Board for Health and Social Care, and the Community Health Index Advisory Group for UK Biobank.

Consenting participants provided baseline information, answered questions, had measurements and biological samples collected. This research has been conducted using the UK Biobank Resource under application number 41975.

ALSPAC

Participants provided informed written consent to participate in the study. Ethics approval for the study was obtained from the ALSPAC Ethics and Law Committee and the local research ethics committees (a full list of the ethics committees that approved different aspects of the ALSPAC studies is available at http://www.bristol.ac.uk/alspac/researchers/research-ethics/). Consent for biological samples has been collected in accordance with the Human Tissue Act (2004).

Consent for publication was obtained from UK Biobank and ALSPAC management teams. The use of these datasets was locally approved by the Centre intégré universitaire de santé et de services sociaux de l’Ouest-de-l’Île-de-Montréal Research Ethics Board under application number IUSMD-21-73.

Identification of the striatal SLC6A3 co-expression gene network and ePGS calculation

Figure 1A shows the steps involved in the identification of the gene co-expression networks and the calculation of the ePGS score. The ePGS was calculated considering genes co-expressed with the SLC6A3 gene in the striatum. As described previously41–48, we began by using brain region-specific RNA sequencing data from mice available at GeneNetwork (http://genenetwork.org/, HBP Rosen Striatum M430V2 (Apr05) RMA Clean)49 to identify SLC6A3 co-expressed genes (absolute value of co-expression correlation with SLC6A3 gene greater or equal to r= 0.5). GeneNetwork was used to obtain gene expression from rodents since our previous findings demonstrated multiple effects of early life adversities, especially poor fetal growth, on dopaminergic mesocorticolimbic system in rodents28–30,50–54. We then converted SLC6A3 co-expressed genes to human orthologs by using the biomaRt package55. Since we were interested in gene networks that were active during the early developmental period in which adversity occurred and when the brain is still undergoing core maturational processes in humans, we used BrainSpan to select autosomal transcripts expressed at least 1.5-fold more during fetal and childhood periods (0–60 months after birth) in comparison to adulthood (20–40 years of age). This process resulted in a list of striatal SLC6A3 co-expressed genes. We then mapped all the existent SNPs in the human ortholog genes comprising the striatum SLC6A3 gene network using biomaRt package55 in R and gathered all gene-SNP pairs from the GTEx dataset in human striatum. These lists were merged with the genotyping data from UK Biobank and ALSPAC cohorts, respectively, retaining only common SNPs and subjecting the final SNP lists to linkage disequilibrium clumping (r2 < 0.2) within 500kb radious to eliminate redundant SNPs. The process resulted in 1532 independent functional SNPs retained in UK Biobank and 1663 SNPs in ALSPAC. The final score included 67 genes in our discovery sample (UKB) (Supplementary information, Supplementary Table 1).

Figure 1.
Figure 1. Construction and characterization of the striatum SLC6A3 gene network.

A, calculation of the expression based polygenic score (ePGS) from the genes co-expressed with the SLC6A3 gene in striatum. GeneNetwork was used to generate a list of genes co-expressed with SLC6A3 in striatum in mice, which were then converted to human orthologs. BrainSpan was used to identify genes overexpressed within striatum in fetal samples and up to 5 years of age in comparison to adult samples. All SNPs from these genes, common between the study sample and GTEx databases, were retained and included in the final list of SNPs. This final list was subjected to linkage disequilibrium clumping, with removal of highly correlated SNPs. Next, for each SNP, a number of alleles at a given SNP (rs1, rs2…) was multiplied by the estimated effect of the genotype-gene expression association from GTEx. The sum of these values over all SNPs provides the striatum SLC6A3 ePGS. B, striatum SLC6A3 ePGS co-expression gene network. Co-expression pattern was mined from GeneMANIA61. The color of the node border represents the correlation sign with the SLC6A3 gene according to GeneNetwork co-expression matrix (dark purple represents negative and light purple positive correlation). Node color intensity represents betweenness (number of times a node acts as a bridge between nodes). Node border width represents the number of connections a node has with other nodes (total degree). C, co-expression of genes included in the striatum SLC6A3 gene network in humans at different ages according to BrainSpan. D, topological properties of the striatum SLC6A3 gene network, showing hubs (with degrees higher than +1SD above mean), bottlenecks (betweenness higher than +1SD above the mean), and hub-bottlenecks. Lines in black indicate mean + 1 SD for degrees and betweenness. Hub and hub-bottleneck genes are related to ribosomal structure. Among the bottleneck genes, HNRNPA1 is involved in the packaging of pre-mRNA into particles and transport from the nucleus to the cytoplasm, as well as splicing. The SDC3 gene may play a role in cell shape organization and has been associated with obesity74.

To calculate the striatal SLC6A3 ePGS, the number of effect alleles (genotype information from the study samples) at a given SNP was weighted using the estimated brain-region-specific effect of the SNP on gene expression from the GTEx data56. We also accounted for the direction of the co-expression of each gene with SLC6A3 in the network, by multiplying the weight by -1 in case the expression of a gene was negatively correlated with the expression of the SLC6A3 gene in the network – therefore, the higher the score, the higher the expression of the genes that compose the network. The sum of the weighted values from all SNPs for each individual in the cohorts resulted in the region-specific striatal scores. The striatal SLC6A3 expression-based polygenic score (ePGS) is a continuous measure that reflects variation of gene expression of the genes co-expressed with the SLC6A3 gene in the striatum.

Comparison between Polygenic risk scores and ePGS

To compare the results obtained with the striatum SLC6A3 ePGS, we generated traditional polygenic risk scores (PRS) using our accelerated pipeline (https://github.com/MeaneyLab/PRSoS)57. A traditional PRS is a cumulative score calculated based on a relevant GWAS that represents a risk for a certain health outcome or trait58. The sum of the allele count weighted by the effect size across all SNPs in GWAS at a specified threshold was used to calculate type 2 diabetes59 and major depression disorder60 PRSs. The number of SNPs included was defined based on the number of SNPs present in our striatum SLC6A3 ePGS calculated in the discovery cohort. For MDD PRS we used the GWAS results that were obtained without UK Biobank or 23andMe subjects.

Functional enrichment analysis

Enrichment analysis was performed using MetaCore® software from Clarivate Analytics (https://portal.genego.com) to characterize the putative biological functions associated with the striatal SLC6A3 co-expression gene network. Genes that comprise the striatal SLC6A3 ePGS were used in the analysis and the whole genome was used as a background. The significance was considered for the false discovery rate (FDR) adjusted p-value <0.05. To investigate network centrality measures, co-expression patterns were mined from geneMANIA61.The gene interactions were then visualized using the Cytoscape® software62. The nodes are the elements of a network (genes) and edges are the connection between these elements. Bottleneck genes are defined as those having a high betweenness (the extent to which genes act as ‘bridges’ between other genes in a network), hub genes are defined as those with a high degree (genes with more connections to other genes). To analyze the topological properties associated with this gene network, the CentiScaPe app in Cytoscape® was used to calculate the degree and betweenness of each gene. We used this information to define the “hub genes” within the network, characterized as nodes with degrees higher than +1SD above the mean; and the “bottlenecks” characterized as nodes with betweenness higher than +1SD above the mean. A gene that is both bottleneck and hub was considered as a central node of the network63. We also mined protein-protein (PPI) network interactions using the STRING database (https://string-db.org)64 and the striatum SLC6A3 ePGS genes, with the objective to query the physical interactions of the genes that compose our genetic score. Although we mapped the mice co-expressed gene list to human orthologs, not necessarily the co-expression features would be recapitulated in humans. In order to confirm if the genes that comprised the striatal SLC6A3 ePGS are co-expressed in humans and to analyze their patterns of co-expression during different life periods in humans, we used the gene expression data from human postmortem samples from the BrainSpan database65 (see Supplementary information, Supplementary Methods).

Outcome measures

Adult cohort

Psychiatric disorder diagnosis was defined based on the primary or secondary diagnosis of a mental, mood, schizophrenia and neurotic disorders according to participants hospital inpatient records, coded according to the International Classification of Diseases version 10 (ICD-10)66 (UK Biobank field 41270; ICD10 codes: F10-F19 Mental and behavioural disorders due to psychoactive substance use, F20-F29 Schizophrenia, schizotypal and delusional disorders, F30-F39 Mood [affective] disorders, F40-F48 Neurotic, stress-related and somatoform disorders). Cardiometabolic disorders diagnosis was defined by the ICD-10 codes from chapter IV Endocrine, nutritional and metabolic diseases and chapter IX Diseases of the circulatory system (UK Biobank fields: 41270; ICD10 codes: E11-Non-insulin-dependent diabetes, I70-Atherosclerosis, I63-Cerebral infarction, I20-I25 Ischaemic heart diseases). The presence of at least one mental disorder diagnosis and at least one cardiometabolic diagnosis was considered a comorbidity case. Comorbidity variable was coded as a binary variable (1 = “yes” or 0 = “no”). T1 structural brain MRI pre-processed imaging data were generated by an image-processing pipeline developed and run on behalf of the UK Biobank67 (Supplementary information, Supplementary Methods).

Adolescent cohort

No diagnoses for the psychiatric and cardiometabolic disorders noted above were available in ALSPAC. As recommended by the American Academy of Pediatrics (AAP)68, we defined disease risk in adolescents using continuous measures of Total difficulties score measured by the Strengths and Difficulties Questionnaire, depression and anxiety scores measured by Computerized Interview Schedule – Revised (CIS-R), Homeostatic Model Assessment of Insulin Resistance (HOMA-IR), and waist circumference (cm) (Supplementary information, Supplementary Methods, Supplementary Table 2, Supplementary Figure 3). We then characterized two groups of children: low and high cardiometabolic and psychiatric comorbidity risk (see Statistical Analysis).

Statistical Analysis

Statistical analysis were performed using R69. For the descriptive statistics, the ePGS groups were defined by median split, and a comparison between low and high ePGS groups was done using Student t test for continuous variables and a chi-square test for categorical variables (Table 1 and Table 2). Significance levels for all tests were set at p < 0.05.

We performed cluster analysis using the mclust package to construct the comorbidity risk variable in ALSPAC adolescent sample 70. This algorithm applies a model-based classification and density estimation of the z-standardized variables based on finite Gaussian mixture modelling. The method assumes that predictors can be explained by an underlying latent categorical variable (cluster) that represents distinct profiles within the sample, both in a qualitative and quantitative manner. We defined a priori a cluster size solution of two (lower and higher risk for comorbidity). All predictors were z-transformed and adjusted for sex prior to entering the clustering procedure. Regression analysis was carried out to demonstrate the difference between the two clusters in the means for each variable used in the cluster analysis (Supplementary information, Supplementary Table 2). The resulting cluster membership, which represented comorbidity risk, was coded as a binary variable (1= “yes comorbidity” or 0= “no comorbidity”).

The gene by environment (G-E) interaction effect on binary outcomes was explored by logistic regression analysis. Birth weight as a continuous variable was used as a proxy for early life environment exposure in UK Biobank (variable ID20022) and ALSPAC. Early life adversity (E), striatal SLC6A3 ePGS (G) and the interaction term between them were included in the model as main predictors for both cohorts. UK Biobank analyses were also adjusted by sex, age, the first forty genetic principal components, genotyping array, and assessment center, and ALSPAC analyses were adjusted by sex and the first ten genetic principal components. In case of a significant gene by environment (G-E) interaction effect, post hoc simple slope analysis was performed to investigate how the environment effect varies as a function of the genetic background71. The directionality of the G-E effect was explored in the UK Biobank, our discovery cohort, using a two-tailed P-value threshold. The directionality of the G-E effect on ALSPAC was anticipated based on the finding from UK Biobank, thus a one-tailed P-value threshold was considered.

The relation between early life adversity, ePGS and gray matter density in UK Biobank was analyzed in a multivariate parallel independent component analysis (pICA). This analysis was applied to identify the effect of early life adversity on the relation between two different data modalities (genetic and gray matter density) in a data-driven manner72. This analysis separately estimates the maximum independent components within each data modality while also maximizing the association between modalities using an entropy term based on information theory72. Each SNP that composes the striatal SLC6A3 ePGS weighted by striatal GTEx data (genotype * GTEx striatum gene expression slope for each SNP) and whole brain voxel based gray matter density were used in the analysis. Weighted SNPs were adjusted for the genetic principal components (ancestry). The subjects were split in two groups according to the birth weight (low birth weight group: subjects with birth weight <= 2.5kg, n = 953) and a randomly selected group of non-low birth weight individuals (subjects with birth weight > 2.5kg, n = 953, please see https://www.who.int/data/nutrition/nlis/info/low-birth-weight), since there was a large discrepancy between cases and controls sample size within the subsample of individuals with T1 structural brain MRI available. Comparison between low birth weight and randomly selected non-low birth weight individuals on main descriptive variables can be seen on Supplementary information (Supplementary information, Supplementary Table 5). Comparison between the randomly selected group and the full sample of non-low birth weight individuals with MRI available can be seen on Supplementary information (Supplementary information, Supplementary Table 6). T1 structural brain MRI pre-processed images were adjusted by age and sex (See Supplementary information, Supplementary Methods). The Fusion ICA Toolbox (http://mialab.mrn.org/software/fit/) within MATLAB® R2019 was used to run the analysis. The number of independent components was estimated using minimum description length criteria72 for the MRI modality and SNP dimensionality inside the toolbox for the genetic modality. Components for both modalities were converted to z-scores and a threshold at |Z| > 2.5 was used to identify significant brain regions and SNPs that contributed the most for the component overall pattern72. Loading coefficients, which describe the presence of the identified component across subjects72, were extracted for each component, modality, and subject. The mean subject-specific loading coefficients of these components from low birth weight and non-low birth weight groups were compared using Student’s t-test. Talairach coordinates were used to identify the anatomical classification of brain areas included in the identified MRI component73. The significant SNPs (|Z| > 2.5) from the identified genetic component were analyzed using MetaCore®, to identify associated gene ontology processes terms (See Figure 3A for graphical representation of pICA analysis).

Results

Characteristics of the striatal SLC6A3 gene network

We developed a polygenic score to explore the genetic moderation of early life adversity on psychiatric and cardiometabolic comorbidities focusing on a specific gene network (Figure 1A and 1B). We first used brain region-specific RNA sequencing data from mice available at GeneNetwork (http://genenetwork.org/)49 to identify genes co-expressed with the SLC6A3 gene in the striatum. These genes were then converted to human orthologs (Supplementary Table 1 and Figure 1B). This list was used to inform the calculation of the expression-based polygenic score (Striatum SLC6A3 ePGS) in UK Biobank and ALSPAC participants as described in the Methods.

To investigate if mouse-generated SLC6A3 gene network was co-expressed in humans, we queried the gene co-expression patterns of the striatum SLC6A3 gene network throughout human development using gene expression data from human postmortem samples65. A high co-expression was expected in childhood/adolescence, as the striatum SLC6A3 gene network was enriched for genes overexpressed in this period of life (see Figure 1C and Supplementary information, Supplementary Methods). Prominent gene co-expression clusters were also seen in adults (Figure 1C). These findings confirm that the striatal SLC6A3 gene network, originally from murine data, is also co-expressed in humans, and that co-expression is observed at different ages. When visualizing and exploring the network properties (Figure 1D), we observed that the central gene (hub) and the hub-bottleneck genes are related to ribosomal structure. Among the bottleneck genes, which are important connectors between groups of genes, we observed HNRNPA1, which is involved in the packaging of pre-mRNA into particles and transportation from the nucleus to the cytoplasm, as well as splicing. We also observed SDC3 gene, which plays a role in cell shape organization and has been associated with obesity74. Protein-protein interactions of the striatum SLC6A3 co-expression network mined from STRING revealed that the network has significantly more interactions than expected by chance (P<1.0e-16), suggesting that a significant number of the genes in this co-expression network also have physical interactions at the protein level.

The main gene ontology processes terms associated with the network (Supplementary information, Supplementary Figure S4) include: insulin signaling and response terms (Insulin receptor signaling pathway via phosphatidylinositol 3-kinase; Insulin receptor signaling pathway; Cellular response to insulin stimulus; Response to insulin), ribosome production related terms (Ribosome biogenesis; Ribosomal large subunit biogenesis), dopamine receptor signaling pathway (Adenylate cyclase-activating dopamine receptor signaling pathway) and inflammatory response related terms (Regulation of cytokine production involved in inflammatory response; Negative regulation of cytokine production involved in inflammatory response).

Striatum SLC6A3 ePGS moderates the association between birth weight and the risk for psychiatric and cardiometabolic comorbidities in adults

Lower birth weight was associated with the presence of comorbidities in UK Biobank (b= -0.206, Odds ratio (OR)= 0.814, 95% confidence interval (95% CI): 0.781 – 0.847, P <0.001). However, in ALSPAC this association was not significant (b= -0.091, OR= 0.913, 95% CI:0.680 – 1.228, P = 0.548).

For the UK Biobank there was no significant main association of the ePGS with comorbidity (b= 0.006, OR= 1.006, 95% CI: 0.977 – 1.036, P=0.678). In contrast, and consistent with our anticipated hypothesis, there was a significant interaction effect between the striatum SLC6A3 ePGS and birth weight on the presence of psychiatric and cardiometabolic comorbidities in UK Biobank adults (b= 0.042, OR = 1.043, 95% CI: 1.001 – 1.086, P = 0.044). The risk for comorbidity increased as birth weight decreased, especially at lower ePGS values (Low ePGS: b = -0.247, OR = 0.781, P< 0.001, 95% CI 0.738 – 0.826; High ePGS: b = -0.166, OR = 0.847, P< 0.001, 95% CI 0.800 – 0.897).

(Figure 2A). (As we considered birth weight as a continuous variable ranging from low to high values and comorbidity as a dichotomous variable, with the presence of comorbidity computed as 1 and absence as 0, the odds ratio represents the negative association between birth weight and the probability of having comorbidity. Results are presented from the perspective of low birth weight).

Figure 2.
Figure 2. Striatum SLC6A3 ePGS moderates the effect of early life adversity on the risk for mental health disorder and cardiometabolic comorbidity.

Probability of having comorbidity in individuals with high and low striatum SLC6A3 ePGS as a function of birth weight. A, UK Biobank cohort, N=225,972. The risk for comorbidity increases as birth weight decreases, especially at lower ePGS values (Low ePGS: b = - 0.247, OR = 0.781, P< 0.001, 95% CI 0.738 – 0.826; High ePGS: b = -0.166, OR = 0.847, P< 0.001, 95% CI 0.800 – 0.897). B, ALSPAC cohort, N= 1,188. The risk for comorbidity increases as birth weight decreases, especially at lower ePGS values (Low ePGS: b = -0.373, OR = 0.688, 95% CI 0.447 – 1.066, P=0.090; High ePGS: b = 0.176, OR = 1.192, P = 0.419, 95% CI 0.779 – 1.825).

In ALSPAC adolescents the G-E model revealed a significant interaction effect between the striatum SLC6A3 ePGS and birth weight on the probability of belonging to the high comorbidity risk cluster (b = 0.271, OR= 1.311, 95% CI: 1.015 – Inf, P = 0.041, n=1,188). The risk for of belonging to the high comorbidity risk cluster increased as birth weight decreased, especially at lower ePGS values (Low ePGS: b = -0.373, OR = 0.688, 95% CI 0.447 – 1.066, P=0.090; High ePGS: b = 0.176, OR = 1.192, P = 0.419, 95% CI 0.779 – 1.825). (Figure 2B). Similar to the findings in adults, there was no significant main effect association of the striatum SLC6A3 ePGS on the comorbidity risk (b= 0.086, OR= 1.090, 95% CI: 0.957 – 1.240, P= 0.195). These results indicate a developmental trajectory, in which early indicators of risk to develop psychiatric and metabolic comorbidities in adulthood can be seen in adolescents as a function of the interaction of the striatum SLC6A3 co-expression gene network and birth weight.

To benchmark our method against the classical polygenic risk score derived from a GWAS, we performed the same G-E interaction analysis using birth weight and PRSs based on GWAS for major depressive disorder60 and type 2 diabetes59. These are phenotypes related to our main outcome, psychiatric and cardiometabolic comorbidity.

We found significant main effects of Type 2 diabetes and MDD PRSs on the comorbidity outcome in UK Biobank, but not in ALSPAC (Supplementary information, Supplementary Table 3) and no significant G-E interaction on comorbidity using these PRSs in the UK Biobank or ALSPAC (Supplementary information, Supplementary Table 4).

SNPs from the striatum SLC6A3 ePGS are related to gray matter variations in the frontal cortex

We then explored the neuroanatomo-functional relevance of the relation between the striatal SLC6A3 gene network and early adversity. Functional refers to the variation in gene expression represented by the weight attributed to the SNPs that compose the striatum SLC6A3 ePGS. We used a multivariate parallel independent component analysis (pICA)72 (Figure 3A and Supplementary information, Supplementary Methods) and investigated correlations between the SNPs from the striatum SLC6A3 ePGS and voxel-based gray matter density in UK Biobank participants from low birth weight and non-low birth weight groups. This analysis identifies independent components within each data modality separately (SNPs and MRI) while also maximizing the association between these two modalities. The estimated number of components for the MRI modality was 28 and for the genetic modality was 34. Only the most significantly linked pair of components that resulted from the multivariate analysis with higher correlation index value was selected to be further explored: the pair combining the genetic component 13 and MRI component 18 (r=-0.201, p=6.779e-19). A statistically significant difference between birth weight groups was observed for both the genetic component 13 (t=2,214, p=0.026) as well as the MRI component 18 (t=-3,318, p<0.001). These differences between the adversity groups suggest that the relations between data pattern variations (i.e., the relationships between SNPs and gray matter) within this pair of components are significantly different between the two birth weight groups. We then explored the content of genetic component 13 and MRI component 18. The subset of significant SNPs within component 13 is related to variations in gray matter density in the frontal cortex, including the prefrontal cortex, and also more specifically the orbitofrontal cortex, part of the prefrontal cortex, cingulate cortex and temporal cortex (Figure 3B). Enrichment analysis of this subset of significant SNPs (Supplementary information Supplementary Table 7) using MetaCore® (FDR<0.05) showed that the most significant gene ontology enrichment terms are related to regulation of dendrite development, regulation of neuron remodeling, positive regulation of nervous system development, pyruvate biosynthetic process, ATP metabolic process and response to epinephrine (Figure 3C).

Figure 3.
Figure 3. Parallel ICA analysis.

A, schematic representation of parallel ICA method. Two different data modalities (SNPs and voxel-based gray matter) were used to establish anatomical-functional correlations between the striatum SLC6A3 ePGS and brain features from UK Biobank participants (N=11,167). Participants were separated into low birth weight and normal birth weight groups. The analysis estimates the maximum independent components within each data modality separately while also maximizing the association between modalities using an entropy term based on information theory. B, significant brain regions associated with SNPs from the striatum SLC6A3 were frontal cortex, including the orbitofrontal and prefrontal cortex, cingulate cortex and temporal cortex.

Color scheme represents the amount of volume variation (cubic centimeter) significantly associated with the subset of SNPs. C, summary of significant gene ontology processes related to SNPs from the striatum SLC6A3 ePGS associated with gray matter.

Discussion

Our study suggests that being born with lower birth weight increases the risk for later comorbidities between cardiometabolic and psychiatric conditions in adulthood. In fact, being born at low birth weight, which reflects prenatal adversity13, independently associates with increased risk for developing both cardiometabolic15–18 and psychiatric disorders19–24 corroborating our findings. Our functional genomics approach provides evidence for the striatal SLC6A3 co-expression gene network as a salient mechanism moderating this association. This finding is aligned with the critical role of the dopaminergic system in environmental responsivity42,75.

Although lower birth weight is associated with an increased risk for co-morbidity in both ePGS groups in the UK Biobank, low ePGS participants have significantly more risk than high ePGS individuals. In the low ePGS group in adolescence, there is a suggestion of increased risk of being part of a high comorbidity risk as birth weight decreases (p=0.09). Although the simple slope for the high ePGS group in adolescence shows a positive inclination between birth weight and risk for comorbidity, this slope is not significant therefore the risk for comorbidity does not vary according to birth weight in the high ePGS group. No information on gestational age was available in the UK Biobank cohort, to maintain consistency, birth weight as a continuous variable and not corrected for gestational age was used in both cohorts. The lack of information about gestational age in our study may be especially affecting the adolescent analysis and this may potentially explain why the simple slope for the low ePRS group does not reach statistical significance, although the interaction between ePRS and birth weight is statistically significant in ALSPAC.

Our enrichment analysis showed that the striatal SLC6A3 gene network is co-expressed in humans across childhood/adolescence and adulthood (Figure 1C), which is aligned with our interaction between striatum SLC6A3 ePGS and birth weight been observed both in adolescents and in adults. Our results therefore demonstrate that striatum SLC6A3 ePGS is able to detect individual differences in response to early adversity at both ages. Based on the comparisons observed in this study, only the striatum SLC6A3 ePGS was capable of capturing gene by environment interaction effects, while different PRSs did not significantly interact with early life adversity to predict the main outcome. GWAS-derived PRS reflect main genetic effects and thus are unlikely to capture individual differences in response to environmental variation. Indeed we found significant main effects of the PRSs of Major depressive disorder60 and Type 2 diabetes 59 on comorbidity in UK Biobank. Overall, these results align with the well-known capacity of PRS to detect main genetic effects, as well as demonstrate the ability of our ePGS technique in identifying responsivity to environmental change as compared to traditional GWAS-based PRS76. PRS main effects were not observed in adolescents from ALSPAC probably due to the specificity of the GWAS to the features of the original discovery sample; the majority of GWASs are generated based on adult samples77, thus limiting the extrapolation of the effects in different ages.

We also identified putative biological mechanisms underlying the moderating effect of the striatum SLC6A3 ePGS on the association between early adversity and cardiometabolic – psychiatric comorbidities (Figure 1D and Supplementary Figure 4). Central genes of the striatum SLC6A3 network are related to ribosomal structure and the entire gene network is significantly enriched for ribosome production related terms (Ribosome biogenesis; Ribosomal large subunit biogenesis) as seen in our gene-ontology analysis. Ribosomes functioning is highly related to cell growth and proliferation and protein synthesis in all cells. The ribosome biogenesis is a critical process to form mature ribosomes78. Dysfunction of ribosomal gene expression has been seen in animals’ models of depression79 and the use of anti-ribosomal P antibodies, that targets phosphorylated protein (P) components of ribosomes, was able to induce depression like behavior in mice80. Interestingly, the same antibody is used to detect systemic lupus erythematosus in humans80. This autoimmune disease is related to high levels of systemic inflammation. Our SLC6A3 gene network is significantly enriched for inflammatory response related terms, especially cytokine production. One of the central genes of our network, GNB1, has been implicated in regulation of inflammasomes, multiple proteins complexes responsible for activation of inflammatory responses81. In fact, many patients with major depressive disorder have elevated levels of inflammatory cytokines 82 and there is evidence linking cardiometabolic syndrome to higher levels of circulating cytokines83. Not surprisingly, our SLC6A3 gene network is also significantly enriched for dopamine receptor signaling pathway. HNRNPA1, a central node of our network, is involved in mRNA transport and synthesis84. mRNA axonal transport and protein synthesis at the terminal is an important mechanism for regulation of neurotransmitter synthesis and reuptake85. Although SLC6A3 gene expression occurs at the level of the ventral tegmental area, having HNRNPA1 as a central node of our network could explain the presence of SLC6A3 mRNA at the striatal terminal, consistent with numerous human post-mortem studies86–89. SLC6A3 mRNA transport to terminals may be a key mechanistic feature of our SLC6A3 striatum gene network. It has been shown that SLC6A3 protein vesicular traffic has a limited contribution to SLC6A3 concentration in synapses90, hence other forms of regulation of SLC6A3 availability in terminals – for instance via mRNA axonal transport and terminal protein synthesis – are likely in place, in agreement to our findings.

Another important gene in the network is the SDC3, that may play a role in cell shape organization and has been associated with obesity74. Obesity is related to increased risk for cardiovascular disease91, high levels of inflammation and insulin resistance92. Our striatum SLC6A3 gene network is significantly enriched for insulin signaling and response related terms. Being born with low birth weight is associated with insulin resistance in children and adolescents93 and insulin resistance is a risk factor for cardiometabolic and brain-based disorders, including type II diabetes, cardiovascular disease, Alzheimer disease, and major depressive disorder94,95. Metformin, a medication to treat insulin resistance, has shown beneficial psychotropic effects in psychiatric conditions96,97. Evidence shows that insulin has a role in modulating mesocorticolimbic DA neurotransmission through different mechanisms, one of which is increasing DA reuptake by activating the phosphatidylinositol (PI) 3-kinase36,98. Insulin also reduces DA release in rodent nucleus accumbens and medial prefrontal cortex slices99. Our significant gene by environment results, using birth weight as our environmental proxy, corroborate with the literature showing elevated risk for developing psychiatric and cardiometabolic disorders among individuals born with low birth weight. Our genetic enrichment analysis results indicate that insulin signaling disturbances may be a potential mechanism involved in the interaction effect between birth weight and the striatum SLC6A3 co-expression gene network on the risk for cardiometabolic and psychiatric comorbidities. This is aligned with many other studies suggesting that altered insulin function is an important mechanism linking early adversity to later disease100,101.

The subset of SNPs from the striatum SLC6A3 ePGS that is related to gray matter density variations in our neuroanatomical-functional correlation analysis is associated with regulation of dendrite morphogenesis, neuron remodeling and positive regulation of nervous system development. These are important processes linked to the prolonged maturation of the mesocorticolimbic dopamine system during the life-course102, which makes striatal dopaminergic axons especially vulnerable to environmental effects during development103. According to these findings, we recently showed that both rodent poor fetal growth and insulin treatment affect the expression of the Netrin-1/DCC axonal guidance cue system, which is involved in the maturation of the mesocorticolimbic DA circuitry104. Pyruvate biosynthetic process and ATP metabolic process also emerged as significant enrichment terms for the subset of significant SNPs related to gray matter density. Both processes have connections with insulin secretion: regulation of insulin secretion in pancreatic β cells is modulated by ATP synthesis and release in mitochondria105 and by pyruvate transport through mitochondrial pyruvate carriers106. The frontal, prefrontal and orbitofrontal cortices were related to the significant subset of SNPs identified by the pICA analysis. This is aligned with evidence demonstrating that resting state functional connectivity between the orbitofrontal cortex and dorsolateral prefrontal cortex is altered in human individuals born small for gestational age, at different ages during development107.

The cingulate and temporal cortices also emerged as a significant brain regions in the pICA analysis. The anterior cingulate has been implicated in affective abnormalities in mood disorders and volume reduction in patients with major depressive disorder108. Abnormal posterior cingulate functional connectivity has also been reported in major depression109. Temporal lobe alterations are related to insulin resistance pathophysiology in different imaging modalities110. Interestingly, the relationship between genetic and MRI components was significantly different between the two early life adversity groups, suggesting that the biological mechanisms represented by the genetic component and the brain regions highlighted by the MRI component are relevant for the effects of early adversity on adult disease.

Taken together, the evidence suggests that ribosomal function, inflammation, and insulin modulation of dopamine function may be underlying mechanisms by which the striatum SLC6A3 gene network moderates the risk for developing psychiatric and cardiometabolic comorbidities in response to early life adversity. These mechanisms might be especially important in brain areas involving the prefrontal and orbitofrontal cortices, cingulate and temporal cortices.

Our study is limited by the fact that ePGS does not consider intronic regions, potentially ignoring other important regulatory elements. Moreover, our developmental results are based on cross sectional studies, and further longitudinal data are needed to better describe this trajectory.

In sum, we observed that the association between environmental and genetic factors can place individuals at risk for adult comorbid chronic conditions from an early age, and that a striatal dopamine transporter gene network expression has a central role in moderating the association of the early environment with the risk for these diseases. These findings open opportunities for the exploration of the understudied field of precision prevention in pediatrics, and the potential design of more effective interventions and primary care strategies.

Data Availability

For UK Biobank, data can be purchased; the study website contains details of all the data that is available at https://biobank.ndph.ox.ac.uk/showcase/. For ALSPAC, data can be purchased; the study website contains details of all the data that is available through a fully searchable data dictionary and variable search tool at http://www.bristol.ac.uk/alspac/researchers/our-data/.

Acknowledgments

We are extremely grateful to all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC and UK Biobank teams, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses. The UK Medical Research Council and Wellcome (Grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This publication is the work of the authors and Patricia Pelufo Silveira and Barbara Barth will serve as guarantors for the contents of this paper. A comprehensive list of grants funding is available on the ALSPAC website (http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf). This research using ALSPAC was specifically funded by the Wellcome Trust and MRC (Grant ref: 076467/Z/05/Z) and Wellcome Trust (Grant ref: 08426812/Z/07/Z). This project was supported by The JPB Foundation through a grant to The JPB Research Network on Toxic Stress: A Project of the Center on the Developing Child at Harvard University. This work was also funded by the Fonds de recherche du Québec – Santé and Canadian Institutes of Health Research (CIHR, PJT-166066, PI Silveira PP).

Footnotes

  • Update on polygenic scores and analyses following revisions. Changes are reflected in the results, figures, tables, and article title.

References

  1. 1.↵
    Prince MJ, Wu F, Guo Y, et al. The burden of disease in older people and implications for health policy and practice. Lancet 2015; 385(9967): 549-62.
  2. 2.↵
    Uijen AA, van de Lisdonk EH. Multimorbidity in primary care: prevalence and trend over the last 20 years. Eur J Gen Pract 2008; 14 Suppl 1: 28–32.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Egede LE. Major depression in individuals with chronic medical disorders: prevalence, correlates and association with health resource utilization, lost productivity and functional disability. Gen Hosp Psychiatry 2007; 29(5): 409–16.
    OpenUrlCrossRefPubMedWeb of Science
  4. 4.↵
    Tian YE, Di Biase MA, Mosley PE, et al. Evaluation of Brain-Body Health in Individuals With Common Neuropsychiatric Disorders. JAMA psychiatry 2023.
  5. 5.↵
    Souama C, Lamers F, Milaneschi Y, et al. Depression, cardiometabolic disease, and their co-occurrence after childhood maltreatment: an individual participant data meta-analysis including over 200,000 participants. BMC medicine 2023; 21(1): 93.
    OpenUrl
  6. 6.↵
    Vos T, Allen C, Arora M, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. The lancet 2016; 388(10053): 1545-602.
  7. 7.↵
    Kyu HH, Abate D, Abate KH, et al. Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. The Lancet 2018; 392(10159): 1859-922.
  8. 8.↵
    Golden SH, Lazo M, Carnethon M, et al. Examining a bidirectional association between depressive symptoms and diabetes. JAMA 2008; 299(23): 2751–9.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    Moulton CD, Pickup JC, Ismail K. The link between depression and diabetes: the search for shared mechanisms. The Lancet Diabetes & Endocrinology 2015; 3(6): 461–71.
    OpenUrl
  10. 10.↵
    Lotfaliany M, Bowe SJ, Kowal P, Orellana L, Berk M, Mohebbi M. Depression and chronic diseases: Co-occurrence and communality of risk factors. J Affect Disord 2018; 241: 461–8.
    OpenUrlCrossRefPubMed
  11. 11.↵
    De Hert M, Detraux J, Vancampfort D. The intriguing relationship between coronary heart disease and mental disorders. Dialogues in Clinical Neuroscience 2018; 20(1): 31–40.
    OpenUrl
  12. 12.↵
    Regenold WT, Thapar RK, Marano C, Gavirneni S, Kondapavuluru PV. Increased prevalence of type 2 diabetes mellitus among psychiatric inpatients with bipolar I affective and schizoaffective disorders independent of psychotropic drug use. J Affect Disord 2002; 70(1): 19–26.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Lawn JE, Ohuma EO, Bradley E, et al. Small babies, big risks: global estimates of prevalence and mortality for vulnerable newborns to accelerate change and improve counting. The Lancet 2023; 401(10389): 1707–19.
    OpenUrl
  14. 14.↵
    Hanson Ma, Gluckman P. Early developmental conditioning of later health and disease: physiology or pathophysiology? Physiological reviews 2014.
  15. 15.↵
    de Mendonça E, de Lima Macêna M, Bueno NB, de Oliveira ACM, Mello CS. Premature birth, low birth weight, small for gestational age and chronic non-communicable diseases in adult life: A systematic review with meta-analysis. Early Hum Dev 2020; 149: 105154.
  16. 16.
    Wang SF, Shu L, Sheng J, et al. Birth weight and risk of coronary heart disease in adults: a meta-analysis of prospective cohort studies. J Dev Orig Health Dis 2014; 5(6): 408–19.
    OpenUrlCrossRefPubMed
  17. 17.
    Hales CN, Barker DJ. Type 2 (non-insulin-dependent) diabetes mellitus: the thrifty phenotype hypothesis. Diabetologia 1992; 35(7): 595–601.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Rich-Edwards JW, Kleinman K, Michels KB, et al. Longitudinal study of birth weight and adult body mass index in predicting risk of coronary heart disease and stroke in women. Bmj 2005; 330(7500): 1115.
  19. 19.↵
    Betts KS, Williams GM, Najman JM, Scott J, Alati R. The association between lower birth weight and comorbid generalised anxiety and major depressive disorder. J Affect Disord 2013; 146(2): 231–7.
    OpenUrl
  20. 20.
    Wiles NJ, Peters TJ, Leon DA, Lewis G. Birth weight and psychological distress at age 45-51 years: results from the Aberdeen Children of the 1950s cohort study. Br J Psychiatry 2005; 187: 21–8.
    OpenUrlAbstract/FREE Full Text
  21. 21.
    Lahti M, Eriksson J, Heinonen K, et al. Late preterm birth, post-term birth, and abnormal fetal growth as risk factors for severe mental disorders from early to late adulthood. Psychological medicine 2015; 45(5): 985–99.
    OpenUrlPubMed
  22. 22.
    Abel KM, Wicks S, Susser ES, et al. Birth weight, schizophrenia, and adult mental disorder: is risk confined to the smallest babies? Archives of general psychiatry 2010; 67(9): 923–30.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.
    Loret de Mola C, de França GV, Quevedo Lde A, Horta BL. Low birth weight, preterm birth and small for gestational age association with adult depression: systematic review and meta-analysis. Br J Psychiatry 2014; 205(5): 340–7.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    Colman I, Ataullahjan A, Naicker K, Van Lieshout RJ. Birth weight, stress, and symptoms of depression in adolescence: evidence of fetal programming in a national Canadian cohort. Can J Psychiatry 2012; 57(7): 422–8.
    OpenUrlCrossRefPubMed
  25. 25.↵
    Gatzke-Kopp LM. The canary in the coalmine: The sensitivity of mesolimbic dopamine to environmental adversity during development. Neuroscience & Biobehavioral Reviews 2011; 35(3): 794–803.
    OpenUrl
  26. 26.↵
    Gragnoli C, Reeves G, Reazer J, Postolache T. Dopamine–prolactin pathway potentially contributes to the schizophrenia and type 2 diabetes comorbidity. Translational psychiatry 2016; 6(4): e785-e.
    OpenUrl
  27. 27.↵
    Nousen EK, Franco JG, Sullivan EL. Unraveling the mechanisms responsible for the comorbidity between metabolic syndrome and mental health disorders. Neuroendocrinology 2013; 98(4): 254–66.
    OpenUrlCrossRefPubMed
  28. 28.↵
    Alves MB, Laureano DP, Dalle Molle R, et al. Intrauterine growth restriction increases impulsive behavior and is associated with altered dopamine transmission in both medial prefrontal and orbitofrontal cortex in female rats. Physiology & behavior 2019; 204: 336–46.
    OpenUrl
  29. 29.↵
    Alves MB, Dalle Molle R, Desai M, Ross MG, Silveira PP. Increased palatable food intake and response to food cues in intrauterine growth-restricted rats are related to tyrosine hydroxylase content in the orbitofrontal cortex and nucleus accumbens. Behavioural brain research 2015; 287: 73–81.
    OpenUrl
  30. 30.↵
    Laureano DP, Alves MB, Miguel PM, et al. Intrauterine growth restriction modifies the accumbal dopaminergic response to palatable food intake. Neuroscience 2018.
  31. 31.↵
    Pizzagalli DA, Berretta S, Wooten D, et al. Assessment of striatal dopamine transporter binding in individuals with major depressive disorder: in vivo positron emission tomography and postmortem evidence. JAMA psychiatry 2019; 76(8): 854–61.
    OpenUrl
  32. 32.↵
    Jones KT, Woods C, Zhen J, Antonio T, Carr KD, Reith ME. Effects of diet and insulin on dopamine transporter activity and expression in rat caudate-putamen, nucleus accumbens, and midbrain. Journal of neurochemistry 2017; 140(5): 728–40.
    OpenUrlCrossRefPubMed
  33. 33.↵
    ter Horst KW, Lammers NM, Trinko R, et al. Striatal dopamine regulates systemic glucose metabolism in humans and mice. Science Translational Medicine 2018; 10(442): eaar3752.
  34. 34.↵
    Huot P, Parent A. Dopaminergic neurons intrinsic to the striatum. Journal of neurochemistry 2007; 101(6): 1441–7.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    Threlfell S, Mohammadi AS, Ryan BJ, et al. Striatal Dopamine Transporter Function Is Facilitated by Converging Biology of α-Synuclein and Cholesterol. Front Cell Neurosci 2021; 15: 658244.
  36. 36.↵
    Garcia BG, Wei Y, Moron JA, Lin RZ, Javitch JA, Galli A. Akt is essential for insulin modulation of amphetamine-induced human dopamine transporter cell-surface redistribution. Molecular pharmacology 2005; 68(1): 102–9.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    Sudlow C, Gallacher J, Allen N, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med 2015; 12(3): e1001779.
    OpenUrlCrossRefPubMed
  38. 38.↵
    Fraser A, Macdonald-Wallis C, Tilling K, et al. Cohort profile: the Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. International journal of epidemiology 2013; 42(1): 97–110.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.
    Boyd A, Golding J, Macleod J, et al. Cohort profile: the ‘children of the 90s’—the index offspring of the Avon Longitudinal Study of Parents and Children. International journal of epidemiology 2013; 42(1): 111–27.
    OpenUrlCrossRefPubMedWeb of Science
  40. 40.↵
    Golding J, Pembrey M, Jones R. ALSPAC--the Avon Longitudinal Study of Parents and Children. I. Study methodology. Paediatric and perinatal epidemiology 2001; 15(1): 74–87.
    OpenUrlCrossRef
  41. 41.↵
    Hari Dass SA, McCracken K, Pokhvisneva I, et al. A biologically-informed polygenic score identifies endophenotypes and clinical conditions associated with the insulin receptor function on specific brain regions. EBioMedicine 2019; 42: 13.
  42. 42.↵
    Miguel PM, Pereira LO, Barth B, et al. Prefrontal cortex dopamine transporter gene network moderates the effect of perinatal hypoxic-ischemic conditions on cognitive flexibility and brain gray matter density in children. Biological psychiatry 2019; 86(8): 621–30.
    OpenUrl
  43. 43.
    Restrepo-Lozano JM, Pokhvisneva I, Wang Z, et al. Corticolimbic DCC gene co-expression networks as predictors of impulsivity in children. Molecular Psychiatry 2022: 1–9.
  44. 44.
    Silveira PP, Pokhvisneva I, Parent C, et al. Cumulative prenatal exposure to adversity reveals associations with a broad range of neurodevelopmental outcomes that are moderated by a novel, biologically informed polygenetic score based on the serotonin transporter solute carrier family C6, member 4 (SLC6A4) gene expression. Development psychopathology 2017; 29(5): 1601–17.
    OpenUrl
  45. 45.
    de Lima RMS, Barth B, Arcego DM, et al. Amygdala 5-HTT Gene Network Moderates the Effects of Postnatal Adversity on Attention Problems: Anatomo-Functional Correlation and Epigenetic Changes. Front Neurosci 2020; 14: 198.
  46. 46.
    de Lima RMS, Barth B, Mar Arcego D, et al. Leptin receptor co-expression gene network moderates the effect of early life adversity on eating behavior in children. Communications Biology 2022; 5(1): 1092.
    OpenUrl
  47. 47.
    de Mendonça Filho EJ, Barth B, Bandeira DR, et al. Cognitive development and brain gray matter susceptibility to prenatal adversities: Moderation by the prefrontal cortex brain-derived neurotrophic factor gene co-expression network. Frontiers in Neuroscience 2021; 15: 744743.
  48. 48.↵
    Dalmaz C, Barth B, Pokhvisneva I, et al. Prefrontal cortex VAMP1 gene network moderates the effect of the early environment on cognitive flexibility in children. Neurobiology of Learning and Memory 2021: 107509.
  49. 49.↵
    Mulligan MK, Mozhui K, Prins P, Williams RW. GeneNetwork: A Toolbox for Systems Genetics. Methods Mol Biol 2017; 1488: 75–120.
  50. 50.↵
    Miguel PM, Deniz BF, Deckmann I, et al. Prefrontal cortex dysfunction in hypoxic-ischaemic encephalopathy contributes to executive function impairments in rats: potential contribution for attention-deficit/hyperactivity disorder. The World Journal of Biological Psychiatry 2018; 19(7): 547–60.
    OpenUrl
  51. 51.
    Dalle Molle R, Laureano DP, Alves MB, et al. Intrauterine growth restriction increases the preference for palatable foods and affects sensitivity to food rewards in male and female adult rats. brain research 2015; 1618: 41–9.
  52. 52.
    Cunha FdS, Dalle Molle R, Portella AK, et al. Both food restriction and high-fat diet during gestation induce low birth weight and altered physical activity in adult rat offspring: the “Similarities in the Inequalities” model. PLoS One 2015; 10(3): e0118586.
    OpenUrlCrossRefPubMed
  53. 53.
    Portella AK, Silveira PP, Laureano DP, et al. Litter size reduction alters insulin signaling in the ventral tegmental area and influences dopamine-related behaviors in adult rats. Behavioural brain research 2015; 278: 66–73.
    OpenUrl
  54. 54.↵
    Silveira PP, Portella AK, Assis SACdN, et al. Early life experience alters behavioral responses to sweet food and accumbal dopamine metabolism. International journal of developmental neuroscience 2010; 28(1): 111–8.
    OpenUrlCrossRefPubMed
  55. 55.↵
    Durinck S, Spellman PT, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nature protocols 2009; 4(8): 1184.
    OpenUrl
  56. 56.↵
    Consortium GT. The Genotype-Tissue Expression (GTEx) project. Nat Genet 2013; 45(6): 580–5.
    OpenUrlCrossRefPubMed
  57. 57.↵
    Chen LM, Yao N, Garg E, et al. PRS-on-Spark (PRSoS): a novel, efficient and flexible approach for generating polygenic risk scores. BMC Bioinformatics 2018; 19(1): 295.
    OpenUrl
  58. 58.↵
    Wray NR, Goddard ME, Visscher PM. Prediction of individual genetic risk to disease from genome-wide association studies. Genome Res 2007; 17(10): 1520–8.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    Cai L, Wheeler E, Kerrison ND, et al. Genome-wide association analysis of type 2 diabetes in the EPIC-InterAct study. Scientific Data 2020; 7(1): 393.
    OpenUrl
  60. 60.↵
    Wray NR, Ripke S, Mattheisen M, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature genetics 2018; 50(5): 668.
    OpenUrlCrossRefPubMed
  61. 61.↵
    Warde-Farley D, Donaldson SL, Comes O, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Research 2010; 38(suppl_2): W214-W20.
  62. 62.↵
    Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome research 2003; 13(11): 2498–504.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    Neves de Oliveira B-H, Dalmaz C, Zeidán-Chuliá F. Network-based identification of altered stem cell pluripotency and calcium signaling pathways in metastatic melanoma. Medical Sciences 2018; 6(1): 23.
    OpenUrl
  64. 64.↵
    Franceschini A, Szklarczyk D, Frankild S, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res 2013; 41(Database issue): D808-15.
  65. 65.↵
    Miller JA, Ding S-L, Sunkin SM, et al. Transcriptional landscape of the prenatal human brain. Nature 2014; 508(7495): 199.
  66. 66.↵
    Organization WH. ICD-10 : international statistical classification of diseases and related health problems : tenth revision. 2 ed; 2004.
  67. 67.↵
    Alfaro-Almagro F, Jenkinson M, Bangerter NK, et al. Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank. Neuroimage 2018; 166: 400–24.
    OpenUrlCrossRefPubMed
  68. 68.↵
    Magge SN, Goodman E, Armstrong SC, et al. The metabolic syndrome in children and adolescents: shifting the focus to cardiometabolic risk factor clustering. Pediatrics 2017; 140(2).
  69. 69.↵
    R CT. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014.
  70. 70.↵
    Scrucca L, Fop M, Murphy TB, Raftery AE. mclust 5: clustering, classification and density estimation using Gaussian finite mixture models. The R journal 2016; 8(1): 289.
    OpenUrl
  71. 71.↵
    Bauer DJ, Curran PJ. Probing interactions in fixed and multilevel regression: Inferential and graphical techniques. Multivariate behavioral research 2005; 40(3): 373–400.
    OpenUrlCrossRefPubMedWeb of Science
  72. 72.↵
    Pearlson GD, Liu J, Calhoun VD. An introductory review of parallel independent component analysis (p-ICA) and a guide to applying p-ICA to genetic data and imaging phenotypes to identify disease-associated biological pathways and systems in common complex disorders. Front Genet 2015; 6: 276.
  73. 73.↵
    Talairach J, Tournoux P. A coplanar stereotaxic atlas of the human brain: An approach to medical cerebral imaging. Thieme Medical Publishers, Inc., New York.; 1988.
  74. 74.↵
    Schüring AN, Lutz F, Tüttelmann F, Gromoll J, Kiesel L, Götte M. Role of syndecan-3 polymorphisms in obesity and female hyperandrogenism. J Mol Med (Berl) 2009; 87(12): 1241–50.
    OpenUrlCrossRefPubMed
  75. 75.↵
    Dahoun T, Nour MM, McCutcheon RA, Adams RA, Bloomfield MA, Howes OD. The relationship between childhood trauma, dopamine release and dexamphetamine-induced positive psychotic symptoms: a [11 C]-(+)-PHNO PET study. Translational psychiatry 2019; 9(1): 1–12.
    OpenUrl
  76. 76.↵
    Silveira PP, Meaney MJ. Examining the biological mechanisms of human mental disorders resulting from gene-environment interdependence using novel functional genomic approaches. Neurobiology of Disease 2023; 178: 106008.
  77. 77.↵
    Sirugo G, Williams SM, Tishkoff SA. The missing diversity in human genetic studies. Cell 2019; 177(1): 26–31.
    OpenUrlPubMed
  78. 78.↵
    Kang J, Brajanovski N, Chan KT, Xuan J, Pearson RB, Sanij E. Ribosomal proteins and human diseases: molecular mechanisms and targeted therapy. Signal Transduct Target Ther 2021; 6(1): 323.
    OpenUrl
  79. 79.↵
    Smagin DA, Kovalenko IL, Galyamina AG, Bragin AO, Orlov YL, Kudryavtseva NN. Dysfunction in Ribosomal Gene Expression in the Hypothalamus and Hippocampus following Chronic Social Defeat Stress in Male Mice as Revealed by RNA-Seq. Neural Plast 2016; 2016: 3289187.
    OpenUrl
  80. 80.↵
    Katzav A, Solodeev I, Brodsky O, et al. Induction of autoimmune depression in mice by anti-ribosomal P antibodies via the limbic system. Arthritis Rheum 2007; 56(3): 938–48.
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    Murakami T, Ruengsinpinya L, Nakamura E, et al. Cutting Edge: G Protein Subunit β 1 Negatively Regulates NLRP3 Inflammasome Activation. J Immunol 2019; 202(7): 1942–7.
    OpenUrlAbstract/FREE Full Text
  82. 82.↵
    Roohi E, Jaafari N, Hashemian F. On inflammatory hypothesis of depression: what is the role of IL-6 in the middle of the chaos? Journal of Neuroinflammation 2021; 18(1): 45.
    OpenUrl
  83. 83.↵
    Karhunen V, Gill D, Huang J, et al. The interplay between inflammatory cytokines and cardiometabolic disease: bi-directional mendelian randomisation study. BMJ Medicine 2023; 2(1): e000157.
    OpenUrlAbstract/FREE Full Text
  84. 84.↵
    Dreyfuss G, Kim VN, Kataoka N. Messenger-RNA-binding proteins and the messages they carry. Nature reviews Molecular cell biology 2002; 3(3): 195–205.
    OpenUrlCrossRefPubMedWeb of Science
  85. 85.↵
    Dalla Costa I, Buchanan CN, Zdradzinski MD, et al. The functional organization of axonal mRNA transport and translation. Nature Reviews Neuroscience 2020: 1–15.
  86. 86.↵
    Choi CS, Hong M, Kim KC, et al. Effects of atomoxetine on hyper-locomotive activity of the prenatally valproate-exposed rat offspring. Biomolecules & therapeutics 2014; 22(5): 406.
    OpenUrl
  87. 87.
    Jijun L, Zaiwang L, Anyuan L, et al. Abnormal expression of dopamine and serotonin transporters associated with the pathophysiologic mechanism of Tourette syndrome. Neurology India 2010; 58(4): 523.
    OpenUrlPubMed
  88. 88.
    Sun Z, Ma Y, Xie L, et al. Behavioral changes and neuronal damage in rhesus monkeys after 10 weeks of ketamine administration involve prefrontal cortex dopamine D2 receptor and dopamine transporter. Neuroscience 2019; 415: 97–106.
    OpenUrl
  89. 89.↵
    Ciliax BJ, Drash GW, Staley JK, et al. Immunocytochemical localization of the dopamine transporter in human brain. Journal of Comparative Neurology 1999; 409(1): 38–56.
    OpenUrlCrossRefPubMedWeb of Science
  90. 90.↵
    Bagalkot TR, Block ER, Bucchin K, et al. Dopamine transporter localization in medial forebrain bundle axons indicates its long-range transport primarily by membrane diffusion with a limited contribution of vesicular traffic on retromer-positive compartments. Journal of Neuroscience 2021; 41(2): 234–50.
    OpenUrlAbstract/FREE Full Text
  91. 91.↵
    Mathieu P, Lemieux I, Després JP. Obesity, inflammation, and cardiovascular risk. Clinical pharmacology & therapeutics 2010; 87(4): 407–16.
    OpenUrl
  92. 92.↵
    Shoelson SE, Herrero L, Naaz A. Obesity, inflammation, and insulin resistance. Gastroenterology 2007; 132(6): 2169–80.
    OpenUrlCrossRefPubMedWeb of Science
  93. 93.↵
    Martín-Calvo N, Goni L, Tur JA, Martínez JA. Low birth weight and small for gestational age are associated with complications of childhood and adolescence obesity: Systematic review and meta-analysis. Obes Rev 2022; 23 Suppl 1: e13380.
    OpenUrl
  94. 94.↵
    Watson KT, Simard JF, Henderson VW, et al. Association of Insulin Resistance With Depression Severity and Remission Status: Defining a Metabolic Endophenotype of Depression. JAMA psychiatry 2021; 78(4): 439–41.
    OpenUrl
  95. 95.↵
    Watson K, Nasca C, Aasly L, McEwen B, Rasgon N. Insulin resistance, an unmasked culprit in depressive disorders: Promises for interventions. Neuropharmacology 2018; 136: 327–34.
    OpenUrlCrossRefPubMed
  96. 96.↵
    Calkin CV, Chengappa KNR, Cairns K, et al. Treating Insulin Resistance With Metformin as a Strategy to Improve Clinical Outcomes in Treatment-Resistant Bipolar Depression (the TRIO-BD Study): A Randomized, Quadruple-Masked, Placebo-Controlled Clinical Trial. J Clin Psychiatry 2022; 83(2).
  97. 97.↵
    Erensoy H, Niafar M, Ghafarzadeh S, Aghamohammadzadeh N, Nader ND. A pilot trial of metformin for insulin resistance and mood disturbances in adolescent and adult women with polycystic ovary syndrome. Gynecol Endocrinol 2019; 35(1): 72–5.
    OpenUrl
  98. 98.↵
    Carvelli L, Morón JA, Kahlig KM, et al. PI 3-kinase regulation of dopamine uptake. Journal of neurochemistry 2002; 81(4): 859–69.
    OpenUrlCrossRefPubMedWeb of Science
  99. 99.↵
    Schoffelmeer AN, Drukarch B, De Vries TJ, Hogenboom F, Schetters D, Pattij T. Insulin modulates cocaine-sensitive monoamine transporter function and impulsive behavior. Journal of Neuroscience 2011; 31(4): 1284–91.
    OpenUrlAbstract/FREE Full Text
  100. 100.↵
    Batra A, Latsko M, Portella AK, Silveira PP. Early adversity and insulin: neuroendocrine programming beyond glucocorticoids. Trends Endocrinol Metab 2021; 32(12): 1031–43.
    OpenUrl
  101. 101.↵
    Alberry B, Silveira PP. Brain insulin signaling as a potential mediator of early life adversity effects on physical and mental health. Neuroscience & Biobehavioral Reviews 2023; 153: 105350.
  102. 102.↵
    Vosberg DE, Leyton M, Flores C. The Netrin-1/DCC guidance system: dopamine pathway maturation and psychiatric disorders emerging in adolescence. Mol Psychiatry 2020; 25(2): 297–307.
    OpenUrlCrossRefPubMed
  103. 103.↵
    Hoops D, Flores C. Making dopamine connections in adolescence. Trends in neurosciences 2017; 40(12): 709–19.
    OpenUrlCrossRefPubMed
  104. 104.↵
    Batra A, Cuesta S, Alves MB, et al. Relationship between insulin and Netrin-1/DCC guidance cue pathway regulation in the prefrontal cortex of rodents exposed to prenatal dietary restriction. Journal of Developmental Origins of Health and Disease 2023; 14(4): 501–7.
    OpenUrl
  105. 105.↵
    Wang C, Geng B, Cui Q, Guan Y, Yang J. Intracellular and extracellular adenosine triphosphate in regulation of insulin secretion from pancreatic β cells (β). J Diabetes 2014; 6(2): 113–9.
    OpenUrl
  106. 106.↵
    Patterson JN, Cousteils K, Lou JW, Fox JEM, MacDonald PE, Joseph JW. Mitochondrial metabolism of pyruvate is essential for regulating glucose-stimulated insulin secretion. Journal of Biological Chemistry 2014; 289(19): 13335–46.
    OpenUrlAbstract/FREE Full Text
  107. 107.↵
    Dalle Molle R, de Mendonça Filho EJ, Minuzzi L, et al. Thrifty-Eating Behavior Phenotype at the Food Court - Programming Goes Beyond Food Preferences. Front Endocrinol (Lausanne) 2022; 13: 882532.
  108. 108.↵
    Yucel K, McKinnon MC, Chahal R, et al. Anterior Cingulate Volumes in Never-Treated Patients with Major Depressive Disorder. Neuropsychopharmacology 2008; 33(13): 3157–63.
    OpenUrlCrossRefPubMedWeb of Science
  109. 109.↵
    Leech R, Sharp DJ. The role of the posterior cingulate cortex in cognition and disease. Brain 2013; 137(1): 12-32.
  110. 110.↵
    Cui Y, Tang TY, Lu CQ, Ju S. Insulin resistance and cognitive impairment: evidence from neuroimaging. Journal of Magnetic Resonance Imaging 2022; 56(6): 1621–49.
    OpenUrl
Back to top
Posted December 15, 2023.
Download PDF

Supplementary Material

Data/Code
Email
Striatal dopamine gene network moderates the effect of early adversity on the risk for adult psychiatric and cardiometabolic comorbidity
Barbara Barth, Danusa Mar Arcego, Euclides José de Mendonça Filho, Randriely Merscher Sobreira de Lima, Carine Parent, Carla Dalmaz, André Krumel Portella, Irina Pokhvisneva, Michael J Meaney, Patricia Pelufo Silveira
medRxiv 2022.04.23.22274209; doi: https://doi.org/10.1101/2022.04.23.22274209
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area