Skip to main content

Advertisement

MEF2A alters the proliferation, inflammation-related gene expression profiles and its silencing induces cellular senescence in human coronary endothelial cells

Article metrics

Abstract

Background

Myocyte enhancer factor 2A (MEF2A) plays an important role in cell proliferation, differentiation and survival. Functional deletion or mutation in MEF2A predisposes individuals to cardiovascular disease mainly caused by vascular endothelial dysfunction. However, the effect of the inhibition of MEF2A expression on human coronary artery endothelial cells (HCAECs) is unclear. In this study, expression of MEF2A was inhibited by specific small interference RNA (siRNA), and changes in mRNA profiles in response to MEF2A knockdown were analyzed using an Agilent human mRNA array.

Results

Silencing of MEF2A in HCAECs accelerated cell senescence and suppressed cell proliferation. Microarray analysis identified 962 differentially expressed genes (DEGs) between the MEF2A knockdown group and the negative control group. Annotation clustering analysis showed that the DEGs were preferentially enriched in gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to proliferation, development, survival, and inflammation. Furthermore, 61 of the 578 downregulated DEGs have at least one potential MEF2A binding site in the proximal promoter and were mostly enriched in the GO terms “reproduction” and “cardiovascular.” The protein–protein interaction network analyzed for the downregulated DEGs and the DEGs in the GO terms “cardiovascular” and “aging” revealed that PIK3CG, IL1B, IL8, and PRKCB were included in hot nodes, and the regulation of the longevity-associated gene PIK3CG by MEF2A has been verified at the protein level, suggesting that PIK3CG might play a key role in MEF2A knockdown induced HCAEC senescence.

Conclusions

MEF2A knockdown accelerates HCAEC senescence, and the underlying molecular mechanism may be involved in down-regulation of the genes related with cell proliferation, development, inflammation and survival, in which PIK3CG may play a key role.

Background

MEF2A is one of the members of the MEF2 family that belongs to the MADS-box superfamily and shares homology in a 58-amino acid domain that mediates DNA binding and dimerization [1]. An MEF2A-encoded protein can act as a homodimer or heterodimer and is involved in several cellular processes, including neuronal differentiation, muscle development, cell growth control, and apoptosis [2]. MEF2A has an essential DNA-binding site in the MyoDa gene control region and is a transcription factor that activates many muscle-specific, growth factor-induced, and stress-induced genes [3, 4]. The critical roles of MEF2A in myocardial development have been confirmed in many studies. Naya et al. observed that mice lacking Mef2a die suddenly within the first week of life and exhibit pronounced right ventricle dilation, myofibrillar fragmentation, mitochondrial disorganization, and fetal cardiac gene program activation [5]. Chen et al. found that MEF2A participates in the pathway mediating cardiomyocyte development by AKT2 [6]. The association of genetic variations in the coding sequence of MEF2A with coronary artery disease (CAD) has also been broadly investigated. The deletion of seven amino acids disrupts the nuclear localization of MEF2A and reduces the MEF2A-mediated activation of transcription, and this deletion is co-segregated with premature CAD in a family lineage [7]. Another study has reported that a novel 6 bp deletion in exon 11 of MEF2A is closely associated with premature CAD [8]. Nonsynonymous mutation or variation in the number of CAG or CCG repeats is correlated with susceptibility to CAD in some populations [9, 10]. However, the roles of MEF2A variants in CAD are controversial because not all studies can show the association between these variants and CAD [11,12,13,14,15]. The inflammation and dysfunction of vascular endothelial cells or smooth muscle cells are implicated in the pathogenesis of CAD. MEF2A promotes vascular smooth muscle cell (VSMC) senescence and regulates H2O2-induced VSMC senescence by interacting with miR-143, and MEF2A inhibition induces the phenotypic switching of VSMCs, leading to increased cell proliferation and migration [16, 17]. Kim et al. demonstrated that the activity of MEF2 (MEF2A and MEF2C) is remarkably impaired in pulmonary arterial endothelial cells derived from the subjects with pulmonary arterial hypertension (PAH), and the restoration of the impaired function of endothelial MEF2 rescues PAH [18]. Silencing of MEF2A in apoE−/− mice by using lentiviral shRNA considerably reduces plaque collagen content and fibrous cap thickness and increases plaque area; conversely, silencing of MEF2A has no obvious effect on plaque lipid content [19]. These studies have suggested that MEF2A plays important roles in cardiomyocytes, VSMCs, and pulmonary arterial endothelial cells, but its functions in coronary artery endothelial cells are unclear. In this study, we found that the inhibition of MEF2A by siRNA promoted HCAEC senescence. To explore the transcriptional change in the downstream molecular network in response to MEF2A knockdown, we conducted mRNA microarray analysis and revealed that many specific pathways might be regulated by MEF2A in HCAECs. We also described a downstream molecular network controlled by MEF2A.

Results

