Human archetypal pluripotent stem cells differentiate into trophoblast stem cells via endogenous BMP5/7 induction without transitioning through naive state.

Primary human trophoblast stem cells (TSCs) and TSCs derived from human pluripotent stem cells (hPSCs) can potentially model placental processes in vitro. Yet, the pluripotent states and factors involved in the differentiation of hPSCs to TSCs remain poorly understood. In this study, we demonstrate that the primed pluripotent state can generate TSCs by activating pathways such as Epidermal Growth Factor (EGF) and Wingless-related integration site (WNT), and by suppressing tumor growth factor beta (TGFß), histone deacetylases (HDAC), and Rho-associated protein kinase (ROCK) signaling pathways, all without the addition of exogenous Bone morphogenetic protein 4 (BMP4)-a condition we refer to as the TS condition. We characterized this process using temporal single-cell RNA sequencing to compare TS conditions with differentiation protocols involving BMP4 activation alone or BMP4 activation in conjunction with WNT inhibition. The TS condition consistently produced a stable, proliferative cell type that closely mimics first-trimester placental cytotrophoblasts, marked by the activation of endogenous retroviral genes and the absence of amnion expression. This was observed across multiple cell lines, including various primed induced pluripotent stem cell (iPSC) and embryonic stem cell (ESC) lines. Primed-derived TSCs can proliferate for over 30 passages and further specify into multinucleated syncytiotrophoblasts and extravillous trophoblast cells. Our research establishes that the differentiation of primed hPSCs to TSC under TS conditions triggers the induction of TMSB4X, BMP5/7, GATA3, and TFAP2A without progressing through a naive state. These findings propose that the primed hPSC state is part of a continuum of potency with the capacity to differentiate into TSCs through multiple routes.

A miR-137-Related Biological Pathway of Risk for Schizophrenia Is Associated With Human Brain Emotion Processing.

BACKGROUND: miR-137 is a microRNA involved in brain development, regulating neurogenesis and neuronal maturation. Genome-wide association studies have implicated miR-137 in schizophrenia risk but do not explain its involvement in brain function and underlying biology. Polygenic risk for schizophrenia mediated by miR-137 targets is associated with working memory, although other evidence points to emotion processing. We characterized the functional brain correlates of miR-137 target genes associated with schizophrenia while disentangling previously reported associations of miR-137 targets with working memory and emotion processing. METHODS: Using RNA sequencing data from postmortem prefrontal cortex (N = 522), we identified a coexpression gene set enriched for miR-137 targets and schizophrenia risk genes. We validated the relationship of this set to miR-137 in vitro by manipulating miR-137 expression in neuroblastoma cells. We translated this gene set into polygenic scores of coexpression prediction and associated them with functional magnetic resonance imaging activation in healthy volunteers (n1 = 214; n2 = 136; n3 = 2075; n4 = 1800) and with short-term treatment response in patients with schizophrenia (N = 427). RESULTS: In 4652 human participants, we found that 1) schizophrenia risk genes were coexpressed in a biologically validated set enriched for miR-137 targets; 2) increased expression of miR-137 target risk genes was mediated by low prefrontal miR-137 expression; 3) alleles that predict greater gene set coexpression were associated with greater prefrontal activation during emotion processing in 3 independent healthy cohorts (n1, n2, n3) in interaction with age (n4); and 4) these alleles predicted less improvement in negative symptoms following antipsychotic treatment in patients with schizophrenia. CONCLUSIONS: The functional translation of miR-137 target gene expression linked with schizophrenia involves the neural substrates of emotion processing.

Investigating trait variability of gene co-expression network architecture in brain by controlling for genomic risk of schizophrenia.

