Skip to main content
Download PDF

STAU2 protein level is controlled by caspases and the CHK1 pathway and regulates cell cycle progression in the non-transformed hTERT-RPE1 cells

BMC Molecular and Cell Biology volume 22, Article number: 16 (2021) Cite this article

Abstract

Background

Staufen2 (STAU2) is an RNA binding protein involved in the posttranscriptional regulation of gene expression. In neurons, STAU2 is required to maintain the balance between differentiation and proliferation of neural stem cells through asymmetric cell division. However, the importance of controlling STAU2 expression for cell cycle progression is not clear in non-neuronal dividing cells. We recently showed that STAU2 transcription is inhibited in response to DNA-damage due to E2F1 displacement from the STAU2 gene promoter. We now study the regulation of STAU2 steady-state levels in unstressed cells and its consequence for cell proliferation.

Results

CRISPR/Cas9-mediated and RNAi-dependent STAU2 depletion in the non-transformed hTERT-RPE1 cells both facilitate cell proliferation suggesting that STAU2 expression influences pathway(s) linked to cell cycle controls. Such effects are not observed in the CRISPR STAU2-KO cancer HCT116 cells nor in the STAU2-RNAi-depleted HeLa cells. Interestingly, a physiological decrease in the steady-state level of STAU2 is controlled by caspases. This effect of peptidases is counterbalanced by the activity of the CHK1 pathway suggesting that STAU2 partial degradation/stabilization fines tune cell cycle progression in unstressed cells. A large-scale proteomic analysis using STAU2/biotinylase fusion protein identifies known STAU2 interactors involved in RNA translation, localization, splicing, or decay confirming the role of STAU2 in the posttranscriptional regulation of gene expression. In addition, several proteins found in the nucleolus, including proteins of the ribosome biogenesis pathway and of the DNA damage response, are found in close proximity to STAU2. Strikingly, many of these proteins are linked to the kinase CHK1 pathway, reinforcing the link between STAU2 functions and the CHK1 pathway. Indeed, inhibition of the CHK1 pathway for 4 h dissociates STAU2 from proteins involved in translation and RNA metabolism.

Conclusions

These results indicate that STAU2 is involved in pathway(s) that control(s) cell proliferation, likely via mechanisms of posttranscriptional regulation, ribonucleoprotein complex assembly, genome integrity and/or checkpoint controls. The mechanism by which STAU2 regulates cell growth likely involves caspases and the kinase CHK1 pathway.

Background

RNA-binding proteins (RBPs) and the posttranscriptional regulation of gene expression that they impose on their associated RNAs play crucial role in cells [1,2,3]. They group mRNAs within regulons coding for functionally-related proteins and allow a timely control of protein synthesis in response to changing cellular environments [4]. RBPs are involved in diverse biological processes that include cell proliferation, DNA damage response and metabolism [5]. Given the broad spectrum of biological functions that depend on RBP activity, it is not surprising that their overexpression or depletion result in strong cellular phenotypes, indicating that their expression levels have to be tightly controlled. Moreover, mutations in genes that code for RBPs are linked to numerous diseases including neurological disorders and various types of cancer [6, 7]. However, these RBPs are part of intricate pathways in which several RBPs and RNAs associated to form functional RNA regulons involved in biological functions. Understanding the role of each molecule in these complexes and how their mutation affects specific biological systems to cause diseases is crucial to design eventual appropriate therapy.

Staufen2 (STAU2) is an RNA-binding protein [8, 9] that binds mRNAs coding for proteins involved in multiple cellular processes including cell cycle regulation [10]. The STAU2 gene, through differential splicing, generates several isoforms, the major ones having molecular masses of 52, 59 and 62 kDa [9]. STAU2 isoforms are mostly cytoplasmic, localizing near the endoplasmic reticulum [9], but can also be found in the nucleus and nucleolus [11]. STAU2 regulates mRNA expression through several posttranscriptional molecular processes such as mRNA localization, differential splicing, regulation of translation, and mRNA decay [12,13,14,15,16]. The physiological consequences of STAU2 downregulation was studied in several animal models. In zebrafish, Stau2 is required for the migration of primordial germ cells and for survival of neurons in the central nervous system [17]. In Xenopus, Stau2 controls the anterior endodermal organ formation [18]. In mouse oocytes, STAU2 is needed for meiosis progression and spindle integrity [19]. In chicken, STAU2 downregulation induces small eye development as a consequence of reduced cell proliferation [20]. Likewise, in developing mouse cortex, STAU2 regulates the balance between neural stem cell maintenance and differentiation [21, 22]: STAU2 downregulation induces cell differentiation while its overexpression produces periventricular neuronal masses. STAU2 depletion in mouse and rat brains impairs hippocampal spatial working memory, spatial novelty detection and/or associative learning and memory [23, 24]. This is consistent with the importance of STAU2 for dendritic spine morphogenesis and for long-term synaptic depression [13, 15]. STAU2 is also linked to the DNA damage response and the apoptotic pathway since STAU2 depletion causes an accumulation of DNA damage and facilitates apoptosis in the HCT116 cancer cell line [25]. In addition, induction of single-stranded break causes the inhibition of STAU2 transcription as a consequence of delocalization of the transcription factor E2F1 from the STAU2 promoter [25]. However, its biological importance and its regulation in unstressed cells are still uncharacterized.

In this paper, we now show that STAU2 knockout facilitates cell proliferation in the non-transformed hTERT-RPE1 cells. We further show that STAU2 protein steady-state level is stabilized by the CHK1 pathway and is decreased by the activity of caspases in unstressed cells. Finally, a genome wide approach reveals that STAU2 is in close proximity to proteins involved in the posttranscriptional regulation of gene expression and to proteins of the nucleolus, including proteins linked to ribosome biogenesis and to DNA repair. Altogether, our results identify novel STAU2 functions at the crossroad of cell cycle regulation and DNA damage response.

Results

STAU2 depletion facilitates cell growth of non-transformed hTERT-RPE1 cells

To study the role of STAU2 in unstressed cells, we first generated STAU2 knockout (STAU2-KO) cells using the non-transformed hTERT-RPE1 cells and the CRISPR/Cas9 technology (Supp Fig. S1). Cells were transfected with a plasmid expressing CRISPR/Cas9 and an RNA guide targeting exon 6 of the STAU2 gene. This exon is just downstream of the AUG initiation codon and common to all STAU2 isoforms. Individual clones were isolated and tested for STAU2 expression by dot blotting (Supp Fig. S1B) revealing that 26% of the selected clones were negative for STAU2 expression. DNA sequencing of STAU2 exon 6 of two STAU2-KO clones showed the presence of short deletions that introduced premature stop codons (Supp Fig. S1A). RT-qPCR quantification revealed that STAU2 mRNA levels were also decreased in all tested clones compared to that in wild type (WT) cells (Supp Fig. S1C).

To determine the impact of STAU2 knockout on cell proliferation, we studied the growth of several STAU2-KO clones compared to that of WT hTERT-RPE1 cells (Fig. 1a). Using the growth curve assay, we first showed that the STAU2-KO clone A4 proliferated faster than wild type cells, indicating that STAU2 depletion interferes with cell cycle progression and potentiates cell growth (Fig. 1a). To rule out a putative off-target effect, we used the colony growth assay to monitor proliferation of additional STAU2-KO clones. Cell growth was significantly increased in the three tested clones (Fig. 1a) compared to WT hTERT-RPE1 cells. A similar increase in cell growth was observed in shRNA-mediated STAU2 depleted cells (Fig. 1b; Supp Fig. S2). Interestingly, STAU2 depletion in HeLa tumor cells and STAU2-KO in the HCT116 cancer cell line (Supp Fig. S3) had no effect on cell proliferation (Fig. 1c,d. Supp Fig. S2) indicating that cancer cells can bypass the consequences of STAU2 depletion that are observed in non-transformed cells.

Fig. 1
figure 1

STAU2 depletion facilitates cell growth of hTERT-RPE1 cells. Growth curves (left panels) and colony growth assays (right panels) were used to monitor cell proliferation. The relative growth of wild-type cells was arbitrary fixed to 1 for the colony assays. *** p-value ≤0.001; ** p-value ≤0.01; * p-value ≤0.05. Two-tailed Student’s t-test for growth curves analysis. One sample t-test for the colony assay. (Middle) Representative Western blots showing STAU2 expression. a STAU2-KO clones (A4, B5, C6) and wild-type (WT) hTERT-RPE1 cells. n = 4. b hTERT-RPE1 cells were infected with viruses expressing shRNA control (shNT) or shRNA against STAU2 (shSTAU2). n = 5. c HeLa cells were infected with viruses expressing shRNA control (shNT) or shRNA against STAU2 (shSTAU2). n = 5. d CRISPR-mediated STAU2-KO HCT116 (STAU2-KO) and CRISPR-treated cells that still express STAU2 as controls (Ctrl). n = 5

Full size image

Caspases are involved in STAU2 degradation

