The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ∼30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.

The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) project received funding from the European Commission for 41/2 years. The objectives of the project revolved around technological developments in nucleic acid analysis. The project partners have discovered, created and developed a huge body of insights into nucleic acid analysis, ranging from improvements and implementation of current technologies to the most promising sequencing technologies that constitute a 3(rd) and 4(th) generation of sequencing methods with nanopores and in situ sequencing, respectively.

Exome sequencing has become a popular method to evaluate undirected mutagenesis experiments in mice. However, the most suitable mouse strain for the biological model may be relatively distant from the standard mouse reference genome. For pinpointing causative variants, a matching reference with gene annotations is essential, but not always readily available.

Cactophilic Drosophila species provide a valuable model to study gene-environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii-D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii-D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

T-cell acute lymphoblastic leukaemias (T-ALL) are aggressive malignant proliferations characterized by high relapse rates and great genetic heterogeneity. TAL1 is amongst the most frequently deregulated oncogenes. Yet, over half of the TAL1(+) cases lack TAL1 lesions, suggesting unrecognized (epi)genetic deregulation mechanisms. Here we show that TAL1 is normally silenced in the T-cell lineage, and that the polycomb H3K27me3-repressive mark is focally diminished in TAL1(+) T-ALLs. Sequencing reveals that >20% of monoallelic TAL1(+) patients without previously known alterations display microinsertions or RAG1/2-mediated episomal reintegration in a single site 5' to TAL1. Using 'allelic-ChIP' and CrispR assays, we demonstrate that such insertions induce a selective switch from H3K27me3 to H3K27ac at the inserted but not the germline allele. We also show that, despite a considerable mechanistic diversity, the mode of oncogenic TAL1 activation, rather than expression levels, impact on clinical outcome. Altogether, these studies establish site-specific epigenetic desilencing as a mechanism of oncogenic activation.

Here, we report DNA methylation and hydroxymethylation dynamics at nucleotide resolution using C/EBPα-enhanced reprogramming of B cells into induced pluripotent cells (iPSCs). We observed successive waves of hydroxymethylation at enhancers, concomitant with a decrease in DNA methylation, suggesting active demethylation. Consistent with this finding, ablation of the DNA demethylase Tet2 almost completely abolishes reprogramming. C/EBPα, Klf4, and Tfcp2l1 each interact with Tet2 and recruit the enzyme to specific DNA sites. During reprogramming, some of these sites maintain high levels of 5hmC, and enhancers and promoters of key pluripotency factors become demethylated as early as 1 day after Yamanaka factor induction. Surprisingly, methylation changes precede chromatin opening in distinct chromatin regions, including Klf4 bound sites, revealing a pioneer factor activity associated with alternation in DNA methylation. Rapid changes in hydroxymethylation similar to those in B cells were also observed during compound-accelerated reprogramming of fibroblasts into iPSCs, highlighting the generality of our observations.

The molecular mechanisms underlying specification from embryonic stem cells (ESCs) and maintenance of neural progenitor cells (NPCs) are largely unknown. Recently, we reported that the Zuotin-related factor 1 (Zrf1) is necessary for chromatin displacement of the Polycomb-repressive complex 1 (PRC1). We found that Zrf1 is required for NPC specification from ESCs and that it promotes the expression of NPC markers, including the key regulator Pax6. Moreover, Zrf1 is essential to establish and maintain Wnt ligand expression levels, which are necessary for NPC self-renewal. Reactivation of proper Wnt signaling in Zrf1-depleted NPCs restores Pax6 expression and the self-renewal capacity. ESC-derived NPCs in vitro resemble most of the characteristics of the self-renewing NPCs located in the developing embryonic cortex, which are termed radial glial cells (RGCs). Depletion of Zrf1 in vivo impairs the expression of key self-renewal regulators and Wnt ligand genes in RGCs. Thus, we demonstrate that Zrf1 plays an essential role in NPC generation and maintenance.

