Direct neuronal reprogramming of mouse astrocytes is associated with multiscale epigenome remodeling and requires Yy1.

Direct neuronal reprogramming is a promising approach to regenerate neurons from local glial cells. However, mechanisms of epigenome remodeling and co-factors facilitating this process are unclear. In this study, we combined single-cell multiomics with genome-wide profiling of three-dimensional nuclear architecture and DNA methylation in mouse astrocyte-to-neuron reprogramming mediated by Neurogenin2 (Ngn2) and its phosphorylation-resistant form (PmutNgn2), respectively. We show that Ngn2 drives multilayered chromatin remodeling at dynamic enhancer-gene interaction sites. PmutNgn2 leads to higher reprogramming efficiency and enhances epigenetic remodeling associated with neuronal maturation. However, the differences in binding sites or downstream gene activation cannot fully explain this effect. Instead, we identified Yy1, a transcriptional co-factor recruited by direct interaction with Ngn2 to its target sites. Upon deletion of Yy1, activation of neuronal enhancers, genes and ultimately reprogramming are impaired without affecting Ngn2 binding. Thus, our work highlights the key role of interactors of proneural factors in direct neuronal reprogramming.

Single-cell metabolic profiling reveals subgroups of primary human hepatocytes with heterogeneous responses to drug challenge.

BACKGROUND: Xenobiotics are primarily metabolized by hepatocytes in the liver, and primary human hepatocytes are the gold standard model for the assessment of drug efficacy, safety, and toxicity in the early phases of drug development. Recent advances in single-cell genomics demonstrate liver zonation and ploidy as main drivers of cellular heterogeneity. However, little is known about the impact of hepatocyte specialization on liver function upon metabolic challenge, including hepatic metabolism, detoxification, and protein synthesis. RESULTS: Here, we investigate the metabolic capacity of individual human hepatocytes in vitro. We assess how chronic accumulation of lipids enhances cellular heterogeneity and impairs the metabolisms of drugs. Using a phenotyping five-probe cocktail, we identify four functional subgroups of hepatocytes responding differently to drug challenge and fatty acid accumulation. These four subgroups display differential gene expression profiles upon cocktail treatment and xenobiotic metabolism-related specialization. Notably, intracellular fat accumulation leads to increased transcriptional variability and diminishes the drug-related metabolic capacity of hepatocytes. CONCLUSIONS: Our results demonstrate that, upon a metabolic challenge such as exposure to drugs or intracellular fat accumulation, hepatocyte subgroups display different and heterogeneous transcriptional responses.

CSF3R T618I Collaborates With RUNX1-RUNX1T1 to Expand Hematopoietic Progenitors and Sensitizes to GLI Inhibition.

Activating colony-stimulating factor-3 receptor gene (CSF3R) mutations are recurrent in acute myeloid leukemia (AML) with t(8;21) translocation. However, the nature of oncogenic collaboration between alterations of CSF3R and the t(8;21) associated RUNX1-RUNX1T1 fusion remains unclear. In CD34+ hematopoietic stem and progenitor cells from healthy donors, double oncogene expression led to a clonal advantage, increased self-renewal potential, and blast-like morphology and distinct immunophenotype. Gene expression profiling revealed hedgehog signaling as a potential mechanism, with upregulation of GLI2 constituting a putative pharmacological target. Both primary hematopoietic cells and the t(8;21) positive AML cell line SKNO-1 showed increased sensitivity to the GLI inhibitor GANT61 when expressing CSF3R T618I. Our findings suggest that during leukemogenesis, the RUNX1-RUNXT1 fusion and CSF3R mutation act in a synergistic manner to alter hedgehog signaling, which can be exploited therapeutically.

epiAneufinder identifies copy number alterations from single-cell ATAC-seq data.

