Single-cell m6A mapping in vivo using picoMeRIP-seq.

Current N6-methyladenosine (m6A) mapping methods need large amounts of RNA or are limited to cultured cells. Through optimized sample recovery and signal-to-noise ratio, we developed picogram-scale m6A RNA immunoprecipitation and sequencing (picoMeRIP-seq) for studying m6A in vivo in single cells and scarce cell types using standard laboratory equipment. We benchmark m6A mapping on titrations of poly(A) RNA and embryonic stem cells and in single zebrafish zygotes, mouse oocytes and embryos.

The RNA m6A landscape of mouse oocytes and preimplantation embryos.

Despite the significance of N6-methyladenosine (m6A) in gene regulation, the requirement for large amounts of RNA has hindered m6A profiling in mammalian early embryos. Here we apply low-input methyl RNA immunoprecipitation and sequencing to map m6A in mouse oocytes and preimplantation embryos. We define the landscape of m6A during the maternal-to-zygotic transition, including stage-specifically expressed transcription factors essential for cell fate determination. Both the maternally inherited transcripts to be degraded post fertilization and the zygotically activated genes during zygotic genome activation are widely marked by m6A. In contrast to m6A-marked zygotic ally-activated genes, m6A-marked maternally inherited transcripts have a higher tendency to be targeted by microRNAs. Moreover, RNAs derived from retrotransposons, such as MTA that is maternally expressed and MERVL that is transcriptionally activated at the two-cell stage, are largely marked by m6A. Our results provide a foundation for future studies exploring the regulatory roles of m6A in mammalian early embryonic development.

RNA m6A modifications in mammalian gametogenesis and pregnancy.

In brief: RNA modifications play key roles in regulating various biological processes. This article discusses and summarizes the recent advances of RNA m6A modifications related to mammalian gametogenesis, early embryonic development, and miscarriage. Abstract: The epitranscriptome is defined as the collection of post-transcriptional chemical modifications of RNA in a cell. RNA methylation refers to the chemical post-transcriptional modification of RNA by selectively adding methyl groups under the catalysis of a methyltransferase. The N6 methyladenosine (m6A) is one of the most common of the more than 100 known RNA modifications. Recent research has revealed that RNA m6A modifications are reversible. Additionally, m6A containing RNA can be selectively identified by immunoprecipitation followed by high-throughput sequencing (MeRIP-SEQ). These two developments have inspired a tremendous effort to unravel the biological role of m6A. The role of RNA m6A modifications in immune regulation, cell division, stem cell renewal, gametogenesis, embryonic development, and placental function has gradually emerged, which is of great significance for the study of post-transcriptional regulation of gene expression in reproductive biology. This review summarizes the current knowledge about RNA m6A modification in a variety of mammalian reproductive events.

Depletion of the m1A writer TRMT6/TRMT61A reduces proliferation and resistance against cellular stress in bladder cancer.

Background: Bladder cancer (BLCA) is a common and deadly disease that results in a reduced quality of life for the patients and a significant economic burden on society. A better understanding of tumorigenesis is needed to improve clinical outcomes. Recent evidence places the RNA modification m1A and its regulatory proteins TRMT6/TRMT61A and ALKBH3 in BLCA pathogenesis. Methods: TRMT6/TRMT61A, ALKBH1, and ALKBH3 expression was examined in human BLCA cell lines and a normal urinary tract epithelium cell line through qRT-PCR and western blot analysis. Prestoblue Cell Viability Reagent, wound-healing assay, and live-cell imaging-based cell displacement analysis, were conducted to assess proliferation, migration, and displacement of this BLCA cell line panel. Cell survival was assessed after inducing cellular stress and activating the unfolded protein response (UPR) with tunicamycin. Moreover, siRNA-mediated gene silencing in two BLCA cell lines (5637 and HT1197) was conducted to investigate the biological roles of TRMT6/TRMT61A. Results: Heterogeneous morphology, proliferation, displacement, tunicamycin sensitivity, and expression levels of m1A regulators were observed among the panel of cell lines examined. In general, TRMT61A expression was increased in BLCA cell lines when compared to SV-HUC-1. Depletion of TRMT6/TRMT61A reduced proliferation capacity in both 5637 and HT1197 cell lines. The average cell displacement of 5637 was also reduced upon TRMT6/TRMT61A depletion. Interestingly, TRMT6/TRMT61A depletion decreased mRNA expression of targets associated with the ATF6-branch of the UPR in 5637 but not in HT1197. Moreover, cell survival after induction of cellular stress was compromised after TRMT6/TRMT61A knockdown in 5637 but not in HT1197 cells. Conclusion: The findings suggest that TRMT6/TRMT61A plays an oncogenic role in BLCA and is involved in desensitizing BLCA cells against cellular stress. Further investigation into the regulation of TRMT6/TRMT61A expression and its impact on cellular stress tolerance may provide insights for future BLCA treatment.

