Biochemical and structural basis for YTH domain of human YTHDC1 binding to methylated adenine in DNA

YTHDF1 facilitates PRC1-mediated H2AK119ub in human ES cells

Polycomb repressive complexes (PRCs) play critical roles in cell fate decisions during normal development as well as disease progression through mediating histone modifications such as H3K27me3 and H2AK119ub. How exactly PRCs recruited to chromatin remains to be fully illuminated. Here, we report that YTHDF1, the N6-methyladenine (m6 A) RNA reader that was previously known to be mainly cytoplasmic, associates with RNF2, a PRC1 protein that mediates H2AK119ub in human embryonic stem cells (hESCs). A portion of YTHDF1 localizes in the nuclei and associates with RNF2/H2AK119ub on a subset of gene loci related to neural development functions. Knock-down YTHDF1 attenuates H2AK119ub modification on these genes and promotes neural differentiation in hESCs. Our findings provide a noncanonical mechanism that YTHDF1 participates in PRC1 functions in hESCs.

Methylation across the central dogma in health and diseases: new therapeutic strategies

The proper transfer of genetic information from DNA to RNA to protein is essential for cell-fate control, development, and health. Methylation of DNA, RNAs, histones, and non-histone proteins is a reversible post-synthesis modification that finetunes gene expression and function in diverse physiological processes. Aberrant methylation caused by genetic mutations or environmental stimuli promotes various diseases and accelerates aging, necessitating the development of therapies to correct the disease-driver methylation imbalance. In this Review, we summarize the operating system of methylation across the central dogma, which includes writers, erasers, readers, and reader-independent outputs. We then discuss how dysregulation of the system contributes to neurological disorders, cancer, and aging. Current small-molecule compounds that target the modifiers show modest success in certain cancers. The methylome-wide action and lack of specificity lead to undesirable biological effects and cytotoxicity, limiting their therapeutic application, especially for diseases with a monogenic cause or different directions of methylation changes. Emerging tools capable of site-specific methylation manipulation hold great promise to solve this dilemma. With the refinement of delivery vehicles, these new tools are well positioned to advance the basic research and clinical translation of the methylation field.

Everything OLD is new again: How structural, functional, and bioinformatic advances have redefined a neglected nuclease family

Overcoming lysogenization defect (OLD) proteins are a conserved family of ATP-powered nucleases that function in anti-phage defense. Recent bioinformatic, genetic, and crystallographic studies have yielded new insights into the structure, function, and evolution of these enzymes. Here we review these developments and propose a new classification scheme to categorize OLD homologs that relies on gene neighborhoods, biochemical properties, domain organization, and catalytic machinery. This taxonomy reveals important similarities and differences between family members and provides a blueprint to contextualize future in vivo and in vitro findings. We also detail how OLD nucleases are related to PARIS and Septu anti-phage defense systems and discuss important mechanistic questions that remain unanswered.

Biochemical and structural characterization of the first-discovered metazoan DNA cytosine-N4 methyltransferase from the bdelloid rotifer Adineta vaga

Much is known about the generation, removal, and roles of 5-methylcytosine (5mC) in eukaryote DNA, and there is a growing body of evidence regarding N6-methyladenine, but very little is known about N4-methylcytosine (4mC) in the DNA of eukaryotes. The gene for the first metazoan DNA methyltransferase generating 4mC (N4CMT) was reported and characterized recently by others, in tiny freshwater invertebrates called bdelloid rotifers. Bdelloid rotifers are ancient, apparently asexual animals, and lack canonical 5mC DNA methyltransferases. Here, we characterize the kinetic properties and structural features of the catalytic domain of the N4CMT protein from the bdelloid rotifer Adineta vaga. We find that N4CMT generates high-level methylation at preferred sites, (a/c)CG(t/c/a), and low-level methylation at disfavored sites, exemplified by ACGG. Like the mammalian de novo 5mC DNA methyltransferase 3A/3B (DNMT3A/3B), N4CMT methylates CpG dinucleotides on both DNA strands, generating hemimethylated intermediates and eventually fully methylated CpG sites, particularly in the context of favored symmetric sites. In addition, like DNMT3A/3B, N4CMT methylates non-CpG sites, mainly CpA/TpG, though at a lower rate. Both N4CMT and DNMT3A/3B even prefer similar CpG-flanking sequences. Structurally, the catalytic domain of N4CMT closely resembles the Caulobacter crescentus cell cycle-regulated DNA methyltransferase. The symmetric methylation of CpG, and similarity to a cell cycle-regulated DNA methyltransferase, together suggest that N4CMT might also carry out DNA synthesis-dependent methylation following DNA replication.

