Aging compromises human islet beta cell function and identity by decreasing transcription factor activity and inducing ER stress.

Pancreatic islet beta cells are essential for maintaining glucose homeostasis. To understand the impact of aging on beta cells, we performed meta-analysis of single-cell RNA sequencing datasets, transcription factor (TF) regulon analysis, high-resolution confocal microscopy, and measured insulin secretion from nondiabetic donors spanning most of the human life span. This revealed the range of molecular and functional changes that occur during beta cell aging, including the transcriptional deregulation that associates with cellular immaturity and reorganization of beta cell TF networks, increased gene transcription rates, and reduced glucose-stimulated insulin release. These alterations associate with activation of endoplasmic reticulum (ER) stress and autophagy pathways. We propose that a chronic state of ER stress undermines old beta cell structure function to increase the risk of beta cell failure and type 2 diabetes onset as humans age.

Activity-dependent modulation of synapse-regulating genes in astrocytes.

Astrocytes regulate the formation and function of neuronal synapses via multiple signals; however, what controls regional and temporal expression of these signals during development is unknown. We determined the expression profile of astrocyte synapse-regulating genes in the developing mouse visual cortex, identifying astrocyte signals that show differential temporal and layer-enriched expression. These patterns are not intrinsic to astrocytes, but regulated by visually evoked neuronal activity, as they are absent in mice lacking glutamate release from thalamocortical terminals. Consequently, synapses remain immature. Expression of synapse-regulating genes and synaptic development is also altered when astrocyte signaling is blunted by diminishing calcium release from astrocyte stores. Single-nucleus RNA sequencing identified groups of astrocytic genes regulated by neuronal and astrocyte activity, and a cassette of genes that show layer-specific enrichment. Thus, the development of cortical circuits requires coordinated signaling between astrocytes and neurons, highlighting astrocytes as a target to manipulate in neurodevelopmental disorders.

Discovering single nucleotide variants and indels from bulk and single-cell ATAC-seq.

Genetic variants and de novo mutations in regulatory regions of the genome are typically discovered by whole-genome sequencing (WGS), however WGS is expensive and most WGS reads come from non-regulatory regions. The Assay for Transposase-Accessible Chromatin (ATAC-seq) generates reads from regulatory sequences and could potentially be used as a low-cost 'capture' method for regulatory variant discovery, but its use for this purpose has not been systematically evaluated. Here we apply seven variant callers to bulk and single-cell ATAC-seq data and evaluate their ability to identify single nucleotide variants (SNVs) and insertions/deletions (indels). In addition, we develop an ensemble classifier, VarCA, which combines features from individual variant callers to predict variants. The Genome Analysis Toolkit (GATK) is the best-performing individual caller with precision/recall on a bulk ATAC test dataset of 0.92/0.97 for SNVs and 0.87/0.82 for indels within ATAC-seq peak regions with at least 10 reads. On bulk ATAC-seq reads, VarCA achieves superior performance with precision/recall of 0.99/0.95 for SNVs and 0.93/0.80 for indels. On single-cell ATAC-seq reads, VarCA attains precision/recall of 0.98/0.94 for SNVs and 0.82/0.82 for indels. In summary, ATAC-seq reads can be used to accurately discover non-coding regulatory variants in the absence of whole-genome sequencing data and our ensemble method, VarCA, has the best overall performance.

Complete loss of H3K9 methylation dissolves mouse heterochromatin organization.

Histone H3 lysine 9 (H3K9) methylation is a central epigenetic modification that defines heterochromatin from unicellular to multicellular organisms. In mammalian cells, H3K9 methylation can be catalyzed by at least six distinct SET domain enzymes: Suv39h1/Suv39h2, Eset1/Eset2 and G9a/Glp. We used mouse embryonic fibroblasts (MEFs) with a conditional mutation for Eset1 and introduced progressive deletions for the other SET domain genes by CRISPR/Cas9 technology. Compound mutant MEFs for all six SET domain lysine methyltransferase (KMT) genes lack all H3K9 methylation states, derepress nearly all families of repeat elements and display genomic instabilities. Strikingly, the 6KO H3K9 KMT MEF cells no longer maintain heterochromatin organization and have lost electron-dense heterochromatin. This is a compelling analysis of H3K9 methylation-deficient mammalian chromatin and reveals a definitive function for H3K9 methylation in protecting heterochromatin organization and genome integrity.