Inhibition of MEF2A in HCAECs promoted cellular senescence and suppressed proliferation

In this study, we attempted to understand the effect of inhibiting the MEF2A expression on HCAECs. Four pairs of MEF2A-specific siRNA were designed to suppress the MEF2A expression, and MEF2A-1527 showed the highest inhibitory activity on the expression of MEF2A (Fig. 1a). MEF2A-1527 had a significant inhibitory effect on MEF2A when it was transfected at different concentrations in HCAEC (Fig. 1b), and a final concentration of 40 nM MEF2A-1527 remarkably inhibited MEF2A at both mRNA (Fig. 1c) and protein (Fig. 1d) levels in HCAEC. The results of senescence-associated β-galactosidase (SA-β-Gal) staining and EdU staining showed that MEF2A knockdown accelerates HCAEC senescence (Fig. 2a) and inhibites cell proliferation (Fig. 2b).

Fig. 1
figure1

MEF2A was inhibited by siRNA on mRNA and protein level. a The inhibition efficiency of four pairs of MEF2A specific siRNA. b The inhibition efficiency of different concentration of MEF2A-1527. c, d Variation of the expression of MEF2A at mRNA and protein levels in HCAEC transfected with MEF2A-1527 at a final concentration of 40 nM

Fig. 2
figure2

Alteration in phenotype of HCAEC in response to MEF2A knockdown. NC: negative control; MEF2A-SI: MEF2A silencing. a Senescent cells were detected using SA-β-Gal staining at 96 h after transfection and senescent cells were stained blue. b The proliferative capability of cells were tested by using EdU staining at 96 h after transfection and the cells stained purple have proliferative activity

Identification of DEGs in response to MEF2A knockdown in HCAECs

The MEF2A-specific siRNA group and the negative control group were subjected to genome-wide expression microarray analyses to determine the underlying regulatory network in response to the inhibition of MEF2A. The raw microarray data obtained in this experiment were deposited in Gene Expression Omnibus (GEO) [20] with an accession number of GSE114114. A total of 1363 genes had more than twofold changes in the expression between the MEF2A knockdown group and the NC group, and the difference in the gene expression patterns could be seen in the heat map (Fig. 3a). The DEGs, namely, MEF2A, SERPINP, IL1B, ENPP1, ZAK, PIK3CG, FOS, MMP1, HLA-DQB1, GUCY1A3, MYL9, GJA5, and CSF2, were selected for real-time q-PCR validation. Up-regulation of MMP1 and MYL9 and downregulation of the other genes except GJA5 in the MEF2A-siRNA group were observed in qPCR and microarray. The microarray result showed that GJA5 was significantly down-regulated, but the qPCR result revealed that GJA5 was up-regulated in the MEF2A-siRNA group (Fig. 3b).

Fig. 3
figure3

DEGs between MEF2A silencing group and NC-group. a The heatmap of the DEGs. Three groups of HCAEC transfected with MEF2A-1527 were respectively indicated with 1527-1, 1527-2 and 1527-3. The mixed sample with three groups of HCAEC transfected with negative control siRNA was defined as NC. b, c The selected DEGs were validated with real-time quantitative PCR. The relative expression of each gene was calculated according to 2−((Cte1−Cte0) −(Ctn1−Ctn0)), here, Cte1 and Cte0 respectively indicate the Ct value of the target gene and the internal reference gene (beta-actin) in MEF2A siRNA-groups; Ctn1 and Ctn0 respectively indicate the Ct value of the target gene and beta-actin in the NC groups. The relative expression of all selected genes was normalized to “1” in the NC groups. c These genes have at least one potential MEF2A binding site in its proximal promoter. *p < 0.05, **p < 0.01, ***p < 0.001 are significant

Of the down-regulated DEGs, 61 were predicted with at least one MEF2A binding site in the proximal promoter (500 bp region upstream from the transcription start site) (Additional file 1: Table S1). Database for Annotation, Visualization, and Integrated Discovery (DAVID) analysis of the 61 genes indicated that many genes were significantly enriched in the terms “reproduction” and “cardiovascular” in the GAD_DISEASE_CLASS category. The terms “secreted,” “disulfide bond,” “signal peptide,” and “zinc finger region:B box-type” were the most significantly enriched (Table 1). A few of these DEGs were selected for qPCR validation, and the expression of TRIM7, SGK196, BCL2A1, NR5A2, and ANK3 significantly differed between the two groups (Fig. 3c). These findings suggested for the first time that the expression of these genes may be directly mediated by MEF2A.

Table 1 GO annotation for the down-regulated DEGs with potential MEF2A binding sites in the proximal promoter

DEGs in response to MEF2A knockdown were mostly enriched in the DAVID GO terms related to growth, signal transduction and inflammation