These results suggest that physiological modulation of the steady-state levels of STAU2 may be advantageous for the fine tuning of cell proliferation in changing cell environment. Therefore, to identify endogenous pathway(s) that control(s) the physiological stability/degradation of STAU2, we first tested pharmacological inhibitors known to target the major protein degradation pathways. hTERT-RPE1 and HCT116 cells were incubated with the proteasome inhibitor MG132 (Fig. 2a) for 6 h or the lysosome inhibitor iLYS (Fig. 2b) for 10 h. Western blotting experiments showed that neither the proteasome inhibitor nor the lysosome inhibitor changed the steady-state levels of STAU2.

Fig. 2
figure 2

STAU2 is a target of caspase(s). Cells were incubated in the presence of different peptidase inhibitors as indicated or DMSO (vehicle). Cell extracts were analyzed by Western blotting for STAU2 expression. All Western blots are representative of at least three independently performed experiments that gave similar results. Quantification of STAU2 protein levels is indicated below the gels. The ratio of STAU2 on the loading control in DMSO-treated cells was fixed to 1. * p-value ≤0.05; ** p-value ≤0.01. One sample t-test. PARP1 was used as a measure of caspase inhibition. a, b hTERT-RPE1 and HCT116 cells were incubated in the presence of the proteasome inhibitor MG132 (10 μM) for 6 h (a) or the lysosome inhibitor iLYS (10 μM) for 10 h (b). c hTERT-RPE1 (top) and HCT116 (bottom) cells were incubated in the presence of caspase inhibitors: ZVAD-FMK (50 μM for 0, 8, 12, 24 h) or emricasan (40 μM for 24 h). d HCT116 cells were transfected with a plasmid coding for STAU259-FLAG3. Cells were then incubated in the presence of DMSO (−) or of the caspase inhibitor Emricasan (40 μM) (+) for 24 h. Full length and degradation fragments were visualized by anti-STAU2 and anti-FLAG antibodies to detect the N- and C-terminal ends of STAU2, respectively. Western blots are representative of two independently performed experiments. β-actin was used as a loading control. (bottom) Schematic representation of STAU2 protein. Predicted caspase cleavage sites are shown. Red and blue boxes: major and minor RNA-binding domains, respectively. Black boxes, regions with RNA-binding consensus sequence but lacking RNA-binding activity in vitro. Hatched boxes, tubulin-binding domain

Full size image

Western blots of HCT116 cells using anti-STAU2 antibody often revealed the presence of a specific band of around 42 kDa (Fig. 2a) that was not observed in STAU2-KO cells (Fig. 1d), suggesting a partial cleavage of STAU2 by specific peptidases in unstressed cells. Therefore, as STAU2 functions are linked to cell proliferation and apoptosis [25], we asked whether caspases, whose activities are, among others, linked to these processes [26], can modulate the steady-state levels of STAU2. Therefore, hTERT-RPE1 and HCT116 cells were treated with the pan-caspase inhibitor ZVAD-FMK for increasing periods of time (Fig. 2c). Western blot analysis indicated that the steady-state levels of STAU2 increased, indicating that caspases contribute to STAU2 degradation in unstressed cells (Fig. 2c). In hTERT-RPE1 cells, STAU2 stabilization can be observed after 8 h of treatment. In HCT116 cells, although STAU2 stabilization was observed at 24 h, the disappearance of the 42 kDa fragment began at 8 h of ZVAD treatment. Incubation of these cells with another caspase inhibitor (emricasan) also protected STAU2 from degradation (Fig. 2c). Inhibition of caspases was confirmed by the stabilization of PARP1, a known target of caspase 3. Concomitant with the stabilization of the full-length STAU2 protein, the 42 kDa band recognized with anti-STAU2 antibody in untreated HCT116 cells totally disappeared upon caspases inhibition. As the anti-STAU2 antibody recognized degradation products containing the N-terminal end of STAU2, we also expressed STAU259-FLAG3 in HCT116 cells and used the anti-FLAG antibody in an attempt to detect C-terminal fragments. In the presence of pan-caspase inhibitor, specific degradation products of around 42 and 15 kDa disappeared when using anti-STAU2 and anti-FLAG antibodies, respectively (Fig. 2d). These results identify at least one putative caspase-cleavage site in the C-terminal end of STAU2. Computer prediction identified three caspase sites that are compatible with the size of the degradation products (Fig. 2d). Altogether, these results indicate that STAU2 is a substrate of caspases in unstressed cells.

STAU2 is in close proximity to proteins involved in RNA post-transcriptional regulation

As a means to identify the pathway(s) that benefits from STAU2 stabilization/degradation, we used BioID2, a genome wide approach that identifies proteins in close proximity to STAU2 (within a range of 10 nm). hTERT-RPE1 cells were infected with a virus expressing either STAU252-biotinylase-HA or STAU252-myc (as control) fusion proteins. We first showed by fluorescence microscopy (Fig. 3a) that the two fusion proteins localized in the cytoplasm as expected from the known distribution of STAU2 in these cells. Then, we confirmed that, in the presence of biotin, STAU252-biotinylase-HA fusion protein can covalently add biotin on proteins (Fig. 3b). Then, we incubated cells with biotin for 16 h. Cells were collected, and biotinylated proteins were pulled down with streptavidin-coupled magnetic beads. Labeled proteins were identified by mass spectrometry (Supp Table S1). Interaction probability was analyzed with SAINT (Significance Analysis of INTeractome) [27] and potential off-targets were rejected using CRAPome (Contaminant Repository for Affinity Purification) [28]. We identified 325 peptides using a score of 0.7 or better (Supp Table S2). Known STAU2 interactors were found in the list of biotinylated proteins, including UPF1 [16], STAU1 [29], CDK1 [30], RSL1D1 [31], ZFR [32], and PABPC1 [33], confirming the efficacy and specificity of the assay. Proteins showing the highest enrichment (fold-changes) are presented along with a heat map representing their average spectral counts (Fig. 3c). In addition, the enriched GO terms are related to known STAU2 functions, such as translation, RNA metabolism, RNA localization, ribonucleoprotein complex biogenesis, RNA decay, splicing and stress granule formation (Supp Table S3). It is thus likely that a major consequence of STAU2 degradation will be on the posttranscriptional regulation of its targeted mRNAs. Interestingly, in addition to proteins of the posttranscriptional pathways, this approach also identified numerous proteins localized in the nucleolus including proteins involved in ribosome biogenesis, response to DNA damage, and DNA replication (Fig. 3d). STAU2-biotinylase often labeled all the proteins found within specific complexes indicating that STAU2 is in close proximity to functional complexes and not only to single proteins. For examples, not only the screen identified PRKDC (DNA-PK) but also its cofactors XRCC5 (KU80) and XRCC6 (KU70) that are recruited together at DNA damage foci [34] and all the components of the HEMIN1-DNA-PK-paraspeckel ribonucleoprotein complex [35]. Strikingly, a large percentage of the labeled proteins involved in DNA replication and/or DNA repair are linked to the CHK1 pathway (Supp Table S4), suggesting that STAU2 expression could also be modulated by this pathway. CHK1 is known to be involved in cell cycle checkpoints, DNA damage repair and DNA replication in unstressed cells [36].

Fig. 3
figure 3

Identification of proteins in proximity to STAU2. a Control (empty vector), STAU252-myc3 and STAU252-BioID2-HA infected hTERT-RPE1 cells were analyzed by immunofluorescence using anti-Myc and anti-HA antibodies to compare the subcellular localization of STAU2-tagged proteins. Both fusion proteins are cytoplasmic and mostly excluded from the nucleus, as expected. DAPI was used to stain the nucleus. b hTERT-RPE1 cells were infected with viruses expressing an empty vector or STAU252-BioID2-HA. Two days post selection, cells were incubated in the absence (−) or presence (+) of 50 μM biotin for 16 h. Biotinylated proteins were visualized by Western blotting using streptavidin-coupled HRP. c 60 proteins with the best interaction probability scores (decreasing fold-change from top to bottom) are listed along with a heat map (right) representing the average spectral counts of STAU2 interactors. d Pie chart of the percentage of STAU2 interactors associated with different biological processes. Percentages represent the number of specific interactors in each molecular process/total number of specific interactors

Full size image

STAU2 protein level is regulated by the CHK1 pathway