m6A RNA methylation of major satellite repeat transcripts facilitates chromatin association and RNA:DNA hybrid formation in mouse heterochromatin.

Heterochromatin has essential functions in maintaining chromosome structure, in protecting genome integrity and in stabilizing gene expression programs. Heterochromatin is often nucleated by underlying DNA repeat sequences, such as major satellite repeats (MSR) and long interspersed nuclear elements (LINE). In order to establish heterochromatin, MSR and LINE elements need to be transcriptionally competent and generate non-coding repeat RNA that remain chromatin associated. We explored whether these heterochromatic RNA, similar to DNA and histones, may be methylated, particularly for 5-methylcytosine (5mC) or methyl-6-adenosine (m6A). Our analysis in mouse ES cells identifies only background level of 5mC but significant enrichment for m6A on heterochromatic RNA. Moreover, MSR transcripts are a novel target for m6A RNA modification, and their m6A RNA enrichment is decreased in ES cells that are mutant for Mettl3 or Mettl14, which encode components of a central RNA methyltransferase complex. Importantly, MSR transcripts that are partially deficient in m6A RNA methylation display impaired chromatin association and have a reduced potential to form RNA:DNA hybrids. We propose that m6A modification of MSR RNA will enhance the functions of MSR repeat transcripts to stabilize mouse heterochromatin.

Altered Neuronal Support and Inflammatory Response in Bipolar Disorder Patient-Derived Astrocytes.

Bipolar disorder (BD) is characterized by cyclical mood shifts. Studies indicate that BD patients have a peripheral pro-inflammatory state and alterations in glial populations in the brain. We utilized an in vitro model to study inflammation-related phenotypes of astrocytes derived from induced pluripotent stem cells (iPSCs) generated from BD patients and healthy controls. BD astrocytes showed changes in transcriptome and induced a reduction in neuronal activity when co-cultured with neurons. IL-1ß-stimulated BD astrocytes displayed a unique inflammatory gene expression signature and increased secretion of IL-6. Conditioned medium from stimulated BD astrocytes reduced neuronal activity, and this effect was partially blocked by IL-6 inactivating antibody. Our results suggest that BD astrocytes are functionally less supportive of neuronal excitability and this effect is partially mediated by IL-6. We confirmed higher IL-6 in blood in a distinct cohort of BD patients, highlighting the potential role of astrocyte-mediated inflammatory signaling in BD neuropathology.

Deficient LEF1 expression is associated with lithium resistance and hyperexcitability in neurons derived from bipolar disorder patients.

Bipolar disorder (BD) is a psychiatric condition characterized by depressive and manic episodes that affect 2% of the world population. The first-line long-term treatment for mood stabilization is lithium (Li). Induced pluripotent stem cell modeling of BD using hippocampal dentate gyrus-like neurons derived from Li-responsive (LR) and Li-non-responsive (NR) patients previously showed neuronal hyperexcitability. Li treatment reversed hyperexcitability only on the LR neurons. In this study we searched for specific targets of Li resistance in NR neurons and found that the activity of Wnt/ß-catenin signaling pathway was severely affected, with a significant decrease in expression of LEF1. Li targets the Wnt/ß-catenin signaling pathway by inhibiting GSK-3ß and releasing ß-catenin that forms a nuclear complex with TCF/LEF1, activating the Wnt/ß-catenin transcription program. Therefore, we propose that downregulation of LEF1 may account for Li resistance in NR neurons. Our results show that valproic acid (VPA), a drug used to treat NR patients that also acts downstream of GSK-3ß, upregulated LEF1 and Wnt/ß-catenin gene targets, increased transcriptional activity of complex ß-catenin/TCF/LEF1, and reduced excitability in NR neurons. In addition, decreasing LEF1 expression in control neurons using shLEF1 caused hyperexcitability, confirming that the impact of VPA on excitability in NR neurons was connected to changes in LEF1 and in the Wnt/ß-catenin pathway. Our results suggest that LEF1 may be a useful target for the discovery of new drugs for BD treatment.

