Depletion of the m1A writer TRMT6/TRMT61A reduces proliferation and resistance against cellular stress in bladder cancer

Background Bladder cancer (BLCA) is a common and deadly disease that results in a reduced quality of life for the patients and a significant economic burden on society. A better understanding of tumorigenesis is needed to improve clinical outcomes. Recent evidence places the RNA modification m1A and its regulatory proteins TRMT6/TRMT61A and ALKBH3 in BLCA pathogenesis. Methods TRMT6/TRMT61A, ALKBH1, and ALKBH3 expression was examined in human BLCA cell lines and a normal urinary tract epithelium cell line through qRT-PCR and western blot analysis. Prestoblue Cell Viability Reagent, wound-healing assay, and live-cell imaging-based cell displacement analysis, were conducted to assess proliferation, migration, and displacement of this BLCA cell line panel. Cell survival was assessed after inducing cellular stress and activating the unfolded protein response (UPR) with tunicamycin. Moreover, siRNA-mediated gene silencing in two BLCA cell lines (5637 and HT1197) was conducted to investigate the biological roles of TRMT6/TRMT61A. Results Heterogeneous morphology, proliferation, displacement, tunicamycin sensitivity, and expression levels of m1A regulators were observed among the panel of cell lines examined. In general, TRMT61A expression was increased in BLCA cell lines when compared to SV-HUC-1. Depletion of TRMT6/TRMT61A reduced proliferation capacity in both 5637 and HT1197 cell lines. The average cell displacement of 5637 was also reduced upon TRMT6/TRMT61A depletion. Interestingly, TRMT6/TRMT61A depletion decreased mRNA expression of targets associated with the ATF6-branch of the UPR in 5637 but not in HT1197. Moreover, cell survival after induction of cellular stress was compromised after TRMT6/TRMT61A knockdown in 5637 but not in HT1197 cells. Conclusion The findings suggest that TRMT6/TRMT61A plays an oncogenic role in BLCA and is involved in desensitizing BLCA cells against cellular stress. Further investigation into the regulation of TRMT6/TRMT61A expression and its impact on cellular stress tolerance may provide insights for future BLCA treatment.


Introduction
Bladder cancer (BLCA) is the 10 th most frequently diagnosed cancer worldwide, accounting for 573,000 new annual diagnoses and 213,000 deaths.Due to gender discrepancy, BLCA is the 6 th most common cancer in men and the 9 th leading cause of cancer death (1).Meanwhile, there is an increasing incidence in women.Thus, BLCA poses a significant public health concern at a global level.
The most frequent type of BLCA is non-muscle invasive bladder cancer (NMIBC), accounting for about 75% of newly diagnosed cases (2).Despite being characterized by a favorable prognosis and overall survival, relapses are frequent, leading to lifelong control schemes, costly treatment, and poor quality of life.Patients with muscle-invasive bladder cancer (MIBC) face poor prognoses and challenges related to major surgery and urinary diversions as well as chemo-and radiotherapy resistance and risk of metastasis (3).The treatment scheme for NMIBC is transurethral resection of the bladder tumor (TURBT), often combined with intravesical instillation of Bacillus Calmette-Gueŕin (BCG).For MIBC, radical cystectomy combined with urinary diversion is currently the primary treatment strategy (4).However, treatment strategies for BLCA are advancing.For NMIBC, this includes strategies to improve the delivery and uptake of current therapeutic agents and the development of new immunotherapies.For MIBC, immune checkpoint inhibitors (ICIs) are promising candidates in maintenance therapy.For more details on advances in BLCA therapeutic strategies, the reader is referred to excellent reviews (5,6).
In BLCA, increased hTrm6p/hTrm61p expression and elevated m1A levels promote disease progression (16).Recently, we reported increased TRMT6/TRMT61A expression accompanied by elevated m1A levels on a tRF class in BLCA and suggested that this promotes the unfolded protein response (UPR), most likely through regulation of m1A base modification on these tRFs (17).UPR is a cellular response activated when cells face an imbalance in the cellular environment, resulting in the accumulation of misfolded or unfolded proteins due to intrinsic or extrinsic stressors.Emerging evidence suggests that by enabling cancer cells to cope with surrounding stress stimuli, the UPR is a central player in malignant transformation and tumorigenesis by regulating cancer cell survival, angiogenesis, metastasis, and chemoresistance (18).Moreover, elevated ALKBH3 expression is found in several cancers, including BLCA (19).Thus, evidence suggests that dysregulation of m1A regulators is involved in BLCA tumorigenesis, but the underlying mechanisms remain elusive.
In this study, a panel of BLCA cell lines and a bladder epithelial control cell line were characterized with regard to expression levels of the m1A methyltransferase complex TRMT6/TRMT61A and demethylases ALKBH1 and ALKBH3.The biological role of TRMT6/TRMT61A on proliferation, migration, and displacement was investigated in two BLCA cell lines.Moreover, the effect of TRMT6/TRMT61A on UPR activation and stress tolerance was further investigated in these two BLCA cell lines.
2 Methods and materials

