Abstract

Background. The transcriptional networks of Cyr61 and its function in cell injury are poorly understood. The present study depicted the lncRNA and mRNA profiles and the involvement in angiotensin II-induced injury after Cyr61 knockdown mediated by CRISPR/Cas9 in HEK293T cells. Methods. HEK293T cells were cultured, and Cyr61 knockdown was achieved by transfection of the CRISPR/Cas9 KO plasmid. lncRNA and mRNA microarrays were used to identify differentially expressed genes (DEGs). Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to determine biofunctions and signaling pathways. RT-PCR was used to validate the microarray results. Cells were divided into four groups: control, Cyr61 knockdown, angiotensin II (Ang II) without Cyr61 knockdown, and Ang II with Cyr61 knockdown. CCK8, western blotting, and flow cytometry analysis were carried out to dissect cellular function. Results. A total of 23184 lncRNAs and 28264 mRNAs were normalized. 26 lncRNAs and 212 mRNAs were upregulated, and 74 lncRNAs and 233 mRNAs were downregulated after Cyr61 knockdown. Analysis of cellular components, molecular functions, biological processes, and regulatory pathways associated with the differentially expressed mRNAs revealed downstream mechanisms of the Cyr61 gene. The differentially expressed genes were affected for small cell lung cancer, axon guidance, Fc gamma R-mediated phagocytosis, MAPK signaling pathway, focal adhesion, insulin resistance, and metabolic pathways. In addition, Cyr61 expression was increased in accordance with induction of cell cycle arrest and apoptosis and inhibition of cell proliferation induced by Ang II. Knockdown of Cyr61 in HEK293T cells promoted cell cycle procession, decreased apoptosis, and promoted cell proliferation. Conclusions. The Cyr61 gene is involved in Ang II-induced injury in HEK293T cells. Functional mechanisms of the differentially expressed lncRNAs and mRNAs as well as identification of metabolic pathways will provide new therapeutic targets for Cyr61-realated diseases.

1. Introduction

Cysteine-rich protein 61 (Cyr61), also known as CCN1, was identified as secreted protein members of the CCN family. Cyr61 which is produced and secreted by endothelial cells, fibroblasts, and smooth muscle cells is a component of the extracellular matrix [1]. Cyr61 has the ability to regulate multiple signaling pathways and play an important role in different biological processes. In vitro, CCN proteins regulate the functions of cells, such as adhesion, migration, proliferation, differentiation, and survival, as well as induce cell apoptosis and extracellular matrix remodeling of different cell types. In vivo, it also plays a vital role in vascular and bone development and angiogenesis [2]. Furthermore, Cyr61 has highlighted the pivotal role this molecule can play in regulating the immune-surveillance process and it has emerged as an important partner when targeting components of the infectious or chronic inflammatory disease processes such as atherosclerosis or rheumatoid arthritis [35]. Additionally, several lines of evidence have shown that overexpression of Cyr61 is involved in the cancer process such as osteosarcoma [6], ovarian tumor [7], and breast adenocarcinoma [8]. Interestingly, Cyr61 was also found to be downregulated in prostate cancer and non-small cell lung carcinoma [9, 10]. This can be ascribed to its direct binding to different integrins in different cell types and a variety of environments [11]. Hence, Cyr61 may play a pivotal role and serve as a potential prognostic indicator or therapeutic target in a variety of diseases.

The mammalian genome is pervasively transcribed, producing a large number of noncoding RNAs (ncRNAs), including long noncoding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides [12, 13]. Studies have found that only one-fifth of transcription across the human genome is associated with protein-coding genes, indicating that there are at least four times more long noncoding than coding RNA sequences [14]. lncRNAs play an important role in the regulation of thousands of protein-encoding genes via various mechanisms, although these are not translated into proteins [15]. Transcriptome profiling has been carried out in various genes and conditions in order to characterize their functions and underlying mechanisms [1619]. Unfortunately, studies of Cyr61 functional recovery have not yet been able to highlight its transcriptional networks. No study so far has reported using the gene chip array to find out about the differentially expressed lncRNAs and mRNAs associated to Cyr61.

