Single-cell Analysis Reveals Inter- and Intratumour Heterogeneity in Metastatic Breast Cancer.

Metastasis is the leading cause of cancer-related deaths of breast cancer patients. Some cancer cells in a tumour go through successive steps, referred to as the metastatic cascade, and give rise to metastases at a distant site. We know that the plasticity and heterogeneity of cancer cells play critical roles in metastasis but the precise underlying molecular mechanisms remain elusive. Here we aimed to identify molecular mechanisms of metastasis during colonization, one of the most important yet poorly understood steps of the cascade. We performed single-cell RNA-Seq (scRNA-Seq) on tumours and matched lung macrometastases of patient-derived xenografts of breast cancer. After correcting for confounding factors such as the cell cycle and the percentage of detected genes (PDG), we identified cells in three states in both tumours and metastases. Gene-set enrichment analysis revealed biological processes specific to proliferation and invasion in two states. Our findings suggest that these states are a balance between epithelial-to-mesenchymal (EMT) and mesenchymal-to-epithelial transitions (MET) traits that results in so-called partial EMT phenotypes. Analysis of the top differentially expressed genes (DEGs) between these cell states revealed a common set of partial EMT transcription factors (TFs) controlling gene expression, including ZNF750, OVOL2, TP63, TFAP2C and HEY2. Our data suggest that the TFs related to EMT delineate different cell states in tumours and metastases. The results highlight the marked interpatient heterogeneity of breast cancer but identify common features of single cells from five models of metastatic breast cancer.

Epigenetics, cryptorchidism, and infertility.

BACKGROUND: Cryptorchid boys with defective mini-puberty and impaired differentiation of Ad spermatogonia (high infertility risk) have altered expression of several genes encoding histone methyltransferases compared to patients with intact differentiation of gonocytes into Ad spermatogonia (low infertility risk). RESULTS: High infertility risk cryptorchid boys display hypogonadotropic hypogonadism, which, together with the diminished expression of histone deacetylases and increased expression of HDAC8 decrotonylase, indicates altered histone marks and, thus, a perturbed histone code. Curative GnRHa treatment induces normalization of histone methyltransferase, chromatin remodeling, and histone deacetylase gene expression. As a result, histone changes induce differentiation of Ad spermatogonia from their precursors and, thus, fertility. In this short report, we describe key functions of histone lysine methyltransferases, chromatin remodeling proteins, and long-noncoding RNAs, and discuss their potential roles in processes leading to infertility. CONCLUSION: Our findings suggest that epigenetic mechanisms are critical to better understanding the root causes underlying male infertility related to cryptorchidism and its possible transgenerational transmission.

