Search

Abstract In brain, the striatum is a heterogenous region involved in reward and goal-directed behaviors. Striatal dysfunction is linked to psychiatric disorders, including opioid use disorder (OUD). Striatal subregions are divided based on neuroanatomy, each with unique roles in OUD. In OUD, the dorsal striatum is involved in altered reward processing, formation of habits, and development of negative affect during withdrawal. Using single nuclei RNA-sequencing, we identified both canonical (e.g., dopamine receptor subtype) and less abundant cell populations (e.g., interneurons) in human dorsal striatum. Pathways related to neurodegeneration, interferon response, and DNA damage were significantly enriched in striatal neurons of individuals with OUD. DNA damage markers were also elevated in striatal neurons of opioid-exposed rhesus macaques. Sex-specific molecular differences in glial cell subtypes associated with chronic stress were found in OUD, particularly female individuals. Together, we describe different cell types in human dorsal striatum and identify cell type-specific alterations in OUD.

Abstract Antipsychotic (AP)-naive first-episode psychosis (FEP) patients display early dysglycemia, including insulin resistance and prediabetes. Metabolic dysregulation may therefore be intrinsic to psychosis spectrum disorders (PSDs), independent of the metabolic effects of APs. However, the potential biological pathways that overlap between PSDs and dysglycemic states remain to be identified. Using meta-analytic approaches of transcriptomic datasets, we investigated whether AP-naive FEP patients share overlapping gene expression signatures with non-psychiatrically ill early dysglycemia individuals. We meta-analyzed peripheral transcriptomic datasets of AP-naive FEP patients and non-psychiatrically ill early dysglycemia subjects to identify common gene expression signatures. Common signatures underwent pathway enrichment analysis and were then used to identify potential new pharmacological compounds via Integrative Library of Integrated Network-Based Cellular Signatures (iLINCS). Our search results yielded 5 AP-naive FEP studies and 4 early dysglycemia studies which met inclusion criteria. We discovered that AP-naive FEP and non-psychiatrically ill subjects exhibiting early dysglycemia shared 221 common signatures, which were enriched for pathways related to endoplasmic reticulum stress and abnormal brain energetics. Nine FDA-approved drugs were identified as potential drug treatments, of which the antidiabetic metformin, the first-line treatment for type 2 diabetes, has evidence to attenuate metabolic dysfunction in PSDs. Taken together, our findings support shared gene expression changes and biological pathways associating PSDs with dysglycemic disorders. These data suggest that the pathobiology of PSDs overlaps and potentially contributes to dysglycemia. Finally, we find that metformin may be a potential treatment for early metabolic dysfunction intrinsic to PSDs.

MicroRNAs (miRNAs) bind on the 3′ untranslated region (3′UTR) of messenger RNAs (mRNAs) and regulate mRNA expression in physiological and pathological conditions, including cancer. Thus, studies have identified miRNAs as potential biomarkers by correlating the miRNA expression with the expression of important mRNAs and/or clinical outcomes in cancers. However, tumors undergo pervasive 3′UTR shortening/lengthening events through alternative polyadenylation (APA), which varies the number of miRNA target sites in mRNA, raising the number of miRNA target sites (numTS) as another important regulatory axis of the miRNA binding effects. In this study, we developed the first statistical method, BIOMATA-APA, to identify predictive miRNAs based on numTS features. Running BIOMATA-APA on The Cancer Genome Atlas (TCGA) and independent cohort data both with immunotherapy and no immunotherapy, we demonstrated for the first time that the numTS feature 1) distinguishes different cancer types, 2) predicts tumor proliferation and immune infiltration status, 3) explains more variation in the proportion of tumor-infiltrating immune cells, 4) predicts response to immune checkpoint blockade (ICB) therapy, and 5) adds prognostic power beyond clinical and miRNA expression. To the best of our knowledge, this is the first pan-cancer study to systematically demonstrate numTS as a novel type of biomarker representing the miRNA binding effects underlying tumorigenesis and pave the way to incorporate miRNA target sites for miRNA biomarker identification. Another advantage of examining the miRNA binding effect using numTS is that it requires only RNA-Seq data, not miRNAs, thus resulting in high power in the miRNA biomarker identification.