Characterization of novel small non-coding RNAs and their modifications in bladder cancer using an updated small RNA-seq workflow.

Background: Bladder cancer (BLCA) is one of the most common cancer types worldwide. The disease is responsible for about 200,000 deaths annually, thus improved diagnostics and therapy is needed. A large body of evidence reveal that small RNAs of less than 40 nucleotides may act as tumor suppressors, oncogenes, and disease biomarkers, with a major focus on microRNAs. However, the role of other families of small RNAs is not yet deciphered. Recent results suggest that small RNAs and their modification status, play a role in BLCA development and are promising biomarkers due to their high abundance in the exomes and body fluids (including urine). Moreover, free modified nucleosides have been detected at elevated levels from the urine of BLCA patients. A genome-wide view of small RNAs, and their modifications, will help pinpoint the molecules that could be used as biomarker or has important biology in BLCA development. Methods: BLCA tumor tissue specimens were obtained from 12 patients undergoing transurethral resection of non-muscle invasive papillary urothelial carcinomas. Genome-wide profiling of small RNAs less than 40 bases long was performed by a modified protocol with TGIRT (thermostable group II reverse transcriptase) to identify novel small RNAs and their modification status. Results: Comprehensive analysis identified not only microRNAs. Intriguingly, 57 ± 15% (mean ± S.D.) of sequencing reads mapped to non-microRNA-small RNAs including tRNA-derived fragments (tRFs), ribosomal RNA-derived fragments (rRFs) and YRNA-derived fragments (YRFs). Misincorporation (mismatch) sites identified potential base modification positions on the small RNAs, especially on tRFs, corresponding to m1A (N1-methyladenosine), m1G (N1-methylguanosine) and m2 2G (N2, N2-dimethylguanosine). We also detected mismatch sites on rRFs corresponding to known modifications on 28 and 18S rRNA. Conclusion: We found abundant non-microRNA-small RNAs in BLCA tumor samples. Small RNAs, especially tRFs and rRFs, contain modifications that can be captured as mismatch by TGIRT sequencing. Both the modifications and the non-microRNA-small RNAs should be explored as a biomarker for BLCA detection or follow-up.

Mechanical coupling of supracellular stress amplification and tissue fluidization during exit from quiescence.

Cellular quiescence is a state of reversible cell cycle arrest that is associated with tissue dormancy. Timely regulated entry into and exit from quiescence is important for processes such as tissue homeostasis, tissue repair, stem cell maintenance, developmental processes, and immunity. However, little is known about processes that control the mechanical adaption to cell behavior changes during the transition from quiescence to proliferation. Here, we show that quiescent human keratinocyte monolayers sustain an actinomyosin-based system that facilitates global cell sheet displacements upon serum-stimulated exit from quiescence. Mechanistically, exposure of quiescent cells to serum-borne mitogens leads to rapid amplification of preexisting contractile sites, leading to a burst in monolayer tension that subsequently drives large-scale displacements of otherwise motility-restricted monolayers. The stress level after quiescence exit correlates with the level of quiescence depth at the time of activation, and a critical stress magnitude must be reached to overcome the cell sheet displacement barrier. The study shows that static quiescent cell monolayers are mechanically poised for motility, and it identifies global stress amplification as a mechanism for overcoming motility restrictions in confined confluent cell monolayers.

Studies on Protein-RNA:DNA Hybrid Interactions by Microscale Thermophoresis (MST).

