Transgenerational inheritance of acquired epigenetic signatures at CpG islands in mice.

Transgenerational epigenetic inheritance in mammals remains a debated subject. Here, we demonstrate that DNA methylation of promoter-associated CpG islands (CGIs) can be transmitted from parents to their offspring in mice. We generated DNA methylation-edited mouse embryonic stem cells (ESCs), in which CGIs of two metabolism-related genes, the Ankyrin repeat domain 26 and the low-density lipoprotein receptor, were specifically methylated and silenced. DNA methylation-edited mice generated by microinjection of the methylated ESCs exhibited abnormal metabolic phenotypes. Acquired methylation of the targeted CGI and the phenotypic traits were maintained and transmitted across multiple generations. The heritable CGI methylation was subjected to reprogramming in parental PGCs and subsequently reestablished in the next generation at post-implantation stages. These observations provide a concrete step toward demonstrating transgenerational epigenetic inheritance in mammals, which may have implications in our understanding of evolutionary biology as well as the etiology, diagnosis, and prevention of non-genetically inherited human diseases.

An Activity-Guided Map of Electrophile-Cysteine Interactions in Primary Human T Cells.

Electrophilic compounds originating from nature or chemical synthesis have profound effects on immune cells. These compounds are thought to act by cysteine modification to alter the functions of immune-relevant proteins; however, our understanding of electrophile-sensitive cysteines in the human immune proteome remains limited. Here, we present a global map of cysteines in primary human T cells that are susceptible to covalent modification by electrophilic small molecules. More than 3,000 covalently liganded cysteines were found on functionally and structurally diverse proteins, including many that play fundamental roles in immunology. We further show that electrophilic compounds can impair T cell activation by distinct mechanisms involving the direct functional perturbation and/or degradation of proteins. Our findings reveal a rich content of ligandable cysteines in human T cells and point to electrophilic small molecules as a fertile source for chemical probes and ultimately therapeutics that modulate immunological processes and their associated disorders.

Glucocorticoid signaling and regulatory T cells cooperate to maintain the hair-follicle stem-cell niche.

Maintenance of tissue homeostasis is dependent on the communication between stem cells and supporting cells in the same niche. Regulatory T cells (Treg cells) are emerging as a critical component of the stem-cell niche for supporting their differentiation. How Treg cells sense dynamic signals in this microenvironment and communicate with stem cells is mostly unknown. In the present study, by using hair follicles (HFs) to study Treg cell-stem cell crosstalk, we show an unrecognized function of the steroid hormone glucocorticoid in instructing skin-resident Treg cells to facilitate HF stem-cell (HFSC) activation and HF regeneration. Ablation of the glucocorticoid receptor (GR) in Treg cells blocks hair regeneration without affecting immune homeostasis. Mechanistically, GR and Foxp3 cooperate in Treg cells to induce transforming growth factor ß3 (TGF-ß3), which activates Smad2/3 in HFSCs and facilitates HFSC proliferation. The present study identifies crosstalk between Treg cells and HFSCs mediated by the GR-TGF-ß3 axis, highlighting a possible means of manipulating Treg cells to support tissue regeneration.

Transcriptome Engineering with RNA-Targeting Type VI-D CRISPR Effectors.

Class 2 CRISPR-Cas systems endow microbes with diverse mechanisms for adaptive immunity. Here, we analyzed prokaryotic genome and metagenome sequences to identify an uncharacterized family of RNA-guided, RNA-targeting CRISPR systems that we classify as type VI-D. Biochemical characterization and protein engineering of seven distinct orthologs generated a ribonuclease effector derived from Ruminococcus flavefaciens XPD3002 (CasRx) with robust activity in human cells. CasRx-mediated knockdown exhibits high efficiency and specificity relative to RNA interference across diverse endogenous transcripts. As one of the most compact single-effector Cas enzymes, CasRx can also be flexibly packaged into adeno-associated virus. We target virally encoded, catalytically inactive CasRx to cis elements of pre-mRNA to manipulate alternative splicing, alleviating dysregulated tau isoform ratios in a neuronal model of frontotemporal dementia. Our results present CasRx as a programmable RNA-binding module for efficient targeting of cellular RNA, enabling a general platform for transcriptome engineering and future therapeutic development.

Cooperative Metabolic Adaptations in the Host Can Favor Asymptomatic Infection and Select for Attenuated Virulence in an Enteric Pathogen.

