Co-transcriptional splicing facilitates transcription of gigantic genes.

Although introns are typically tens to thousands of nucleotides, there are notable exceptions. In flies as well as humans, a small number of genes contain introns that are more than 1000 times larger than typical introns, exceeding hundreds of kilobases (kb) to megabases (Mb). It remains unknown why gigantic introns exist and how cells overcome the challenges associated with their transcription and RNA processing. The Drosophila Y chromosome contains some of the largest genes identified to date: multiple genes exceed 4Mb, with introns accounting for over 99% of the gene span. Here we demonstrate that co-transcriptional splicing of these gigantic Y-linked genes is important to ensure successful transcription: perturbation of splicing led to the attenuation of transcription, leading to a failure to produce mature mRNA. Cytologically, defective splicing of the Y-linked gigantic genes resulted in disorganization of transcripts within the nucleus suggestive of entanglement of transcripts, likely resulting from unspliced long RNAs. We propose that co-transcriptional splicing maintains the length of nascent transcripts of gigantic genes under a critical threshold, preventing their entanglement and ensuring proper gene expression. Our study reveals a novel biological significance of co-transcriptional splicing.

The nuclear receptor THRB facilitates differentiation of human PSCs into more mature hepatocytes.

To understand the mechanisms regulating the in vitro maturation of hPSC-derived hepatocytes, we developed a 3D differentiation system and compared gene regulatory elements in human primary hepatocytes with those in hPSC-hepatocytes that were differentiated in 2D or 3D conditions by RNA-seq, ATAC-seq, and H3K27Ac ChIP-seq. Regulome comparisons showed a reduced enrichment of thyroid receptor THRB motifs in accessible chromatin and active enhancers without a reduced transcription of THRB. The addition of thyroid hormone T3 increased the binding of THRB to the CYP3A4 proximal enhancer, restored the super-enhancer status and gene expression of NFIC, and reduced the expression of AFP. The resultant hPSC-hepatocytes showed gene expression, epigenetic status, and super-enhancer landscape closer to primary hepatocytes and activated regulatory regions including non-coding SNPs associated with liver-related diseases. Transplanting the hPSC-hepatocytes resulted in the engraftment of human hepatocytes into the mouse liver without disrupting normal liver histology. This work implicates the environmental factor-nuclear receptor axis in regulating the maturation of hPSC-hepatocytes.

Genome-wide CRISPR screen identifies PRC2 and KMT2D-COMPASS as regulators of distinct EMT trajectories that contribute differentially to metastasis.

Epithelial-mesenchymal transition (EMT) programs operate within carcinoma cells, where they generate phenotypes associated with malignant progression. In their various manifestations, EMT programs enable epithelial cells to enter into a series of intermediate states arrayed along the E-M phenotypic spectrum. At present, we lack a coherent understanding of how carcinoma cells control their entrance into and continued residence in these various states, and which of these states favour the process of metastasis. Here we characterize a layer of EMT-regulating machinery that governs E-M plasticity (EMP). This machinery consists of two chromatin-modifying complexes, PRC2 and KMT2D-COMPASS, which operate as critical regulators to maintain a stable epithelial state. Interestingly, loss of these two complexes unlocks two distinct EMT trajectories. Dysfunction of PRC2, but not KMT2D-COMPASS, yields a quasi-mesenchymal state that is associated with highly metastatic capabilities and poor survival of patients with breast cancer, suggesting that great caution should be applied when PRC2 inhibitors are evaluated clinically in certain patient cohorts. These observations identify epigenetic factors that regulate EMP, determine specific intermediate EMT states and, as a direct consequence, govern the metastatic ability of carcinoma cells.

Direct and Indirect Regulators of Epithelial-Mesenchymal Transition-Mediated Immunosuppression in Breast Carcinomas.

The epithelial-to-mesenchymal transition, which conveys epithelial (E) carcinoma cells to quasi-mesenchymal (qM) states, enables them to metastasize and acquire resistance to certain treatments. Murine tumors composed of qM mammary carcinoma cells assemble an immunosuppressive tumor microenvironment (TME) and develop resistance to anti-CTLA4 immune-checkpoint blockade (ICB) therapy, unlike their E counterparts. Importantly, minority populations of qM cells within a tumor can cross-protect their more E neighbors from immune attack. The underlying mechanisms of immunosuppression and cross-protection have been unclear. We demonstrate that abrogation of qM carcinoma cell-derived factors (CD73, CSF1, or SPP1) prevents the assembly of an immunosuppressive TME and sensitizes otherwise refractory qM tumors partially or completely to anti-CTLA4 ICB. Most strikingly, mixed tumors in which minority populations of carcinoma cells no longer express CD73 are now sensitized to anti-CTLA4 ICB. Finally, loss of CD73 also enhances the efficacy of anti-CTLA4 ICB during the process of metastatic colonization. SIGNIFICANCE: Minority populations of qM carcinoma cells, which likely reside in human breast carcinomas, can cross-protect their E neighbors from immune attack. Understanding the mechanisms by which qM carcinoma cells resist antitumor immune attack can help identify signaling channels that can be interrupted to potentiate the efficacy of checkpoint blockade immunotherapies.This article is highlighted in the In This Issue feature, p. 995.