Excluding the DEGs that lacked the ENTREZ GENE ID information, we found that 384 and 578 of 962 DEGs were up-regulated and down-regulated, respectively (Additional file 2: Table S2). The annotation clusters containing GO terms, including glycoprotein, glycosylation site:N-linked (GlcNAc), signal peptide, disulfide bond, secreted, and signal, were significantly enriched in the DEGs (Additional file 3: Table S3). Functional GO term enrichment analysis revealed that preferential changes in the expression of the genes in response to MEF2A knockdown were involved in diverse biological processes, including “female pregnancy,” “cell–cell signaling,” and “metabolic process.” The DEGs were likely implicated in many molecular functions, including “phosphatidylinositol-4,5-bisphosphate 3-kinase activity,” “cytokine activity,” “serine-type endopeptidase inhibitor activity,” and “growth factor activity” (Table 2). The most significantly enriched GO terms in cellular component included extracellular space, extracellular region, and plasma membrane (Table 2). Pathway enrichment analysis showed that most of the DEGs down-regulated by the MEF2A knockdown were enriched in immune-related pathways (Table 3). These GO terms are closely associated with growth, development, cell survival, and inflammation, suggesting the extensive involvement of MEF2A in mediating cell functions.

Table 2 The top 10 enriched biological process, cell component and molecular function GO terms for DEGs, which were sorted respectively by p-value in ascending
Table 3 The top 10 significant KEGG pathway enriched in the down-regulated DEGs

PIK3CG and several chemotactic factors were included in the hot nodes in protein–protein interaction (PPI) networks

The PPI network of the down-regulated DEGs contained 502 nodes and 281 edges, and the p-value of PPI enrichment was < 1.0 × 10−16. Most of the downregulated DEGs formed a complicated interaction network (Fig. 4a). GNG13, GNG2, GABBR2, CXCR7, CXCL3, IL8, CNR1, IL1B, PRKCB, PMCH, P2RY14, and PIK3CG were presented at the hot nodes and may play important roles in the biological process (Fig. 4a). The majority of the genes in the network were significantly enriched (FDR < 0.01) in the categories of extracellular space and extracellular region in the GO term of cellular component, and many of these genes were involved in signal transduction (Table 4). The KEGG pathways “amoebiasis, Staphylococcus aureus infections, influenza A, rheumatoid arthritis, and olfactory transduction” were significantly enriched (FDR < 0.05) for the DEGs in the PPI network (Table 4). However, most of the up-regulated genes are independent of each other and only a small number of them can form an interaction network (Fig. 4b).

Fig. 4
figure4

The functional association protein networks were constructed with STRING for the DEGs. Settings for the interaction search as follows: minimum required interaction score is ≥ 0.7 (high confidence); active interaction sources are Textmining, Experiments, Databases, Co-expression, Neighborhood, Gene Fusion and Co-occurrence; the max number of interactors to show is query proteins only for the first shell and none for the second shell; hide the disconnected nodes in the network. The pies indicated nodes and the full line indicated the edges of the PPI network. a Network for down-regulated genes. b Network for UP-regulated genes

Table 4 GO annotation for the down-regulated DEGs in the PPI network

In order to discover the genes possibly involved in related diseases, the DEGs were annotated to the category “GAD_DISEASE_CLASS” and the terms “metabolic, chemdependency, cardiovascular, hematological, psych, neurological, renal, vision, pharmacogenomic, reproduction, normal variation, immune, aging, and developmental” were preferentially enriched (Fig. 5a). Among the DEGs in the GO term “aging,” 23 were up-regulated and 33 were down-regulated. The PPI network of the DEGs in the GO term “aging” showed that 22 DEGs, including PIK3CG, TXK, HDAC9, PPARG, IL1B, IL8, and PCK1, were in a correlative interaction network (Fig. 5b). A PPI network analysis performed for 244 DEGs in the GO term “cardiovascular” showed that these genes had more interactions among themselves than the expected observation (PPI enrichment p < 1 × 10−16), suggesting that the genes were biologically connected as a group (Fig. 5c). PIK3CG, IL8, IL1B, and CSF were included in the hot nodes with multiple interactions, implying their major role in the GO terms “aging” and “cardiovascular.” Further validation experiments showed that the mRNA and protein levels of PIK3CG were significantly down-regulated when MEF2A was inhibited (Fig. 6a), while the mRNA and protein levels of PIK3CG were significantly up-regulated when overexpressing MEF2A (Fig. 6b).

Fig. 5
figure5

Enrichment of the GO terms in GAD_DISEASE_CLASS for the DEGs and PPI network for certain GO term ‘aging’ and ‘cardiovascular’ related DEGs. In the PPI network, each node represents all the proteins produced by a single, protein-coding gene locus. Edges represent protein–protein associations, and associations are meant to be specific and meaningful. a The Go terms significantly enriched for the DEGs. b Interaction network for the DEGs enriched in the GO term ‘aging’. A “d” indicates a gene with down-regulatory signal on the microarray. A “u” indicates a gene with up-regulatory signal on the microarray. c Interaction network for the DEGs enriched in the GO term ‘cardiovascular’