The effect of schizophrenia (SCZ) genetic risk on gene expression in brain remains elusive. A popular approach to this problem has been the application of gene co-expression network algorithms (e.g., WGCNA). To improve reliability with this method it is critical to remove unwanted sources of variance while also preserving biological signals of interest. In this WCGNA study of RNA-Seq data from postmortem prefrontal cortex (78 neurotypical donors, EUR ancestry), we tested the effects of SCZ genetic risk on co-expression networks. Specifically, we implemented a novel design in which gene expression was adjusted by linear regression models to preserve or remove variance explained by biological signal of interest (GWAS genomic scores for SCZ risk-(GS-SCZ), and genomic scores- GS of height (GS-Ht) as a negative control), while removing variance explained by covariates of non-interest. We calculated co-expression networks from adjusted expression (GS-SCZ and GS-Ht preserved or removed), and consensus between them (representative of a "background" network free of genomic scores effects). We then tested the overlap between GS-SCZ preserved modules and background networks reasoning that modules with reduced overlap would be most affected by GS-SCZ biology. Additionally, we tested these modules for convergence of SCZ risk (i.e., enrichment in PGC3 SCZ GWAS priority genes, enrichment in SCZ risk heritability and relevant biological ontologies. Our results highlight key aspects of GS-SCZ effects on brain co-expression networks, specifically: 1) preserving/removing SCZ genetic risk alters the co-expression modules; 2) biological pathways enriched in modules affected by GS-SCZ implicate processes of transcription, translation and metabolism that converge to influence synaptic transmission; 3) priority PGC3 SCZ GWAS genes and SCZ risk heritability are enriched in modules associated with GS-SCZ effects. Overall, our results indicate that gene co-expression networks that selectively integrate information about genetic risk can reveal novel combinations of biological pathways involved in schizophrenia.

Heritability of amygdala reactivity to angry faces and its replicable association with the schizophrenia risk locus of miR-137.

BACKGROUND: Among healthy participants, the interindividual variability of brain response to facial emotions is associated with genetic variation, including common risk variants for schizophrenia, a heritable brain disorder characterized by anomalies in emotion processing. We aimed to identify genetic variants associated with heritable brain activity during processing of facial emotions among healthy participants and to explore the impact of these identified variants among patients with schizophrenia. METHODS: We conducted a data-driven stepwise study including samples of healthy twins, unrelated healthy participants and patients with schizophrenia. Participants approached or avoided pictures of faces with negative emotional valence during functional magnetic resonance imaging (fMRI). RESULTS: We investigated 3 samples of healthy participants - including 28 healthy twin pairs, 289 unrelated healthy participants (genome-wide association study [GWAS] discovery sample) and 90 unrelated healthy participants (replication sample) - and 1 sample of 48 patients with schizophrenia. Among healthy twins, we identified the amygdala as the brain region with the highest heritability during processing of angry faces (heritability estimate 0.54, p < 0.001). Subsequent GWAS in both discovery and replication samples of healthy non-twins indicated that amygdala activity was associated with a polymorphism in the miR-137 locus (rs1198575), a micro-RNA strongly involved in risk for schizophrenia. A significant effect in the same direction was found among patients with schizophrenia (p = 0.03). LIMITATIONS: The limited sample size available for GWAS analyses may require further replication of results. CONCLUSION: Our data-driven approach shows preliminary evidence that amygdala activity, as evaluated with our task, is heritable. Our genetic associations preliminarily suggest a role for miR-137 in brain activity during explicit processing of facial emotions among healthy participants and patients with schizophrenia, pointing to the amygdala as a brain region whose activity is related to miR-137.

Prioritization of potential causative genes for schizophrenia in placenta.

Our earlier work has shown that genomic risk for schizophrenia converges with early life complications in affecting risk for the disorder and sex-biased neurodevelopmental trajectories. Here, we identify specific genes and potential mechanisms that, in placenta, may mediate such outcomes. We performed TWAS in healthy term placentae (N = 147) to derive candidate placental causal genes that we confirmed with SMR; to search for placenta and schizophrenia-specific associations, we performed an analogous analysis in fetal brain (N = 166) and additional placenta TWAS for other disorders/traits. The analyses in the whole sample and stratifying by sex ultimately highlight 139 placenta and schizophrenia-specific risk genes, many being sex-biased; the candidate molecular mechanisms converge on the nutrient-sensing capabilities of placenta and trophoblast invasiveness. These genes also implicate the Coronavirus-pathogenesis pathway and showed increased expression in placentae from a small sample of SARS-CoV-2-positive pregnancies. Investigating placental risk genes for schizophrenia and candidate mechanisms may lead to opportunities for prevention that would not be suggested by study of the brain alone.