Author Correction: Mitochondrial metabolism promotes adaptation to proteotoxic stress.

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

Mitochondrial metabolism promotes adaptation to proteotoxic stress.

The mechanisms by which cells adapt to proteotoxic stress are largely unknown, but are key to understanding how tumor cells, particularly in vivo, are largely resistant to proteasome inhibitors. Analysis of cancer cell lines, mouse xenografts and patient-derived tumor samples all showed an association between mitochondrial metabolism and proteasome inhibitor sensitivity. When cells were forced to use oxidative phosphorylation rather than glycolysis, they became proteasome-inhibitor resistant. This mitochondrial state, however, creates a unique vulnerability: sensitivity to the small molecule compound elesclomol. Genome-wide CRISPR-Cas9 screening showed that a single gene, encoding the mitochondrial reductase FDX1, could rescue elesclomol-induced cell death. Enzymatic function and nuclear-magnetic-resonance-based analyses further showed that FDX1 is the direct target of elesclomol, which promotes a unique form of copper-dependent cell death. These studies explain a fundamental mechanism by which cells adapt to proteotoxic stress and suggest strategies to mitigate proteasome inhibitor resistance.

Acquisition of a hybrid E/M state is essential for tumorigenicity of basal breast cancer cells.

Carcinoma cells residing in an intermediate phenotypic state along the epithelial-mesenchymal (E-M) spectrum are associated with malignant phenotypes, such as invasiveness, tumor-initiating ability, and metastatic dissemination. Using the recently described CD104+/CD44hi antigen marker combination, we isolated highly tumorigenic breast cancer cells residing stably-both in vitro and in vivo-in an intermediate phenotypic state and coexpressing both epithelial (E) and mesenchymal (M) markers. We demonstrate that tumorigenicity depends on individual cells residing in this E/M hybrid state and cannot be phenocopied by mixing two cell populations that reside stably at the two ends of the spectrum, i.e., in the E and in the M state. Hence, residence in a specific intermediate state along the E-M spectrum rather than phenotypic plasticity appears critical to the expression of tumor-initiating capacity. Acquisition of this E/M hybrid state is facilitated by the differential expression of EMT-inducing transcription factors (EMT-TFs) and is accompanied by the expression of adult stem cell programs, notably, active canonical Wnt signaling. Furthermore, transition from the highly tumorigenic E/M state to a fully mesenchymal phenotype, achieved by constitutive ectopic expression of Zeb1, is sufficient to drive cells out of the E/M hybrid state into a highly mesenchymal state, which is accompanied by a substantial loss of tumorigenicity and a switch from canonical to noncanonical Wnt signaling. Identifying the gatekeepers of the various phenotypic states arrayed along the E-M spectrum is likely to prove useful in developing therapeutic approaches that operate by shifting cancer cells between distinct states along this spectrum.

Predicting microRNA targeting efficacy in Drosophila.

BACKGROUND: MicroRNAs (miRNAs) are short regulatory RNAs that derive from hairpin precursors. Important for understanding the functional roles of miRNAs is the ability to predict the messenger RNA (mRNA) targets most responsive to each miRNA. Progress towards developing quantitative models of miRNA targeting in Drosophila and other invertebrate species has lagged behind that of mammals due to the paucity of datasets measuring the effects of miRNAs on mRNA levels. RESULTS: We acquired datasets suitable for the quantitative study of miRNA targeting in Drosophila. Analyses of these data expanded the types of regulatory sites known to be effective in flies, expanded the mRNA regions with detectable targeting to include 5' untranslated regions, and identified features of site context that correlate with targeting efficacy in fly cells. Updated evolutionary analyses evaluated the probability of conserved targeting for each predicted site and indicated that more than a third of the Drosophila genes are preferentially conserved targets of miRNAs. Based on these results, a quantitative model was developed to predict targeting efficacy in insects. This model performed better than existing models, and it drives the most recent version, v7, of TargetScanFly. CONCLUSIONS: Our evolutionary and functional analyses expand the known scope of miRNA targeting in flies and other insects. The existence of a quantitative model that has been developed and trained using Drosophila data will provide a valuable resource for placing miRNAs into gene regulatory networks of this important experimental organism.