Mutations of ESR1, the gene encoding the estrogen receptor alpha, are associated with acquired resistance to therapy in luminalbreast cancer. Here the authors show that ESR1 mutant tumors gain basal-like features with increased expression of basal cytokeratines and immune activation.

Hepatocellular carcinoma (HCC) is one of the most lethal human cancers. Liver transplantation has been an effective approach to treat liver cancer. However, significant numbers of patients with HCC experience cancer recurrence, and the selection of suitable candidates for liver transplant remains a challenge. We developed a model to predict the likelihood of HCC recurrence after liver transplantation based on transcriptome and whole‐exome sequencing analyses. We used a training cohort and a subsequent testing cohort based on liver transplantation performed before or after the first half of 2012. We found that the combination of transcriptome and mutation pathway analyses using a random forest machine learning correctly predicted HCC recurrence in 86.8% of the training set. The same algorithm yielded a correct prediction of HCC recurrence of 76.9% in the testing set. When the cohorts were combined, the prediction rate reached 84.4% in the leave‐one‐out cross‐validation analysis. When the transcriptome analysis was combined with Milan criteria using the k‐top scoring pairs (k‐TSP) method, the testing cohort prediction rate improved to 80.8%, whereas the training cohort and the combined cohort prediction rates were 79% and 84.4%, respectively. Application of the transcriptome/mutation pathways RF model on eight tumor nodules from 3 patients with HCC yielded 8/8 consistency, suggesting a robust prediction despite the heterogeneity of HCC. Conclusion: The genome prediction model may hold promise as an alternative in selecting patients with HCC for liver transplant.

Abstract Severe and persistent disruptions to sleep and circadian rhythms are common in people with opioid use disorder (OUD). Preclinical evidence suggests altered molecular rhythms in the brain modulate opioid reward and relapse. However, whether molecular rhythms are disrupted in the brains of people with OUD remained an open question, critical to understanding the role of circadian rhythms in opioid addiction. Using subjects’ times of death as a marker of time of day, we investigated transcriptional rhythms in the brains of subjects with OUD compared to unaffected comparison subjects. We discovered rhythmic transcripts in both the dorsolateral prefrontal cortex (DLPFC) and nucleus accumbens (NAc), key brain areas involved in OUD, that were largely distinct between OUD and unaffected subjects. Fewer rhythmic transcripts were identified in DLPFC of subjects with OUD compared to unaffected subjects, whereas in the NAc, nearly double the number of rhythmic transcripts was identified in subjects with OUD. In NAc of subjects with OUD, rhythmic transcripts peaked either in the evening or near sunrise, and were associated with an opioid, dopamine, and GABAergic neurotransmission. Associations with altered neurotransmission in NAc were further supported by co-expression network analysis which identified OUD-specific modules enriched for transcripts involved in dopamine, GABA, and glutamatergic synaptic functions. Additionally, rhythmic transcripts in DLPFC and NAc of subjects with OUD were enriched for genomic loci associated with sleep-related GWAS traits, including sleep duration and insomnia. Collectively, our findings connect transcriptional rhythm changes in opioidergic, dopaminergic, GABAergic signaling in the human brain to sleep-related traits in opioid addiction.

Twelve-hour (12 h) ultradian rhythms are a well-known phenomenon in coastal marine organisms. While 12 h cycles are observed in human behavior and physiology, no study has measured 12 h rhythms in the human brain. Here, we identify 12 h rhythms in transcripts that either peak at sleep/wake transitions (approximately 9 AM/PM) or static times (approximately 3 PM/AM) in the dorsolateral prefrontal cortex, a region involved in cognition. Subjects with schizophrenia (SZ) lose 12 h rhythms in genes associated with the unfolded protein response and neuronal structural maintenance. Moreover, genes involved in mitochondrial function and protein translation, which normally peak at sleep/wake transitions, peak instead at static times in SZ, suggesting suboptimal timing of these essential processes.