Characterization of epitranscriptome reader proteins experimentally and in silico: Current knowledge and future perspectives beyond the YTH domain

To date, over 150 chemical modifications to the four canonical RNA bases have been discovered, known collectively as the epitranscriptome. Many of these modifications have been implicated in a variety of cellular processes and disease states. Additional work has been done to identify proteins known as "readers" that selectively interact with RNAs that contain specific chemical modifications. Protein interactomes with N6-methyladenosine (m⁶A), N1-methyladenosine (m¹A), N5-methylcytosine (m⁵C), and 8-oxo-7,8-dihydroguanosine (8-oxoG) have been determined, mainly through experimental advances in proteomics techniques. However, relatively few proteins have been confirmed to bind directly to RNA containing these modifications. Furthermore, for many of these protein readers, the exact binding mechanisms as well as the exclusivity for recognition of modified RNA species remain elusive, leading to questions regarding their roles within different cellular processes. In the case of the YT-521B homology (YTH) family of proteins, both experimental and in silico techniques have been leveraged to provide valuable biophysical insights into the mechanisms of m⁶A recognition at atomic resolution. To date, the YTH family is one of the best characterized classes of readers. Here, we review current knowledge about epitranscriptome recognition of the YTH domain proteins from previously published experimental and computational studies. We additionally outline knowledge gaps for proteins beyond the well-studied human YTH domains and the current in silico techniques and resources that can enable investigation of protein interactions with modified RNA outside of the YTH-m⁶A context.

The potential role of m6A reader YTHDF1 as diagnostic biomarker and the signaling pathways in tumorigenesis and metastasis in pan-cancer

m6A is an important RNA methylation in progression of various human cancers. As the m6A reader protein, YTHDF1 is reported to accelerate m6A-modified mRNAs translation in cytoplasm. It is highly expressed in various human cancers and contributes to the progression and metastasis of cancers. YTHDF1 was closely associated with poor prognosis and also used as a molecular marker for clinical diagnosis or therapy in human cancers. It has been reported to promote chemoresistance to Adriamycin, Cisplatin and Olaparib by increasing mRNA stability of its target molecule. Moreover, it contributes to CSC-like characteristic of tumor cells and inducing the antitumor immune microenvironment. Here, we reviewed the clinical diagnostic and prognostic values of YTHDF1, as well as the molecular mechanisms of YTHDF1 in progression and metastasis of human cancers.

High-precision mapping reveals rare N6-deoxyadenosine methylation in the mammalian genome

N⁶-deoxyadenosine methylation (6mA) is the most widespread type of DNA modification in prokaryotes and is also abundantly distributed in some unicellular eukaryotes. However, 6mA levels are remarkably low in mammals. The lack of a precise and comprehensive mapping method has hindered more advanced investigations of 6mA. Here, we report a new method MM-seq (modification-induced mismatch sequencing) for genome-wide 6mA mapping based on a novel detection principle. We found that modified DNA bases are prone to form a local open region that allows capture by antibody, for example, via a DNA breathing or base-flipping mechanism. Specified endonuclease or exonuclease can recognize the antibody-stabilized mismatch-like structure and mark the exact modified sites for sequencing readout. Using this method, we examined the genomic positions of 6mA in bacteria (E. coli), green algae (C. reinhardtii), and mammalian cells (HEK239T, Huh7, and HeLa cells). In contrast to bacteria and green algae, human cells possess a very limited number of 6mA sites which are sporadically distributed across the genome of different cell types. After knocking out the RNA m⁶A methyltransferase METTL3 in mouse ES cells, 6mA becomes mostly diminished. Our results imply that rare 6mA in the mammalian genome is introduced by RNA m⁶A machinery via a non-targeted mechanism.

Advances in detecting N6-methyladenosine modification in circRNAs