Genetics and Brain Transcriptomics of Completed Suicide.

OBJECTIVE: The authors sought to study the transcriptomic and genomic features of completed suicide by parsing the method chosen, to capture molecular correlates of the distinctive frame of mind of individuals who die by suicide, while reducing heterogeneity. METHODS: The authors analyzed gene expression (RNA sequencing) from postmortem dorsolateral prefrontal cortex of patients who died by suicide with violent compared with nonviolent means, nonsuicide patients with the same psychiatric disorders, and a neurotypical group (total N=329). They then examined genomic risk scores (GRSs) for each psychiatric disorder included, and GRSs for cognition (IQ) and for suicide attempt, testing how they predict diagnosis or traits (total N=888). RESULTS: Patients who died by suicide by violent means showed a transcriptomic pattern remarkably divergent from each of the other patient groups but less from the neurotypical group; consistently, their genomic profile of risk was relatively low for their diagnosed illness as well as for suicide attempt, and relatively high for IQ: they were more similar to the neurotypical group than to other patients. Differentially expressed genes (DEGs) associated with patients who died by suicide by violent means pointed to purinergic signaling in microglia, showing similarities to a genome-wide association study of Drosophila aggression. Weighted gene coexpression network analysis revealed that these DEGs were coexpressed in a context of mitochondrial metabolic activation unique to suicide by violent means. CONCLUSIONS: These findings suggest that patients who die by suicide by violent means are in part biologically separable from other patients with the same diagnoses, and their behavioral outcome may be less dependent on genetic risk for conventional psychiatric disorders and be associated with an alteration of purinergic signaling and mitochondrial metabolism.

A multimodal study of a first episode psychosis cohort: potential markers of antipsychotic treatment resistance.

Treatment resistant (TR) psychosis is considered to be a significant cause of disability and functional impairment. Numerous efforts have been made to identify the clinical predictors of TR. However, the exploration of molecular and biological markers is still at an early stage. To understand the TR condition and identify potential molecular and biological markers, we analyzed demographic information, clinical data, structural brain imaging data, and molecular brain imaging data in 7 Tesla magnetic resonance spectroscopy from a first episode psychosis cohort that includes 136 patients. Age, gender, race, smoking status, duration of illness, and antipsychotic dosages were controlled in the analyses. We found that TR patients had a younger age at onset, more hospitalizations, more severe negative symptoms, a reduction in the volumes of the hippocampus (HP) and superior frontal gyrus (SFG), and a reduction in glutathione (GSH) levels in the anterior cingulate cortex (ACC), when compared to non-TR patients. The combination of multiple markers provided a better classification between TR and non-TR patients compared to any individual marker. Our study shows that ACC-GSH, HP and SFG volumes, and age at onset, could potentially be biomarkers for TR diagnosis, while hospitalization and negative symptoms could be used to evaluate the progression of the disease. Multimodal cohorts are essential in obtaining a comprehensive understanding of brain disorders.

The interaction between cannabis use and a CB1-related polygenic co-expression index modulates dorsolateral prefrontal activity during working memory processing.