Microglia are resident macrophages of the brain, performing roles related to brain homeostasis, including modulation of synapses, trophic support, phagocytosis of apoptotic cells and debris, as well as brain protection and repair. Studies assessing morphological and transcriptional features of microglia found regional differences as well as sex differences in some investigated brain regions. However, markers used to isolate microglia in many previous studies are not expressed exclusively by microglia or cannot be used to identify and isolate microglia in all contexts. Here, fluorescent activated cell sorting was used to isolate cells expressing the microglia-specific marker TMEM119 from prefrontal cortex (PFC), striatum, and midbrain in mice. RNA-sequencing was used to assess the transcriptional profile of microglia, focusing on brain region and sex differences. We found striking brain region differences in microglia-specific transcript expression. Most notable was the distinct transcriptional profile of midbrain microglia, with enrichment for pathways related to immune function; these midbrain microglia exhibited a profile similar to disease-associated or immune-surveillant microglia. Transcripts more highly expressed in PFC isolated microglia were enriched for synapse-related pathways while microglia isolated from the striatum were enriched for pathways related to microtubule polymerization. We also found evidence for a gradient of expression of microglia-specific transcripts across the rostral-to-caudal axes of the brain, with microglia extracted from the striatum exhibiting a transcriptional profile intermediate between that of the PFC and midbrain. We also found sex differences in expression of microglia-specific transcripts in all 3 brain regions, with many selenium-related transcripts more highly expressed in females across brain regions. These results suggest that the transcriptional profile of microglia varies between brain regions under homeostatic conditions, suggesting that microglia perform diverse roles in different brain regions and even based on sex.

Abstract The interplay between the immune system and tumor progression is well recognized. However, current human breast cancer immunophenotyping studies are mostly focused on primary tumors with metastatic breast cancer lesions remaining largely understudied. To address this gap, we examined exome-capture RNA sequencing data from 50 primary breast tumors (PBTs) and their patient-matched metastatic tumors (METs) in brain, ovary, bone and gastrointestinal tract. We used gene expression signatures as surrogates for tumor infiltrating lymphocytes (TILs) and compared TIL patterns in PBTs and METs. Enrichment analysis and deconvolution methods both revealed that METs had a significantly lower abundance of total immune cells, including CD8+ T cells, regulatory T cells and dendritic cells. An exception was M2-like macrophages, which were significantly higher in METs across the organ sites examined. Multiplex immunohistochemistry results were consistent with data from the in-silico analysis and showed increased macrophages in METs. We confirmed the finding of a significant reduction in immune cells in brain METs (BRMs) by pathologic assessment of TILs in a set of 49 patient-matched pairs of PBT/BRMs. These findings indicate that METs have an overall lower infiltration of immune cells relative to their matched PBTs, possibly due to immune escape. RNAseq analysis suggests that the relative levels of M2-like macrophages are increased in METs, and their potential role in promoting breast cancer metastasis warrants further study.

Abstract Recent advances in large-scale gene expression profiling necessitate concurrent development of biostatistical approaches to reveal meaningful biological relationships. Most analyses rely on significance thresholds for identifying differentially expressed genes. We use an approach to compare gene expression datasets using ‘threshold-free’ comparisons. Significance cut-offs to identify genes shared between datasets may be too stringent and may miss concordant patterns of gene expression with potential biological relevance. A threshold-free approach gaining popularity in several research areas, including neuroscience, is Rank–Rank Hypergeometric Overlap (RRHO). Genes are ranked by their p-value and effect size direction, and ranked lists are compared to identify significantly overlapping genes across a continuous significance gradient rather than at a single arbitrary cut-off. We have updated the previous RRHO analysis by accurately detecting overlap of genes changed in the same and opposite directions between two datasets. Here, we use simulated and real data to show the drawbacks of the previous algorithm as well as the utility of our new algorithm. For example, we show the power of detecting discordant transcriptional patterns in the postmortem brain of subjects with psychiatric disorders. The new R package, RRHO2, offers a new, more intuitive visualization of concordant and discordant gene overlap.

Abstract Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer following invasive ductal carcinoma (IDC). ILC differs from IDC in a number of histological and clinical features, such as single strand growth, difficulty in detection, and frequent late recurrences. To understand the molecular pathways involved in the clinical characteristics of ILC, we compared the gene expression profiles of luminal A ILC and luminal A IDC using data from TCGA and utilized samples from METABRIC as a validation data set. Top pathways that were significantly enriched in ILC were related to immune response. ILC exhibited a higher activity of almost all types of immune cells based on cell type-specific signatures compared to IDC. Conversely, pathways that were less enriched in ILC were related to protein translation and metabolism, which we functionally validated in cell lines. The higher immune activity uncovered in our study highlights the currently unexplored potential of a response to immunotherapy in a subset of patients with ILC. Furthermore, the lower rates of protein translation and metabolism - known features of tumor dormancy - may play a role in the late recurrences of ILC and lower detection rate in mammography and PET scanning.