Single-cell open chromatin profiling via scATAC-seq has become a mainstream measurement of open chromatin in single-cells. Here we present epiAneufinder, an algorithm that exploits the read count information from scATAC-seq data to extract genome-wide copy number alterations (CNAs) for individual cells, allowing the study of CNA heterogeneity present in a sample at the single-cell level. Using different cancer scATAC-seq datasets, we show that epiAneufinder can identify intratumor clonal heterogeneity in populations of single cells based on their CNA profiles. We demonstrate that these profiles are concordant with the ones inferred from single-cell whole genome sequencing data for the same samples. EpiAneufinder allows the inference of single-cell CNA information from scATAC-seq data, without the need of additional experiments, unlocking a layer of genomic variation which is otherwise unexplored.

CIARA: a cluster-independent algorithm for identifying markers of rare cell types from single-cell sequencing data.

A powerful feature of single-cell genomics is the possibility of identifying cell types from their molecular profiles. In particular, identifying novel rare cell types and their marker genes is a key potential of single-cell RNA sequencing. Standard clustering approaches perform well in identifying relatively abundant cell types, but tend to miss rarer cell types. Here, we have developed CIARA (Cluster Independent Algorithm for the identification of markers of RAre cell types), a cluster-independent computational tool designed to select genes that are likely to be markers of rare cell types. Genes selected by CIARA are subsequently integrated with common clustering algorithms to single out groups of rare cell types. CIARA outperforms existing methods for rare cell type detection, and we use it to find previously uncharacterized rare populations of cells in a human gastrula and among mouse embryonic stem cells treated with retinoic acid. Moreover, CIARA can be applied more generally to any type of single-cell omic data, thus allowing the identification of rare cells across multiple data modalities. We provide implementations of CIARA in user-friendly packages available in R and Python.

EpiScanpy: integrated single-cell epigenomic analysis.

EpiScanpy is a toolkit for the analysis of single-cell epigenomic data, namely single-cell DNA methylation and single-cell ATAC-seq data. To address the modality specific challenges from epigenomics data, epiScanpy quantifies the epigenome using multiple feature space constructions and builds a nearest neighbour graph using epigenomic distance between cells. EpiScanpy makes the many existing scRNA-seq workflows from scanpy available to large-scale single-cell data from other -omics modalities, including methods for common clustering, dimension reduction, cell type identification and trajectory learning techniques, as well as an atlas integration tool for scATAC-seq datasets. The toolkit also features numerous useful downstream functions, such as differential methylation and differential openness calling, mapping epigenomic features of interest to their nearest gene, or constructing gene activity matrices using chromatin openness. We successfully benchmark epiScanpy against other scATAC-seq analysis tools and show its outperformance at discriminating cell types.

Brief homogeneous TCR signals instruct common iNKT progenitors whose effector diversification is characterized by subsequent cytokine signaling.

Innate-like T cell populations expressing conserved TCRs play critical roles in immunity through diverse developmentally acquired effector functions. Focusing on the prototypical lineage of invariant natural killer T (iNKT) cells, we sought to dissect the mechanisms and timing of fate decisions and functional effector differentiation. Utilizing induced expression of the semi-invariant NKT cell TCR on double positive thymocytes, an initially highly synchronous wave of iNKT cell development was triggered by brief homogeneous TCR signaling. After reaching a uniform progenitor state characterized by IL-4 production potential and proliferation, effector subsets emerged simultaneously, but then diverged toward different fates. While NKT17 specification was quickly completed, NKT1 cells slowly differentiated and expanded. NKT2 cells resembled maturing progenitors, which gradually diminished in numbers. Thus, iNKT subset diversification occurs in dividing progenitor cells without acute TCR input but utilizes multiple active cytokine signaling pathways. These data imply a two-step model of iNKT effector differentiation.

Region-level epimutation rates in Arabidopsis thaliana.