To determine if STAU2 is a downstream factor in the CHK1 pathway as suggested by the BioID2 assay, we treated the non-transformed hTERT-RPE1 and the cancer HCT116 cells with the CHK1 inhibitor PF-477736 (PF47) for 8 h (Fig. 4a). Our results indicated that STAU2 steady-state level significantly decreased following CHK1 inhibition. As control, we showed that the amount of CHK1 protein was downregulated in PF47-treated cells compared to untreated cells. Inhibition of CHK1 by two additional inhibitors (iCHK1 and CHIR124) also resulted in STAU2 depletion (Supp Fig. S4A), ruling out an off-target effect. Interestingly, the decrease of STAU2 protein levels coincided with a cleavage of PARP1, suggesting that CHK1 inhibition could activate caspases. To further eliminate putative off-target effect, we also treated cells with lower inhibitor concentrations (1 μM vs 20 μM) using longer incubation times. hTERT-RPE1 and HCT116 cells were incubated in the presence of CHK1 inhibitor PF47 for 48 hours, and STAU2 protein levels were quantified by Western blotting (Supp Fig. 4B). As observed above, inhibition of CHK1 with lower inhibitor concentrations still caused a decrease in STAU2 protein levels. To study the mechanism of downregulation, STAU2 protein and mRNA levels were quantified by Western blotting and RT-qPCR (Fig. 4b), respectively. Our results indicated that STAU2 steady-state level significantly decreased following CHK1 inhibition while mRNA levels were stable. Therefore, the decrease in STAU2 protein levels did not rely on transcriptional downregulation. These results indicate that STAU2 steady-state protein level is regulated by the CHK1 pathway, likely through an active mechanism of protein degradation. To determine if STAU2 is part of a feedback loop with CHK1, we quantified the amount of CHK1 protein in WT and STAU2-KO hTERT-RPE1 cells. Western blotting experiments indicated that the amount of CHK1 was decreased in STAU2-KO cell line (Fig. 4c), suggesting that STAU2 expression may have a positive feedback effect on CHK1 expression or activity.

Fig. 4
figure 4

CHK1 inhibition causes a decrease in the steady-state levels of STAU2 protein. a hTERT-RPE1 and HCT116 cells were incubated in the presence of the CHK1 inhibitor PF47 (20 μM) for 8.5 h. Cell extracts were analyzed by Western blotting. The vehicle DMSO was used as control and β-actin as a loading control. PARP1 cleavage was used as a measure of apoptosis. Quantification of STAU2 protein levels is indicated below the blots. Western blots are representative of at least three independently performed experiments that gave similar results. b Cells were incubated in the presence of CHK1 inhibitor (PF47 20 μM) for 6.5 h. STAU2 protein expression was analyzed by western blotting, while STAU2 mRNA levels were quantified by RT-qPCR. The ratio of STAU2 (protein or mRNA) on actin (protein or mRNA, respectively) in DMSO-treated cells was fixed to 1. n = 3. ** p-value ≤0.01. One sample t-test. c WT and STAU2-KO A4 hTERT-RPE1 cells were analyzed by Western blotting for expression of STAU2, CHK1 and PARP1. Actin was used as a loading control. d WT and STAU2-KO A4 hTERT-RPE1 cells were treated with the pan-caspase inhibitor emricasan. Cell proliferation using growth curve assays was quantified every day using the crystal violet retention assay. Statistic: Dunnett’s multiple comparisons (e) Dotplot representation of protein abundance in proximity to STAU2 in the presence or absence of the CHK1 inhibitor PF47. Color range indicates peptide abundance and the size of circles their relative abundance

Full size image

It is well known that CHK1 activation activates checkpoint controls, and therefore its inactivation may be involved in the enhanced cell growth observed in STAU2-depleted cells. To determine if caspase activation observed in STAU2-depleted cells may also be involved in cell growth, growth curve (Fig. 4d) and colony formation (Supp Fig. S5) assays were performed in the absence or presence of the pan-caspase inhibitor emricasan. Our results showed that caspase inhibition slightly reduced cell proliferation, suggesting that the activation of caspases observed in STAU2-KO cells might be beneficial for cell growth.

Inhibition of the CHK1 pathway dissociates STAU2 from RNA metabolism

To determine the immediate consequence of CHK1 inhibition on STAU2 proximal interactome, we identified proteins in proximity to STAU2 during the first 4 h following inhibition of CHK1 by PF47 using the Turbo-biotinylase. In contrast to BioID2 which necessitates labeling with biotin for 16–24 h, the high activity of turboID allows short incubation times and thus higher temporal resolution of protein interaction. A virus expressing STAU252/turbo-biotinylase was infected in hTERT-RPE1 cells and the activity of the fusion protein was induced by addition of biotin in the presence or absence of the CHK1 inhibitor PF47. Biotinylated proteins were isolated 4 h later and analyzed by mass spectrometry (Supp Table S5). In the absence of PF47, STAU2 is in close proximity to proteins involved in translation and RNA metabolism as described above (Supp Tables S6, S7). However, in the presence of PF47, 57 of the previously identified proximity partners were absent or less abundant (Supp Table S8). Proteins showing the highest fold change are presented in Fig. 4e. Proteins enriched in GO terms related to translation, miRNA metabolic processes and mRNA catabolism were lost (Supp Table S9). Interestingly, this change in STAU2 interactome occurs prior to STAU2 degradation, indicating that STAU2 functions are rapidly affected in response to CHK1 inhibition, before its degradation. Interestingly, only two proteins (SMEG1 and SMEG2) increased their proximity to STAU2 in these conditions. The role of these proteins is poorly understood. They were nevertheless previously described as modulators of zinc-dependent proteases via their high capacity to binding zinc [37] and shown to co-immunoprecipitate with several proteins including proteins involved in checkpoint controls and DNA damage repair [38].

Discussion

In this paper, we investigated the consequence of STAU2 depletion in non-transformed and cancer cells. We show that STAU2 depletion accelerates cell cycle progression in non-transformed cells but not in cancer cells. Interestingly, under physiological conditions, the steady state level of STAU2 protein is controlled by caspases and by the activity of the CHK1 pathway (Fig. 5). This loop likely contributes to the fine-tuning of cell proliferation in changing cell environment. The regulation of STAU2 at the protein level is in addition to that already observed at the transcriptional level by the transcription factor E2F1 [25]. Regulation of STAU2 by the E2F1 and CHK1 pathways that induce and facilitate S phase progression and DNA surveillance suggests a role for STAU2 in pathways that merge with mechanisms of DNA replication and/or DNA maintenance during S phase.

Fig. 5
figure 5

Schematic representation of the STAU2/CHK1 regulatory loop. CHK1 is well known for its role in cell cycle arrest and DNA repair following single-stranded DNA damage. In unstressed cells, CHK1 is involved in DNA replication, DNA repair and cell cycle arrest. Our results suggest that CHK1 positively controls the amount of STAU2 via inhibition of caspases. In turn, STAU2 expression positively regulates the amount of CHK1 protein, likely via the posttranscriptional regulation of CHK1 mRNA [10], linking STAU2 to cell cycle arrest and DNA repair. In addition to its role on CHK1 regulation, STAU2 may control the posttranscriptional regulation of multiple mRNAs coding for proteins involved in DNA repair, cell proliferation and DNA replication pathways. This figure was created with PowerPoint (Microsoft Office Professionnel Plus 2016)

Full size image

STAU2 regulates cell growth

The mechanism by which STAU2 depletion facilitates cell growth is not clear. STAU2, as an RNA-binding protein involved in the posttranscriptional regulation of gene expression, was previously shown to bind mRNAs coding for proteins involved in cell cycle regulation and cell cycle proliferation, including CHK1 mRNA [10]. We hypothesize that in unstressed cells, STAU2 expression may stabilize CHK1 mRNA, thus the amount of CHK1 protein, causing activation of checkpoint controls. Upon STAU2 depletion, the amount of the CHK1 protein decreased (Fig. 4c) likely via altered posttranscriptional regulation of its mRNA. The physiological decrease of CHK1 activity that releases checkpoint controls may be part of the mechanism that leads to enhanced cell proliferation upon STAU2 depletion. This mechanism contributes to the fine-tuning of both STAU2 and CHK1 expression.

Similarly, STAU2 was shown to bind 26 mRNAs involved in the TGF-β pathway, including the central protein TGF-β receptor I [10]. It was recently reported that dysregulation of the TGF-β pathway in RPE1 cells leads to cell cycle dysregulation and tumoral transformation [39]. The phenotypes observed upon STAU2 depletion are indeed compatible with oncogenic transformation. Oncogenic transformation was previously observed when RNA-binding proteins inappropriately misregulated proto-oncogenes and/or tumor suppressor genes [3, 6, 7, 40, 41]. The hyper-proliferative phenotype that results from the misregulation of these RNA-binding proteins is often followed by an accumulation of DNA damage due to the loss of cell cycle checkpoints [42, 43]. Accumulation of DNA damages are indeed observed in STAU2-KO hTERT-RPE1 cells (Supp Fig. S6) and in RNAi-mediated STAU2 knockdown in HCT116 cells [25]. Decrease in the amount of CHK1 upon STAU2 depletion may be part of the mechanism that leads to increase DNA damages. Oncogenic transformation may also explain why tumor cells are not affected by STAU2 depletion. Transformed cell lines are known to be permanently stimulated by oncogenic stimulus [44, 45], and therefore additional oncogenic stimulus may not affect their proliferation.

STAU2 protein is degraded by caspases and stabilized by the CHK1 pathway

