Search

Although metastatic disease is the leading cause of cancer-related deaths, its tumor microenvironment remains poorly characterized due to technical and biospecimen limitations. In this study, we assembled a multi-modal spatial and cellular map of 67 tumor biopsies from 60 patients with metastatic breast cancer across diverse clinicopathological features and nine anatomic sites with detailed clinical annotations. We combined single-cell or single-nucleus RNA sequencing for all biopsies with a panel of four spatial expression assays (Slide-seq, MERFISH, ExSeq and CODEX) and H&E staining of consecutive serial sections from up to 15 of these biopsies. We leveraged the coupled measurements to provide reference points for the utility and integration of different experimental techniques and used them to assess variability in cell type composition and expression as well as emerging spatial expression characteristics across clinicopathological and methodological diversity. Finally, we assessed spatial expression and co-localization features of macrophage populations, characterized three distinct spatial phenotypes of epithelial-to-mesenchymal transition and identified expression programs associated with local T cell infiltration versus exclusion, showcasing the potential of clinically relevant discovery in such maps., Competing Interests: Competing interests A. Regev, J.K. and D.L.A. are co-inventors on patent application number 17/156,392, filed by the Broad Institute, for inventions relating to work in this manuscript. A. Regev is a co-founder and equity holder of Celsius Therapeutics, an equity holder in Immunitas and, until 31 July 2020, a scientific advisory board (SAB) member of Thermo Fisher Scientific, Syros Pharmaceuticals, Neogene Therapeutics and Asimov. From 1 August 2020, A. Regev is an employee of Genentech and has equity in Roche. A. Rotem is an employee of AstraZeneca since September 2020 and has equity in AstraZeneca. E.S.B. co-founded a company that is exploring commercial applications of expansion microscopy. G.N. holds equity in and consults for Akoya Biosciences. J.G.T.Z. owns stocks in the biotechnology exchange-traded funds CNCR, IDNA, IBB and XBI, owns stock in Novo Nordisk and owns stock in Adaptive Biotechnologies, 2seventy bio and bluebird bio. J.J-V. is a contractor at Genentech since 2021. M.S. is a contractor at Genentech since November 2020. O.R.-R. is a co-inventor on patent applications filed by the Broad Institute for inventions related to single-cell genomics. She has given numerous lectures on the subject of single-cell genomics to a wide variety of audiences and, in some cases, has received remuneration to cover time and costs. O.R.-R. is an employee of Genentech since 19 October 2020 and has equity in Roche. S.J.R. is a member of the SAB of Immunitas Therapeutics and receives research funding from Bristol Myers Squibb and Kite/Gilead. X.Z. is a co-founder of and consultant for Vizgen. The other authors declare no competing interests., (© 2024. The Author(s).)

Spatio-molecular data and microscopy images provide complementary information, essential to study structure and function of spatially organised multicellular systems such as healthy or diseased tissues. However, aligning these two types of data can be challenging due to distortions and differences in resolution, orientation, and position. Manual registration is tedious but may be necessary for challenging samples as well as for the generation of ground-truth data sets that enable benchmarking of existing and emerging automated alignment tools. To make the process of manual registration more convenient, efficient, and integrated, we created BoReMi, a python-based, Jupyter-integrated, visual tool that offers all the relevant functionalities for aligning and registering spatio-molecular data and associated microscopy images. We showcase BoReMi's utility using publicly available data and images and make BoReMi as well as an interactive demo available on GitHub., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Ishar et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

UNC93B1 is critical for trafficking and function of nucleic acid-sensing Toll-like receptors (TLRs) TLR3, TLR7, TLR8, and TLR9, which are essential for antiviral immunity. Overactive TLR7 signaling induced by recognition of self-nucleic acids has been implicated in systemic lupus erythematosus (SLE). Here, we report UNC93B1 variants (E92G and R336L) in four patients with early-onset SLE. Patient cells or mouse macrophages carrying the UNC93B1 variants produced high amounts of TNF-α and IL-6 and upon stimulation with TLR7/TLR8 agonist, but not with TLR3 or TLR9 agonists. E92G causes UNC93B1 protein instability and reduced interaction with TLR7, leading to selective TLR7 hyperactivation with constitutive type I IFN signaling. Thus, UNC93B1 regulates TLR subtype-specific mechanisms of ligand recognition. Our findings establish a pivotal role for UNC93B1 in TLR7-dependent autoimmunity and highlight the therapeutic potential of targeting TLR7 in SLE.