Failure to maintain DNA methylation patterns during plant development can occasionally give rise to so-called "spontaneous epimutations". These stochastic methylation changes are sometimes heritable across generations and thus accumulate in plant genomes over time. Recent evidence indicates that spontaneous epimutations have a major role in shaping patterns of methylation diversity in plant populations. Using single CG dinucleotides as units of analysis, previous work has shown that the epimutation rate is several orders of magnitude higher than the genetic mutation rate. While these large rate differences have obvious implications for understanding genome-methylome co-evolution, the functional relevance of single CG methylation changes remains questionable. In contrast to single CG, solid experimental evidence has linked methylation gains and losses in larger genomic regions with transcriptional variation and heritable phenotypic effects. Here we show that such region-level changes arise stochastically at about the same rate as those at individual CG sites, are only marginal dependent on region size and cytosine density, but strongly dependent on chromosomal location. We also find consistent evidence that region-level epimutations are not restricted to CG contexts but also frequently occur in non-CG regions at the genome-wide scale. Taken together, our results support the view that many differentially methylated regions (DMRs) in natural populations originate from epimutation events and may not be effectively tagged by proximal SNPs. This possibility reinforces the need for epigenome-wide association studies (EWAS) in plants as a way to identify the epigenetic basis of complex traits.

AlphaBeta: computational inference of epimutation rates and spectra from high-throughput DNA methylation data in plants.

Stochastic changes in DNA methylation (i.e., spontaneous epimutations) contribute to methylome diversity in plants. Here, we describe AlphaBeta, a computational method for estimating the precise rate of such stochastic events using pedigree-based DNA methylation data as input. We demonstrate how AlphaBeta can be employed to study transgenerationally heritable epimutations in clonal or sexually derived mutation accumulation lines, as well as somatic epimutations in long-lived perennials. Application of our method to published and new data reveals that spontaneous epimutations accumulate neutrally at the genome-wide scale, originate mainly during somatic development and that they can be used as a molecular clock for age-dating trees.

A genome assembly and the somatic genetic and epigenetic mutation rate in a wild long-lived perennial Populus trichocarpa.

BACKGROUND: Plants can transmit somatic mutations and epimutations to offspring, which in turn can affect fitness. Knowledge of the rate at which these variations arise is necessary to understand how plant development contributes to local adaption in an ecoevolutionary context, particularly in long-lived perennials. RESULTS: Here, we generate a new high-quality reference genome from the oldest branch of a wild Populus trichocarpa tree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we use a multi-omics approach to quantify age-related somatic changes at the genetic, epigenetic, and transcriptional level. We show that the per-year somatic mutation and epimutation rates are lower than in annuals and that transcriptional variation is mainly independent of age divergence and cytosine methylation. Furthermore, a detailed analysis of the somatic epimutation spectrum indicates that transgenerationally heritable epimutations originate mainly from DNA methylation maintenance errors during mitotic rather than during meiotic cell divisions. CONCLUSION: Taken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

Author Correction: Widespread long-range cis-regulatory elements in the maize genome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Altering microtubule dynamics is synergistically toxic with spindle assembly checkpoint inhibition.

Chromosomal instability (CIN) and aneuploidy are hallmarks of cancer. As most cancers are aneuploid, targeting aneuploidy or CIN may be an effective way to target a broad spectrum of cancers. Here, we perform two small molecule compound screens to identify drugs that selectively target cells that are aneuploid or exhibit a CIN phenotype. We find that aneuploid cells are much more sensitive to the energy metabolism regulating drug ZLN005 than their euploid counterparts. Furthermore, cells with an ongoing CIN phenotype, induced by spindle assembly checkpoint (SAC) alleviation, are significantly more sensitive to the Src kinase inhibitor SKI606. We show that inhibiting Src kinase increases microtubule polymerization rates and, more generally, that deregulating microtubule polymerization rates is particularly toxic to cells with a defective SAC. Our findings, therefore, suggest that tumors with a dysfunctional SAC are particularly sensitive to microtubule poisons and, vice versa, that compounds alleviating the SAC provide a powerful means to treat tumors with deregulated microtubule dynamics.

breakpointR: an R/Bioconductor package to localize strand state changes in Strand-seq data.