Circular RNAs (circRNAs) are a class of noncoding RNAs with covalently single-stranded closed loop structures derived from back-splicing event of linear precursor mRNAs (pre-mRNAs). N6-methyladenosine (m⁶A), the most abundant epigenetic modification in eukaryotic RNAs, has been shown to play a crucial role in regulating the fate and biological function of circRNAs, and thus affecting various physiological and pathological processes. Accurate identification of m⁶A modification in circRNAs is an essential step to fully elucidate the crosstalk between m⁶A and circRNAs. In recent years, the rapid development of high-throughput sequencing technology and bioinformatic methodology has propelled the establishment of a multitude of approaches to detect circRNAs and m⁶A modification, including in vitro-based and in silico methods. Based on this, the research community has started on a new journey to develop methods for identification of m⁶A modification in circRNAs. In this review, we provide a comprehensive review and evaluation of the existing methods responsible for detecting circRNAs, m⁶A modification, and especially, m⁶A modification in circRNAs, which mainly focused on those developed based on high-throughput technologies and methodology of bioinformatics. This handy reference can help researchers figure out towards which direction this field will go.

Enzymatic Characterization of In Vitro Activity of RNA Methyltransferase PCIF1 on DNA

PCIF1 and FTO are a pair of human mRNA cap-specific modification enzymes that have opposing activities. PCIF1 adds a methyl group to the N6-position of 2'O-methyladenosine (Am), generating N6, 2'O-dimethyladenosine (m6Am), when Am is the cap-proximal nucleotide. FTO removes the N6-methyl group from m6Am. In addition, FTO has a demethylase activity on a broad spectrum of various RNA substrates, as well as on DNA N6-methyldeoxyadenosine (m6dA). While the existence of m6dA in mammalian DNA remains controversial, we show here that PCIF1 has significant methylation activity on single stranded DNA deoxyadenosine, double stranded RNA/DNA hybrids, and double stranded DNA, though with lower catalytic efficiency than that on its preferred RNA substrate. PCIF1 has activities in the order ssRNA > RNA/DNA hybrid > ssDNA > dsDNA. We discuss the implications of PCIF1 generation, and FTO removal, of DNA adenine methylation.

Genome-Wide Identification, Classification and Expression Analysis of m6A Gene Family in Solanum lycopersicum

Advanced knowledge of messenger RNA (mRNA) N⁶-methyladenosine (m⁶A) and DNA N⁶-methyldeoxyadenosine (6 mA) redefine our understanding of these epigenetic modifications. Both m⁶A and 6mA carry important information for gene regulation, and the corresponding catalytic enzymes sometimes belong to the same gene family and need to be distinguished. However, a comprehensive analysis of the m⁶A gene family in tomato remains obscure. Here, 24 putative m⁶A genes and their family genes in tomato were identified and renamed according to BLASTP and phylogenetic analysis. Chromosomal location, synteny, phylogenetic, and structural analyses were performed, unravelling distinct evolutionary relationships between the MT-A70, ALKBH, and YTH protein families, respectively. Most of the 24 genes had extensive tissue expression, and 9 genes could be clustered in a similar expression trend. Besides, SlYTH1 and SlYTH3A showed a different expression pattern in leaf and fruit development. Additionally, qPCR data revealed the expression variation under multiple abiotic stresses, and LC-MS/MS determination exhibited that the cold stress decreased the level of N⁶ 2′-O dimethyladenosine (m⁶Am). Notably, the orthologs of newly identified single-strand DNA (ssDNA) 6mA writer–eraser–reader also existed in the tomato genome. Our study provides comprehensive information on m⁶A components and their family proteins in tomato and will facilitate further functional analysis of the tomato N⁶-methyladenosine modification genes.

Evolutionary History of RNA Modifications at N6-Adenosine Originating from the R-M System in Eukaryotes and Prokaryotes

Simple Summary