Background: Deep learning has revolutionized medical image analysis in cancer pathology, where it had a substantial clinical impact by supporting the diagnosis and prognostic rating of cancer. Among the first available digital resources in the field of brain cancer is glioblastoma, the most common and fatal brain cancer. At the histologic level, glioblastoma is characterized by abundant phenotypic variability that is poorly linked with patient prognosis. At the transcriptional level, 3 molecular subtypes are distinguished with mesenchymal-subtype tumors being associated with increased immune cell infiltration and worse outcome., Results: We address genotype-phenotype correlations by applying an Xception convolutional neural network to a discovery set of 276 digital hematozylin and eosin (H&E) slides with molecular subtype annotation and an independent The Cancer Genome Atlas-based validation cohort of 178 cases. Using this approach, we achieve high accuracy in H&E-based mapping of molecular subtypes (area under the curve for classical, mesenchymal, and proneural = 0.84, 0.81, and 0.71, respectively; P < 0.001) and regions associated with worse outcome (univariable survival model P < 0.001, multivariable P = 0.01). The latter were characterized by higher tumor cell density (P < 0.001), phenotypic variability of tumor cells (P < 0.001), and decreased T-cell infiltration (P = 0.017)., Conclusions: We modify a well-known convolutional neural network architecture for glioblastoma digital slides to accurately map the spatial distribution of transcriptional subtypes and regions predictive of worse outcome, thereby showcasing the relevance of artificial intelligence-enabled image mining in brain cancer., (© The Author(s) 2024. Published by Oxford University Press GigaScience.)

Background: The prognostic roles of clinical and laboratory markers have been exploited to model risk in patients with primary CNS lymphoma, but these approaches do not fully explain the observed variation in outcome. To date, neuroimaging or molecular information is not used. The aim of this study was to determine the utility of radiomic features to capture clinically relevant phenotypes, and to link those to molecular profiles for enhanced risk stratification., Methods: In this retrospective study, we investigated 133 patients across 9 sites in Austria (2005-2018) and an external validation site in South Korea (44 patients, 2013-2016). We used T1-weighted contrast-enhanced MRI and an L1-norm regularized Cox proportional hazard model to derive a radiomic risk score. We integrated radiomic features with DNA methylation profiles using machine learning-based prediction, and validated the most relevant biological associations in tissues and cell lines., Results: The radiomic risk score, consisting of 20 mostly textural features, was a strong and independent predictor of survival (multivariate hazard ratio = 6.56 [3.64-11.81]) that remained valid in the external validation cohort. Radiomic features captured gene regulatory differences such as in BCL6 binding activity, which was put forth as testable treatment target for a subset of patients., Conclusions: The radiomic risk score was a robust and complementary predictor of survival and reflected characteristics in underlying DNA methylation patterns. Leveraging imaging phenotypes to assess risk and inform epigenetic treatment targets provides a concept on which to advance prognostic modeling and precision therapy for this aggressive cancer., Competing Interests: M.P. has received honoraria for lectures, consultation, or advisory board participation from the following for-profit companies: Bayer, Bristol-Myers Squibb, Novartis, Gerson Lehrman Group (GLG), CMC Contrast, GlaxoSmithKline, Mundipharma, Roche, BMJ Journals, MedMedia, Astra Zeneca, AbbVie, Lilly, Medahead, Daiichi Sankyo, Sanofi, Merck Sharp & Dome, Tocagen, Adastra, Gan & Lee Pharmaceuticals, but declares no nonfinancial competing interests. G.L. holds shares in the company contextflow and has received honoraria for lectures from the following for-profit companies: Boehringer Ingelheim, Novartis, and declares no nonfinancial competing interests. All other authors declare no financial or nonfinancial conflicts of interest., (© The Author(s) 2023. Published by Oxford University Press, the Society for Neuro-Oncology and the European Association of Neuro-Oncology.)