Convergent findings indicate that cannabis use and variation in the cannabinoid CB1 receptor coding gene (CNR1) modulate prefrontal function during working memory (WM). Other results also suggest that cannabis modifies the physiological relationship between genetically induced expression of CNR1 and prefrontal WM processing. However, it is possible that cannabis exerts its modifying effect on prefrontal physiology by interacting with complex molecular ensembles co-regulated with CB1. Since co-regulated genes are likely co-expressed, we investigated how genetically predicted co-expression of a molecular network including CNR1 interacts with cannabis use in modulating WM processing in humans. Using post-mortem human prefrontal data, we first computed a polygenic score (CNR1-PCI), combining the effects of single nucleotide polymorphisms (SNPs) on co-expression of a cohesive gene set including CNR1, and positively correlated with such co-expression. Then, in an in vivo study, we computed CNR1-PCI in 88 cannabis users and 147 non-users and investigated its interaction with cannabis use on brain activity during WM. Results revealed an interaction between cannabis use and CNR1-PCI in the dorsolateral prefrontal cortex (DLPFC), with a positive relationship between CNR1-PCI and DLPFC activity in cannabis users and a negative relationship in non-users. Furthermore, DLPFC activity in cannabis users was positively correlated with the frequency of cannabis use. Taken together, our results suggest that co-expression of a CNR1-related network predicts WM-related prefrontal activation as a function of cannabis use. Furthermore, they offer novel insights into the biological mechanisms associated with the use of cannabis.

The genetic architecture of human brainstem structures and their involvement in common brain disorders.

Brainstem regions support vital bodily functions, yet their genetic architectures and involvement in common brain disorders remain understudied. Here, using imaging-genetics data from a discovery sample of 27,034 individuals, we identify 45 brainstem-associated genetic loci, including the first linked to midbrain, pons, and medulla oblongata volumes, and map them to 305 genes. In a replication sample of 7432 participants most of the loci show the same effect direction and are significant at a nominal threshold. We detect genetic overlap between brainstem volumes and eight psychiatric and neurological disorders. In additional clinical data from 5062 individuals with common brain disorders and 11,257 healthy controls, we observe differential volume alterations in schizophrenia, bipolar disorder, multiple sclerosis, mild cognitive impairment, dementia, and Parkinson's disease, supporting the relevance of brainstem regions and their genetic architectures in common brain disorders.

Multivariate patterns of gray matter volume in thalamic nuclei are associated with positive schizotypy in healthy individuals.

BACKGROUND: Previous models suggest biological and behavioral continua among healthy individuals (HC), at-risk condition, and full-blown schizophrenia (SCZ). Part of these continua may be captured by schizotypy, which shares subclinical traits and biological phenotypes with SCZ, including thalamic structural abnormalities. In this regard, previous findings have suggested that multivariate volumetric patterns of individual thalamic nuclei discriminate HC from SCZ. These results were obtained using machine learning, which allows case-control classification at the single-subject level. However, machine learning accuracy is usually unsatisfactory possibly due to phenotype heterogeneity. Indeed, a source of misclassification may be related to thalamic structural characteristics of those HC with high schizotypy, which may resemble structural abnormalities of SCZ. We hypothesized that thalamic structural heterogeneity is related to schizotypy, such that high schizotypal burden would implicate misclassification of those HC whose thalamic patterns resemble SCZ abnormalities. METHODS: Following a previous report, we used Random Forests to predict diagnosis in a case-control sample (SCZ = 131, HC = 255) based on thalamic nuclei gray matter volumes estimates. Then, we investigated whether the likelihood to be classified as SCZ (π-SCZ) was associated with schizotypy in 174 HC, evaluated with the Schizotypal Personality Questionnaire. RESULTS: Prediction accuracy was 72.5%. Misclassified HC had higher positive schizotypy scores, which were correlated with π-SCZ. Results were specific to thalamic rather than whole-brain structural features. CONCLUSIONS: These findings strengthen the relevance of thalamic structural abnormalities to SCZ and suggest that multivariate thalamic patterns are correlates of the continuum between schizotypy in HC and the full-blown disease.

NURR1 and ERR1 Modulate the Expression of Genes of a DRD2 Coexpression Network Enriched for Schizophrenia Risk.