Tuft Cells Inhibit Pancreatic Tumorigenesis in Mice by Producing Prostaglandin D2.

BACKGROUND & AIMS: Development of pancreatic ductal adenocarcinoma (PDA) involves acinar to ductal metaplasia and genesis of tuft cells. It has been a challenge to study these rare cells because of the lack of animal models. We investigated the role of tuft cells in pancreatic tumorigenesis. METHODS: We performed studies with LSL-KrasG12D/+;Ptf1aCre/+ mice (KC; develop pancreatic tumors), KC mice crossed with mice with pancreatic disruption of Pou2f3 (KPouC mice; do not develop tuft cells), or mice with pancreatic disruption of the hematopoietic prostaglandin D synthase gene (Hpgds, KHC mice) and wild-type mice. Mice were allowed to age or were given caerulein to induce pancreatitis; pancreata were collected and analyzed by histology, immunohistochemistry, RNA sequencing, ultrastructural microscopy, and metabolic profiling. We performed laser-capture dissection and RNA-sequencing analysis of pancreatic tissues from 26 patients with pancreatic intraepithelial neoplasia (PanIN), 19 patients with intraductal papillary mucinous neoplasms (IPMNs), and 197 patients with PDA. RESULTS: Pancreata from KC mice had increased formation of tuft cells and higher levels of prostaglandin D2 than wild-type mice. Pancreas-specific deletion of POU2F3 in KC mice (KPouC mice) resulted in a loss of tuft cells and accelerated tumorigenesis. KPouC mice had increased fibrosis and activation of immune cells after administration of caerulein. Pancreata from KPouC and KHC mice had significantly lower levels of prostaglandin D2, compared with KC mice, and significantly increased numbers of PanINs and PDAs. KPouC and KHC mice had increased pancreatic injury after administration of caerulein, significantly less normal tissue, more extracellular matrix deposition, and higher PanIN grade than KC mice. Human PanIN and intraductal papillary mucinous neoplasm had gene expression signatures associated with tuft cells and increased expression of Hpgds messenger RNA compared with PDA. CONCLUSIONS: In mice with KRAS-induced pancreatic tumorigenesis, loss of tuft cells accelerates tumorigenesis and increases the severity of caerulein-induced pancreatic injury, via decreased production of prostaglandin D2. These data are consistent with the hypothesis that tuft cells are a metaplasia-induced tumor attenuating cell type.

Dynamical Electrical Complexity Is Reduced during Neuronal Differentiation in Autism Spectrum Disorder.

Neuronal activity can be modeled as a nonlinear dynamical system to yield measures of neuronal state and dysfunction. The electrical recordings of stem cell-derived neurons from individuals with autism spectrum disorder (ASD) and controls were analyzed using minimum embedding dimension (MED) analysis to characterize their dynamical complexity. MED analysis revealed a significant reduction in dynamical complexity in ASD neurons during differentiation, which was correlated to bursting and spike interval measures. MED was associated with clinical endpoints, such as nonverbal intelligence, and was correlated with 53 differentially expressed genes, which were overrepresented with ASD risk genes related to neurodevelopment, cell morphology, and cell migration. Spatiotemporal analysis also showed a prenatal temporal enrichment in cortical and deep brain structures. Together, we present dynamical analysis as a paradigm that can be used to distinguish disease-associated cellular electrophysiological and transcriptional signatures, while taking into account patient variability in neuropsychiatric disorders.

Targeting LIF-mediated paracrine interaction for pancreatic cancer therapy and monitoring.