Under physiological conditions, STAU2 protein level is decreased by caspase activity (Fig. 2) and stabilized by the CHK1 pathway (Fig. 4). As CHK1 activation inhibits caspase activity [46] and CHK1 inhibition increases PARP1 degradation via caspases (Fig. 4), it is tempting to propose that CHK1 indirectly controls STAU2 steady-state levels via the activation/inhibition of caspase activity (Fig. 5). It is known that the expression of STAU2 [25] and of CHK1 [47] is upregulated at the G1/S phase of the cell cycle by the transcription factor E2F1. Our results now support the hypothesis that the amounts of STAU2 and CHK1 proteins are then positively regulated by CHK1-mediated caspase inhibition and by STAU2-mediated posttranscriptional regulation of CHK1 mRNA, respectively. Although well characterized for their roles in response to massive DNA damages and induction of apoptosis [48,49,50,51,52], caspases also play other important roles in unstressed cells related to cell proliferation, differentiation and cellular reprogramming [53,54,55]. In turn, CHK1 is involved in cell cycle control, especially during the S phase, where it is required for firing late origin of DNA replication [56]. CHK1 regulates checkpoint controls and is essential for cell survival in dividing cells [57] and its inhibition induces an accumulation of DNA damage and apoptosis [58, 59].

Alternatively, CHK1 may control STAU2 stability via activation of downstream kinases as previously described for the CHK1-mediated control of the RNA-binding protein HuR via the regulation of CDK1 [60]. Interestingly, STAU2 is a target of CDK1 [30] and CDK1 was found in the list of proteins in proximity to STAU2 (Supp Table S1). In contrast, we have no evidence that STAU2 is a direct target of CHK1. CHK1 was not found in the list of proteins in proximity to STAU2 and, using an in-vitro kinase assay, we did not observe phosphorylation of STAU2 by a purified CHK1 kinase (data not shown).

STAU2 is in close proximity to proteins involved in RNA posttranscriptional regulation and proteins of the nucleolus

Degradation/stabilization of STAU2 and/or modulation of STAU2 functions via proximal partners should influence the pathway(s) in which STAU2 is involved. The BioID2 and TurboID experiments that detect proteins in close proximity to STAU2 link STAU2 to its well recognized functions in RNA posttranscriptional regulation. Indeed, STAU2 is well known for its roles in translation [15], RNA localization [9, 12, 15], splicing [14], mRNA decay [16], and stress granule formation [61]. It is likely that STAU2 regulation by caspases and/or the CHK1 pathway fine tunes the expression of mRNA regulons, resulting in well-ordered cell proliferation.

Interestingly, these genome-wide experiments also identify many proteins that can be found in the nucleolus. The nucleolus is well known for its role in ribosome biogenesis and for its involvement in non-ribosomal ribonucleoprotein complex formation [62]. Labeling of proteins found in the nucleolus compartment is consistent with earlier observation showing that STAU2 can migrate in the nucleolus [11]. The role of STAU2 in the nucleolus is unclear but it was proposed that STAU2 may assemble ribonucleoprotein complexes in the nucleolus [63]. The presence of ribosomes in STAU2-containing ribonucleoprotein complexes [9] further suggests a functional link between ribosome biogenesis and RNP formation in the nucleolus. STAU2 indeed co-immunoprecipitates with the ribosome biogenesis factor RSL1D1 [31]. The nucleolus might function as a checkpoint to verify the potential functional integrity of RNP and RNP-ribosome complexes or to tightly regulate their activity and/or release. Depletion of STAU2 may thus cause the formation of less-functional ribonucleoprotein complexes that could be released ahead of time or in inappropriate cellular compartments, impairing cell proliferation.

Other proteins in proximity to STAU2 are members of the DNA damage response. Most of these proteins can be found in the nucleoplasm and in the nucleolus and, therefore, could be non-specifically labeled via their co-localization with STAU2 in the nucleolus. However, STAU2 depletion causes an accumulation of DNA damages and therefore could be somehow involved in processes of DNA repair. Several processes could link STAU2 to DNA repair pathways. First, rDNA genes are the most transcribed genes and thus collisions between the transcription machinery and the replication fork are frequent [64, 65]. Consistently, several proteins involved in DNA replication are found in proximity to STAU2. The clashes between the machineries lead to replication fork stalling and single-stranded breaks [62]. This, in turn, activates the CHK1 kinase pathway [66,67,68]. In addition, R-Loops forms between nascent mRNAs and template DNA strands leading to genomic instability [69,70,71]. Remarkably, dysregulation of post-transcriptional processes increases R-loops formation [69, 72]. Depletion of STAU2 may thus facilitate the formation of R-loops whereas STAU2 stabilization upon CHK1 activation may help to attenuate the effects of R-loop formation. Alternatively, several studies have revealed the importance of lncRNAs in the DNA damage repair processes [73]. lncRNAs act as scaffolds for DNA repair protein recruitment. For example, lncRNA LINP1 participates in DNA-PK recruitment to DSB [74] and lncRNA NEAT1 is needed for the formation of DNA-PK/HEXIM1 paraspeckle complexes [35]. In other cases, lncRNAs play a role in chromatin modification, which is an indispensable step for the DNA repair process [75]. It is conceivable that STAU2 may be required to localize lncRNAs at the DNA damage sites to support the repair process.

Conlusions

In this paper, we describe a complex network of proteins in proximity to STAU2 that influences the roles of STAU2 in mRNA translation and RNA metabolism. This study also opens new research fields in post-transcriptional regulation of gene expression related to the fine-tuning of cell proliferation, ribosome biogenesis and DNA damage response.

Methods

Plasmids and cloning strategies