Shank3/PROSAP2 gene mutations are associated with cognitive impairment ranging from mental retardation to autism. Shank3 is a large scaffold postsynaptic density protein implicated in dendritic spines and synapse formation; however, its specific functions have not been clearly demonstrated. We have used RNAi to knockdown Shank3 expression in neuronal cultures and showed that this treatment specifically reduced the synaptic expression of the metabotropic glutamate receptor 5 (mGluR5), but did not affect the expression of other major synaptic proteins. The functional consequence of Shank3 RNAi knockdown was impaired signaling via mGluR5, as shown by reduction in ERK1/2 and CREB phosphorylation induced by stimulation with (S)-3,5-dihydroxyphenylglycine (DHPG) as the agonist of mGluR5 receptors, impaired mGluR5-dependent synaptic plasticity (DHPG-induced long-term depression), and impaired mGluR5-dependent modulation of neural network activity. We also found morphological abnormalities in the structure of synapses (spine number, width, and length) and impaired glutamatergic synaptic transmission, as shown by reduction in the frequency of miniature excitatory postsynaptic currents (mEPSC). Notably, pharmacological augmentation of mGluR5 activity using 3-cyano-N-(1,3-diphenyl-1H-pyrazol-5-yl)-benzamide as the positive allosteric modulator of these receptors restored mGluR5-dependent signaling (DHPG-induced phosphorylation of ERK1/2) and normalized the frequency of mEPSCs in Shank3-knocked down neurons. These data demonstrate that a deficit in mGluR5-mediated intracellular signaling in Shank3 knockdown neurons can be compensated by 3-cyano-N-(1,3-diphenyl-1H-pyrazol-5-yl)-benzamide; this raises the possibility that pharmacological augmentation of mGluR5 activity represents a possible new therapeutic approach for patients with Shank3 mutations.

Six extant species of non-human great apes are currently recognized: Sumatran and Bornean orangutans, eastern and western gorillas, and chimpanzees and bonobos [1]. However, large gaps remain in our knowledge of fine-scale variation in hominoid morphology, behavior, and genetics, and aspects of great ape taxonomy remain in flux. This is particularly true for orangutans (genus: Pongo), the only Asian great apes and phylogenetically our most distant relatives among extant hominids [1]. Designation of Bornean and Sumatran orangutans, P. pygmaeus (Linnaeus 1760) and P. abelii (Lesson 1827), as distinct species occurred in 2001 [1, 2]. Here, we show that an isolated population from Batang Toru, at the southernmost range limit of extant Sumatran orangutans south of Lake Toba, is distinct from other northern Sumatran and Bornean populations. By comparing cranio-mandibular and dental characters of an orangutan killed in a human-animal conflict to those of 33 adult male orangutans of a similar developmental stage, we found consistent differences between the Batang Toru individual and other extant Ponginae. Our analyses of 37 orangutan genomes provided a second line of evidence. Model-based approaches revealed that the deepest split in the evolutionary history of extant orangutans occurred ∼3.38 mya between the Batang Toru population and those to the north of Lake Toba, whereas both currently recognized species separated much later, about 674 kya. Our combined analyses support a new classification of orangutans into three extant species. The new species, Pongo tapanuliensis, encompasses the Batang Toru population, of which fewer than 800 individuals survive. VIDEO ABSTRACT.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .

Somatic cell reprogramming is the process that allows differentiated cells to revert to a pluripotent state. In contrast to the extensively studied rewiring of epigenetic and transcriptional programs required for reprogramming, the dynamics of post-transcriptional changes and their associated regulatory mechanisms remain poorly understood. Here we study the dynamics of alternative splicing changes occurring during efficient reprogramming of mouse B cells into induced pluripotent stem (iPS) cells and compare them to those occurring during reprogramming of mouse embryonic fibroblasts.