Pancreatic ductal adenocarcinoma (PDAC) has a dismal prognosis largely owing to inefficient diagnosis and tenacious drug resistance. Activation of pancreatic stellate cells (PSCs) and consequent development of dense stroma are prominent features accounting for this aggressive biology1,2. The reciprocal interplay between PSCs and pancreatic cancer cells (PCCs) not only enhances tumour progression and metastasis but also sustains their own activation, facilitating a vicious cycle to exacerbate tumorigenesis and drug resistance3-7. Furthermore, PSC activation occurs very early during PDAC tumorigenesis8-10, and activated PSCs comprise a substantial fraction of the tumour mass, providing a rich source of readily detectable factors. Therefore, we hypothesized that the communication between PSCs and PCCs could be an exploitable target to develop effective strategies for PDAC therapy and diagnosis. Here, starting with a systematic proteomic investigation of secreted disease mediators and underlying molecular mechanisms, we reveal that leukaemia inhibitory factor (LIF) is a key paracrine factor from activated PSCs acting on cancer cells. Both pharmacologic LIF blockade and genetic Lifr deletion markedly slow tumour progression and augment the efficacy of chemotherapy to prolong survival of PDAC mouse models, mainly by modulating cancer cell differentiation and epithelial-mesenchymal transition status. Moreover, in both mouse models and human PDAC, aberrant production of LIF in the pancreas is restricted to pathological conditions and correlates with PDAC pathogenesis, and changes in the levels of circulating LIF correlate well with tumour response to therapy. Collectively, these findings reveal a function of LIF in PDAC tumorigenesis, and suggest its translational potential as an attractive therapeutic target and circulating marker. Our studies underscore how a better understanding of cell-cell communication within the tumour microenvironment can suggest novel strategies for cancer therapy.

Histone serotonylation is a permissive modification that enhances TFIID binding to H3K4me3.

Chemical modifications of histones can mediate diverse DNA-templated processes, including gene transcription1-3. Here we provide evidence for a class of histone post-translational modification, serotonylation of glutamine, which occurs at position 5 (Q5ser) on histone H3 in organisms that produce serotonin (also known as 5-hydroxytryptamine (5-HT)). We demonstrate that tissue transglutaminase 2 can serotonylate histone H3 tri-methylated lysine 4 (H3K4me3)-marked nucleosomes, resulting in the presence of combinatorial H3K4me3Q5ser in vivo. H3K4me3Q5ser displays a ubiquitous pattern of tissue expression in mammals, with enrichment observed in brain and gut, two organ systems responsible for the bulk of 5-HT production. Genome-wide analyses of human serotonergic neurons, developing mouse brain and cultured serotonergic cells indicate that H3K4me3Q5ser nucleosomes are enriched in euchromatin, are sensitive to cellular differentiation and correlate with permissive gene expression, phenomena that are linked to the potentiation of TFIID4-6 interactions with H3K4me3. Cells that ectopically express a H3 mutant that cannot be serotonylated display significantly altered expression of H3K4me3Q5ser-target loci, which leads to deficits in differentiation. Taken together, these data identify a direct role for 5-HT, independent from its contributions to neurotransmission and cellular signalling, in the mediation of permissive gene expression.

Single-dose CRISPR-Cas9 therapy extends lifespan of mice with Hutchinson-Gilford progeria syndrome.

Hutchinson-Gilford progeria syndrome (HGPS) is a rare lethal genetic disorder characterized by symptoms reminiscent of accelerated aging. The major underlying genetic cause is a substitution mutation in the gene coding for lamin A, causing the production of a toxic isoform called progerin. Here we show that reduction of lamin A/progerin by a single-dose systemic administration of adeno-associated virus-delivered CRISPR-Cas9 components suppresses HGPS in a mouse model.

BART: bioinformatics array research tool.