As one of the earliest and most extensively studied hormonal systems, the renin-angiotensin system (RAS) which is involved in many physiological and pathological processes is an atypical hormonal system in several ways [20, 21], and the new components and functions of the RAS have still not been unraveled [22]. The RAS has powerful effects including the control of blood pressure and sodium homeostasis as well as fluid through integrated actions in the kidney, the cardiovascular system, and the central nervous system. Along with its impact on blood pressure, the RAS also affects others pathological processes including inflammation and immune responses [23]. Research on the RAS has contributed significantly to advances in the understanding and treatment of cardiovascular diseases [24]. As one of the major bioactive products of the RAS, angiotensin II (Ang II) has a well-known role in cardiovascular regulation and is a key regulator of renal inflammation and fibrosis as well as blood pressure and renal hemodynamics [2527]. Several pathological effects are caused by the Ang II-induced TGF-β pathway which is a potent inflammatory as well as fibrotic and apoptotic cytokine [28, 29]. In addition, Ang II participates in cell proliferation, matrix degradation, inflammation, and apoptosis by activating a multiplicity of signaling pathways [30, 31].

In this study, the lncRNA and mRNA profiles were determined after Cyr61 knockdown mediated by CRISPR/Cas9 in HEK293T cells. The signal transduction and metabolic pathways associated with these differentially expressed lncRNA and mRNA were analyzed. The effect of Cry61 on Ang II-induced cell injury in HEK293T cells was also analyzed. The results provide the functional mechanisms of Cry61 as well as potential new targets for Cyr61-associated diseases.

2. Materials and Methods

2.1. Cell Lines and Cell Culture

The HEK293T cell line (which was purchased from the Molecular Microbiology and Immunology, Keck School of Medicine, University of Southern California) was cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, Su Zhou, USA) in the presence of 10% fetal bovine serum (FBS, Lanzhou Minhai Bio-Engineering, Gansu, China). All of the cells were maintained at 37°C in a humidified atmosphere with 5% CO2.

2.2. Knockdown of Cyr61 Gene Expression

The guide RNA sequences used are shown in Table 1. Cyr61 CRISPR/Cas9 KO plasmid, homology-directed repair (HDR) transfection plasmid, and UltraCruz Transfection Reagent were purchased from Santa Cruz Biotechnology (USA), and gene knockdowns were performed according to the manufacturer’s protocol. The HEK293T cells were seeded in 6-well plates at cells/3 mL per well, and the serum-free medium was changed after 70% confluence was achieved by routine culture. 1 μg of each of Cyr61 CRISPR/Cas9 KO plasmid and HDR transfection plasmid and 10 μL of UltraCruz reagents were added to 300 μL of serum-free and antibiotic-free medium and in 6-well plates for incubation at room temperature for 10 min, followed by addition of antibiotic-free 10% FBS-containing media to a final volume of 2 mL. Wild-type HEK293T cells were treated with different concentrations of puromycin, and the results showed that the lowest inhibitory concentration to the HEK293T cells after 7 d of treatment was 8 μg/mL and this concentration was subsequently used for screening. At 48 h posttransfection, the culture medium was replaced with 8 μg/mL puromycin-containing DMEM medium for screening the cells.

2.3. Classification of Cells

The HEK293T cells were divided into four groups: (1) control group, nontransfected cells, (2) Cyr61-downregulated group, the third generation incubated after transfection, (3) Ang II group, nontransfected cells treated with Ang II (10-7 mol/L), and (4) Cyr61 downregulated + Ang II group, transfected cells treated with Ang II (10-7 mol/L). Cyr61 expression in the control group and Cyr61-downregulated group were identified by western blotting, and these were subjected to microarray analysis. Cell proliferation, apoptosis, and cell cycle assays were performed for all the groups.

2.4. RNA Extraction and Quality Control

Total RNA was extracted from each sample by soaking the samples in TRIzol Reagent (Invitrogen, Grand Island, NY, USA) in accordance with the manufacturer’s instructions. RNA quantity and quality were measured using a NanoDrop ND-1000, and RNA integrity was assessed by standard denaturing agarose gel electrophoresis. RNA quality test results are shown as Supplemental Material (Table 2). For spectrophotometer analysis, only samples with OD260/OD280 and OD260/OD230 were acceptable for further analysis.

2.5. Identification of Differentially Expressed lncRNAs and mRNAs

The mRNA-lncRNA Affymetrix microarrays (raw data) were used in this study. Limma (linear models for microarray data) package in R was used to identify differentially expressed lncRNAs and mRNAs between the 2 groups through use of the -test. Fold and were regarded as the criteria for differential expression.

2.6. Gene Ontology (GO) and Pathway Enrichment Analysis