Genetic and epigenetic determinants establish a continuum of Hsf1 occupancy and activity across the yeast genome.

Heat shock factor 1 is the master transcriptional regulator of molecular chaperones and binds to the same cis-acting heat shock element (HSE) across the eukaryotic lineage. In budding yeast, Hsf1 drives the transcription of ∼20 genes essential to maintain proteostasis under basal conditions, yet its specific targets and extent of inducible binding during heat shock remain unclear. Here we combine Hsf1 chromatin immunoprecipitation sequencing (seq), nascent RNA-seq, and Hsf1 nuclear depletion to quantify Hsf1 binding and transcription across the yeast genome. We find that Hsf1 binds 74 loci during acute heat shock, and these are linked to 46 genes with strong Hsf1-dependent expression. Notably, Hsf1's induced DNA binding leads to a disproportionate (∼7.5-fold) increase in nascent transcription. Promoters with high basal Hsf1 occupancy have nucleosome-depleted regions due to the presence of "pioneer factors." These accessible sites are likely critical for Hsf1 occupancy as the activator is incapable of binding HSEs within a stably positioned, reconstituted nucleosome. In response to heat shock, however, Hsf1 accesses nucleosomal sites and promotes chromatin disassembly in concert with the Remodels Structure of Chromatin (RSC) complex. Our data suggest that the interplay between nucleosome positioning, HSE strength, and active Hsf1 levels allows cells to precisely tune expression of the proteostasis network.

SYK kinase mediates brown fat differentiation and activation.

Brown adipose tissue (BAT) metabolism influences glucose homeostasis and metabolic health in mice and humans. Sympathetic stimulation of ß-adrenergic receptors in response to cold induces proliferation, differentiation, and UCP1 expression in pre-adipocytes and mature brown adipocytes. Here we show that spleen tyrosine kinase (SYK) is upregulated during brown adipocyte differentiation and activated by ß-adrenergic stimulation. Deletion or inhibition of SYK, a kinase known for its essential roles in the immune system, blocks brown and white pre-adipocyte proliferation and differentiation in vitro, and results in diminished expression of Ucp1 and other genes regulating brown adipocyte function in response to ß-adrenergic stimulation. Adipocyte-specific SYK deletion in mice reduces BAT mass and BAT that developed consisted of SYK-expressing brown adipocytes that had escaped homozygous Syk deletion. SYK inhibition in vivo represses ß-agonist-induced thermogenesis and oxygen consumption. These results establish SYK as an essential mediator of brown fat formation and function.

Integrin-ß4 identifies cancer stem cell-enriched populations of partially mesenchymal carcinoma cells.

Neoplastic cells within individual carcinomas often exhibit considerable phenotypic heterogeneity in their epithelial versus mesenchymal-like cell states. Because carcinoma cells with mesenchymal features are often more resistant to therapy and may serve as a source of relapse, we sought to determine whether such cells could be further stratified into functionally distinct subtypes. Indeed, we find that a basal epithelial marker, integrin-ß4 (ITGB4), can be used to enable stratification of mesenchymal-like triple-negative breast cancer (TNBC) cells that differ from one another in their relative tumorigenic abilities. Notably, we demonstrate that ITGB4+ cancer stem cell (CSC)-enriched mesenchymal cells reside in an intermediate epithelial/mesenchymal phenotypic state. Among patients with TNBC who received chemotherapy, elevated ITGB4 expression was associated with a worse 5-year probability of relapse-free survival. Mechanistically, we find that the ZEB1 (zinc finger E-box binding homeobox 1) transcription factor activity in highly mesenchymal SUM159 TNBC cells can repress expression of the epithelial transcription factor TAp63α (tumor protein 63 isoform 1), a protein that promotes ITGB4 expression. In addition, we demonstrate that ZEB1 and ITGB4 are important in modulating the histopathological phenotypes of tumors derived from mesenchymal TNBC cells. Hence, mesenchymal carcinoma cell populations are internally heterogeneous, and ITGB4 is a mechanistically driven prognostic biomarker that can be used to identify the more aggressive subtypes of mesenchymal carcinoma cells in TNBC. The ability to rapidly isolate and mechanistically interrogate the CSC-enriched, partially mesenchymal carcinoma cells should further enable identification of novel therapeutic opportunities to improve the prognosis for high-risk patients with TNBC.