Multiple schizophrenia (SCZ) risk loci may be involved in gene co-regulation mechanisms, and analysis of coexpressed gene networks may help to clarify SCZ molecular basis. We have previously identified a dopamine D2 receptor (DRD2) coexpression module enriched for SCZ risk genes and associated with cognitive and neuroimaging phenotypes of SCZ, as well as with response to treatment with antipsychotics. Here we aimed to identify regulatory factors modulating this coexpression module and their relevance to SCZ. We performed motif enrichment analysis to identify transcription factor (TF) binding sites in human promoters of genes coexpressed with DRD2. Then, we measured transcript levels of a group of these genes in primary mouse cortical neurons in basal conditions and upon overexpression and knockdown of predicted TFs. Finally, we analyzed expression levels of these TFs in dorsolateral prefrontal cortex (DLPFC) of SCZ patients. Our in silico analysis revealed enrichment for NURR1 and ERR1 binding sites. In neuronal cultures, the expression of genes either relevant to SCZ risk (Drd2, Gatad2a, Slc28a1, Cnr1) or indexing coexpression in our module (Btg4, Chit1, Osr1, Gpld1) was significantly modified by gain and loss of Nurr1 and Err1. Postmortem DLPFC expression data analysis showed decreased expression levels of NURR1 and ERR1 in patients with SCZ. For NURR1 such decreased expression is associated with treatment with antipsychotics. Our results show that NURR1 and ERR1 modulate the transcription of DRD2 coexpression partners and support the hypothesis that NURR1 is involved in the response to SCZ treatment.SIGNIFICANCE STATEMENT In the present study, we provide in silico and experimental evidence for a role of the TFs NURR1 and ERR1 in modulating the expression pattern of genes coexpressed with DRD2 in human DLPFC. Notably, genetic variations in these genes is associated with SCZ risk and behavioral and neuroimaging phenotypes of the disease, as well as with response to treatment. Furthermore, this study presents novel findings on a possible interplay between D2 receptor-mediated dopamine signaling involved in treatment with antipsychotics and the transcriptional regulation mechanisms exerted by NURR1. Our results suggest that coexpression and co-regulation mechanisms may help to explain some of the complex biology of genetic associations with SCZ.

Publisher Correction: Common brain disorders are associated with heritable patterns of apparent aging of the brain.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The 2019 Schizophrenia International Research Society Conference, 10-14 April, Orlando, Florida: A summary of topics and trends.

The Schizophrenia International Research Society (SIRS) recently held its first North American congress, which took place in Orlando, Florida from 10-14 April 2019. The overall theme of this year's congress was United in Progress - with the aim of cultivating a collaborative effort towards advancing the field of schizophrenia research. Student travel awardees provided reports of the oral sessions and concurrent symposia that took place during the congress. A collection of these reports is summarized and presented below and highlights the main themes and topics that emerged during the congress. In summary, the congress covered a broad range of topics relevant to the field of psychiatry today.

Common brain disorders are associated with heritable patterns of apparent aging of the brain.

Common risk factors for psychiatric and other brain disorders are likely to converge on biological pathways influencing the development and maintenance of brain structure and function across life. Using structural MRI data from 45,615 individuals aged 3-96 years, we demonstrate distinct patterns of apparent brain aging in several brain disorders and reveal genetic pleiotropy between apparent brain aging in healthy individuals and common brain disorders.

Brain-derived neurotrophic factor and schizophrenia.

The brain-derived neurotrophic factor (BDNF) is a secretory growth factor that promotes neuronal proliferation and survival, synaptic plasticity and long-term potentiation in the central nervous system. Brain-derived neurotrophic factor biosynthesis and secretion are chrono-topically regulated processes at the cellular level, accounting for specific localizations and functions. Given its role in regulating brain development and activity, BDNF represents a potentially relevant gene for schizophrenia, and indeed BDNF and its non-synonymous functional variant, rs6265 (C → T, Val → Met) have been widely studied in psychiatric genetics. Human and animal studies have indicated that brain-derived neurotrophic factor is relevant for schizophrenia-related phenotypes, and that: (1) fine-tuned regulation of brain-derived neurotrophic factor secretion and activity is necessary to guarantee brain optimal development and functioning; (2) the Val → Met substitution is associated with impaired activity-dependent secretion of brain-derived neurotrophic factor; (3) disruption of brain-derived neurotrophic factor signaling is associated with altered synaptic plasticity and neurodevelopment. However, genome-wide association studies failed to associate the BDNF locus with schizophrenia, even though a sub-threshold association exists. Here, we will review studies focused on the relationship between the genetic variation of BDNF and schizophrenia, trying to fill the gap between genetic risk per se and insights from molecular biology. A deeper understanding of brain-derived neurotrophic factor biology and of the epigenetic regulation of brain-derived neurotrophic factor and its interactome during development may help clarifying the potential role of this gene in schizophrenia, thus informing development of brain-derived neurotrophic factor-based strategies of prevention and treatment of this disorder.