GO analysis is frequently used in functional enrichment studies of large-scale genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to analyze the biological pathways, involving the differentially expressed mRNAs. In the current study, DAVID (the database for Annotation Visualization and Integrated Discovery) software was used to investigate the functional enrichment condition for the up- and downregulated differentially expressed mRNAs. was selected as the threshold.

2.7. Quantitative Real-Time PCR Validation

To validate the microarray data, we selected the top 3 most significant upregulated (FRMD1, SERPINF1, and FEZ1) and the top 3 most significant downregulated (PRR21, REG3G, and ACAT2) mRNAs as well as the top 3 most significant upregulated (RP11-659F24.1, RP11-966I7.4, and LAMB2P1) and the top 3 most significant downregulated (ANKRD30BL, CH17-360D5.2, and SOX2-OT) lncRNAs from mRNAs and lncRNAs that were expressed aberrantly. Briefly, 2 μg of total RNA from each sample was used for the synthesis of first-strand cDNA using a FastKing RT Kit (Tiangen, Beijing, China) according to the manufacturer’s protocol. Following first-strand cDNA synthesis, a PCR reaction was carried out in a 20 μL reaction volume containing ×1 SuperReal PreMix Plus, 0.6 μL of each specific forward and reverse primer, and 2 μL of cDNA template. PCR assays were set at an initial denaturation step at 95°C for 3 min, followed by 40 cycles of 95°C for 10 s and 60°C for 30 s in a LightCycler® 96 system (Roche, Switzerland). Relative changes in gene expression were calculated using the 2-∆∆Ct method with glyceraldehyde 3-phosphate dehydrogenase (GAPDH) as a reference gene as described previously [32, 33]. Each qPCR was carried out in triplicate, with each triplicate data point repeated 3 times. For the triplicates, only samples that differed by <0.3 Ct values were used in the final calculations. Less than about 5% of samples fell outside this range. In addition, standard curves were generated for each gene using a 50-fold dilution range. The data shown represented the means of three experiments. The primers for each gene are listed in Table 3.

2.8. Western Blotting Analysis

The HEK293T cells were harvested and washed twice with PBS, lysed in ice-cold radio-immunoprecipitation assay buffer (RIPA, Beyotime Biotechnology, Shanghai, China) with a freshly added 0.01% protease inhibitor cocktail (CWBiotech, Beijing, China), and incubated on ice for 30 minutes. Protein concentration was measured by bicinchoninic acid (BCA) protein assay kit (Beyotime Biotechnology, Shanghai, China). A total of 50 μg of protein was subjected to electrophoresis using SDS-PAGE and transferred by electrophoresis onto a nitrocellulose membrane (GE Healthcare, Freiburg, Germany). Blots were visualized using enhanced chemiluminescence (ECL, Millipore, Billerica, MA, USA) after antibody binding. The antibody against Cyr61 was purchased from Santa Cruz Biotechnology (USA); antibodies against GAPDH were purchased from Cell Signaling Technology Biotech (Danvers, MA, USA). The bands were quantified by their densitometry with ImageJ software (NIH, USA).

2.9. Cell Proliferation Assay

Cell proliferation was performed using Cell Counting Kit-8 (CCK-8, Solarbio Life Sciences, Beijing, China) assay. The HEK293T cells were seeded in 96-well plates at a density of cells per well and cultured at 37°C in 5% CO2 for 12 hours and treated with Ang II (10-7 mol/L) [34]. After incubation for 0, 12, 24, and 48 h, 10 μL CCK-8 solution was added to each well and incubated for another 1 to 4 hours. The absorbance was measured in the multimode microplate reader TriStar LB 941 (Berthold Technologies, Bad Wildbad, Germany) at 450 nm. Each experiment was performed in triplicate.

2.10. Cell Apoptosis Assay

Apoptosis was determined by using a flow cytometer (FACSCanto II, BD Biosciences, San Jose, CA, USA), and an Annexin V fluorescein isothiocyanate (FITC)/propidium iodide (PI) double-stain assay was performed in accordance with the manufacturer’s protocol (FITC Annexin V/PI, BD Biosciences, San Diego, USA). After incubation for 48 hours, each supernatant was collected in a centrifuge tube and each group of cells was washed three times with phosphate-buffered saline (PBS), trypsinized, centrifuged ( at room temperature) for 5 min, adjusted to cells/mL, and suspended in binding buffer containing Annexin V-FITC and PI. After incubation for a further 15 min at room temperature in the dark, the fluorescent intensity was measured using a flow cytometer. Each experiment was performed in triplicate.