LACTB is a tumour suppressor that modulates lipid metabolism and cell state.

Post-mitotic, differentiated cells exhibit a variety of characteristics that contrast with those of actively growing neoplastic cells, such as the expression of cell-cycle inhibitors and differentiation factors. We hypothesized that the gene expression profiles of these differentiated cells could reveal the identities of genes that may function as tumour suppressors. Here we show, using in vitro and in vivo studies in mice and humans, that the mitochondrial protein LACTB potently inhibits the proliferation of breast cancer cells. Its mechanism of action involves alteration of mitochondrial lipid metabolism and differentiation of breast cancer cells. This is achieved, at least in part, through reduction of the levels of mitochondrial phosphatidylserine decarboxylase, which is involved in the synthesis of mitochondrial phosphatidylethanolamine. These observations uncover a novel mitochondrial tumour suppressor and demonstrate a connection between mitochondrial lipid metabolism and the differentiation program of breast cancer cells, thereby revealing a previously undescribed mechanism of tumour suppression.

Suppression of 19S proteasome subunits marks emergence of an altered cell state in diverse cancers.

The use of proteasome inhibitors to target cancer's dependence on altered protein homeostasis has been greatly limited by intrinsic and acquired resistance. Analyzing data from thousands of cancer lines and tumors, we find that those with suppressed expression of one or more 19S proteasome subunits show intrinsic proteasome inhibitor resistance. Moreover, such proteasome subunit suppression is associated with poor outcome in myeloma patients, where proteasome inhibitors are a mainstay of treatment. Beyond conferring resistance to proteasome inhibitors, proteasome subunit suppression also serves as a sentinel of a more global remodeling of the transcriptome. This remodeling produces a distinct gene signature and new vulnerabilities to the proapoptotic drug, ABT-263. This frequent, naturally arising imbalance in 19S regulatory complex composition is achieved through a variety of mechanisms, including DNA methylation, and marks the emergence of a heritably altered and therapeutically relevant state in diverse cancers.

A Genome-wide CRISPR Screen in Toxoplasma Identifies Essential Apicomplexan Genes.

Apicomplexan parasites are leading causes of human and livestock diseases such as malaria and toxoplasmosis, yet most of their genes remain uncharacterized. Here, we present the first genome-wide genetic screen of an apicomplexan. We adapted CRISPR/Cas9 to assess the contribution of each gene from the parasite Toxoplasma gondii during infection of human fibroblasts. Our analysis defines ∼200 previously uncharacterized, fitness-conferring genes unique to the phylum, from which 16 were investigated, revealing essential functions during infection of human cells. Secondary screens identify as an invasion factor the claudin-like apicomplexan microneme protein (CLAMP), which resembles mammalian tight-junction proteins and localizes to secretory organelles, making it critical to the initiation of infection. CLAMP is present throughout sequenced apicomplexan genomes and is essential during the asexual stages of the malaria parasite Plasmodium falciparum. These results provide broad-based functional information on T. gondii genes and will facilitate future approaches to expand the horizon of antiparasitic interventions.

Activation of PKA leads to mesenchymal-to-epithelial transition and loss of tumor-initiating ability.

The epithelial-to-mesenchymal transition enables carcinoma cells to acquire malignancy-associated traits and the properties of tumor-initiating cells (TICs). TICs have emerged in recent years as important targets for cancer therapy, owing to their ability to drive clinical relapse and enable metastasis. Here, we propose a strategy to eliminate mesenchymal TICs by inducing their conversion to more epithelial counterparts that have lost tumor-initiating ability. We report that increases in intracellular levels of the second messenger, adenosine 3',5'-monophosphate, and the subsequent activation of protein kinase A (PKA) induce a mesenchymal-to-epithelial transition (MET) in mesenchymal human mammary epithelial cells. PKA activation triggers epigenetic reprogramming of TICs by the histone demethylase PHF2, which promotes their differentiation and loss of tumor-initiating ability. This study provides proof-of-principle for inducing an MET as differentiation therapy for TICs and uncovers a role for PKA in enforcing and maintaining the epithelial state.

MERAV: a tool for comparing gene expression across human tissues and cell types.