Transferring annotations of single-cell-, spatial- and multi-omics data is often challenging owing both to technical limitations, such as low spatial resolution or high dropout fraction, and to biological variations, such as continuous spectra of cell states. Based on the concept that these data are often best described as continuous mixtures of cells or molecules, we present a computational framework for the transfer of annotations to cells and their combinations (TACCO), which consists of an optimal transport model extended with different wrappers to annotate a wide variety of data. We apply TACCO to identify cell types and states, decipher spatiomolecular tissue structure at the cell and molecular level and resolve differentiation trajectories using synthetic and biological datasets. While matching or exceeding the accuracy of specialized tools for the individual tasks, TACCO reduces the computational requirements by up to an order of magnitude and scales to larger datasets (for example, considering the runtime of annotation transfer for 1 M simulated dropout observations)., (© 2023. The Author(s).)

Competing Interests: Declaration of interests M.L. receives research funding from Biogen Inc. unrelated to the current manuscript. E.J.F. is on the scientific advisory board of Viosera Therapeutics and is a paid consultant for Merck and Mestag Therapeutics. J.K. is an inventor on patent application number 17156392 entitled “Molecular spatial mapping of metastatic tumor microenvironment.”

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution., (© 2023. The Author(s).)

Single-cell RNA sequencing has revealed extensive transcriptional cell state diversity in cancer, often observed independently of genetic heterogeneity, raising the central question of how malignant cell states are encoded epigenetically. To address this, here we performed multiomics single-cell profiling-integrating DNA methylation, transcriptome and genotype within the same cells-of diffuse gliomas, tumors characterized by defined transcriptional cell state diversity. Direct comparison of the epigenetic profiles of distinct cell states revealed key switches for state transitions recapitulating neurodevelopmental trajectories and highlighted dysregulated epigenetic mechanisms underlying gliomagenesis. We further developed a quantitative framework to directly measure cell state heritability and transition dynamics based on high-resolution lineage trees in human samples. We demonstrated heritability of malignant cell states, with key differences in hierarchal and plastic cell state architectures in IDH-mutant glioma versus IDH-wild-type glioblastoma, respectively. This work provides a framework anchoring transcriptional cancer cell states in their epigenetic encoding, inheritance and transition dynamics., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Methods for highly multiplexed RNA imaging are limited in spatial resolution and thus in their ability to localize transcripts to nanoscale and subcellular compartments. We adapt expansion microscopy, which physically expands biological specimens, for long-read untargeted and targeted in situ RNA sequencing. We applied untargeted expansion sequencing (ExSeq) to the mouse brain, which yielded the readout of thousands of genes, including splice variants. Targeted ExSeq yielded nanoscale-resolution maps of RNAs throughout dendrites and spines in the neurons of the mouse hippocampus, revealing patterns across multiple cell types, layer-specific cell types across the mouse visual cortex, and the organization and position-dependent states of tumor and immune cells in a human metastatic breast cancer biopsy. Thus, ExSeq enables highly multiplexed mapping of RNAs from nanoscale to system scale., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Single-cell genomics is essential to chart tumor ecosystems. Although single-cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumors, single-nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each requires customization to different tissue and tumor types, posing a barrier to adoption. Here, we have developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We analyzed 216,490 cells and nuclei from 40 samples across 23 specimens spanning eight tumor types of varying tissue and sample characteristics. We evaluated protocols by cell and nucleus quality, recovery rate and cellular composition. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types, but at different proportions. Our work provides guidance for studies in a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.

Spatial and molecular characteristics determine tissue function, yet high-resolution methods to capture both concurrently are lacking. Here, we developed high-definition spatial transcriptomics, which captures RNA from histological tissue sections on a dense, spatially barcoded bead array. Each experiment recovers several hundred thousand transcript-coupled spatial barcodes at 2-μm resolution, as demonstrated in mouse brain and primary breast cancer. This opens the way to high-resolution spatial analysis of cells and tissues.

The marsupial Tasmanian devil (Sarcophilus harrisii) faces extinction due to transmissible devil facial tumor disease (DFTD). To unveil the molecular underpinnings of this transmissible cancer, we combined pharmacological screens with an integrated systems-biology characterization. Sensitivity to inhibitors of ERBB tyrosine kinases correlated with their overexpression. Proteomic and DNA methylation analyses revealed tumor-specific signatures linked to the evolutionary conserved oncogenic STAT3. ERBB inhibition blocked phosphorylation of STAT3 and arrested cancer cells. Pharmacological blockade of ERBB or STAT3 prevented tumor growth in xenograft models and restored MHC class I expression. This link between the hyperactive ERBB-STAT3 axis and major histocompatibility complex class I-mediated tumor immunosurveillance provides mechanistic insights into horizontal transmissibility and puts forward a dual chemo-immunotherapeutic strategy to save Tasmanian devils from DFTD. VIDEO ABSTRACT., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

