High resolution mapping of the tumor microenvironment using integrated single-cell, spatial and in situ analysis.

Single-cell and spatial technologies that profile gene expression across a whole tissue are revolutionizing the resolution of molecular states in clinical samples. Current commercially available technologies provide whole transcriptome single-cell, whole transcriptome spatial, or targeted in situ gene expression analysis. Here, we combine these technologies to explore tissue heterogeneity in large, FFPE human breast cancer sections. This integrative approach allowed us to explore molecular differences that exist between distinct tumor regions and to identify biomarkers involved in the progression towards invasive carcinoma. Further, we study cell neighborhoods and identify rare boundary cells that sit at the critical myoepithelial border confining the spread of malignant cells. Here, we demonstrate that each technology alone provides information about molecular signatures relevant to understanding cancer heterogeneity; however, it is the integration of these technologies that leads to deeper insights, ushering in discoveries that will progress oncology research and the development of diagnostics and therapeutics.

Synapse-Enriched m6A-Modified Malat1 Interacts with the Novel m6A Reader, DPYSL2, and Is Required for Fear-Extinction Memory.

The RNA modification N6-methyladenosine (m6A) regulates the interaction between RNA and various RNA binding proteins within the nucleus and other subcellular compartments and has recently been shown to be involved in experience-dependent plasticity, learning, and memory. Using m6A RNA-sequencing, we have discovered a distinct population of learning-related m6A- modified RNAs at the synapse, which includes the long noncoding RNA metastasis-associated lung adenocarcinoma transcript 1 (Malat1). RNA immunoprecipitation and mass spectrometry revealed 12 new synapse-specific learning-induced m6A readers in the mPFC of male C57/BL6 mice, with m6A-modified Malat1 binding to a subset of these, including CYFIP2 and DPYSL2. In addition, a cell type- and synapse-specific, and state-dependent, reduction of m6A on Malat1 impairs fear-extinction memory; an effect that likely occurs through a disruption in the interaction between Malat1 and DPYSL2 and an associated decrease in dendritic spine formation. These findings highlight the critical role of m6A in regulating the functional state of RNA during the consolidation of fear-extinction memory, and expand the repertoire of experience-dependent m6A readers in the synaptic compartment.SIGNIFICANCE STATEMENT We have discovered that learning-induced m6A-modified RNA (including the long noncoding RNA, Malat1) accumulates in the synaptic compartment. We have identified several new m6A readers that are associated with fear extinction learning and demonstrate a causal relationship between m6A-modified Malat1 and the formation of fear-extinction memory. These findings highlight the role of m6A in regulating the functional state of an RNA during memory formation and expand the repertoire of experience-dependent m6A readers in the synaptic compartment.

Compartment-specific dendritic information processing in striatal cholinergic interneurons is reconfigured by peptide neuromodulation.

Cholinergic interneurons are central hubs of the striatal neuronal network, controlling information processing in a behavioral-state-dependent manner. It remains unknown, however, how such state transitions influence the integrative properties of these neurons. To address this, we made simultaneous somato-dendritic recordings from identified rodent cholinergic interneurons, revealing that action potentials are initiated at dendritic sites because of a dendritic axonal origin. Functionally, this anatomical arrangement ensured that the action potential initiation threshold was lowest at axon-bearing dendritic sites, a privilege efficacy powerfully accentuated at the hyperpolarized membrane potentials achieved in cholinergic interneurons following salient behavioral stimuli. Experimental analysis revealed the voltage-dependent attenuation of the efficacy of non-axon-bearing dendritic excitatory input was mediated by the recruitment of dendritic potassium channels, a regulatory mechanism that, in turn, was controlled by the pharmacological activation of neurokinin receptors. Together, these results indicate that the neuropeptide microenvironment dynamically controls state- and compartment-dependent dendritic information processing in striatal cholinergic interneurons.

Uncovering the spatial landscape of molecular interactions within the tumor microenvironment through latent spaces.

