1. Introduction
Rhabdomyosarcoma (RMS), the most prevalent soft tissue malignancy in children and adolescents, arises from abnormal differentiation of skeletal muscle precursor cells [
1]. Traditionally classified on morphological features, RMS comprises three primary subtypes (embryonal, alveolar, and pleomorphic), driven by distinct molecular mechanisms, presenting unique challenges for therapeutic interventions. Hence, recent advancements in molecular biology and genetics have facilitated improvements in RMS diagnostics, classification, patient risk stratification, and management strategies [
2]. Moreover, elucidating the molecular biology of RMS holds promise for identifying novel therapeutic targets [
3].
The diagnostics of RMS relies on immunohistochemical detection of muscle markers like desmin, myoglobin, myogenin, and MyoD. However, there is a need for new markers as the conventional ones are not always sensitive and specific enough to ensure an accurate diagnosis [
4]. Among potential muscle-specific markers, the ankyrin repeat domain 1 (ANKRD1) has emerged for differential diagnostics of RMS, due to its higher expression in RMS compared to other sarcomas [
5]. Despite that, its role and regulation in RMS remain unclear.
ANKRD1, (also known as cardiac ankyrin repeat protein, CARP), a muscle ankyrin repeat protein (MARP) family member, is a highly conserved, multifunctional protein [
6]. This mechanical stress-inducible protein is involved in cardiac and skeletal muscle remodeling [
7]. It participates in cardiogenesis, mechanosensing, regulation of gene expression, and intracellular signaling [6,8-11]. ANKRD1 is mostly expressed in cardiomyocytes, both in the nucleus (role of transcriptional regulator) and at the sarcomeric I-band (role in sarcomere organization). It is able to translocate from the cytoplasm to the nucleus in response to mechanical stimuli, contributing to the transformation of biomechanical stress signals to the regulation of gene expression.
Changes in ANKRD1 expression have been observed in various tumors, including RMS and ovarian carcinoma [
12,
13]. Altered ANKRD1 expression has prognostic implications in ovarian and renal clear cell carcinoma and has been associated with chemotherapy resistance in ovarian and lung carcinoma [14-16]. Furthermore, ANKRD1 inhibition sensitizes ovarian cancer cells to endoplasmic reticulum stress-induced apoptosis [
13], while ANKRD1 upregulation promotes osteosarcoma cell proliferation and invasion [
17]. In contrast, ANKRD1 upregulation suppresses colony formation in hepatocellular, lung, colon, and prostate cancer cell lines [
18,
19]. The YAP/TAZ oncogene complex and tumor suppressor p53 are involved in the transcriptional regulation of ANKRD1 expression [
20,
21,
22].
Considering the significance of ANKRD1 in carcinogenesis, its association with chemotherapy resistance in various tumor types, and its role in normal muscle cell differentiation, which is disrupted in RMS, we aimed to further investigate this potential diagnostic marker by examining its expression at the mRNA and protein level, as well as its intracellular localization and protein turnover, in cell lines representing three main subtypes of RMS.
2. Materials and Methods
2.1. Cell culture, treatments, and transfection
Human rhabdomyosarcoma cell lines: RD (ATCC CCL-136), HS-729 (ATCC HTB-153) and SJRH30 (ATCC CRL-2061), and HeLa (ATCC CCL-2) cells were used in this study. Cells were cultured in Dulbecco’s Modified Eagle’s Medium - high glucose (Sigma-Aldrich, USA) supplemented with 10% of heat-inactivated fetal bovine serum (Gibco, Thermo Fisher Scientific, USA) and antibiotic-antimycotic mixture (Sigma-Aldrich, USA). Cells were maintained in a humidified environment with 5% CO2 at 37°C and subcultured twice a week. Cellular treatments included incubation with 10 µM proteasome inhibitor, MG132 (Merck KGaA, Germany) for 4 and 8 h. Transient transfections were performed at 90% confluence, with the pCMV-Tag2B expression vector containing ANKRD1 or ANKRD2 open reading frame, using Lipofectamine 2000 Reagent (Invitrogen, Thermo Fisher Scientific, USA) according to the manufacturer’s instructions.
2.2. RT-PCR and quantitative PCR
Total RNA was extracted from cells using the TRIzol reagent (Thermo Fisher Scientific, USA), according to the manufacturer’s protocol. RNA concentration and purity were determined by spectrophotometry using the BioSpec-nano (Shimadzu, Japan). Complementary DNA (cDNA) was synthesized from 2 μg of total RNA using random hexamer primers and High Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific, USA), according to the manufacturer’s instructions.
ANKRD1 and
GAPDH fragments were amplified using Native Taq DNA polymerase (EURx, Poland), cDNA as a template, and primers, listed in
Table 1. The quantity of cDNA in each reaction corresponded to 100 ng of total RNA. The PCR reaction was performed in the 2720 Thermal Cycler (Applied Biosystems, USA).
GAPDH was used as a cDNA quality control.
Quantitative PCR reactions were carried out in technical triplicate for each sample and for three biological replicates of each cell line on a 7900HT Fast Real-Time System (Applied Biosystems, USA) using TaqMan Gene Expression Assay (Applied Biosystems, USA). The GAPDH transcript was used as an internal reference to normalize the mRNA levels in different samples. The data were analyzed using the 2−∆Ct method. The results are presented as mean ± standard deviation.
2.3. Sequencing of the ANKRD1 open reading frame
ANKRD1 open reading frame was amplified using the cDNA templates, Phusion High-Fidelity PCR Master Mix (Thermo Fisher Scientific, USA), and primers given in
Table 1. The PCR reaction was performed in the Mastercycler gradient (Eppendorf, Germany). Generated fragments were sequenced using the BigDye
TM Terminator Version 3.1 Ready Reaction Kit (Applied Biosystems, USA) according to the manufacturer’s protocol, on a 3130 Genetic Analyzer (Applied Biosystems, USA) with primers listed in
Table 1.
2.4. Western blot
Total cellular proteins were extracted as previously described [
23], separated by SDS-PAGE, and blotted onto the PVDF membrane (Immobilon P, Millipore, USA). The anti-ANKRD1 antibody against full-length ANKRD1, generated by the procedure reported in Kojic et al. [
24] was used at 1:750 dilution; rabbit polyclonal anti-GAPDH antibody (Guangzhou Dingguo Biology, China) was used at 1:2000 dilution; mouse monoclonal anti-FLAG antibody (Sigma-Aldrich, USA), and rabbit polyclonal anti-p21 antibody (Santa Cruz Biotechnology, USA) were used at 1:1000 dilution.
2.5. Immunoprecipitation
For ANKRD1 protein enrichment, RMS cell extracts containing 0.5 or 1 mg of total proteins were incubated overnight at 4°C with anti-ANKRD1 antibody. The following day, the antibody-antigen complex was pulled out of the mixture using the Pierce Protein A/G UltraLink Resin (Thermo Fisher Scientific, USA). Resins were then washed with phosphate-buffered saline (PBS), resuspended in 30 µl of coIP buffer (125mM Tris-HCl pH 6.8, 4% SDS, 20% glycerol, 100mM DTT, bromophenol blue), incubated for 5 min at 95°C and centrifuged for 1 min at 2000 rpm. The supernatants were analyzed by Western blot.
2.6. Immunocytochemistry and microscopy
After seeding on coverslips and reaching 80% confluence, cells were fixed in 4% paraformaldehyde, permeabilized with 0.5% TritonX in PBS, and incubated overnight at 4°C in blocking solution (3% bovine serum albumin in PBS). Coverslips were then incubated with primary antibodies, diluted in 1% bovine serum albumin in PBS, for 3 h at room temperature, washed with PBS, and incubated with secondary antibodies labeled with Alexa Fluor 488 or Alexa Fluor 568 (Invitrogen, Thermo Fisher Scientific, USA) for 1 h. Primary antibodies were: anti-ANKRD1 (mouse, 1:100 dilution), anti-coilin (rabbit, 1:300 dilution, Cell Signaling Technology, USA), and anti-PML (goat, 1:100 dilution, Santa Cruz Biotechnology, USA). Nuclei were stained with DAPI (100 ng/ml) before mounting with Dako fluorescence medium (Agilent Technologies, Denmark). Confocal imaging and image manipulation were performed using the Leica SP8 system and Leica Application Suite X (Leica, Germany).
2.7. Statistical analysis
Statistical analysis was conducted using SPSS software v.21.0 (IBM SPSS Statistics, USA). The normal distribution and variance homogeneity of the data were verified by the Shapiro–Wilk and the Levene’s test, respectively. One-way ANOVA was employed for the comparison of the ANKRD1 transcript levels between the RMS cell lines. Results are presented as mean±SD, and p-value ≤ 0.05 was considered significant.
4. Discussion
ANKRD1 had been essentially investigated in muscle cells, where numerous reports had suggested its involvement in muscle differentiation under physiological and stress conditions, as well as in mechanosignaling pathways in mature muscle cells [6,8-10]. Parallel studies have reported an increased ANKRD1 expression in cells and tissues originating from certain primary tumors, such as rhabdomyosarcoma [
10] and ovarian carcinoma [
14]. Transcriptomic data from the Human Protein Atlas states increased ANKRD1 transcript level in cell lines derived from several tumor types, including rhabdomyosarcoma, testis, liver, ovarian, and pancreatic cancer. Additionally, ANKRD1 was suggested to play a role in cancer development since it is implicated in the progression of cancers, including those affecting ovaries, breast, and pancreas [
13,
26,
27], and in the chemoresistance of ovarian and lung cancer [
14,
15]. In RMS, ANKRD1 was previously investigated as a potential diagnostic marker. Ishiguro et al. [
5] reported upregulation of ANKRD1 in RMS compared to other sarcoma types, indicating its potential utility in differential diagnostics.
Driven by many findings suggesting the role of ANKRD1 in molecular mechanisms of tumorigenesis, we performed this in vitro study to examine its expression profile and regulation in cell lines derived from three RMS subtypes: RD (embryonal), SJRH30 (alveolar), and HS-729 (pleomorphic).
In RD and SJRH30 cell lines, we verified the ANKRD1 transcript expression as reported in Human Protein Atlas, and additionally detected a similar ANKRD1 mRNA level in HS-729 cells. In contrast to the similar mRNA levels, the ANKRD1 protein exhibited differential expression among these cell lines, suggesting different regulatory mechanisms controlling its protein level in RMS subtypes. Western blot analysis revealed detectable ANKRD1 protein levels in SJRH30 and RD cells (SJRH30 > RD), while in HS-729 cells the corresponding signal was not observed, even after immunoprecipitation. However, ANKRD1 was detected in HS-729 cells by immunofluorescence. This could be attributed to the differences in sensitivity or antibody affinities in two methods.
The fast protein turnover rate by proteasomal degradation could explain the discrepancy between ANKRD1 mRNA and protein levels and different amounts of endogenous ANKRD1 detected in protein extracts derived from different batches of cells, since proteasomal inhibition led to an increase in the ANKRD1 protein level in RMS cell lines. Similarly to other proteins containing PEST sequences that serve as proteolytic signals, ANKRD1 possesses a short intracellular half-life due to rapid degradation in 26S proteasomes [
28,
29]. This mechanism could also contribute to our inability to achieve ANKRD1 overexpression in RMS cells. However, in line with previous evidence [
30], the overexpression of ANKRD2, a closely related protein sharing structural motifs with ANKRD1, including PEST sequences, was achieved. Additionally, we attained ANKRD1 overexpression in HeLa cells, despite its short half-life due to the proteasomal degradation in this cell line [
28]. Our results imply that resistance to excessive ANKRD1 levels is a specific feature of RMS cells and mechanisms for regulation of its protein level in RMS cells warrant further exploration. For example, it would be worth investigating if ANKRD1 is implicated in growth inhibition of RMS cells triggered by the loss of proteasomal activity induced by bortezomib treatment [
31,
32].
Similar cases of cell lines exhibiting intolerance to gene overexpression have been reported. Lee et al. [
33] demonstrated that mouse embryonic cell lines do not tolerate SUMO1 overexpression, adapting by reducing SUMO1 levels. Analogous observations include findings on glucocorticoid receptor expression in small-cell lung cancer cell lines [
34] and research on SIRT6 overexpression in cancer cell lines [
35]. The list goes on, including many pro-apoptotic and tumor suppressor proteins, like Bax, caspase-9, Beclin1, and p21 [
36,
37,
38,
39]. Additionally, the role of ANKRD1 in apoptosis modulation in various cancer types with its differential pro- and anti-apoptotic effects [,
13,
17,
18,
19,
27], further confounds the understanding of its function in RMS cells. The function of ANKRD1 in tumor cells, particularly in apoptosis, should be considered in regard to p53. In addition to their interaction at protein level, we have previously shown that ANKRD1 upregulates p53 activity [
20]. Furthermore, overexpressed ANKRD1 triggers the transcriptional activation of p53, inducing cellular apoptosis in cardiomyocytes. [
40]. Our speculation that overexpression of ANKRD1 might be involved in induction of apoptotic death of RMS cells via activation of p53 remains to be clarified.
In muscle cells, ANKRD1 has dual, cytoplasmic and nuclear localization. It participates in sarcomere mechanosensing by translocation to the nucleus upon mechanical stimuli where it participates in the regulation of stress response genes [
6,
11]. In RMS cells we also observed both cytoplasmic and nuclear localization of ANKRD1 in RMS cells, which is in alignment with the findings of Ishiguro et al. [
10]. However, inconsistently with their results, we detected higher ANKRD1 abundance in the nucleus compared to the cytoplasm. This could arise from differences in sample types (tissue sections vs. cell lines) or in the use of different antibodies (with different epitopes and sensitivity). Nuclear localization of ANKRD1 in RMS and RMS cell lines suggests its regulatory role. Accordingly, the next step would be to identify transcriptional targets downstream of ANKRD1 specific to RMS.
We detected ANKRD1 in structures resembling nuclear bodies in nuclei of RD and HS-729 cells. To identify these structures, the first hint was to explore the localization of ANKRD1 in PML bodies, since another MARP family member ANKRD2 is known to co-localize with PML protein in the PML bodies of human myoblasts [
41]. However, there was no co-localization of ANKRD1 with a marker of PML bodies. Then we tested if ANKRD1 is localized to Cajal bodies. There was co-localization with coilin, resident protein of Cajal bodies, membraneless organelles implicated in the assembly of spliceosomal snRNPs [
42], and in the maintenance of telomeres notably in cancer cells [
43]. Although having many known protein partners, ANKRD1’s interactions with Cajal bodies-resident proteins like fibrillarin, SMN, and coilin have not been reported. Furthermore, there is no evidence that ANKRD1 could bind to RNA molecules. However, the reported interaction of ANKRD1 with YB-1 [
44], an RNA-binding protein suggested as a target for RMS therapy [
45], introduces intriguing prospects. Although rare, YB-1's presence in Cajal bodies [
46] and its role as a component of spliceosomes [
47,
48] further emphasize the need for further research into the interactions of ANKRD1 with YB-1 and related proteins in Cajal bodies. Another possible direction would be to investigate if ANKRD1, as a downstream target of the Hippo pathway [
19] involved in the responses to oxidative stress [
49], has a role in the redox homeostasis driven by Cajal bodies in cancer cells [
50].
Although the presence of ANKRD1 within the nucleoli fibrillar center is documented in various cell lines (proteinatlas.org), our study does not report ANKRD1 in nuclear structures other than Cajal bodies. Given the dynamic relationship between Cajal bodies and nucleoli in ribosome biogenesis, telomerase localization, RNA interference, and stress responses [
51], translocation of ANKRD1 between Cajal bodies and nucleoli might present a topic for further investigation.