BACKGROUND: Microarray experiments comprise more than half of all series in the Gene Expression Omnibus (GEO). However, downloading and analyzing raw or semi-processed microarray data from GEO is not intuitive and requires manual error-prone analysis and a bioinformatics background. This is due to a lack of standardization in array platform fabrication as well as the lack of a simple interactive tool for clustering, plotting, differential expression testing, and testing for functional enrichment. RESULTS: We introduce the Bioinformatics Array Research Tool (BART), an R Shiny web application that automates the microarray download and analysis process across diverse microarray platforms. It provides an intuitive interface, automatically downloads and parses data from GEO, suggests groupings of samples for differential expression testing, performs batch effect correction, outputs quality control plots, converts probe IDs, generates full lists of differentially expressed genes, and performs functional enrichment analysis. We show that BART enables a more comprehensive analysis of a wider range of microarray datasets on GEO by comparing it to four leading online microarray analysis tools. CONCLUSIONS: BART allows a scientist with no bioinformatics background to extract knowledge from their own microarray data or microarray experiments available from GEO. BART is functional on more microarray experiments and provides more comprehensive analyses than extant microarray analysis tools. BART is hosted on bart.salk.edu , includes a user tutorial, and is available for download from https://bitbucket.org/Luisa_amaral/bart .

The Aging Astrocyte Transcriptome from Multiple Regions of the Mouse Brain.

Aging brains undergo cognitive decline, associated with decreased neuronal synapse number and function and altered metabolism. Astrocytes regulate neuronal synapse formation and function in development and adulthood, but whether these properties change during aging, contributing to neuronal dysfunction, is unknown. We addressed this by generating aged and adult astrocyte transcriptomes from multiple mouse brain regions. These data provide a comprehensive RNA-seq database of adult and aged astrocyte gene expression, available online as a resource. We identify astrocyte genes altered by aging across brain regions and regionally unique aging changes. Aging astrocytes show minimal alteration of homeostatic and neurotransmission-regulating genes. However, aging astrocytes upregulate genes that eliminate synapses and partially resemble reactive astrocytes. We further identified heterogeneous expression of synapse-regulating genes between astrocytes from different cortical regions. We find that alterations to astrocytes in aging create an environment permissive to synapse elimination and neuronal damage, potentially contributing to aging-associated cognitive decline.

The RNA-Binding Protein NONO Coordinates Hepatic Adaptation to Feeding.

The mechanisms by which feeding and fasting drive rhythmic gene expression for physiological adaptation to daily rhythm in nutrient availability are not well understood. Here we show that, upon feeding, the RNA-binding protein NONO accumulates within speckle-like structures in liver cell nuclei. Combining RNA-immunoprecipitation and sequencing (RIP-seq), we find that an increased number of RNAs are bound by NONO after feeding. We further show that NONO binds and regulates the rhythmicity of genes involved in nutrient metabolism post-transcriptionally. Finally, we show that disrupted rhythmicity of NONO target genes has profound metabolic impact. Indeed, NONO-deficient mice exhibit impaired glucose tolerance and lower hepatic glycogen and lipids. Accordingly, these mice shift from glucose storage to fat oxidation, and therefore remain lean throughout adulthood. In conclusion, our study demonstrates that NONO post-transcriptionally coordinates circadian mRNA expression of metabolic genes with the feeding/fasting cycle, thereby playing a critical role in energy homeostasis.

