Gastric cancer is a major contributor to global cancer burden, being the third leading cause of cancer death worldwide and in both sexes. This type of cancer is thought to be consequence of a multi-step process, resulting from different genetic and epigenetic changes. Specifically, dysfunction of oncogenes and tumor suppressor genes contributes to this malignant disease, and many candidate genes have been implicated to serve as gastric cancer biomarkers.
In this context, the piwi like RNA-mediated gene silencing 1 (PIWIL1) gene, located in 12q24.33 and having 22 exons, became an attractive target for gastric cancer treatment. PIWIL1 protein is expressed at increased levels in cancer tissues, stem cells and germ cells, but it has been shown to be absent in normal somatic tissues. This means that it could be a potential target for therapy, since most non-cancer cells would not be affected by cytotoxic effects[3-7].
PIWIL1 plays a key role in tumor cell viability, migration and invasion, and its expression is associated with the maintenance of stem-like characteristics of tumors, which in turn contribute to more severe histological grade, advanced stage and worse clinical outcome[8-10].
Wang et al showed that expression of PIWIL1 in gastric cancer tissue was significantly higher than in adjacent-to-tumor tissue (tumor front). They also demonstrated that patients with a lower expression of PIWIL1 presented a significant better overall survival rate compared to patients with a higher expression levels. Additionally, the 5-year survival rate of patients with a higher expression level of PIWIL1 was significantly lower (36.5% vs 67.6%).
Liu et al reported that expression of PIWIL1 progressively increases during cancer development. The expression ratio in normal gastric tissues, atrophic gastritis, intestinal metaplasia and gastric cancers varied from 10% to 76%.
To further investigate the potential functions of the PIWIL1 gene, Liu et al also silenced PIWIL1 by antisense or short hairpin RNA and noted that suppression of this gene inhibited the growth of gastric cancer cells and induced G2/M arrest. Although relevant information regarding the possible role of PIWIL1 in gastric cancer carcinogenesis is provided by the current literature, the exact molecular mechanisms involved in this carcinogenic process remain unclear.
A recently introduced technology, based on the adaptive immune system of prokaryotes and known as type II clustered, regularly interspaced, short palindromic repeats (CRISPR)/associated protein (Cas), has been demonstrated to cleave double-stranded DNA and has emerged as a relevant genome editing tool[13-15].
This technology can be used both to perform permanent gene knockouts and the site-specific integration of a gene (knock-in)[16-19]. Importantly, it allows for the permanent silencing of the target gene, and it also creates a stable and permanent cell line with the desired modification[14,16].
Here we applied CRISPR/Cas9 technology for the first time to knockout PIWIL1 gene in a gastric cancer cell line and analyzed its phenotypic modifications.
MATERIALS AND METHODS
The human gastric cancer cell line AGP01 was maintained in DMEM-F12 medium supplemented with 10% fetal bovine serum. The cell culture grew attached to a plastic flask in a monolayer in a humidified incubator maintained at 37 °C and 5% CO2.
The AGP01 cell line was established by our research group in 2009 from cancer cells present in the ascitic fluid of a female individual with intestinal gastric cancer, located at the antrum and the body region of the stomach, and staged as T3N2M1. The cell line was tested and authenticated by conventional cytogenetics. Recently, the AGP01 cell line was tested by multicolor-fluorescence in situ hybridization (FISH), and results are presented here.
24-color-FISH using all human whole chromosome painting probes
24-color-FISH using simultaneous all human whole chromosome painting (WCP) probes was done as previously reported[21,22]. A total of 20 metaphases was analyzed, using a fluorescence microscope (Axio Imager Z1 mot; Carl Zeiss AG, Oberkochen, Germany) equipped with appropriate filter sets to discriminate between a maximum of five fluorochromes and the counterstain DAPI; the latter was used to induce a GTG-like banding pattern. Image capturing and processing were carried out using ISIS imaging system (MetaSystems, Altlussheim, Germany).
Targeted knockout of PIWIL1 using the CRISPR-Cas9 system
The CRISPR-Cas9 system used was purchased from Dharmacon GE Life Sciences (Lafayette, CO, United States). First, 1 × 104 AGP01 cells/well were seeded in DMEM-F12 medium to a 96-well plate for 24 h. Subsequently, transfection was performed using CR-0046-03-005 (Dharmacon GE Life Sciences) for 48 h. For the transfection procedure, a solution containing 1 μL of the CRISPR RNAs (crRNAs) mixed with the trans-activating small RNA, 2 μL of Cas9 and 7 μL of DMEM-F12 medium (for each well) was prepared first. In another tube, 0.4 μL of DharmafecDUO and 9.6 μL of DMEM-F12 medium (to each well) were added.
The solutions were then incubated for 5 min at room temperature before being combined. After combining, the solution was incubated for 20 min at room temperature, and finally 80 μL of DMEM-F12 medium/10% fetal calf serum (FCS) (per well) was added. At the end of the transfection, all contents (from each well) were transferred to a 24-well plate containing DMEM-F12 medium/10% FCS/1% pen-strep.
After 24 h, samples were treated with 6 μg/mL of puromycin for 72 h to select the resistant clones. Next, 40 cells were plated per well in a 6-well plate to isolate the clones by the filter paper method.
This method consists of using a cut and autoclaved piece of filter paper so that it, after being soaked in trypsin, can be positioned above a single colony of cells that has grown from an isolated cell, allowing for the collection of this clone. Subsequently, each clone grew in a separate well in a 6-well plate, so it could reach the confluence needed to perform DNA extraction and sequencing.
DNA was extracted using the Wizard Genomic DNA Purification Kit (Promega Corporation, Madison, WI, United States) according to the manufacturer’s instructions. For PCR, specific primers targeting the binding region of the purchased crRNA were constructed using the online program Primer3 (Supplementary Table 1).
The quantities of the reagents used in the PCR for a final volume of 12 μL were as follows: 6.25 μL of nuclease-free H2O, 0.5 μL of forward primer (10 ng/μL), 0.5 μL of reverse primer (ng/μL), 4.25 μL of GoTaq Colorless Master Mix 2× (Promega Corporation) and 1 μL of DNA (10 ng/μL).
The conditions using the MasterCycler Gradient thermal cycler (Eppendorf, Hamburg, Germany) were 1 cycle at 95 °C for 3 min for initial denaturation followed by 35 cycles consisting of denaturation at 94 °C for 2 min, primer annealing at 59 °C for 1 min and extension at 70 °C for 2 min, ending with 1 cycle at 70 °C for a final extension for 30 min.
For direct sequencing of the PCR product, the quantities of the reagents used for a final volume of 20 μL were as follows: 15 μL of nuclease-free water, 0.5 μL of forward or reverse primer (10 ng/μL), 0.5 μL of Big Dye, 3 μL of Save Money and 1 μL of the PCR reaction. For this reaction, the ABI PRISM Big Dye Terminator 3.1 Cycle Sequencing Kit (Applied Biosystems, Hercules, CA, United States) was used.
The sequencing was performed using the MasterCycler Gradient (Eppendorf) thermal cycler according to the following 25-cycle thermocycling conditions: denaturation at 96 °C for 50 s, primer annealing at 59 °C for 1 min and extension at 60 °C for 4 min, ending with 1 cycle at 4 °C for 5 min.
After this procedure, a precipitation step was carried out in order to purify the product of the reaction before continuing. For this step, samples were washed with 70% isopropanol and 70% ethanol. Subsequently, sequencing was performed using the ABI Prism 3500 DNA Sequencer (Applied Biosystems). The methodology used was based on the biochemical synthesis of the DNA strand by the Sanger method.
Reference sequence for exon 15 in the PIWIL1 gene was obtained from the National Center for Biotechnology Information (NCBI) and compared with the DNA sequence of the modified cell line as well as the negative control (cell line without the gene knockout). To infer the effect of changes in protein synthesis, we applied Gene Runner v.3.05 (Hastings Software Inc., Hastings, NY, United States; http://www.generunner.com).
Cells were grown in 12-well plates at a density of 2 × 105 cells/well and maintained for 24 h in 5% CO2 at 37 °C. After this period, cells were injured with a 10 μL tip in the center of each well. The medium was then removed to eliminate suspended cells, and wells were washed with 1 × phosphate buffered saline before fresh DMEM-F12 medium/10% FCS/1% pen-strep was added again.
The behavior of cells was observed and photographed immediately after injury and at 6 h, 12 h and 24 h after injury. All experiments were performed in triplicate.
Boyden chamber invasion assay
Boyden inserts (8 μm pores) (BD BiosciencesTM, Franklin Lakes, NJ, United States) were coated with 200 μL of Matrigel (10-13 mg/mL) in 12-well plates. Cells (2 × 105) were seeded in the upper chamber in 1 mL of DMEM without fetal bovine serum. In the lower chamber, DMEM-F12 medium/10% FCS/1% pen-strep was added, functioning as a chemoattractant for the cells present in the upper chamber.
After 48 h, the remaining cells above the filter were removed by scraping with a sterile swab. The cells at the bottom of the filter were fixed with 4% paraformaldehyde and stained with Giemsa. Cells were photographed and analyzed using a light microscope and counted in optical fields (100 ×). All experiments were performed in triplicate.
Total RNA extraction
The mRNA extraction was performed using Promega’s Total RNA Isolation System kit, according to the manufacturer’s specifications. AGP01 and AGP01 PIWIL1 knockout cells were prepared for mRNA extraction. Samples were lysed with lysis buffer containing beta-mercaptoethanol and then diluted in RNA dilution buffer. The samples were centrifuged for 10 min at maximum speed. Then, 95% ethanol was added to ensure adequate membrane binding conditions.
The samples were then transferred to centrifuge columns where the RNA could bind the membrane of the column, facilitating washing to eliminate possible contaminants as well as favoring the extraction of high quality of total RNA.
At the end of the procedure, the RNA was diluted in 60 μL of nuclease-free water. The total RNA was quantified using a Nanodrop spectrophotometer ND-1000 UV-VIS version 3.2.1 (Nanodrop Technologies, Wilmington, DE, United States). The RNA quality was also evaluated by analyzing the A260/A280 ratio according to the manufacturer’s specifications. Purified RNA was stored at -80 °C for the microarray expression assay.
For the microarray assay, a one-color microarray-based gene expression analysis kit (Agilent Technologies, Santa Clara, CA, United States) was used according to the protocol provided by the manufacturer. The gene expression profile was evaluated in both cell lines (AGP01 with and without PIWIL1 knockout).
The total RNA obtained during the extraction phase was used as the template for the synthesis of the first cDNA strand by reverse transcription using T7 RNA polymerase. Synthesis of the second cDNA strand was used as the template for the in vitro transcription reaction for cRNA production. The cRNA was then incorporated into the fluorochrome 3-cyanine (Cy-3) using the Low Input Quick Amp Labeling kit (Agilent Technologies) according to the protocol provided by the manufacturer. Thereafter, the cRNA purification process was carried out.
The cRNA was quantified by a spectrophotometer (pmol/L), by which it was possible to analyze the absorbance ratio (260 nm/280 nm) and the cRNA (ng/μL) concentration in each sample. After, hybridization was performed for 17 h in a hybridization chamber at 65 °C at 10 rpm. After this period, the slide was washed and immediately scanned in the Agilent G4900DA SureScan Microarray System.
The following setup was used to scan the microarray slides for one color: scan region of 61 mm × 21.6 mm, 5 μm scan resolution, dye channel of green. Next, the images were obtained by using Feature Extraction v10.10 software, and the data were analyzed with GeneSpring GX 9.0 and IPathwayguide (Advaita Bioinformatics Company, Plymouth, MI, United States) programs. Gene identification followed a restriction criterion with a fold-change of > 2.
Differential expression and gene ontology enrichment analysis
To identify differentially expressed (DE) mRNAs, we compared the probes’ expression profiles before and after PIWIL1 knockout. Probes with a mean fold-change < 0.5 and a mean fold-change of > 2 [|Log2(fold-change)| > 1] were selected for differential analysis. Student’s t test was performed in 1222 selected probes, false discovery rate adjustment was performed, and genes with an adjusted P-value of < 0.05 were tagged as DE. For the volcano plot, the fold-change and the P-value of all probes were used, and probes were tagged as DE following the previous criteria. Gene ontology enrichment was performed using org.Hs.eg.db and Gostats R libraries. All graphical and statistical analyses were performed in the R platform (R Core Team, 2017, Vienna, Austria; https://http://www.r-project.org/).
24-color-FISH using all human WCP probes
24-color-FISH using all human WCP probes revealed in AGP01 cell line a complex karyotype as follows: 63,XX,inv(1)(p12q43),der(1)(1pter->1p12::1q43->1p12::9p12->9pter),+der(1)t(1;12)(q21;q12),+del(1)(p12),+del(2)(p12),+der(2)t(2;8)(p12;q11.2),inv(3)(p21q13)x2,+inv(3)(p21q13),+der(3)t(3;5)(p14;q13),t(4;14)(p12;q11.2),dic(4;12)(p15;q12),+der(4)t(2;4)(p or q?;q12),del(5)(p13)x2,+der(5)x2,del(6)(q12)x2,+del(6)(p21),inv(7)(p12;q11.2),der(7),der(9)t(9;acro)(p21;p12),der(9)t(6;9)(p12;p12),t(13;13)(p10),der(15)(:q13->p11::p11->qter),-15,der(16)t(X;16)(q or p?;q23),+der(16)t(12;16)(q12;q23),21p+,t(22;22)(p10) (Figure 1).
Figure 1 AGP01 cell line multicolor-FISH showing several chromosomal changes, including the monosomy of chromosome 12, where the PIWIL1 gene is located.
From the analysis, we verified that the unique chromosome 12 (where PIWIL1 gene is located) remains intact, without translocations or derivative chromosomes. It is important to note that the monosomy of chromosome 12 agrees with the sequencing result, since the 7 bp insert sequence was observed in hemizygous status.
Targeted knockout of the PIWIL1 gene using the CRISPR-Cas9 system
PIWIL1 gene knockout was successful, as determined by Sanger sequencing. The latter revealed an insertion of 7 adenines in the PIWIL1 gene sequence (Figure 2), which caused a frameshift mutation that impaired protein synthesis.
Figure 2 Electropherograms showing the insertion of seven adenines in the PIWIL1 gene sequence after using the CRISPR-Cas9 system (B) in comparison with the reference sequence of the negative control (A).
Prediction of the encoded protein indicated a premature stop codon (Figure 3), suggesting that this insertion generates a truncated protein consisting of 573 amino acids (the wild-type contains 861) with a loss-of-function phenotype, which means that knockout was efficient.
Figure 3 Amino acid sequences of the wild-type PIWIL1 protein and the PIWIL1 protein after the insertion of seven adenines by the CRISPR-Cas9 system.
Notably, the PIWIL1 knockout cell line remained viable and could be used for further experiments.
The PIWIL1 gene knockout caused a significant decrease in AGP01 migration capacity after 24 h (P < 0.01; Figure 4), which is consistent with the fact that this protein is related to various pathways that regulate cell motility.
Figure 4 Analysis of the migration capacity of the AGP01 cell line with and without PIWIL1 gene knockout.
T0: Immediately after injury; T6: 6 h after injury; T12: 12 h after injury; T24: 24 h after injury. NC: Negative control. aP < 0.01. Two-way ANOVA, Bonferroni post-test. Photomicrography of AGP01 cell migration. A: Immediately after injury; B: 6 h after injury; C: 12 h after injury; D: 24 h after injury. The black lines represent approximation of the edges over time, demonstrating the migration capacity of the cells.
Boyden chamber invasion assay
Also, PIWIL1 gene knockout caused a significant decrease in AGP01 invasiveness (P < 0.001; Figure 5), which is also consistent with the fact that this protein is related to various pathways that regulate cell motility.
Figure 5 Analysis of the invasion capacity of the AGP01 cell line with and without the PIWIL1 gene knockout.
Statistically significant difference between groups was shown by the Student’s t test (cP < 0.001). Photomicrography of the cell invasion assay demonstrating the decrease in the number of cells that invaded when PIWIL1 was knocked out. NC: Negative control.
mRNA array and gene ontology enrichment
Differential analysis: By comparing expression profiles after PIWIL1 permanent knockout in the AGP01 cell line, a total of 251 mRNA were found to be DE [adjusted P-value of < 0.05 and |Log2(Fold-Change)| > 1]. A total of 43 mRNA probes were up-regulated and 208 were down-regulated (Figure 6). The DE genes are described in Supplementary Table 2.
Figure 6 Volcano plot comparing gene expression after PIWIL1 permanent knockout in the gastric cancer AGP01 cell line.
Differentially expressed probes [adjusted P-value of < 0.05 and |Log2(Fold-Change)| > 1] are on superior left and right areas (red). mRNAs involved in invasion and migration processes are in green.
We performed a functional analysis by grouping all DE mRNAs. This approach revealed that a total of 35 genes (five up-regulated and 30 down-regulated) encoded proteins involved in invasion and migration cellular processes (Supplementary Figure 1; Supplementary Tables 3 and 4).
According to the current literature, 9 of these 35 genes (DOCK2, ZNF503, PDE4D, ABL1, ABL2, LPAR1, SMAD2, WASF3 and DACH1) are possibly related to the mechanisms used by PIWIL1 to promote carcinogenic effects related to migration and invasion, since their functions are consistent with the changes observed (being up- or down-regulated after knockout).
Evidence indicating that PIWIL1 is an oncogene that regulates several cellular mechanisms important for the carcinogenic process of many types of cancers is rising in recent years. PIWIL1 seems to be especially implicated in proliferating activity of cancer cells[8,10,26,27].
Here, we introduce a new gastric cell line, i.e., AGP01, and provide a first description of its karyotype. As is usual for cancer cell lines, there were slight changes from cell to cell concerning number of chromosomes or single cell chromosomal rearrangements. The here-given karyotype was the most frequently observed one. Interestingly, inversions, dicentrics and reciprocal translocations of homologous acrocentric (#13, #15, #22) was observed. Overall, gain of the following regions was present: #1, #3, large parts of #2 and #5, 8q11.2 to 8qter, 12pter to 12q12. Besides, the following regions were under-represented: 2pter to 2p12, 9pter to 9p12, 15q12 to 15qter. These imbalances are in concordance with the literature.
In AGP01, we performed for the first time an in vitro knockout experiment of PIWIL1 gene using the CRISPR-Cas9 system. It could be shown that absence of this gene significantly impairs migration and invasion capacity of AGP01 cells. Thus, the AGP01 cell line behaves like that previously reported for gastric cancer cells or lung adenocarcinoma. Together, these studies suggest that PIWIL1 expression is strongly associated with an increased aggressiveness of cancer cells.
According to Wang et al, one of the mechanisms by which PIWIL1 regulates the migration and invasion of cancer cells is by promoting the expression of MMP2 and MMP9, two important metalloproteinases involved in the degradation of the extracellular matrix, thereby creating paths for the locomotion of cancer cells.
Additionally, Amaar et al demonstrated that the over-expression of PIWIL1 down-regulates the tumor suppressor gene IGFBP5, a member of the insulin-like growth factor binding protein family and whose expression is implicated in suppressing epithelial-mesenchymal transition and reducing the expression of E-cadherin and HIF1α, indicating that is it critically related to cancer progression[33,34].
Data obtained from our gene expression experiments also provided corroborating evidence that the PIWIL1 gene plays a key role in cancer cell migration and invasion because several genes involved in these cellular processes were observed as DE when the cell lines were compared before and after PIWIL1 knockout.
Many studies have demonstrated the oncogenic activities of the DOCK2, ZNF503, PDE4D, ABL1, ABL2, LPAR1, SMAD2 and WASF3 genes and their relation to tumor aggressiveness in several types of cancer, including gastric cancer[35-58]. Interestingly, PIWIL1 knockout led to a decreased expression of these genes as well as an increased expression of the tumor suppressor gene DACH1, demonstrating that PIWIL1 plays a crucial role in the pathway of development and progression of gastric cancer, and is likely a promising candidate for therapeutic intervention.
Zhu et al reported that the over-expression of DACH1 impaired the proliferation and invasion ability of lung adenocarcinoma cells in vitro via the down-regulation of PRX3, an oncoprotein required for the maintenance of mitochondrial function and tumorigenesis. DACH1 expression also inhibited epithelial-mesenchymal transition and metastasis by affecting TGF-β signaling and decreased proliferation of cancer cells by inducing cell cycle arrest at the G2/M phase[60,61].
Regarding the reported oncogenes, Rahrmann et al observed that PDE4D is over-expressed in human prostate cancer and demonstrated that the knockdown of this gene reduced the growth and migration of prostate cancer cells in vitro as well as the growth and proliferation rate of prostate cancer xenografts in vivo. Delyon et al demonstrated that PDE4D is also over-expressed in melanoma cell lines and pinpointed this gene as a regulator of cell invasion by interacting with FAK through RACK1, constituting a signaling pathway that when activated promotes tumor progression and metastasis.
Recent studies have determined the role of ABL members from the tyrosine kinase family, ABL1 and ABL2, in the development of many types of solid tumors. These proteins induce the activation of actin polymerization machinery by modulating the expression of several MMPs to promote morphological changes, including the formation of membrane protrusions and altered cell adhesion. Consequently, activation of ABL1 and ABL2 in cancer cells promote enhanced proliferation, migration and invasion, as well as drug resistance[39-42].
DOCK2 (dedicator of cytokinesis) belongs to the DOCK family of proteins and is expressed in hematopoietic cells. According to Kulkarni et al, DOCK2 has been reported to activate Rac, which is known to regulate several crucial processes, including lymphocyte migration, activation and differentiation of T cells. Wang et al knocked out DOCK2 in a B-cell lymphoma cell line and observed a decrease in Rac1 expression. Additionally, analysis of the growth curves of both cell lines demonstrated that the DOCK2 knockout grew less than DOCK2, as evidenced by the lower cell proliferation.
ZNF503 is expressed in the mammary gland and other tissues, and there is a high incidence of association between this gene deregulation and tumor aggressiveness in several kinds of tissues, such as lung, kidney and intestine. Shahi et al performed scratch and 3D Matrigel culture assays in two mammary epithelial cell lines to analyze the cell motility and migration. Both assays demonstrated that cell lines with ZNF503 knockout migrated less, did not close the gaps, and inhibited invasiveness when compared to the control cells. These data indicate that ZNF503 promotes cellular invasion and migration, and high levels of this gene are closely related with poor patient survival, breast cancer progression and increased metastasis.
Yu et al demonstrated that the lysophosphatidic acid receptor 1 (LPAR1) gene is related to migration and invasion in ascites from ovarian cancer and is expressed at higher levels in metastatic cell lines, when compared to non-metastatic cell lines. They also observed that the presence of high levels of lysophosphatidic acid are directly connected to cell migration stimulation, and LPAR1 silencing reduced lysophosphatidic acid-induced invasion. Additionally, in breast tumors, a higher expression of LPAR1 is related to a worse lung metastasis-free survival rate.
Wiskott-Aldrich syndrome protein family 3 (WASF3) is an important gene, which has C-terminal domains that are responsible for actin polymerization activation, playing a role in cell proliferation and migration. The WASF3 gene is normally over-expressed in several types of tumors, such as breast cancer, osteosarcoma and prostate cancer[52,53]. In gastric cancer, little is known about this gene; however, since micro (mi)RNAs and their targets are considered potential biomarkers for gastric cancer, Wang et al performed a luciferase assay and western blotting to investigate the relationship between miR-218 and WASF3. Their results demonstrated that over-expression of WASF3 harms miR-218 and results in the inhibition of cell proliferation and migration, suggesting that WASF3 is over-expressed in gastric cancer and induces cell proliferation and migration. Additionally, in qRT-PCR, WASF3 mRNA expression levels were higher when compared to normal gastric cell lines.
Smad2 is the first intracellular protein in the signaling cascade of the TGF-β1 signaling pathway, which is involved in the progression of gastric cancer. In advanced stages of cancer, TGF-β1 acts as an oncogene, regulating multiple cellular functions, including stimulation of proliferation, differentiation and the inhibition of apoptosis[55,57].
Interestingly, Lv et al observed that TGF-β1 levels in peritoneal lavage fluid are directly connected to peritoneal metastasis. Corroborating evidence was provided by Shinto et al, whose experiments demonstrated that p-Smad2 expression was higher in diffuse-type tumors and in peritoneal metastasis cases. Notably, the AGP01 cell line used to perform the PIWIL1 knockout in our study was obtained from a patient with peritoneal metastasis, and we found SMAD2 was over-expressed.
Taken together, these data reinforce the idea that PIWIL1 plays a crucial role in the signaling pathway of gastric cancer, regulating several genes involved in migration and invasion processes; therefore, its use as a therapeutic target may generate promising results in the treatment of gastric cancer, mainly in patients with peritoneal carcinomatosis, which is a condition associated with poor prognosis and a decreased overall survival.
Gastric cancer (GC) remains a major public health problem, having the third highest incidence of death worldwide. Piwi like RNA-mediated gene silencing 1 (PIWIL1) is involved in regulation of widespread biological processes, including stem cell proliferation, embryogenesis, growth, and development, and has been found to be frequently over-expressed in various tumor types, including GC. Previous studies have demonstrated that PIWIL1 is implicated in improving tumor malignant behavior. PIWIL1 expression has been shown to be absent in normal somatic tissues, making it a very intriguing target for therapy. We attempted to investigate the role of PIWIL1 on the migration and invasion capacity of metastatic GC cells, using the AGP01 cell line, as well as checking the expression status of genes and proteins involved in these cellular processes, in order to elucidate the mechanisms by which PIWIL1 provokes tumorigenic effects and to shed light on potential new strategies to target PIWIL1
Many aspects of gastric carcinogenesis remain elusive, and much effort has been made to improve patient prognosis. The PIWIL1 has been identified as a novel extremely highly expressed gene in many types of cancer and its expression in GC tissue is related to poorer overall survival, suggesting that high expression of PIWIL1 is associated with poor prognosis and that it could be used as a predictive marker or even a target for therapy. Although PIWIL1 has been correlated with worse outcome, the involved mechanisms remain unclear, and many hypotheses are being tested. Once the upstream and downstream signaling pathways of PIWIL1 are elucidated, it will be possible to create new therapeutic strategies for gastric carcinogenesis, in order to improve the overall health of patients affected by this disease.
We performed permanent knockout of the PIWL1 gene to verify phenotypic modifications in the AGP01 metastatic GC cell line, as well as alterations in expression level of mRNA and protein, in an attempt to better understand the mechanisms by which PIWIL1 promotes tumor malignant behavior. This research demonstrates the importance of studying PIWIL1 in GC, since data obtained through the achievement of our objectives showed that this protein has a crucial role in gastric carcinogenesis, promoting molecular and phenotypic alterations compatible with enhanced tumor aggressiveness. The elucidation of the role of PIWIL1 protein in cancer cell invasion and migration will pave the way for developing potential clinical interventions, aiming to control GC dissemination.
We applied CRISPR/Cas9 technology to knockout the PIWIL1 gene in a metastatic GC cell line, and analyzed its phenotypic modifications, as well as alterations in gene and protein expression. CRISPR-Cas9 technology was considered in 2015 as one of the most important technological advances of science. Mainly, it allows permanent silencing of the target gene and also creates a stable and permanent cell line with the desired modification. By this way, multiple experiments can be carried on, including long term evaluation of the downstream events caused by the molecular alteration, as well as discovering potential pathways influenced by the studied gene. Therefore, after permanent knockout of PIWIL1 in the AGP01 cell line, we analyzed phenotypic modifications by performing wound-healing and Boyden chamber invasion assays, to assess migration and invasion, respectively. Moreover, aiming to shed light on the molecular mechanisms used by PIWIL1 to make changes in the migration and invasion capability of cells, we carried out proteomic and microarray assays, using multidimensional protein identification technology (commonly known as MudPIT) and a one-color microarray-based gene expression analysis kit, respectively.
PIWIL1 gene knockout was successfully performed and confirmed by Sanger sequencing, which revealed an insertion of seven adenines in the PIWIL1 gene sequence. In silico prediction of the encoded protein pointed to the appearance of a premature termination codon, suggesting that this insertion generates a truncated protein with a loss-of-function phenotype. PIWIL1 knockout promoted a significant decrease in cell migration and invasion capacity (P < 0.01 and P < 0.001, respectively), which is consistent with data present in the literature demonstrating that this protein is implicated in several signaling pathways that regulate cell motility. By comparing expression profiles after PIWIL1 knockout, a total of 251 mRNA were found to be differentially expressed, with 43 up-regulated and 208 down-regulated mRNA. A functional analysis grouping all differentially expressed mRNAs demonstrated that 35 genes encoded proteins were involved in invasion and migration cellular processes. After extensive review of data presented in the literature, we selected 9 of these 35 genes (DOCK2, ZNF503, PDE4D, ABL1, ABL2, LPAR1, SMAD2, WASF3 and DACH1) as possibly related to the mechanisms used by PIWIL1 to promote carcinogenic effects related to migration and invasion, since their functions are consistent with the changes observed (being up- or down-regulated after knockout). Additionally, the analysis of proteomic data revealed that PIWIL1 knockout caused modification in the expression of 27 proteins involved in epithelial-mesenchymal transition (EMT). Twenty-two oncoproteins related to EMT promotion, including FGFR1, PCNA, ACTN4, GSN and TUBB3, were expressed in the AGP01 cell line before knockout and reduced to a level that were not detectable by the technique after knockout. On the other hand, PIWIL1 knockout caused an increase in the expression of six proteins implicated in EMT suppression, such as ACSM3, ADGRG1 and ANPEP, that were absent in AGP01 before knockout. To the best of our knowledge, this is the first report describing molecular alteration compatible with phenotypic alterations after permanent knockout of PIWL1 in GC. Detailed mechanisms leading to PIWIL1 over-expression in cancer as well as the pathways by which this protein improves the malignant phenotype should be further investigated.
In the current study, we pioneered the performance of an in vitro knockout of the PIWIL1 gene by using the CRISPR-Cas9 system, and found that absence of this gene significantly impaired the migration and invasion capacity of the AGP01 cell line, besides modifying mRNA and protein expression of potential molecular targets involved in the EMT process. The results of such experiments contributed to understanding of the mechanisms used by PIWIL1 to promote alteration in migration and invasion capacity of gastric cells during tumorigenesis, and also revealed the participation of new players related to PIWIL1 expression, such as FGFR1, PCNA, ACTN4, PDE4D and SMAD2. Our results demonstrated that knockout of PIWIL1 promotes several changes in cell phenotype, suggesting the critical role of the PIWIL1 oncogene in GC, and confirmed the hypothesis that PIWIL1 expression provokes migration, invasiveness and EMT as potential mechanisms of improved tumor aggressiveness. The presented findings open new perspectives for molecular interventions in GC
Definite silencing of PIWIL1 by the CRISPR-Cas9 system resulted in robust findings favoring the discovery of new mechanisms involved in gastric carcinogenesis. The presented results must be validated by other researchers, and if confirmed, might lead to innovative interventions aiming to treat GC.