Recent advances in spatial transcriptomics (STs) enable gene expression measurements from a tissue sample while retaining its spatial context. This technology enables unprecedented in situ resolution of the regulatory pathways that underlie the heterogeneity in the tumor as well as the tumor microenvironment (TME). The direct characterization of cellular co-localization with spatial technologies facilities quantification of the molecular changes resulting from direct cell-cell interaction, as it occurs in tumor-immune interactions. We present SpaceMarkers, a bioinformatics algorithm to infer molecular changes from cell-cell interactions from latent space analysis of ST data. We apply this approach to infer the molecular changes from tumor-immune interactions in Visium spatial transcriptomics data of metastasis, invasive and precursor lesions, and immunotherapy treatment. Further transfer learning in matched scRNA-seq data enabled further quantification of the specific cell types in which SpaceMarkers are enriched. Altogether, SpaceMarkers can identify the location and context-specific molecular interactions within the TME from ST data.

A single-cell and spatially resolved atlas of human breast cancers.

Breast cancers are complex cellular ecosystems where heterotypic interactions play central roles in disease progression and response to therapy. However, our knowledge of their cellular composition and organization is limited. Here we present a single-cell and spatially resolved transcriptomics analysis of human breast cancers. We developed a single-cell method of intrinsic subtype classification (SCSubtype) to reveal recurrent neoplastic cell heterogeneity. Immunophenotyping using cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) provides high-resolution immune profiles, including new PD-L1/PD-L2+ macrophage populations associated with clinical outcome. Mesenchymal cells displayed diverse functions and cell-surface protein expression through differentiation within three major lineages. Stromal-immune niches were spatially organized in tumors, offering insights into antitumor immune regulation. Using single-cell signatures, we deconvoluted large breast cancer cohorts to stratify them into nine clusters, termed 'ecotypes', with unique cellular compositions and clinical outcomes. This study provides a comprehensive transcriptional atlas of the cellular architecture of breast cancer.

Spatial transcriptomics at subspot resolution with BayesSpace.

Recent spatial gene expression technologies enable comprehensive measurement of transcriptomic profiles while retaining spatial context. However, existing analysis methods do not address the limited resolution of the technology or use the spatial information efficiently. Here, we introduce BayesSpace, a fully Bayesian statistical method that uses the information from spatial neighborhoods for resolution enhancement of spatial transcriptomic data and for clustering analysis. We benchmark BayesSpace against current methods for spatial and non-spatial clustering and show that it improves identification of distinct intra-tissue transcriptional profiles from samples of the brain, melanoma, invasive ductal carcinoma and ovarian adenocarcinoma. Using immunohistochemistry and an in silico dataset constructed from scRNA-seq data, we show that BayesSpace resolves tissue structure that is not detectable at the original resolution and identifies transcriptional heterogeneity inaccessible to histological analysis. Our results illustrate BayesSpace's utility in facilitating the discovery of biological insights from spatial transcriptomic datasets.

Multi-omic analysis of stroke recurrence in African Americans from the Vitamin Intervention for Stroke Prevention (VISP) clinical trial.

African Americans endure a nearly two-fold greater risk of suffering a stroke and are 2-3 times more likely to die from stroke compared to those of European ancestry. African Americans also have a greater risk of recurrent stroke and vascular events, which are deadlier and more disabling than incident stroke. Stroke is a multifactorial disease with both heritable and environmental risk factors. We conducted an integrative, multi-omic study on 922 plasma metabolites, 473,864 DNA methylation loci, and 556 variants from 50 African American participants of the Vitamin Intervention for Stroke Prevention clinical trial to help elucidate biomarkers contributing to recurrent stroke rates in this high risk population. Sixteen metabolites, including cotinine, N-delta-acetylornithine, and sphingomyelin (d17:1/24:1) were identified in t-tests of recurrent stroke outcome or baseline smoking status. Serum tricosanoyl sphingomyelin (d18:1/23:0) levels were significantly associated with recurrent stroke after adjusting for covariates in Cox Proportional Hazards models. Weighted Gene Co-expression Network Analysis identified moderate correlations between sphingolipid markers and clinical traits including days to recurrent stroke. Integrative analyses between genetic variants in sphingolipid pathway genes identified 29 nominal associations with metabolite levels in a one-way analysis of variance, while epigenomic analyses identified xenobiotics, predominately smoking-associated metabolites and pharmaceutical drugs, associated with methylation profiles. Taken together, our results suggest that metabolites, specifically those associated with sphingolipid metabolism, are potential plasma biomarkers for stroke recurrence in African Americans. Furthermore, genetic variation and DNA methylation may play a role in the regulation of these metabolites.