Fig. 6
figure6

The expression of PIK3CG varies with the level of MEF2A. a, b Changes in the expression of PIK3CG at the mRNA and protein levels when MEF2A is inhibited; NC: transfected with negative control siRNA; MEF2A-SI: transfected with MEF2A specific siRNA (MEF2A-1527). c, d Changes in the expression of PIK3CG at mRNA and protein levels when MEF2A was over-expressed

Discussion

In the present study, inhibition of MEF2A by siRNA promoted HCAEC senescence, suggesting that it has an anti-aging functional role. This protective role appears to be different from its role in VSMC, which exhibits a detrimental aging-promoting role [16, 17]. However, an in vivo research has shown that MEF2A may have an anti-atherosclerotic function in ApoE−/− mice [19], and functional deletion in MEF2A gene is co-segregated with premature CAD in several family lineages [7, 8], suggesting the protective role of MEF2A in vessels. The above findings imply that a decrease in MEF2A expression in vascular endothelial cells may result in vessel dysfunction. Therefore, it is important to understand the profile and the interaction network of the DEGs affected by knockdown of MEF2A in HCAEC.

In our study, the DEGs in response to MEF2A knockdown were significantly enriched in the GO terms associated with growth, proliferation, development, survival, and inflammation. The down-regulated genes in MEF2A-knocked down HCAECs were involved in a complicated regulatory network. Among the hot nodes, PIK3CG, PRKCB, and IL8 have been extensively investigated. PIK3CG plays important roles in immunity, hypertension, longevity, and circulation [21, 22]. p110γ, the encoding protein of PIK3CG, is principally involved in the control of leukocyte chemotaxis and inflammatory reactions in the immune system, affects multiple aspects of adrenergic receptor signaling and cardiac function in cardiomyocytes [21], and participates in angiotensin II induced calcium influx and governing the vascular contractile response [22]. Thus, p110γ has been considered as a promising target for drug design in the treatment of heart failure and hypertension, whereas the function of p110γ in vascular endothelial cells is unclear. In our study, the GO term analysis based on the GAD_disease indicated that PIK3CG participated in the aging pathway and interacted with many inflammatory factors in the PPI network. In practice, mass evidences have indicated that PI3K pathway is involved in cellular senescence and proliferation [23,24,25,26]. In the present study, multiple experiments showed that the expression of PIK3CG varied with the expression level of MEF2A, suggesting that the expression of PIK3CG is regulated by MEF2A. The expression of PIK3CG was significantly down-regulated in MEF2A-knocked down HCAECs, suggesting that PIK3CG may be involved in the molecular mechanism of MEF2A silencing induced senescence. In addition, many inflammatory factors were overtly down-regulated simultaneously with PIK3CG, indicating that the down-regulation of inflammatory factors in response to MEF2A silencing might be mediated by PIK3CG.

PRKCB might promote the removal of LDLs from the blood by increasing the expression of LDLR and SR-A and reducing foam cell formation, which might mitigate atherogenesis. However, the characteristics of PRKCB that plays a pro-inflammatory role in endothelial cells and promotes VSMC proliferation and migration are associated with pro-atherosclerosis as observed in PRKCB knockout mice [27].

IL8 is another extensively studied gene that serves as a risk factor of cardiovascular disease [28]. The PPI network shows a large number of interactions between IL8 and other genes, suggesting that MEf2A-mediated inflammatory pathway might greatly depend on the participation of IL-8. Suzuki et al. [29] revealed that MEF2A may mediate vascular inflammation by stimulating MCP-1 expression in VSMCs and macrophage infiltration. MEF2A and MEF2D may function as activators or repressors of gene transcription in human macrophages [30]. In our study, many genes involved in macrophage infiltration were downregulated in MEF2A-knocked down HCAECs, suggesting the role of MEF2A in regulating inflammation. Furthermore, 8 of the top 10 KEGG pathways enriched by the down-regulated DEGs were involved in immune-related diseases, indicating that MEF2A participated in mediating the immune system.

Guanine nucleotide binding protein subunit gamma 2 and 13 (GNG2 and GNG13) play key roles in signal transduction pathways. In our study, they were significantly down-regulated in MEF2A-knocked down HCAECs. GNG13 is colocalized with alpha-gustducin in taste receptor cells, and antibodies against Ggamma13 block the denatonium-induced increase in inositol trisphosphate (IP3) in taste tissue [31]. In our study, most olfactory receptors were distinctively down-regulated in MEF2A-knocked down HCAECs, and almost all of them interacted with GNG13 in the PPI network. These findings suggested that the expression of olfactory receptors might be regulated by MEF2A through GNG13 and PI3K.