2.11. Cell Cycle Assay

For cell cycle analysis, each group of cells was harvested after 48 h. Cell cycle analysis was conducted with the cell cycle detection kit (KeyGen, China) according to the manufacturer’s instructions. Briefly, cells were harvested by regular trypsin digestion and rinsed with PBS. After fixation in 70% ethanol overnight at -20C, cells were rinsed again in PBS and stained in a staining solution containing PI and RNase A. Flow cytometry (Canto II, BD Biosciences, USA) was employed to analyze the cellular DNA content. Each treatment was performed in triplicate.

2.12. Statistical Analysis

The results are presented as the mean using SPSS23. One-way analysis of variance (ANOVA) was used among diverse groups, and independent samples were analyzed by Student’s -test when appropriate. A value of was regarded as statistically significant.

To determine if DEGs overlapped with or were significantly enriched with a specific gene set, 28264 expressed genes in our samples were used as background for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses assessed by DAVID.

3. Results

3.1. The Expression Levels of Cyr61 Protein in HEK293T Cells

To examine the biological function of Cyr61 in HEK293T cells, the expression of Cyr61 was knocked down in HEK293T cells by the CRISPR/Cas9 KO plasmid, and three separate CRISPR gRNA sequences targeted to the Cyr61 gene were designed. The HEK293T cells were transfected with CRISPR/Cas9 vectors containing each of the three target sequences. After screening by puromycin, one homozygous was found by western blotting analysis. As shown in Figure 1, the Cyr61 protein expression levels were significantly decreased in HEK293T cells from to and from to (Cyr61/GAPDH, respectively) after knockdown from control levels.

3.2. Differential Gene Expression Analysis between Experimental and Control Groups

100 differentially expressed lncRNAs with a fold change cutoff of 1.0 (26 upregulated and 74 downregulated, ) were identified from the lncRNA microarrays of the Cyr61 experimental and control groups. Volcano plots were constructed to identify differences in the lncRNAs. The top 10 most significant up- and downregulated lncRNAs are shown in Table 4. RP11-659F24.1 and ANKRD30BL were the most up- and downregulated lncRNA transcript with a fold change of 2.87 and 2.13, respectively. The volcano plots display the expression profiles for all detected lncRNAs (Figure 2(a)).

446 differentially expressed mRNAs were also identified by using mRNA microarrays with a fold change cutoff of 1.0 (212 upregulated and 223 downregulated, ) in Cyr61 knocked down HEK293T cells when compared with controls. The top 10 most significant up- and downregulated mRNAs are shown in Table 5. FRMD1and PRR21 were the most up- and downregulated mRNA transcripts with the fold change of 2.7 and 2.02, respectively. The volcano plots display the expression profiles for all detected mRNAs (Figure 2(b)). The heat map is shown in Figure 2(c).

3.3. Gene Ontology Analysis

The 446 differentially expressed mRNAs underwent Gene Ontology (GO) analyses by using the DAVID software. GO analysis of the differentially expressed genes was performed to determine the gene product attributes in biological processes, cellular components, and molecular functions. denoted the significance of GO term enrichment in the differentially expressed genes, and the lower the value, the more significant the GO term. The top 10 biological processes with the highest values are listed in Table 6 and Figure 3(a), and the most significant process is related to “regulation of cell migration.” The top 10 cellular components with the highest values are shown in Table 7 and Figure 3(b), and the most significant one is the “cytosol.” The top 10 molecular functions terms with the highest values are displayed in Table 8 and Figure 3(c), and the most significant one is “poly (A) RNA binding.”

3.4. KEGG Pathway Analysis

Pathway analysis is a functional analysis that maps genes to the KEGG pathways. The target genes for differentially expressed mRNAs were mapped onto signaling pathways for KEGG pathway analysis (Table 9 and Figure 3(d)). The differentially expressed genes were affected for small cell lung cancer, axon guidance, Fc gamma R-mediated phagocytosis, MAPK signaling pathway, focal adhesion, insulin resistance, and metabolic pathways. The results revealed that the differentially expressed gene annotation for metabolic pathway and MAPK signaling pathway are the top 2 signaling pathways, while the small lung cancer is the most significant.

3.5. Validation of lncRNA Expression by RT-qPCR