RéSUMé: CONTEXTE: Chez les garçons cryptorchides qui présentent une minipuberté défectueuse et une différenciation altérée des spermatogonies Ad (risque élevé d'infertilité), l'expression de plusieurs gènes codant pour les histone méthyltransférases est altérée par rapport aux garçons présentant une différenciation intacte des gonocytes en spermatogonie Ad (faible risque d'infertilité). RéSULTATS: Les garçons cryptorchides à risque élevé d'infertilité présentent un hypogonadisme hypogonadotrope, qui, avec la diminution de l'expression des histone désacétylases et l'augmentation de l'expression de la décrotonyase HDAC8, indiquent une altération des marques d'histones et, par conséquent, un code d'histones perturbé. Le traitement curatif par la GnRHa induit une normalisation de l'histone-méthyltransférase, du remodelage de la chromatine et l'expression du gène de l'histone-désacétylase. En conséquence, les changements d'histones induisent la différenciation des spermatogonies Ad à partir de leurs précurseurs, et donc la fertilité. Dans le court rapport qui suit, nous décrivons les fonctions clés des histones lysine méthyltransférases, des protéines de remodelage de la chromatine et des ARN longs non codants; leur rôles potentiels dans les processus menant à l'infertilité sont discutés. CONCLUSION: Nos résultats suggèrent que les mécanismes épigénétiques sont un élément critique pour une meilleure compréhension des causes sous-jacentes de l'infertilité masculine liée à la cryptorchidie et sa transmission transgénérationnelle.

Convergence of multiple RNA-silencing pathways on GW182/TNRC6.

The RNA-binding protein TRIM71/LIN-41 is a phylogenetically conserved developmental regulator that functions in mammalian stem cell reprogramming, brain development, and cancer. TRIM71 recognizes target mRNAs through hairpin motifs and silences them through molecular mechanisms that await identification. Here, we uncover that TRIM71 represses its targets through RNA-supported interaction with TNRC6/GW182, a core component of the miRNA-induced silencing complex (miRISC). We demonstrate that AGO2, TRIM71, and UPF1 each recruit TNRC6 to specific sets of transcripts to silence them. As cellular TNRC6 levels are limiting, competition occurs among the silencing pathways, such that the loss of AGO proteins or of AGO binding to TNRC6 enhances the activities of the other pathways. We conclude that a miRNA-like silencing activity is shared among different mRNA silencing pathways and that the use of TNRC6 as a central hub provides a means to integrate their activities.

mutscan-a flexible R package for efficient end-to-end analysis of multiplexed assays of variant effect data.

Multiplexed assays of variant effect (MAVE) experimentally measure the effect of large numbers of sequence variants by selective enrichment of sequences with desirable properties followed by quantification by sequencing. mutscan is an R package for flexible analysis of such experiments, covering the entire workflow from raw reads up to statistical analysis and visualization. The core components are implemented in C++ for efficiency. Various experimental designs are supported, including single or paired reads with optional unique molecular identifiers. To find variants with changed relative abundance, mutscan employs established statistical models provided in the edgeR and limma packages. mutscan is available from https://github.com/fmicompbio/mutscan .

A multiple super-enhancer region establishes inter-TAD interactions and controls Hoxa function in cranial neural crest.

Enhancer-promoter interactions preferentially occur within boundary-insulated topologically associating domains (TADs), limiting inter-TAD interactions. Enhancer clusters in linear proximity, termed super-enhancers (SEs), ensure high target gene expression levels. Little is known about SE topological regulatory impact during craniofacial development. Here, we identify 2232 genome-wide putative SEs in mouse cranial neural crest cells (CNCCs), 147 of which target genes establishing CNCC positional identity during face formation. In second pharyngeal arch (PA2) CNCCs, a multiple SE-containing region, partitioned into Hoxa Inter-TAD Regulatory Element 1 and 2 (HIRE1 and HIRE2), establishes long-range inter-TAD interactions selectively with Hoxa2, that is required for external and middle ear structures. HIRE2 deletion in a Hoxa2 haploinsufficient background results in microtia. HIRE1 deletion phenocopies the full homeotic Hoxa2 knockout phenotype and induces PA3 and PA4 CNCC abnormalities correlating with Hoxa2 and Hoxa3 transcriptional downregulation. Thus, SEs can overcome TAD insulation and regulate anterior Hoxa gene collinear expression in a CNCC subpopulation-specific manner during craniofacial development.

Multimodal characterization of murine gastruloid development.

Gastruloids are 3D structures generated from pluripotent stem cells recapitulating fundamental principles of embryonic pattern formation. Using single-cell genomic analysis, we provide a resource mapping cell states and types during gastruloid development and compare them with the in vivo embryo. We developed a high-throughput handling and imaging pipeline to spatially monitor symmetry breaking during gastruloid development and report an early spatial variability in pluripotency determining a binary response to Wnt activation. Although cells in the gastruloid-core revert to pluripotency, peripheral cells become primitive streak-like. These two populations subsequently break radial symmetry and initiate axial elongation. By performing a compound screen, perturbing thousands of gastruloids, we derive a phenotypic landscape and infer networks of genetic interactions. Finally, using a dual Wnt modulation, we improve the formation of anterior structures in the existing gastruloid model. This work provides a resource to understand how gastruloids develop and generate complex patterns in vitro.

N-acetylcysteine overcomes NF1 loss-driven resistance to PI3Kα inhibition in breast cancer.

A genome-wide PiggyBac transposon-mediated screen and a resistance screen in a PIK3CAH1047R-mutated murine tumor model reveal NF1 loss in mammary tumors resistant to the phosphatidylinositol 3-kinase α (PI3Kα)-selective inhibitor alpelisib. Depletion of NF1 in PIK3CAH1047R breast cancer cell lines and a patient-derived organoid model shows that NF1 loss reduces sensitivity to PI3Kα inhibition and correlates with enhanced glycolysis and lower levels of reactive oxygen species (ROS). Unexpectedly, the antioxidant N-acetylcysteine (NAC) sensitizes NF1 knockout cells to PI3Kα inhibition and reverts their glycolytic phenotype. Global phospho-proteomics indicates that combination with NAC enhances the inhibitory effect of alpelisib on mTOR signaling. In public datasets of human breast cancer, we find that NF1 is frequently mutated and that such mutations are enriched in metastases, an indication for which use of PI3Kα inhibitors has been approved. Our results raise the attractive possibility of combining PI3Kα inhibition with NAC supplementation, especially in patients with drug-resistant metastases associated with NF1 loss.

Evidence that direct inhibition of transcription factor binding is the prevailing mode of gene and repeat repression by DNA methylation.

Cytosine methylation efficiently silences CpG-rich regulatory regions of genes and repeats in mammalian genomes. To what extent this entails direct inhibition of transcription factor (TF) binding versus indirect inhibition via recruitment of methyl-CpG-binding domain (MBD) proteins is unclear. Here we show that combinatorial genetic deletions of all four proteins with functional MBDs in mouse embryonic stem cells, derived neurons or a human cell line do not reactivate genes or repeats with methylated promoters. These do, however, become activated by methylation-restricted TFs if DNA methylation is removed. We identify several causal TFs in neurons, including ONECUT1, which is methylation sensitive only at a motif variant. Rampantly upregulated retrotransposons in methylation-free neurons feature a CRE motif, which activates them in the absence of DNA methylation via methylation-sensitive binding of CREB1. Our study reveals methylation-sensitive TFs in vivo and argues that direct inhibition, rather than indirect repression by the tested MBD proteins, is the prevailing mechanism of methylation-mediated repression at regulatory regions and repeats.

monaLisa: an R/Bioconductor package for identifying regulatory motifs.

SUMMARY: Proteins binding to specific nucleotide sequences, such as transcription factors, play key roles in the regulation of gene expression. Their binding can be indirectly observed via associated changes in transcription, chromatin accessibility, DNA methylation and histone modifications. Identifying candidate factors that are responsible for these observed experimental changes is critical to understand the underlying biological processes. Here, we present monaLisa, an R/Bioconductor package that implements approaches to identify relevant transcription factors from experimental data. The package can be easily integrated with other Bioconductor packages and enables seamless motif analyses without any software dependencies outside of R. AVAILABILITY AND IMPLEMENTATION: monaLisa is implemented in R and available on Bioconductor at https://bioconductor.org/packages/monaLisa with the development version hosted on GitHub at https://github.com/fmicompbio/monaLisa. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cooperation between HDAC3 and DAX1 mediates lineage restriction of embryonic stem cells.

Mouse embryonic stem cells (mESCs) are biased toward producing embryonic rather than extraembryonic endoderm fates. Here, we identify the mechanism of this barrier and report that the histone deacetylase Hdac3 and the transcriptional corepressor Dax1 cooperatively limit the lineage repertoire of mESCs by silencing an enhancer of the extraembryonic endoderm-specifying transcription factor Gata6. This restriction is opposed by the pluripotency transcription factors Nr5a2 and Esrrb, which promote cell type conversion. Perturbation of the barrier extends mESC potency and allows formation of 3D spheroids that mimic the spatial segregation of embryonic epiblast and extraembryonic endoderm in early embryos. Overall, this study shows that transcriptional repressors stabilize pluripotency by biasing the equilibrium between embryonic and extraembryonic lineages that is hardwired into the mESC transcriptional network.

The NFIB-ERO1A axis promotes breast cancer metastatic colonization of disseminated tumour cells.

Metastasis is the main cause of deaths related to solid cancers. Active transcriptional programmes are known to regulate the metastatic cascade but the molecular determinants of metastatic colonization remain elusive. Using an inducible piggyBac (PB) transposon mutagenesis screen, we have shown that overexpression of the transcription factor nuclear factor IB (NFIB) alone is sufficient to enhance primary mammary tumour growth and lung metastatic colonization. Mechanistically and functionally, NFIB directly increases expression of the oxidoreductase ERO1A, which enhances HIF1α-VEGFA-mediated angiogenesis and colonization, the last and fatal step of the metastatic cascade. NFIB is thus clinically relevant: it is preferentially expressed in the poor-prognostic group of basal-like breast cancers, and high expression of the NFIB/ERO1A/VEGFA pathway correlates with reduced breast cancer patient survival.

Molecular clues in the regulation of mini-puberty involve neuronal DNA binding transcription factor NHLH2.

Gonadotropin releasing hormone agonist (GnRHa) treatment following surgery to correct cryptorchidism restores mini-puberty via endocrinological and transcriptional effects and prevents adult infertility in most cases. Several genes are important for central hypogonadotropic hypogonadism in mammals, including many that are transcribed in both the brain and testis. However, the expression of these genes in prepubertal gonads has not been studied systematically, and little is known about the effect of hormone therapy on their testicular and neuronal expression levels. In this review, we interpret histological sections, data on hormone levels, and RNA profiling data from adult normal testes compared to pre-pubertal low infertility risk (LIR) and high infertility risk (HIR) patients randomly treated with surgery in combination with GnRHa or only surgery. We organize 31 target genes relevant for idiopathic hypogonadotropic hypogonadism and cryptorchidism into five classes depending on their expression levels in HIR versus LIR samples and their response to GnRHa treatment. Nescient-helix-loop-helix 2 (NHLH2) was the only gene showing a decreased mRNA level in HIR patients and an increase after GnRHa treatment. This phenomenon may reflect a broader effect of hormone treatment on gene expression in both testicular and central nervous system tissues, which could explain why the hypothalamus-pituitary-testicular axis is permanently restored by the administration of GnRHa.

RéSUMé: Le traitement par l'agoniste de l'hormone de libération des gonadotrophines (GnRHa) suite à une intervention chirurgicale pour cryptorchidie rétablit la mini-puberté par des effets endocrinologiques et transcriptionnels et prévient l'infertilité adulte dans la plupart des cas. Plusieurs gènes jouent un rôle important dans l'hypogonadisme hypogonadotrope central chez les mammifères, dont certains sont transcrits à la fois dans le cerveau et les testicules. Cependant, l'expression de ces gènes dans les gonades prépubères n'a pas été étudiée systématiquement et l'effet de l'hormonothérapie sur leurs niveaux d'expression testiculaire et neuronale n'est pas connu. Dans cette revue, nous interprétons des coupes histologiques, des données sur les taux d'hormones et des données de profilage d'ARN provenant de testicules normaux adultes et des patients prépubères à faible risque d'infertilité (LIR) et à haut risque d'infertilité (HIR) traités par chirurgie en association avec la GnRHa ou seulement la chirurgie dans le cadre d'une étude randomisée. Nous organisons 31 gènes cibles pertinents pour l'hypogonadisme hypogonadotrope idiopathique et la cryptorchidie en cinq classes en fonction de leurs niveaux d'expression dans les échantillons HIR et LIR et de leur réponse au traitement par GnRHa. Nescient-helix-loop-helix 2 (NHLH2) était l'unique gène dont le niveau d'ARNm diminue chez les patients HIR par rapport aux LIR et augmente suite au traitement par GnRHa. Ce phénomène pourrait être révélateur d'un effet généralisé du traitement hormonal sur l'expression des gènes dans les tissus testiculaires et du système nerveux central. Cela pourrait expliquer pourquoi l'axe hypothalamo-hypophyso-gonadique est définitivement rétablie par l'administration de la GnRHa.

A unique bipartite Polycomb signature regulates stimulus-response transcription during development.

Rapid cellular responses to environmental stimuli are fundamental for development and maturation. Immediate early genes can be transcriptionally induced within minutes in response to a variety of signals. How their induction levels are regulated and their untimely activation by spurious signals prevented during development is poorly understood. We found that in developing sensory neurons, before perinatal sensory-activity-dependent induction, immediate early genes are embedded into a unique bipartite Polycomb chromatin signature, carrying active H3K27ac on promoters but repressive Ezh2-dependent H3K27me3 on gene bodies. This bipartite signature is widely present in developing cell types, including embryonic stem cells. Polycomb marking of gene bodies inhibits mRNA elongation, dampening productive transcription, while still allowing for fast stimulus-dependent mark removal and bipartite gene induction. We reveal a developmental epigenetic mechanism regulating the rapidity and amplitude of the transcriptional response to relevant stimuli, while preventing inappropriate activation of stimulus-response genes.

Mammalian SWI/SNF continuously restores local accessibility to chromatin.

Chromatin accessibility is a hallmark of regulatory regions, entails transcription factor (TF) binding and requires nucleosomal reorganization. However, it remains unclear how dynamic this process is. In the present study, we use small-molecule inhibition of the catalytic subunit of the mouse SWI/SNF remodeler complex to show that accessibility and reduced nucleosome presence at TF-binding sites rely on persistent activity of nucleosome remodelers. Within minutes of remodeler inhibition, accessibility and TF binding decrease. Although this is irrespective of TF function, we show that the activating TF OCT4 (POU5F1) exhibits a faster response than the repressive TF REST. Accessibility, nucleosome depletion and gene expression are rapidly restored on inhibitor removal, suggesting that accessible chromatin is regenerated continuously and in a largely cell-autonomous fashion. We postulate that TF binding to chromatin and remodeler-mediated nucleosomal removal do not represent a stable situation, but instead accessible chromatin reflects an average of a dynamic process under continued renewal.

Preprocessing choices affect RNA velocity results for droplet scRNA-seq data.

Experimental single-cell approaches are becoming widely used for many purposes, including investigation of the dynamic behaviour of developing biological systems. Consequently, a large number of computational methods for extracting dynamic information from such data have been developed. One example is RNA velocity analysis, in which spliced and unspliced RNA abundances are jointly modeled in order to infer a 'direction of change' and thereby a future state for each cell in the gene expression space. Naturally, the accuracy and interpretability of the inferred RNA velocities depend crucially on the correctness of the estimated abundances. Here, we systematically compare five widely used quantification tools, in total yielding thirteen different quantification approaches, in terms of their estimates of spliced and unspliced RNA abundances in five experimental droplet scRNA-seq data sets. We show that there are substantial differences between the quantifications obtained from different tools, and identify typical genes for which such discrepancies are observed. We further show that these abundance differences propagate to the downstream analysis, and can have a large effect on estimated velocities as well as the biological interpretation. Our results highlight that abundance quantification is a crucial aspect of the RNA velocity analysis workflow, and that both the definition of the genomic features of interest and the quantification algorithm itself require careful consideration.

OBF1 and Oct factors control the germinal center transcriptional program.

OBF1 is a specific coactivator of the POU family transcription factors OCT1 and OCT2. OBF1 and OCT2 are B cell-specific and indispensable for germinal center (GC) formation, but their mechanism of action is unclear. Here, we show by chromatin immunoprecipitation-sequencing that OBF1 extensively colocalizes with OCT1 and OCT2. We found that these factors also often colocalize with transcription factors of the ETS family. Furthermore, we showed that OBF1, OCT2, and OCT1 bind widely to the promoters or enhancers of genes involved in GC formation in mouse and human GC B cells. Short hairpin RNA knockdown experiments demonstrated that OCT1, OCT2, and OBF1 regulate each other and are essential for proliferation of GC-derived lymphoma cell lines. OBF1 downregulation disrupts the GC transcriptional program: genes involved in GC maintenance, such as BCL6, are downregulated, whereas genes related to exit from the GC program, such as IRF4, are upregulated. Ectopic expression of BCL6 does not restore the proliferation of GC-derived lymphoma cells depleted of OBF1 unless IRF4 is also depleted, indicating that OBF1 controls an essential regulatory node in GC differentiation.

Phenotypic landscape of intestinal organoid regeneration.

The development of intestinal organoids from single adult intestinal stem cells in vitro recapitulates the regenerative capacity of the intestinal epithelium1,2. Here we unravel the mechanisms that orchestrate both organoid formation and the regeneration of intestinal tissue, using an image-based screen to assay an annotated library of compounds. We generate multivariate feature profiles for hundreds of thousands of organoids to quantitatively describe their phenotypic landscape. We then use these phenotypic fingerprints to infer regulatory genetic interactions, establishing a new approach to the mapping of genetic interactions in an emergent system. This allows us to identify genes that regulate cell-fate transitions and maintain the balance between regeneration and homeostasis, unravelling previously unknown roles for several pathways, among them retinoic acid signalling. We then characterize a crucial role for retinoic acid nuclear receptors in controlling exit from the regenerative state and driving enterocyte differentiation. By combining quantitative imaging with RNA sequencing, we show the role of endogenous retinoic acid metabolism in initiating transcriptional programs that guide the cell-fate transitions of intestinal epithelium, and we identify an inhibitor of the retinoid X receptor that improves intestinal regeneration in vivo.

ExploreModelMatrix: Interactive exploration for improved understanding of design matrices and linear models in R.

Linear and generalized linear models are used extensively in many scientific fields, to model observed data and as the basis for hypothesis tests. The use of such models requires specification of a design matrix, and subsequent formulation of contrasts representing scientific hypotheses of interest. Proper execution of these steps requires a thorough understanding of the meaning of the individual coefficients, and is a frequent source of uncertainty for end users. Here, we present an R/Bioconductor package, ExploreModelMatrix, which enables interactive exploration of design matrices and linear model diagnostics. Given a sample data table and a desired design formula, the package displays how the model coefficients are combined to give the fitted values for each combination of predictor variables, which allows users to both extract the interpretation of each individual coefficient, and formulate desired linear contrasts. In addition, the interactive interface displays informative characteristics for the regular linear model corresponding to the provided design, such as variance inflation factors and the pseudoinverse of the design matrix. We envision the package and the built-in collection of common types of linear model designs to be useful for teaching and self-learning purposes, as well as for assisting more experienced users in the interpretation of complex model designs.

Cell Types of the Human Retina and Its Organoids at Single-Cell Resolution.

Human organoids recapitulating the cell-type diversity and function of their target organ are valuable for basic and translational research. We developed light-sensitive human retinal organoids with multiple nuclear and synaptic layers and functional synapses. We sequenced the RNA of 285,441 single cells from these organoids at seven developmental time points and from the periphery, fovea, pigment epithelium and choroid of light-responsive adult human retinas, and performed histochemistry. Cell types in organoids matured in vitro to a stable "developed" state at a rate similar to human retina development in vivo. Transcriptomes of organoid cell types converged toward the transcriptomes of adult peripheral retinal cell types. Expression of disease-associated genes was cell-type-specific in adult retina, and cell-type specificity was retained in organoids. We implicate unexpected cell types in diseases such as macular degeneration. This resource identifies cellular targets for studying disease mechanisms in organoids and for targeted repair in human retinas.

Zfp281 orchestrates interconversion of pluripotent states by engaging Ehmt1 and Zic2.

Developmental cell fate specification is a unidirectional process that can be reverted in response to injury or experimental reprogramming. Whether differentiation and de-differentiation trajectories intersect mechanistically is unclear. Here, we performed comparative screening in lineage-related mouse naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), and identified the constitutively expressed zinc finger transcription factor (TF) Zfp281 as a bidirectional regulator of cell state interconversion. We showed that subtle chromatin binding changes in differentiated cells translate into activation of the histone H3 lysine 9 (H3K9) methyltransferase Ehmt1 and stabilization of the zinc finger TF Zic2 at enhancers and promoters. Genetic gain-of-function and loss-of-function experiments confirmed a critical role of Ehmt1 and Zic2 downstream of Zfp281 both in driving exit from the ESC state and in restricting reprogramming of EpiSCs. Our study reveals that cell type-invariant chromatin association of Zfp281 provides an interaction platform for remodeling the cis-regulatory network underlying cellular plasticity.