DAVID GO functional annotation enrichment based on GAD revealed that the cardiovascular system was one of the most prominently enriched terms of the DEGs in response to MEF2A knockdown. More than 200 DEGs were related to the cardiovascular system and exhibited complicated interactions with one another. Many hot nodes in the PPI network are inflammatory factors, including IL8, IL1B, CXCL3, CXCL12, and HLA genes, which are associated with cardiovascular diseases [28, 32,33,34]. MEF2A functional deletion carriers were subjected to premature CAD [7, 8], strongly supporting the important roles of MEF2A in providing protection against cardiovascular disease. Our findings show that numerous genes in inflammatory pathways were down-regulated as MEF2A was silenced, indicating that MEF2A may be an important regulator to mediate the expression of inflammatory factors, and that this role of MEF2A appears to be inconsistent with its protective effect against CAD. Alternative splicing of the pre-mRNA of MEF2A resulted in different isoforms which have distinct activity and are differentially distributed in various tissues [35]. Aberrant alternative splicing isoforms were observed to be associated with myotonic dystrophy and neuromuscular disorders [36]. More than 20 MEF2A isoforms have been submitted to NCBI database so far. Yet the functional diversity of the isoforms is unclear. Silencing of MEF2A by siRNA may alter the proportion of different splicing isoforms compared to in the normal physiological status of the cells. Therefore, when we attempt to explore the underlying molecular mechanism that MEF2A participate in a certain biological or pathological process, which MEF2A alternative splicing isoforms exert major effect on the process should be taken into account. We found that the functional role of MEF2A in vascular endothelial cells is inconsistent with its roles in VSMC reported by other studies, possibly due to differential splicing of MEF2A in the two kinds of cells.

Among the down-regulated DEGs, 61 genes were found with at least one potential MEF2A binding site in the proximal promoter. If extended the search region in the promoters of the DEGs, there would be more genes found with MEF2A binding sites. Therefore, many DEGs may be directly regulated by MEF2A via binding to their cis-elements in the promoters. TRIM7, SGK196, BCL2A1, NR5A2, ANK3 and PYHIM9 have MEF2A binding site in their proximal promoter, which were down-regulated in MEF2A-knocked down HCAEC via verifying by qPCR. These genes are respectively involved in metabolism, development, differentiation, apoptosis and survival [37, 38].

Conclusion

MEF2A knockdown in HCAECs promotes cell senescence and remarkably alters mRNA profile, primarily related to signal transduction, proliferation, development, inflammation, and other processes, and may lead to cell dysfunction by promoting cellular senescence and inflammation or by inhibiting cell proliferation. PIK3CG has been shown to be associated with longevity, and its expression varies with the expression level of MEF2A in HCAEC, suggesting that PIK3CG may play an important role in cell senescence induced by MEF2A silencing.

Materials and methods

Cell culture

The primary HCAECs were purchased from Cell Applications and cultured in EGM-2 medium (Lonza, USA) supplemented with growth factors. The cells were further grown at 37 °C in a humidified atmosphere with 5% (v/v) CO2 and passaged at 1:3 when they reached approximately 85–95% confluency.

Construction of MEF2A overexpression plasmid

The complete MEF2A coding region fragment was amplified using high fidelity PCR polymerase and cDNA template derived from HCAEC, and the amplified fragment was ligated into the mammalian expression vector pc3-gab at the sites EcoRI and BglII. Then, transformation and plasmid extraction were performed, and the inserted encoding sequence was verified by Sanger sequencing to be same to the MEF2A alternative splice variant XM_005254915.2_V_X15.

Transfection

siRNA was diluted to a 20 μM solution with RNAse-free and sterile deionized water and stored below − 20 °C. HCAECs were seeded in a plate at a density of 7500 cells/cm2 18–24 h before transfection and cultured at 37 °C in a humidified atmosphere with 5% (v/v) CO2. Lipofectamine® RNAiMAX and Lipofectamine® 3000 reagent were respectively used to guide siRNA and plasmids into the cells, and the procedure was performed in accordance with the manufacturer’s instructions. The next treatments were performed at 72 h post-transfection.

Detection of cell senescence and proliferation

The cells were subjected to SA-β-Gal staining by using a β-galactosidase staining kit (Cell Signaling Technology, USA) in accordance with the manufacturer’s instructions. The stained cells were examined with a high-resolution microscope at 400× magnification, and blue stain corresponded to senescent cells. Cell proliferation was detected with an EdU staining kit (Ribobio, China) in accordance with the manufacturer’s instructions. The stained cells were imaged with a fluorescence inversion microscope system at 400× magnification, and purple stain indicated proliferating cells. The senescence rate and the proliferation rate were respectively determined by calculating the proportion of the senescent cells and the proliferating cells in the total number of cells in 15 visions per well. The results were expressed as the mean of triplicates ± standard error of the mean (SEM).

RNA isolation

Total RNA of cells was isolated with Trizol reagent (Invitrogen, Carlsbad, USA) in accordance with the manufacturer’s protocol. RNA concentration was determined by measuring the absorbance at 260 nm with a NanoVue spectrophotometer (GE Healthcare, UK). RNA integrity was assessed through electrophoresis on agarose gel.