The m⁶A is the most abundant and well-studied modification of mRNA, and plays an important role in transcription and translation. It is known to be evolutionarily conserved machinery present in the last eukaryotic common ancestor (LECA). The writers and erasers responsible for adding or removing m⁶A belong to specific protein families, respectively, suggesting that these members are evolutionarily related. However, only some of these mRNA m⁶A modification-associated proteins have been studied from an evolutionary perspective, while there has been no comprehensive and systematic analysis of the distributions and evolutionary history of N6mA-associated proteins in the three kingdoms of life. In this study, we identified orthologues of all the reported N6mA-associated proteins in 88 organisms from three kingdoms of life and comprehensively reconstructed the evolutionary history of the RNA N6mA modification machinery. The results demonstrate that RNA N6mA-MTases are derived from at least two different types of prokaryotic DNA MTases (class α and β MTases). As the m⁶A reader, YTH proteins may be acquired by LECA from an individual prokaryotic YTH-domain protein that evolved from the N-terminals of an R-M system endonuclease. In addition, the origin of eukaryotic ALKBH family proteins is inferred to be driven by at least two occasions of independent HTG from the bacterial ALKB family.

Abstract

Methylation at the N6-position of adenosine (N6mA) on mRNA (m⁶A) is one of the most widespread, highly selective and dynamically regulated RNA modifications and plays an important role in transcription and translation. In the present study, a comprehensive analysis of phylogenetic relationships, conserved domain sequence characteristics and protein structure comparisons were employed to explore the distribution of RNA N6mA modification (m⁶A, m^6,6A, m⁶Am, m^{6, 6}Am and m⁶t⁶A)-associated proteins (writers, readers and erasers) in three kingdoms of life and reveal the evolutionary history of these modifications. These findings further confirmed that the restriction-modification (R-M) system is the origin of DNA and RNA N6mA modifications. Among them, the existing mRNA m⁶A modification system derived from the last eukaryotic common ancestor (LECA) is the evolutionary product of elements from the last universal common ancestor (LUCA) or driven by horizontal gene transfer (HGT) from bacterial elements. The subsequent massive gene gains and losses contribute to the development of unique and diverse functions in distinct species. Particularly, RNA methyltransferases (MTases) as the writer responsible for adding N6mA marks on mRNA and ncRNAs may have evolved from class α and β prokaryotic “orphan” MTases originating from the R-M system. The reader, YTH proteins that specifically recognize the m⁶A deposit, may be acquired by LECA from an individual prokaryotic YTH-domain protein that evolved from N-terminals of an R-M system endonuclease. The eraser, which emerged from the ALKB family (ALKBH5 and FTO) in eukaryotes, may be driven by independent HTG from bacterial ALKB proteins. The evolutionary history of RNA N6mA modifications was inferred in the present study, which will deepen our understanding of these modifications in different species.

m6A: An Emerging Role in Programmed Cell Death

Programmed cell death is an active extinction process, including autophagy, ferroptosis, pyroptosis, apoptosis, and necroptosis. m6A is a reversible RNA modification which undergoes methylation under the action of methylases (writers), and is demethylated under the action of demethylases (erasers). The RNA base site at which m6A is modified is recognized by specialized enzymes (readers) which regulate downstream RNA translation, decay, and stability. m6A affects many aspects of mRNA metabolism, and also plays an important role in promoting the maturation of miRNA, the translation and degradation of circRNA, and the stability of lncRNA. The regulatory factors including writers, erasers and readers promote or inhibit programmed cell death via up-regulating or down-regulating downstream targets in a m6A-dependent manner to participate in the process of disease. In this review, we summarize the functions of m6A with particular reference to its role in programmed cell death.

A plant-like mechanism coupling m6A reading to polyadenylation safeguards transcriptome integrity and developmental gene partitioning in Toxoplasma

Correct 3’end processing of mRNAs is one of the regulatory cornerstones of gene expression. In a parasite that must adapt to the regulatory requirements of its multi-host life style, there is a need to adopt additional means to partition the distinct transcriptional signatures of the closely and tandemly arranged stage-specific genes. In this study, we report our findings in T. gondii of an m6A-dependent 3’end polyadenylation serving as a transcriptional barrier at these loci. We identify the core polyadenylation complex within T. gondii and establish CPSF4 as a reader for m6A-modified mRNAs, via a YTH domain within its C-terminus, a feature which is shared with plants. We bring evidence of the specificity of this interaction both biochemically, and by determining the crystal structure at high resolution of the T. gondii CPSF4-YTH in complex with an m6A-modified RNA. We show that the loss of m6A, both at the level of its deposition or its recognition is associated with an increase in aberrantly elongated chimeric mRNAs emanating from impaired transcriptional termination, a phenotype previously noticed in the plant model Arabidopsis thaliana. Nanopore direct RNA sequencing shows the occurrence of transcriptional read-through breaching into downstream repressed stage-specific genes, in the absence of either CPSF4 or the m6A RNA methylase components in both T. gondii and A. thaliana. Taken together, our results shed light on an essential regulatory mechanism coupling the pathways of m6A metabolism directly to the cleavage and polyadenylation processes, one that interestingly seem to serve, in both T. gondii and A. thaliana, as a guardian against aberrant transcriptional read-throughs.

