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.

Enhanced copper-resistance gene repertoire in Alteromonas macleodii strains isolated from copper-treated marine coatings.

Copper is prevalent in coastal ecosystems due to its use as an algaecide and as an anti-fouling agent on ship hulls. Alteromonas spp. have previously been shown to be some of the early colonizers of copper-based anti-fouling paint but little is known about the mechanisms they use to overcome this initial copper challenge. The main models of copper resistance include the Escherichia coli chromosome-based Cue and Cus systems; the plasmid-based E. coli Pco system; and the plasmid-based Pseudomonas syringae Cop system. These were all elucidated from strains isolated from copper-rich environments of agricultural and/or enteric origin. In this work, copper resistance assays demonstrated the ability of Alteromonas macleodii strains CUKW and KCC02 to grow at levels lethal to other marine bacterial species. A custom database of Hidden Markov Models was designed based on proteins from the Cue, Cus, and Cop/Pco systems and used to identify potential copper resistance genes in CUKW and KCC02. Comparative genomic analyses with marine bacterial species and bacterial species isolated from copper-rich environments demonstrated that CUKW and KCC02 possess genetic elements of all systems, oftentimes with multiple copies, distributed throughout the chromosome and mega-plasmids. In particular, two copies of copA (the key player in cytoplasmic detoxification), each with its own apparent MerR-like transcriptional regulator, occur on a mega-plasmid, along with multiple copies of Pco homologs. Genes from both systems were induced upon exposure to elevated copper levels (100 µM- 3 mM). Genomic analysis identified one of the merR-copA clusters occurs on a genomic island (GI) within the plasmid, and comparative genomic analysis found that either of the merR-copA clusters, which also includes genes coding for a cupredoxin domain-containing protein and an isoprenylcysteine methyltransferase, occurs on a GI across diverse bacterial species. These genomic findings combined with the ability of CUKW and KCC02 to grow in copper-challenged conditions are couched within the context of the genome flexibility of the Alteromonas genus.

Retinoic acid signaling is critical during the totipotency window in early mammalian development.

Totipotent cells hold enormous potential for regenerative medicine. Thus, the development of cellular models recapitulating totipotent-like features is of paramount importance. Cells resembling the totipotent cells of early embryos arise spontaneously in mouse embryonic stem (ES) cell cultures. Such '2-cell-like-cells' (2CLCs) recapitulate 2-cell-stage features and display expanded cell potential. Here, we used 2CLCs to perform a small-molecule screen to identify new pathways regulating the 2-cell-stage program. We identified retinoids as robust inducers of 2CLCs and the retinoic acid (RA)-signaling pathway as a key component of the regulatory circuitry of totipotent cells in embryos. Using single-cell RNA-seq, we reveal the transcriptional dynamics of 2CLC reprogramming and show that ES cells undergo distinct cellular trajectories in response to RA. Importantly, endogenous RA activity in early embryos is essential for zygotic genome activation and developmental progression. Overall, our data shed light on the gene regulatory networks controlling cellular plasticity and the totipotency program.

A cell in hand is worth two in the embryo: recent advances in 2-cell like cell reprogramming.

Multicellular organisms develop from a single cell, the zygote. This feature is referred to as totipotency. In the mouse, only the zygote and the 2-cell stage embryo display this attribute. Cells resembling the embryonic 2-cell stage blastomeres were identified in embryonic stem (ES) cell cultures as '2-cell-like cells' (2CLCs). This discovery brought the first cellular model with the possibility to investigate some features of the totipotent embryo and the molecular mechanisms regulating totipotency in vitro. In this article, we discuss the latest advancements on the research on 2CLCs, which have uncovered an intricate reprogramming process regulated by proteins as well as metabolites and ncRNAs. These recent findings have shed light on the combinatorial regulation of 2-cell-like cell emergence and the nature of their unique attributes.

A distinct metabolic state arises during the emergence of 2-cell-like cells.

