Postmortem Gene Expression Research Articles

Postmortem communication.

The phenomenon of near death and dying experiences has been both of popular interest and of scientific speculation. However, the reality of mental perception at the point of death is currently a subjective experience and has not been formally evaluated. While postmortem gene expression, even in humans, has been evaluated, restoration of postmortem brain activity has heretofore only been attempted in animal models, at the molecular and cellular levels. Meanwhile, progress has been made to translate brain activity of living humans into speech and images. This paper proposes two inter-related thought experiments. First, assuming progress and refinement of the technology of translating human brain activity into interpretable speech and images, can an objective analysis of death experiences be obtained by utilizing these technologies on dying humans? Second, can human brain function be revived postmortem and, if so, can the relevant technologies be utilized for communication with (recently) deceased individuals? In this paper, these questions are considered and possible implications explored.

Multi‐scale brain attributes contribute to the distribution of diffuse glioma subtypes

AbstractGliomas are primary brain tumors and are among the most malignant types. Adult‐type diffuse gliomas can be classified based on their histological and molecular signatures as IDH‐wildtype glioblastoma, IDH‐mutant astrocytoma, and IDH‐mutant and 1p/19q‐codeleted oligodendroglioma. Recent studies have shown that each subtype of glioma has its own specific distribution pattern. However, the mechanisms underlying the specific distributions of glioma subtypes are not entirely clear despite partial explanations such as cell origin. To investigate the impact of multi‐scale brain attributes on glioma distribution, we constructed cumulative frequency maps for diffuse glioma subtypes based on T1w structural images and evaluated the spatial correlation between tumor frequency and diverse brain attributes, including postmortem gene expression, functional connectivity metrics, cerebral perfusion, glucose metabolism, and neurotransmitter signaling. Regression models were constructed to evaluate the contribution of these factors to the anatomic distribution of different glioma subtypes. Our findings revealed that the three different subtypes of gliomas had distinct distribution patterns, showing spatial preferences toward different brain environmental attributes. Glioblastomas were especially likely to occur in regions enriched with synapse‐related pathways and diverse neurotransmitter receptors. Astrocytomas and oligodendrogliomas preferentially occurred in areas enriched with genes associated with neutrophil‐mediated immune responses. The functional network characteristics and neurotransmitter distribution also contributed to oligodendroglioma distribution. Our results suggest that different brain transcriptomic, neurotransmitter, and connectomic attributes are the factors that determine the specific distributions of glioma subtypes. These findings highlight the importance of bridging diverse scales of biological organization when studying neurological dysfunction.

Genetic influences on brain and cognitive health and their interactions with cardiovascular conditions and depression

Approximately 40% of dementia cases could be prevented or delayed by modifiable risk factors related to lifestyle and environment. These risk factors, such as depression and vascular disease, do not affect all individuals in the same way, likely due to inter-individual differences in genetics. However, the precise nature of how genetic risk profiles interact with modifiable risk factors to affect brain health is poorly understood. Here we combine multiple data resources, including genotyping and postmortem gene expression, to map the genetic landscape of brain structure and identify 367 loci associated with cortical thickness and 13 loci associated with white matter hyperintensities (P < 5×10−8), with several loci also showing a significant association with cognitive function. We show that among 220 unique genetic loci associated with cortical thickness in our genome-wide association studies (GWAS), 95 also showed evidence of interaction with depression or cardiovascular conditions. Polygenic risk scores based on our GWAS of inferior frontal thickness also interacted with hypertension in predicting executive function in the Canadian Longitudinal Study on Aging. These findings advance our understanding of the genetic underpinning of brain structure and show that genetic risk for brain and cognitive health is in part moderated by treatable mid-life factors.

Global-brain functional connectivity related with trait anxiety and its association with neurotransmitters and gene expression profiles