Microarray analysis

Three independent groups of HCAECs transfected with MEF2A-specific siRNA (MEF2A-Homo-1527) were separately analyzed with an Agilent human mRNA array (AHMA), and three independent groups of HCAECs transfected with negative control siRNA were collected together into a tube. RNA was extracted for AHMA analysis. Three MEF2A-specific siRNA groups were indicated as 1527-1, 1527-2, and 1527-3, and the mixed control group was designated as NC. The AHMA was designed with eight identical arrays/slide (8 × 60 K format), and each array contained probes examining about 27,958 Entrez gene RNAs. The whole microarray analyses, including sample preparation, labeling, hybridization, scanning, and treatment of original data, were processed by CapitalBio Corporation (Beijing, China). DEGs were analyzed by using the DAVID v6.8, and a protein–protein interaction (PPI) network was searched using the STRING tool [39]. The DEGs containing a potential MEF2A binding site in their promoter were predicted with LASAGNA-Search 2.0 [40].

Real-time fluorescence quantitative PCR (q-PCR)

Total RNA was reverse transcribed by using a ReverTra Ace qPCR RT kit (Toyobo, Japan), and q-PCR was performed with SYBR® Green mix (Takara, China) in accordance with the manufacturer’s instructions. Alteration in the gene expression relative to the control group was analyzed with a comparative CT method by using beta-actin for endogenous control and NC cell for calibration. The primers used in this study are shown in Table 5.

Table 5 The primer sequences for real time quantitative PCR

Statistical analysis

Data were expressed as mean ± SEM with more than three independent experiments. The statistical significance of the comparison between groups was examined with Student’s t test, and p < 0.05 indicates significant differences.

Abbreviations

MEF2A:

myocyte enhancer factor 2A

HCAEC:

human coronary artery endothelial cells

DEG:

differentially expressed genes

GO:

gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

CAD:

coronary artery disease

SA-β-Gal:

senescence-associated β-galactosidase

VSMC:

vascular smooth muscle cell

PAH:

pulmonary arterial hypertension

LDL:

low density lipoprotein cholesterol

LDLR:

low density lipoprotein cholesterol receptor

PPI:

protein–protein interaction

siRNA:

small interference RNA

AHMA:

Agilent human mRNA array

DAVID:

Database for Annotation, Visualization, and Integrated Discovery