Glioblastoma is characterized by widespread genetic and transcriptional heterogeneity, yet little is known about the role of the epigenome in glioblastoma disease progression. Here, we present genome-scale maps of DNA methylation in matched primary and recurring glioblastoma tumors, using data from a highly annotated clinical cohort that was selected through a national patient registry. We demonstrate the feasibility of DNA methylation mapping in a large set of routinely collected FFPE samples, and we validate bisulfite sequencing as a multipurpose assay that allowed us to infer a range of different genetic, epigenetic, and transcriptional characteristics of the profiled tumor samples. On the basis of these data, we identified subtle differences between primary and recurring tumors, links between DNA methylation and the tumor microenvironment, and an association of epigenetic tumor heterogeneity with patient survival. In summary, this study establishes an open resource for dissecting DNA methylation heterogeneity in a genetically diverse and heterogeneous cancer, and it demonstrates the feasibility of integrating epigenomics, radiology, and digital pathology for a national cohort, thereby leveraging existing samples and data collected as part of routine clinical practice.

Micronutrient status of parents can affect long term health of their progeny. Around 2 billion humans are affected by chronic micronutrient deficiency. In this study we use zebrafish as a model system to examine morphological, molecular and epigenetic changes in mature offspring of parents that experienced a one-carbon (1-C) micronutrient deficiency. Zebrafish were fed a diet sufficient, or marginally deficient in 1-C nutrients (folate, vitamin B12, vitamin B6, methionine, choline), and then mated. Offspring livers underwent histological examination, RNA sequencing and genome-wide DNA methylation analysis. Parental 1-C micronutrient deficiency resulted in increased lipid inclusion and we identified 686 differentially expressed genes in offspring liver, the majority of which were downregulated. Downregulated genes were enriched for functional categories related to sterol, steroid and lipid biosynthesis, as well as mitochondrial protein synthesis. Differential DNA methylation was found at 2869 CpG sites, enriched in promoter regions and permutation analyses confirmed the association with parental feed. Our data indicate that parental 1-C nutrient status can persist as locus specific DNA methylation marks in descendants and suggest an effect on lipid utilization and mitochondrial protein translation in F 1 livers. This points toward parental micronutrients status as an important factor for offspring health and welfare.

CRISPR-based genetic screens are accelerating biological discovery, but current methods have inherent limitations. Widely used pooled screens are restricted to simple readouts including cell proliferation and sortable marker proteins. Arrayed screens allow for comprehensive molecular readouts such as transcriptome profiling, but at much lower throughput. Here we combine pooled CRISPR screening with single-cell RNA sequencing into a broadly applicable workflow, directly linking guide RNA expression to transcriptome responses in thousands of individual cells. Our method for CRISPR droplet sequencing (CROP-seq) enables pooled CRISPR screens with single-cell transcriptome resolution, which will facilitate high-throughput functional dissection of complex regulatory mechanisms and heterogeneous cell populations.

Developmental tumors in children and young adults carry few genetic alterations, yet they have diverse clinical presentation. Focusing on Ewing sarcoma, we sought to establish the prevalence and characteristics of epigenetic heterogeneity in genetically homogeneous cancers. We performed genome-scale DNA methylation sequencing for a large cohort of Ewing sarcoma tumors and analyzed epigenetic heterogeneity on three levels: between cancers, between tumors, and within tumors. We observed consistent DNA hypomethylation at enhancers regulated by the disease-defining EWS-FLI1 fusion protein, thus establishing epigenomic enhancer reprogramming as a ubiquitous and characteristic feature of Ewing sarcoma. DNA methylation differences between tumors identified a continuous disease spectrum underlying Ewing sarcoma, which reflected the strength of an EWS-FLI1 regulatory signature and a continuum between mesenchymal and stem cell signatures. There was substantial epigenetic heterogeneity within tumors, particularly in patients with metastatic disease. In summary, our study provides a comprehensive assessment of epigenetic heterogeneity in Ewing sarcoma and thereby highlights the importance of considering nongenetic aspects of tumor heterogeneity in the context of cancer biology and personalized medicine.