To validate the up- and downregulation of Cyr61 changes of mRNA and lncRNA expression detected by microarray, six mRNAs (FRMD1, SERPINF1, FEZ1, PRR21, REG3G, and ACAT2) and six lncRNAs (RP11-659F24.1, RP11-966I7.4, LAMB2P1, ANKRD30BL, CH17-360D5.2, and SOX2-OT) were selected. Their expression was examined with RT-qPCR. We observed a good agreement across the two methods. The results of the RT-qPCR were similar to those obtained from the microarray. In both microarrays and RT-qPCR, the expression levels of the mRNAs FRMD1, SERPINF1, and FEZ1 and lncRNAs RP11-659F24.1, RP11-966I7.4, and LAMB2P1 were upregulated, and the expression levels of the mRNAs PRR21, REG3G, and ACAT2 and lncRNAs ANKRD30BL, CH17-360D5.2, and SOX2-OT were downregulated after knocking down the Cyr61 gene. This validation indicated good reproducibility and reliability of the observed changes in expression detected by mRNA and lncRNA microarrays (Figure 4).

3.6. Downregulation Cyr61 Promotes Proliferation of HEK293T Cells

To investigate the role of Cyr61 knockdown on the growth of HEK293T cell lines, we treated the HEK293T cells with Ang II at various time points (0, 12, 24, and 48 h) after knocking down the Cyr61 gene and then performed the CCK-8 assay to examine the proliferation of HEK293T cells. Knockdown of the Cyr61 gene increased the cell proliferation of HEK293T cells significantly, and Ang II inhibited proliferation when compared with control groups. The proliferation rate of Ang II with Cyr61 knockdown groups was significantly increased compared to cells treated with Ang II without Cyr61 knockdown. There was no significant difference between controls and Ang II with Cyr61 knockdown (Figure 5). Taken together, these data suggest a pro-proliferative role of downregulation of the Cyr61 gene and an antiproliferative role of Ang II in HEK293T cells.

3.7. Downregulation of Cyr61 Inhibits Apoptosis of HEK293T Cells

In order to investigate the role of Cyr61 on the apoptosis of HEK293T cells, we treated the HEK293T cells with Ang II after knocking down the Cyr61 gene for 48 h, and then cells were stained with Annexin V-FITC and PI, followed by flow cytometry analysis. Under both conditions, the apoptosis rate was decreased significantly in Cyr61 knockdown groups () and significant cell apoptosis was observed in the Ang II groups () compared with controls (). The percentage of Annexin V-FITC for Ang II with Cyr61 knockdown groups were significantly decreased () compared with Ang II without Cyr61 knockdown (). However, there was no significant difference between controls () and Ang II with Cyr61 knockdown (; Figure 6). Taken together, these data suggest a proliferative role after downregulation of the Cyr61 gene and an antiproliferative role of Ang II in HEK293T cells.

3.8. Downregulation of Cyr61 Promotes Cell Cycle Progression of HEK293T Cells

Having found the interaction between Cyr61 and Ang II on the apoptosis and proliferation in HEK293T cells, the interaction between these on the cell cycle was investigated by treating HEK293T cells with Ang II for 48 h after knocking down the Cyr61 gene. Cells were stained with PI and collected to be examined on a flow cytometer. An increase in the percentage of cells in the G2/M phase () and a marked decrease in the percentage of cells in the S phase () in Cyr61 knockdown groups compared with the other groups were seen. In addition, a decrease in the percentage of cells in the G2/M phase () and a marked increase in the percentage of cells in the S phase () in Ang II groups were also seen. There was no significant difference between controls (G2/M phase: , S phase: ) and Ang II with Cyr61 knockdown groups (G2/M phase: , S phase: ; Figure 7). Downregulation of Cyr61 promoted the cell cycle progression, and treatment with Ang II inhibited the cell cycle progression.

3.9. Ang II Promotes Cyr61 Expression in HEK293 T Cells

To investigate whether Ang II induces Cyr61 expression in 293 cells, we incubated HEK293T cells for 48 hours in DMEM with or without Ang II. Cyr61 transcript levels as determined by western blots were elevated (, Cyr61/GAPDH) in HEK293T cells exposed to Ang II compared with controls (, Cyr61/GAPDH). Furthermore, the Cyr61 protein transcript levels of HEK293T cells which had the Cyr61 gene knocked down was elevated significantly (, Cyr61/GAPDH) after treating with Ang II compared to cells not treated with Ang II (, Cyr61/GAPDH; Figure 8).

4. Discussion