BackgroundNumerous studies have explored the neural correlates of trait anxiety, a predisposing factor for several stress-related disorders. However, the findings from previous studies are inconsistent, which might be due to the limited regions of interest (ROI). A recent approach, named global-brain functional connectivity (GBC), has been demonstrated to address the shortcomings of ROI-based analysis. Furthermore, research on the transcriptome–connectome association has provided an approach to link the microlevel transcriptome profile with the macroscale brain network. In this paper, we aim to explore the neurobiology of trait anxiety with an imaging transcriptomic approach using GBC, biological neurotransmitters, and transcriptome profiles. MethodsUsing a sample of resting-state fMRI data, we investigated trait anxiety-related alteration in GBC. We further used behavioral analysis, spatial correlation analysis, and postmortem gene expression to separately assess the cognitive functions, neurotransmitters, and transcriptional profiles related to alteration in GBC in individuals with trait anxiety. ResultsGBC values in the ventromedial prefrontal cortex and the precuneus were negatively correlated with levels of trait anxiety. This alteration was correlated with behavioral terms including social cognition, emotion, and memory. A strong association was revealed between trait anxiety-related alteration in GBC and neurotransmitters, including dopaminergic, serotonergic, GABAergic, and glutamatergic systems in the ventromedial prefrontal cortex and the precuneus. The transcriptional profiles explained the functional connectivity, with correlated genes enriched in transmembrane signaling. LimitationsSeveral limitations should be taken into account in this research. For example, future research should consider using some different approaches based on dynamic or task-based FC analysis, include more neurotransmitter receptors, additional gene expression data from different samples or more genes related to other stress-related disorders. Meanwhile, it is of great significance to include a larger sample size of individuals with a diagnosis of MDD or other disorders for analysis and comparison and apply stricter multiple-comparison correction and threshold settings in future research. ConclusionsOur research employed multimodal data to investigate GBC in the context of trait anxiety and to establish its associations with neurotransmitters and transcriptome profiles. This approach may improve understanding of the neural mechanism, together with the biological and molecular genetic foundations of GBC in trait anxiety.

In Vivo Characterization of Intracortical Probes with Focused Ion Beam-Etched Nanopatterned Topographies.

(1) Background: Intracortical microelectrodes (IMEs) are an important part of interfacing with the central nervous system (CNS) and recording neural signals. However, recording electrodes have shown a characteristic steady decline in recording performance owing to chronic neuroinflammation. The topography of implanted devices has been explored to mimic the nanoscale three-dimensional architecture of the extracellular matrix. Our previous work used histology to study the implant sites of non-recording probes and showed that a nanoscale topography at the probe surface mitigated the neuroinflammatory response compared to probes with smooth surfaces. Here, we hypothesized that the improvement in the neuroinflammatory response for probes with nanoscale surface topography would extend to improved recording performance. (2) Methods: A novel design modification was implemented on planar silicon-based neural probes by etching nanopatterned grooves (with a 500 nm pitch) into the probe shank. To assess the hypothesis, two groups of rats were implanted with either nanopatterned (n = 6) or smooth control (n = 6) probes, and their recording performance was evaluated over 4 weeks. Postmortem gene expression analysis was performed to compare the neuroinflammatory response from the two groups. (3) Results: Nanopatterned probes demonstrated an increased impedance and noise floor compared to controls. However, the recording performances of the nanopatterned and smooth probes were similar, with active electrode yields for control probes and nanopatterned probes being approximately 50% and 45%, respectively, by 4 weeks post-implantation. Gene expression analysis showed one gene, Sirt1, differentially expressed out of 152 in the panel. (4) Conclusions: this study provides a foundation for investigating novel nanoscale topographies on neural probes.

Functional Connectome Hierarchy in Schizotypy and Its Associations With Expression of Schizophrenia-Related Genes.