Pathogen virulence exists on a continuum. The strategies that drive symptomatic or asymptomatic infections remain largely unknown. We took advantage of the concept of lethal dose 50 (LD50) to ask which component of individual non-genetic variation between hosts defines whether they survive or succumb to infection. Using the enteric pathogen Citrobacter, we found no difference in pathogen burdens between healthy and symptomatic populations. Iron metabolism-related genes were induced in asymptomatic hosts compared to symptomatic or naive mice. Dietary iron conferred complete protection without influencing pathogen burdens, even at 1000× the lethal dose of Citrobacter. Dietary iron induced insulin resistance, increasing glucose levels in the intestine that were necessary and sufficient to suppress pathogen virulence. A short course of dietary iron drove the selection of attenuated Citrobacter strains that can transmit and asymptomatically colonize naive hosts, demonstrating that environmental factors and cooperative metabolic strategies can drive conversion of pathogens toward commensalism.

Uptake of oxidized lipids by the scavenger receptor CD36 promotes lipid peroxidation and dysfunction in CD8+ T cells in tumors.

A common metabolic alteration in the tumor microenvironment (TME) is lipid accumulation, a feature associated with immune dysfunction. Here, we examined how CD8+ tumor infiltrating lymphocytes (TILs) respond to lipids within the TME. We found elevated concentrations of several classes of lipids in the TME and accumulation of these in CD8+ TILs. Lipid accumulation was associated with increased expression of CD36, a scavenger receptor for oxidized lipids, on CD8+ TILs, which also correlated with progressive T cell dysfunction. Cd36-/- T cells retained effector functions in the TME, as compared to WT counterparts. Mechanistically, CD36 promoted uptake of oxidized low-density lipoproteins (OxLDL) into T cells, and this induced lipid peroxidation and downstream activation of p38 kinase. Inhibition of p38 restored effector T cell functions in vitro, and resolution of lipid peroxidation by overexpression of glutathione peroxidase 4 restored functionalities in CD8+ TILs in vivo. Thus, an oxidized lipid-CD36 axis promotes intratumoral CD8+ T cell dysfunction and serves as a therapeutic avenue for immunotherapies.

Telomere-to-mitochondria signalling by ZBP1 mediates replicative crisis.

Cancers arise through the accumulation of genetic and epigenetic alterations that enable cells to evade telomere-based proliferative barriers and achieve immortality. One such barrier is replicative crisis-an autophagy-dependent program that eliminates checkpoint-deficient cells with unstable telomeres and other cancer-relevant chromosomal aberrations1,2. However, little is known about the molecular events that regulate the onset of this important tumour-suppressive barrier. Here we identified the innate immune sensor Z-DNA binding protein 1 (ZBP1) as a regulator of the crisis program. A crisis-associated isoform of ZBP1 is induced by the cGAS-STING DNA-sensing pathway, but reaches full activation only when associated with telomeric-repeat-containing RNA (TERRA) transcripts that are synthesized from dysfunctional telomeres. TERRA-bound ZBP1 oligomerizes into filaments on the outer mitochondrial membrane of a subset of mitochondria, where it activates the innate immune adapter protein mitochondrial antiviral-signalling protein (MAVS). We propose that these oligomerization properties of ZBP1 serve as a signal amplification mechanism, where few TERRA-ZBP1 interactions are sufficient to launch a detrimental MAVS-dependent interferon response. Our study reveals a mechanism for telomere-mediated tumour suppression, whereby dysfunctional telomeres activate innate immune responses through mitochondrial TERRA-ZBP1 complexes to eliminate cells destined for neoplastic transformation.

Dynamic regulation of histone modifications and long-range chromosomal interactions during postmitotic transcriptional reactivation.