MOTIVATION: Strand-seq is a specialized single-cell DNA sequencing technique centered around the directionality of single-stranded DNA. Computational tools for Strand-seq analyses must capture the strand-specific information embedded in these data. RESULTS: Here we introduce breakpointR, an R/Bioconductor package specifically tailored to process and interpret single-cell strand-specific sequencing data obtained from Strand-seq. We developed breakpointR to detect local changes in strand directionality of aligned Strand-seq data, to enable fine-mapping of sister chromatid exchanges, germline inversion and to support global haplotype assembly. Given the broad spectrum of Strand-seq applications we expect breakpointR to be an important addition to currently available tools and extend the accessibility of this novel sequencing technique. AVAILABILITY AND IMPLEMENTATION: R/Bioconductor package https://bioconductor.org/packages/breakpointR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gene networks in cancer are biased by aneuploidies and sample impurities.

Gene regulatory network inference is a standard technique for obtaining structured regulatory information from, for instance, gene expression measurements. Methods performing this task have been extensively evaluated on synthetic, and to a lesser extent real data sets. In contrast to these test evaluations, applications to gene expression data of human cancers are often limited by fewer samples and more potential regulatory links, and are biased by copy number aberrations as well as cell mixtures and sample impurities. Here, we take networks inferred from TCGA cohorts as an example to show that (1) transcription factor annotations are essential to obtain reliable networks, and (2) even for state of the art methods, we expect that between 20 and 80% of edges are caused by copy number changes and cell mixtures rather than transcription factor regulation.

Widespread long-range cis-regulatory elements in the maize genome.

Genetic mapping studies on crops suggest that agronomic traits can be controlled by gene-distal intergenic loci. Despite the biological importance and the potential agronomic utility of these loci, they remain virtually uncharacterized in all crop species to date. Here, we provide genetic, epigenomic and functional molecular evidence to support the widespread existence of gene-distal (hereafter, distal) loci that act as long-range transcriptional cis-regulatory elements (CREs) in the maize genome. Such loci are enriched for euchromatic features that suggest their regulatory functions. Chromatin loops link together putative CREs with genes and recapitulate genetic interactions. Putative CREs also display elevated transcriptional enhancer activities, as measured by self-transcribing active regulatory region sequencing. These results provide functional support for the widespread existence of CREs that act over large genomic distances to control gene expression.

Endogenous retroviral insertions drive non-canonical imprinting in extra-embryonic tissues.

BACKGROUND: Genomic imprinting is an epigenetic phenomenon that allows a subset of genes to be expressed mono-allelically based on the parent of origin and is typically regulated by differential DNA methylation inherited from gametes. Imprinting is pervasive in murine extra-embryonic lineages, and uniquely, the imprinting of several genes has been found to be conferred non-canonically through maternally inherited repressive histone modification H3K27me3. However, the underlying regulatory mechanisms of non-canonical imprinting in post-implantation development remain unexplored. RESULTS: We identify imprinted regions in post-implantation epiblast and extra-embryonic ectoderm (ExE) by assaying allelic histone modifications (H3K4me3, H3K36me3, H3K27me3), gene expression, and DNA methylation in reciprocal C57BL/6 and CAST hybrid embryos. We distinguish loci with DNA methylation-dependent (canonical) and independent (non-canonical) imprinting by assaying hybrid embryos with ablated maternally inherited DNA methylation. We find that non-canonical imprints are localized to endogenous retrovirus-K (ERVK) long terminal repeats (LTRs), which act as imprinted promoters specifically in extra-embryonic lineages. Transcribed ERVK LTRs are CpG-rich and located in close proximity to gene promoters, and imprinting status is determined by their epigenetic patterning in the oocyte. Finally, we show that oocyte-derived H3K27me3 associated with non-canonical imprints is not maintained beyond pre-implantation development at these elements and is replaced by secondary imprinted DNA methylation on the maternal allele in post-implantation ExE, while being completely silenced by bi-allelic DNA methylation in the epiblast. CONCLUSIONS: This study reveals distinct epigenetic mechanisms regulating non-canonical imprinted gene expression between embryonic and extra-embryonic development and identifies an integral role for ERVK LTR repetitive elements.

Quantification of Aneuploidy in Mammalian Systems.