Microscale thermophoresis (MST) is a technology that allows for quantitative analysis of interactions between biomolecules with low sample consumption. MST uses localized temperature fields to measure the diffusion rates of the free and bound states of a fluorescently labeled protein, and to determine the dissociation constant KD by fitting of the binding isotherm with a 1:1 binding model. Here, we describe the use of MST for quantitative analysis of the interaction of the N-terminal his-tagged 6-methyladenine (m6A) reader protein YTHDF2 with m6A modified and unmodified RNA, in single-strand configuration or with RNA:DNA hybrid substrates. The described protocol is also suitable for studies of interactions with proteins binding to double-stranded RNA or DNA substrates.

ALKBH5 regulates somatic cell reprogramming in a phase-specific manner.

Establishment of the pluripotency regulatory network in somatic cells by introducing four transcription factors [octamer binding transcription factor 4 (OCT4; also known as POU5F1), sex determining region Y (SRY)-box 2 (SOX2), Kruppel-like factor 4 (KLF4) and cellular myelocytomatosis (c-MYC)] provides a promising tool for cell-based therapies in regenerative medicine. Nevertheless, the mechanisms at play when generating induced pluripotent stem cells from somatic cells are only partly understood. Here, we show that the RNA-specific N6-methyladenosine (m6A) demethylase ALKBH5 regulates somatic cell reprogramming in a stage-specific manner through its catalytic activity. Knockdown or knockout of Alkbh5 in the early reprogramming phase impairs reprogramming efficiency by reducing the proliferation rate through arresting the cells at G2/M phase and decreasing the upregulation of epithelial markers. On the other hand, ALKBH5 overexpression at the early reprogramming phase has no significant impact on reprogramming efficiency, whereas overexpression at the late phase enhances reprogramming by stabilizing Nanog transcripts, resulting in upregulated Nanog expression. Our study provides mechanistic insight into the crucial dynamic role of ALKBH5, mediated through its catalytic activity, in regulating somatic cell reprogramming at the post-transcriptional level. This article has an associated First Person interview with the first author of the paper.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer.

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

Waves of sumoylation support transcription dynamics during adipocyte differentiation.

Tight control of gene expression networks required for adipose tissue formation and plasticity is essential for adaptation to energy needs and environmental cues. However, the mechanisms that orchestrate the global and dramatic transcriptional changes leading to adipocyte differentiation remain to be fully unraveled. We investigated the regulation of nascent transcription by the sumoylation pathway during adipocyte differentiation using SLAMseq and ChIPseq. We discovered that the sumoylation pathway has a dual function in differentiation; it supports the initial downregulation of pre-adipocyte-specific genes, while it promotes the establishment of the mature adipocyte transcriptional program. By characterizing endogenous sumoylome dynamics in differentiating adipocytes by mass spectrometry, we found that sumoylation of specific transcription factors like PPARγ/RXR and their co-factors are associated with the transcription of adipogenic genes. Finally, using RXR as a model, we found that sumoylation may regulate adipogenic transcription by supporting the chromatin occurrence of transcription factors. Our data demonstrate that the sumoylation pathway supports the rewiring of transcriptional networks required for formation of functional adipocytes. This study also provides the scientists in the field of cellular differentiation and development with an in-depth resource of the dynamics of the SUMO-chromatin landscape, SUMO-regulated transcription and endogenous sumoylation sites during adipocyte differentiation.

N6-methyladenosine (m6A) depletion regulates pluripotency exit by activating signaling pathways in embryonic stem cells.