A general approach for detecting expressed mutations in AML cells using single cell RNA-sequencing.

Virtually all tumors are genetically heterogeneous, containing mutationally-defined subclonal cell populations that often have distinct phenotypes. Single-cell RNA-sequencing has revealed that a variety of tumors are also transcriptionally heterogeneous, but the relationship between expression heterogeneity and subclonal architecture is unclear. Here, we address this question in the context of Acute Myeloid Leukemia (AML) by integrating whole genome sequencing with single-cell RNA-sequencing (using the 10x Genomics Chromium Single Cell 5' Gene Expression workflow). Applying this approach to five cryopreserved AML samples, we identify hundreds to thousands of cells containing tumor-specific mutations in each case, and use the results to distinguish AML cells (including normal-karyotype AML cells) from normal cells, identify expression signatures associated with subclonal mutations, and find cell surface markers that could be used to purify subclones for further study. This integrative approach for connecting genotype to phenotype is broadly applicable to any sample that is phenotypically and genetically heterogeneous.

Quantification of glutathione in the human brain by MR spectroscopy at 3 Tesla: Comparison of PRESS and MEGA-PRESS.

PURPOSE: Glutathione (GSH) is an important intracellular antioxidant in the brain. A number of studies report its measurement by localized 1 H spectroscopy using PRESS and STEAM. This study evaluates the reliability and accuracy of GSH measurements from PRESS at 3 Tesla (T) and compares the results to those obtained with MEGA-PRESS. METHODS: Phantoms containing brain metabolites, identical except for variable GSH concentration between 0 and 24 mM, were scanned using PRESS (echo time (TE) = 35 ms) and MEGA-PRESS (optimized TE = 130 ms) at 3 T. Spectra of the anterior cingulate cortex and occipital cortex in seven healthy volunteers were also acquired. RESULTS: Phantom GSH concentrations from 0 to 3mM were unreliably quantified using PRESS, although at 4 mM and above there was a linear relationship between measured and true concentrations (R2 = 0.99). Using MEGA-PRESS, there was no signal detected at 0 mM GSH, plus a linear relationship (R2 = 0.99) over the full range from 0-24 mM. In brain, concentrations calculated from MEGA-PRESS and PRESS were significantly different in occipital cortex (P < 0.001). Moreover, only MEGA-PRESS reported significant differences in GSH between the two brain regions (P = 0.003). CONCLUSION: Due to uncertainties in GSH quantification raised by the study, the authors conclude that physiological concentrations (<4 mM) of GSH cannot be reliably quantified from PRESS (TE = 35 ms) spectra at 3 T. Magn Reson Med 78:1257-1266, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Glutathione in the human brain: Review of its roles and measurement by magnetic resonance spectroscopy.

We review the transport, synthesis and catabolism of glutathione in the brain as well as its compartmentation and biochemistry in different brain cells. The major reactions involving glutathione are reviewed and the factors limiting its availability in brain cells are discussed. We also describe and critique current methods for measuring glutathione in the human brain using magnetic resonance spectroscopy, and review the literature on glutathione measurements in healthy brains and in neurological, psychiatric, neurodegenerative and neurodevelopmental conditions In summary: Healthy human brain glutathione concentration is ∼1-2 mM, but it varies by brain region, with evidence of gender differences and age effects; in neurological disease glutathione appears reduced in multiple sclerosis, motor neurone disease and epilepsy, while being increased in meningiomas; in psychiatric disease the picture is complex and confounded by methodological differences, regional effects, length of disease and drug-treatment. Both increases and decreases in glutathione have been reported in depression and schizophrenia. In Alzheimer's disease and mild cognitive impairment there is evidence for a decrease in glutathione compared to age-matched healthy controls. Improved methods to measure glutathione in vivo will provide better precision in glutathione determination and help resolve the complex biochemistry of this molecule in health and disease.

Genome-wide meta-analysis of homocysteine and methionine metabolism identifies five one carbon metabolism loci and a novel association of ALDH1L1 with ischemic stroke.

Circulating homocysteine levels (tHcy), a product of the folate one carbon metabolism pathway (FOCM) through the demethylation of methionine, are heritable and are associated with an increased risk of common diseases such as stroke, cardiovascular disease (CVD), cancer and dementia. The FOCM is the sole source of de novo methyl group synthesis, impacting many biological and epigenetic pathways. However, the genetic determinants of elevated tHcy (hyperhomocysteinemia), dysregulation of methionine metabolism and the underlying biological processes remain unclear. We conducted independent genome-wide association studies and a meta-analysis of methionine metabolism, characterized by post-methionine load test tHcy, in 2,710 participants from the Framingham Heart Study (FHS) and 2,100 participants from the Vitamin Intervention for Stroke Prevention (VISP) clinical trial, and then examined the association of the identified loci with incident stroke in FHS. Five genes in the FOCM pathway (GNMT [p = 1.60 × 10(-63)], CBS [p = 3.15 × 10(-26)], CPS1 [p = 9.10 × 10(-13)], ALDH1L1 [p = 7.3 × 10(-13)] and PSPH [p = 1.17 × 10(-16)]) were strongly associated with the difference between pre- and post-methionine load test tHcy levels (ΔPOST). Of these, one variant in the ALDH1L1 locus, rs2364368, was associated with incident ischemic stroke. Promoter analyses reveal genetic and epigenetic differences that may explain a direct effect on GNMT transcription and a downstream affect on methionine metabolism. Additionally, a genetic-score consisting of the five significant loci explains 13% of the variance of ΔPOST in FHS and 6% of the variance in VISP. Association between variants in FOCM genes with ΔPOST suggest novel mechanisms that lead to differences in methionine metabolism, and possibly the epigenome, impacting disease risk. These data emphasize the importance of a concerted effort to understand regulators of one carbon metabolism as potential therapeutic targets.

Cyclic AMP signaling control of action potential firing rate and molecular circadian pacemaking in the suprachiasmatic nucleus.

Circadian pacemaking in suprachiasmatic nucleus (SCN) neurons revolves around transcriptional/posttranslational feedback loops, driven by protein products of "clock" genes. These loops are synchronized and sustained by intercellular signaling, involving vasoactive intestinal peptide (VIP) via its VPAC2 receptor, which positively regulates cAMP synthesis. In turn, SCN cells communicate circadian time to the brain via a daily rhythm in electrophysiological activity. To investigate the mechanisms whereby VIP/VPAC2/cAMP signaling controls SCN molecular and electrical pacemaking, we combined bioluminescent imaging of circadian gene expression and whole-cell electrophysiology in organotypic SCN slices. As a potential direct target of cAMP, we focused on hyperpolarization-activated, cyclic nucleotide-gated (HCN) channels. Mutations of VIP-ergic signaling compromised the SCN molecular pacemaker, diminishing the amplitude and intercellular synchrony of circadian gene expression. These deficits were transiently reversed by elevation of cAMP. Similarly, cellular synchrony in electrical firing rates was lost in SCN slices lacking the VPAC2 receptor for VIP. Whole-cell current-clamp recordings in wild-type (WT) slices revealed voltage responses shaped by the conductance I(h), which is mediated by HCN channel activity. The influence of I(h) on voltage responses showed a modest peak in early circadian day, identifying HCN channels as a putative mediator of cAMP-dependent circadian effects on firing rate. I(h), however, was unaffected by loss of VIP-ergic signaling in VPAC2-null slices, and inhibition of cAMP synthesis had no discernible effect on I(h) but did suppress gene expression and SCN firing rates. Moreover, only sustained but not acute, pharmacological blockade of HCN channels reduced action potential (AP) firing. Thus, our evidence suggests that in the SCN, cAMP-mediated signaling is not a principal regulator of HCN channel function and that HCN is not a determinant of AP firing rate. VIP/cAMP-dependent signaling sustains the SCN molecular oscillator and action potential firing via mechanisms yet to be identified.

Rai1 haploinsufficiency causes reduced Bdnf expression resulting in hyperphagia, obesity and altered fat distribution in mice and humans with no evidence of metabolic syndrome.

Smith-Magenis syndrome (SMS) is a genetic disorder caused by haploinsufficiency of the retinoic acid induced 1 (RAI1) gene. In addition to intellectual disabilities, behavioral abnormalities and sleep disturbances, a majority of children with SMS also have significant early-onset obesity. To study the role of RAI1 in obesity, we investigated the growth and obesity phenotype in a mouse model haploinsufficient for Rai1. Data show that Rai1(+/-) mice are hyperphagic, have an impaired satiety response and have altered abdominal and subcutaneous fat distribution, with Rai1(+/-) female mice having a higher proportion of abdominal fat when compared with wild-type female mice. Expression analyses revealed that Bdnf (brain-derived neurotrophic factor), a gene previously associated with hyperphagia and obesity, is downregulated in the Rai1(+/-) mouse hypothalamus, and reporter studies show that RAI1 directly regulates the expression of BDNF. Even though the Rai1(+/-) mice are significantly obese, serum analyses do not reveal any evidence of metabolic syndrome. Supporting these findings, a caregiver survey revealed that even though a high incidence of abdominal obesity is observed in females with SMS, they did not exhibit a higher incidence of indicators of metabolic syndrome above the general population. We conclude that Rai1 haploinsufficiency represents a single-gene model of obesity with hyperphagia, abnormal fat distribution and altered hypothalamic gene expression associated with satiety, food intake, behavior and obesity. Linking RAI1 and BDNF provides a more thorough understanding of the role of Rai1 in growth and obesity and insight into the complex pathogenicity of obesity, behavior and sex-specific differences in adiposity.

On violations of Le Chatelier's principle for a temperature change in small systems observed for short times.

Le Chatelier's principle states that when a system is disturbed, it will shift its equilibrium to counteract the disturbance. However for a chemical reaction in a small, confined system, the probability of observing it proceed in the opposite direction to that predicted by Le Chatelier's principle, can be significant. This work gives a molecular level proof of Le Chatelier's principle for the case of a temperature change. Moreover, a new, exact mathematical expression is derived that is valid for arbitrary system sizes and gives the relative probability that a single experiment will proceed in the endothermic or exothermic direction, in terms of a microscopic phase function. We show that the average of the time integral of this function is the maximum possible value of the purely irreversible entropy production for the thermal relaxation process. Our result is tested against computer simulations of the unfolding of a polypeptide. We prove that any equilibrium reaction mixture on average responds to a temperature increase by shifting its point of equilibrium in the endothermic direction.

Heart energy metabolism after intestinal ischaemia and reperfusion.

BACKGROUND/PURPOSE: Multiple organ failure subsequent to intestinal ischaemia and reperfusion (I/R) includes cardiac failure, but little is known about heart energy metabolism in this setting. This study investigates the effects of intestinal I/R on heart energy metabolism and evaluates the effects of moderate hypothermia. METHODS: Adult rats underwent intestinal ischaemia for 60 minutes followed by 120 minutes of reperfusion. Animals were maintained at either normothermia (36 degrees to 38 degrees C) or moderate hypothermia (30 degrees to 32 degrees C). In experiment A, 2 groups were studied: (1) sham at normothermia; (2) I/R at normothermia. After death, the heart was removed. Cardiac phosphoenergetics were assessed by 31P magnetic resonance spectroscopy; data are expressed as micromoles per gram. In experiment B, 4 groups were studied: (1) sham at normothermia, (2) I/R at normothermia, (3) sham at hypothermia, (4) I/R at hypothermia. At the end of the experiment, the heart was harvested. The activity of carnitine palmitoyl transferase I (CPT I), an important enzyme in the control of fatty acid oxidation, was measured; data are expressed as nanomoles per minute per unit citrate synthase. Results are expressed as mean +/- SEM. RESULTS: In experiment A, there were no differences between the 2 study groups in cardiac phosphocreatine, inorganic phosphate, adenosine triphosphate (ATP), or in the ratio of inorganic phosphate to ATP. In experiment B, CPT I activity was decreased significantly after I/R at normothermia compared with normothermic sham, but this enzyme inhibition was prevented by hypothermia (3.9 +/- 0.2; v I/R). CONCLUSIONS: These results suggest that although cardiac ATP supply was maintained during intestinal I/R at normothermia, the balance of substrate utilisation was shifted from fatty acid oxidation to carbohydrate utilisation. However, moderate hypothermia modified these changes. The beneficial effect of moderate hypothermia on cardiac metabolism during intestinal I/R has potential clinical application in various surgical conditions.