To date, the CRISPR/Cas9 system is known as a molecular scissor and is widely used in various studies, including cancer research, drug discovery, treatment of mental disease, and applications in plant husbandry [35]. A limitation of this method is that the Cas9 construct is either permanently integrated or must be later removed with a subsequent reagent delivery and/or clonal selection step to achieve editing without scarring [36, 37]. In this study, we knocked down Cyr61 via CRISPR/Cas9 plasmid cotransfection with the Cyr61 HDR plasmid. Cyr61 HDR utilizes template DNA (usually the homologous chromosome) to repair DNA in a precise or “error-free” manner [38]. Cyr61 HDR plasmid cotransfected with Cyr61 CRISPR/Cas9 KO plasmid and was designed for repair of the site-specific Cas9-induced DNA cleavage within the Cyr61 gene. During repair, the Cyr61 HDR plasmid incorporates a puromycin resistance gene to enable selection of stable knockout cells and an RFP (red fluorescent protein) gene to visually confirm that the transfection has occurred. Using nonintegrated plasmid vectors that express a puromycin N-acetyltransferase (PAC) gene, whose expression and translation is linked to that of Cas9, we selected for cells based on their early expression levels of Cas9 protein. Edited cells isolated using this method did not contain any detectable off-target mutations and displayed expected functional phenotypes after directed differentiation [38].

Cyr61 is a secreted, extracellular matrix- (ECM-) associated signaling protein of the CCN family [11]. Cyr61 is capable of regulating a broad range of cellular activities, including cell adhesion, migration, proliferation, differentiation, apoptosis, and senescence through interaction with cell surface integrin receptors and heparan sulfate proteoglycans. During embryonic development, Cyr61 is critical for cardiac septal morphogenesis, blood vessel formation in placenta, and vascular integrity. In adulthood, Cyr61 plays important roles in inflammation and tissue repair and is associated with diseases related to chronic inflammation, including rheumatoid arthritis, atherosclerosis, diabetes-related nephropathy and retinopathy, and many different forms of cancers.

lncRNAs are defined as transcripts longer than 200 nucleotides [39]. It is estimated that at least 90% of RNAs transcribed by the human genome are lncRNAs [40] and was previously considered as “junk gene” or transcriptional “noise” [41]. Although the function of lncRNA remains unclear, in recent years, numerous studies have found that lncRNAs have important biological functions; regulate gene expression via epigenetic regulation, transcriptional regulation, and transcriptional regulation [42]; and participate in embryonic development, cell proliferation, differentiation, apoptosis, and several other biological processes. lncRNAs also play important roles in modulating the innate and adaptive immune responses and immune cell development [43, 44]. Moreover, emerging evidence suggests that lncRNAs play a decisive role in such autoimmune diseases as systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), type 1 diabetes mellitus (T1DM), multiple sclerosis (MS) [45], autism spectrum disorders [46], tumor biology [47, 48], and thoracic and abdominal aortic aneurysms [49]. Future research that deepens our understanding of the biological roles and properties of lncRNAs will enable researchers to harness their potential as valuable diagnostic and therapeutic targets for human diseases.

Unfortunately, until now, Cyr61 functional recovery highlighting its transcriptional networks has not been reported and lncRNAs may play an important regulatory role and possibly mediate Cyr61-involved diseases. To explore the molecular mechanisms involved in Cyr61-induced changes of cell functions, in this study, high-throughput microarray technology designed for genome-wide identification was performed to detect the profiles of aberrantly expressed lncRNAs and mRNAs in downregulated Cyr61 HEK293 cells. A total of 23184 lncRNAs and 28264 mRNAs were normalized. 26 lncRNAs and 212 mRNAs were upregulated, and 74 lncRNAs and 233 mRNAs were downregulated after Cyr61 knockdown. Analysis of the cellular components, molecular functions, biological processes, and regulatory pathways associated with the differentially expressed mRNAs revealed the downstream mechanisms of the Cyr61 gene. The differentially expressed genes were affected for small cell lung cancer, axon guidance, Fc gamma R-mediated phagocytosis, MAPK signaling pathway, focal adhesion, insulin resistance, and metabolic pathways.