Plasmids coding for STAU252-FLAG3 and STAU259-FLAG3 were previously described [25]. To generate the retroviral pMSCVpuro-STAU252-BioID-HA construction, STAU252 was PCR-amplified using Phusion polymerase (NEB) and STAU252-FLAG3 as template (FW: 5′ TAAGCAGCTAGCATGCTTCAAATAAATCAGATGTTCT 3′; RV: 5′ TGCTTAACCGGTCTACCTGAAAGCCTTGAATCCT 3′). The PCR product was digested by NheI and AgeI and cloned into pcDNA3.1 BioID2-HA vector [MCS-BioID2-HA was a gift from Kyle Roux (Addgene plasmid # 74224; http://n2t.net/addgene:74224; RRID:Addgene_74224)] [76]. The resulting product was PCR-amplified and cloned into retroviral pMSCV puromycin vector after EcoRV digestion (FW: 5′.

GCTAGCATGCTTCAAATAAAT 3; RV: 5′ GTTTAAACTTAAGCTTCTATGCG 3′). Three copies of the myc sequence were introduced at the NotI site.

To generate the retroviral pMSCVpuro-STAU252-TurboID-HA and pMSCVpuro-YFP-TurboID-HA vectors, STAU252-FLAG and YFP were PCR-amplified using Phusion polymerase (NEB) and oligonucleotide primers (STAU2: FW: 5′ TAAGCAGCGGCCGCATGCTTCAAATAAATCAGATGTTCTCAG 3′; RV: 5’TGCTTAGCTAGCGGATCCGAATTCGAATCCGGAGACGTACGACCGGTCTACCTGAAAGCCTTGAATCCTTG 3′) (YFP: FW: 5′ TAAGCAGCGGCCGCGCCACCATGGTGAGCAAG 3′; RV: 5′ TGCTTAGCTAGCGGATCCGAATTCGAATCCGGAGACGTACGACCGGTCCTTGTACAGCTCGTCCATGC 3′). The PCR products were digested by NotI and NheI and cloned into pcDNA3.1 TurboID-HA vector [V5-TurboID-NES_pCDNA3 was a gift from Alice Ting (Addgene plasmid # 107169; http://n2t.net/addgene:107169; RRID: Addgene_107169)]. The resulting products were PCR-amplified and cloned into retroviral pMSCV puromycin vector after NotI digestion (STAU2: FW: 5′ TAAGCAGCGGCCGCGCCACCATGCTTCAAATAAATCAG 3′; RV: 5′ TGCTTAGCGGCCGCCTATGCGTAATCCGGTACATCGTAAGGGTATCCCTTTTCGGCAGACCGCAG 3′) (YFP: FW: 5′ TAAGCAGCGGCCGCGCCACCATGGTGAGCAAG 3′; RV: 5′ TGCTTAGCGGCCGCTTACTATGCGTAATCCGGTACATCGTAAGGGTATCCGTCCAGGGTCAGGCGCTC 3′).

Generation of STAU2-KO hTERT-RPE1 cell lines and genomic sequencing

hTERT-RPE1 cells were transfected with a plasmid coding for GFP, Cas9, and a sgRNA targeting exon 6 of the STAU2 gene (Horizon Discovery). Forty-eight hours post-transfection, GFP positive cells were sorted by FACS and plated into a 10 cm dish. Forty-eight hours after the first sorting, cells were once again isolated by FACS, and individual cells were grown into 96-well plates until colonies formed. The loss of STAU2 expression was monitored by dot blotting using anti-STAU2 antibody. Genomic DNA of several clones was isolated (Genomic DNA Miniprep Kit, Bio Basic) and PCR-amplified with specifics primers flanking exon 6 of the STAU2 gene (FW: 5′ AGCAGAATTCTTGGATAGGATAGAACAGAATTTGG 3′; RV: 5′ ATTAGGATCCACACACATAGCAGACAACATAAC 3′). PCR products were cloned into a pBluescript SK (+) vector (Stratagene) and sequenced (Sanger Sequencing Services Genome Québec).

Antibodies and reagent

Primary Antibodies against CHK1 (2G1D5), PARP1 (46D11), H2AX (D17A3), and ɣH2AX (20E3) were purchased from Cell Signaling; against GAPDH (0411), HA (12CA5) and Myc (9E10) were purchased from Santa Cruz Biotechnology. Anti-STAU2 (HPA019155), anti-β-Actin (A5441), anti-Flag (F3165), and HRP-Streptavidin were obtained from Sigma. All primary antibodies were used at 1:1000 dilution. MG132 (C2211), iCHK1 (681637) and DMSO were purchased from Millipore-Sigma; ZVAD-FMK (S7023), Emricasan (S7775), PF47 (PF-477736, S2904), and CHIR124 (S2683) were obtained from Selleckchem.

Cell culture

hTERT-RPE1 (non-transformed retinal pigment epithelial cell line immortalized with hTERT telomerase), HCT116 (colorectal carcinoma cell line) and HeLa (cervical cancer cell line) cells were cultured in Dulbecco modified Eagle’s medium (DMEM, Wisent) supplemented with 10% fetal bovine serum (Wisent), 100 μg/ml streptomycin and 100 units/ml penicillin (Wisent) under 5% CO2 atmosphere.

DNA transfection and infection

For transient expression, cells were transfected with lipofectamine 2000 (Invitrogen) or Mirus X2 (Mirus Bio. LLC) at approximatively 60% confluency. For infections, Phoenix cells were transfected at approximatively 50% confluency using lipofectamine 2000, with 10 μg of retroviral plasmids (pMSCV-puromycin) and 5 μg of packaging plasmid. Virus-containing supernatant were collected, filtered (0.45 μm), and added to target cells (HCT116, hTERT-RPE1) with polybrene (8 μg/ml). Infected cells were selected on puromycin (HCT116, 2 μg/ml; hTERT-RPE1, 5 μg/ml).

Protein electrophoresis and Western blot analysis

Total cell extracts were prepared in lysis buffer (50 mM Tris–HCl pH 7.5, 15 mM EDTA, 0.5% Triton X-100, 100 mM NaCl, 1 mM dithiothreitol [DTT] and a protease inhibitor cocktail [Roche]), resolved by SDS–polyacrylamide gel electrophoresis and transferred to nitrocellulose membrane. Membranes were blocked in 5% milk at room temperature and incubated with primary antibodies overnight at 4 °C prior to incubation with HRP-conjugated secondary antibody [polyclonal anti-mouse (1/3000), Dako: P0447; polyclonal anti-rabbit (1/5000), Dako: P0448] for 1 h at room temperature. Membranes were processed using Perkin Elmer Western Lightning Plus-ECL. Data were collected either on X-ray films (Fujifilm) or with the ChemiDoc MP Imaging System (Bio-Rad Laboratories). Western blot signals were quantified with the ImageLab (Bio-Rad Laboratories) software or with ImageJ software (NIH.gov). For the detection of biotinylated proteins by western blot, membranes were blocked in PBS-BSA 5% and incubated with HRP-streptavidin 0.3% in PBS-BSA 3% for 45 min at room temperature.

RNA isolation and RT-qPCR

Total RNA was isolated from cell extracts using the Geneaid extraction kit. Purified RNA was resuspended in 40 μl of water and digested with DNase using the TURBO DNA-free kit (Ambion). Reverse transcription reactions were done with 1 μg of RNA using RevertAid H Minus First Strand cDNA Synthesis kit (Thermo Scientific. qPCR was performed with Luna® Universal qPCR Master Mix (NEB) on a LightCycler 96 instrument (Roche). Samples were run in triplicates and normalized to actin expression.

Cell growth assays

To determine cell proliferation rates, cells were plated at the same density (Growth curve assay: HeLa: 5000 cells, hTERT-RPE1: 8000 cells and HCT116 5000 cells. Colony assay: HeLa 15,000 cells, hTERT-RPE1: 15000 or 25,000 cells, HCT116 15,000 cells) and allowed to grow for 9 days. Cell proliferation was quantified using a crystal violet retention assay (39). For growth curve assay, cells were harvested every day and the number of cells was counted with an automated hemacytometer.

BioID2 and TurboID sample preparations

One (TurboID) or two (BioID2) biological replicates were processed independently as described [77]. Briefly, stable cell line expressing STAU252-BioID-HA, STAU252-TurboID-HA, or YFP-TurboID-HA were generated by retroviral infection. After puromycin selection, cells were grown to 70% confluency and were incubated in the presence of 50 μM biotin for 16 h (BioID2) or 4 h (TurboID) and lysed in 700 μL RIPA lysis buffer (50 mM Tris pH 8, 150 mM NaCl, 0.1% SDS, 0.5% sodium deoxycholate, 1% Triton X-100, 1× protease inhibitor cocktail, Sigma-Aldrich). Cell extracts were sonicated and biotinylated proteins pulled down with 250 μL streptavidin-coated magnetic beads (Dynabeads MyOne Streptavidin T1). Beads were then washed twice with 1 mL of RIPA lysis buffer, once with 1 mL of 1 M KCl, once with 1 mL of 0.1 M Na2CO3, once with 1 mL of 2 M urea in 10 mM Tris-HCl (pH 8.0), and three time with 1 mL RIPA lysis buffer. Beads were then resuspended in 100 μL 50 mM ammonium bicarbonate solution.

Protein digestion and LC-MS/MS

Proteins were identified at the proteomic core facility for LC-MS/MS analysis. Proteins were digested on beads with 1 μg Sequencing Grade Modified Trypsin (Promega) overnight at 37 °C. Peptides in the supernatant were reduced with 9 mM dithiothreitol at 37 °C for 30 min and alkylated with 17 mM iodoacetamide at room temperature for 20 min. After desalting, peptides were eluted in 10% ammonium hydroxide/90% methanol (v/v) and resuspended in 5% FA. Peptides were loaded into a 75 μm i.d. × 150 mm Self-Pack C18 column installed in the Easy-nLC II system (Proxeon Biosystems). The buffers used for chromatography were 0.2% formic acid (buffer A) and 90% acetonitrile/0.2% formic acid (buffer B). Peptides were eluted with a two slope gradient at a flowrate of 250 nL/min. Solvent B first increased from 1 to 35% in 105 min and then from 35 to 84% B in 15 min. The HPLC system was coupled to a Q Exactive mass spectrometer (Thermo Scientific) through a Nanospray Flex Ion Source. Nanospray and S-lens voltages were set to 1.3–1.8 kV and 50 V, respectively. Capillary temperature was set to 225 °C. Full scan MS survey spectra (m/z 360–2000) in profile mode were acquired in the Orbitrap with a resolution of 70,000. The 16 most intense peptide ions were fragmented in the collision cell and MS/MS spectra were analyzed in the Orbitrap.

Protein identification and enrichment analysis

The peak list files were generated with Proteome Discoverer (version 2.1) using the following parameters: minimum mass set to 500 Da, maximum mass set to 6000 Da, no grouping of MS/MS spectra, precursor charge set to auto, and minimum number of fragment ions set to 5. Protein database searching was performed with Mascot 2.6 (Matrix Science) against the Uniprot human protein database (May 16th, 2018). The mass tolerances for precursor and fragment ions were set to 10 ppm and 0.6 Da, respectively. Trypsin was used as the enzyme allowing for up to 1 missed cleavage. Cysteine carbamidomethylation was specified as a fixed modification, and methionine oxidation as variable modifications. Data interpretation was performed using Scaffold (version 4.8.9). Protein lists were submitted to the CRAPome software to calculate SAINT score and the fold change scores (FC-A, FC-B). P Proteins with SAINT probability (SP) greater than 0.7 were selected. REPRINT tools were used to generate visual analysis of our mass spectrometry data.

Epifluorescence microscopy

Cells were seeded onto 12 × 12 mm glass coverslips in 6-well plates and incubated overnight at 37 °C. Cells were fixed in 4% PFA for 15 min at room temperature and permeabilized in 10 mM Tris (pH 7.5), 0.5% Triton-X-100 (v/v) and 0.1% BSA (w/v) for 10 min. Cells were then stained with anti-myc (1600) or anti-HA: (1500) for 1 h at room temperature. Secondary antibody (Alexa Fluor 488 goat) was added for 1 h at room temperature. Coverslips were washed once with PBS and incubated for 10 min with PBS containing DAPI (0.5 μg/ml) at room temperature. Coverslips were mounted using Dako mounting medium (Dako Faramount Aqueous Mounting Medium; S3025). Images were acquired with an ECLIPSE TE2000U microscope using 100X/1.40 Oil objective and chroma filter 41,001 (GFP) or chroma filter 31,000 (DAPI) (NIKON). Images processing was performed using ImageJ software. The data are representative of at least 15 fields of view.

Availability of data and materials

All data generated or analyzed during this study are included in this published article [and its supplementary information files].

Abbreviations

STAU2:

Staufen2

CRISPR:

Clustered Regularly Interspaced Short Palindromic Repeats

CHK1:

checkpoint kinase 1

DDR:

DNA damage response

RBP:

RNA-binding protein

KO:

knockout

WT:

wild type

shRNA:

short hairpin RNA

GO:

gene ontology

lncRNA:

long non-coding RNA

CRAPome:

Contaminant Repository for Affinity Purification

REPRINT:

Resource for Evaluation of Protein Interaction Networks

SAINT:

Significance Analysis of INTeractome

References

  1. Hong S. RNA binding protein as an emerging therapeutic target for Cancer prevention and treatment. J Cancer Prev. 2017;22(4):203–10.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Moore S, Jarvelin AI, Davis I, Bond GL, Castello A. Expanding horizons: new roles for non-canonical RNA-binding proteins in cancer. Curr Opin Genet Dev. 2017;48:112–20.

    Article  PubMed  CAS  Google Scholar 

  3. Pereira B, Billaud M, Almeida R. RNA-binding proteins in Cancer: old players and new actors. Trends Cancer. 2017;3(7):506–28.

    Article  CAS  PubMed  Google Scholar 

  4. Keene JD. RNA regulons: coordination of post-transcriptional events. Nature Reviews. 2007;8(7):533–43.

    Article  CAS  PubMed  Google Scholar 

  5. Huang R, Han M, Meng L, Chen X. Transcriptome-wide discovery of coding and noncoding RNA-binding proteins. Proc Natl Acad Sci U S A. 2018;115(17):E3879–E87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Lukong KE, Chang KW, Khandjian EW, Richard S. RNA-binding proteins in human genetic disease. Trends Genet. 2008;24(8):416–25.

    Article  CAS  PubMed  Google Scholar 

  7. Castello A, Fischer B, Hentze MW, Preiss T. RNA-binding proteins in Mendelian disease. Trends Genet. 2013;29(5):318–27.

    Article  CAS  PubMed  Google Scholar 

  8. Buchner G, Bassi MT, Andolfi G, Ballabio A, Franco B. Identification of a novel homolog of the Drosophila staufen protein in the chromosome 8q13-q21.1 region. Genomics. 1999;62(1):113–8.

    Article  CAS  PubMed  Google Scholar 

  9. Duchaine TF, Hemraj I, Furic L, Deitinghoff A, Kiebler MA, DesGroseillers L. Staufen2 isoforms localize to the somatodendritic domain of neurons and interact with different organelles. J Cell Sci. 2002;115(Pt 16):3285–95.

    CAS  PubMed  Google Scholar 

  10. Furic L, Maher-Laporte M, DesGroseillers L. A genome-wide approach identifies distinct but overlapping subsets of cellular mRNAs associated with Staufen1- and Staufen2-containing ribonucleoprotein complexes. RNA. 2008;14(2):324–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Macchi P, Brownawell AM, Grunewald B, DesGroseillers L, Macara IG, Kiebler MA. The brain-specific double-stranded RNA-binding protein Staufen2: nucleolar accumulation and isoform-specific exportin-5-dependent export. J Biol Chem. 2004;279(30):31440–4.

    Article  CAS  PubMed  Google Scholar 

  12. Tang SJ, Meulemans D, Vazquez L, Colaco N, Schuman E. A role for a rat homolog of staufen in the transport of RNA to neuronal dendrites. Neuron. 2001;32(3):463–75.

    Article  CAS  PubMed  Google Scholar 

  13. Goetze B, Tuebing F, Xie Y, Dorostkar MM, Thomas S, Pehl U, et al. The brain-specific double-stranded RNA-binding protein Staufen2 is required for dendritic spine morphogenesis. J Cell Biol. 2006;172(2):221–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. O'Leary DA, Sharif O, Anderson P, Tu B, Welch G, Zhou Y, et al. Identification of small molecule and genetic modulators of AON-induced dystrophin exon skipping by high-throughput screening. PLoS One. 2009;4(12):e8348.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Lebeau G, Miller LC, Tartas M, McAdam R, Laplante I, Badeaux F, et al. Staufen 2 regulates mGluR long-term depression and map 1b mRNA distribution in hippocampal neurons. Learning Memory. 2011;18(5):314–26.

    Article  CAS  PubMed  Google Scholar 

  16. Park E, Gleghorn ML, Maquat LE. Staufen2 functions in Staufen1-mediated mRNA decay by binding to itself and its paralog and promoting UPF1 helicase but not ATPase activity. Proc Natl Acad Sci U S A. 2013;110(2):405–12.

    Article  CAS  PubMed  Google Scholar 

  17. Ramasamy S, Wang H, Quach HN, Sampath K. Zebrafish Staufen1 and Staufen2 are required for the survival and migration of primordial germ cells. Dev Biol. 2006;292(2):393–406.

    Article  CAS  PubMed  Google Scholar 

  18. Bilogan CK, Horb ME. Xenopus staufen2 is required for anterior endodermal organ formation. Genesis. 2012;50(3):251–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Cao Y, Du J, Chen D, Wang Q, Zhang N, Liu X, et al. RNA- binding protein Stau2 is important for spindle integrity and meiosis progression in mouse oocytes. Cell Cycle. 2016;15(19):2608–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Cockburn DM, Charish J, Tassew NG, Eubanks J, Bremner R, Macchi P, et al. The double-stranded RNA-binding protein Staufen 2 regulates eye size. Mol Cell Neurosci. 2012;51(3–4):101–11.

    Article  CAS  PubMed  Google Scholar 

  21. Kusek G, Campbell M, Doyle F, Tenenbaum SA, Kiebler M, Temple S. Asymmetric segregation of the double-stranded RNA binding protein Staufen2 during mammalian neural stem cell divisions promotes lineage progression. Cell Stem Cell. 2012;11(4):505–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Vessey JP, Amadei G, Burns SE, Kiebler MA, Kaplan DR, Miller FD. An asymmetrically localized Staufen2-dependent RNA complex regulates maintenance of mammalian neural stem cells. Cell Stem Cell. 2012;11(4):517–28.

    Article  CAS  PubMed  Google Scholar 

  23. Berger SM, Fernandez-Lamo I, Schonig K, Moya SMF, Ehses J, Schieweck R, et al. Forebrain-specific, conditional silencing of Staufen2 alters synaptic plasticity, learning, and memory in rats. Genome Biol. 2017;18(1):222.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Popper B, Demleitner A, Bolivar VJ, Kusek G, Snyder-Keller A, Schieweck R, et al. Staufen2 deficiency leads to impaired response to novelty in mice. Neurobiol Learn Mem. 2018;150:107–15.

    Article  CAS  PubMed  Google Scholar 

  25. Zhang X, Trepanier V, Beaujois R, Viranaicken W, Drobetsky E, DesGroseillers L. The downregulation of the RNA-binding protein Staufen2 in response to DNA damage promotes apoptosis. Nucleic Acids Res. 2016;44(8):3695–712.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. McArthur K, Kile BT. Apoptotic Caspases: multiple or mistaken identities? Trends Cell Biol. 2018;28(6):475–93.

    Article  CAS  PubMed  Google Scholar 

  27. Choi H, Larsen B, Lin ZY, Breitkreutz A, Mellacheruvu D, Fermin D, et al. SAINT: probabilistic scoring of affinity purification-mass spectrometry data. Nat Methods. 2011;8(1):70–3.

    Article  CAS  PubMed  Google Scholar 

  28. Mellacheruvu D, Wright Z, Couzens AL, Lambert JP, St-Denis NA, Li T, et al. The CRAPome: a contaminant repository for affinity purification-mass spectrometry data. Nature Methods. 2013;10(8):730.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Elbarbary RA, Li W, Tian B, Maquat LE. STAU1 binding 3′ UTR IRAlus complements nuclear retention to protect cells from PKR-mediated translational shutdown. Genes Dev. 2013;27(13):1495–510.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Beaujois R, Ottoni E, Zhang X, Gagnon C, HSine S, Mollet S, et al. The M-phase specific hyperphosphorylation of Staufen2 involved the cyclin-dependent kinase CDK1. BMC Cell Biol. 2017;18(1):25.

  31. Lessard F, Igelmann S, Trahan C, Huot G, Saint-Germain E, Mignacca L, et al. Senescence-associated ribosome biogenesis defects contributes to cell cycle arrest through the Rb pathway. Nat Cell Biol. 2018;20(7):789–99.

    Article  CAS  PubMed  Google Scholar 

  32. Elvira G, Massie B, DesGroseillers L. The zinc-finger protein ZFR is critical for Staufen 2 isoform specific nucleocytoplasmic shuttling in neurons. J Neurochem. 2006;96(1):105–17.

    Article  CAS  PubMed  Google Scholar 

  33. Maher-Laporte M, Berthiaume F, Moreau M, Julien LA, Lapointe G, Mourez M, et al. Molecular composition of staufen2-containing ribonucleoproteins in embryonic rat brain. PLoS One. 2010;5(6):e11350.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Gottlieb TM, Jackson SP. The DNA-dependent protein kinase: requirement for DNA ends and association with Ku antigen. Cell. 1993;72(1):131–42.

    Article  CAS  PubMed  Google Scholar 

  35. Morchikh M, Cribier A, Raffel R, Amraoui S, Cau J, Severac D, et al. HEXIM1 and NEAT1 long non-coding RNA form a multi-subunit complex that regulates DNA-mediated innate immune response. Mol Cell. 2017;67(3):387–99 e5.

    Article  CAS  PubMed  Google Scholar 

  36. Zhang Y, Hunter T. Roles of Chk1 in cell biology and cancer therapy. Int J Cancer. 2014;134(5):1013–23.

    Article  CAS  PubMed  Google Scholar 

  37. Jonsson M, Lundwall A, Malm J. The semenogelins: proteins with functions beyond reproduction? Cell Mol Life Sci. 2006;63(24):2886–8.

    Article  CAS  PubMed  Google Scholar 

  38. Huttlin EL, Ting L, Bruckner RJ, Gebreab F, Gygi MP, Szpyt J, et al. The BioPlex network: a systematic exploration of the human Interactome. Cell. 2015;162(2):425–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Parapuram SK, Chang B, Li L, Hartung RA, Chalam KV, Nair-Menon JU, et al. Differential effects of TGFbeta and vitreous on the transformation of retinal pigment epithelial cells. Invest Ophthalmol Vis Sci. 2009;50(12):5965–74.

    Article  PubMed  Google Scholar 

  40. Sommer G, Dittmann J, Kuehnert J, Reumann K, Schwartz PE, Will H, et al. The RNA-binding protein La contributes to cell proliferation and CCND1 expression. Oncogene. 2011;30(4):434–44.

    Article  CAS  PubMed  Google Scholar 

  41. Cookson MR. RNA-binding proteins implicated in neurodegenerative diseases. Wiley Interdiscip Rev RNA. 2017;8(1). https://doi.org/10.1002/wrna.1397.

  42. Shortt J, Johnstone RW. Oncogenes in cell survival and cell death. Cold Spring Harb Perspect Biol. 2012;4(12).

  43. Hills SA, Diffley JF. DNA replication and oncogene-induced replicative stress. Curr Biol. 2014;24(10):R435–44.

    Article  CAS  PubMed  Google Scholar 

  44. Ahmed D, Eide PW, Eilertsen IA, Danielsen SA, Eknaes M, Hektoen M, et al. Epigenetic and genetic features of 24 colon cancer cell lines. Oncogenesis. 2013;2:e71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Mittelman D, Wilson JH. The fractured genome of HeLa cells. Genome Biol. 2013;14(4):111.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Myers K, Gagou ME, Zuazua-Villar P, Rodriguez R, Meuth M. ATR and Chk1 suppress a caspase-3-dependent apoptotic response following DNA replication stress. PLoS Genet. 2009;5(1):e1000324.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Carrassa L, Broggini M, Vikhanskaya F, Damia G. Characterization of the 5'flanking region of the human Chk1 gene: identification of E2F1 functional sites. Cell Cycle. 2003;2(6):604–9.

    Article  CAS  PubMed  Google Scholar 

  48. Zhao H, Piwnica-Worms H. ATR-mediated checkpoint pathways regulate phosphorylation and activation of human Chk1. Mol Cell Biol. 2001;21(13):4129–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Jackson SP, Bartek J. The DNA-damage response in human biology and disease. Nature. 2009;461(7267):1071–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Blasius M, Forment JV, Thakkar N, Wagner SA, Choudhary C, Jackson SP. A phospho-proteomic screen identifies substrates of the checkpoint kinase Chk1. Genome Biol. 2011;12(8):R78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. McIlwain DR, Berger T, Mak TW. Caspase functions in cell death and disease. Cold Spring Harb Perspect Biol. 2013;5(4):a008656.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Patil M, Pabla N, Dong Z. Checkpoint kinase 1 in DNA damage response and cell cycle regulation. Cell Mol Life Sci. 2013;70(21):4009–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Li J, Yuan J. Caspases in apoptosis and beyond. Oncogene. 2008;27(48):6194–206.

    Article  CAS  PubMed  Google Scholar 

  54. Kuranaga E. Beyond apoptosis: caspase regulatory mechanisms and functions in vivo. Genes Cells. 2012;17(2):83–97.

    Article  CAS  PubMed  Google Scholar 

  55. Baena-Lopez LA, Arthurton L, Xu DC, Galasso A. Non-apoptotic Caspase regulation of stem cell properties. Semin Cell Dev Biol. 2018;82:118–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Petermann E, Woodcock M, Helleday T. Chk1 promotes replication fork progression by controlling replication initiation. Proc Natl Acad Sci U S A. 2010;107(37):16090–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Wilsker D, Petermann E, Helleday T, Bunz F. Essential function of Chk1 can be uncoupled from DNA damage checkpoint and replication control. Proc Natl Acad Sci U S A. 2008;105(52):20752–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Rodriguez R, Meuth M. Chk1 and p21 cooperate to prevent apoptosis during DNA replication fork stress. Mol Biol Cell. 2006;17(1):402–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Meuth M. Chk1 suppressed cell death. Cell Div. 2010;5:21.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  60. Kim HH, Abdelmohsen K, Gorospe M. Regulation of HuR by DNA damage response kinases. J Nucleic Acids. 2010;2010.

  61. Thomas MG, Martinez Tosar LJ, Loschi M, Pasquini JM, Correale J, Kindler S, et al. Staufen recruitment into stress granules does not affect early mRNA transport in oligodendrocytes. Mol Biol Cell. 2005;16(1):405–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Lindstrom MS, Jurada D, Bursac S, Orsolic I, Bartek J, Volarevic S. Nucleolus as an emerging hub in maintenance of genome stability and cancer pathogenesis. Oncogene. 2018;37(18):2351–66.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Kiebler MA, Jansen RP, Dahm R, Macchi P. A putative nuclear function for mammalian Staufen. Trends Biochem Sci. 2005;30(5):228–31.

    Article  CAS  PubMed  Google Scholar 

  64. Takeuchi Y, Horiuchi T, Kobayashi T. Transcription-dependent recombination and the role of fork collision in yeast rDNA. Genes Dev. 2003;17(12):1497–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Garcia-Muse T, Aguilera A. Transcription-replication conflicts: how they occur and how they are resolved. Nat Rev Mol Cell Biol. 2016;17(9):553–63.

    Article  CAS  PubMed  Google Scholar 

  66. Zou L, Elledge SJ. Sensing DNA Damage through ATRIP recognition of RPA-ssDNA complexes. Science. 2003;300(5625):1542–8.

    Article  CAS  PubMed  Google Scholar 

  67. Shechter D, Costanzo V, Gautier J. ATR and ATM regulate the timing of DNA replication origin firing. Nat Cell Biol. 2004;6(7):648–55.

    Article  CAS  PubMed  Google Scholar 

  68. Jazayeri A, Falck J, Lukas C, Bartek J, Smith GC, Lukas J, et al. ATM- and cell cycle-dependent regulation of ATR in response to DNA double-strand breaks. Nat Cell Biol. 2006;8(1):37–45.

    Article  CAS  PubMed  Google Scholar 

  69. Li X, Manley JL. Cotranscriptional processes and their influence on genome stability. Genes Dev. 2006;20(14):1838–47.

    Article  CAS  PubMed  Google Scholar 

  70. Aguilera A, Gomez-Gonzalez B. Genome instability: a mechanistic view of its causes and consequences. Nature Reviews. 2008;9(3):204–17.

    Article  CAS  PubMed  Google Scholar 

  71. Hamperl S, Cimprich KA. The contribution of co-transcriptional RNA:DNA hybrid structures to DNA damage and genome instability. DNA Repair (Amst). 2014;19:84–94.

    Article  CAS  Google Scholar 

  72. Paulsen RD, Soni DV, Wollman R, Hahn AT, Yee MC, Guan A, et al. A genome-wide siRNA screen reveals diverse cellular processes and pathways that mediate genome stability. Mol Cell. 2009;35(2):228–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Su M, Wang H, Wang W, Wang Y, Ouyang L, Pan C, et al. LncRNAs in DNA damage response and repair in cancer cells. Acta Biochim Biophys Sin Shanghai. 2018;50(5):433–9.

    Article  CAS  PubMed  Google Scholar 

  74. Lees-Miller SP, Beattie TL, Tainer JA. Noncoding RNA joins Ku and DNA-PKcs for DNA-break resistance in breast cancer. Nat Struct Mol Biol. 2016;23(6):509–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Dianatpour A, Ghafouri-Fard S. The role of long non coding RNAs in the repair of DNA double Strand breaks. Int J Mol Cell Med. 2017;6(1):1–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  76. Kim DI, Jensen SC, Noble KA, Kc B, Roux KH, Motamedchaboki K, et al. An improved smaller biotin ligase for BioID proximity labeling. Mol Biol Cell. 2016;27(8):1188–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Branon TC, Bosch JA, Sanchez AD, Udeshi ND, Svinkina T, Carr SA, et al. Efficient proximity labeling in living cells and organisms with TurboID. Nature Biotechnology. 2018;36(9):880.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Tripathi S, Pohl MO, Zhou YY, Rodriguez-Frandsen A, Wang GJ, Stein DA, et al. Meta- and orthogonal integration of influenza "OMICs” data defines a role for UBR4 in virus budding. Cell Host Microbe 2015;18(6):723–735.

Download references

Acknowledgements

We thank Dr. Xin Zhang for useful suggestions, Louise Cournoyer for help with cell culture, Alicia Pellerin-Viger, Bellastrid Doran, for technical assistance, and Dr. Armelle Le Campion for her support with flow cytometry. We also thank the proteomic core facilities at the Institut de Recherches en Immunologie et Cancer (IRIC) and at the Institut de Recherches Cliniques de Montréal (IRCM) for the proteomics analysis.

Funding

This work was supported by a grant (RGPIN-2019-05027) from the Natural Sciences and Engineering Research Council of Canada (NSERC) to LDG. The funding agency had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

  1. Département de Biochimie et Médecine Moléculaire, Faculté de médecine, Université de Montréal, 2900 Édouard Montpetit, Montréal, QC, H3T 1J4, Canada

    Lionel Condé, Yulemi Gonzalez Quesada, Florence Bonnet-Magnaval, Rémy Beaujois & Luc DesGroseillers

Authors
  1. Lionel Condé
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yulemi Gonzalez Quesada
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Florence Bonnet-Magnaval
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Rémy Beaujois
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Luc DesGroseillers
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

LC, YGQ, FBM, and RB acquired data. LC, YGQ, FBM, RB and LDG designed the work, analyzed and interpreted data. LC and LDG wrote the manuscript. LC, YGQ, FBM, RB and LDG approved the final version, including author’s contributions to the study.

Corresponding author

Correspondence to Luc DesGroseillers.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Knockout of STAU2 in the nontransformed hTERT-RPE1 cell line by the CRISPR/Cas9 technique. (A) Schematic representation of the STAU2 gene. Exons are indicated as well as the differential splicing of exon 5 and the differential choice of stop codons that contribute to the expression of several protein isoforms. The 3’end sequence of STAU2 exon 6 of the human genome (WT) is shown. The target sequence of the RNA guide used for CRISPR gene editing is underlined. The sequences of two different CRISPR/Cas9-derived STAU2-KO clones are shown below. Dashed lines indicate deleted nucleotides; the premature stop codons are indicated. (B) hTERT-RPE1 cells were transfected with plasmids expressing a Cas9/sgRNA complex targeting exon 6 of the STAU2 gene. Colonies grown from single cells were screened for STAU2 protein expression by dot blotting. (C) Four different STAU2-KO clones and one CRISPR-derived clones that express STAU2 (I2) were analyzed by RT-qPCR for STAU2 expression. The ratio of STAU2 mRNAs on actin mRNA in CRISPR-control cells (I2) cells was arbitrary fixed to 1. The graph represents the means and standard deviation of three independently performed experiments. *** p-value ≤0.001; ** p-value ≤0.01. One sample t-test. Note that clones A4, B5 and C6 were originally named F3, F4 and G4, respectively.

Additional file 2: Figure S2.

STAU2 depletion facilitates cell growth. Representative images of the colony growth assays (quantified in Fig. 1). The pictures are representative of three independently performed experiments. (A) WT and STAU2-KO (A4, B5) hTERT-RPE1 cells. (B,C) hTERT-RPE1 (B) and HeLa (C) cells infected with non-targeting shRNA (shNT) or shRNA against STAU2 (shSTAU2). (D) CRISPR-infected cells that still express STAU2 (Ctrl) and STAU2-KO (STAU2-KO) HCT116 cells.

Additional file 3: Figure S3.

Knockout of STAU2 in HCT116 cancer cells. (A) Schematic representation of the STAU2 gene, including the 5’end sequence of STAU2 exon 6 of the human genome and the position of the RNA guide RNA (underlined). (B) Western blot of CRISPR-transfected HCT116 cells grown from single cells to monitor STAU2 protein expression. 35% of the selected clones were negative for STAU2 expression.

Additional file 4: Figure S4.

CHK1 inhibition causes a decrease in the steady-state levels of STAU2 protein. (A) HCT116 cells were incubated in the presence of CHK1 inhibitors (PF47 20 μM, iCHK1 20 μM for 8.5 h and CHIR124 200 nM for 24 h). (B) hTERT-RPE1 and HCT116 cells were incubated in the presence of low concentration of the CHK1 inhibitor PF47 (1 μM) for 48 h. Cell extracts were analyzed by Western blotting. The vehicle DMSO was used as control and β-actin as a loading control. PARP1 cleavage was used as a measure of apoptosis. Quantification of STAU2 protein levels is indicated below the blots. Western blots are representative of at least three independently performed experiments that gave similar results.

Additional file 5: Figure S5.

Caspases inhibition alters cell growth. WT and STAU2-KO A4 hTERT-RPE1 cells were treated with the pan-caspase inhibitor emricasan and allow to grow for 7 days. Colony growth assays were used to monitor cell proliferation. Left: representative growth of cells plated in triplicates. Right: Quantification of cell growth from three independently performed experiments. The relative growth of wild-type cells was arbitrary fixed to 1. ** p-value ≤0.01; * p-value ≤0.05. One sample t-test.

Additional file 6: Figure S6.

STAU2 depleted cells accumulate DNA damages. (A) Protein extracts isolated from WT and STAU2-KO hTERT-RPE1 cells were analyzed by Western Blotting for STAU2 and H2AX expression. Phosphorylated H2AX (γH2AX), a marker of DNA damages, was also revealed. The western blot is representative of three independently performed experiments. (B) Quantification of ɣH2AX protein levels in four different hTERT-RPE1 STAU2-KO cells compared to WT cells. Protein quantification represents the means and standard deviation of three independently performed experiments. The ratio of ɣH2AX on H2AX in wild type (WT) cells was arbitrary fixed to 1. ** p-value ≤0.01. One sample t-test.

Additional file 7: Figure S7.

Uncropped images of immunoblots used in figures.

Additional file 8: Table S1.

List of all biotinylated proteins in the BioID2 assay.

Additional file 9: Table S2.

List of STAU2 interactors with an interaction probability score of 0.7 or better as calculated with the SAINT and CRAPome softwares.

Additional file 10: Table S3.

Gene ontology (GO) of STAU2 interactors in hTERT-RPE1 cells. List of significant shared GO terms that are over-represented in BioID2, using Metascape (A Gene Annotation & Analysis Resources [78]).

Additional file 11: Table S4.

List of STAU2 interactors involved in DNA repair and/or replication pathways. Proteins linked to the CHK1 pathway are indicated.

Additional file 12: Table S5.

List of all biotinylated peptides in the TurboID assay. The TurboID assay was performed in standard conditions or in the presence of the CHK1 inhibitor PF47.

Additional file 13: Table S6.

List of STAU2 interactors in the TurboID assay with an interaction probability score of 0.7 or better, using the SAINT and CRAPome softwares.

Additional file 14: Table S7.

Gene ontology (GO) of STAU2 interactors in hTERT-RPE1 cells. List of significant shared GO terms that are over-represented in TurboID, using Metascape (A Gene Annotation & Analysis Resources [78]).

Additional file 15: Table S8.

List of STAU2 interactors that are lost in PF47 conditions.

Additional file 16: Table S9.

Gene ontology (GO) of STAU2 interactors that show reduced proximity to STAU2 in hTERT-RPE1 cells treated with the CHK1 inhibitor PF47 compared to untreated cells. List of significant shared GO terms that are over-represented in TurboID, using Metascape (A Gene Annotation & Analysis Resources [78]).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Condé, L., Gonzalez Quesada, Y., Bonnet-Magnaval, F. et al. STAU2 protein level is controlled by caspases and the CHK1 pathway and regulates cell cycle progression in the non-transformed hTERT-RPE1 cells. BMC Mol and Cell Biol 22, 16 (2021). https://doi.org/10.1186/s12860-021-00352-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12860-021-00352-y

Keywords

BMC Molecular and Cell Biology

ISSN: 2661-8850

Contact us