Prefrontal Coexpression of Schizophrenia Risk Genes Is Associated With Treatment Response in Patients.

BACKGROUND: Gene coexpression networks are relevant to functional and clinical translation of schizophrenia risk genes. We hypothesized that schizophrenia risk genes converge into coexpression pathways that may be associated with gene regulation mechanisms and with response to treatment in patients with schizophrenia. METHODS: We identified gene coexpression networks in two prefrontal cortex postmortem RNA sequencing datasets (n = 688) and replicated them in four more datasets (n = 1295). We identified and replicated (p values < .001) a single module enriched for schizophrenia risk loci (13 risk genes in 10 loci). In silico screening of potential regulators of the schizophrenia risk module via bioinformatic analyses identified two transcription factors and three microRNAs associated with the risk module. To translate postmortem information into clinical phenotypes, we identified polymorphisms predicting coexpression and combined them to obtain an index approximating module coexpression (Polygenic Coexpression Index [PCI]). RESULTS: The PCI-coexpression association was successfully replicated in two independent brain transcriptome datasets (n = 131; p values < .05). Finally, we tested the association between the PCI and short-term treatment response in two independent samples of patients with schizophrenia treated with olanzapine (n = 167). The PCI was associated with treatment response in the positive symptom domain in both clinical cohorts (p values < .05). CONCLUSIONS: In summary, our findings in 1983 samples of human postmortem prefrontal cortex show that coexpression of a set of genes enriched for schizophrenia risk genes is relevant to treatment response. This coexpression pathway may be coregulated by transcription factors and microRNA associated with it.

Thalamic connectivity measured with fMRI is associated with a polygenic index predicting thalamo-prefrontal gene co-expression.

The functional connectivity between thalamic medio-dorsal nucleus (MD) and cortical regions, especially the dorsolateral prefrontal cortex (DLPFC), is implicated in attentional processing and is anomalous in schizophrenia, a brain disease associated with polygenic risk and attentional deficits. However, the molecular and genetic underpinnings of thalamic connectivity anomalies are unclear. Given that gene co-expression across brain areas promotes synchronous interregional activity, our aim was to investigate whether coordinated expression of genes relevant to schizophrenia in MD and DLPFC may reflect thalamic connectivity anomalies in an attention-related network including the DLPFC. With this aim, we identified in datasets of post-mortem prefrontal mRNA expression from healthy controls a gene module with robust overrepresentation of genes with coordinated MD-DLPFC expression and enriched for schizophrenia genes according to the largest genome-wide association study to date. To link this gene cluster with imaging phenotypes, we computed a Polygenic Co-Expression Index (PCI) combining single-nucleotide polymorphisms predicting module co-expression. Finally, we investigated the association between PCI and thalamic functional connectivity during attention through fMRI Independent Component Analysis in 265 healthy participants. We found that PCI was positively associated with connectivity strength of a thalamic region overlapping with the MD within an attention brain circuit. These findings identify a novel association between schizophrenia-related genes and thalamic functional connectivity. Furthermore, they highlight the association between gene expression co-regulation and brain connectivity, such that genes with coordinated MD-DLPFC expression are associated with coordinated activity between the same brain regions. We suggest that gene co-expression is a plausible mechanism underlying biological phenotypes of schizophrenia.

Reproducible grey matter patterns index a multivariate, global alteration of brain structure in schizophrenia and bipolar disorder.