Based on pathway analysis, thirteen major mRNAs were found to be involved in the MAPK signaling pathway. The MAPK signaling pathway is essential in regulating numerous cellular processes including inflammation, cell stress response, cell differentiation, cell division, cell proliferation, metabolism, motility, and apoptosis. The role of the MAPK pathway in cancer, immune disorders, and neurodegenerative diseases has been well recognized. The results of gene microarray analysis suggest that downregulation of Cyr61 regulates the cellular function either by regulating the signaling pathway or by regulating the expression of some others genes. This warrants further study.

The mechanisms by which overexpression of Cyr61 increases cancer risk are likely to be multifactorial. Part of the effects seen appears to be due to the direct interaction between Cyr61 protein with the MAPK signaling pathway and p53 as well as an effect on cell apoptosis [7, 50, 51]. An effect mediated by the MAPK signaling pathway is given additional support by our transcriptome analysis. In addition, based on our findings, Cyr61 may also contribute to malignant transformation indirectly through its effects on other genes, which were found to be differentially expressed in this study and have been implicated in malignant transformation. These genes including SERPINF1, CRABP2, and TXNIP [5254]. REG3G, ACAT2, and ATF5 have also been found in various cancers [5557]. Several of the differentially expressed lncRNAs we show in Table 4 have also been found to be associated with cancer development, including CH17-360D5.2, SOX2-OT, and HOXA11 [5860]. These findings suggest that the effect of Cyr61 on malignant transformation is multifactorial and not simply due to a direct interaction with the MAPK signaling pathway.

Additionally, we also found that the differently expressed genes are involved in other pathophysiological processes unrelated to cancer. For mRNAs, for example, PRR21 may be involved in stress responses that are related to phosphorylation of mitochondrial proteins [61]. Mutation of PIGW is associated with West syndrome and hyperphosphatasia with mental retardation syndrome [62]. Rbm19 plays an important role in the development of intestinal epithelium cells [63]. Thus, the genetic changes associated with Cyr61 may be involved in the occurrence and development of many diseases.

GO analysis predicted that up- and downregulated mRNAs were associated with several biological processes, cellular components, and molecular functions in HEK 293 cells with downregulated Cyr61. The GO annotation indicated that these gene products were attributed to the progress of a variety of biochemical reactions, such as cell-cell adhesion, negative regulation of cell migration, protein binding, mRNA processing, cell-cell adherens junction and poly(A) RNA binding. Cell adhesion and migration are the central processes in the development and maintenance of multicellular organisms. Tissue formation during embryonic development, wound healing, and immune responses all require the orchestrated movement of cells in particular directions to specific locations. Errors during these processes may have serious consequences, including intellectual disability, vascular disease, tumor formation, and metastasis [64].

It has been reported that Cyr61 mediates numerous cellular activities, including cell migration, adhesion, apoptosis, and proliferation [65]. As a cell adhesion protein, through binding to particular integrins, Cyr61 associates with the adhesion activities of endothelial cells, fibroblast cells, smooth muscle cells, and monocytes [6669]. Numerous studies have shown that Cyr61 is involved in cell migration, including smooth muscle cell migration, fibroblast-like synoviocyte invasion, and cancerous adenocarcinoma [7072], and is seen as a promising therapeutic target. Nevertheless, the mechanism by which Cyr61 participates in cell adhesion and migration is still not completely clear because both cell adhesion and migration are complex biological processes. Our research shows that a large number of abnormal genes are involved in cell adhesion and migration. It is anticipated that this can provide additional research directions for Cyr61-involved diseases. Subsequently, exploration of the correlations between abnormal genes and signal pathways is critical for developing new biomarkers for the early diagnosis and therapeutic surveillance of Cyr61-involved diseases.

We also studied the correlation between Ang II and Cyr61 in HEK293 cells. It is well recognized that the RAS is an important hormonal system in humans and exists not only as a circulating and paracrine but also as an intracrine system [20, 73]. There have been many well-documented pathophysiologic functions of the intrarenal RAS. Increasing lines of evidence have indicated the RAS was involved in a variety of diseases including tumors, depression, glaucomatous, cardiovascular diseases, kidney diseases, and Parkinson’s disease [7478]. The pathological processes include fibrosis, inflammation, and immune responses, and it also induces inflammation, renal cell growth, mitogenesis, apoptosis, migration, and differentiation [79], and the adverse effects of the RAS on these diseases are due mainly to increased levels of Ang II, which exerts its effects through the classical ACE-AngII-AT1R axis [80].