During mitosis, transcription of genomic DNA is dramatically reduced, before it is reactivated during nuclear reformation in anaphase/telophase. Many aspects of the underlying principles that mediate transcriptional memory and reactivation in the daughter cells remain unclear. Here, we used ChIP-seq on synchronized cells at different stages after mitosis to generate genome-wide maps of histone modifications. Combined with EU-RNA-seq and Hi-C analyses, we found that during prometaphase, promoters, enhancers, and insulators retain H3K4me3 and H3K4me1, while losing H3K27ac. Enhancers globally retaining mitotic H3K4me1 or locally retaining mitotic H3K27ac are associated with cell type-specific genes and their transcription factors for rapid transcriptional activation. As cells exit mitosis, promoters regain H3K27ac, which correlates with transcriptional reactivation. Insulators also gain H3K27ac and CCCTC-binding factor (CTCF) in anaphase/telophase. This increase of H3K27ac in anaphase/telophase is required for posttranscriptional activation and may play a role in the establishment of topologically associating domains (TADs). Together, our results suggest that the genome is reorganized in a sequential order, in which histone methylations occur first in prometaphase, histone acetylation, and CTCF in anaphase/telophase, transcription in cytokinesis, and long-range chromatin interactions in early G1. We thus provide insights into the histone modification landscape that allows faithful reestablishment of the transcriptional program and TADs during cell division.

AMPK regulation of Raptor and TSC2 mediate metformin effects on transcriptional control of anabolism and inflammation.

Despite being the frontline therapy for type 2 diabetes, the mechanisms of action of the biguanide drug metformin are still being discovered. In particular, the detailed molecular interplays between the AMPK and the mTORC1 pathway in the hepatic benefits of metformin are still ill defined. Metformin-dependent activation of AMPK classically inhibits mTORC1 via TSC/RHEB, but several lines of evidence suggest additional mechanisms at play in metformin inhibition of mTORC1. Here we investigated the role of direct AMPK-mediated serine phosphorylation of RAPTOR in a new RaptorAA mouse model, in which AMPK phospho-serine sites Ser722 and Ser792 of RAPTOR were mutated to alanine. Metformin treatment of primary hepatocytes and intact murine liver requires AMPK regulation of both RAPTOR and TSC2 to fully inhibit mTORC1, and this regulation is critical for both the translational and transcriptional response to metformin. Transcriptionally, AMPK and mTORC1 were both important for regulation of anabolic metabolism and inflammatory programs triggered by metformin treatment. The hepatic transcriptional response in mice on high-fat diet treated with metformin was largely ablated by AMPK deficiency under the conditions examined, indicating the essential role of this kinase and its targets in metformin action in vivo.

Paternal RLIM/Rnf12 is a survival factor for milk-producing alveolar cells.

In female mouse embryos, somatic cells undergo a random form of X chromosome inactivation (XCI), whereas extraembryonic trophoblast cells in the placenta undergo imprinted XCI, silencing exclusively the paternal X chromosome. Initiation of imprinted XCI requires a functional maternal allele of the X-linked gene Rnf12, which encodes the ubiquitin ligase Rnf12/RLIM. We find that knockout (KO) of Rnf12 in female mammary glands inhibits alveolar differentiation and milk production upon pregnancy, with alveolar cells that lack RLIM undergoing apoptosis as they begin to differentiate. Genetic analyses demonstrate that these functions are mediated primarily by the paternal Rnf12 allele due to nonrandom maternal XCI in mammary epithelial cells. These results identify paternal Rnf12/RLIM as a critical survival factor for milk-producing alveolar cells and, together with population models, reveal implications of transgenerational epigenetic inheritance.

Heterozygous Mutations in SMARCA2 Reprogram the Enhancer Landscape by Global Retargeting of SMARCA4.

Mammalian SWI/SNF complexes are multi-subunit chromatin remodeling complexes associated with an ATPase (either SMARCA4 or SMARCA2). Heterozygous mutations in the SMARCA2 ATPase cause Nicolaides-Baraitser syndrome (NCBRS), an intellectual disability syndrome associated with delayed speech onset. We engineered human embryonic stem cells (hESCs) to carry NCBRS-associated heterozygous SMARCA2 K755R or R1159Q mutations. While SMARCA2 mutant hESCs were phenotypically normal, differentiation to neural progenitors cells (NPCs) was severely impaired. We find that SMARCA2 mutations cause enhancer reorganization with loss of SOX3-dependent neural enhancers and prominent emergence of astrocyte-specific de novo enhancers. Changes in chromatin accessibility at enhancers were associated with an increase in SMARCA2 binding and retargeting of SMARCA4. We show that the AP-1 family member FRA2 is aberrantly overexpressed in SMARCA2 mutant NPCs, where it functions as a pioneer factor at de novo enhancers. Together, our results demonstrate that SMARCA2 mutations cause impaired differentiation through enhancer reprogramming via inappropriate targeting of SMARCA4.

CDK12 phosphorylates 4E-BP1 to enable mTORC1-dependent translation and mitotic genome stability.