Genome-wide expression profiling of the human brain has revealed genes that are differentially expressed across the lifespan. Characterizing these genes adds to our understanding of both normal functions and pathological conditions. Additionally, the specific cell-types that contribute to the motor, sensory and cognitive declines during aging are unclear. Here we test if age-related genes show higher expression in specific neural cell types. Our study leverages data from two sources of murine single-cell expression data and two sources of age-associations from large gene expression studies of postmortem human brain. We used nonparametric gene set analysis to test for age-related enrichment of genes associated with specific cell-types; we also restricted our analyses to specific gene ontology groups. Our analyses focused on a primary pair of single-cell expression data from the mouse visual cortex and age-related human post-mortem gene expression information from the orbitofrontal cortex. Additional pairings that used data from the hippocampus, prefrontal cortex, somatosensory cortex and blood were used to validate and test specificity of our findings. We found robust age-related up-regulation of genes that are highly expressed in oligodendrocytes and astrocytes, while genes highly expressed in layer 2/3 glutamatergic neurons were down-regulated across age. Genes not specific to any neural cell type were also down-regulated, possibly due to the bulk tissue source of the age-related genes. A gene ontology-driven dissection of the cell-type enriched genes highlighted the strong down-regulation of genes involved in synaptic transmission and cell-cell signaling in the Somatostatin (Sst) neuron subtype that expresses the cyclin dependent kinase 6 (Cdk6) and in the vasoactive intestinal peptide (Vip) neuron subtype expressing myosin binding protein C, slow type (Mybpc1). These findings provide new insights into cell specific susceptibility to normal aging, and suggest age-related synaptic changes in specific inhibitory neuron subtypes.

Black versus white older Americans are more likely to experience frailty, a condition associated with adverse health outcomes. To reduce racial disparities in health, a complete understanding of the pathophysiology of frailty is needed. Metabolomics may further our understanding by characterizing differences in the body during a vigorous versus frail state. We sought to identify metabolites and biological pathways associated with vigor to frailty among 287 black men ages 70−81 from the Health, Aging, and Body Composition study. Using liquid chromatography-mass spectrometry, 350 metabolites were measured in overnight-fasting plasma. The Scale of Aging Vigor in Epidemiology (SAVE) measured vigor to frailty based on weight change, strength, energy, gait speed, and physical activity. Thirty-seven metabolites correlated with SAVE scores (p < 0.05), while adjusting for age and site. Fourteen metabolites remained significant after multiple comparisons adjustment (false discovery rate < 0.30). Lower values of tryptophan, methionine, tyrosine, asparagine, C14:0 sphingomyelin, and 1-methylnicotinamide, and higher values of glucoronate, N-carbamoyl-beta-alanine, isocitrate, creatinine, C4-OH carnitine, cystathionine, hydroxyphenylacetate, and putrescine were associated with frailer SAVE scores. Pathway analyses identified nitrogen metabolism, aminoacyl-tRNA biosynthesis, and the citric acid cycle. Future studies need to confirm these SAVE-associated metabolites and pathways that may indicate novel mechanisms involved in the frailty syndrome.

Microarray technology has been widely applied to the analysis of many malignancies, however, integrative analyses across multiple studies are rarely investigated. In this study we performed a meta-analysis on the expression profiles of four published studies analyzing organ donor, benign tissues adjacent to tumor and tumor tissues from liver, prostate, lung and bladder samples. We identified 99 distinct multi-cancer biomarkers in the comparison of all three tissues in liver and prostate and 44 in the comparison of normal versus tumor in liver, prostate and lung. The bladder samples appeared to have a different list of biomarkers from the other three cancer types. The identified multi-cancer biomarkers achieved high accuracy similar to using whole genome in the within-cancer-type prediction. They also performed superior than the one using whole genome in inter-cancer-type prediction. To test the validity of the multi-cancer biomarkers, 23 independent prostate cancer samples were evaluated and 96% accuracy was achieved in inter-study prediction from the original prostate, liver and lung cancer data sets respectively. The result suggests that the compact lists of multi-cancer biomarkers are important in cancer development and represent the common signatures of malignancies of multiple cancer types. Pathway analysis revealed important tumorogenesis functional categories.