Recognition of N6-Methyladenosine by the YTHDC1 YTH Domain Studied by Molecular Dynamics and NMR Spectroscopy: The Role of Hydration

The YTH domain of YTHDC1 belongs to a class of protein "readers", recognizing the N6-methyladenosine (m⁶A) chemical modification in mRNA. Static ensemble-averaged structures revealed details of N6-methyl recognition via a conserved aromatic cage. Here, we performed molecular dynamics (MD) simulations along with nuclear magnetic resonance (NMR) and isothermal titration calorimetry (ITC) to examine how dynamics and solvent interactions contribute to the m⁶A recognition and negative selectivity toward an unmethylated substrate. The structured water molecules surrounding the bound RNA and the methylated substrate's ability to exclude bulk water molecules contribute to the YTH domain's preference for m⁶A. Intrusions of bulk water deep into the binding pocket disrupt binding of unmethylated adenosine. The YTHDC1's preference for the 5'-Gm⁶A-3' motif is partially facilitated by a network of water-mediated interactions between the 2-amino group of the guanosine and residues in the m⁶A binding pocket. The 5'-Im⁶A-3' (where I is inosine) motif can be recognized too, but disruption of the water network lowers affinity. The D479A mutant also disrupts the water network and destabilizes m⁶A binding. Our interdisciplinary study of the YTHDC1 protein-RNA complex reveals an unusual physical mechanism by which solvent interactions contribute toward m⁶A recognition.

The language of chromatin modification in human cancers

The genetic information of human cells is stored in the context of chromatin, which is subjected to DNA methylation and various histone modifications. Such a 'language' of chromatin modification constitutes a fundamental means of gene and (epi)genome regulation, underlying a myriad of cellular and developmental processes. In recent years, mounting evidence has demonstrated that miswriting, misreading or mis-erasing of the modification language embedded in chromatin represents a common, sometimes early and pivotal, event across a wide range of human cancers, contributing to oncogenesis through the induction of epigenetic, transcriptomic and phenotypic alterations. It is increasingly clear that cancer-related metabolic perturbations and oncohistone mutations also directly impact chromatin modification, thereby promoting cancerous transformation. Phase separation-based deregulation of chromatin modulators and chromatin structure is also emerging to be an important underpinning of tumorigenesis. Understanding the various molecular pathways that underscore a misregulated chromatin language in cancer, together with discovery and development of more effective drugs to target these chromatin-related vulnerabilities, will enhance treatment of human malignancies.

Human MettL3-MettL14 RNA adenine methyltransferase complex is active on double-stranded DNA containing lesions

MettL3-MettL14 methyltransferase complex has been studied widely for its role in RNA adenine methylation. This complex is also recruited to UV- and X-ray exposed DNA damaged sites, and its methyltransfer activity is required for subsequent DNA repair, though in theory this could result from RNA methylation of short transcripts made at the site of damage. We report here that MettL3-MettL14 is active in vitro on double-stranded DNA containing a cyclopyrimidine dimer – a major lesion of UV radiation-induced products – or an abasic site or mismatches. Furthermore, N6-methyladenine (N6mA) decreases misincorporation of 8-oxo-guanine (8-oxoG) opposite to N6mA by repair DNA polymerases. When 8-oxoG is nevertheless incorporated opposite N6mA, the methylation inhibits N6mA excision from the template (correct) strand by the adenine DNA glycosylase (MYH), implying that the methylation decreases inappropriate misrepair. Finally, we observed that the N6mA reader domain of YTHDC1, which is also recruited to sites of DNA damage, binds N6mA that is located across from a single-base gap between two canonical DNA helices. This YTHDC1 complex with a gapped duplex is structurally similar to DNA complexes with FEN1 and GEN1 – two members of the nuclease family that act in nucleotide excision repair, mismatch repair and homologous recombination, and which incise distinct non-B DNA structures. Together, the parts of our study provide a plausible mechanism for N6mA writer and reader proteins acting directly on lesion-containing DNA, and suggest in vivo experiments to test the mechanisms involving methylation of adenine.