The RNA polymerase II (RNAPII) C-terminal domain kinase, CDK12, regulates genome stability, expression of DNA repair genes, and cancer cell resistance to chemotherapy and immunotherapy. In addition to its role in mRNA biosynthesis of DNA repair genes, we show here that CDK12 phosphorylates the mRNA 5' cap-binding repressor, 4E-BP1, to promote translation of mTORC1-dependent mRNAs. In particular, we found that phosphorylation of 4E-BP1 by mTORC1 (T37 and T46) facilitates subsequent CDK12 phosphorylation at two Ser-Pro sites (S65 and T70) that control the exchange of 4E-BP1 with eIF4G at the 5' cap of CHK1 and other target mRNAs. RNA immunoprecipitation coupled with deep sequencing (RIP-seq) revealed that CDK12 regulates release of 4E-BP1, and binding of eIF4G, to many mTORC1 target mRNAs, including those needed for MYC transformation. Genome-wide ribosome profiling (Ribo-seq) further identified specific CDK12 "translation-only" target mRNAs, including many mTORC1 target mRNAs as well as many subunits of mitotic and centromere/centrosome complexes. Accordingly, confocal imaging analyses revealed severe chromosome misalignment, bridging, and segregation defects in cells deprived of CDK12 or CCNK. We conclude that the nuclear RNAPII-CTD kinase CDK12 cooperates with mTORC1, and controls a specialized translation network that is essential for mitotic chromosome stability.

In vivo reprogramming of wound-resident cells generates skin epithelial tissue.

Large cutaneous ulcers are, in severe cases, life threatening1,2. As the global population ages, non-healing ulcers are becoming increasingly common1,2. Treatment currently requires the transplantation of pre-existing epithelial components, such as skin grafts, or therapy using cultured cells2. Here we develop alternative supplies of epidermal coverage for the treatment of these kinds of wounds. We generated expandable epithelial tissues using in vivo reprogramming of wound-resident mesenchymal cells. Transduction of four transcription factors that specify the skin-cell lineage enabled efficient and rapid de novo epithelialization from the surface of cutaneous ulcers in mice. Our findings may provide a new therapeutic avenue for treating skin wounds and could be extended to other disease situations in which tissue homeostasis and repair are impaired.

Dynamic regulation of CTCF stability and sub-nuclear localization in response to stress.

The nuclear protein CCCTC-binding factor (CTCF) has diverse roles in chromatin architecture and gene regulation. Functionally, CTCF associates with thousands of genomic sites and interacts with proteins, such as cohesin, or non-coding RNAs to facilitate specific transcriptional programming. In this study, we examined CTCF during the cellular stress response in human primary cells using immune-blotting, quantitative real time-PCR, chromatin immunoprecipitation-sequence (ChIP-seq) analysis, mass spectrometry, RNA immunoprecipitation-sequence analysis (RIP-seq), and Airyscan confocal microscopy. Unexpectedly, we found that CTCF is exquisitely sensitive to diverse forms of stress in normal patient-derived human mammary epithelial cells (HMECs). In HMECs, a subset of CTCF protein forms complexes that localize to Serine/arginine-rich splicing factor (SC-35)-containing nuclear speckles. Upon stress, this species of CTCF protein is rapidly downregulated by changes in protein stability, resulting in loss of CTCF from SC-35 nuclear speckles and changes in CTCF-RNA interactions. Our ChIP-seq analysis indicated that CTCF binding to genomic DNA is largely unchanged. Restoration of the stress-sensitive pool of CTCF protein abundance and re-localization to nuclear speckles can be achieved by inhibition of proteasome-mediated degradation. Surprisingly, we observed the same characteristics of the stress response during neuronal differentiation of human pluripotent stem cells (hPSCs). CTCF forms stress-sensitive complexes that localize to SC-35 nuclear speckles during a specific stage of neuronal commitment/development but not in differentiated neurons. We speculate that these particular CTCF complexes serve a role in RNA processing that may be intimately linked with specific genes in the vicinity of nuclear speckles, potentially to maintain cells in a certain differentiation state, that is dynamically regulated by environmental signals. The stress-regulated activity of CTCF is uncoupled in persistently stressed, epigenetically re-programmed "variant" HMECs and certain cancer cell lines. These results reveal new insights into CTCF function in cell differentiation and the stress-response with implications for oxidative damage-induced cancer initiation and neuro-degenerative diseases.

