N6-methyladenosine methyltransferase METTL3 affects the phenotype of cerebral arteriovenous malformation via modulating Notch signaling pathway
Journal of Biomedical Science volume 27, Article number: 62 (2020)
Cerebral arteriovenous malformation (AVM) is a serious life-threatening congenital cerebrovascular disease. Specific anatomical features, such as nidus size, location, and venous drainage, have been validated to affect treatment outcomes. Until recently, molecular biomarkers and corresponding molecular mechanism related to anatomical features and treatment outcomes remain unknown.
RNA N6-methyladenosine (m6A) Methyltransferase METTL3 was identified as a differentially expressed gene in groups with different lesion sizes by analyzing the transcriptome sequencing (RNA-seq) data. Tube formation and wound healing assays were performed to investigate the effect of METTL3 on angiogenesis. In addition, Methylated RNA Immunoprecipitation Sequencing technology (MeRIP-seq) was performed to screen downstream targets of METTL3 in endothelial cells and to fully clarify the specific underlying molecular mechanisms affecting the phenotype of cerebral AVM.
In the current study, we found that the expression level of METTL3 was reduced in the larger pathological tissues of cerebral AVMs. Moreover, knockdown of METTL3 significantly affected angiogenesis of the human endothelial cells. Mechanistically, down-regulation of METTL3 reduced the level of heterodimeric Notch E3 ubiquitin ligase formed by DTX1 and DTX3L, thereby continuously activating the Notch signaling pathway. Ultimately, the up-regulated downstream genes of Notch signaling pathway dramatically affected the angiogenesis of endothelial cells. In addition, we demonstrated that blocking Notch pathway with DAPT could restore the phenotype of METTL3 deficient endothelial cells.
Our findings revealed the mechanism by which m6A modification regulated the angiogenesis and might provide potential biomarkers to predict the outcome of treatment, as well as provide suitable pharmacological targets for preventing the formation and progression of cerebral AVM.
Cerebral arteriovenous malformation (AVM) is a fatal congenital vascular disease. Maldevelopment or lack of capillary network in the lesion regions leads to direct connection between cerebral arteries and veins, resulting in a series of brain hemodynamic disorders . In addition to conservative treatment, microsurgical resection, embolization, stereotactic radiosurgery (SRS), or a combination of these modalities have been evolved to treat cerebral AVMs . Specific anatomical features of cerebral AVMs form the basis of several commonly used grading scales that have been validated to predict treatment outcomes. The most widely used Spetzler-Martin grading of cerebral AVMs were categorized based on nidus size (maximum diameter, MD), location (in or not in the eloquent cortex, which refers to the regions that control motor, sensory, visual, and language functions), and the presence or absence of deep venous drainage. It can be used not only to predict the outcome of microsurgical treatment, but also to predict the SRS outcome [3, 4]. Among these anatomical features, it is particularly noteworthy that nidus size was also significantly associated with seizures and other clinical symptoms. However, genetic and molecular biomarkers associated with anatomical features and risks of treatment outcomes have not yet been identified, and a large amount of research is still needed to elucidate the progression and pathogenesis of cerebral AVM.
Recently, it is well established that N6-methyladenosine (m6A) is the most common and abundant RNA molecular modification in eukaryotes . And the catalytic subunit METTL3 and METTL14, along with the regulate subunit WTAP and KIAA1429, form the core methyltransferase complex to catalyze the m6A modification of adenosine on RNA [6,7,8,9,10]. Accumulating evidences have shown that m6A modifications exert great influence on the multiple metabolic processes of RNA, such as shearing, nuclear transport, translation ability, stability, and transcription [11,12,13,14,15,16,17,18,19], and play critical roles in many bioprocesses, including circadian rhythm, DNA damage response, sex determination, developmental arrest, neuronal disorder, infectious diseases and tumorigenesis [20, 21]. Especially noteworthy is that METTL3 deficiency results in continuous activation of Notch signaling in dorsal aorta endothelial cells, eventually blocking the endothelial-to-haematopoietic transition . In addition, Notch signaling pathway has complex and context-dependent effects on angiogenesis [23, 24], and has been reported to be activated in human cerebral AVMs [25, 26]. More interestingly, both gain and loss of Notch signaling can result in AVM formation [27,28,29,30,31]. However, it has not yet been determined whether METTL3 is implicated in the formation and progression of human cerebral AVMs by regulating m6A modification.
Here, we reported that the expression level of METTL3 was down-regulated in larger cerebral AVMs lesions. Moreover, knockdown of METTL3 significantly affected the tube formation and migration of endothelial cells, suggesting the properties of vascular endothelial cell can be modulated by m6A modification. Therefore, METTL3 might be used as a biomarker linked to the size of AVM lesions to predict the treatment outcomes. In order to fully clarify the specific underlying molecular mechanism, we performed RNA transcriptome sequencing (RNA-seq) and Methylated RNA Immunoprecipitation Sequencing technology (MeRIP-seq) to profile the changes of transcriptome-wide m6A modification sites and mRNA expression level in METTL3 deficient endothelial cells. Finally, we determined that DTX3L, the direct downstream target of METTL3, synergized with DTX1 to modulate the Notch signaling pathway to regulate angiogenesis in vitro. Our findings might contribute to understanding in more detail how arteriovenous malformations develop, and also might provide suitable targets for pharmacological treatment to prevent the progression and enlargement of cerebral AVMs in the future.
Materials and methods
Patients and samples
Detailed information on patient recruitment and sample preparation could be found in our previous study . In Brief, cerebral AVMs samples were collected from consecutive patients undergoing surgical treatment. The clinical diagnosis of cerebral AVM was confirmed using digital subtraction angiography. The size and location (in or not in the eloquent cortex) of the lesion, the presence or absence of deep venous drainage and Spetzler-Martin grade were determined based on the imaging results. Finally, from September 2016 to November 2017, a total of 57 patients were enrolled in this study. Informed consents were obtained from all patients, and this study was approved by the institutional review board of Beijing Tiantan Hospital, Capital Medical University. The baseline characteristics of all patients are summarized in Table S1.
The Human Umbilical Vein Endothelial Cells (HUVECs) were purchased from ScienCell (Carlsbad, CA), and cultured in Endothelial Cell Medium (ECM, ScienCell) supplemented with 5% fetal bovine serum (FBS, Gibco), 100 U/ml penicillin and 100 μg/ml streptomycin.
The siRNAs were designed and synthesized by RIBOBIO (Table S2). Endothelial cells were transfected with siRNAs using Lipofectamine RNAiMAX (Invitrogen). After transfection for 48 h, cells were harvested for subsequent mRNA or protein expression analysis.
Plasmid construction and transfection
The plasmids expressing Flag-tagged Homo sapiens METTL3 or Flag-tagged Homo sapiens HEY2 were synthesized by Shanghai Genechem Co., Ltd. Constructed plasmids were transfected into the endothelial cells according to the manufacturer’s instructions of jetPRIME kit (Polyplus-transfection).
RNA extraction, cDNA synthesis and quantitative real-time PCR (qRT-PCR)
Total RNA was extracted from endothelial cells according to the manufacturer’s instructions for TRIzol reagent (Invitrogen, USA) and then dissolved in RNase-free water. The cDNA was synthesized using PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa Co. Dalian, China). Quantitative Real-Time PCR was performed using the SYBR® Premix Ex Taq™ II (Tli RNaseH Plus) (TaKaRa) on the QuantStudio™ real-time PCR system (Applied Biosystems, Foster City, CA, USA). The specific primers used in this paper were listed in Table S3.
Cells were harvested and lysed in RIPA lysis buffer with protease and phosphatase inhibitors. Protein samples were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to 0.45 μm polyvinylidene difluoride (PVDF) membrane (Merck Millipore). After blocking for 1 h in 5% bovine serum albumin (BSA), the membranes were incubated with the specific primary antibodies (Table S4). And then, the PVDF membrane was incubated with Horseradish peroxidase-conjugated secondary antibodies and immunoreactive bands were visualized with enhanced chemiluminescence reagent (Merck Millipore) according to the manufacturer’s instructions.
Methylated RNA immunoprecipitation (MeRIP)
MeRIP were performing using Magna MeRIP m6A Kit (Millipore, 17–10,499) according to the manufacturer’s instructions. Briefly, 300 μg of total RNA was chemically fragmented into about 100 nucleotides in length by incubation at 94 °C for 1 min in fragmentation buffer, followed by magnetic immunoprecipitation with the monoclonal antibody toward m6A. Methylated RNA was eluted by competition with free m6A, and extracted with RNeasy kit (Qiagen). Both the IP samples and the input samples without immunoprecipitation were used for RNA-seq library generation with NEBNext® Ultra II Directional RNA Library Prep Kit (New England Biolabs, Inc., USA). Thereafter, sequencing was performed by Cloud-Seq Biotech Ltd. Co. (Shanghai, China). The raw data of this study have been deposited in GEO database (GSE142386). For m6A-IP-qPCR, total RNA was chemically fragmented into about 300 nucleotides in length. One-tenth of the fragmented RNA was saved as input control, and further analysed by qPCR along with eluted methylated RNA.
mRNA stability assays
Control and genes silencing endothelial cells were treated with 10 μg per ml actinomycin D (MCE, HY-17559). The total RNAs were then extracted by TRIzol (Invitrogen) at indicated time points and analyzed by qRT-PCR. The turnover rate and half-life of mRNA was estimated according to a previously published paper .
Cell migration experiments were conducted according to the manufacturer’s protocol, 70 μl cell suspension at the 5 × 105 cells/ml concentration was seeding into the ibidi Culture-Insert 2 Well in μ-Dish 35 mm, which would result in a confluent layer within 24 h. After appropriate cell attachment (incubated at 37 °C and 5% CO2 for 24 h), the Culture-Insert 2 Well was gently removed by using sterile tweezer. Then, the used wells were filled with cell free medium, and the original images were obtained by using the Fluorescence Inversion Microscope System. 24 h later, the migrated images were acquired and analyzed with the original images by using the Image J software.
Tube formation assays
Gel matrix was prepared according to the manufacturer’s protocol or reference. The inner well of μ-Slide was filled with 10 μl liquid gel, which could make the gel polymerize under appropriate conditions. Cell suspension at 3 × 105 cells/ml concentration was prepared, and 50 μl of the cell suspension was applied into the upper well. Then the μ–Slide was covered with the supplied lid and incubated at 37 °C and 5% CO2 as usual. About 24 h later, the tube formation images were taken under the Fluorescence Inversion Microscope System and analyzed by using the Image J software.
All experiments were performed at least three independent replicates and statistical analyses were performed with GraphPad Prism6 software. Statistical significance was calculated by two-way ANOVA, one-way ANOVA and/or unpaired Student’s t-test. A p value less than 0.05 was considered statistically significant.
METTL3 is related to the nidus size of cerebral AVMs
A total of 57 patients with cerebral AVMs were included in this study. The baseline characteristics of all patients were summarized in Table S1. Specific anatomical features and Spetzler-Martin grade were determined according to the imaging results of digital subtraction angiography. There were 13 patients (22.8%) with deep vein drainage. The lesions of 25 patients (43.9%) were located in eloquent brain areas. The maximum diameter of lesions ranged from 1.1 cm to 5.9 cm with a mean of 3.3 cm. According to Spetzler-Martin grade, 18 patients (31.6%) were classified as grade 1, 16 patients (28.1%) were classified as grade 2, 15 patients (26.3%) were classified as grade 3, and 8 patients (14.0%) were classified as grade 4. In addition, 21 patients (36.8%) had intracranial hemorrhage.
Previously, we have performed transcriptome sequencing and identified that WTAP was down-regulated in the cerebral AVMs lesions compared to the normal intracranial vascular tissues . Here, we analyzed the differential expression profiles between specific groups of anatomical features, with particular attention to the expression of m6A-related molecules. Finally, we found that the METTL3 expression levels were reduced in the larger nidus compared to smaller nidus (Fig. 1a). Meanwhile, no significant differences were found between different hemorrhagic presentation, different locations, or different Spetzler-Martin grade groups (Fig. 1b-e). However, in the hemorrhagic group, we found that the METTL3 expression levels were not only related to the nidus size of cerebral AVMs, but also related to the age, locations, and Spetzler-Martin grade (Fig. S1). Therefore, METTL3 might be a biomarker related to the nidus size of cerebral AVMs and might be involved in the progression of cerebral AVM.
METTL3 is required for angiogenesis of endothelial cells
In cerebral AVMs, abnormalities in angiogenesis such as tube formation and migration of vascular endothelial cells have been reported. To fully determined the effects of METTL3 on these characteristics, we knocked down METTL3 in endothelial cells by transfecting specific siRNA, and verified knockdown efficiency with qRT-PCR and western blot (Fig. 2a and b). The data of wound healing and migration assays indicated the migratory ability of METTL3 deficient endothelial cells was evidently decreased (Fig. 2c). Moreover, the tube formation results showed that knockdown of METTL3 significantly inhibited formation of capillary-like structures (Fig. 2d). Overall, these results revealed that the expression level of METTL3 could significantly affect angiogenesis of vascular endothelial cells and might therefore affect the size of the lesions.
Notch pathway is activated in METTL3 deficient endothelial cells
In order to identify the target genes involved in angiogenesis, we performed mRNA transcriptome sequencing to interrogate the expression changes between the control and METTL3 deficient endothelial cells. The data showed 706 and 680 genes were significantly up-regulated and down-regulated in METTL3 deficient endothelial cells, respectively (Fig. 3a). Intriguingly, similar to METTL3 deficiency in zebrafish embryos, pathway analysis showed that highly up-regulated genes in METTL3 deficient human endothelial cells also significantly enriched in the Notch signaling pathway (Fig. 3b). Notch is activated by a unique process that includes ligand binding and multistep proteolytic processing, which induces the transcription of Notch target genes, such as HES1, HEY1, and the like bHLH transcription factors [34, 35]. Naturally, we examined the expression of the Notch downstream genes in METTL3 deficient and overexpressing endothelial cells. The results showed that knockdown of METTL3 significantly up-regulated the mRNA and protein levels of the putative arterial endothelial marker HEY2 in endothelial cells (Fig. 3c and d). And, correspondingly, overexpression of METTL3 significantly decreased the expression level of HEY2 (Fig. 3c and d). Therefore, we concluded that the Notch signaling pathway was regulated by METTL3 in endothelial cells. Subsequently, we verified the effects of HEY2 on tube formation and migration in endothelial cells. Same as the phenotype of endothelial cells with knockdown of METTL3, overexpression of HEY2 significantly affected the formation of capillary-like tubes and migration of endothelial cells (Fig. 3e and f).
BMPs belong to the TGF-β superfamily of secreted growth factors and bind receptors to induce phosphorylation and nuclear translocation of SMAD transcription factors (SMAD1/5/9) and regulate vessel growth. However, previous study has proved that Notch directly regulates SMAD6 expression to affect BMP responsiveness of endothelial cells and new vessel branch formation . Therefore, we also detected the expression level of SMAD6 in METTL3 deficient endothelial cells. Indeed, with the activation of the Notch signaling pathway, the expression level of SMAD6 also increased significantly in METTL3 deficient endothelial cells, while the phosphorylation modification of SMAD1/5/9 and SAMD2/3 significantly decreased (Fig. S2). Simultaneously, we found that overexpression of METTL3 suppressed the expression of SMAD6 and up-regulated BMP responsiveness (Fig. S2). Taken together, we confirmed that a decrease in METTL3 expression levels could activate the Notch pathway to affect angiogenesis in endothelial cells.
METTL3 activates Notch pathway via reducing the expression of DTX3L in endothelial cells
METTL3 is well known to be the most important m6A methyltransferase. Knockdown of METTL3 by specific siRNA significantly reduced the m6A modification level of RNAs in endothelial cells (Fig. S3A). To further clarify the downstream targets of METTL3, we performed Methylated RNA Immunoprecipitation Sequencing technology (MeRIP-seq) to profile the changes of transcriptome-wide m6A modification sites between METTL3 deficient and control endothelial cells. Finally, 1576 m6A modified transcripts were identified (Fig. 4a). Among them, 730 genes with significantly decreased m6A modification were identified as the potential targets of METTL3 (Fig. 4a). Pathway analysis revealed that the low levels of m6A modified transcripts were enriched in many biological processes such as MAPK signaling pathway, cAMP signaling pathway, and other immune-related signaling pathways (Fig. S3B). Previous study has showed that Notch1a is a directed target of METTL3 . Surprisingly, although Notch1 were significantly up-regulated and down-regulated in METTL3 deficient and overexpressing endothelial cells, respectively, knockdown of METTL3 did not result in reduction of the m6A modification in Notch1 mRNA (Fig. S3C and D). In addition, knockdown of YTHDF2 did not up-regulated the mRNA and protein levels of Notch1 in human endothelial cells (Fig. S3E and F). A hypothesis raised by these results was that METTL3 could active the Notch pathway by regulating another molecular via m6A modification, rather than directly modulating Notch1 expression.
Then, we further analysed the expression levels of 730 potential targets with reduced m6A modification in METTL3 deficient endothelial cells. Two hundred thirty seven targets were shown to be differentially expressed genes, including 230 significantly down-regulated and 7 up-regulated genes (Fig. 4b). Previous study has demonstrated that DTX1 and DTX3L, functioning as a heterodimeric Notch E3 ligase, concertedly down-regulate Notch activity . Interestingly, MeRIP-seq and RNA-seq results showed that the m6A modification and expression levels of DTX3L mRNA were decreased dramatically in METTL3 deficient endothelial cells (Fig. 4c). Moreover, our results of qRT-PCR and western blot showed that mRNA and protein levels of DTX3L were significantly down-regulated and up-regulated in METTL3 deficient and over-expressed endothelial cells (Fig. 4d and e). m6A-IP-qPCR assays also confirmed that DTX3L mRNA exhibited noticeably decreased m6A modification level in METTL3 deficient endothelial cells (Fig. 4f). Furthermore, knockdown of METTL3 markedly shortened half-live of DTX3L transcript, indicating that METTL3 regulated DTX3L expression through modulating mRNA stability (Fig. 4g). Therefore, we identified that DTX3L was a direct downstream target of METTL3.
To investigate whether DTX3L affected the Notch signaling pathway, we then knocked it down in endothelial cells. Finally, we found that the endogenous NICD level and the expression level of HEY2 were both up-regulated in DTX3L deficient endothelial cells (Fig. 4h and i). In addition, the endogenous NICD level was increased or decreased in METTL3 deficient or overexpressing endothelial cells (Fig. 4e). Therefore, these results suggested that METTL3 could activate the Notch pathway by m6A-dependent regulating its downstream target DTX3L. Tube formation and wound healing assays also showed that knockdown of DTX3L significantly affected angiogenesis of endothelial cells (Fig. 4j and k). In addition to DTX3L, the expression level of DTX1 was also decreased or increased in METTL3 deficient or overexpressing endothelial cells, respectively (Fig. S4). Previously, DTX3L and the human family of DTX proteins (DTX1, DTX2, and DTX3) have been found to function as E3 ligases based on their capacity for self-ubiquitination, while heterodimerization of DTX3L and DTXs enhances this activity . Correspondingly, we speculated that down-regulated METTL3 in cerebral AVMs attenuated DTX3L and DTX1 to synergistically activate Notch signaling pathway, ultimately affecting angiogenesis of endothelial cells.
IGF2BPs enhances DTX3L mRNA stability via an m6A-dependent manner
Previous studies have identified the IGF2BP family members, as m6A readers, promote the stability of their target mRNAs in an m6A-dependent manner [39, 40]. To elucidate whether IGF2BPs involved in the regulation of DTX3L, IGF2BP1, IGF2BP2 and IGF2BP3 were knocked down in endothelial cells, respectively. The RNA stability assays showed that the DTX3L mRNA half-lives were significantly affected by the inhibition of IGF2BP1 and IGF2BP3 (Fig. 5a and b). Owing to knockdown of IGF2BP3 mostly shortened the DTX3L mRNA half-lives, we choose it for further analysis (Fig. 5b). As expected, the mRNA and protein levels of DTX3L decreased dramatically after interfering the expression of IGF2BP3 in endothelial cells (Fig. 5c and e). Therefore, these results suggested that DTX3L was regulated via an m6A-IGF2BPs-dependent manner. Moreover, the endogenous NICD level and the expression level of HEY2 were up-regulated in IGF2BP3 deficient endothelial cells (Fig. 5d and e). Functionally, interference of IGF2BP3 significantly affected the tube formation and migration of endothelial cells (Fig. 5f and g). In summary, we concluded that METTL3 modulated the expression level of DTX3L depending on IGF2BPs.
Notch pathway inhibitor DAPT rescues defective angiogenesis caused by knockdown of METTL3
Considering that the γ-secretase inhibitor DAPT can suppress Notch signaling pathway, we hypothesized that blocking Notch pathway with DAPT could restore the phenotype of METTL3 deficient endothelial cells. To address this hypothesis, the control and METTL3 deficient endothelial cells were treaded with DMSO and DAPT, respectively. As we speculated, compared to DMSO treatment, inhibition of Notch signaling by DAPT significantly reduced the mRNA and protein level of HEY2 in METTL3 deficient endothelial cells (Fig. 6a and b). Additionally, DAPT treatment enhanced tube formation and migration of METTL3 deficient endothelial cells (Fig. 6c and d). Therefore, DAPT may be an important potential therapeutic intervention for cerebral AVMs.
The clinical symptoms such as cerebral hemorrhage, epilepsy, and stroke, caused by cerebral AVM, seriously threaten the life safety of patients. Currently, the pathogenesis of cerebral arteriovenous malformation remains unclear. It is well known epigenetic mechanisms, including DNA, chromatin and histone modifications, play an important role in regulating gene expression . More recently, emerging epitranscriptomics have revealed the extremely important role of RNA modifications in shaping the transcriptomic landscapes [42, 43]. In particular, m6A modifications are dynamic and reversible in mammalian cells, regulating all stages of the RNA life cycle . Abnormal m6A modification levels lead to RNA metabolic disorders and contribute to the development of diseases . Surprisingly, in this study, we found that the METTL3 expression levels were related to the nidus size of cerebral AVMs (Fig. 1a). In addition, depletion of METTL3 resulted in significantly lower m6A modification levels in RNA and affected angiogenesis of endothelial cells (Fig. 2c, d, and S3A). Therefore, m6A modification on RNA, which have been proposed as another layer of epigenetic regulation similar to DNA and histone methylation, may be associated with the development of cerebral AVMs.
Due to the lack of an intervening capillary bed, cerebral arteries and veins are directly connected, leading to hemodynamic changes in cerebral AVMs . In addition, abnormal angiogenesis usually results in structural and functional defects in new blood vessels . Typical angioarchitectural features, hemodynamic alterations and abnormal angiogenesis may result in disordered vascular hyperplasia in cerebral AVMs. Abnormal tangles of dilated vascular structures, including tortuous arteries and dilated veins, form nidus in cerebral AVMs. Therefore, disregulated angiogenesis may affect the size of the lesions in cerebral AVMs. The most widely used Spetzler-Martin grading of cerebral AVMs can be used not only to predict the outcome of microsurgical treatment, but also to predict the SRS outcome [3, 4]. Excellent results were seen with the resection of small AVMs in noneloquent brain with superficial venous drainage. Damage to the regions that control motor, sensory, visual, and language functions (termed eloquent areas of the cortex), as well as impairing deep white-matter pathways and basal ganglia structures are associated with a poor clinical outcome . Larger cerebral AVMs surrounded by a larger volume of parenchyma may result in a higher likelihood of affecting these zones. Taking these into account, METTL3 might be identified as a molecular biomarker linked to treatment outcomes. Moreover, a better understanding of how METTL3 affects lesion size will help elucidate the progression of cerebral AVMs.
Similar to METTL3 deficiency in zebrafish embryos , knockdown of METTL3 significantly activated the Notch pathway in endothelial cells (Fig. 3a and b). However, m6A modifications on Notch1 mRNA were not altered in METTL3 deficient endothelial cells (Fig. S3C). Intriguingly, DTX1 and DTX3L, which synergistically blocked Notch1 activity through functioning as a heterodimeric Notch E3 ligase , were down-regulated in the METTL3 deficient endothelial cells (Fig. 4d, e, and S4). Moreover, DTX3L was identified as the directed downstream gene of METTL3 and was maintained by an m6A-IGF2BPs-dependent mechanism in endothelial cells (Fig. 4 and Fig. 5). Thus, Notch pathway was continuously activated due to the inhibition of ubiquitination degradation, suggesting another layer of dysregulation of Notch signaling that may be involved in cerebral AVMs.
The TGF-β superfamily, consisting of various soluble ligands and their cognate membrane-bound receptors, plays an important role in a diverse range of processes during vessel development . The cross-talk between Notch signaling pathway and BMP pathway has been shown to modulate endothelial cell behavior during vascular development. For example, HEY2 acts as a downstream of Notch to repress BMP2 expression [48,49,50]. However, along with Notch activation, the up-regulated target HEY2 did not inhibited BMP2 in METTL3 deficient endothelial cells (data not shown). Multiple lines of evidences suggest that HEY proteins primarily act as direct transcriptional repressors. Therefore, additional targets of HEY2 might be involved in the regulation of angiogenesis in endothelial cells. Moreover, Notch can directly regulate transcription of SMAD6 to repress BMP pathway responsiveness , which was also confirmed in this study (Fig. S2). Silencing of METTL3 dramatically up-regulated the expression of SMAD6 and inhibited the phosphorylation modification of the SMAD1/5/9 and SMAD2/3 (Fig. S2B). These findings suggested that METTL3 could affect angiogenesis of endothelial cells by negatively regulating TGF-β signaling pathway.
Except for anticonvulsants or acute symptomatic treatment to control associated seizures, so far, there are no drugs available to prevent bleeding or to effectively treat cerebral AVMs. Therefore, it is necessary to reveal the underlying mechanisms of cerebral AVMs and develop drugs accordingly. In this study, we found that down-regulated METTL3 could up-regulate Notch signaling pathway activity in cerebral AVMs. Moreover, Notch signaling pathway had also been shown to be abnormally activated in cerebral AVMs , so blocking Notch signaling may be a potential therapeutic pathway. Here, we demonstrated that DAPT could restore angiogenesis defects caused by reduced METTL3 expression levels and subsequent continuous activation of Notch signaling. Therefore, our results showed the potential therapeutic value of DAPT, and might lay a solid foundation for pharmacological treatment of cerebral AVMs.
Collectively, our work suggests that METTL3 may be a potential molecular biomarker related to nidus size of cerebral AVMs and is likely to be used to predict the outcomes of treatment. Moreover, our results reveal a new paradigm for m6A-Notch axis regulating angiogenesis in cerebral AVM, which will provide suitable pharmacological targets for more effective treatment of patients.
Availability of data and materials
The datasets used and analysed during the current study are available from the corresponding author on reasonable request.
quantitative Real-Time PCR
RNA transcriptome sequencing
Methylated RNA Immunoprecipitation sequencing
Human Umbilical Vein Endothelial Cell
Friedlander RM. Clinical practice. Arteriovenous malformations of the brain. N Engl J Med. 2007;356(26):2704–12. https://doi.org/10.1056/NEJMcp067192.
Colby GP, Coon AL, Huang J, Tamargo RJ. Historical perspective of treatments of cranial arteriovenous malformations and dural arteriovenous fistulas. Neurosurg Clin N Am. 2012;23(1):15–25. https://doi.org/10.1016/j.nec.2011.10.001.
Spetzler RF, Martin NA. A proposed grading system for arteriovenous malformations. J Neurosurg. 1986;65(4):476–83. https://doi.org/10.3171/jns.1986.65.4.0476.
Solomon RA, Connolly ES Jr. Arteriovenous malformations of the brain. N Engl J Med. 2017;376(19):1859–66. https://doi.org/10.1056/NEJMra1607407.
Roundtree IA, Evans ME, Pan T, He C. Dynamic RNA modifications in gene expression regulation. Cell. 2017;169(7):1187–200. https://doi.org/10.1016/j.cell.2017.05.045.
Bokar JA, Shambaugh ME, Polayes D, Matera AG, Rottman FM. Purification and cDNA cloning of the AdoMet-binding subunit of the human mRNA (N6-adenosine)-methyltransferase. Rna. 1997;3(11):1233–47.
Bokar JA, Rath-Shambaugh ME, Ludwiczak R, Narayan P, Rottman F. Characterization and partial purification of mRNA N6-adenosine methyltransferase from HeLa cell nuclei. Internal mRNA methylation requires a multisubunit complex. J Biol Chem. 1994;269(26):17697–704.
Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, et al. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014;10(2):93–5. https://doi.org/10.1038/nchembio.1432.
Ping XL, Sun BF, Wang L, Xiao W, Yang X, Wang WJ, et al. Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase. Cell Res. 2014;24(2):177–89. https://doi.org/10.1038/cr.2014.3.
Schwartz S, Mumbach MR, Jovanovic M, Wang T, Maciag K, Bushkin GG, et al. Perturbation of m6A writers reveals two distinct classes of mRNA methylation at internal and 5′ sites. Cell Rep. 2014;8(1):284–96. https://doi.org/10.1016/j.celrep.2014.05.048.
Xiao W, Adhikari S, Dahal U, Chen YS, Hao YJ, Sun BF, et al. Nuclear m (6) a reader YTHDC1 regulates mRNA splicing. Mol Cell. 2016;61(4):507–19. https://doi.org/10.1016/j.molcel.2016.01.012.
Zhao X, Yang Y, Sun BF, Shi Y, Yang X, Xiao W, et al. FTO-dependent demethylation of N6-methyladenosine regulates mRNA splicing and is required for adipogenesis. Cell Res. 2014;24(12):1403–19. https://doi.org/10.1038/cr.2014.151.
Fustin JM, Doi M, Yamaguchi Y, Hida H, Nishimura S, Yoshida M, et al. RNA-methylation-dependent RNA processing controls the speed of the circadian clock. Cell. 2013;155(4):793–806. https://doi.org/10.1016/j.cell.2013.10.026.
Roundtree IA, Luo GZ, Zhang Z, Wang X, Zhou T, Cui Y, et al. YTHDC1 mediates nuclear export of N(6)-methyladenosine methylated mRNAs. eLife. 2017;6. https://doi.org/10.7554/eLife.31311.
Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, et al. N(6)-methyladenosine modulates messenger RNA translation efficiency. Cell. 2015;161(6):1388–99. https://doi.org/10.1016/j.cell.2015.05.014.
Li A, Chen YS, Ping XL, Yang X, Xiao W, Yang Y, et al. Cytoplasmic m(6) a reader YTHDF3 promotes mRNA translation. Cell Res. 2017;27(3):444–7. https://doi.org/10.1038/cr.2017.10.
Wang X, Lu Z, Gomez A, Hon GC, Yue Y, Han D, et al. N6-methyladenosine-dependent regulation of messenger RNA stability. Nature. 2014;505(7481):117–20. https://doi.org/10.1038/nature12730.
Du H, Zhao Y, He J, Zhang Y, Xi H, Liu M, et al. YTHDF2 destabilizes m(6)A-containing RNA through direct recruitment of the CCR4-NOT deadenylase complex. Nat Commun. 2016;7:12626. https://doi.org/10.1038/ncomms12626.
Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18(1):31–42. https://doi.org/10.1038/nrm.2016.132.
Roignant JY, Soller M. M(6) a in mRNA: an ancient mechanism for fine-tuning gene expression. Trends Genet. 2017;33(6):380–90. https://doi.org/10.1016/j.tig.2017.04.003.
Deng X, Su R, Feng X, Wei M, Chen J. Role of N(6)-methyladenosine modification in cancer. Curr Opin Genet Dev. 2018;48:1–7. https://doi.org/10.1016/j.gde.2017.10.005.
Zhang C, Chen Y, Sun B, Wang L, Yang Y, Ma D, et al. M(6) a modulates haematopoietic stem and progenitor cell specification. Nature. 2017;549(7671):273–6. https://doi.org/10.1038/nature23883.
Gridley T. Notch signaling during vascular development. Proc Natl Acad Sci U S A. 2001;98(10):5377–8. https://doi.org/10.1073/pnas.101138098.
Gridley T. Notch signaling in vascular development and physiology. Development. 2007;134(15):2709–18. https://doi.org/10.1242/dev.004184.
ZhuGe Q, Zhong M, Zheng W, Yang GY, Mao X, Xie L, et al. Notch-1 signalling is activated in brain arteriovenous malformations in humans. Brain. 2009;132(Pt 12):3231–41. https://doi.org/10.1093/brain/awp246.
Li S, Wang R, Wang Y, Li H, Zheng J, Duan R, et al. Receptors of the Notch signaling pathway are associated with hemorrhage of brain arteriovenous malformations. Mol Med Rep. 2014;9(6):2233–8. https://doi.org/10.3892/mmr.2014.2061.
Gale NW, Dominguez MG, Noguera I, Pan L, Hughes V, Valenzuela DM, et al. Haploinsufficiency of delta-like 4 ligand results in embryonic lethality due to major defects in arterial and vascular development. Proc Natl Acad Sci U S A. 2004;101(45):15949–54. https://doi.org/10.1073/pnas.0407290101.
Krebs LT, Shutter JR, Tanigaki K, Honjo T, Stark KL, Gridley T. Haploinsufficient lethality and formation of arteriovenous malformations in Notch pathway mutants. Genes Dev. 2004;18(20):2469–73. https://doi.org/10.1101/gad.1239204.
Krebs LT, Starling C, Chervonsky AV, Gridley T. Notch1 activation in mice causes arteriovenous malformations phenocopied by ephrinB2 and EphB4 mutants. Genesis. 2010;48(3):146–50. https://doi.org/10.1002/dvg.20599.
Murphy PA, Kim TN, Lu G, Bollen AW, Schaffer CB, Wang RA. Notch4 normalization reduces blood vessel size in arteriovenous malformations. Sci Transl Med. 2012;4(117):117ra118. https://doi.org/10.1126/scitranslmed.3002670.
Murphy PA, Lu G, Shiah S, Bollen AW, Wang RA. Endothelial Notch signaling is upregulated in human brain arteriovenous malformations and a mouse model of the disease. Laboratory Invest J Techn Methods Pathol. 2009;89(9):971–82. https://doi.org/10.1038/labinvest.2009.62.
Huo R, Fu W, Li H, Jiao Y, Yan Z, Wang L, et al. RNA sequencing reveals the activation of Wnt signaling in low flow rate brain arteriovenous malformations. J Am Heart Assoc. 2019;8(12):e012746. https://doi.org/10.1161/JAHA.119.012746.
Chen CY, Ezzeddine N, Shyu AB. Messenger RNA half-life measurements in mammalian cells. Methods Enzymol. 2008;448:335–57. https://doi.org/10.1016/S0076-6879(08)02617-7.
Briot A, Iruela-Arispe ML. Blockade of specific NOTCH ligands: a new promising approach in cancer therapy. Cancer Discov. 2015;5(2):112–4. https://doi.org/10.1158/2159-8290.CD-14-1501.
Kopan R, Ilagan MX. The canonical Notch signaling pathway: unfolding the activation mechanism. Cell. 2009;137(2):216–33. https://doi.org/10.1016/j.cell.2009.03.045.
Mouillesseaux KP, Wiley DS, Saunders LM, Wylie LA, Kushner EJ, Chong DC, et al. Notch regulates BMP responsiveness and lateral branching in vessel networks via SMAD6. Nat Commun. 2016;7:13247. https://doi.org/10.1038/ncomms13247.
Choi D, Park E, Jung E, Seong YJ, Yoo J, Lee E, et al. Laminar flow downregulates Notch activity to promote lymphatic sprouting. J Clin Invest. 2017;127(4):1225–40. https://doi.org/10.1172/JCI87442.
Takeyama K, Aguiar RC, Gu L, He C, Freeman GJ, Kutok JL, et al. The BAL-binding protein BBAP and related Deltex family members exhibit ubiquitin-protein isopeptide ligase activity. J Biol Chem. 2003;278(24):21930–7. https://doi.org/10.1074/jbc.%20M301157200.
Huang H, Weng H, Sun W, Qin X, Shi H, Wu H, et al. Recognition of RNA N(6)-methyladenosine by IGF2BP proteins enhances mRNA stability and translation. Nat Cell Biol. 2018;20(3):285–95. https://doi.org/10.1038/s41556-018-0045-z.
Li T, Hu PS, Zuo Z, Lin JF, Li X, Wu QN, et al. METTL3 facilitates tumor progression via an m(6)A-IGF2BP2-dependent mechanism in colorectal carcinoma. Mol Cancer. 2019;18(1):112. https://doi.org/10.1186/s12943-019-1038-7.
Goldberg AD, Allis CD, Bernstein E. Epigenetics: a landscape takes shape. Cell. 2007;128(4):635–8. https://doi.org/10.1016/j.cell.2007.02.006.
Saletore Y, Meyer K, Korlach J, Vilfan ID, Jaffrey S, Mason CE. The birth of the Epitranscriptome: deciphering the function of RNA modifications. Genome Biol. 2012;13(10):175. https://doi.org/10.1186/gb-2012-13-10-175.
Widagdo J, Anggono V. The m6A-epitranscriptomic signature in neurobiology: from neurodevelopment to brain plasticity. J Neurochem. 2018;147(2):137–52. https://doi.org/10.1111/jnc.14481.
Wei W, Ji X, Guo X, Ji S. Regulatory role of N(6) -methyladenosine (m(6) a) methylation in RNA processing and human diseases. J Cell Biochem. 2017;118(9):2534–43. https://doi.org/10.1002/jcb.25967.
MS J, HM D, TF M. Computational network modeling of Intranidal hemodynamic compartmentalization in a theoretical three-dimensional brain arteriovenous malformation. Front Physiol. 2019;10:1250. https://doi.org/10.3389/fphys.2019.01250.
J F. Angiogenesis in cancer, vascular, rheumatoid and other disease. Nat Med. 1995;1(1):27–31. https://doi.org/10.1038/nm0195-27.
Fang J, Hirschi K. Molecular regulation of arteriovenous endothelial cell specification. F1000Research. 2019;8. https://doi.org/10.12688/f1000research.16701.1.
Luna-Zurita L, Prados B, Grego-Bessa J, Luxan G, del Monte G, Benguria A, et al. Integration of a Notch-dependent mesenchymal gene program and Bmp2-driven cell invasiveness regulates murine cardiac valve formation. J Clin Invest. 2010;120(10):3493–507. https://doi.org/10.1172/JCI42666.
Rutenberg JB, Fischer A, Jia H, Gessler M, Zhong TP, Mercola M. Developmental patterning of the cardiac atrioventricular canal by Notch and hairy-related transcription factors. Development. 2006;133(21):4381–90. https://doi.org/10.1242/dev.02607.
Heisig J, Weber D, Englberger E, Winkler A, Kneitz S, Sung WK, et al. Target gene analysis by microarrays and chromatin immunoprecipitation identifies HEY proteins as highly redundant bHLH repressors. PLoS Genet. 2012;8(5):e1002728. https://doi.org/10.1371/journal.pgen.1002728.
We thank Cloud-Seq Biotech Ltd. Co. (Shanghai, China) for the MeRIP-Seq service and the subsequent bioinformatics analysis.
This study was supported by “13th Five-Year Plan” National Science and Technology Supporting Plan (2015BAI12B04), Beijing Science and Technology Supporting PlanD16110000381605 and the Beijing Municipal Administration of Hospitals’ Mission Plan (SML20150501).
Ethics approval and consent to participate
This study was performed in accordance with ethical standards, according to the Declaration of Helsinki, and according to national and international guidelines. Informed consents were obtained from all patients, and this study was approved by the institutional review board of Beijing Tiantan Hospital, Capital Medical University.
Consent for publication
The authors read and approved the manuscript, and declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Baseline characteristics of the samples.
siRNA used in this paper.
Primers used in this paper.
Antibody used in this paper.
The expression levels of METTL3 in different groups of bleeding cerebral AVMs. (A) The expression levels of METTL3 in different size groups of bleeding cerebral AVMs. Small, 0 cm < MD ≤ 3 cm, n = 10; Medium, 3 cm < MD ≤ 6 cm, n = 11. (B) The expression levels of METTL3 in different age groups of bleeding cerebral AVMs. 0 < age ≤ 20, n = 10; 20 < age ≤ 40, n = 7; 40 < age ≤ 60, n = 4. (C) The expression levels of METTL3 in different location groups of bleeding cerebral AVMs. Eloquent, n = 11; Noneloquent, n = 10. (D) The expression levels of METTL3 in different Spetzler-Martin grading groups of bleeding cerebral AVMs. Grade 1, n = 7; Grade 2, n = 3; Grade 3, n = 7; Grade 4, n = 4. P values were calculated using one-way ANOVA and/or Student’s t-test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
SMAD6 is up-regulated in METTL3 deficient endothelial cells. (A) qRT-PCR and (B) western blot analysis of the expression level of indicated genes in METTL3 silencing cells. Data are shown as mean ± SEM of three independent experiments. P values were calculated using Student’s t-test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
Knockdown of METTL3 does not reduce m6A enrichment in Notch1 mRNA. (A) Dot blot showing the m6A modification level in control and METTL3 deficient endothelial cells. (B) GO analysis of the low m6A modification genes in METTL3 knockdown endothelial cells. (C) Integrative Genomics Viewer (IGV) tracks displaying MeRIP-seq read distribution in Notch1 mRNA of control and METTL3 knockdown endothelial cells. (D) qRT-PCR analysis of the expression level of Notch1 after silencing or overexpressing METTL3. (E) qRT-PCR and (F) western blot analysis of the expression level of Notch1 in YTHDF2 silencing cells. Data are shown as mean ± SEM of three independent experiments. P values were calculated using Student’s t-test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
DTX1 is regulated by METTL3 in endothelial cells. (A) qRT-PCR and (B) western blot analysis of the expression level of DTX1 in METTL3 knockdown endothelial cells. Data are shown as mean ± SEM of three independent experiments. P values were calculated using Student’s t-test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
About this article
Cite this article
Wang, Lj., Xue, Y., Huo, R. et al. N6-methyladenosine methyltransferase METTL3 affects the phenotype of cerebral arteriovenous malformation via modulating Notch signaling pathway. J Biomed Sci 27, 62 (2020). https://doi.org/10.1186/s12929-020-00655-w
- Cerebral arteriovenous malformation
- Nidus size
- Notch signaling pathway