References

  1. 1.

    Santelli E, Richmond TJ. Crystal structure of MEF2A core bound to DNA at 1.5 A resolution. J Mol Biol. 2000;297(2):437–49.

  2. 2.

    Ornatsky OI, McDermott JC. MEF2 protein expression, DNA binding specificity and complex composition, and transcriptional activity in muscle and non-muscle cells. J Biol Chem. 1996;271(40):24927–33.

  3. 3.

    Kaushal S, Schneider JW, Nadal-Ginard B, Mahdavi V. Activation of the myogenic lineage by MEF2A, a factor that induces and cooperates with MyoD. Science. 1994;266(5188):1236–40.

  4. 4.

    Erck C, Seidl K. MyoD and MEF2A mediate activation and repression of the p75NGFR gene during muscle development. Biochem Biophys Res Commun. 1998;245(3):871–7.

  5. 5.

    Naya FJ, Black BL, Wu H, Bassel-Duby R, Richardson JA, Hill JA, Olson EN. Mitochondrial deficiency and cardiac sudden death in mice lacking the MEF2A transcription factor. Nat Med. 2002;8(11):1303–9.

  6. 6.

    Chen D, Chen F, Xu Y, Zhang Y, Li Z, Zhang H, Pan T, Su Y, Wan M, Wang X, et al. AKT2 deficiency induces retardation of myocyte development through EndoG-MEF2A signaling in mouse heart. Biochem Biophys Res Commun. 2017;493(4):1410–7.

  7. 7.

    Wang L, Fan C, Topol SE, Topol EJ, Wang Q. Mutation of MEF2A in an inherited disorder with features of coronary artery disease. Science. 2003;302(5650):1578–81.

  8. 8.

    Xu DL, Tian HL, Cai WL, Zheng J, Gao M, Zhang MX, Zheng ZT, Lu QH. Novel 6-bp deletion in MEF2A linked to premature coronary artery disease in a large Chinese family. Mol Med Rep. 2016;14(1):649–54.

  9. 9.

    Maddhuri S, Gudala S, Lakkaraju C, Malempati AR, Nallari NP, Mundluru HP. Association of Mef2a gene polymorphisms in early onset of coronary artery disease of south Indian cohort. Meta Gene. 2018;15:16–22.

  10. 10.

    Foroughmand AM, Shahbazi Z, Galehdari H, Purmahdi Borujeni M, Dinarvand P, Golabgirkhademi K. Association of MEF2A gene polymorphisms with coronary artery disease. Iran Red Crescent Med J. 2014;16(8):e13533.

  11. 11.

    Altshuler D, Hirschhorn JN. MEF2A sequence variants and coronary artery disease: a change of heart? J Clin Invest. 2005;115(4):831–3.

  12. 12.

    Bhagavatula MR, Fan C, Shen GQ, Cassano J, Plow EF, Topol EJ, Wang Q. Transcription factor MEF2A mutations in patients with coronary artery disease. Hum Mol Genet. 2004;13(24):3181–8.

  13. 13.

    Liu Y, Niu W, Wu Z, Su X, Chen Q, Lu L, Jin W. Variants in exon 11 of MEF2A gene and coronary artery disease: evidence from a case-control study, systematic review, and meta-analysis. PLoS ONE. 2012;7(2):e31406.

  14. 14.

    Dai DP, Zhou XY, Xiao Y, Xu F, Sun FC, Ji FS, Zhang ZX, Hu JH, Guo J, Zheng JD, et al. Structural changes in exon 11 of MEF2A are not related to sporadic coronary artery disease in Han Chinese population. Eur J Clin Invest. 2010;40(8):669–77.

  15. 15.

    Jiang Y, Liu HB. Myocyte enhancer factor-2A gene mutation and coronary artery disease. Chin Med J. 2015;128(19):2688–91.

  16. 16.

    Zhao W, Zheng XL, Peng DQ, Zhao SP. Myocyte enhancer factor 2A regulates hydrogen peroxide-induced senescence of vascular smooth muscle cells via microRNA-143. J Cell Physiol. 2015;230(9):2202–11.

  17. 17.

    Zhao W, Zhao SP, Peng DQ. The effects of myocyte enhancer factor 2A gene on the proliferation, migration and phenotype of vascular smooth muscle cells. Cell Biochem Funct. 2012;30(2):108–13.

  18. 18.

    Kim J, Hwangbo C, Hu X, Kang Y, Papangeli I, Mehrotra D, Park H, Ju H, McLean DL, Comhair SA, et al. Restoration of impaired endothelial myocyte enhancer factor 2 function rescues pulmonary arterial hypertension. Circulation. 2015;131(2):190–9.

  19. 19.

    Zhou WP, Zhang H, Zhao YX, Liu GQ, Zhang JY. RNA interference of myocyte enhancer factor 2A accelerates atherosclerosis in apolipoprotein E-deficient mice. PLoS ONE. 2015;10(3):e0121823.

  20. 20.

    Gene Expression Omnibus. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE114114. Accessed 12 Mar 2019.

  21. 21.

    Perino A, Ghigo A, Hirsch E. Leukocyte and cardiac phosphoinositide 3-kinase gamma activity in pressure overload-induced cardiac failure. Trends Cardiovasc Med. 2010;20(8):273–6.

  22. 22.

    Carnevale D, Lembo G. PI3Kgamma in hypertension: a novel therapeutic target controlling vascular myogenic tone and target organ damage. Cardiovasc Res. 2012;95(4):403–8.

  23. 23.

    Xu Y, Li N, Xiang R, Sun P. Emerging roles of the p38 MAPK and PI3K/AKT/mTOR pathways in oncogene-induced senescence. Trends Biochem Sci. 2014;39(6):268–76.

  24. 24.

    Liu S, Liu S, Wang X, Zhou J, Cao Y, Wang F, Duan E. The PI3K-Akt pathway inhibits senescence and promotes self-renewal of human skin-derived precursors in vitro. Aging Cell. 2011;10(4):661–74.

  25. 25.

    Vredeveld LC, Possik PA, Smit MA, Meissl K, Michaloglou C, Horlings HM, Ajouaou A, Kortman PC, Dankort D, McMahon M, et al. Abrogation of BRAFV600E-induced senescence by PI3K pathway activation contributes to melanomagenesis. Genes Dev. 2012;26(10):1055–69.

  26. 26.

    Li Y, Shan Z, Liu C, Yang D, Wu J, Men C, Xu Y. MicroRNA-294 promotes cellular proliferation and motility through the PI3K/AKT and JAK/STAT pathways by upregulation of NRAS in bladder cancer. Biochemistry. 2017;82(4):474–82.

  27. 27.

    Fan HC, Fernandez-Hernando C, Lai JH. Protein kinase C isoforms in atherosclerosis: pro- or anti-inflammatory? Biochem Pharmacol. 2014;88(2):139–49.

  28. 28.

    Apostolakis S, Vogiatzi K, Amanatidou V, Spandidos DA. Interleukin 8 and cardiovascular disease. Cardiovasc Res. 2009;84(3):353–60.

  29. 29.

    Suzuki E, Satonaka H, Nishimatsu H, Oba S, Takeda R, Omata M, Fujita T, Nagai R, Hirata Y. Myocyte enhancer factor 2 mediates vascular inflammation via the p38-dependent pathway. Circ Res. 2004;95(1):42–9.

  30. 30.

    Aude-Garcia C, Collin-Faure V, Bausinger H, Hanau D, Rabilloud T, Lemercier C. Dual roles for MEF2A and MEF2D during human macrophage terminal differentiation and c-Jun expression. Biochem J. 2010;430(2):237–44.

  31. 31.

    Huang L, Shanker YG, Dubauskaite J, Zheng JZ, Yan W, Rosenzweig S, Spielman AI, Max M, Margolskee RF. Ggamma13 colocalizes with gustducin in taste receptor cells and mediates IP3 responses to bitter denatonium. Nat Neurosci. 1999;2(12):1055–62.

  32. 32.

    Aukrust P, Yndestad A, Smith C, Ueland T, Gullestad L, Damas JK. Chemokines in cardiovascular risk prediction. Thromb Haemost. 2007;97(5):748–54.

  33. 33.

    Liu B, Xiong L, Tian C, Zhou Q, Zhong Y, Li A, Liu SM. HLA-DRB1*12:02:01 plays a protective role against coronary artery disease in women of southern Han Chinese descent. Hum Immunol. 2012;73(1):122–6.

  34. 34.

    Xiong Y, Wang L, Mo P, Huang G, Li A, Chai R, Lin X, Zhong Y, Liu B, Ou W, et al. Association between HLA-DQB1 alleles and susceptibility to coronary artery disease in Southern Han Chinese. Hum Immunol. 2017;78(9):540–6.

  35. 35.

    Zhu B, Ramachandran B, Gulick T. Alternative pre-mRNA splicing governs expression of a conserved acidic transactivation domain in myocyte enhancer factor 2 factors of striated muscle and brain. J Biol Chem. 2005;280(31):28749–60.

  36. 36.

    Bachinski LL, Sirito M, Bohme M, Baggerly KA, Udd B, Krahe R. Altered MEF2 isoforms in myotonic dystrophy and other neuromuscular disorders. Muscle Nerve. 2010;42(6):856–63.

  37. 37.

    Di Costanzo S, Balasubramanian A, Pond HL, Rozkalne A, Pantaleoni C, Saredi S, Gupta VA, Sunu CM, Yu TW, Kang PB, et al. POMK mutations disrupt muscle development leading to a spectrum of neuromuscular presentations. Hum Mol Genet. 2014;23(21):5781–92.

  38. 38.

    Medina-Echeverz J, Haile LA, Zhao F, Gamrekelashvili J, Ma C, Metais JY, Dunbar CE, Kapoor V, Manns MP, Korangy F, et al. IFN-gamma regulates survival and function of tumor-induced CD11b+ Gr-1high myeloid derived suppressor cells by modulating the anti-apoptotic molecule Bcl2a1. Eur J Immunol. 2014;44(8):2457–67.

  39. 39.

    STRING. http://string-db.org. Accessed 12 Mar 2019.

  40. 40.

    Lee C, Huang CH. LASAGNA-search: an integrated web tool for transcription factor binding site search and visualization. Biotechniques. 2013;54(3):141–53.