In this study, we knocked down the Cyr61 gene of HEK293T cells by using the Cyr61 CRISPR/Cas9 KO plasmid and Cyr61 HDR transfection plasmid and then treated the HEK293T cells with Ang II. The cell functions of these HEK293T cells were then observed including the cell cycle, cell proliferation, and cell apoptosis. Our results showed that the expression of Cyr61 protein was significantly decreased after knocking down the Cyr61 gene by Cyr61 CRISPR/Cas9 KO plasmid and HDR transfection plasmid in HEK293T cells compared with the controls, and knockdown of the Cry61 gene can promote cell proliferation and inhibit apoptosis significantly. These results may be related to promoting the cell cycle. In this study, cell-cycle distribution analysis showed significantly decreased cell populations in the S phase and increased cell populations in the G2/M phase when the Cyr61 gene was knocked down. Our results are similar to those of some previous studies. For example, Xiong et al. reported that Ang II significantly inhibited cell survival, induced cell cycle arrest, and enhanced cell apoptosis in HK-2 cells [81].

Interestingly, the HEK293T cells with knockdown of the Cyr61 gene treated with Ang II have a higher cell proliferation rate and G2/M phase and lower apoptosis rate than the Ang II-added group of cells. We further examined the expression of Cyr61 after treating with Ang II. Consistent with the findings of Hilfiker et al. [82], we found that Cyr61 transcript levels were elevated in HEK293T cells exposed to Ang II compared with controls. Furthermore, the Cyr61 protein transcript levels of HEK293T cells with the Cyr61 gene knocked down were elevated significantly after treating with Ang II compared to cells without Ang II treatment. In contrast, the proliferation and apoptosis of Cyr61-downregulated + Ang II cells were not different to control cells, although the Cyr61 level of downregulated + Ang II cells was lower than that of the controls was. This probably means that Cyr61 is not the only factor which is involved in Ang II-induced apoptosis and antiproliferation of these cells and suggests that some others factors and signaling pathways, such as AT2 receptor (angiotensin type 2 receptor), p38 signaling pathway, and IGF-IR-PI3K-Akt signaling are also involved in these processes [8385]. Accordingly, further research into the mechanisms of Ang II-induced apoptosis and antiproliferation is needed. Taken together, downregulation of Cyr61 is of benefit to the cell cycle and is antiapoptotic, and Ang II induces HEK293T cell injury, at least in part, by upregulating the expression of Cyr61 protein.

There are some limitations in our study. The cells transfected with the vector that expresses puromycin N-transferase (PAC) were not involved, and the unspecific effects of PAC on HEK293T cells were not observed in the present study. The puromycin selection may enhance Cas-9 expression in the transfected cells so unspecific effects by Cas-9 cannot be ruled out. lncRNAs can play their biological functions by interacting with mRNAs. The biological meanings of differentially expressed lncRNAs will be more convinced if the lncRNAs and mRNA coexpression network can be performed between the validated lncRNAs and their related mRNA based on the correlation analysis.

5. Conclusions

Cyr61 is mutated in a variety of diseases. We report for the first time that many lncRNAs and mRNAs are significantly upregulated or downregulated and that specific metabolic pathways may play important roles in the downstream events of the Cyr61 gene. The significantly altered lncRNAs and mRNAs levels may be related to many pathological processes. It appears that Cyr61 is involved in Ang II-induced injury in HEK293T cells. It is envisaged that studies of the functional mechanisms of these differentially expressed lncRNAs and mRNAs as well as exploration of the pathways involved will provide novel targets for Cyr61-involved diseases.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Ethical Approval

The study was approved by the Committee of Animal Care and Use of Youjiang Medical University for Nationalities, and all procedures were performed according to the National Institute of Health Guidelines.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Authors’ Contributions

Junjie Wang, Dongdong Fu, and Soulixay Senouthai carried out the experimental work and conceived of the study and participated in its design and coordination. Yan Jiang and Rentong Hu participated in the design of the study and together with Yanwu You performed the statistical analysis. Junjie Wang and Yanwu You drafted the manuscript. All authors read and approved the final manuscript. Junjie Wang, Dongdong Fu, and Soulixay Senouthai are co-first authors.

Acknowledgments

This study was supported by grants from the National Natural Science Foundation of China (No. 81560271), the Key Project of Scientific Research of the Guangxi Colleges and Universities (No. KY2015ZD092), and the Program of Natural Science Foundation of Guangxi (Nos. 2014GXNSFAA118253 and 2017GXNSFDA198005). The authors would like to thank Dr. Dev Sooranna, Imperial College London, for editing the manuscript.