N6-methyladenosine (m6A) deposition on messenger RNA (mRNA) controls embryonic stem cell (ESC) fate by regulating the mRNA stabilities of pluripotency and lineage transcription factors (TFs) [P. J. Batista et al., Cell Stem Cell 15, 707-719 (2014); Y. Wang et al., Nat. Cell Biol. 16, 191-198 (2014); and S. Geula et al., Science 347, 1002-1006 (2015)]. If the mRNAs of these two TF groups become stabilized, it remains unclear how the pluripotency or lineage commitment decision is implemented. We performed noninvasive quantification of Nanog and Oct4 TF protein levels in reporter ESCs to define cell-state dynamics at single-cell resolution. Long-term single-cell tracking shows that immediate m6A depletion by Mettl3 knock-down in serum/leukemia inhibitory factor supports both pluripotency maintenance and its departure. This is mediated by differential and opposing signaling pathways. Increased FGF5 mRNA stability activates pErk, leading to Nanog down-regulation. FGF5-mediated coactivation of pAkt reenforces Nanog expression. In formative stem cells poised toward differentiation, m6A depletion activates both pErk and pAkt, increasing the propensity for mesendodermal lineage induction. Stable m6A depletion by Mettl3 knock-out also promotes pErk activation. Higher pErk counteracts the pluripotency exit delay exhibited by stably m6A-depleted cells upon differentiation. At single-cell resolution, we illustrate that decreasing m6A abundances activates pErk and pAkt-signaling, regulating pluripotency departure.

NEIL1 and NEIL2 DNA glycosylases modulate anxiety and learning in a cooperative manner in mice.

Oxidative DNA damage in the brain has been implicated in neurodegeneration and cognitive decline. DNA glycosylases initiate base excision repair (BER), the main pathway for oxidative DNA base lesion repair. NEIL1 and NEIL3 DNA glycosylases affect cognition in mice, while the role of NEIL2 remains unclear. Here, we investigate the impact of NEIL2 and its potential overlap with NEIL1 on behavior in knockout mouse models. Neil1-/-Neil2-/- mice display hyperactivity, reduced anxiety and improved learning. Hippocampal oxidative DNA base lesion levels are comparable between genotypes and no mutator phenotype is found. Thus, impaired canonical repair is not likely to explain the altered behavior. Electrophysiology suggests reduced axonal activation in the hippocampal CA1 region in Neil1-/-Neil2-/- mice and lack of NEIL1 and NEIL2 causes dysregulation of genes in CA1 relevant for synaptic function. We postulate a cooperative function of NEIL1 and NEIL2 in genome regulation, beyond canonical BER, modulating behavior in mice.

Alleviation of C⋠C Mismatches in DNA by the Escherichia coli Fpg Protein.

DNA polymerase III mis-insertion may, where not corrected by its 3'→ 5' exonuclease or the mismatch repair (MMR) function, result in all possible non-cognate base pairs in DNA generating base substitutions. The most thermodynamically unstable base pair, the cytosine (C)⋅C mismatch, destabilizes adjacent base pairs, is resistant to correction by MMR in Escherichia coli, and its repair mechanism remains elusive. We present here in vitro evidence that C⋅C mismatch can be processed by base excision repair initiated by the E. coli formamidopyrimidine-DNA glycosylase (Fpg) protein. The k cat for C⋅C is, however, 2.5 to 10 times lower than for its primary substrate 8-oxoguanine (oxo8G)⋅C, but approaches those for 5,6-dihydrothymine (dHT)⋅C and thymine glycol (Tg)⋅C. The K M values are all in the same range, which indicates efficient recognition of C⋅C mismatches in DNA. Fpg activity was also exhibited for the thymine (T)⋅T mismatch and for N 4- and/or 5-methylated C opposite C or T, Fpg activity being enabled on a broad spectrum of DNA lesions and mismatches by the flexibility of the active site loop. We hypothesize that Fpg plays a role in resolving C⋅C in particular, but also other pyrimidine⋅pyrimidine mismatches, which increases survival at the cost of some mutagenesis.

ALKBH7-mediated demethylation regulates mitochondrial polycistronic RNA processing.

Members of the mammalian AlkB family are known to mediate nucleic acid demethylation1,2. ALKBH7, a mammalian AlkB homologue, localizes in mitochondria and affects metabolism3, but its function and mechanism of action are unknown. Here we report an approach to site-specifically detect N1-methyladenosine (m1A), N3-methylcytidine (m3C), N1-methylguanosine (m1G) and N2,N2-dimethylguanosine (m22G) modifications simultaneously within all cellular RNAs, and discovered that human ALKBH7 demethylates m22G and m1A within mitochondrial Ile and Leu1 pre-tRNA regions, respectively, in nascent polycistronic mitochondrial RNA4-6. We further show that ALKBH7 regulates the processing and structural dynamics of polycistronic mitochondrial RNAs. Depletion of ALKBH7 leads to increased polycistronic mitochondrial RNA processing, reduced steady-state mitochondria-encoded tRNA levels and protein translation, and notably decreased mitochondrial activity. Thus, we identify ALKBH7 as an RNA demethylase that controls nascent mitochondrial RNA processing and mitochondrial activity.

