Overviews of the MeRIP-sequencing and RNA-sequencing profiles in MCAO/R Rats
The schematic protocol of this part was shown (Fig. 1a). The 2,3,5-triphenyltetrazolium chloride (TTC) staining was performd to roughly estimate infarct volume (Fig. 1b). All of the MeRIP-seq data were deposited into NCBI-SRA database under the project number of PRJNA73198. According to the MeRIP-peak region situations, there were 23,121 and 23,000 m6A peaks identified in the sham and MCAO/R tissues respectively of which 3,319 peaks within the same mRNAs overlapped in this two groups (Fig. 1c). To clarify the preferential distribution of m6A modification in transcripts, we investigated the metagene profiles of all identified m6A peaks in the entire transcriptome which showed that m6A peaks were mainly enriched in the 3-UTR and Exon regions in both groups (Fig. 1d). To further elucidate the m6A methylation differences between this two groups, a differential m6A methylation level analysis was carried out. It was demonstrated that when compared with the sham group, a total of 9,192 differential m6A-methylated peaks were identified among which 2,686 hypermethylated and 7,506 hypomethylated m6A sites were discovered respectively (fold changes ≥ 1.5, p < 0.05 vs Sham) in the MCAO/R group (Fig. 1e). In addition, transcriptome profiles of altered mRNAs were detectd through RNA-seq. 1,497 differentially expressed genes (DEGs) were uncovered (MCAO/R vs Sham, fold changes ≥ 2 and p < 0.05) including 989 upregulated and 508 downregulated DEGs (Fig. 1e, right)
TFAP2B standed out based on combinant interactive bioinformatics analyses of the MeRIP and RNA-sequencing data
In order to excavate key DEGs which mRNA levels changed due to the m6A modification as much as possible, combinant interactive bioinformatics analyses of the MeRIP-Seq and RNA-Seq data were carried out (Fig. 2a). It was identified 83 DEGs with hypermethylated m6A peaks among which 74 DEGs were significantly upregulated (hyper-up) and 9 were downregulated (hyper-down). On the other hand, 359 DEGs with hypomethylated m6A peaks, consisting of 78 upregulated (hypo-up) and 281 downregulated (hypo-down) (Fig. 2b). The results suggested that DEGs with altered m6A-methylated peaks were mainly enriched in the “hyper-up” or “hypo-down” which is consistent with the previous literature [4]. Then, clustering was made again according to the following parameters: the significance of DEG (p < 0.05, |log FC| > 1.0) and with significant alterations in m6A abundance (p < 0.01, |log FC| > 1.0), based the results, top20 DEGs were picked out again (Fig. 2c). Among these DEGs, transcription factor TFAP2B attracted our interest which has never been studied in ischemic cerebrosis. In order to make the visual contrast more intuitive and clear, the peaks of IP and the corresponding input samples were drawn overlapped of TFAP2B (Fig. 2d). Ordinary QPCR and m6A-specific qPCR were implemented to verify the mRNA levels (Fig. 2e). and m6A abundance (Fig. 2f) differences between sham and MCAO/R tissues respectivly, and the result showed both the relative mRNA levels and the m6A-modified abundance of TFAP2B were significantly lower in MCAO/R when compared with sham, which were consistent with the RNA-seq and MeRIP-Seq result
TFAP2B was inhibited in primary neurons under the OGD/R context
OGD/R induced neurons are widely considered as classic cell models of CIRI. In order to further explore the expression of TFAP2B in OGD/R induced neurons, OGD/R neuron models were established, in which OGD lasted for 1,2,4 hrs respectively followed by 12 hr of reoxygenation (Fig. 3a). The primary neurons of newborn rats were isolated and cultured (Fig. 3b) . The mRNA levels of HIF-α were detected to be a preliminary indicator of the successful construction of the model. It was shown that when OGD lasted for 2 hr, the relative mRNA level of HIF-α began to increase. when OGD lasted for up to 4 hr, the relative mRNA level of HIF-α showed a very significant upregulation (vs NC, Fig. 3c). Next, the mRNA levels of TFAP2B were detected accordingly. It was shown that when OGD lasted for both 2 hrs and 4 hrs, the relative mRNA levels of TFAP2B reduced significantly compared with NC group (Fig. 3d).Then western blot was carried out and the result showed TFAP2B expression was significantly reduced under OGD/R induction compared to the NC (Fig. 3e). Besides, m6A-specific QPCR was conducted and the results showed the m6A abundance of TFAP2B in OGD/R neurons was obviously reduced when compared with NC group which was in line with that in MCAO/R tissues (Fig. 3f).
TFAP2B effectively improved the survival of OGD/R induced primary neurons
The experimental design of this part is shown in Fig. 4a. TFAP2B-overexpressed lentivirus was infected into primary rat neuronal cells to construct the gain-of-function neuronal cell model which was verified at both mRNA and protein levels (Fig. 4b and c). Annexin-V/FITC assay showed that compared with the negative control group (Lv-NC-Ov), the apoptosis of TFAP2B-overexpressed neurons decreased from 14.5–6.8% after OGD/R treatment (Fig. 4d). In order to confirm this result again, tunel staining was performed and the result demonstrated that after OGD/R treatment, the proportion of neurons with red apoptotic signal in Lv-TFAP2B-Ov group was significantly lower than that in negative control (Lv-NC-Ov) group (Fig. 4e). As neurons hardly proliferate, a Edu assay was conducted to reflec cell vitality. It could be observed that after OGD/R, the number of edu-positive neurons showing red signal increased obviously in Lv-TFAP2B-Ov neurons when compared with the Lv-NC-Ov neurons (Fig. 4f). The above results demonstrated that TFAP2B could effectively improve the survival of primary neurons under OGD.
RNA-sequencing of primary neurons with TFAP2B overexpression.
As a transcription factor, exploring its transcription targets is crucial for revealing the potential molecular mechanisms of its biological function. To explore the mechanism through which TFAP2B improve the survival of OGD/R induced primary neurons, we performed RNA sequencing of cell extracts from TFAP2B-overexpressing primary neurons and controlled primary neurons, treated with OGD/R, to seek candidate target genes (RNA-sequencing data were deposited into NCBI-SRA database under the project number of PRJNA73199). It was found that 1624 mRNA were up-regulated and 575 mRNA were down regulated significantly in TFAP2B overexpressed neurons (Lv-TFAP2B-Ov vs Lv-NC-Ov, p < 0.05, |log FC| > 1.0) (Fig. 5a). Further, t-test was carried out, and DEGs with significant differences were screened out (p < 0.05); Next, DEGs with FPKM ≥ 10.0 were screened out and were sorted according to their coefficient of variation (CV), and the top 25 genes with the upper and lower levels were selected separately and the cluster heat maps were generated (Fig. 5b). Among these genes, BNIP3 is an important mediator of mitophagy [23,26] and became the focus of our study, although we cannot rule out the potential function of the other genes. Then western blot was used to empirically confirm the trends of BNIP3 and validated it successfully, when TFAP2B was overexpressed, both the mRNA level and the protein expression of BNIP3 were increased correspondingly (Fig. 5c and d). Subsequently, Jaspar (https://jaspar.genereg.net/) was used to predict the potential binding segments of TFAP2B in the BNIP3 promoter region. It was found that there are two peridcted putative sites with relative profile score threshold 90% in the + strand (-1452-1442 nt;. -1267-1257 nt) (Fig. 5e). Furthermore, luciferase reporter assays were used to validate whether TFAP2B activating transcription of BNIP3 promoter. BNIP3 promoter sequences containing wild-type (WT-P1 and P2) or mutant-type (MUT-P1 and P2) were inserted in the luciferase reporter plasmid promoter. It was found that luciferase activity from the WT-P1, but not the MUT-P1, was markedly increased through TFAP2B overexpressing which demonstrated that the P1 segments(-1452-1442 nt) is one of the binding site of TFAP2B in the BNIP3 promoter (Fig. 5f and g).
TFAP2B protects against the OGD/R induced apoptosis through mediating mitophagy.
BNIP3 serves as mitochondria associated protein that can promote mitophagy. Due to the transcriptional activation of BNIP3 by TFAP2B, it was speculated that TFAP2B alleviated neuronal apoptosis through promoting mitophagy in damaged mitochondria under the OGD/R context and experiments were performed (Fig. 6a). To verify the hypothesis, firstly, western blots were implemented to evaluate autophagy markers (LC3II/I, P62) and mitochondrial mass markers ( TOMM20). The results demonstrated that the LC3II/I was enhanced and the autophagy substrate P62 decreased in TFAP2B-Overexpressed neuron cells (vs.Lv-NC-Ov) under OGD/R treatment suggesting autophagy promotion. Correspondingly, there was a certain degree of Tomm20 recovery in TFAP2B-Overexpressed neuron cells (vs.Lv-NC-Ov) under OGD/R treatment implying the recovery of mitochondrial mass to some extent (Fig. 6b and c). Besides, the expression of apoptotic executing protein c-caspase3 was also significantly reduced in TFAP2B-Overexpressed neuron cells (vs.Lv-NC-Ov) under OGD/R treatment (vs.Lv-NC-Ov) (Fig. 6b and c). The involvement of TFAP2B was further evaluated by mitophagosome formation through the co-localization of LC3B and Mito-Tracker Green. As indicated by colocalization of RFP-LC3B and MitoTracker, mitophagy was increased a little in OGD/R-treated cells (Lv-NC-Ov + OGD/R vs. Lv-NC-Ov). As expected, TFAP2B overexpression significantly promoted the mitophagy formation as indicated by the increased co-localization of LC3B and MitoTracker Green. These results suggested TFAP2B mediates OGD/R-induced mitophagosome formation in neurons (Fig. 6d). To explore whether the regulation of TFAP2B on mitophagy was partially mediated by BNIP3, BNIP3 was silenced (Lv-BNIP3-si) in neurons overexpressing TFAP2B using lentivirus infection again. Subsequently, the mitophagy of neurons under the OGD/R was evaluated again. It was found that the mitophagy formation was decreased obviously as indicated by the reduced co-localization of LC3B and MitoTracker Green when BNIP3 was knocked-down in TFAP2B overexpressed neurons (Fig. 6e). On the other hand, QPCR assays showed that, under OGD/R, the relative mRNA levels of mitochondrial mass markers Tommy20 and CoxⅣ which had originally recovered with TFAP2B overexpression, decreased again with BNIP3 knockdown, while the apoptotic protein BAX, which had significantly decreased with TFAPB overexpression, increased again with BNIP3 KD (Lv-BNIP3-si vs. Lv-NC-si) (Fig. 6f).
TFAP2B overexpression alleviated the MCAO/R damage to a certain extent in vivo
In vivo, Nissl staining and NeuN + tunel IF double staining were performed (Fig. 7a). The Nissl body is the main site for protein synthesis in neurons. When neurons are damaged, Nissl bodies in the neurons significantly decrease or even disappear, ultimately leading to neuronal necrosis. The Nissl dying showed that when compared with AAv-NC-Ov rats, the neurons in the sham groups both (AAv-TFAP2B-Ov and AAv-NC-Ov were plump in shape, with obvious nuclei and abundant Nissl bodies, suggesting the normal viability of nerurons. However, after MCAO/R, without TFAP2B overexpression, neurons were wrinkled and necrotic but the with TFAP2B overexpression, neuronal injury was alleviated obviously when compared with the AAv-NC-Ov group (Fig. 7b). In order to show the apoptosis of nerve cells more intuitively, immunofluorescence mulitiple staining was performed using Cy3 labeled neuron-specific nuclear protein (NeuN, red) ,tunel (green) and DAPI (blue). By global immunofluorescence staining, it can be seen that no obvious TUNEL positive areas of nerve cells (green fluorescence indicates apoptosis) were found of the sham group.After MCAO/R, the right half of the brains in both groups showed obvious neuronal cell apoptosis (green). However, in the TFAP2B overexpression group, the apoptosis areas was significantly narrowed compared with the AAv-NC-Ov group (Fig. 7c). At the same time, the local results also showed that after MCAO/R, the AAv-TFAP2B-Ov group showed a lower neuronal apoptosis rate when compared with the AAv-NC-Ov group (Fig. 7d). These results proved that TFAP2B overexpression alleviated the MCAO/R damage to a certain extent in vivo