Determining the effect of DNA methylation on chromatin structure and function in higher organisms is challenging due to the extreme complexity of epigenetic regulation. We studied a simpler model system, budding yeast, that lacks DNA methylation machinery making it a perfect model system to study the intrinsic role of DNA methylation in chromatin structure and function. We expressed the murine DNA methyltransferases in Saccharomyces cerevisiae and analyzed the correlation between DNA methylation, nucleosome positioning, gene expression and 3D genome organization. Despite lacking the machinery for positioning and reading methylation marks, induced DNA methylation follows a conserved pattern with low methylation levels at the 5' end of the gene increasing gradually toward the 3' end, with concentration of methylated DNA in linkers and nucleosome free regions, and with actively expressed genes showing low and high levels of methylation at transcription start and terminating sites respectively, mimicking the patterns seen in mammals. We also see that DNA methylation increases chromatin condensation in peri-centromeric regions, decreases overall DNA flexibility, and favors the heterochromatin state. Taken together, these results demonstrate that methylation intrinsically modulates chromatin structure and function even in the absence of cellular machinery evolved to recognize and process the methylation signal.

Here, we describe how the speed of C/EBPα-induced B cell to macrophage transdifferentiation (BMT) can be regulated, using both mouse and human models. The identification of a mutant of C/EBPα (C/EBPαR35A) that greatly accelerates BMT helped to illuminate the mechanism. Thus, incoming C/EBPα binds to PU.1, an obligate partner expressed in B cells, leading to the release of PU.1 from B cell enhancers, chromatin closing and silencing of the B cell program. Released PU.1 redistributes to macrophage enhancers newly occupied by C/EBPα, causing chromatin opening and activation of macrophage genes. All these steps are accelerated by C/EBPαR35A, initiated by its increased affinity for PU.1. Wild-type C/EBPα is methylated by Carm1 at arginine 35 and the enzyme's perturbations modulate BMT velocity as predicted from the observations with the mutant. Increasing the proportion of unmethylated C/EBPα in granulocyte/macrophage progenitors by inhibiting Carm1 biases the cell's differentiation toward macrophages, suggesting that cell fate decision velocity and lineage directionality are closely linked processes.

Patient-derived 3D cell culture systems are currently advancing cancer research since they potentiate the molecular analysis of tissue-like properties and drug response under well-defined conditions. However, our understanding of the relationship between the heterogeneity of morphological phenotypes and the underlying transcriptome is still limited. To address this issue, we here introduce "pheno-seq" to directly link visual features of 3D cell culture systems with profiling their transcriptome. As prototypic applications breast and colorectal cancer (CRC) spheroids were analyzed by pheno-seq. We identified characteristic gene expression signatures of epithelial-to-mesenchymal transition that are associated with invasive growth behavior of clonal breast cancer spheroids. Furthermore, we linked long-term proliferative capacity in a patient-derived model of CRC to a lowly abundant PROX1-positive cancer stem cell subtype. We anticipate that the ability to integrate transcriptome analysis and morphological patho-phenotypes of cancer cells will provide novel insight on the molecular origins of intratumor heterogeneity.

TWIST1 is thought to be a novel oncogene. Understanding the molecular mechanisms regulating the TWIST1 gene expression profiles in tumor cells may give new insights regarding prognostic factors and novel therapeutic targets in veterinary oncology. In the present study we partially isolated the TWIST1 gene in Felis catus and performed comparative studies. Several primer combinations were used based on the alignments of homologous DNA sequences. After PCR amplification, three bands were obtained, purified and sequenced. Several bioinformatic tools were utilized to carry out the comparative studies. Higher similarity was found between the isolated TWIST1 gene in Felis catus and Homo sapiens (86%) than between Homo sapiens and Rattus norvegicus or Mus musculus (75%). Partial amino acid sequence showed no change in the four species analyzed. This confirmed that coding sequences presented high similarity (~96%) between man and cat. These results give the first insights regarding the TWIST1 gene in cat but further studies are required in order to establish, or not, its role in tumor formation and progression in veterinary oncology.