The m6A reader YTHDC2 inhibits lung adenocarcinoma tumorigenesis by suppressing SLC7A11-dependent antioxidant function

The biological functions of N6-methyladenosine (m⁶A) RNA methylation are mainly dependent on the reader; however, its role in lung tumorigenesis remains unclear. Here, we have demonstrated that the m⁶A reader YT521-B homology domain containing 2 (YTHDC2) is frequently suppressed in lung adenocarcinoma (LUAD). Downregulation of YTHDC2 was associated with poor clinical outcome of LUAD. YTHDC2 decreased tumorigenesis in a spontaneous LUAD mouse model. Moreover, YTHDC2 exhibited antitumor activity in human LUAD cells. Mechanistically, YTHDC2, via its m⁶A-recognizing YTH domain, suppressed cystine uptake and blocked the downstream antioxidant program. Administration of cystine downstream antioxidants to pulmonary YTHDC2-overexpressing mice rescued lung tumorigenesis. Furthermore, solute carrier 7A11 (SLC7A11), the catalytic subunit of system X_C⁻, was identified to be the direct target of YTHDC2. YTHDC2 destabilized SLC7A11 mRNA in an m⁶A-dependent manner because YTHDC2 preferentially bound to m⁶A-modified SLC7A11 mRNA and thereafter promoted its decay. Clinically, a large proportion of acinar LUAD subtype cases exhibited simultaneous YTHDC2 downregulation and SLC7A11 elevation. Patient-derived xenograft (PDX) mouse models generated from acinar LUAD showed sensitivity to system X_C⁻ inhibitors. Collectively, the promotion of cystine uptake via the suppression of YTHDC2 is critical for LUAD tumorigenesis, and blocking this process may benefit future treatment.

•

The m⁶A reader YTHDC2 is frequently suppressed in LUAD and indicates poor prognosis.

•

YTHDC2 suppresses the antioxidant function of system X_C⁻ via its m⁶A reading domain.

•

The mRNA encoding SLC7A11 is a direct target of YTHDC2.

•

YTHDC2 preferentially accelerates the decay of m⁶A-methylated SLC7A11 mRNA.

•

LUAD with YTHDC2 suppression is sensitive to system X_C⁻ inhibitors.

A Role for N6-Methyladenine in DNA Damage Repair

The leading cause of mutation due to oxidative damage is 8-oxo-2'-deoxyguanosine (8-oxoG) mispairing with adenine (Ade), which can occur in two ways. First, guanine of a G:C DNA base pair can be oxidized. If not repaired in time, DNA polymerases can mispair Ade with 8-oxoG in the template. This 8-oxoG:A can be repaired by enzymes that remove Ade opposite to template 8-oxoG, or 8-oxoG opposite to Cyt. Second, free 8-oxo-dGTP can be misincorporated by DNA polymerases into DNA opposite template Ade. However, there is no known repair activity that removes 8-oxoG opposite to template Ade. We suggest that a major role of N6-methyladenine in mammalian DNA is minimizing incorporation of 8-oxoG opposite to Ade by DNA polymerases following adduct formation.

The Potential Roles of RNA N6-Methyladenosine in Urological Tumors

N6-methyladenosine (m⁶A) is regarded as the most abundant, prevalent and conserved internal mRNA modification in mammalian cells. M⁶A can be catalyzed by m⁶A methyltransferases METTL3, METTL14 and WTAP (writers), reverted by demethylases ALKBH5 and FTO (erasers), and recognized by m⁶A -binding proteins such as YTHDF1/2/3, IGF2BP1/2/3 and HNRNPA2B1 (readers). Emerging evidence suggests that m⁶A modification is significant for regulating many biological and cellular processes and participates in the pathological development of various diseases, including tumors. This article reviews recent studies on the biological function of m⁶A modification and the methylation modification of m⁶A in urological tumors.