Reprogramming pancreatic stellate cells via p53 activation: A putative target for pancreatic cancer therapy.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by an extremely dense fibrotic stroma, which contributes to tumor growth, metastasis, and drug resistance. During tumorigenesis, quiescent pancreatic stellate cells (PSCs) are activated and become major contributors to fibrosis, by increasing growth factor signaling and extracellular matrix deposition. The p53 tumor suppressor is known to restrict tumor initiation and progression through cell autonomous mechanisms including apoptosis, cell cycle arrest, and senescence. There is growing evidence that stromal p53 also exerts anti-tumor activity by paracrine mechanisms, though a role for stromal p53 in PDAC has not yet been described. Here, we demonstrate that activation of stromal p53 exerts anti-tumor effects in PDAC. We show that primary cancer-associated PSCs (caPSCs) isolated from human PDAC express wild-type p53, which can be activated by the Mdm2 antagonist Nutlin-3a. Our work reveals that p53 acts as a major regulator of PSC activation and as a modulator of PDAC fibrosis. In vitro, p53 activation by Nutlin-3a induces profound transcriptional changes, which reprogram activated PSCs to quiescence. Using immunofluorescence and lipidomics, we have also found that p53 activation induces lipid droplet accumulation in both normal and tumor-associated fibroblasts, revealing a previously undescribed role for p53 in lipid storage. In vivo, treatment of tumor-bearing mice with the clinical form of Nutlin-3a induces stromal p53 activation, reverses caPSCs activation, and decreases fibrosis. All together our work uncovers new functions for stromal p53 in PDAC.

Differentiation of Inflammation-Responsive Astrocytes from Glial Progenitors Generated from Human Induced Pluripotent Stem Cells.

Astrocyte dysfunction and neuroinflammation are detrimental features in multiple pathologies of the CNS. Therefore, the development of methods that produce functional human astrocytes represents an advance in the study of neurological diseases. Here we report an efficient method for inflammation-responsive astrocyte generation from induced pluripotent stem cells (iPSCs) and embryonic stem cells. This protocol uses an intermediate glial progenitor stage and generates functional astrocytes that show levels of glutamate uptake and calcium activation comparable with those observed in human primary astrocytes. Stimulation of stem cell-derived astrocytes with interleukin-1ß or tumor necrosis factor α elicits a strong and rapid pro-inflammatory response. RNA-sequencing transcriptome profiling confirmed that similar gene expression changes occurred in iPSC-derived and primary astrocytes upon stimulation with interleukin-1ß. This protocol represents an important tool for modeling in-a-dish neurological diseases with an inflammatory component, allowing for the investigation of the role of diseased astrocytes in neuronal degeneration.

Influence of donor age on induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) are being pursued as a source of cells for autologous therapies, many of which will be aimed at aged patients. To explore the impact of age on iPSC quality, we produced iPSCs from blood cells of 16 donors aged 21-100. We find that iPSCs from older donors retain an epigenetic signature of age, which can be reduced through passaging. Clonal expansion via reprogramming also enables the discovery of somatic mutations present in individual donor cells, which are missed by bulk sequencing methods. We show that exomic mutations in iPSCs increase linearly with age, and all iPSC lines analyzed carry at least one gene-disrupting mutation, several of which have been associated with cancer or dysfunction. Unexpectedly, elderly donors (>90 yrs) harbor fewer mutations than predicted, likely due to a contracted blood progenitor pool. These studies establish that donor age is associated with an increased risk of abnormalities in iPSCs and will inform clinical development of reprogramming technology.

A Highly Polymorphic Copy Number Variant in the NSF Gene is Associated with Cocaine Dependence.

Cocaine dependence is a complex psychiatric disorder involving both genetic and environmental factors. Several neurotransmitter systems mediate cocaine's effects, dependence and relapse, being the components of the neurotransmitter release machinery good candidates for the disorder. Previously, we identified a risk haplotype for cocaine dependence in the NSF gene, encoding the protein N-Ethylmaleimide-Sensitive Factor essential for synaptic vesicle turnover. Here we examined the possible contribution to cocaine dependence of a large copy number variant (CNV) that encompasses part of the NSF gene. We performed a case-control association study in a discovery sample (359 cases and 356 controls) and identified an association between cocaine dependence and the CNV (P = 0.013), that was confirmed in the replication sample (508 cases and 569 controls, P = 7.1e-03) and in a pooled analysis (P = 1.8e-04), with an over-representation of low number of copies in cases. Subsequently, we studied the functional impact of the CNV on gene expression and found that the levels of two NSF transcripts were significantly increased in peripheral blood mononuclear cells (PBMC) along with the number of copies of the CNV. These results, together with a previous study from our group, support the role of NSF in the susceptibility to cocaine dependence.