Forced transcription factor expression can transdifferentiate somatic cells into other specialised cell types or reprogram them into induced pluripotent stem cells (iPSCs) with variable efficiency. To better understand the heterogeneity of these processes, we used single-cell RNA sequencing to follow the transdifferentation of murine pre-B cells into macrophages as well as their reprogramming into iPSCs. Even in these highly efficient systems, there was substantial variation in the speed and path of fate conversion. We predicted and validated that these differences are inversely coupled and arise in the starting cell population, with Mychigh large pre-BII cells transdifferentiating slowly but reprogramming efficiently and Myclow small pre-BII cells transdifferentiating rapidly but failing to reprogram. Strikingly, differences in Myc activity predict the efficiency of reprogramming across a wide range of somatic cell types. These results illustrate how single cell expression and computational analyses can identify the origins of heterogeneity in cell fate conversion processes.

The concept of tissue-specific gene expression posits that lineage-determining transcription factors (LDTFs) determine the open chromatin profile of a cell via collaborative binding, providing molecular beacons to signal-dependent transcription factors (SDTFs). However, the guiding principles of LDTF binding, chromatin accessibility and enhancer activity have not yet been systematically evaluated. We sought to study these features of the macrophage genome by the combination of experimental (ChIP-seq, ATAC-seq and GRO-seq) and computational approaches. We show that Random Forest and Support Vector Regression machine learning methods can accurately predict chromatin accessibility using the binding patterns of the LDTF PU.1 and four other key TFs of macrophages (IRF8, JUNB, CEBPA and RUNX1). Any of these TFs alone were not sufficient to predict open chromatin, indicating that TF binding is widespread at closed or weakly opened chromatin regions. Analysis of the PU.1 cistrome revealed that two-thirds of PU.1 binding occurs at low accessible chromatin. We termed these sites labelled regulatory elements (LREs), which may represent a dormant state of a future enhancer and contribute to macrophage cellular plasticity. Collectively, our work demonstrates the existence of LREs occupied by various key TFs, regulating specific gene expression programs triggered by divergent macrophage polarizing stimuli.

Forced expression of reprogramming factors can convert somatic cells into induced pluripotent stem cells (iPSCs). Here we studied genome topology dynamics during reprogramming of different somatic cell types with highly distinct genome conformations. We find large-scale topologically associated domain (TAD) repositioning and alterations of tissue-restricted genomic neighborhoods and chromatin loops, effectively erasing the somatic-cell-specific genome structures while establishing an embryonic stem-cell-like 3D genome. Yet, early passage iPSCs carry topological hallmarks that enable recognition of their cell of origin. These hallmarks are not remnants of somatic chromosome topologies. Instead, the distinguishing topological features are acquired during reprogramming, as we also find for cell-of-origin-dependent gene expression patterns.

While analyzing the DNA methylome of multiple myeloma (MM), a plasma cell neoplasm, by whole-genome bisulfite sequencing and high-density arrays, we observed a highly heterogeneous pattern globally characterized by regional DNA hypermethylation embedded in extensive hypomethylation. In contrast to the widely reported DNA hypermethylation of promoter-associated CpG islands (CGIs) in cancer, hypermethylated sites in MM, as opposed to normal plasma cells, were located outside CpG islands and were unexpectedly associated with intronic enhancer regions defined in normal B cells and plasma cells. Both RNA-seq and in vitro reporter assays indicated that enhancer hypermethylation is globally associated with down-regulation of its host genes. ChIP-seq and DNase-seq further revealed that DNA hypermethylation in these regions is related to enhancer decommissioning. Hypermethylated enhancer regions overlapped with binding sites of B cell-specific transcription factors (TFs) and the degree of enhancer methylation inversely correlated with expression levels of these TFs in MM. Furthermore, hypermethylated regions in MM were methylated in stem cells and gradually became demethylated during normal B-cell differentiation, suggesting that MM cells either reacquire epigenetic features of undifferentiated cells or maintain an epigenetic signature of a putative myeloma stem cell progenitor. Overall, we have identified DNA hypermethylation of developmentally regulated enhancers as a new type of epigenetic modification associated with the pathogenesis of MM.

In conformational disorders, it is not evident which amyloid aggregates affect specific molecular mechanisms or cellular pathways, which cause disease because of their quantity and mechanical features and which states in aggregate formation are pathogenic. Due to the increasing consensus that prefibrillar oligomers play a major role in conformational diseases, there is a growing interest in understanding the characteristics of metastable polypeptide associations.