Next Article in Journal
Control of the Epithelial-to-Mesenchymal Transition and Cancer Metastasis by Autophagy-Dependent SNAI1 Degradation
Next Article in Special Issue
Nevoid Basal Cell Carcinoma Syndrome: PTCH1 Mutation Profile and Expression of Genes Involved in the Hedgehog Pathway in Argentinian Patients
Previous Article in Journal
Epigenetically Altered T Cells Contribute to Lupus Flares
Previous Article in Special Issue
Targeting Hedgehog Signalling through the Ubiquitylation Process: The Multiple Roles of the HECT-E3 Ligase Itch
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 8
Issue 2
10.3390/cells8020128
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regulation of Survivin Isoform Expression by GLI Proteins in Ovarian Cancer

by
Diana Trnski
1,†,
Maja Gregorić
2,†,
Sonja Levanat
1,
Petar Ozretić
1,
Nikolina Rinčić
1,
Tajana Majić Vidaković
1,3,
Držislav Kalafatić
4,5,
Ivana Maurac
4,
Slavko Orešković
4,5,
Maja Sabol
1,* and
Vesna Musani
1,*
1
Ruđer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
2
Zagreb Health School, Medvedgradska 55, 10000 Zagreb, Croatia
3
PP Orahovica, Pustara 1, 33513 Zdenci, Croatia
4
Department of Obstetrics and Gynaecology, University Hospital Centre Zagreb, Petrova 13, 10000 Zagreb, Croatia
5
School of Medicine, University of Zagreb, Petrova 13, 10000 Zagreb, Croatia
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cells 2019, 8(2), 128; https://doi.org/10.3390/cells8020128
Submission received: 31 December 2018 / Revised: 1 February 2019 / Accepted: 4 February 2019 / Published: 6 February 2019
(This article belongs to the Collection Hedgehog Signal Transduction in Physiology and Disease)
Browse Figures
Review Reports Versions Notes

Abstract

:
Ovarian cancer (OC) is the most lethal female gynecological malignancy, mostly due to diagnosis in late stages when treatment options are limited. Hedgehog-GLI (HH-GLI) signaling is a major developmental pathway involved in organogenesis and stem cell maintenance, and is activated in OC. One of its targets is survivin (BIRC5), an inhibitor of apoptosis protein (IAP) that plays a role in multiple processes, including proliferation and cell survival. We wanted to investigate the role of different GLI proteins in the regulation of survivin isoform expression (WT, 2α, 2B, 3B, and Δex3) in the SKOV-3 OC cell line. We demonstrated that survivin isoforms are downregulated in GLI1 and GLI2 knock-out cell lines, but not in the GLI3 knock-out. Treatment of GLI1 knock-out cells with GANT-61 shows an additional inhibitory effect on several isoforms. Additionally, we examined the expression of survivin isoforms in OC samples and the potential role of BIRC5 polymorphisms in isoform expression. Clinical samples showed the same pattern of survivin isoform expression as in the cell line, and several BIRC5 polymorphisms showed the correlation with isoform expression. Our results showed that survivin isoforms are regulated both by different GLI proteins and BIRC5 polymorphisms in OC.

1. Introduction

Most malignant ovarian tumors in adult women are epithelial ovarian tumors. These neoplasms are classified into different morphological categories according to the appearance of epithelial cells in serous, mucinous, endometrial, clear cell carcinomas, mixed, squamous, transient, and non-differentiated tumors [1]. Serous adenocarcinomas are the most common and make up 75% of all ovarian cancer. Each histological type of ovarian tumor and the degree of differentiation of the cells is associated with certain molecular-genetic changes.
The Hedgehog signaling pathway is a developmental pathway involved in formation of various tissues and organs, including the ovary. In mammals, canonical signal transduction is triggered by binding of the ligand Hedgehog (HH) to the transmembrane receptor Patched (PTCH). This leads to receptor internalization and exposure of the protein Smoothened (SMO) on the cell surface. Activation of SMO triggers the release of transcription factors GLI (GLI1-3) from Suppressor of Fused (SUFU) and the translocation of GLI to the nucleus. GLI, then, activates transcription of HH target genes which are involved in cell cycle regulation, proliferation, adhesion, epithelial-mesenchymal transition, self-renewal, and pathway autoregulation [2,3]. Recently, survivin has also been identified as a novel HH target gene. Vlčková et al demonstrated that expression of survivin is directly regulated by GLI2, more specifically by the GLI2ΔN isoform [4]. A study by Brun et al showed that survivin is overexpressed in HH-driven medulloblastoma. They suggest it may represent a novel therapeutic target for this disease [5].
Survivin is the smallest member of the inhibitors of apoptosis (IAP) family of proteins in mammals [6]. The gene for survivin, BIRC5 (baculoviral IAP repeat containing 5), is located on human chromosome 17q25 [7]. The BIRC5 gene has four dominant exons (1, 2, 3, and 4) and two cryptic exons (2B and 3B). Alternative splicing of its pre-mRNA produces at least five different mRNAs encoding five different proteins with different intracellular distribution and anti-apoptotic properties: wild type survivin (S WT), survivin 2α (S 2α), survivin 2B (S 2B), survivin 3B (S 3B), and survivin-ΔEx3 (S ΔEx3) (Figure 1) [8].
S WT, discovered in 1997, contains a BIR domain necessary for its anti-apoptotic function [9]. It is expressed during development but is not expressed in most differentiated adult tissues. The overexpression of survivin is common in almost all tumors and is indicative of decreased overall survival, increased rate of recurrence, and resistance to therapy [6].
S 2α is the smallest isoform with a truncated BIR domain [10]. The reports about its function in cancer is conflicting, while some studies report that it binds to and attenuates anti-apoptotic activity of WT survivin or correlates with expression in lower stages of the disease [8,10], other studies report its association with resistance to treatment [11,12].
S 2B is the longest survivin isoform, but the insertion of cryptic exon 2B interrupts the BIR domain [13]. Its function remains unclear. While some studies report that it promotes cell death, or that its expression is inversely correlated with the tumor stage [8,11,14,15,16,17], others report its expression being correlated with treatment resistant cancers [18].
S 3B lacks the carboxyl-terminal coiled-coil domain of WT survivin [19]. It is considered cytoprotective [20,21] and its overexpression has been correlated with shorter overall survival and resistance to therapy [11,22].
S Δex3 contains a bipartite nuclear localization signal (NLS) and localizes to nucleus in cancer cells [23,24]. Its expression in tumors is associated with aggressive disease and unfavorable prognosis [24,25].
Survivin exhibits cell-cycle-dependent expression that is mostly controlled at transcription level. Its accumulation during mitosis is also influenced by posttranslational modifications that affect its stability. When expressed during mitosis, it is located in various sections of the mitotic apparatus such as centrosomes, microtubules, and anaphase spindles, and remains of the mitotic apparatus [6]. The homologous deletion of survivin results in early embryo death, which shows its crucial role in cellular development, differentiation, and homeostasis. It is selectively expressed in cancer cells, but not in healthy tissues. Excess accumulation of survivin is associated with the development of disease, disease recovery, and prognosis in various cancers, including bladder cancer, cervix, head and neck, prostate, skin, and ovarian [7]. A global deregulation of the BIRC5 gene mediated by oncogenes (including STAT3, E2F or mutated RAS) or by the loss of tumor suppressors such as p53 or APC, accounts for the selective expression of survivin in cancer [26]. Since survivin is expressed exclusively in cancer cells, it is an interesting target for targeted therapy and new methods for detection of survivin, as well as new inhibitors, are constantly being developed [18,26,27,28].
Several BIRC5 polymorphisms have been studied and have been associated with susceptibility to lung [29], gastric [30], bladder [31], oral [32], and liver cancer [33] as well as age of onset in ovarian cancer [7] and survival in colorectal cancer [8] and breast cancer [34].
X-ray crystallography has shown that survivin is organized as a dimer [6]. Since its isoforms exhibit various apoptosis-related properties, it is believed that the formation of survivin heterodimers with its isoforms may be important for regulating the function of survivin [8]. Since survivin isoforms can affect the activity of wild type survivin, and the transcription of this gene is regulated by HH-GLI signaling, it is important to investigate which factors contribute to the expression of certain isoforms. In this paper our goal was to examine the role of individual GLI transcription factors in the transcriptional regulation of survivin isoforms.