Type 1 diabetes is characterized by the destruction of pancreatic β cells, and generating new insulin-producing cells from other cell types is a major aim of regenerative medicine. One promising approach is transdifferentiation of developmentally related pancreatic cell types, including glucagon-producing α cells. In a genetic model, loss of the master regulatory transcription factor Arx is sufficient to induce the conversion of α cells to functional β-like cells. Here, we identify artemisinins as small molecules that functionally repress Arx by causing its translocation to the cytoplasm. We show that the protein gephyrin is the mammalian target of these antimalarial drugs and that the mechanism of action of these molecules depends on the enhancement of GABA A receptor signaling. Our results in zebrafish, rodents, and primary human pancreatic islets identify gephyrin as a druggable target for the regeneration of pancreatic β cell mass from α cells., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Some members of the physiological human microbiome occasionally cause life-threatening disease even in immunocompetent individuals. A prime example of such a commensal pathogen is Neisseria meningitidis, which normally resides in the human nasopharynx but is also a leading cause of sepsis and epidemic meningitis. Using N. meningitidis as model organism, we tested the hypothesis that virulence of commensal pathogens is a consequence of within host evolution and selection of invasive variants due to mutations at contingency genes, a mechanism called phase variation. In line with the hypothesis that phase variation evolved as an adaptation to colonize diverse hosts, computational comparisons of all 27 to date completely sequenced and annotated meningococcal genomes retrieved from public databases showed that contingency genes are indeed enriched for genes involved in host interactions. To assess within-host genetic changes in meningococci, we further used ultra-deep whole-genome sequencing of throat-blood strain pairs isolated from four patients suffering from invasive meningococcal disease. We detected up to three mutations per strain pair, affecting predominantly contingency genes involved in type IV pilus biogenesis. However, there was not a single (set) of mutation(s) that could invariably be found in all four pairs of strains. Phenotypic assays further showed that these genetic changes were generally not associated with increased serum resistance, higher fitness in human blood ex vivo or differences in the interaction with human epithelial and endothelial cells in vitro. In conclusion, we hypothesize that virulence of meningococci results from accidental emergence of invasive variants during carriage and without within host evolution of invasive phenotypes during disease progression in vivo., Competing Interests: The authors have declared that no competing interests exist.

Hematopoietic stem cells give rise to all blood cells in a differentiation process that involves widespread epigenome remodeling. Here we present genome-wide reference maps of the associated DNA methylation dynamics. We used a meta-epigenomic approach that combines DNA methylation profiles across many small pools of cells and performed single-cell methylome sequencing to assess cell-to-cell heterogeneity. The resulting dataset identified characteristic differences between HSCs derived from fetal liver, cord blood, bone marrow, and peripheral blood. We also observed lineage-specific DNA methylation between myeloid and lymphoid progenitors, characterized immature multi-lymphoid progenitors, and detected progressive DNA methylation differences in maturing megakaryocytes. We linked these patterns to gene expression, histone modifications, and chromatin accessibility, and we used machine learning to derive a model of human hematopoietic differentiation directly from DNA methylation data. Our results contribute to a better understanding of human hematopoietic stem cell differentiation and provide a framework for studying blood-linked diseases., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

Tissue-resident macrophages support embryonic development and tissue homeostasis and repair. The mechanisms that control their differentiation remain unclear. We report here that erythro-myeloid progenitors in mice generate premacrophages (pMacs) that simultaneously colonize the whole embryo from embryonic day 9.5 in a chemokine-receptor-dependent manner. The core macrophage program initiated in pMacs is rapidly diversified as expression of transcriptional regulators becomes tissue-specific in early macrophages. This process appears essential for macrophage specification and maintenance, as inactivation of Id3 impairs the development of liver macrophages and results in selective Kupffer cell deficiency in adults. We propose that macrophage differentiation is an integral part of organogenesis, as colonization of organ anlagen by pMacs is followed by their specification into tissue macrophages, hereby generating the macrophage diversity observed in postnatal tissues., (Copyright © 2016, American Association for the Advancement of Science.)

Inflammatory bowel disease (IBD) arises by unknown environmental triggers in genetically susceptible individuals. Epigenetic regulation of gene expression may integrate internal and external influences and may thereby modulate disease susceptibility. Epigenetic modification may also affect the germ-line and in certain contexts can be inherited to offspring. This study investigates epigenetic alterations consequent to experimental murine colitis induced by dextran sodium sulphate (DSS), and their paternal transmission to offspring. Genome-wide methylome- and transcriptome-profiling of intestinal epithelial cells (IECs) and sperm cells of males of the F0 generation, which received either DSS and consequently developed colitis (F0(DSS)), or non-supplemented tap water (F0(Ctrl)) and hence remained healthy, and of their F1 offspring was performed using reduced representation bisulfite sequencing (RRBS) and RNA-sequencing (RNA-Seq), respectively. Offspring of F0(DSS) males exhibited aberrant methylation and expression patterns of multiple genes, including Igf1r and Nr4a2, which are involved in energy metabolism. Importantly, DSS colitis in F0(DSS) mice was associated with decreased body weight at baseline of their F1 offspring, and these F1 mice exhibited increased susceptibility to DSS-induced colitis compared to offspring from F0(Ctrl) males. This study hence demonstrates epigenetic transmissibility of metabolic and inflammatory traits resulting from experimental colitis.

Pancreatic islets of Langerhans contain several specialized endocrine cell types, which are commonly identified by the expression of single marker genes. However, the established marker genes cannot capture the complete spectrum of cellular heterogeneity in human pancreatic islets, and existing bulk transcriptome datasets provide averages across several cell populations. To dissect the cellular composition of the human pancreatic islet and to establish transcriptomes for all major cell types, we performed single-cell RNA sequencing on 70 cells sorted from human primary tissue. We used this dataset to validate previously described marker genes at the single-cell level and to identify specifically expressed transcription factors for all islet cell subtypes. All data are available for browsing and download, thus establishing a useful resource of single-cell expression profiles for endocrine cells in human pancreatic islets., (© 2015 The Authors. Published under the terms of the CC BY NC ND 4.0 license.)

Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for DNA methylation analysis in non-model organisms, we developed an integrated approach for studying DNA methylation differences independent of a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, carp, sea bass, and zebrafish). Bioinformatically, we developed the RefFreeDMA software to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions between samples or groups of individuals (http://RefFreeDMA.computational-epigenetics.org). The identified regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell-type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

Methods for single-cell genome and transcriptome sequencing have contributed to our understanding of cellular heterogeneity, whereas methods for single-cell epigenomics are much less established. Here, we describe a whole-genome bisulfite sequencing (WGBS) assay that enables DNA methylation mapping in very small cell populations (μWGBS) and single cells (scWGBS). Our assay is optimized for profiling many samples at low coverage, and we describe a bioinformatic method that analyzes collections of single-cell methylomes to infer cell-state dynamics. Using these technological advances, we studied epigenomic cell-state dynamics in three in vitro models of cellular differentiation and pluripotency, where we observed characteristic patterns of epigenome remodeling and cell-to-cell heterogeneity. The described method enables single-cell analysis of DNA methylation in a broad range of biological systems, including embryonic development, stem cell differentiation, and cancer. It can also be used to establish composite methylomes that account for cell-to-cell heterogeneity in complex tissue samples., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

Genetic control of male or female gonad development displays between different groups of organisms a remarkable diversity of "master sex-determining genes" at the top of the genetic hierarchies, whereas downstream components surprisingly appear to be evolutionarily more conserved. Without much further studies, conservation of sequence has been equalized to conservation of function. We have used the medaka fish to investigate the generality of this paradigm. In medaka, the master male sex-determining gene is dmrt1bY, a highly conserved downstream regulator of sex determination in vertebrates. To understand its function in orchestrating the complex gene regulatory network, we have identified targets genes and regulated pathways of Dmrt1bY. Monitoring gene expression and interactions by transgenic fluorescent reporter fish lines, in vivo tissue-chromatin immunoprecipitation and in vitro gene regulation assays revealed concordance but also major discrepancies between mammals and medaka, notably amongst spatial, temporal expression patterns and regulations of the canonical Hedgehog and R-spondin/Wnt/Follistatin signaling pathways. Examination of Foxl2 protein distribution in the medaka ovary defined a new subpopulation of theca cells, where ovarian-type aromatase transcriptional regulation appears to be independent of Foxl2. In summary, these data show that the regulation of the downstream regulatory network of sex determination is less conserved than previously thought.