It has been widely accepted that 5-methylcytosine is the only form of DNA methylation in mammalian genomes. Here we identify N6-methyladenine as another form of DNA modification in mouse embryonic stem cells. Alkbh1 encodes a demethylase for N6-methyladenine. An increase of N6-methyladenine levels in Alkbh1-deficient cells leads to transcriptional silencing. N6-methyladenine deposition is inversely correlated with the evolutionary age of LINE-1 transposons; its deposition is strongly enriched at young (<1.5 million years old) but not old (>6 million years old) L1 elements. The deposition of N6-methyladenine correlates with epigenetic silencing of such LINE-1 transposons, together with their neighbouring enhancers and genes, thereby resisting the gene activation signals during embryonic stem cell differentiation. As young full-length LINE-1 transposons are strongly enriched on the X chromosome, genes located on the X chromosome are also silenced. Thus, N6-methyladenine developed a new role in epigenetic silencing in mammalian evolution distinct from its role in gene activation in other organisms. Our results demonstrate that N6-methyladenine constitutes a crucial component of the epigenetic regulation repertoire in mammalian genomes.

Peroxisomes are components of virtually all eukaryotic cells. While much is known about peroxisomal matrix protein import, our understanding of how peroxisomal membrane proteins (PMPs) are targeted and inserted into the peroxisome membrane is extremely limited. Here, we show that PEX19 binds a broad spectrum of PMPs, displays saturable PMP binding, and interacts with regions of PMPs required for their targeting to peroxisomes. Furthermore, mislocalization of PEX19 to the nucleus leads to nuclear accumulation of newly synthesized PMPs. At steady state, PEX19 is bimodally distributed between the cytoplasm and peroxisome, with most of the protein in the cytoplasm. We propose that PEX19 may bind newly synthesized PMPs and facilitate their insertion into the peroxisome membrane. This hypothesis is supported by the observation that the loss of PEX19 results in degradation of PMPs and/or mislocalization of PMPs to the mitochondrion.

Liver metastasis, the leading cause of colorectal cancer mortality, exhibits a highly heterogeneous and suppressive immune microenvironment. Here, we sequenced 97 matched samples by using single-cell RNA sequencing and spatial transcriptomics. Strikingly, the metastatic microenvironment underwent remarkable spatial reprogramming of immunosuppressive cells such as MRC1+ CCL18+ M2-like macrophages. We further developed scMetabolism, a computational pipeline for quantifying single-cell metabolism, and observed that those macrophages harbored enhanced metabolic activity. Interestingly, neoadjuvant chemotherapy could block this status and restore the antitumor immune balance in responsive patients, whereas the nonresponsive patients deteriorated into a more suppressive one. Our work described the immune evolution of metastasis and uncovered the black box of how tumors respond to neoadjuvant chemotherapy. Significance: We present a single-cell and spatial atlas of colorectal liver metastasis and found the highly metabolically activated MRC1+ CCL18+ M2-like macrophages in metastatic sites. Efficient neoadjuvant chemotherapy can slow down such metabolic activation, raising the possibility to target metabolism pathways in metastasis. This article is highlighted in the In This Issue feature, p. 1

BackgroundMost published genome sequences are drafts, and most are dominated by computational gene prediction. Draft genomes typically incorporate considerable sequence data that are not assigned to chromosomes, and predicted genes without quality confidence measures. The current Actinidia chinensis (kiwifruit) ‘Hongyang’ draft genome has 164 Mb of sequences unassigned to pseudo-chromosomes, and omissions have been identified in the gene models.ResultsA second genome of an A. chinensis (genotype Red5) was fully sequenced. This new sequence resulted in a 554.0 Mb assembly with all but 6 Mb assigned to pseudo-chromosomes. Pseudo-chromosomal comparisons showed a considerable number of translocation events have occurred following a whole genome duplication (WGD) event some consistent with centromeric Robertsonian-like translocations. RNA sequencing data from 12 tissues and ab initio analysis informed a genome-wide manual annotation, using the WebApollo tool. In total, 33,044 gene loci represented by 33,123 isoforms were identified, named and tagged for quality of evidential support. Of these 3114 (9.4%) were identical to a protein within ‘Hongyang’ The Kiwifruit Information Resource (KIR v2). Some proportion of the differences will be varietal polymorphisms. However, as most computationally predicted Red5 models required manual re-annotation this proportion is expected to be small. The quality of the new gene models was tested by fully sequencing 550 cloned ‘Hort16A’ cDNAs and comparing with the predicted protein models for Red5 and both the original ‘Hongyang’ assembly and the revised annotation from KIR v2. Only 48.9% and 63.5% of the cDNAs had a match with 90% identity or better to the original and revised ‘Hongyang’ annotation, respectively, compared with 90.9% to the Red5 models.ConclusionsOur study highlights the need to take a cautious approach to draft genomes and computationally predicted genes. Our use of the manual annotation tool WebApollo facilitated manual checking and correction of gene models enabling improvement of computational prediction. This utility was especially relevant for certain types of gene families such as the EXPANSIN like genes. Finally, this high quality gene set will supply the kiwifruit and general plant community with a new tool for genomics and other comparative analysis.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-4656-3) contains supplementary material, which is available to authorized users.