Pluripotent stem cells are thought of as a surrogate of early developmental stages that sustain the capacity to generate all cell types in the body, thereby constituting an invaluable tool to address the mechanisms underlying cellular plasticity. In the mouse, cells resembling totipotent 2-cell-stage embryos (2-cell-like cells) arise at a very low frequency in embryonic stem cell (ESC) cultures. However, the extent to which these early-embryonic-like cells recapitulate the molecular features of the early embryo is unclear. Here, we have undertaken a characterization of some of the metabolic features of early-embryonic-like cells in culture. Our data indicate that early-embryonic-like cells exhibit decreased glycolytic and respiratory activity, lower levels of reactive oxygen species and increased glucose uptake, suggesting a shift of the metabolic programme during 2-cell-like cell reprogramming. Accordingly, we find that 2-cell-like cells can be induced by defined metabolites. Thus, in addition to their transcriptional and chromatin features, 2-cell-like cells recapitulate some of the metabolic features of their in vivo counterpart. Altogether, our work underscores a distinct metabolic state of early-embryonic-like cells and identifies compounds that can induce their emergence in vitro.

SUMO Safeguards Somatic and Pluripotent Cell Identities by Enforcing Distinct Chromatin States.

Understanding general principles that safeguard cellular identity should reveal critical insights into common mechanisms underlying specification of varied cell types. Here, we show that SUMO modification acts to stabilize cell fate in a variety of contexts. Hyposumoylation enhances pluripotency reprogramming in vitro and in vivo, increases lineage transdifferentiation, and facilitates leukemic cell differentiation. Suppressing sumoylation in embryonic stem cells (ESCs) promotes their conversion into 2-cell-embryo-like (2C-like) cells. During reprogramming to pluripotency, SUMO functions on fibroblastic enhancers to retain somatic transcription factors together with Oct4, Sox2, and Klf4, thus impeding somatic enhancer inactivation. In contrast, in ESCs, SUMO functions on heterochromatin to silence the 2C program, maintaining both proper H3K9me3 levels genome-wide and repression of the Dux locus by triggering recruitment of the sumoylated PRC1.6 and Kap/Setdb1 repressive complexes. Together, these studies show that SUMO acts on chromatin as a glue to stabilize key determinants of somatic and pluripotent states.

Lamin B1 mapping reveals the existence of dynamic and functional euchromatin lamin B1 domains.

Lamins (A/C and B) are major constituents of the nuclear lamina (NL). Structurally conserved lamina-associated domains (LADs) are formed by genomic regions that contact the NL. Lamins are also found in the nucleoplasm, with a yet unknown function. Here we map the genome-wide localization of lamin B1 in an euchromatin-enriched fraction of the mouse genome and follow its dynamics during the epithelial-to-mesenchymal transition (EMT). Lamin B1 associates with actively expressed and open euchromatin regions, forming dynamic euchromatin lamin B1-associated domains (eLADs) of about 0.3 Mb. Hi-C data link eLADs to the 3D organization of the mouse genome during EMT and correlate lamin B1 enrichment at topologically associating domain (TAD) borders with increased border strength. Having reduced levels of lamin B1 alters the EMT transcriptional signature and compromises the acquisition of mesenchymal traits. Thus, during EMT, the process of genome reorganization in mouse involves dynamic changes in eLADs.

Starting embryonic transcription for the first time.

Three studies highlight DUX proteins as key transcription factors regulating embryonic genome activation in early mammalian development.

Lysyl oxidase-like 2 (LOXL2) oxidizes trimethylated lysine 4 in histone H3.

Methylation of histone H3 lysine 4 is linked to active transcription and can be removed by LSD1 or the JmjC domain-containing proteins by amino-oxidation or hydroxylation, respectively. Here we describe that its deamination can be catalyzed by lysyl oxidase-like 2 protein (LOXL2), presenting an unconventional chemical mechanism for H3K4 modification. Infrared spectroscopy and mass spectrometry analyses demonstrated that recombinant LOXL2 specifically deaminates trimethylated H3K4. Moreover, by regulating H3K4me3 deamination, LOXL2 activity is linked with the transcriptional control of the CDH1 gene. These results reveal the existence of further H3 modification as well as a novel mechanism for H3K4me3 demethylation. DATABASE: The GEO accession number for the data referred to this paper is GSE35600.