AMPK governs lineage specification through Tfeb-dependent regulation of lysosomes.

Faithful execution of developmental programs relies on the acquisition of unique cell identities from pluripotent progenitors, a process governed by combinatorial inputs from numerous signaling cascades that ultimately dictate lineage-specific transcriptional outputs. Despite growing evidence that metabolism is integrated with many molecular networks, how pathways that control energy homeostasis may affect cell fate decisions is largely unknown. Here, we show that AMP-activated protein kinase (AMPK), a central metabolic regulator, plays critical roles in lineage specification. Although AMPK-deficient embryonic stem cells (ESCs) were normal in the pluripotent state, these cells displayed profound defects upon differentiation, failing to generate chimeric embryos and preferentially adopting an ectodermal fate at the expense of the endoderm during embryoid body (EB) formation. AMPK(-/-) EBs exhibited reduced levels of Tfeb, a master transcriptional regulator of lysosomes, leading to diminished endolysosomal function. Remarkably, genetic loss of Tfeb also yielded endodermal defects, while AMPK-null ESCs overexpressing this transcription factor normalized their differential potential, revealing an intimate connection between Tfeb/lysosomes and germ layer specification. The compromised endolysosomal system resulting from AMPK or Tfeb inactivation blunted Wnt signaling, while up-regulating this pathway restored expression of endodermal markers. Collectively, these results uncover the AMPK pathway as a novel regulator of cell fate determination during differentiation.

An RNA polymerase ribozyme that synthesizes its own ancestor.

The RNA-based organisms from which modern life is thought to have descended would have depended on an RNA polymerase ribozyme to copy functional RNA molecules, including copying the polymerase itself. Such a polymerase must have been capable of copying structured RNAs with high efficiency and high fidelity to maintain genetic information across successive generations. Here the class I RNA polymerase ribozyme was evolved in vitro for the ability to synthesize functional ribozymes, resulting in the markedly improved ability to synthesize complex RNAs using nucleoside 5'-triphosphate (NTP) substrates. The polymerase is descended from the class I ligase, which contains the same catalytic core as the polymerase. The class I ligase can be synthesized by the improved polymerase as three separate RNA strands that assemble to form a functional ligase. The polymerase also can synthesize the complement of each of these three strands. Despite this remarkable level of activity, only a very small fraction of the assembled ligases retain catalytic activity due to the presence of disabling mutations. Thus, the fidelity of RNA polymerization should be considered a major impediment to the construction of a self-sustained, RNA-based evolving system. The propagation of heritable information requires both efficient and accurate synthesis of genetic molecules, a requirement relevant to both laboratory systems and the early history of life on Earth.

Accurate annotation of human protein-coding small open reading frames.

Functional protein-coding small open reading frames (smORFs) are emerging as an important class of genes. However, the number of translated smORFs in the human genome is unclear because proteogenomic methods are not sensitive enough, and, as we show, Ribo-seq strategies require additional measures to ensure comprehensive and accurate smORF annotation. Here, we integrate de novo transcriptome assembly and Ribo-seq into an improved workflow that overcomes obstacles with previous methods, to more confidently annotate thousands of smORFs. Evolutionary conservation analyses suggest that hundreds of smORF-encoded microproteins are likely functional. Additionally, many smORFs are regulated during fundamental biological processes, such as cell stress. Peptides derived from smORFs are also detectable on human leukocyte antigen complexes, revealing smORFs as a source of antigens. Thus, by including additional validation into our smORF annotation workflow, we accurately identify thousands of unannotated translated smORFs that will provide a rich pool of unexplored, functional human genes.

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.

SUMMER, a shiny utility for metabolomics and multiomics exploratory research.

INTRODUCTION: Cellular metabolites are generated by a complex network of biochemical reactions. This makes interpreting changes in metabolites exceptionally challenging. OBJECTIVES: To develop a computational tool that integrates multiomics data at the level of reactions. METHODS: Changes in metabolic reactions are modeled with input from transcriptomics/proteomics measurements of enzymes and metabolomic measurements of metabolites. RESULTS: We developed SUMMER, which identified more relevant signals, key metabolic reactions, and relevant underlying biological pathways in a real-world case study. CONCLUSION: SUMMER performs integrative analysis for data interpretation and exploration. SUMMER is freely accessible at http://summer.salk.edu and the code is available at https://bitbucket.org/salkigc/summer .