Schizotypy has been conceptualized as a continuum of symptoms with marked genetic, neurobiological, and sensory-cognitive overlaps to schizophrenia. Hierarchical organization represents a general organizing principle for both the cortical connectome supporting sensation-to-cognition continuum and gene expression variability across the cortex. However, a mapping of connectome hierarchy to schizotypy remains to be established. Importantly, the underlying changes of the cortical connectome hierarchy that mechanistically link gene expressions to schizotypy are unclear. The present study applied novel connectome gradient on resting-state fMRI data from 1013 healthy young adults to investigate schizotypy-associated sensorimotor-to-transmodal connectome hierarchy and assessed its similarity with the connectome hierarchy of schizophrenia. Furthermore, normative and differential postmortem gene expression data were utilized to examine transcriptional profiles linked to schizotypy-associated connectome hierarchy. We found that schizotypy was associated with a compressed functional connectome hierarchy. Moreover, the pattern of schizotypy-related hierarchy exhibited a positive correlation with the connectome hierarchy observed in schizophrenia. This pattern was closely colocated with the expression of schizophrenia-related genes, with the correlated genes being enriched in transsynaptic, receptor signaling and calcium ion binding. The compressed connectome hierarchy suggests diminished functional system differentiation, providing a novel and holistic system-level basis for various sensory-cognition deficits in schizotypy. Importantly, its linkage with schizophrenia-altered hierarchy and schizophrenia-related gene expression yields new insights into the neurobiological continuum of psychosis. It also provides mechanistic insight into how gene variation may drive alterations in functional hierarchy, mediating biological vulnerability of schizotypy to schizophrenia.

Association of protein distribution and gene expression revealed by positron emission tomography and postmortem gene expression in the dopaminergic system of the human brain.

The topological distribution of dopamine-related proteins is determined by gene transcription and subsequent regulations. Recent research strategies integrating positron emission tomography with a transcriptome atlas have opened new opportunities to understand the influence of regulation after transcription on protein distribution. Previous studies have reported that messenger (m)-RNA expression levels spatially correlate with the density maps of serotonin receptors but not with those of transporters. This discrepancy may be due to differences in regulation after transcription between presynaptic and postsynaptic proteins, which have not been studied in the dopaminergic system. Here, we focused on dopamine D1 and D2/D3 receptors and dopamine transporters and investigated their region-wise relationship between mRNA expression and protein distribution. We examined the region-wise correlation between regional binding potentials of the target region relative to that of non-displaceable tissue (BPND) values of 11C-SCH-23390 and mRNA expression levels of dopamine D1 receptors (D1R); regional BPND values of 11C-FLB-457 and mRNA expression levels of dopamine D2/D3 receptors (D2/D3R); and regional total distribution volume (VT) values of 18F-FE-PE2I and mRNA expression levels of dopamine transporters (DAT) using Spearman's rank correlation. We found significant positive correlations between regional BPND values of 11C-SCH-23390 and the mRNA expression levels of D1R (r = 0.769, p = 0.0021). Similar to D1R, regional BPND values of 11C-FLB-457 positively correlated with the mRNA expression levels of D2R (r = 0.809, p = 0.0151) but not with those of D3R (r = 0.413, p = 0.3095). In contrast to D1R and D2R, no significant correlation between VT values of 18F-FE-PE2I and mRNA expression levels of DAT was observed (r = -0.5934, p = 0.140). We found a region-wise correlation between the mRNA expression levels of dopamine D1 and D2 receptors and their respective protein distributions. However, we found no region-wise correlation between the mRNA expression levels of dopamine transporters and their protein distributions, indicating different regulatory mechanisms for the localization of pre- and postsynaptic proteins. These results provide a broader understanding of the application of the transcriptome atlas to neuroimaging studies of the dopaminergic nervous system.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Connectome gradient dysfunction in major depression and its association with gene expression profiles and treatment outcomes.