Altered oncomodules underlie chromatin regulatory factors driver mutations.

Chromatin regulatory factors (CRFs), are known to be involved in tumorigenesis in several cancer types. Nevertheless, the molecular mechanisms through which driver alterations of CRFs cause tumorigenesis remain unknown. Here, we developed a CRFs Oncomodules Discovery approach, which mines several sources of cancer genomics and perturbaomics data. The approach prioritizes sets of genes significantly miss-regulated in primary tumors (oncomodules) bearing mutations of driver CRFs. We applied the approach to eleven TCGA tumor cohorts and uncovered oncomodules potentially associated to mutations of five driver CRFs in three cancer types. Our results revealed, for example, the potential involvement of the mTOR pathway in the development of tumors with loss-of-function mutations of MLL2 in head and neck squamous cell carcinomas. The experimental validation that MLL2 loss-of-function increases the sensitivity of cancer cell lines to mTOR inhibition lends further support to the validity of our approach. The potential oncogenic modules detected by our approach may guide experiments proposing ways to indirectly target driver mutations of CRFs.

LOXL2 Oxidizes Methylated TAF10 and Controls TFIID-Dependent Genes during Neural Progenitor Differentiation.

Protein function is often regulated and controlled by posttranslational modifications, such as oxidation. Although oxidation has been mainly considered to be uncontrolled and nonenzymatic, many enzymatic oxidations occur on enzyme-selected lysine residues; for instance, LOXL2 oxidizes lysines by converting the ε-amino groups into aldehyde groups. Using an unbiased proteomic approach, we have identified methylated TAF10, a member of the TFIID complex, as a LOXL2 substrate. LOXL2 oxidation of TAF10 induces its release from its promoters, leading to a block in TFIID-dependent gene transcription. In embryonic stem cells, this results in the inactivation of the pluripotency genes and loss of the pluripotent capacity. During zebrafish development, the absence of LOXL2 resulted in the aberrant overexpression of the neural progenitor gene Sox2 and impaired neural differentiation. Thus, lysine oxidation of the transcription factor TAF10 is a controlled protein modification and demonstrates a role for protein oxidation in regulating pluripotency genes.

A new role for LOX and LOXL2 proteins in transcription regulation.

The lysyl oxidase (LOX) family of proteins (LOX and LOXL1-LOXL4) oxidize amino groups located in the ε-position in lysines to generate an aldehyde group. In general, they are considered as extracellular proteins and have elastin and collagen as their main substrates. However, recent findings suggest a critical intracellular role for LOX and LOXL2 in transcriptional regulation. In this review, we highlight what is known about the transcriptional role of these two members of the family. Intriguingly, both the intracellular localization of these proteins and the fact that histones have been revealed to be their substrates place this family of proteins within the epigenetic field.

Regulation of heterochromatin transcription by Snail1/LOXL2 during epithelial-to-mesenchymal transition.

Although heterochromatin is enriched with repressive traits, it is also actively transcribed, giving rise to large amounts of noncoding RNAs. Although these RNAs are responsible for the formation and maintenance of heterochromatin, little is known about how their transcription is regulated. Here, we show that the Snail1 transcription factor represses mouse pericentromeric transcription, acting through the H3K4 deaminase LOXL2. Since Snail1 plays a key role in the epithelial-to-mesenchymal transition (EMT), we analyzed the regulation of heterochromatin transcription in this process. At the onset of EMT, one of the major structural heterochromatin proteins, HP1α, is transiently released from heterochromatin foci in a Snail1/LOXL2-dependent manner, concomitantly with a downregulation of major satellite transcription. Moreover, preventing the downregulation of major satellite transcripts compromised the migratory and invasive behavior of mesenchymal cells. We propose that Snail1 regulates heterochromatin transcription through LOXL2, thus creating the favorable transcriptional state necessary for completing EMT.