Download references

Authors’ contributions

YX and BL conceptualized the study. LW, WJ and LP performed the experiment. YX, WL, AL and YZ helped analyze and interpreted the data. WO, BL and SL prepared the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Dunzhou Mo for giving many suggestions on the microarray data analysis.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and analysed during the current study are available in the Gene Expression Omnibus (GEO) repository with an accession number of GSE114114 [20].

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

This study was supported by the innovation project from the Department of Education of Guangdong Province (2014KTSCX098), Guangdong Provincial Science and Technology plan project (2016A020215142) and the General Programs of the National Natural Science Foundation of China (81370463 to Liu B and 81570259 to Liu SM).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Correspondence to Benrong Liu or Shi-ming Liu.

Additional files

12867_2019_125_MOESM1_ESM.xlsx

Additional file 1: Table S1. The information of the predicted MEF2A binding sites on the proximal promoter of the DEGs.

12867_2019_125_MOESM2_ESM.xlsx

Additional file 2: Table S2. The genes with more than twofold change in expression in the siRNA group compared to the negative control. Up, up-regulated in MEF2A specific siRNA group; down, down-regulated in MEF2A specific siRNA group.

12867_2019_125_MOESM3_ESM.xlsx

Additional file 3: Table S3. Annotation for the differentially expressed genes.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xiong, Y., Wang, L., Jiang, W. et al. MEF2A alters the proliferation, inflammation-related gene expression profiles and its silencing induces cellular senescence in human coronary endothelial cells. BMC Molecular Biol 20, 8 (2019) doi:10.1186/s12867-019-0125-z

Download citation

Keywords

  • Myocyte enhancer factor 2A
  • Human coronary artery endothelial cells
  • Vascular endothelial dysfunction
  • Senescence
  • Cardiovascular disease
  • PIK3CG