Patients with major depressive disorder (MDD) exhibit concurrent deficits in both sensory and higher-order cognitive processing. Connectome studies have suggested a principal primary-to-transmodal gradient in functional brain networks, supporting the spectrum from sensation to cognition. However, whether this gradient structure is disrupted in patients with MDD and how this disruption associates with gene expression profiles and treatment outcome remain unknown. Using a large cohort of resting-state fMRI data from 2227 participants (1148 MDD patients and 1079 healthy controls) recruited at nine sites, we investigated MDD-related alterations in the principal connectome gradient. We further used Neurosynth, postmortem gene expression, and an 8-week antidepressant treatment (20 MDD patients) data to assess the meta-analytic cognitive functions, transcriptional profiles, and treatment outcomes related to MDD gradient alterations, respectively. Relative to the controls, MDD patients exhibited global topographic alterations in the principal primary-to-transmodal gradient, including reduced explanation ratio, gradient range, and gradient variation (Cohen's d = 0.16-0.21), and focal alterations mainly in the primary and transmodal systems (d = 0.18-0.25). These gradient alterations were significantly correlated with meta-analytic terms involving sensory processing and higher-order cognition. The transcriptional profiles explained 53.9% variance of the altered gradient pattern, with the most correlated genes enriched in transsynaptic signaling and calcium ion binding. The baseline gradient maps of patients significantly predicted symptomatic improvement after treatment. These results highlight the connectome gradient dysfunction in MDD and its linkage with gene expression profiles and clinical management, providing insight into the neurobiological underpinnings and potential biomarkers for treatment evaluation in this disorder.

Acute sleep deprivation upregulates serotonin 2A receptors in the frontal cortex of mice via the immediate early gene Egr3.

Serotonin 2A receptors (5-HT2ARs) mediate the hallucinogenic effects of psychedelic drugs and are a key target of the leading class of medications used to treat psychotic disorders. These findings suggest that dysfunction of 5-HT2ARs may contribute to the symptoms of schizophrenia, a mental illness characterized by perceptual and cognitive disturbances. Indeed, numerous studies have found that 5-HT2ARs are reduced in the brains of individuals with schizophrenia. However, the mechanisms that regulate 5-HT2AR expression remain poorly understood. Here we show that a physiologic environmental stimulus, sleep deprivation, significantly upregulates 5-HT2ARs levels in the mouse frontal cortex in as little as 6–8 hours (for mRNA and protein, respectively). This induction requires the activity-dependent immediate early gene transcription factor early growth response 3 (Egr3) as it does not occur in Egr3 deficient (−/−) mice. Using chromatin immunoprecipitation, we show that EGR3 protein binds to the promoter of Htr2a, the gene that encodes the 5-HT2ARs, in the frontal cortex in vivo, and drives expression of in vitro reporter constructs via two EGR3 binding sites in the Htr2a promoter. These results suggest that EGR3 directly regulates Htr2a expression, and 5-HT2ARs levels, in the frontal cortex in response to physiologic stimuli. Analysis of publicly available post-mortem gene expression data revealed that both EGR3 and HTR2A mRNA are reduced in the prefrontal cortex of schizophrenia patients compared to controls. Together these findings suggest a mechanism by which environmental stimuli alter levels of a brain receptor that may mediate the symptoms, and treatment, of mental illness.

A review and meta-analysis of gene expression profiles in suicide

Suicide claims over 800,000 deaths worldwide, making it a serious public health problem. The etiopathophysiology of suicide remains unclear and is highly complex, and postmortem gene expression studies can offer insights into the molecular biological mechanism underlying suicide. In the current study, we conducted a meta-analysis of postmortem brain gene expression in relation to suicide. We identified five gene expression datasets for postmortem orbitofrontal, prefrontal, or dorsolateral prefrontal cortical brain regions from the Gene Expression Omnibus repository. After quality control, the total sample size was 380 (141 suicide deaths and 239 deaths from other causes). We performed the analyses using two meta-analytic approaches. We further performed pathway and cell-set enrichment analyses. We found reduced expression of the KCNJ2 (Potassium Inwardly Rectifying Channel Subfamily J Member 2), A2M (Alpha-2-Macroglobulin), AGT (Angiotensinogen), PMP2 (Peripheral Myelin Protein 2), and VEZF1 (Vascular Endothelial Zinc Finger 1) genes (FDR p<0.05). Our findings support the involvement of astrocytes, stress response, immune system, and microglia in suicide. These findings will require further validation in additional large datasets.

Elucidating cellular contributions to tau-mediated neurodegeneration using drosophila and single-cell transcriptomics.

Alzheimer's Dementia (AD) is neuropathologically defined by aggregations of amyloid-beta plaques and tau neurofibrillary tangles in the brain, with tau pathology being most strongly correlated with cognitive decline. Systems biology approaches profiling post-mortem human AD brain tissue have highlighted involvement of key processes such as innate immunity and translational regulation. However, interpretation of cross-sectional post-mortem gene expression data is confounded by co-morbidities and mixed brain pathologies at autopsy. Moreover, characterizing longitudinal, age-related changes in Alzheimer's dementia (AD) is of particular importance because clinical manifestation of AD is a decades-long process, and age is the strongest known risk factor. We previously completed RNA-sequencing and proteomics of several highly tractable, longitudinal, Drosophila models of tauopathy to characterize tau- and age-specific perturbations. However, cellular heterogeneity confounds interpretation of bulk tissue data. We thus leverage single cell RNA sequencing (scRNA-seq) to explore cell-specific contributions to tau- and age-mediated gene expression. Whole brains of adult Drosophila pan-neuronally expressing a mutant variant of human Tau (R406W) and controls were dissociated for scRNA-seq (10x Genomics Chromium) at 1-, 10-, and 20- days of age. Perturbations in cellular composition and gene expression were quantified. To address cross-species conservation of tau-mediated scRNAseq signatures, we leverage existing human AD brain scRNA-seq data. As expected, we highlight evidence for tau-mediated degeneration of Kenyon cells, a neuronal population in the Drosophila "mushroom body" that is required for learning and memory, along with a concomitant increase in glial cell numbers (e.g. gliosis). Gene set enrichment analysis reveal overrepresentation of acetylcholine receptor regulation pathways in vulnerable mushroom body neurons. Interestingly, regularized regression analyses highlight innate immune activity as predictive of tau-related neuronal reduction. Cluster-specific transcriptional signatures between flies and humans demonstrate strong correlations among inferred cell identities. Our results comprise a powerful single cell transcriptomic resource for studying tau-mediated disruption of gene expression and dynamic age-dependent changes at cellular resolution in a tractable model system that is amenable for high throughput functional validation.

Cell survival and DNA damage repair are promoted in the human blood thanatotranscriptome shortly after death

RNA analysis of post-mortem tissues, or thanatotranscriptomics, has become a topic of interest in forensic science due to the essential information it can provide in forensic investigations. Several studies have previously investigated the effect of death on gene transcription, but it has never been conducted with samples of the same individual. For the first time, a longitudinal mRNA expression analysis study was performed with post-mortem human blood samples from individuals with a known time of death. The results reveal that, after death, two clearly differentiated groups of up- and down-regulated genes can be detected. Pathway analysis suggests active processes that promote cell survival and DNA damage repair, rather than passive degradation, are the source of early post-mortem changes of gene expression in blood. In addition, a generalized linear model with an elastic net restriction predicted post-mortem interval with a root mean square error of 4.75 h. In conclusion, we demonstrate that post-mortem gene expression data can be used as biomarkers to estimate the post-mortem interval though further validation using independent sample sets is required before use in forensic casework.

Evaluation of time passed since death by examination of oxidative stress markers, histopathological, and molecular changes of major organs in male albino rats.

Recent biochemical, metabolic, and molecular profiles of various body fluids showed more accurate correlation to the postmortem interval than the traditional physical examination. Our study aimed to evaluate time passed since death in relation to oxidative stress markers, HMGB1 genetic expression, histopathological examination, and BCL2 immunohistochemical analysis in major organs (heart, kidney, and testis). Forty-two adult male rats were included and randomly divided into seven equal groups. After sacrification, the rodents were kept at room temperature and major organs were obtained at 0, 12, 24, 48, 72, 96, and 120h. Malonaldehyde (MDA), superoxide dismutase (SOD), reduced glutathione (GSH) tissue levels, High mobility group box 1 protein (HMGB1) gene expression, histopathological, and B cell lymphoma 2 (BCL2) immunohistochemical expressions were analyzed. Postmortem interval was correlated to different tissue levels of MDA, SOD, and GSH. HMGB1 showed enhanced postmortem gene expression with a peak at 48h after death. Obvious time-dependent histopathological changes were observed in all the examined organs. Dilated spaces, extravasation, and fragmentation scores in heart specimens were higher at 96 and 120h compared with the other groups. Renal changes in the form of shrunken glomeruli, loss of tubular epithelium, and hyalinization and testicular findings in the form of epithelial detachment, vacuolation, and loss of sperms started at 72h postmortem. BCL2 expression began to decrease 24h and became negative at 96h after death. In conclusion, HMGB1 gene expression can be used for estimation of time passed since death as it shows time-dependent changes in the form of a progressive increase with a peak at 48h then it begins to decline. Oxidants and antioxidants are correlated to PMI until 120h after death. Histopathological changes in the heart, kidney, and testis are also time-dependent until the 5th day after death. BCL2 immunohistochemical expression begins to decline 24h until 96h after death when it becomes negative.

NIMG-13. GENETIC, CELLULAR, AND CONNECTOMIC CHARACTERIZATION OF THE ADULT HUMAN BRAIN REGIONS COMMONLY PLAGUED BY GLIOMA

Abstract For decades, it has been known that gliomas follow a nonrandom spatial distribution, appearing more often in some brain regions (e.g. the insula) compared to others (e.g. the occipital lobe). A better understanding of the glioma localization patterns could lend clues to the origins of these types of tumors, and consequently inform treatment targets. Following hypotheses derived from prior research into neuropsychiatric disease and cancer, gliomas may be expected to localize to brain regions characterized by functional hubness, stem-like cells, and transcription of genetic drivers of gliomagenesis. We combined neuroimaging data from 335 adult patients with high- and low-grade glioma to form a replicable tumor frequency map. Using this map, we demonstrate that glioma frequency is elevated in association cortex and correlated with multiple graph-theoretical metrics of high functional connectedness. Brain regions populated with putative cells-of-origin for glioma, neural stem cells and oligodendrocyte precursor cells, exhibited a high glioma frequency. Leveraging a human brain atlas of post-mortem gene expression, we found that gliomas were localized to brain regions enriched with expression of genes associated with chromatin organization and synaptic signaling. A set of glioma proto-oncogenes was enriched among the transcriptomic correlates of glioma distribution. Finally, a regression model incorporating connectomic, cellular, and genetic factors explained 58% of the variance in glioma frequency. These results add to previous literature reporting the vulnerability of hub regions to neurological disease, as well as provide support for cancer stem cell theories of glioma. Our findings illustrate how factors of diverse scale, from genetic to connectomic, can independently influence the anatomic localization of brain dysfunction.

Genetic, cellular, and connectomic characterization of the brain regions commonly plagued by glioma.

For decades, it has been known that gliomas follow a non-random spatial distribution, appearing more often in some brain regions (e.g. the insula) compared to others (e.g. the occipital lobe). A better understanding of the localization patterns of gliomas could provide clues to the origins of these types of tumours, and consequently inform treatment targets. Following hypotheses derived from prior research into neuropsychiatric disease and cancer, gliomas may be expected to localize to brain regions characterized by functional hubness, stem-like cells, and transcription of genetic drivers of gliomagenesis. We combined neuroimaging data from 335 adult patients with high- and low-grade glioma to form a replicable tumour frequency map. Using this map, we demonstrated that glioma frequency is elevated in association cortex and correlated with multiple graph-theoretical metrics of high functional connectedness. Brain regions populated with putative cells of origin for glioma, neural stem cells and oligodendrocyte precursor cells, exhibited a high glioma frequency. Leveraging a human brain atlas of post-mortem gene expression, we found that gliomas were localized to brain regions enriched with expression of genes associated with chromatin organization and synaptic signalling. A set of glioma proto-oncogenes was enriched among the transcriptomic correlates of glioma distribution. Finally, a regression model incorporating connectomic, cellular, and genetic factors explained 58% of the variance in glioma frequency. These results add to previous literature reporting the vulnerability of hub regions to neurological disease, as well as provide support for cancer stem cell theories of glioma. Our findings illustrate how factors of diverse scale, from genetic to connectomic, can independently influence the anatomic localization of brain dysfunction.

Integration of postmortem amygdala expression profiling, GWAS, and functional cell culture assays: neuroticism-associated synaptic vesicle glycoprotein 2A (SV2A) gene is regulated by miR-133a and miR-218

Recent genome-wide studies have begun to identify gene variants, expression profiles, and regulators associated with neuroticism, anxiety disorders, and depression. We conducted a set of experimental cell culture studies of gene regulation by micro RNAs (miRNAs), based on genome-wide transcriptome, proteome, and miRNA expression data from twenty postmortem samples of lateral amygdala from donors with known neuroticism scores. Using Ingenuity Pathway Analysis and TargetScan, we identified a list of mRNA–protein–miRNA sets whose expression patterns were consistent with miRNA-based translational repression, as a function of trait anxiety. Here, we focused on one gene from that list, which is of particular translational significance in Psychiatry: synaptic vesicle glycoprotein 2A (SV2A) is the binding site of the anticonvulsant drug levetiracetam ((S)-α-Ethyl-2-oxo-1-pyrrolidineacetamide), which has shown promise in anxiety disorder treatments. We confirmed that SV2A is associated with neuroticism or anxiety using an original GWAS of a community cohort (N = 1,706), and cross-referencing a published GWAS of multiple cohorts (Ns ranging from 340,569 to 390,278). Postmortem amygdala expression profiling implicated three putative regulatory miRNAs to target SV2A: miR-133a, miR-138, and miR-218. Moving from association to experimental causal testing in cell culture, we used a luciferase assay to demonstrate that miR-133a and miR-218, but not miR-138, significantly decreased relative luciferase activity from the SV2A dual-luciferase construct. In human neuroblastoma cells, transfection with miR-133a and miR-218 reduced both endogenous SV2A mRNA and protein levels, confirming miRNA targeting of the SV2A gene. This study illustrates the utility of combining postmortem gene expression data with GWAS to guide experimental cell culture assays examining gene regulatory mechanisms that may contribute to complex human traits. Identifying specific molecular mechanisms of gene regulation may be useful for future clinical applications in anxiety disorders or other forms of psychopathology.

Genetic Architecture and Molecular Neuropathology of Human Cocaine Addiction.

We integrated genomic and bioinformatic analyses, using data from the largest genome-wide association study of cocaine dependence (CD; n = 6546; 82.37% with CD; 57.39% male) and the largest postmortem gene-expression sample of individuals with cocaine use disorder (CUD; n = 36; 51.35% with CUD; 100% male). Our genome-wide analyses identified one novel gene (NDUFB9) associated with the genetic predisposition to CD in African-Americans. The genetic architecture of CD was similar across ancestries. Individual genes associated with CD demonstrated modest overlap across European-Americans and African-Americans, but the genetic liability for CD converged on many similar tissue types (brain, heart, blood, liver) across ancestries. In a separate sample, we investigated the neuronal gene expression associated with CUD by using RNA sequencing of dorsal-lateral prefrontal cortex neurons. We identified 133 genes differentially expressed between CUD case patients and cocaine-free control subjects, including previously implicated candidates for cocaine use/addiction (FOSB, ARC, KCNJ9/GIRK3, NR4A2, JUNB, and MECP2). Differential expression analyses significantly correlated across European-Americans and African-Americans. While genes significantly associated with CD via genome-wide methods were not differentially expressed, two of these genes (NDUFB9 and C1qL2) were part of a robust gene coexpression network associated with CUD involved in neurotransmission (GABA, acetylcholine, serotonin, and dopamine) and drug addiction. We then used a "guilt-by-association" approach to unravel the biological relevance of NDUFB9 and C1qL2 in the context of CD. In sum, our study furthers the understanding of the genetic architecture and molecular neuropathology of human cocaine addiction and provides a framework for translating biological meaning into otherwise obscure genome-wide associations.SIGNIFICANCE STATEMENT Our study further clarifies the genetic and neurobiological contributions to cocaine addiction, provides a rapid approach for generating testable hypotheses for specific candidates identified by genome-wide research, and investigates the cross-ancestral biological contributions to cocaine use disorder/dependence for individuals of European-American and African-American ancestries.

Identification of cadaveric liver tissues using thanatotranscriptome biomarkers

Thanatotranscriptome studies involve the examination of mRNA transcript abundance and gene expression patterns in the internal organs of deceased humans. Postmortem gene expression is indicative of the cellular status of a corpse at the time of death, a portion of which may represent a cascade of molecular events occasioned by death. Specific gene biomarkers identify perceptible transcriptional changes induced by stochastic responses to the cessation of biological functions. Transcriptome analyses of postmortem mRNA from a tissue fragment may determine unique molecular identifiers for specific organs and demonstrate unique patterns of gene expression that can provide essential contextual anatomical information. We evaluated the impact of targeted transcriptome analysis using RNA sequencing to reveal global changes in postmortem gene expression in liver tissues from 27 Italian and United States corpses: 3.5-hour-old to 37-day-old. We found that our single blind study using eight liver tissue-specific gene biomarkers (e.g. AMBP and AHSG) is highly specific, with autopsy-derived organ samples correctly identified as tissues originating from postmortem livers. The results demonstrate that 98–100% of sequencing reads were mapped to these liver biomarkers. Our findings indicate that gene expression signatures of mRNA exposed up to 37 days of autolysis, can be used to validate the putative identity of tissue fragments.

SCN11A mRNA levels in female bipolar disorder PBMCs as tentative biomarker for distinct patient sub-phenotypes.

Bipolar disorder (BD) is a complex neuropsychiatric disorder characterized by recurrent mania and depression episodes and requiring lifelong treatment with mood stabilizing drugs. Several lines of evidence, including with BD patient iPSC-derived neurons, suggest that neuronal hyperexcitability may underlie the key clinical symptoms of BD. Indeed, higher mRNA levels of SCN11A, coding for the voltage-gated sodium channel NaV 1.9 implicated in nociception, were detected in iPSC-derived neurons from BD patients, and were normalized by in vitro lithium. Here we studied SCN11A expression in peripheral blood mononuclear cells (PBMCs) from well-phenotyped female BD patients and controls and evaluated their association with several clinical sub-phenotypes. We observed higher mRNA levels of SCN11A in PBMCs from female BD patients with no records of alcohol dependence (p = .0050), no records of psychosis (p = .0097), or no records of suicide attempts (p = .0409). A trend was observed for higher SCN11A expression (FD = 1.91; p = .052) in BD PBMCs compared with controls. Datamining of published postmortem gene expression datasets indicated higher SCN11A expression in dorsolateral prefrontal cortex and orbitofrontal cortex tissues from BD patients compared with controls. Higher phenotype-associated expression levels in PBMC from BD patients were also observed for ID2 (alcohol dependence, suicide attempts) and HDGFRP3 (seasonal BD pattern). Our findings suggest that higher PBMC SCN11A expression levels may be associated with certain behavioral BD sub-phenotypes, including lack of alcohol dependence and psychosis, among BD patients. The NaV 1.9 voltage-gated sodium channel thus deserves consideration as a tentative phenotype modifier in BD.