Schizophrenia is a severe mental disorder characterized by numerous subtle changes in brain structure and function. Machine learning allows exploring the utility of combining structural and functional brain magnetic resonance imaging (MRI) measures for diagnostic application, but this approach has been hampered by sample size limitations and lack of differential diagnostic data. Here, we performed a multi-site machine learning analysis to explore brain structural patterns of T1 MRI data in 2668 individuals with schizophrenia, bipolar disorder or attention-deficit/ hyperactivity disorder, and healthy controls. We found reproducible changes of structural parameters in schizophrenia that yielded a classification accuracy of up to 76% and provided discrimination from ADHD, through it lacked specificity against bipolar disorder. The observed changes largely indexed distributed grey matter alterations that could be represented through a combination of several global brain-structural parameters. This multi-site machine learning study identified a brain-structural signature that could reproducibly differentiate schizophrenia patients from controls, but lacked specificity against bipolar disorder. While this currently limits the clinical utility of the identified signature, the present study highlights that the underlying alterations index substantial global grey matter changes in psychotic disorders, reflecting the biological similarity of these conditions, and provide a roadmap for future exploration of brain structural alterations in psychiatric patients.

Genetic Variation of a DRD2 Co-expression Network is Associated with Changes in Prefrontal Function After D2 Receptors Stimulation.

Dopamine D2 receptors (D2Rs) contribute to the inverted U-shaped relationship between dopamine signaling and prefrontal function. Genetic networks from post-mortem human brain revealed 84 partner genes co-expressed with DRD2. Moreover, eight functional single nucleotide polymorphisms combined into a polygenic co-expression index (PCI) predicted co-expression of this DRD2 network and were associated with prefrontal function in humans. Here, we investigated the non-linear association of the PCI with behavioral and Working Memory (WM) related brain response to pharmacological D2Rs stimulation. Fifty healthy volunteers took part in a double-blind, placebo-controlled, functional MRI (fMRI) study with bromocriptine and performed the N-Back task. The PCI by drug interaction was significant on both WM behavioral scores (P = 0.046) and related prefrontal activity (all corrected P < 0.05) using a polynomial PCI model. Non-linear responses under placebo were reversed by bromocriptine administration. fMRI results on placebo were replicated in an independent sample of 50 participants who did not receive drug administration (P = 0.034). These results match earlier evidence in non-human primates and confirm the physiological relevance of this DRD2 co-expression network. Results show that in healthy subjects, different alleles evaluated as an ensemble are associated with non-linear prefrontal responses. Therefore, brain response to a dopaminergic drug may depend on a complex system of allelic patterns associated with DRD2 co-expression.

Transcriptomic context of DRD1 is associated with prefrontal activity and behavior during working memory.

Dopamine D1 receptor (D1R) signaling shapes prefrontal cortex (PFC) activity during working memory (WM). Previous reports found higher WM performance associated with alleles linked to greater expression of the gene coding for D1Rs (DRD1). However, there is no evidence on the relationship between genetic modulation of DRD1 expression in PFC and patterns of prefrontal activity during WM. Furthermore, previous studies have not considered that D1Rs are part of a coregulated molecular environment, which may contribute to D1R-related prefrontal WM processing. Thus, we hypothesized a reciprocal link between a coregulated (i.e., coexpressed) molecular network including DRD1 and PFC activity. To explore this relationship, we used three independent postmortem prefrontal mRNA datasets (total n = 404) to characterize a coexpression network including DRD1 Then, we indexed network coexpression using a measure (polygenic coexpression index-DRD1-PCI) combining the effect of single nucleotide polymorphisms (SNPs) on coexpression. Finally, we associated the DRD1-PCI with WM performance and related brain activity in independent samples of healthy participants (total n = 371). We identified and replicated a coexpression network including DRD1, whose coexpression was correlated with DRD1-PCI. We also found that DRD1-PCI was associated with lower PFC activity and higher WM performance. Behavioral and imaging results were replicated in independent samples. These findings suggest that genetically predicted expression of DRD1 and of its coexpression partners stratifies healthy individuals in terms of WM performance and related prefrontal activity. They also highlight genes and SNPs potentially relevant to pharmacological trials aimed to test cognitive enhancers modulating DRD1 signaling.