
RNA sequencing has become the method of choice for comprehensive identification and quantification of gene transcripts. Development of next generation sequencing (NGS) technology probably enables large scale RNA sequencing, thereby analysis of the change in gene expression profile becomes easier for researchers. In most cases, when cells encounter physical, chemical, or biological stimuli, the cells respond by changing expression of the genes encoded in their genome. Gene transcripts are the resulting products of the cell’s response. Therefore, identification of transcripts that are increased or decreased might be an invaluable resource for analysis of the cell’s responses to the stimulus. Statistically significant change in expression of a gene in the cells between two experimental conditions is regarded as a differentially expressed genes (DEGs). Through analysis of readouts and changing amount of DEGs, the signaling networks working in the cells in response to the stimulus can be determined.
According to the World Health Organization (WHO), cancer has long been a highly ranked disease, leading deaths worldwide, accounting for nearly 10 million deaths in 2020. In Korea, gastric cancer was the most common cancer and it remained the second most common cancer in 2020 [1]. Discovery of DEGs that might contribute to formation of cancer must be a valuable foundation for development of cancer therapeutics and even cancer vaccines. Because development of cancer takes time, treatment of normal cells with a strong cancer inducer is a good alternative method for discovery of the cancer-promoting DEGs in a short period of time. A well-known carcinogen, N-methyl-N’-nitro-N-nitrosoguanidine (MNNG), has been broadly used to induce diverse cancer formation in normal cells and animals since its successful induction of stomach cancer formation in treated rats [2]. MNNG-induced cancer models have been used in determining the effects of
In the current study, an analysis of MNNG-induced DEGs in normal stomach cells was performed by comparison with DEGs from established stomach cancer cell lines at approximately change level, protein-protein interactions (PPIs) in the pathway for cancer formation. These analyses resulted in discovery of new DEGs and their PPIs in the pathways for cancer formation.
Normal human stomach cell line, HS738, and human stomach cancer cell lines, AGS and MKN45, were cultured under DMEM containing 10% fetal bovine serum (FBS). DMSO or MNNG dissolved in DMSO was added to the culture medium to be 50 μM. Two incubation periods, 6 hours or 24 hours, were adopted for monitoring direct induction or indirect induction included, respectively. An RNeasy kit from Qiagen was used for purification of total RNA from the treated cells. Sequencing of total mRNA was performed at E-Biogen (Seoul, South Korea) using NGS.
An analysis of the sequenced mRNA information provided in Excel format was performed using Excel based DEG analysis (ExDEGA, ver. 3.0.1) provided by E-Biogen. Collection of DEGs from the MNNG-treated cells was based on
GO were performed using a web-based DAVID Bioinformatics Resources 6.8. The database for annotation, visualization and integrated discovery (DAVID) V6.8 comprises a full knowledgebase update and provides a comprehensive set of functional annotation tools to understand the biological meaning behind the long list of genes. An analysis of the connection between DEG and biological pathways was based on the KEGG database. Analysis of functional protein association networks of the selected genes was performed using the STRING web site, from which a biological database and web resource of known and predicted protein-protein interaction (PPI) are available. In this study, only the known PPIs were used to draw networks, which were visualized using Cytoscape software (version 3.5.1).
DEGs from normal stomach cells treated with MNNG were listed from top to bottom according to their induction or reduction folds (≥3.0 folds, log2 ≥1) and determination of significance was based on their
Functional categorization of the MNNG-regulated DEGs was performed using DAVID analysis. Only the genes with known functions were considered allowing duplication for the analysis. Of the up-regulated DEGs, the highest number participated in “signaling transduction”, and the “cellular response to DNA damage stimulus” was the second most common GO. Given that MNNG causes damage to DNA, it was assumed that a number of up-regulated DEGs belonged to the GO categories of “cellular response to DNA damage” and “signal transduction” (Figure 2). In the case of the reduced genes, a large number of genes belonged to GO of cell proliferation, such as “cell division”, “mitotic nuclear division” or “sister chromatid cohesion” (Figure 2). Results of these analyses indicate that the MNNG-induced DEGs are responsible for signaling pathways required for cancer formation, while the reduced DEGs are involved in cell division.
Next, we determined the number of MNNG-regulated DEGs in common with the DEGs from two stomach cancer cell lines (AGS and MKN45). In spite of its well-known cancer-promoting activity, few MNNG-regulated DEGs were expressed in common with the two cancer cell lines. Of the DEGs from MNNG-treated normal cells, 10 and eight DEGs were commonly up-regulated in AGS or MKN45 cells, respectively, while an increase in seven DEGs was observed across all of the cells (Figure 3). In the case of the MNNG-reduced DEGs, 33 DEGs were commonly down-regulated across all of the cells, while only two DEGs were commonly down-regulated in AGS cells. No common gene was accordingly decreased between MNNG-treated stomach cells and MKN45 stomach cancer cells (Figure 3). The number of contra-regulated DEGs was also determined, where expression was opposite between the compared cell types. Of the MNNG-induced DEGs, 49 DEGs belonged to that category. The relationship of gene expression was much stronger between the two cancer cells compared with that of MNNG-treated normal cells. A large number of DEGs in the two cancer cell lines have the same direction of expression: 746 up-regulated and 1440 down-regulated DEGs were overlapped between each other (Figure 3).
Individual genes in the MNNG-regulated DEGs were mapped for retrieval of a PPI network using the STRING database [8]. The STRING analysis was performed using genes that could be searched from “KEGG pathways in cancer”. Species was limited to “
Some of the MNNG-regulated DEGs were not in common with the DEGs from the cancer cells. An analysis of PPIs of DEGs that belonged exclusively to MNNG-treated normal cells was performed using STRING, and their roles in possible pathways were examined. PPI between TUBA1C and TUBA1B was identified from the DEGs which showed less expression in 24 hours than in 6 hours since MNNG treatment. These two proteins have been found to have signal roles in degenerative neuronal diseases, including “Parkinson disease”, “Amyotrophic lateral sclerosis”, and “Alzheimer disease” (Figure 5A). Next, an analysis of MNNG-regulated DEGs in normal stomach cells at 6 hours after treatment was performed for identification of PPIs. Among the up-regulated DEGs, four groups of PPIs were identified, however, the number of PPIs with identified roles in the KEGG pathway was limited. PPIs between LAMB4 and LAMB1 and between HIST1H2BH and HIST1H2BM may have roles in cytoskeleton rearrangement (Figure 5B). Three groups of PPIs were identified from the down-regulated DEGs at 6 hours after MNNG treatment. Most of the proteins involved in the identified PPIs had roles in the “PI3K-Akt signaling pathway”. In particular, the “PI3K-AKT signaling pathway” was the only known pathway for the PPI between PIAPH2 and PIAPH3 (Figure 5C). These results indicate that MNNG-induced gene expression at early phase shows an increase of genes for cytoskeleton proteins and a decrease of genes for the PI3K-AKT signaling pathway.
Transformation of normal cells to cancer cells essentially relies on changes in expression of genes required during the transformation [9]. Regardless of its manner of expression, whether it is genetically induced or epigenetically induced, it is obvious that a gene expression pattern that is different from normal is prerequisite. Therefore, attaining an understanding of DEGs can provide the underlying mechanism of cancer development. Development of a rapid sequencing technique known as NGS has enabled a systemic approach to analysis of DEGs. Due to the capacity for rapid sequencing, almost all of the mRNA identity in a cell can be identified in a short period of time. In this study, because MNNG is recognized as a cancer-inducer to a broad range of cancers
Previously, study of MNNG-induced gene expression was conducted by others using an oligonucleotide microarray containing ∼8,000 probe sets [11]. According to the results, up-regulated genes belonged to extracellular matrix remodeling (i.e.
Two differentiated methods of analysis were employed in this study. The first analysis was a com-parison between MNNG-regulated DEGs and DEGs from the two cancer cell lines, to determine whether DEGs contributed to cancer at early stage development or later stage maintenance. The second analysis was that of protein interaction networks of the DEGs with roles in the cancer pathway. Simply indicating an increase or decrease of expression of a gene cannot determine the difference between normal cells and cancer cells. Actual functions in cells are performed by proteins that interacted with other molecules and each other, therefore analysis of PPIs of protein products of the DEGs results in a clearer understanding of the function of the gene and its contribution to cancer formation [13, 14]. Here a couple of significant MNNG-induced PPIs that might have the capacity to promote stomach cancer formation were identified. Of the MNNG up-regulated DEGs at 6 hours, PPIs between LAMB4 (laminin beta 4) and LAMB1 (laminin beta 4), and between HIST1H2BH (histone H2B type 1-H) and HIST1H2BM (histone H2B type 1-M) PPIs were exclusively identified (Figure 5B). Therefore, in MNNG-induced early phase, even when not reaching cancer status, cytoskeletal changes and histone modification could occur. Subsequently, PPIs found in cancer cells transformed from normal cells by MNNG treatment were identified by analysis of MNNG-induced DEGs overlapped with DEGs from the two cancer cell lines. From the larger number of PPIs identified, proteins including FGFR1, EGF, GNAI1, GNB4, and ITGB1 had centered roles. FGFRs 1-4, which are transmembrane proteins, have oncogenic or anti-oncogenic roles depending on the ligands and cellular environments [15]. EGF regulates important cellular processes, such as proliferation, differentiation, and survival by binding to its receptor, EGFR [16]. GNAI1 (guanine- nucleotide binding subunit I subunit alpha-1) and GNB4 (guanine-nucleotide binding subunit beta-4) are parts of heterotrimeric signal-transducing proteins [17]. ITGB1 (Integrin subunit beta 1) is a cell surface receptor that works as a collagen receptor [18]. Collectively, these identified proteins indicate that MNNG-induced disturbance at the area of cellular communication with the environment, cell surface, might be critical to inducing cancer formation. In addition to the changes in membrane, the PI3K-AKT signaling pathway was the down-stream signaling pathway that was mainly modified by MNNG treatment. Changes in interaction of membrane proteins together with down-regulation of the PI3K-AKT signaling pathway lead to disruption of proper cellular responses. Recent research papers used MNNG as a proven cancer-inducer for examining anti-cancer effects of the molecules [6, 7]. But these studies had been conducted without considering MNNG-induced genetic changes. This study could enable more precise interpretation about function of the anti-cancer molecules and developing more effective cancer therapeutics by providing target genes or protein networks to be intervened. We expect that these results may shed light on the understanding of gene expression and their protein networks that contributes to neoplasm initiation in normal stomach cells.
암은 전 세계적인 건강문제이다. 암의 종류는 다양하나 암의 발생과정에서의 유사성은 상당히 높다. 모든 암은 유전자의 발현 변화가 있다는 점에서 공통점을 갖는다. 암 발생과정에서 변화되는 유전자의 발현을 이해하는 것은 암치료제나 암백신 개발에 큰 도움을 줄 것이다. 이번 연구에서는 잘 알려진 암발생 물질인 MNNG를 이용하여 정상 위세포에서 암발생 시 변화되는 유전자 발현을 분석하였다. MNNG 처리에 의해 정상 세포와는 다르게 발현되는 유전자인 DEG들을 찾고 이들을 암세포에서 발현되는 DEG들과 비교하였다. 여기에 더해 DEG의 단백질 결과물들의 기능과 단백질들 간의 결합을 분석하여 단순히 유전자의 발현뿐만 아니라 이들 단백질의 신호전달 과정에서의 연관성도 함께 분석하였다. 이 결과 위암 발생과정에서 관여하는 유전자들과 이들 유전자의 단백질 결과물들의 상호 작용을 밝히게 되었다.
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) grant (http://nrf.re.kr) to Dongju Jung (NRF-2019R1I1A3A01061981).
None
Kim TJ, Undergraduate student; Kim MK, Graduate student; Jung D, Professor.