The oncogenic transformation of normal cells into malignant, rapidly proliferating cells requires major alterations in cell physiology. For example, the transformed cells remodel their metabolic processes to supply the additional demand for cellular building blocks. We have recently demonstrated essential metabolic processes in tumor progression through the development of a methodological analysis of gene expression. Here, we present the Metabolic gEne RApid Visualizer (MERAV, http://merav.wi.mit.edu), a web-based tool that can query a database comprising ∼4300 microarrays, representing human gene expression in normal tissues, cancer cell lines and primary tumors. MERAV has been designed as a powerful tool for whole genome analysis which offers multiple advantages: one can search many genes in parallel; compare gene expression among different tissue types as well as between normal and cancer cells; download raw data; and generate heatmaps; and finally, use its internal statistical tool. Most importantly, MERAV has been designed as a unique tool for analyzing metabolic processes as it includes matrixes specifically focused on metabolic genes and is linked to the Kyoto Encyclopedia of Genes and Genomes pathway search.

A respiratory chain controlled signal transduction cascade in the mitochondrial intermembrane space mediates hydrogen peroxide signaling.

Reactive oxygen species (ROS) such as hydrogen peroxide (H2O2) govern cellular homeostasis by inducing signaling. H2O2 modulates the activity of phosphatases and many other signaling molecules through oxidation of critical cysteine residues, which led to the notion that initiation of ROS signaling is broad and nonspecific, and thus fundamentally distinct from other signaling pathways. Here, we report that H2O2 signaling bears hallmarks of a regular signal transduction cascade. It is controlled by hierarchical signaling events resulting in a focused response as the results place the mitochondrial respiratory chain upstream of tyrosine-protein kinase Lyn, Lyn upstream of tyrosine-protein kinase SYK (Syk), and Syk upstream of numerous targets involved in signaling, transcription, translation, metabolism, and cell cycle regulation. The active mediators of H2O2 signaling colocalize as H2O2 induces mitochondria-associated Lyn and Syk phosphorylation, and a pool of Lyn and Syk reside in the mitochondrial intermembrane space. Finally, the same intermediaries control the signaling response in tissues and species responsive to H2O2 as the respiratory chain, Lyn, and Syk were similarly required for H2O2 signaling in mouse B cells, fibroblasts, and chicken DT40 B cells. Consistent with a broad role, the Syk pathway is coexpressed across tissues, is of early metazoan origin, and displays evidence of evolutionary constraint in the human. These results suggest that H2O2 signaling is under control of a signal transduction pathway that links the respiratory chain to the mitochondrial intermembrane space-localized, ubiquitous, and ancient Syk pathway in hematopoietic and nonhematopoietic cells.

Dihydropyrimidine accumulation is required for the epithelial-mesenchymal transition.

It is increasingly appreciated that oncogenic transformation alters cellular metabolism to facilitate cell proliferation, but less is known about the metabolic changes that promote cancer cell aggressiveness. Here, we analyzed metabolic gene expression in cancer cell lines and found that a set of high-grade carcinoma lines expressing mesenchymal markers share a unique 44 gene signature, designated the "mesenchymal metabolic signature" (MMS). A FACS-based shRNA screen identified several MMS genes as essential for the epithelial-mesenchymal transition (EMT), but not for cell proliferation. Dihydropyrimidine dehydrogenase (DPYD), a pyrimidine-degrading enzyme, was highly expressed upon EMT induction and was necessary for cells to acquire mesenchymal characteristics in vitro and for tumorigenic cells to extravasate into the mouse lung. This role of DPYD was mediated through its catalytic activity and enzymatic products, the dihydropyrimidines. Thus, we identify metabolic processes essential for the EMT, a program associated with the acquisition of metastatic and aggressive cancer cell traits.

Altered translation of GATA1 in Diamond-Blackfan anemia.

Ribosomal protein haploinsufficiency occurs in diverse human diseases including Diamond-Blackfan anemia (DBA), congenital asplenia and T cell leukemia. Yet, how mutations in genes encoding ubiquitously expressed proteins such as these result in cell-type- and tissue-specific defects remains unknown. Here, we identify mutations in GATA1, encoding the critical hematopoietic transcription factor GATA-binding protein-1, that reduce levels of full-length GATA1 protein and cause DBA in rare instances. We show that ribosomal protein haploinsufficiency, the more common cause of DBA, can lead to decreased GATA1 mRNA translation, possibly resulting from a higher threshold for initiation of translation of this mRNA in comparison with other mRNAs. In primary hematopoietic cells from patients with mutations in RPS19, encoding ribosomal protein S19, the amplitude of a transcriptional signature of GATA1 target genes was globally and specifically reduced, indicating that the activity, but not the mRNA level, of GATA1 is decreased in patients with DBA associated with mutations affecting ribosomal proteins. Moreover, the defective hematopoiesis observed in patients with DBA associated with ribosomal protein haploinsufficiency could be partially overcome by increasing GATA1 protein levels. Our results provide a paradigm by which selective defects in translation due to mutations affecting ubiquitous ribosomal proteins can result in human disease.