Intrinsic Strand-Incision Activity of Human UNG: Implications for Nick Generation in Immunoglobulin Gene Diversification.

Uracil arises in cellular DNA by cytosine (C) deamination and erroneous replicative incorporation of deoxyuridine monophosphate opposite adenine. The former generates C → thymine transition mutations if uracil is not removed by uracil-DNA glycosylase (UDG) and replaced by C by the base excision repair (BER) pathway. The primary human UDG is hUNG. During immunoglobulin gene diversification in activated B cells, targeted cytosine deamination by activation-induced cytidine deaminase followed by uracil excision by hUNG is important for class switch recombination (CSR) and somatic hypermutation by providing the substrate for DNA double-strand breaks and mutagenesis, respectively. However, considerable uncertainty remains regarding the mechanisms leading to DNA incision following uracil excision: based on the general BER scheme, apurinic/apyrimidinic (AP) endonuclease (APE1 and/or APE2) is believed to generate the strand break by incising the AP site generated by hUNG. We report here that hUNG may incise the DNA backbone subsequent to uracil excision resulting in a 3´-α,ß-unsaturated aldehyde designated uracil-DNA incision product (UIP), and a 5´-phosphate. The formation of UIP accords with an elimination (E2) reaction where deprotonation of C2´ occurs via the formation of a C1´ enolate intermediate. UIP is removed from the 3´-end by hAPE1. This shows that the first two steps in uracil BER can be performed by hUNG, which might explain the significant residual CSR activity in cells deficient in APE1 and APE2.

Modifications and interactions at the R-loop.

R-loops are tripartite structures consisting of an RNA:DNA hybrid and a displaced single-stranded DNA [1]. They are widespread and occupy up to 5 % of the mammalian genomes [2]. R-loops have a key role in genome stability, and known functions associated with gene regulation, DNA replication, chromatin patterning, immunoglobuline gene recombination and DNA Double-strand break repair [3-7]. Novel methodology, including the application of the S9.6 antibody, have more recently led to detailed knowledge on the genome-wide distribution of the R-loops as well as the identification of the R-loop interactome [8-10]. The regulation of R-loops was recently shown to also depend on dynamic RNA-methylation, including METTL3/14 dependent 6-methylAdenines (m6As) and METTL8 dependent 3-methylCytosines (m3Cs) [11-13].

KDM4A regulates the maternal-to-zygotic transition by protecting broad H3K4me3 domains from H3K9me3 invasion in oocytes.

The importance of germline-inherited post-translational histone modifications on priming early mammalian development is just emerging1-4. Histone H3 lysine 9 (H3K9) trimethylation is associated with heterochromatin and gene repression during cell-fate change5, whereas histone H3 lysine 4 (H3K4) trimethylation marks active gene promoters6. Mature oocytes are transcriptionally quiescent and possess remarkably broad domains of H3K4me3 (bdH3K4me3)1,2. It is unknown which factors contribute to the maintenance of the bdH3K4me3 landscape. Lysine-specific demethylase 4A (KDM4A) demethylates H3K9me3 at promoters marked by H3K4me3 in actively transcribing somatic cells7. Here, we report that KDM4A-mediated H3K9me3 demethylation at bdH3K4me3 in oocytes is crucial for normal pre-implantation development and zygotic genome activation after fertilization. The loss of KDM4A in oocytes causes aberrant H3K9me3 spreading over bdH3K4me3, resulting in insufficient transcriptional activation of genes, endogenous retroviral elements and chimeric transcripts initiated from long terminal repeats during zygotic genome activation. The catalytic activity of KDM4A is essential for normal epigenetic reprogramming and pre-implantation development. Hence, KDM4A plays a crucial role in preserving the maternal epigenome integrity required for proper zygotic genome activation and transfer of developmental control to the embryo.