High-throughput next generation sequencing karyotyping has emerged as a powerful tool for the detection of genomic heterogeneity in normal tissues and cancers. Here we describe a single-cell whole genome sequencing (scWGS) platform to assess whole-chromosome aneuploidy, structural aneuploidies involving only chromosome fragments and more local small copy number alterations in individual cells. We provide a detailed protocol for the isolation, library preparation, low coverage sequencing and data analysis of single cells. Since our approach does not involve a whole-genome preamplification step, our method allows for acquisition of reliable high-resolution single-cell copy number profiles. Moreover, the protocol allows multiplexing of 384 single-cell libraries in one sequencing run, thereby significantly reducing sequencing costs and can be completed in 3-4 days starting from single cell isolation to analysis of sequencing data.

METHimpute: imputation-guided construction of complete methylomes from WGBS data.

BACKGROUND: Whole-genome bisulfite sequencing (WGBS) has become the standard method for interrogating plant methylomes at base resolution. However, deep WGBS measurements remain cost prohibitive for large, complex genomes and for population-level studies. As a result, most published plant methylomes are sequenced far below saturation, with a large proportion of cytosines having either missing data or insufficient coverage. RESULTS: Here we present METHimpute, a Hidden Markov Model (HMM) based imputation algorithm for the analysis of WGBS data. Unlike existing methods, METHimpute enables the construction of complete methylomes by inferring the methylation status and level of all cytosines in the genome regardless of coverage. Application of METHimpute to maize, rice and Arabidopsis shows that the algorithm infers cytosine-resolution methylomes with high accuracy from data as low as 6X, compared to data with 60X, thus making it a cost-effective solution for large-scale studies. CONCLUSIONS: METHimpute provides methylation status calls and levels for all cytosines in the genome regardless of coverage, thus yielding complete methylomes even with low-coverage WGBS datasets. The method has been extensively tested in plants, but should also be applicable to other species. An implementation is available on Bioconductor.

MLL2 conveys transcription-independent H3K4 trimethylation in oocytes.

Histone 3 K4 trimethylation (depositing H3K4me3 marks) is typically associated with active promoters yet paradoxically occurs at untranscribed domains. Research to delineate the mechanisms of targeting H3K4 methyltransferases is ongoing. The oocyte provides an attractive system to investigate these mechanisms, because extensive H3K4me3 acquisition occurs in nondividing cells. We developed low-input chromatin immunoprecipitation to interrogate H3K4me3, H3K27ac and H3K27me3 marks throughout oogenesis. In nongrowing oocytes, H3K4me3 was restricted to active promoters, but as oogenesis progressed, H3K4me3 accumulated in a transcription-independent manner and was targeted to intergenic regions, putative enhancers and silent H3K27me3-marked promoters. Ablation of the H3K4 methyltransferase gene Mll2 resulted in loss of transcription-independent H3K4 trimethylation but had limited effects on transcription-coupled H3K4 trimethylation or gene expression. Deletion of Dnmt3a and Dnmt3b showed that DNA methylation protects regions from acquiring H3K4me3. Our findings reveal two independent mechanisms of targeting H3K4me3 to genomic elements, with MLL2 recruited to unmethylated CpG-rich regions independently of transcription.

Ecological plant epigenetics: Evidence from model and non-model species, and the way forward.

Growing evidence shows that epigenetic mechanisms contribute to complex traits, with implications across many fields of biology. In plant ecology, recent studies have attempted to merge ecological experiments with epigenetic analyses to elucidate the contribution of epigenetics to plant phenotypes, stress responses, adaptation to habitat, and range distributions. While there has been some progress in revealing the role of epigenetics in ecological processes, studies with non-model species have so far been limited to describing broad patterns based on anonymous markers of DNA methylation. In contrast, studies with model species have benefited from powerful genomic resources, which contribute to a more mechanistic understanding but have limited ecological realism. Understanding the significance of epigenetics for plant ecology requires increased transfer of knowledge and methods from model species research to genomes of evolutionarily divergent species, and examination of responses to complex natural environments at a more mechanistic level. This requires transforming genomics tools specifically for studying non-model species, which is challenging given the large and often polyploid genomes of plants. Collaboration among molecular geneticists, ecologists and bioinformaticians promises to enhance our understanding of the mutual links between genome function and ecological processes.