Cell proliferation and survival assay
PrestoBlue ™ Cell Viability Reagent (A13262, Invitrogen) was utilized for cell viability, and this assay was carried out according to the manufacturer's instructions.Briefly, cells were seeded onto 96well plates (VWR/Falcon, 353072).Cell proliferation and survival (refer to Section 2.4 Tunicamycin treatment) were assessed by adding 10 µL Prestoblue (Invitrogen, A13262) to each well containing 100 µL medium.After incubating for 10 minutes (15) at 37°C, fluorescence was measured on a plate reader (VICTOR Nivo Multimode plate reader, Perkin Elmer).Between three and five independent experiments were carried out to measure proliferation and survival after tunicamycin treatment.

Tunicamycin treatment
For cell survival after tunicamycin (Merck/Sigma-Aldrich, #T7765) treatment, approximately 6000 cells were plated in 96 well plates.After 4 hours (h), tunicamycin was added to the wells at different concentrations (100, 200, 400, 600 and 800 nM), and cells were allowed to grow for 76 h.Cell viability was analyzed at four timepoints: 4 h after seeding (at the time tunicamycin was added, T=0) and 24, 48 and 72 h (T=24, T=48 and T=72) after tunicamycin addition.To induce ER stress and subsequently activate the UPR, cells were treated with 5 µg/mL tunicamycin for 8 h.

RNA isolation
Total RNA of bladder cancer cell lines was extracted using RNAzol® RT reagent (Molecular Research Center, #MR-RN190-500) according to the manufacturer's instructions.RNA quality and quantity was assessed with Nanodrop (Thermo Scientific ™ NanoDrop ™ One Microvolume UV-Vis Spectrophotometer).

Reverse transcription-quantitative PCR (RT-qPCR assays)
cDNA was obtained using SuperScript ™ IV VILO ™ Master Mix with ezDNase ™ Enzyme (Invitrogen, #11766050) according to manufacturer's instructions.500 ng total RNA (refer to Section 2.5 RNA isolation) was used per cDNA reaction.For mRNA analysis, qRT-PCR was performed with PowerUp ™ SYBR ™ Green Master Mix (Applied Biosystems; Thermo Fisher Scientific, USA, #A25742) using the StepOne Plus and QuantStudio ™ 5 Real-Time PCR Systems (Applied Biosystems; Thermo Fisher Scientific, USA, #A28140).4 ng cDNA was used in 10 µL reaction.All sequences for primers used are listed in Supplementary Table 1.Fold changes in mRNA expression were determined using the 2 -DDCT method and normalized against GAPDH or ACTB as the endogenous control and further normalized relative to SV-HUC1 with StepOne ™ Software v2.3.

Cell migration
For morphology analysis, 5 x 10 5 cells were seeded in 6 well plates, and pictures were taken after 24-36 hours with an Olympus CKX53 inverted microscope.For the scratch wound healing assay, 2 x 10 6 cells were seeded in 24-well plates.After 24 h, a pipette tip (1000 ml) was used to create a scratch wound.Next, the cell cultures were washed using PBS and incubated at 37C and 5% CO 2 .The wounds were observed under an Olympus CKX53 inverted microscope for 24 h in situ to compare the number of cells that migrated into the wound areas.

Analysis of cell displacements
Urinary bladder cancer cell lines (SV-HUC-1, 5637, SW780, T24, 5HT1376, HT1197) were seeded at sub-confluent densities (4000 cells/well) in 96-well glass bottom plates (Greiner Sensoplate, M4187-16EA, Merck).After seeding, the plates were placed in a CO 2 incubator overnight to allow cells to attach to the glass surface.Subsequently, cells were monitored using an ImageXpress Micro Confocal high-content microscope (Molecular Devices), equipped with an environmental chamber maintaining 5% CO 2 and 37°C, and controlled by the MetaXpress 6 software.Time lapse series of phase contrast images were acquired in widefield mode using a 20x 0.45 NA Ph1 air objective, camera binning = 2, a time interval of 3 min between frames and a total imaging period of 20 h.Cell displacement were analyzed by particle tracking using the TrackMate plugin in Fiji ImageJ (20) in combination with in house Python-based script (Python 3.7.6).Average displacement speed was calculated as the average of all tracked cell displacements within a 12 h period of imaging.

Cell displacement analysis after siRNA knockdown
5637 cells treated with siRNAs were seeded at a density of 50000 cells/well in 96-well glass bottom plates (Greiner Sensoplate, M4187-16EA, Merck) and placed in a CO 2 incubator overnight.Prior to imaging, cells were stained with NucBlue solution (Invitrogen) to label the nuclei of live cells.In brief, cells were incubated at 37°C and 5% CO 2 for 30 min with NucBlue solution diluted 1:40 in culture medium.Subsequently, cells were washed twice with pre-warmed culture medium.Image acquisition was carried out in an ImageXpress Micro Confocal high-content microscope (Molecular Devices) in widefield mode.Time lapse series were acquired for a total time period of 4.5 h using a time interval of 4 min between frames, a 20x 0.45 NA Ph1 air objective, filter set for detection of DAPI and phase contrast, and an environmental control gasket that maintain 37°C and 5% CO 2 .Cell displacement was analyzed by particle tracking using the TrackMate plugin in Fiji ImageJ (20) in combination with inhouse Python-based script (Python 3.7.6).Average displacement speed was calculated as the average of all tracked cell displacements within a 4 h period of imaging.

Quantification of RNA modifications using LC-MS/MS
Cells were harvested and stored in RNAlater ™ solution (Invitrogen, #AM7021) prior to total RNA isolation (refer to Section 2.5 RNA isolation).tRNAs were enriched using an Agilent 1260 Infinity II Analytical-Scale LC-UV Purification System with a Bio SEC-3 300 Å, 2.1 x 300 mm column (Agilent Technologies) chromatographed isocratically with 100 mM ammonium acetate pH 7 at 0.280 ml/min and 40°C for 20 min.Chromatograms were recorded at 260 nm and tRNAs were collected and lyophilized and solved in 30 ml of water.The RNA was enzymatically hydrolyzed to ribonucleosides by 20 U benzonase (Santa Cruz Biotech) and 0.2 U nuclease P1 (Sigma) in 10 mM ammonium acetate pH 6.0 and 1 mM magnesium chloride at 40°C for 1 h, then added ammonium bicarbonate to 50 mM, 0.002 U phoshodiesterase I and 0.1 U alkaline phosphatase (Sigma) and incubated further at 37°C for 1 h.The hydrolysates were added 3 volumes of acetonitrile and centrifuged (16,000 g, 20 min, 4°C).The supernatants were lyophilized and dissolved in 50 µl water for LC-MS/MS analysis of modified and canonical ribonucleosides.Chromatographic separation was performed using an Agilent 1290 Infinity II UHPLC system with an ZORBAX RRHD Eclipse Plus C18 150 x 2.1 mm ID (1.8 mm) column protected with an ZORBAX RRHD Eclipse Plus C18 5 x 2.1 mm ID (1.8 µm) guard column (Agilent).The mobile phase consisted of water and methanol (both added 0.1% formic acid) run at 0.23 ml/min, for modifications starting with 5% methanol for 0.5 min followed by a 2.5 min gradient of 5-15% methanol, a 3 min gradient of 15-95% methanol and 4 min re-equilibration with 5% methanol.A portion of each sample was diluted for the analysis of unmodified ribonucleosides which was chromatographed isocratically with 20% methanol.Mass spectrometric detection was performed using an Agilent 6495 Triple Quadrupole system with electrospray ionization, monitoring the mass transitions 268.

Small RNA-seq by TGIRT and data analysis
Cells were treated with 5 µg/mL tunicamycin for 8 h to induce ER stress, and total RNA was subsequently extracted from the samples (refer to Section 2.5 RNA isolation).Small RNA-seq libraries were prepared and analyzed similar to previously (17).500 ng purified total RNA was used as starting material and followed with NEBNext Small RNA Library Prep Set for Illumina (NEB #E7330) with slight modifications: 20-fold diluted adaptors were used, TGIRT-III (InGex #TGIRT50) was used for reverse transcription before final PCR amplification for 15 cycles.PAGE gel purified libraries were pooled and sequenced at UAB Genomics Core Lab with 100-base single-end by Illumina NovaSeq.For data analysis, cutadapt v3.4 (21) was used to trim 3' adaptor sequence and discard trimmed reads that are shorter than 15 nucleotides or contain 5' adaptor sequence.To map small RNAs to human reference RNA sequences, unitas v1.7.3 (22) was used with setting that allows 1 internal mismatch.Reads mapped to tRFs were grouped by parental tRNA amino acid groups and tRF types.Mismatch index (0-100%) is calculated for each nucleotide on tRFs by dividing summed mismatch to summed read coverage.To plot mismatch index across different conditions, mismatch index at the indicated position is scaled across samples and represented as heatplot in Made4 R package (23).For differential analysis, microRNA and tRF counts (at least 10 read counts) were considered as input for DESeq2 wald test (24).DESeq2 results were visualized by EnhancedVolcano R package (https://github.com/kevinblighe/EnhancedVolcano).

Statistical analysis
Statistical analysis was carried out using the software GraphPad Prism 10.1.1For comparisons across different cell lines the two-tailed unequal variances Welch's t-test was applied.When comparing measurements with the same cell line, the two-tailed student's t-test was used.IC50 was determined by regression built in analysis nonlinear regression equation [inhibitor] vs. normalized response.A twoway ANOVA test with Dunnets multiple comparisons test with a single pooled variance was used for comparing each cell lines survival against the SV-HUC-1 control line for different levels of tunicamycin.

Expansion of urothelial cell lines
SV-HUC-1 cells derived from normal human urinary tract epithelium and five human bladder cancer cell lines (5637, T24, SW780, HT1376, and HT1197) were expanded in vitro in their preferred culture medium (Table 1) .The cultured cells presented quite different features.Upon seeding, SV-HUC-1 (Figure 1A), 5637 (Figure 1B), and T24 (Figure 1D) attached as individual cells.In contrast, SW780 (Figure 1C) and HT1376 (Figure 1F) formed wide clusters, while HT1197 (Figure 1E) formed more compact clusters.All cells proliferated to yield confluent cultures.Cell proliferation of all cell lines was examined by using the PrestoBlue Cell Viability Reagent (Figure 2A).The results showed the slowest proliferation for SV-HUC-1 cells compared to the rest of the bladder carcinoma cell lines, with cell lines T24 and SW780 cells proliferating the fastest.There was no correlation between the grade of the cancer cell lines and proliferation.

Cell displacement and migration
Cell migration is a key feature of numerous biological processes such as embryo development, tissue formation, immune defense, inflammation, and cancer progression.We used two different approaches to measure cell migration in the six cell lines.First, the average cell displacement speed was investigated by time-lapse imaging for the cell lines in their preferred media (Figure 2B).Five of the cell lines displayed similar average displacement speed.Interestingly, one cell line, 5637, displayed a remarkably higher average displacement speed.Second, a wound-healing assay was performed, and the results confirmed the previous observations (Figure 2C).After 12 hours, 5637 cells had migrated across the wound.

Cell survival after tunicamycin treatment
In BLCA, ER stress and the UPR have been implicated in affecting various cellular processes, including cell proliferation, apoptosis, and therapy resistance (25).Tunicamycin induces ER stress by inhibiting N-linked glycosylation, and we wanted to explore its effect on the six cell lines.Thus, in the next set of experiments, cells were treated with increasing concentrations of tunicamycin (0, 50, 100, 200, 400, 600 nM) to induce cellular stress and subsequent activation of the UPR.Subsequently, cell survival was investigated 72 hours after treatment (Figure 3).The cell lines showed varying degrees of dosage sensitivity to tunicamycin, and the specific IC50 values were determined by nonlinear regression (Supplementary Figure 3).To assess the relative sensitivity of each cell line to tunicamycin and compare it to the control cell line, we conducted a two-way ANOVA analysis (Figure 3).The results of the ANOVA, displayed in Table 2, demonstrate a significant realtionship between tunicamycin concentration, cell line and cell survival (**** p < 0.001).This indicates that the tunicamycin sensitivity varies depending on the cell line.The control cell line and four BLCA cell lines reflected a dosage-dependent tunicamycin tolerance, with SV-HUC1, 5637, and T24 having the highest sensitivity, followed by SW780 and HT1376.However, only the grade 4 cell line HT1197 survival was relatively unaffected, even with exposure to high dosages of the drug.

Expression levels of m1A regulators in urothelial carcinoma cell lines
The expression of m1A regulators has previously been investigated in various cancers, and upregulation of TRMT6, TRMT61A, and ALKBH3 has been reported in BLCA.We wanted to further investigate the expression levels of m1A methyltransferase subunits TRMT6 and TRMT61A and m1A demethylases ALKBH1 and ALKBH3 at the mRNA (Figure 4A) and protein level (Figure 4B).To our knowledge, this is the first time ALKBH1 expression has been investigated in BLCA.Noteworthy, T24 generally had low expression of all m1A regulators both at the mRNA and protein levels compared to SV-HUC-1.On the other hand, TRMT61A was highly overexpressed at the protein level in all the other BLCA cell lines.Compared to SV-HUC1, upregulated TRMT6 protein expression was observed at moderate levels in HT1376 and SW780 and, more prominently, in 5637.ALKBH1 protein expression levels were comparable to SV-HUC1 for HT1376 and 5637 and downregulated in HT1197 and SW780.ALKBH3 mRNA and protein expression were remarkably low in HT1197.In contrast, HT1376 and 5637 had higher ALKBH3 protein expression than SV-HUC1.In conclusion, there did not seem to be a correlation between expression levels and the stage or grade of the cancer cell lines.Furthermore, TRMT61A and ALKBH3 were upregulated in almost all BLCA cell lines compared to the control SV-HUC-1.

TRMT6/TRMT61A silencing inhibits cell proliferation
We further explored the role of TRMT6/TRMT61A in 5637 and HT1197 cells to investigate the contribution of the m1A methyltransferase complex.These cell lines were chosen based on the observed discrepancies in tunicamycin sensitivity and the malignant behavior of the BLCA cells.Moreover, both cell lines had high TRMT6/TRMT61A protein expression.The cells were transfected with siRNA pools targeting TRMT6 and TRMT61A to decrease the expression of these targets.The mRNA and protein expression levels of TRMT6 and TRMT61A were assessed by qRT-PCR and western blot, respectively (Figure 5A).We subsequently examined the cell proliferation of the two cell lines by the PrestoBlue Cell viability for 72 h to assess the effect of TRMT6/ TRMT61A knockdown on cell proliferation (Figure 5B).Downregulation of TRMT6/TRMT61A reduced cell proliferation in 5637, and a significant inhibitory effect was detected in HT1197.The effect of TRMT6/TRMT61A depletion on the total m1A levels in total tRNA was quantitated by liquid chromatography-tandem mass spectrometry (LC-MS/MS) in TRMT6/TRMT61A knockdown and control samples (Supplementary Figure 1).A modest decrease in m1A level was detected in the 5637 cell line after TRMT6/TRMT61A depletion.However, no noticeable difference was observed for HT1197.

TRMT6/TRMT61A knockdown inhibits BLCA cell migration in vitro
To assess the role of TRMT6/TRMT61A on cell migration, a wound-healing assay was performed after TRMT6/TRMT61A knockdown for 72 hours.For 5637, siTRMT6/TRMT61A cells showed a reduced migration capacity, although not statistically significant (Figure 6A).However, cell displacement analysis revealed that siTRMT6/TRMT61A cells had a significantly lower average displacement speed as compared to the control (Figure 6B).

TRMT6/TRMT61A knockdown influences the UPR and sensitizes cells to tunicamycin
Cells were treated with 5 µg/mL tunicamycin for 8 hours to induce ER stress and subsequent UPR activation (26).Using qRT- Characterization of cell proliferation, migration, and displacement.(A) Cell proliferation of BLCA and control cell lines grown under preferred cell culture conditions.Proliferation was compared for t = 24-, 48-and 72-hours using Welch's t-test, two tailed.Proliferation was significantly higher for T24 at t=24 (** p = 0.016) and HT1197 at t=72 (* p = 0.0317) (B) Cell displacement analysis.Average displacement speed was calculated based on particle tracking data, each well contained 4000 cells.The boxplot represents the average cell displacement speed over a time period of 12 h, n = 8-12 separate microscopic fields of view.(C) Wound healing assay to assess cell invasiveness.Quantitation of cell migration (%) at T=12 and T=24 hours after the wound was made.
PCR, mRNA expression levels of the Site-1 protease (S1P), activating transcription factor 6 (ATF6), and cyclic AMPresponsive element-binding protein 3-like protein 2 (CREB3L2) were investigated 8 hours after tunicamycin treatment and showed increased mRNA expression suggesting UPR activation (Figure 7A).S1P is involved in the cleavage and formation of the active forms of transcription factors ATF6 (ATF6 1-373 ) and CREB3L2.In 5637 cells, significantly lower mRNA levels of CREB3L2 were detected in siTRMT6/TRMT61A samples after 8h tunicamycin treatment compared to the control.In the next set of experiments, we investigated whether reduction of TRMT6/ TRMT61A affected the ability to tolerate cellular stress in 5637 and HT1197 cells.Following knockdown of TRMT6/TRMT61A, cells were treated with increasing concentrations of tunicamycin (0, 50, 100, 200, 400, 600 nM), followed by an assessment of cell survival using the PrestoBlue assay 72 h after drug treatment (Figure 7B).In 5637, depletion of TRMT6/61A decreased the tolerance for tunicamycin, reflected by a significantly lower IC50 for cells treated with siTRMT6/61A compared to the control.In HT1197, cell survival remained unchanged, and as previously observed, these cells were seemingly unaffected by even high dosages of the drug.TGIRT-based small RNA-seq was performed to investigate whether depletion of TRMT6/TRMT61A altered expression levels of small RNAs in the 5637 cell line.No significant changes in small RNA expression in siTRMT6/TRMT61A compared to siCtrl were detected (Supplementary Figures 2A, B).Cells were treated with 5 µg/mL tunicamycin for 8 hours to induce ER stress and UPR activation.However, no significant changes in small RNA expression in siTRMT6/siTRMT61A treated with tunicamycin compared to siCtrl treated with tunicamycin were detected (Supplementary Figures 2C, D).Similar to our previous report, the m1A mismatch rate was relatively decreased in siTRMT6/ TRMT61A samples across a range of tRFs (Figure 8A).On the other hand, tunicamycin treatment did not seem to have a significant effect on the tRF m1A modification status (Figures 8B-D).

Discussion
BLCA is a disease with high incidence and a contributor to mortality worldwide.Clinical management of the disease poses challenges due to a high level of heterogeneity and, thus, variable responses to treatments.Hence, there is a need for improved treatment of BLCA.The posttranscriptional RNA modification m1A is gaining increased attention as the aberrant expression of its regulatory proteins is associated with various cancers, including glioma (14), lung (27), and pancreatic (28,29) cancer.High TRMT6 has also been reported in BLCA patients (16,17), and elevated TRMT61A expression level in the 5637 BLCA cell line is suggested to play a central role in the pathogenesis of BLCA (16).Here, we characterized the expression levels of m1A methyltransferase complex TRMT6 and TRMT61A, as well as the m1A demethylases ALKBH1 and ALKBH3 in five BLCA cell lines of different BLCA grades and stages and in a control cell line derived from a healthy donor.The 5637 cell line showed elevated TRMT6 protein expression in contrast to most of the other BLCA cell lines, which had comparable protein expression levels as the control SV-HUC1.Apart from T24, all BLCA cell lines had upregulated levels of TRMT61A.However, the two subunits are both required for m1A methyltransferase activity (9).Thus, the collective m1A methyltransferase activity is upregulated in 5637 compared to the SV-HUC1 control.It also seems plausible to assume this to be the case for HT1376 and SW780, despite these cell lines having only modestly elevated TRMT6 expression compared to SV-HUC1.
m1A has been linked to various cellular processes, including cell proliferation, cell migration, and invasion, which are all central features of cancer cells.It has been shown that m1A in tRNAs can promote cell proliferation.In glioma, glioblastoma, and hepatocellular carcinoma, overexpression of TRMT6/TRMT61A promotes proliferation and malignant transformation (13)(14)(15)30).In line with this, our results showed that the downregulation of the TRMT6/TRMT61A m1A methyltransferase complex reduced cell proliferation in both 5637 and HT1197 cell lines.The downregulation of TRMT6 was recently found to suppress glioma cell migration and invasion (14).In contrast, ALKBH3 promotes cancer cell invasion through m1A demethylation, causing tRNA destabilization (12).Here, cell tracking and subsequent cell displacement analysis of 5637 cells after TRMT6/TRMT61A knockdown revealed that the average cell displacement speed was significantly reduced compared to the control.Similarly, this trend was seen in wound-healing assays for the investigation of migration properties.
The UPR is an intracellular signaling pathway utilized by eukaryotic cells to restore endoplasmic reticulum (ER) homeostasis during stress-related conditions (31).The pro-survival UPR has been reported as an important factor in the progression of several cancers, including BLCA (32).In the T24 BLCA cell line, we linked TRMT6/TRMT61A activity of mediating m1A modification status on tRF-3s to a role in the maintenance of UPR homeostasis (17).In the current study, downregulation of TRMT6/TRMT61A and subsequent treatment with the UPR stress inducer tunicamycin showed a significantly lower level of mRNA expression of CREB3L2 and slightly reduced mRNA levels of S1P, which have been previously identified as targets for m1A-modified tRF3-mediated silencing (17).Although this requires further investigation, the results could suggest that the ATF6-branch of the UPR is affected by the TRMT6/TRMT61A-m1A-axis in the 5637 cell line.The HT1197 cell line did not seem to be affected by the same mechanism.Notably, a variable degree of Tm sensitivity has been detected across gastric cancer (GC) cell lines, wherein MDR GC cells were more sensitive to Tm than the parental cells, and the increased Tm sensitivity was suggested to correlate with basal ER stress levels (33).
If the UPR cannot restore the ER to normality, ER stress leads to cell dysfunction and, ultimately, cell death (34).Here, we wanted to investigate if depletion of TRMT6/TRMT61A combined with tunicamycin-mediated induction of UPR would influence BLCA cell survival.Upon downregulation of TRMT6/TRMT61A in 5637 but not HT1197 cells had a significantly reduced cell survival rate after treatment, even with low concentrations of tunicamycin.Moreover, TGIRT-based small RNA-seq of 5637 samples revealed that depletion of TRMT6/TRMT61A resulted in decreased m1A mismatch rate across a range of tRFs compared to siCtrl which was expected.However, tunicamycin treatment did not seem to influence the m1A mismatch index of these tRFs.Based on these findings, it seems like TRMT6/TRMT61A may be involved in desensitizing BLCA cells against cellular stress.This could offer an interesting opportunity to investigate if these cells with a A B The Unfolded Protein Response and cell survival.(A) 5637 and HT1197 cells were untreated (T=0) or treated with tunicamycin (5 ug/mL) for 8h (T=8) before qRT-PCR was used to investigate mRNA expression levels of S1P, ATF6, and CREB3L2.b-actin was used as endogenous control, and mRNA expression was normalized to siCtrl T=0.Bar graphs show relative mRNA levels represented as mean ± SD, based on three independent knockdown experiments.* p value = 0.04, Unpaired T-test, two-tailed.(B) Cell survival after treatment with increasing concentration of tunicamycin.Measurements were done 72 h after drug treatment.Nonlinear regression was performed to determine the IC50 for the cell lines depleted of TRMT6/61A compared to the control.Regression line is marked as a solid line on the graph.IC50 *** p-value=0.0003.Unpaired T-test, two-tailed.
compromised stress tolerance would be more susceptible to chemotherapeutic agents commonly used in BLCA treatment to mitigate the risk of recurrence, such as mitomycin C. Notably, further studies would include a TRMT6/TRMT61A knockout model system to give more insight into the role of the TRMT6/ TRMT61A methyltransferase complex in BLCA.
BLCA is broadly categorized into NMIBC and MIBC, which comprise heterogeneous subtypes (2).In this study, the 5637 and HT1197 cell lines were used in our investigations, representing NMIBC and a metastatic BLCA, respectively.The two cell lines do not seem to respond similarly to the assays conducted after downregulating TRMT6/TRMT61A expression, which could reflect that the underlying mechanisms for BLCA pathogenesis differ in the two systems.It further raises the question of whether TRMT6/ TRMT61A is linked to any grade or stage of BLCA development.Provided that the mutation rates in these cell lines have been characterized and found to differ (35), this should not be ruled out.

6
FIGURE 6 Cell migration and displacement of 5637 are inhibited by the knockdown of TRMT6/TRMT61A.(A)Wound-healing assay.Bars represent a relative quantitation of migration distance (%) at 7 hours after the scratch was made.(B) Cell displacement analysis.Average displacement speed was calculated based on particle tracking data.The boxplot represents the average cell displacement speed over a time period of 4 hours, n = 9 separate microscopic fields of view.*** p-value = 8.56•10 -9 .

TABLE 1
Information overview of normal human urinary tract epithelium and five human bladder cancer cell lines used in this study.