Expression profiling of SCN8A and NDUFC2 genes in colorectal carcinoma
Summary. Aim: The expression differences of SCN8A (which encodes type VIII alpha subunit of voltage gated sodium channel) and NDUFC2 (which encodes C2 subunit of Complex I enzyme in oxidative phosphorylation) genes were evaluated in paired colorectal cancer (CRC) tissues which was relied on our partial transcriptome analysis data in cancer cell lines. Materials and Methods: A total of 62 paired tissues of CRC patients (34 male, 28 female) were included in the study. The mRNA levels of SCN8A and NDUFC2 genes were determined by using real-time PCR (qRT-PCR and semiquantitative PCR). Results: SCN8A gene expression level was significantly lower in tumor tissues (p = 0.0128) and in the patients with the age below 45 years (p = 0.0049). There were also meaningful relationships between the gender, grade of CRC, tumor location, histopathological classification, and SCN8A expression. There was no NDUFC2 differential expression. However, the tumors taken from right colon had significantly lower NDUFC2 expression. Conclusion: Although the voltage gated sodium channels (VGSCs) and Complex I (CI) were associated to a number of diseases including different types of cancers, the different subunits of CI and individual members of VGSCs seem to be cancer type-specific in varying proportions.
Submitted: January 20, 2015
*Correspondence: Fax: +90 342 3603928;
Abbreviations used: CI — Complex I; Cp — crossing point; CRC — colorectal cancer; GAPDH — glyceraldehyde-3-phosphate dehydrogenase; Nav1.6 — sodium channel, voltage gated, type VIII, alpha subunit; qRT-PCR — real-time PCR; VGSCs — the voltage gated sodium channels.
Colorectal cancer (CRC) is one of the most common types of cancer worldwide . In the United States, CRC is the third in the leading causes of death due to cancer among men and women . Thirty percent of CRC cases have familial basis while 67% are sporadic, and the rest 3% of the cases have inherited CRC syndromes. The complex etiology of sporadic CRC involves a combined effect of genetic and environmental risk factors .
Various genes are affected in terms of expression during the multistep malignant transformation into CRC . Among those genes, plasma membrane ion channel encoding genes are of special importance since they are involved in intracellular signaling events and they activate specific cellular responses. However, ion channels research is considered relatively novel in cancer studies  and the transcriptional regulation is not well characterized for the sodium channel genes . The evidences suggesting the correlation of ion channels with cancer are accumulating and the voltage gated sodium channels (VGSCs) were associated to invasion and metastasis behaviors of several types of cancers [4, 6–8]. Among those channels, sodium channel, voltage gated, type VIII, alpha subunit (Nav1.6) is encoded by SCN8A gene. SCN8A has specific roles in neurons and skeletal and cardiac muscles . The gene is located on 12q13 and the transcript variants of this gene encode two distinct proteins which are 1980 and 1939 amino acids (aa) in length (http://www.ncbi.nlm.nih.gov/gene/6334). An approximate 40-fold overexpression of SCN8A gene was detected in cervical cancer samples . Another significant increase in SCN8A mRNA expression in prostate cell lines was detected, as well .
Complex I (CI), NADH: ubiquinone oxidoreductase is of critical importance for cell physiology and is the first multiprotein complex involved in the formation of the proton-motive force in mitochondria. Its main function is to transport the electrons by oxidation of NADH followed by reduction of ubiquinone [11, 12]. The genes encoding various CI subunits were associated with diseases including severe neonatal lactic acidosis, Leigh syndrome, neonatal cardiomyopathy with lactic acidosis, leukodystrophy with macrocephaly, and hepatopathy with renal tubulopathy [12, 13] and malignancies including breast, colon and thyroid cancers [14–17]. Human CI is composed of 45 subunits making it the largest and most complicated enzyme in the oxidative phosphorylation system. Thirty-nine of the subunits of mammalian CI are encoded by the nuclear genome and the rest is encoded by the mitochondrial genome . Among those, NDUFC2 gene is located on 11q14.1 and the transcript variants of this gene encode three distinct proteins which are 119, 88, and 96 aa. in length (http://www.ncbi.nlm.nih.gov/gene/4718). The downregulation of NDUFC2 gene in papillary thyroid carcinoma was previously reported . The chromosomal cytoband region 11q14.1 spanning NDUFC2 gene was shown as amplified in estrogen receptor negative breast cancer and this region was associated with worse prognosis of this disease, as well . Additionally, NDUFC2 gene was suggested as being one of the candidate genes involved in microsatellite instability of colon adenomas .
SCN8A and NDUFC2 genes are members of different and almost unrelated pathways which also are of importance in terms of CRC development and progression. However, this pilot study was relied on our partial transcriptome analysis results in which we determined altered expression of SCN8A and NDUFC2 genes in various cancer cell lines. Current study was planned to confirm our previous findings and it is aimed to evaluate the gene expression differences of SCN8A and NDUFC2 genes in CRC tissues.
MATERIALS AND METHODS
Subjects and sample collection. A total of 62 (34 male, mean age — 58.8±11.2; 28 female, mean age — 53.8 ± 12.5) CRC patients who underwent surgery were included in the study. Paired tumor and tumor-surrounding normal tissues of the patients were used. Samples were placed into liquid nitrogen immediate after the surgery and kept in −80 °C until RNA isolation.
Fifty-four of the samples were adenocarcinoma while 8 were mucinous adenocarcinoma, according to the histopathological classification. Six of the tumors were located in cecum, 12 in right colon, 2 in transverse colon, 11 in left colon, 10 in sigmoid colon, and 21 were located in rectum.
This study was approved by the local ethical committee in concordance with the declaration of Helsinki (University of Gaziantep, Clinical research ethical committee, protocol # 219/11.06.2013). Written consent was taken from all of the patients.
RNA isolation and gene expression analysis. RNA isolation was made by using RNA isolation kit (QIAGEN, Germany) according to the instructions of the manufacturer. RNA concentrations were determined by using spectrophotometer (NanoDrop, ND-1000, USA). cDNA synthesis was performed by using AMV first strand cDNA synthesis kit (New England Biolabs, UK) according to the instructions of the manufacturer. Expression primers spanning SCN8A and NDUFC2 mRNAs were designed by using Primer-BLAST (http://www.ncbi.nlm.nih.gov/tools/primerblast/). Primer sequences were presented in Table 1. GAPDH (glyceraldehyde-3-phosphate dehydrogenase) gene was used as the housekeeping gene used for normalization of the expression results of the genes of interest.
Table 1. The characteristics of the primer pairs used for expression analysis
Gene expression levels were determined by using real-time PCR (qRT-PCR) on LightCycler® 480 instrument. PCR conditions were as follows: 2 min at 50 °C, 10 min at 95 °C, 40 cycles of 15 s at 95 °C, 30 s at 60 °C, and 30 s at 72 °C. Semi-quantitative PCR analysis was conducted prior to qRT-PCR analysis in order to confirm the findings (Fig. 1). Thermal conditions were as follows; 5 min at 95 °C, 35 cycles of 30 s at 95 °C, 30 s at appropriate annealing temperature presented in Table 1, 30 s at 72 °C and 5 min at 72 °C.
Fig. 1. Semi-quantitative PCR results of NDUFC2 and SCN8A genes. Top: Two percent agarose gel image of randomly selected 8 samples’ results were shown. GAPDH gene was utilized for normalization of the expression values by using. NormalTarget Gene/Reference Gene-TumorTarget Gene/Reference Gene formula. Bottom: The graphs show the average expression levels of the samples. The intensity of the bands were quantitated and densitometric analysis was made by using ImageJ software (v1.46r). A general decrease in the expression values of NDUFC2 and SCN8A genes in tumor tissues compared to normal tissues were prominent.
NDUFC2; Mean, N:1,24, T:1,16; Std. Deviation, N:0,62, T:0,44; Std. Error, N:0,079, T:0,056
SCN8A; Mean, N:1,29, T:1,14; Std. Deviation, N:0,56, T:0,54; Std. Error, N:0,072, T:0,069
Statistical analysis. After normalization by using GAPDH expression, Cp (Crossing point) values were used for the statistical analysis. Calculation was made by using 2-ΔCp formula as previously described . The expression levels of SCN8A and NDUFC2 in tumor and paired normal tissues were compared. The relationships between age, gender, grade of the tumor, tumor location, histopathological classification, and gene expression levels were also evaluated. Wilcoxon test was used for the statistical analyses by using GraphPad Prism (v.5) software. All statistical tests were two-tailed and p < 0.05 was considered as significant.
SCN8A gene expression level was significantly lower in tumor tissues compared to paired normal tissues (p = 0.0128) (Table 2, Fig. 2). The decrease in expression was much more significant in the patients with the age below 45 years (p = 0.0049). Of the other patient parameters evaluated, there were also meaningful relationships between the gender (female; p = 0.012), grade of CRC (grade II; p = 0.02), tumor location (cecum; p = 0.031 and right colon; p = 0.021), histopathological classification (adenocarcinoma; p = 0.017) and SCN8A expression CRC (Table 2).
There was no NDUFC2 gene expression difference between the tumors and paired normal tissues in overall (Fig. 3). However, the tumors taken from right colon had significantly lower NDUFC2 expression when compared to paired normal tissues (p = 0.042) (Table 2).
Table 2. The patients characteristics, tumor grades, tumor locations, histopathological classifications and the results of the statistical analyses
NS — not significant.
Fig. 2. Calculated 2-∆Cp values of tumor and paired normal tissues for SCN8A gene. N — normal tissue; T — tumor tissue; *denotes statistical significance
Fig. 3. Calculated 2-∆Cp values of tumor and paired normal tissues for NDUFC2 gene. N — normal tissue; T — tumor tissue
As of the other types of cancer, the complexity of CRC is a major obstacle in understanding of the genetics of the disease . Therefore, molecular classification of CRC still remains a challenge that researchers have to overcome in the coming years . The implementation of effective CRC biomarkers may be realized by better understanding of cancer. In our study, we have evaluated the gene expression differences of SCN8A and NDUFC2 genes both of which are involved in different and important cellular processes, in CRC tissues and surrounding healthy colonic mucosa, by using qRT-PCR method. The rationale of the study was based on our preliminary results.
Differential expression and/or genetically alterations of the members of the VGSCs might be involved in the progression of the malignant growth . In our study, the expression of SCN8A gene which encodes Nav1.6 was significantly lower in CRC samples. Interestingly, this decrease was much more significant in the patients with the age below 45 years. Accumulating data suggest that VGSCs, especially Nav1.5, Nav1.6, and Nav1.7, were associated with various types of cancers in terms of oncogenic advantage and invasion capacity [7, 10, 20, 21]. However, the expression profiles of the individual members of VGSCs are highly variable. In a study aiming to investigate functional expression of alpha subunits of VGSCs in human ovarian cancer; only Nav1.1, Nav1.3, Nav1.4, and Nav1.5 mRNAs were overexpressed in ovarian cancers among the 12 members of VGSCs . In the same study, only Nav1.2, Nav1.4, Nav1.5, and Nav1.7 mRNAs were overexpressed in metastatic ovarian cancers. Taken together, the individual members of VGSCs were differentially expressed specific to cancer type and its malignancy potential. Previous studies reported overexpression of SCN8A gene in cervical cancer and prostate cancers [9, 10]. However, in our study there was downregulation of SCN8A mRNA in CRC tissues. This decrease, appears to be CRC-specific, may be due to the aforementioned cancer-type dependent expression of VGSCs. The other alpha subunits of VGSCs could also be associated with CRC. Additionally, the alternatively-spliced members of VGSCs may also involve in CRC progression . The overexpression of VGSCs might be an important factor in triggering metastasis and motility due to the altered electrophysiological status of the cell . The altered expressions of VGSCs, either as a whole or as particular members, might be key in developing potential diagnostic or therapeutic targets . Because this was a pilot study, we did not determine the protein expression levels of SCN8A and we did not determine the expression levels of VGSCs other than SCN8A, either. Further studies are needed.
Except for the samples excised from right colon, there was no significant expressional change for NDUFC2 gene in this study. This was contrary to our expectation. This might be due to the heterogenicity of cancer cells. NDUFC2 gene encodes one of the 45 subunits of CI . CI is a huge complex and it is placed in a complicated physically and functionally interconnected cellular network . In our previous study , we reported downregulation of NDUFC2 in papillary thyroid carcinoma. Additionally, CI was associated to a number of diseases including different types of cancers [14,15]. However, the type of CI subunit and whether they are encoded from nuclear or mitochondrial genome seems to be disease-specific. Therefore, the evaluation of CI subunits as a whole and the protein expressions of these individual subunits are warranted. NDUFC2 protein has three isoforms (http://www.ncbi.nlm.nih.gov/gene/4718). The primer sets used in this study for gene expression analysis of NDUFC2 gene were spanning the shared sequence of all three isoforms. Hence, the expression levels we determined were in overall; rather than isoform-specific. The potential differential expressions of NDUFC2 splice variants should be tested. As the SCN8A gene, NDUFC2 gene expression was decreased for the samples taken out from the right colon. This data should also be verified in larger cohorts.
CRC is one of the most common cancers worldwide and diagnosis of CRC could be conducted in a non-invasive manner by the usage of molecular markers. The discovery of new markers and the combinations of them with the others is thought to increase the diagnostic efficiency. This pilot study represents the expression profiling of SCN8A and NDUFC2 genes each of which encode one of the subunits of two important (and unrelated with each other) complexes in terms of CRC development and progression. Further studies aiming to evaluate the relationship between CRC and other CI subunits and other members of VGSCs in terms of mRNA and protein expressions are warranted.
Conflict of interest
The authors declare that there are no conflicts of interest.
No Comments » Add comments