The Escherichia coli alkA Gene Is Activated to Alleviate Mutagenesis by an Oxidized Deoxynucleoside.

The cellular methyl donor S-adenosylmethionine (SAM) and other endo/exogenous agents methylate DNA bases non-enzymatically into products interfering with replication and transcription. An important product is 3-methyladenine (m3A), which in Escherichia coli is removed by m3A-DNA glycosylase I (Tag) and II (AlkA). The tag gene is constitutively expressed, while alkA is induced by sub-lethal concentrations of methylating agents. We previously found that AlkA exhibits activity for the reactive oxygen-induced thymine (T) lesion 5-formyluracil (fU) in vitro. Here, we provide evidence for AlkA involvement in the repair of oxidized bases by showing that the adenine (A) ⋅ T → guanine (G) ⋅ cytosine (C) mutation rate increased 10-fold in E. coli wild-type and alkA - cells exposed to 0.1 mM 5-formyl-2'-deoxyuridine (fdU) compared to a wild-type specific reduction of the mutation rate at 0.2 mM fdU, which correlated with alkA gene induction. G⋅C → A⋅T alleviation occurred without alkA induction (at 0.1 mM fdU), correlating with a much higher AlkA efficiency for fU opposite to G than for that to A. The common keto form of fU is the AlkA substrate. Mispairing with G by ionized fU is favored by its exclusion from the AlkA active site.

OXR1A, a Coactivator of PRMT5 Regulating Histone Arginine Methylation.

Oxidation resistance gene 1 (OXR1) protects cells against oxidative stress. We find that male mice with brain-specific isoform A knockout (Oxr1A-/-) develop fatty liver. RNA sequencing of male Oxr1A-/- liver indicates decreased growth hormone (GH) signaling, which is known to affect liver metabolism. Indeed, Gh expression is reduced in male mice Oxr1A-/- pituitary gland and in rat Oxr1A-/- pituitary adenoma cell-line GH3. Oxr1A-/- male mice show reduced fasting-blood GH levels. Pull-down and proximity ligation assays reveal that OXR1A is associated with arginine methyl transferase PRMT5. OXR1A-depleted GH3 cells show reduced symmetrical dimethylation of histone H3 arginine 2 (H3R2me2s), a product of PRMT5 catalyzed methylation, and chromatin immunoprecipitation (ChIP) of H3R2me2s shows reduced Gh promoter enrichment. Finally, we demonstrate with purified proteins that OXR1A stimulates PRMT5/MEP50-catalyzed H3R2me2s. Our data suggest that OXR1A is a coactivator of PRMT5, regulating histone arginine methylation and thereby GH production within the pituitary gland.

Deletion of Endonuclease V suppresses chemically induced hepatocellular carcinoma.

Endonuclease V (EndoV) is a conserved inosine-specific ribonuclease with unknown biological function. Here, we present the first mouse model lacking EndoV, which is viable without visible abnormalities. We show that endogenous murine EndoV cleaves inosine-containing RNA in vitro, nevertheless a series of experiments fails to link an in vivo function to processing of such transcripts. As inosine levels and adenosine-to-inosine editing often are dysregulated in hepatocellular carcinoma (HCC), we chemically induced HCC in mice. All mice developed liver cancer, however, EndoV-/- tumors were significantly fewer and smaller than wild type tumors. Opposed to human HCC, adenosine deaminase mRNA expression and site-specific editing were unaltered in our model. Loss of EndoV did not affect editing levels in liver tumors, however mRNA expression of a selection of cancer related genes were reduced. Inosines are also found in certain tRNAs and tRNAs are cleaved during stress to produce signaling entities. tRNA fragmentation was dysregulated in EndoV-/- livers and apparently, inosine-independent. We speculate that the inosine-ribonuclease activity of EndoV is disabled in vivo, but RNA binding allowed to promote stabilization of transcripts or recruitment of proteins to fine-tune gene expression. The EndoV-/- tumor suppressive phenotype calls for related studies in human HCC.