2. Materials and Methods

2.1. Generation of Knockout Cell Lines

SKOV-3 cell line [3] was maintained in DMEM supplemented with 10% FBS and penicillin/streptomycin. The CRISPR/Cas9 system was used to generate the knock-out cell lines. sgRNA molecules targeting the area surrounding the STG site of GLI1, GLI2, and GLI3 genes were generated using the online tool at crispr.mit.edu web site. Five top ranking sgRNA molecules were selected, and the one closest to the ATG site of each gene was selected. The sgRNA oligos were annealed, phosphorylated and cloned into the pX330-U6-Chimeric_BB-CBh-hSpCas9 (plasmid #42230, Addgene, Watertown, MA, USA) [35] into the BbsI restriction site, and the resulting vectors were sequenced on the ABI Prism 310 sequencing machine. The oligo sequences are listed in Supplementary Table S1. The SKOV-3 cell line was transfected with the generated vectors using the Xfect transfection reagent (Clontech, Mountain View, CA, USA). At 24 h post-transfection, the cells were trypsinized and plated into 96-well plates at the density of 1 cell/well. Several cell lines were propagated from single cells, and protein expression for each of the three GLI proteins was examined by Western blot to select the lines with the best knock-out for each protein.

2.2. Cell Culture Experiments

Cells were plated in 6-well plates, treated with two concentrations of GANT-61 (10 and 20 µM). The cells were collected 24 h after treatment. For the transfection of GLI1 (pcDNA4NLSMTGLI1, a kind gift from Prof. F. Aberger, Austria), GLI2 (p4TO6MTGLI2, a kind gift from Prof. M. Stevanovic, Serbia), GLI3 (pcDNA4/TO/GLI3richtig, a kind gift from Prof. M. Stevanovic, Serbia), and GLI3R (EGLI3-PHS, a kind gift from Prof. R. Toftgard, Sweden) the cells were seeded in 6-well plates and transfected the next day with 5 µg of plasmid DNA using the X-fect reagent (Clontech) following the manufacturer’s instructions.

2.3. Clinical Samples

Forty ovarian carcinoma (OC) samples and nine healthy fallopian tubes (FT) tissue samples (excised for reasons other than malignant transformations) were collected at the Department of Obstetrics and Gynaecology, University Hospital Centre Zagreb, School of Medicine, University of Zagreb. Patient blood samples were also collected. Blood samples collected from 74 healthy elderly women with no personal history of ovarian cancer were used as control for comparing polymorphism frequencies.
All patients gave their informed consent before the samples were taken, and samples were collected with the approval of the hospital’s Ethical Committee (number of approval 02/21 AG, issued by the University Hospital Centre Zagreb on March 7th, 2017). The study was carried out following the rules of the Declaration of Helsinki Principles. All tissue samples taken during surgery were immediately placed in a vial containing 1 mL RNALater solution (Invitrogen, Carlsbad, CA, USA), kept at 4 °C overnight, and DNA and RNA were extracted the following day. DNA from tissue samples was extracted by the standard phenol-chloroform method and from blood by the salting out method.

2.4. Western Blot

Proteins were extracted by sonication in RIPA buffer supplemented with protease inhibitors (Roche, Basel, Switzerland). Protein concentration was determined using the Pierce BCA Protein Assay Kit (Thermo Fisher, Waltham, MA, USA). Western blot was performed as previously described [36] using the following primary antibodies: GLI1 (V812, Cell Signaling Technology, 1:1000, Danvers, MA, USA), GLI2 (ARP31885_T100, Aviva Systems Biology, 1:1000, San Diego, CA, USA), GLI3 (19949-1-AP, ProteinTech, 1:600, Rosemont, IL, USA), GFP for detection of GFP tagged GLI3R (Abcam ab6556, 1:1000, Cambridge, UK), and β-actin (60008-1-Ig, ProteinTech, 1:4000) and γ-tubulin (Santa Cruz Biotechnology, sc-7396, 1:500, Dallas, TX, USA) as loading controls.

2.5. Expression Analysis

RNA was extracted using TRIZol reagent (Thermo Fisher) and reverse-transcribed using TaqMan Reverse Transcription Reagents (Applied Biosystems, Foster City, CA, USA). qRT-PCR was performed on a CFX96 machine (Bio-Rad Laboratories, Hercules, CA, USA) using Sso Fast EVAGreen Supermix (Bio-Rad).
The PCR conditions were as follows: initial denaturation at 95 °C for 3 min, 40 cycles of 95 °C for 15 sec, 61 °C for 1 min, and finally melting curve from 70 °C to 95 °C. All experiments were performed at least in duplicate. Primer sequences are listed in Supplementary Table S2. Expression levels were calculated using the 2−ΔΔCt formula, relative to the housekeeping gene TBP [37]. For samples with no expression of a specific target after 40 cycles of qRT-PCR, the Ct was set to 40 to enable the statistical analysis of data [38].

2.6. Genotyping

Whole coding region (including alternate exons S 2α, S 2B, and S 3B) was genotyped, including the six SNPs in the BIRC5 promoter and four SNPs in the 3′UTR region selected from the National Center for Biotechnology Information SNP database (http://www.ncbi.nlm.nih.gov/snp). Primers were designed using the Primer3 online tool (http://bioinfo.ut.ee/primer3/) [39]. Thirteen PCR fragments were analyzed using high resolution melting analysis on High-Resolution Melter (HR-1, Idaho Technology, Salt Lake City, UT, USA) as described in Cvok et al, 2008 [40] followed by Sanger sequencing (ABI PRISM 310 Genetic Analyzer, Applied Biosystems). Due to presence of four different polymorphisms in the PCR product of exon 4 and the beginning of 3′UTR, it was directly sequenced. PCR primer sequences and cycling conditions are listed in Supplementary Table S3. Due to high GC content in the DNA sequence of promoter region, in all PCR fragments located in the promoter CG RICH buffer (Roche) was added.

2.7. Statistical Analysis

D'Agostino–Pearson test was used for testing distribution normality of the expression data. For variables which showed normal distribution after logarithmic transformation, independent sample t-test and one-way analysis of variance were used to test the differences in isoforms’ expressions. Otherwise Mann–Whitney test and Kruskal–Wallis test were used. Nonparametric Spearman rank correlation coefficient (ρ) was used to assess the correlation of expression between various isoforms. For comparing polymorphism frequencies Fisher’s exact test (2 × 2) and χ2 test (3 × 2) were used. Two-tailed p values less than 0.05 were considered statistically significant. Statistical analysis was performed using MedCalc, version 18.2.1 (MedCalc Software bvba, Ostend, Belgium).

3. Results and Discussion

3.1. GLI Regulation of Survivin Isoform Expression

To test the contribution of each GLI protein to the expression levels of five different survivin isoforms, we generated SKOV-3 knock-out lines for GLI1, GLI2, and GLI3 proteins and examined the expression of these isoforms by qRT-PCR. All three KO cell lines exhibit reduced levels of full length GLI proteins, as well as GLI2 and GLI3 repressor forms. However, the band corresponding to the GLI2ΔN isoform is still present in the SKOV-3 GLI2 KO cell line (Supplementary Figure S1). The expression levels of all isoforms were significantly downregulated in the SKOV-3 GLI1 KO cell line. In the SKOV-3 GLI2 KO cell line, only the S 3B and S Δex3 isoforms were downregulated, whereas knocking out GLI3 had no effect on either isoform (Figure 2). Therefore, we decided to focus on the SKOV-3 GLI1 KO line.
SKOV-3 wild type (WT) and the GLI1 KO lines were treated with two doses of the GLI1/2 inhibitor GANT-61 to test the responsiveness of different isoforms. GANT-61 showed a dose-dependent additive effect to knocking out GLI1 on the expression of S WT, S 2α, and S 2B, which suggests that these three isoforms are regulated by both GLI1 and an additional GANT-61 target, probably GLI2. The S 2B isoform shows the strongest response to inhibition with GANT-61, and even though it is affected by the GLI1 KO, it is less pronounced than in the other isoforms, suggesting that this isoform is regulated primarily by GLI2. The lack of effect on expression of these isoforms in the GLI2 KO cell line may be due to the persisting GLI2ΔN isoform, which is known to be a much more potent activator than the GLI2 full length protein [41]. The S 3B and S Δex3 isoforms do not respond to additional GANT-61 inhibition, suggesting that they are primarily regulated through GLI1 (Figure 3).
Our results correspond to those of Vlčkova et al [4], who showed that GLI1 and GLI2 have a moderate effect on survivin promotor activation and GLI3 displayed no promotor activation. In their study, they studied the GLI2ΔN isoform separately and showed that this GLI2 isoform has the strongest activation potential.
To validate the role of specific GLI transcription factors in the expression of survivin isoforms we overexpressed all three GLI transcription factors in the SKOV-3 cell line, including the GLI3 repressor (GLI3R) (Supplementary Figure S2). Transfection of GLI1 and GLI2 expectedly activated the HH-GLI pathway, detected by the upregulation of GLI1 and PTCH1 genes. Overexpression of GLI3 showed no effect on HH-GLI pathway activity, whereas overexpression of GLI3R downregulated the pathway (Figure 4A). Only slight changes in gene expression levels of survivin isoforms were observed upon GLI transcription factor transfection, probably due to an already high level of activity of these genes in cancer cells. On the other hand, GLI3R overexpression, which downregulates the HH-GLI signaling pathway (Figure 4A), leads to a significant decrease in expression of all survivin isoforms, except S 2B (Figure 4B).

3.2. Expression of Splice Variants

To verify that the SKOV-3 cell line is the appropriate model for studying survivin isoforms in OC, we compared basal expression levels of all five isoforms in cells to those in OC tissues. The expression patterns were shown to be the same (p < 0.001), implying that this cell line is a valid model for studying survivin isoforms in OC (Figure 5).
The mRNA expression of S WT, S 2α, S 2B, S 3B, and S Δex3 was analyzed in 29 OC samples and compared to FT. Isoforms S WT and S 2α were detected in 96.5% of OC samples, S 2B in 58.6%, S 3B in 68.9%, and S Δex3 in 89.6% of OC samples. In contrast, the only splice variants detected in healthy FT samples were the weakly expressed S WT and S 2α. No other isoforms were detected in FT. The levels of all survivin isoforms expression differed significantly from each other (p < 0.001) (Table 1). The highest expression was determined for the S WT and S 2α isoforms, while S 2B had the lowest expression. OC samples showed significantly higher expression of these two isoforms compared to healthy fallopian tube tissue (p = 0.0001 and p = 0.0015, respectively) (Figure 6).
The expression levels of survivin isoforms was compared to the levels of GLI1, GLI2, and GLI3 expression which were published previously for the same set of OC tissue samples (data not shown) [36]. The only isoform which was correlated with the GLI genes was the S 2B isoform, which correlated with GLI2 expression (ϱ = 0.55, p = 0.019). Interestingly, in SKOV-3 cells the S 2B isoform showed significant upregulation only after transfection of GLI2 (Figure 4B), indicating that this isoform is regulated by GLI2.

3.3. Genotyping

Fifteen different polymorphisms were detected in the analyzed samples (OC and control) (Table 2). Most of the polymorphisms were located in the 5' and 3′UTR regions including the promoter region. All polymorphisms were in the Hardy–Weinberg equilibrium. There was no significant difference in distribution of genotypes and alleles between OC samples and controls for any of variant. However, c.-235A variant was more frequent in OC but with borderline significance (p = 0.053). Likewise, the frequency of heterozygous genotype c.-235G>A was higher in OC with marginal significance (p = 0.051) (Supplementary Table S4).
The allele frequencies for polymorphisms were in accordance with previous reports [29,34,42,43,44,45,46] or in NCBI database (http://www.ncbi.nlm.nih.gov).

3.4. Linkage Disequilibrium

Analysis showed there is a difference in number of variants in linkage disequilibrium between OC samples and controls. OC cases showed increased non-random association of alleles at multiple loci. In case of c.-1547C>T and c.9386T>C, these two polymorphisms always appear in the same combination in all tested OC samples showing complete linkage disequilibrium (r2 = 1.0) (Figure 7).
LD has often been reported between various BIRC5 polymorphisms. Most of the correlations found here have already been described [29,34,44,47,48]. Our group recently described the LD between c.-644C>T and c.9809T>C in breast cancer (submitted). LD between c.235T>C and c.221+209T>C found in almost total LD, has not been described before.

3.5. Correlation of Splice Variant Expression with Genetic Polymorphisms in OC Patients

Seven polymorphisms showed significant association with expression of survivin isoforms.
c.-1547C>T showed the highest association with splice variant expression (Supplementary Figure S1). Samples with homozygous c.-1547T/T genotype showed significantly higher expression of isoforms S 2α (p = 0.016), S 2B (p = 0.036), and S 3B (p = 0.028). Major T allele was also significantly associated with higher expression of S 2α (p = 0.033), and S 3B (p = 0.045) than C allele.
c.9194G>A also showed significant association with variant expression (Supplementary Figure S2). Samples with heterozygous GA genotype had significantly higher expression of S WT (p = 0.0001) and S 2α (p = 0.0005) than AA genotype. Minor G allele was also significantly associated with higher expression of S WT (p < 0.0001) and S 2α (p < 0.0001) compared to A allele.
c.9386T>C and c.10611C>A showed significant association with the expression of S 2α and S 3B (Supplementary Figure S3). Samples with major c.9386TT genotype had significantly higher expression of S 2α (p = 0.032) and S 3B (p = 0.036) than TC genotype. Samples with major c.10611CC also had significantly higher expression of S 2α (P = 0.040) and S 3B (p = 0.036) than CA genotype.
c.-625G>C, c.-235G>A, and c.221+209T>C also showed some significant associations with isoform expression (Supplementary Figure S4). Samples with minor homozygous c.-625CC genotype had significantly lower expression of S 2α in comparison to heterozygous GC genotype. Samples with major G allele of c.-235G>A (p = 0.047) or major T allele of c.221+209T>C (p = 0.045) had significantly higher expression of S 3B isoform than the minor alleles.
Association of BIRC5 polymorphisms and expression of survivin isoforms have not been explored a lot in the literature. There are some reports about the influence of several polymorphisms on total survivin expression, most notably c.-31G>C [46,48,49,50], but also c.9386T>C [51] and c.9809T>C [47]. Antonacopoulou et al studied the correlation of c.-31G>C and c.9386C>T with the survivin isoform expression in colorectal cancer, and they found the same association of c.9386C>T genotypes with the expression of S 2α [8].

4. Conclusions

We have detected survivin isoform expression (S WT, S 2α, S 2B, S 3B, and S Δex3) in the ovarian cancer cell line SKOV-3, that follows the trend of isoform expression levels in its GLI KO variants and in ovarian cancer samples (Figure 5). The BIRC5 gene has been recognized as one of the transcriptional targets of the HH-GLI signaling pathway [4]. The GLI2ΔN isoform has been indicated as the main transcription factor that regulates BIRC5 gene expression. We have previously detected activity of the HH-GLI pathway in ovarian cancer [36] and therefore we wanted to examine the role of the GLI transcription factors not only in the regulation of survivin wild type expression, but also its isoforms. Because of the different roles of certain isoforms and their ability to regulate survivin function after heterodimer formation it is important to investigate how certain isoforms are regulated. Interestingly, our results confirm that survivin is regulated by the HH-GLI signaling pathway, but not all isoforms are regulated in the same manner. More specifically, all isoforms are regulated by GLI1, but isoforms S WT, S 2α, and S 2B can additionally be inhibited with the GLI1/2 inhibitor GANT-61, suggesting that these isoforms can also be regulated through other GANT-61 targets, most likely, GLI2.
Seven BIRC5 polymorphisms are also associated with the expression of survivin isoforms. Therefore, there is an intricate regulation of expression between genetic influence and transcription factors such as GLI proteins to determine the final expression of survivin isoform expression and their anti- and pro-apoptotic role that could determine the fate of the cells.

Supplementary Materials

The following are available at https://www.mdpi.com/2073-4409/8/2/128/s1, Figure S1: Western blots showing the protein expression of GLI1-3 proteins in SKOV-3 knock-out lines; Figure S2: Western blots showing the overexpression of GLI1-3 proteins in SKOV-3 cell line; Figure S3: Association between c.-1547C>T and survivin isoform expression; Figure S4: Association between c.9194G>A and survivin isoform expression; Figure S5: Association between c.9386T>C and c.10611C>A and survivin isoform expression; Figure S6: Association between c.-625G>C, c.-235G>A and c.221+209T>C and survivin isoform expression. Table S1: Oligos used for cloning; Table S2: Primer sequences for splice variant expression; Table S3: Sequences of genotyping primers; Table S4: Association of genotype and allele frequencies.

Author Contributions

Conceptualization, V.M., M.S., and D.T.; methodology, V.M. and M.S.; validation, P.O., S.L., and V.M.; formal analysis, P.O.; investigation, M.G., N.R., and T.M.V.; resources, D.K., I.M., and S.O.; data curation, P.O. and V.M.; writing—original draft preparation, D.T., and M.G; writing—review and editing, M.S., S.L., and V.M.; visualization, D.T.; supervision, S.L. and S.O.; project administration, V.M and S.L.; funding acquisition, S.L.

Funding

This research was funded by Croatian Science Foundation, grant number IP-2016-06-1268 and donation from Terry Fox Foundation. The APC was funded by Croatian Science Foundation.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Bell, S.M.; Schreiner, C.M.; Goetz, J.A.; Robbins, D.J.; Scott, W.J. Shh signaling in limb bud ectoderm: Potential role in teratogen-induced postaxial ectrodactyly. Dev. Dyn. 2005, 233, 313–325. [Google Scholar] [CrossRef] [Green Version]
  2. Sabol, M.; Trnski, D.; Musani, V.; Ozretić, P.; Levanat, S. Role of GLI Transcription Factors in Pathogenesis and Their Potential as New Therapeutic Targets. Int. J. Mol. Sci. 2018, 19, 2562. [Google Scholar] [CrossRef]
  3. Sabol, M.; Car, D.; Musani, V.; Ozretic, P.; Oreskovic, S.; Weber, I.; Levanat, S. The Hedgehog signaling pathway in ovarian teratoma is stimulated by Sonic Hedgehog which induces internalization of Patched. Int. J. Oncol. 2012, 41, 1411–1418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Vlčková, K.; Ondrušová, L.; Vachtenheim, J.; Réda, J.; Dundr, P.; Zadinová, M.; Žáková, P.; Poučková, P. Survivin, a novel target of the Hedgehog/GLI signaling pathway in human tumor cells. Cell Death Dis. 2016, 7, e2048. [Google Scholar] [CrossRef] [Green Version]
  5. Brun, S.N.; Markant, S.L.; Esparza, L.A.; Garcia, G.; Terry, D.; Huang, J.-M.; Pavlyukov, M.S.; Li, X.-N.; Grant, G.A.; Crawford, J.R.; et al. Survivin as a therapeutic target in Sonic hedgehog-driven medulloblastoma. Oncogene 2015, 34, 3770–3779. [Google Scholar] [CrossRef] [PubMed]
  6. Altieri, D.C. Validating survivin as a cancer therapeutic target. Nat. Rev. Cancer 2003, 3, 46–54. [Google Scholar] [CrossRef] [PubMed]
  7. Han, C.H.; Wei, Q.; Lu, K.K.; Liu, Z.; Mills, G.B.; Wang, L.-E. Polymorphisms in the survivin promoter are associated with age of onset of ovarian cancer. Int. J. Clin. Exp. Med. 2009, 2, 289–299. [Google Scholar]
  8. Antonacopoulou, A.G.; Floratou, K.; Bravou, V.; Kottorou, A.; Dimitrakopoulos, F.-I.; Marousi, S.; Stavropoulos, M.; Koutras, A.K.; Scopa, C.D.; Kalofonos, H.P. The survivin -31 snp in human colorectal cancer correlates with survivin splice variant expression and improved overall survival. Anal. Cell Pathol. 2010, 33, 177–189. [Google Scholar] [CrossRef]
  9. Ambrosini, G.; Adida, C.; Altieri, D.C. A novel anti-apoptosis gene, survivin, expressed in cancer and lymphoma. Nat. Med. 1997, 3, 917–921. [Google Scholar] [CrossRef]
  10. Caldas, H.; Honsey, L.E.; Altura, R.A. Survivin 2α: A novel Survivin splice variant expressed in human malignancies. Mol. Cancer 2005, 4, 11. [Google Scholar] [CrossRef]
  11. Boidot, R.; Vegran, F.; Lizard-Nacol, S. Predictive value of survivin alternative transcript expression in locally advanced breast cancer patients treated with neoadjuvant chemotherapy. Int. J. Mol. Med. 2009, 23, 285–291. [Google Scholar]
  12. Zaki Dizaji, M.; Ghaffari, S.H.; Hosseini, E.; Alizadeh, N.; Rostami, S.; Momeny, M.; Alimoghaddam, K.; Ghavamzadeh, A. Survivin isoform expression in arsenic trioxide-treated acute promyelocytic leukemia cell line and patients: The odd expression pattern of survivin-2α. Asia Pac. J. Clin. Oncol. 2017, 13, e21–e30. [Google Scholar] [CrossRef] [PubMed]
  13. Mahotka, C.; Wenzel, M.; Springer, E.; Gabbert, H.E.; Gerharz, C.D. Survivin-deltaEx3 and survivin-2B: Two novel splice variants of the apoptosis inhibitor survivin with different antiapoptotic properties. Cancer Res. 1999, 59, 6097–6102. [Google Scholar] [PubMed]
  14. Krieg, A.; Mahotka, C.; Krieg, T.; Grabsch, H.; Müller, W.; Takeno, S.; Suschek, C.V.; Heydthausen, M.; Gabbert, H.E.; Gerharz, C.D. Expression of different survivin variants in gastric carcinomas: First clues to a role of survivin-2B in tumour progression. Br. J. Cancer 2002, 86, 737–743. [Google Scholar] [CrossRef] [PubMed]
  15. Ling, X.; Cheng, Q.; Black, J.D.; Li, F. Forced expression of survivin-2B abrogates mitotic cells and induces mitochondria-dependent apoptosis by blockade of tubulin polymerization and modulation of Bcl-2, Bax, and survivin. J. Biol. Chem. 2007, 282, 27204–27214. [Google Scholar] [CrossRef]
  16. Cho, G.-S.; Ahn, T.S.; Jeong, D.; Kim, J.-J.; Kim, C.-J.; Cho, H.-D.; Park, D.-K.; Baek, M.-J. Expression of the survivin-2B splice variant related to the progression of colorectal carcinoma. J. Korean Surg. Soc. 2011, 80, 404–411. [Google Scholar] [CrossRef] [PubMed]
  17. Shi, K.; An, J.; Shan, L.; Jiang, Q.; Li, F.; Ci, Y.; Wu, P.; Duan, J.; Hui, K.; Yang, Y.; et al. Survivin-2B promotes autophagy by accumulating IKK alpha in the nucleus of selenite-treated NB4 cells. Cell Death Dis. 2014, 5, e1071. [Google Scholar] [CrossRef] [PubMed]
  18. Vivas-Mejia, P.E.; Rodriguez-Aguayo, C.; Han, H.-D.; Shahzad, M.M.K.; Valiyeva, F.; Shibayama, M.; Chavez-Reyes, A.; Sood, A.K.; Lopez-Berestein, G. Silencing survivin splice variant 2B leads to antitumor activity in taxane--resistant ovarian cancer. Clin. Cancer Res. 2011, 17, 3716–3726. [Google Scholar] [CrossRef]
  19. Badran, A.; Yoshida, A.; Ishikawa, K.; Goi, T.; Yamaguchi, A.; Ueda, T.; Inuzuka, M. Identification of a novel splice variant of the human anti-apoptopsis gene survivin. Biochem. Biophys. Res. Commun. 2004, 314, 902–907. [Google Scholar] [CrossRef]
  20. Knauer, S.K.; Bier, C.; Schlag, P.; Fritzmann, J.; Dietmaier, W.; Rödel, F.; Klein-Hitpass, L.; Kovács, A.F.; Döring, C.; Hansmann, M.-L.; et al. The survivin isoform survivin-3B is cytoprotective and can function as a chromosomal passenger complex protein. Cell Cycle 2007, 6, 1502–1509. [Google Scholar] [CrossRef]
  21. Sawai, K.; Goi, T.; Hirono, Y.; Katayama, K.; Yamaguchi, A. Survivin-3B gene decreases the invasion-inhibitory effect of colon cancer cells with 5-fluorouracil. Oncol. Res. 2010, 18, 541–547. [Google Scholar] [CrossRef]
  22. Végran, F.; Boidot, R.; Bonnetain, F.; Cadouot, M.; Chevrier, S.; Lizard-Nacol, S. Apoptosis gene signature of Survivin and its splice variant expression in breast carcinoma. Endocr. Relat. Cancer 2011, 18, 783–792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Mahotka, C.; Liebmann, J.; Wenzel, M.; Suschek, C.V.; Schmitt, M.; Gabbert, H.E.; Gerharz, C.D. Differential subcellular localization of functionally divergent survivin splice variants. Cell Death Differ. 2002, 9, 1334–1342. [Google Scholar] [CrossRef] [PubMed]
  24. Lopergolo, A.; Tavecchio, M.; Lisanti, S.; Ghosh, J.C.; Dohi, T.; Faversani, A.; Vaira, V.; Bosari, S.; Tanigawa, N.; Delia, D.; et al. Chk2 phosphorylation of survivin-DeltaEx3 contributes to a DNA damage-sensing checkpoint in cancer. Cancer Res. 2012, 72, 3251–3259. [Google Scholar] [CrossRef] [PubMed]
  25. Vandghanooni, S.; Eskandani, M.; Montazeri, V.; Halimi, M.; Babaei, E.; Feizi, M.A.H. Survivin-deltaEx3: A novel biomarker for diagnosis of papillary thyroid carcinoma. J. Cancer Res. Ther. 2011, 7, 325–330. [Google Scholar] [PubMed]
  26. Altieri, D.C. Targeted therapy by disabling crossroad signaling networks: The survivin paradigm. Mol. Cancer Ther. 2006, 5, 478–482. [Google Scholar] [CrossRef] [PubMed]
  27. Stobiecka, M.; Chalupa, A.; Dworakowska, B. Piezometric biosensors for anti-apoptotic protein survivin based on buried positive-potential barrier and immobilized monoclonal antibodies. Biosens. Bioelectron. 2016, 84, 37–43. [Google Scholar] [CrossRef]
  28. Ratajczak, K.; Krazinski, B.E.; Kowalczyk, A.E.; Dworakowska, B.; Jakiela, S.; Stobiecka, M. Optical Biosensing System for the Detection of Survivin mRNA in Colorectal Cancer Cells Using a Graphene Oxide Carrier-Bound Oligonucleotide Molecular Beacon. Nanomaterials 2018, 8, 510. [Google Scholar] [CrossRef]
  29. Jang, J.S.; Kim, K.M.; Kang, K.H.; Choi, J.E.; Lee, W.K.; Kim, C.H.; Kang, Y.M.; Kam, S.; Kim, I.-S.; Jun, J.E.; et al. Polymorphisms in the survivin gene and the risk of lung cancer. Lung Cancer 2008, 60, 31–39. [Google Scholar] [CrossRef]
  30. Cheng, Z.-J.; Hu, L.-H.; Huang, S.-J. [Correlation of -31G/C polymorphisms of survivin promoter to tumorigenesis of gastric carcinoma]. Ai Zheng 2008, 27, 258–263. [Google Scholar]
  31. Kawata, N.; Tsuchiya, N.; Horikawa, Y.; Inoue, T.; Tsuruta, H.; Maita, S.; Satoh, S.; Mitobe, Y.; Narita, S.; Habuchi, T. Two survivin polymorphisms are cooperatively associated with bladder cancer susceptibility. Int. J. Cancer 2011, 129, 1872–1880. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Weng, C.J.; Hsieh, Y.H.; Chen, M.K.; Tsai, C.M.; Lin, C.W.; Yang, S.F. Survivin SNP-carcinogen interactions in oral cancer. J. Dent. Res. 2012, 91, 358–363. [Google Scholar] [CrossRef] [PubMed]
  33. Hsieh, Y.-S.; Tsai, C.-M.; Yeh, C.-B.; Yang, S.-F.; Hsieh, Y.-H.; Weng, C.-J. Survivin T9809C, an SNP located in 3’-UTR, displays a correlation with the risk and clinicopathological development of hepatocellular carcinoma. Ann. Surg. Oncol. 2012, 19, S625–S633. [Google Scholar] [CrossRef] [PubMed]
  34. Shi, H.; Bevier, M.; Johansson, R.; Enquist-Olsson, K.; Henriksson, R.; Hemminki, K.; Lenner, P.; Försti, A. Prognostic impact of polymorphisms in the MYBL2 interacting genes in breast cancer. Breast Cancer Res. Treat. 2012, 131, 1039–1047. [Google Scholar] [CrossRef] [PubMed]
  35. Cong, L.; Ran, F.A.; Cox, D.; Lin, S.; Barretto, R.; Habib, N.; Hsu, P.D.; Wu, X.; Jiang, W.; Marraffini, L.A.; et al. Multiplex Genome Engineering Using CRISPR/Cas Systems. Science 2013, 339, 819–823. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ozretić, P.; Trnski, D.; Musani, V.; Maurac, I.; Kalafatić, D.; Orešković, S.; Levanat, S.; Sabol, M. Non-canonical Hedgehog signaling activation in ovarian borderline tumors and ovarian carcinomas. Int. J. Oncol. 2017, 51, 1869–1877. [Google Scholar] [CrossRef] [PubMed]
  37. Hofstetter, G.; Berger, A.; Fiegl, H.; Slade, N.; Zorić, A.; Holzer, B.; Schuster, E.; Mobus, V.J.; Reimer, D.; Daxenbichler, G.; et al. Alternative splicing of p53 and p73: The novel p53 splice variant p53delta is an independent prognostic marker in ovarian cancer. Oncogene 2010, 29, 1997–2004. [Google Scholar] [CrossRef]
  38. McCall, M.N.; McMurray, H.R.; Land, H.; Almudevar, A. On non-detects in qPCR data. Bioinformatics 2014, 30, 2310–2316. [Google Scholar] [CrossRef] [Green Version]
  39. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3--new capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef]
  40. Cvok, M.L.; Cretnik, M.; Musani, V.; Ozretic, P.; Levanat, S. New sequence variants in BRCA1 and BRCA2 genes detected by high-resolution melting analysis in an elderly healthy female population in Croatia. Clin. Chem. Lab. Med. 2008, 46, 1376–1383. [Google Scholar] [CrossRef]
  41. Roessler, E.; Ermilov, A.N.; Grange, D.K.; Wang, A.; Grachtchouk, M.; Dlugosz, A.A.; Muenke, M. A previously unidentified amino-terminal domain regulates transcriptional activity of wild-type and disease-associated human GLI2. Hum. Mol. Genet. 2005, 14, 2181–2188. [Google Scholar] [CrossRef] [Green Version]
  42. Dai, J.; Jin, G.; Dong, J.; Chen, Y.; Xu, L.; Hu, Z.; Shen, H. Prognostic significance of survivin polymorphisms on non-small cell lung cancer survival. J. Thorac. Oncol. 2010, 5, 1748–1754. [Google Scholar] [CrossRef] [PubMed]
  43. Aynaci, E.; Coskunpinar, E.; Eren, A.; Kum, O.; Oltulu, Y.M.; Akkaya, N.; Turna, A.; Yaylim, I.; Yildiz, P. Association between survivin gene promoter -31G/C and -644C/T polymorphisms and non-small cell lung cancer. Genet. Mol. Res. 2013, 12, 3975–3982. [Google Scholar] [CrossRef] [PubMed]
  44. Yang, X.; Xiong, G.; Chen, X.; Xu, X.; Wang, K.; Fu, Y.; Yang, K.; Bai, Y. Polymorphisms of survivin promoter are associated with risk of esophageal squamous cell carcinoma. J. Cancer Res. Clin. Oncol. 2009, 135, 1341–1349. [Google Scholar] [CrossRef] [PubMed]
  45. Hmeljak, J.; Erčulj, N.; Dolžan, V.; Kern, I.; Cör, A. BIRC5 promoter SNPs do not affect nuclear survivin expression and survival of malignant pleural mesothelioma patients. J. Cancer Res. Clin. Oncol. 2011, 137, 1641–1651. [Google Scholar] [CrossRef]
  46. Boidot, R.; Vegran, F.; Jacob, D.; Chevrier, S.; Gangneux, N.; Taboureau, J.; Oudin, C.; Rainville, V.; Mercier, L.; Lizard-Nacol, S. The expression of BIRC5 is correlated with loss of specific chromosomal regions in breast carcinomas. Genes Chromosomes Cancer 2008, 47, 299–308. [Google Scholar] [CrossRef]
  47. Pu, F.; Shao, Z.; Yang, S.; Liu, J.; Lin, S.; Ma, X.; Yang, H. Association between functional variants in BIRC5/survivin gene 3’ untranslated region and mRNA expression in lymphoblastoid cell lines. Oncol. Lett. 2015, 10, 2319–2322. [Google Scholar] [CrossRef] [PubMed]
  48. Li, X.-B.; Li, S.-N.; Yang, Z.-H.; Cao, L.; Duan, F.-L.; Sun, X.-W. Polymorphisms of survivin and its protein expression are associated with colorectal cancer susceptibility in Chinese population. DNA Cell Biol. 2013, 32, 236–242. [Google Scholar] [CrossRef]
  49. Rosato, A.; Menin, C.; Boldrin, D.; Dalla Santa, S.; Bonaldi, L.; Scaini, M.C.; Del Bianco, P.; Zardo, D.; Fassan, M.; Cappellesso, R.; et al. Survivin expression impacts prognostically on NSCLC but not SCLC. Lung Cancer 2013, 79, 180–186. [Google Scholar] [CrossRef]
  50. Kafadar, D.; Yaylim, I.; Kafadar, A.M.; Cacina, C.; Ergen, A.; Kaynar, M.Y.; Isbir, T. Investigation of Survivin Gene Polymorphism and Serum Survivin Levels in Patients with Brain Tumors. Anticancer Res. 2018, 38, 5991–5998. [Google Scholar] [CrossRef]
  51. Zu, Y.; Ban, J.; Xia, Z.; Wang, J.; Cai, Y.; Ping, W.; Sun, W. Genetic variation in a miR-335 binding site in BIRC5 alters susceptibility to lung cancer in Chinese Han populations. Biochem. Biophys. Res. Commun. 2013, 430, 529–534. [Google Scholar] [CrossRef] [PubMed]
Cells 08 00128 g001 550
Figure 1. The exon structure of five splice isoforms of survivin. Exons are drawn relative to their size. Vertical bars indicate the site of stop codon of each isoform.
Figure 1. The exon structure of five splice isoforms of survivin. Exons are drawn relative to their size. Vertical bars indicate the site of stop codon of each isoform.
Cells 08 00128 g001
Cells 08 00128 g002 550
Figure 2. Levels of survivin isoform expression in the SKOV-3 GLI KO lines compared with SKOV-3 WT cells. * represents statistically significant downregulation in comparison with the WT cell line (p < 0.05).
Figure 2. Levels of survivin isoform expression in the SKOV-3 GLI KO lines compared with SKOV-3 WT cells. * represents statistically significant downregulation in comparison with the WT cell line (p < 0.05).
Cells 08 00128 g002
Cells 08 00128 g003 550
Figure 3. Survivin isoform expression levels in the SKOV-3 WT cell line and the SKOV-3 GLI1 KO cell line after treatment with the GLI1/2 inhibitor GANT-61 for 24h. NT, non-treated samples; G 10µM, cells treated with 10 µM GANT-61; G 20µM, cells treated with 20 µM GANT-61. * represents statistically significant downregulation (p < 0.05).
Figure 3. Survivin isoform expression levels in the SKOV-3 WT cell line and the SKOV-3 GLI1 KO cell line after treatment with the GLI1/2 inhibitor GANT-61 for 24h. NT, non-treated samples; G 10µM, cells treated with 10 µM GANT-61; G 20µM, cells treated with 20 µM GANT-61. * represents statistically significant downregulation (p < 0.05).
Cells 08 00128 g003
Cells 08 00128 g004 550
Figure 4. The effect of GLI1, GLI2, GLI3, and GLI3R overexpression in the SKOV-3 WT cell line on (A) the expression levels of GLI1, GLI2, GLI3, and PTCH1 and (B) the expression levels of survivin isoforms. * represents statistically significant changes in expression compared with cells transfected with empty vector.
Figure 4. The effect of GLI1, GLI2, GLI3, and GLI3R overexpression in the SKOV-3 WT cell line on (A) the expression levels of GLI1, GLI2, GLI3, and PTCH1 and (B) the expression levels of survivin isoforms. * represents statistically significant changes in expression compared with cells transfected with empty vector.
Cells 08 00128 g004
Cells 08 00128 g005 550
Figure 5. Relative expression of survivin isoforms in the SKOV-3 cell line (circles) and ovarian carcinoma (OC) tissue samples (triangles).
Figure 5. Relative expression of survivin isoforms in the SKOV-3 cell line (circles) and ovarian carcinoma (OC) tissue samples (triangles).
Cells 08 00128 g005
Cells 08 00128 g006 550
Figure 6. Relative expression of S WT (A) and S 2α (B) in fallopian tube (FT) and ovarian carcinoma (OC) samples.
Figure 6. Relative expression of S WT (A) and S 2α (B) in fallopian tube (FT) and ovarian carcinoma (OC) samples.
Cells 08 00128 g006
Cells 08 00128 g007 550
Figure 7. Pairwise linkage disequilibrium (LD) of eleven BIRC5 sequence variants with highest minor allele frequency. The location of each sequence variant along the BIRC5 gene is relative to the real nucleotide position. The number in each diamond indicates the intensity of LD (r2 × 102) between respective pairs of sequence variants. The LD strength is also represented by shades of gray (0 [white] < r2 < 1 [black]). (A) Controls. (B) Cases.
Figure 7. Pairwise linkage disequilibrium (LD) of eleven BIRC5 sequence variants with highest minor allele frequency. The location of each sequence variant along the BIRC5 gene is relative to the real nucleotide position. The number in each diamond indicates the intensity of LD (r2 × 102) between respective pairs of sequence variants. The LD strength is also represented by shades of gray (0 [white] < r2 < 1 [black]). (A) Controls. (B) Cases.
Cells 08 00128 g007
Table
Table 1. The association of survivin isoform expression in OC samples.
Table 1. The association of survivin isoform expression in OC samples.
S 2αS 2BS 3BS Δex3--
0.920.690.810.75S WTϱ
<0.00010.002<0.0001<0.0001p
-0.670.800.81S 2αϱ
-0.005<0.0001<0.0001p
--0.550.67S 2Bϱ
--0.0210.004p
---0.66S 3Bϱ
---0.002p
------
Spearman rank correlation coefficient (ϱ)
Table
Table 2. BIRC5 polymorphisms in OC samples and healthy controls. MAF, minor allele frequency; OC, ovarian carcinoma.
Table 2. BIRC5 polymorphisms in OC samples and healthy controls. MAF, minor allele frequency; OC, ovarian carcinoma.
Gene RegionSNP ID NumberNucleotide ChangeMAF in OC (%)MAF in Controls (%)p-Value
promoterrs3764383c.-1547C>T *23.125.00.871
promoterrs143396310c.-1458C>T1.31.41.000
promoterrs8073903c.-644T>C26.333.10.297
promoterrs8073069c.-625G>C21.322.31.000
promoterrs17878731c.-267G>A1.30.71.000
promoterrs17878467c.-241C>T5.010.80.219
promoterrs17887126c.-235G>A6.31.40.053
5′UTRrs9904341c.-31G>C48.837.20.093
intron 2rs4789551c.221+209T>C6.44.70.756
exon 4rs2071214c.9194G>A **3.83.41.000
3′UTRrs17885521c.9288G>C1.32.01.000
3′UTRrs17882627c.9342G>A2.51.40.614
3′UTRrs2239680c.9386T>C23.823.01.000
3′UTRrs1042489c.9809T>C25.035.80.104
3′UTRrs2661694c.10611C>A26.325.71.000
* C is the minor allele, ** G is the minor allele.

Share and Cite

MDPI and ACS Style

Trnski, D.; Gregorić, M.; Levanat, S.; Ozretić, P.; Rinčić, N.; Vidaković, T.M.; Kalafatić, D.; Maurac, I.; Orešković, S.; Sabol, M.; et al. Regulation of Survivin Isoform Expression by GLI Proteins in Ovarian Cancer. Cells 2019, 8, 128. https://doi.org/10.3390/cells8020128

AMA Style

Trnski D, Gregorić M, Levanat S, Ozretić P, Rinčić N, Vidaković TM, Kalafatić D, Maurac I, Orešković S, Sabol M, et al. Regulation of Survivin Isoform Expression by GLI Proteins in Ovarian Cancer. Cells. 2019; 8(2):128. https://doi.org/10.3390/cells8020128

Chicago/Turabian Style

Trnski, Diana, Maja Gregorić, Sonja Levanat, Petar Ozretić, Nikolina Rinčić, Tajana Majić Vidaković, Držislav Kalafatić, Ivana Maurac, Slavko Orešković, Maja Sabol, and et al. 2019. "Regulation of Survivin Isoform Expression by GLI Proteins in Ovarian Cancer" Cells 8, no. 2: 128. https://doi.org/10.3390/cells8020128

APA Style

Trnski, D., Gregorić, M., Levanat, S., Ozretić, P., Rinčić, N., Vidaković, T. M., Kalafatić, D., Maurac, I., Orešković, S., Sabol, M., & Musani, V. (2019). Regulation of Survivin Isoform Expression by GLI Proteins in Ovarian Cancer. Cells, 8(2), 128. https://doi.org/10.3390/cells8020128

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop