- Research
- Open access
- Published:
Multi-omics profiling reveals dysregulated ribosome biogenesis and impaired cell proliferation following knockout of CDR2L
BMC Cancer volume 24, Article number: 645 (2024)
Abstract
Background
Cerebellar degeneration-related (CDR) proteins are associated with paraneoplastic cerebellar degeneration (PCD) – a rare, neurodegenerative disease caused by tumour-induced autoimmunity against neural antigens resulting in degeneration of Purkinje neurons in the cerebellum. The pathogenesis of PCD is unknown, in large part due to our limited understanding of the functions of CDR proteins. To this end, we performed an extensive, multi-omics analysis of CDR-knockout cells focusing on the CDR2L protein, to gain a deeper understanding of the properties of the CDR proteins in ovarian cancer.
Methods
Ovarian cancer cell lines lacking either CDR1, CDR2, or CDR2L were analysed using RNA sequencing and mass spectrometry-based proteomics to assess changes to the transcriptome, proteome and secretome in the absence of these proteins.
Results
For each knockout cell line, we identified sets of differentially expressed genes and proteins. CDR2L-knockout cells displayed a distinct expression profile compared to CDR1- and CDR2-knockout cells. Knockout of CDR2L caused dysregulation of genes involved in ribosome biogenesis, protein translation, and cell cycle progression, ultimately causing impaired cell proliferation in vitro. Several of these genes showed a concurrent upregulation at the transcript level and downregulation at the protein level.
Conclusions
Our study provides the first integrative multi-omics analysis of the impact of knockout of the CDR genes, providing both new insights into the biological properties of the CDR proteins in ovarian cancer, and a valuable resource for future investigations into the CDR proteins.
Introduction
Paraneoplastic neurological syndromes are rare, immune-mediated disorders triggered by cancer [1]. These syndromes are characterized by circulating antibodies directed against antigens expressed by neurons and cancer cells. Paraneoplastic cerebellar degeneration (PCD) is one of the most common forms of paraneoplastic neurological syndromes [2, 3]. It presents with subacute onset of limb and truncal ataxia, dysarthria, and nystagmus [4]. In patients with PCD, the dominant autoantibody detected in serum and cerebrospinal fluid is anti-Yo which is most frequently seen in patients with breast, ovarian or other gynaecological cancers [2, 4]. Anti-Yo targets three intracellular antigens expressed by Purkinje neurons and cancer cells – the cerebellar degeneration-related (CDR) proteins CDR1, CDR2, and CDR2-like (CDR2L) [5,6,7,8]. The interactions between anti-Yo and CDR proteins have been shown to cause Purkinje neuron death [9, 10].
We demonstrated recently that CDR2L is likely the major antigen of anti-Yo [11]. Increased expression of CDR2L mRNA has been observed in ovarian tumours from PCD patients with anti-Yo antibodies compared to ovarian tumours from patients without PCD or anti-Yo antibodies [12]. Due to our limited understanding of the biological functions of the CDR2L protein, the consequences of this increased expression in tumours are not known. Previous studies showed that CDR2L localizes to the cell cytoplasm in association with membrane-bound and free ribosomes in both Purkinje neurons and cancer cells [13,14,15]. Further, we demonstrated that CDR2L interacts with several ribosomal proteins, suggesting a role in protein synthesis [15].
In the present study, we used a multi-omics approach to explore changes in the transcriptome, proteome, and secretome induced by knockout of CDR1, CDR2, and CDR2L in ovarian cancer cells. Since CDR2L is considered the major antigen in Yo-associated PCD, we characterized the effects of CDR2L knockout in detail. We found that loss of CDR2L had a distinct effect on the transcriptome and proteome compared to loss of CDR1 or of CDR2. Following knockout of CDR2L, genes associated with ribosome biogenesis, protein translation and cell cycle progression were downregulated, and cell proliferation was impaired. This study provides an extensive analysis of CDR-knockout cells and offers a valuable resource for future investigations into the CDR proteins.
Results
Knockout of genes encoding CDR proteins induces expression changes in the transcriptome, proteome, and secretome of ovarian cancer cells
Using differential expression analysis, we investigated changes induced by knockout of CDR1, CDR2, and CDR2L in the transcriptome, and the proteomes of cell lysate and conditioned cell medium. These three omics datasets were separately normalized and subjected to pre-filtering, resulting in 18,719 genes in the transcriptomic dataset (Table S1): the majority of which were protein-coding (80.2%). There were 7,074 proteins in the cell lysate proteomics dataset (referred to as the proteome; Table S2) and 3,588 in the cell medium proteomics dataset (referred to as the secretome; Table S3). Comparison of the transcriptomics and proteomics data revealed that 6,881 (97.3%) and 3,470 (96.7%) proteins in the proteome and secretome datasets, respectively, were matched with protein-coding mRNAs (Fig. 1A). There were significant positive correlations between the transcriptome and proteome (Spearman’s rank correlation coefficient across wild-type and knockout cell lines, ρ = 0.53, p-value < 2.2e-16), between transcriptome and secretome (ρ = 0.38, p-value < 2.2e-16), and between proteome and secretome (ρ = 0.45, p-value < 2.2e-16; Figure S1).
Proteins present in the secretome include secreted proteins, proteins shed by membrane vesicles, and proteins leaked from dying cells. To assess the sources of proteins identified in the secretome, we compared these proteins with proteins in the Plasma Proteome Database [16], which are reportedly present in the plasma and serum, and with those in the Vesiclepedia Database [17], which have been detected in extracellular vesicles. Of the 3,588 proteins in the secretome, 55.8% were found in the Plasma Proteome Database and 86.3% were found in the Vesiclepedia database, suggesting that most of proteins in the secretome were derived from extracellular vesicles.
We performed differential expression analysis on each dataset comparing knockout cells to the wild-type (WT) cells, resulting in differentially expressed (DE) genes and proteins in the transcriptome, proteome, and secretome for each knockout cell line (Table 1). Effect sizes and significance values are listed in Tables S4-S6. There were more DE genes in the transcriptome in all group comparisons than in the proteome. Since mass spectrometry-based proteomics is less sensitive than RNA sequencing, resulting in fewer detected proteins than genes, the lower number of DE genes in the proteome compared to the transcriptome was not unexpected. The proportion of genes that were DE was similar between the transcriptome (mean proportion across all three knockout cell lines = 32.0%) and proteome (35.7%), whereas the secretome had a much higher proportion of DE genes (58.5%) with more than half of secreted proteins affected by the knockout, irrespective of CDR protein (Fig. 1B).
Some genes were DE in all three knockout cell lines. These constituted 15.7% of DE genes in the transcriptome and 15.4% of DE proteins in the proteome (FDR < 0.05 in all three knockout analyses, where the intersection was performed irrespective of direction of change). In the secretome, this proportion was more than twice as high with 34.4% of DE genes common to the three knockout cell lines (Fig. 1C). Over-representation analysis on these common DE genes was performed to identify the biological processes affected by the loss of any of the CDR genes (Table S7). In the transcriptome, cell adhesion, cell motility, extracellular structure organization, and negative regulation of nervous system development were among the enriched pathways. In the proteome, several metabolic processes were enriched, including organic hydroxy compound metabolic process, steroid metabolic process, and fatty acid metabolic process. In the secretome, ribosomes, RNA catabolic process, and protein localization to endoplasmic reticulum were enriched.
Knockout of CDR2L results in a unique transcriptomic and proteomic expression profile
To understand the differential effect of the knockout, we compared the number of DE genes between the three knockout cell lines. In the transcriptome, 43.4% of genes were DE in CDR2L-knockout cells compared to 22.8% in CDR2- and 30.0% in CDR1-knockout cells (Fig. 1C). In line with the higher proportion of DE genes, CDR2L also showed the highest number of uniquely altered genes (i.e., not significantly DE in the other knockout cell lines) with 29.0% of genes in the transcriptome uniquely altered in the CDR2L-knockout cells (Fig. 1C). Similarly, in the proteome, 44.6% of genes were DE in CDR2L-knockout cells, with 27.3% of genes in the proteome uniquely altered in CDR2L-knockout cells. The high number of uniquely altered genes in CDR2L-knockout cells suggests that loss of CDR2L has a distinct effect on the cells compared to loss of the other CDR proteins. In the secretome, however, most DE genes were common to all knockout cell lines, with CDR2L-knockout cells having the lowest number of uniquely altered genes.
Using principal component analysis on the three omic datasets separately, we observed a clustering of samples into the respective knockout groups, confirming a good accordance between biological replicates. The CDR2L-knockout cells were separated from the other cells, including WT, along the first principal component in the transcriptome and proteome but not in the secretome (Fig. 1D), further strengthening our hypothesis that loss of CDR2L has a unique effect inside the cell. We investigated this further using linear regression and found a significant association between the first principal component (representing the main axis of variation in the dataset) and the condition variable CDR2L versus the other cell lines for the transcriptome (p-value = 3.4e-04) and the proteome (p-value = 4.3e-08) but not for the secretome (p-value = 0.357). When comparing the 1000 genes with the highest absolute loadings in the first principal component, only 5.5% were common between the transcriptome and proteome, suggesting that the genes driving the separation between the CDR2L-knockout cell line and the other cell lines are different between the transcriptome and proteome. Based on these findings, we decided to further characterize the likely unique intracellular expression profile of the CDR2L-knockout cells.
Gene set enrichment analysis reveals discordant changes in enriched pathways at transcript and protein level upon knockout of CDR2L
To better understand the biological functions of the genes with expression altered by knockout of CDR2L, we performed gene set enrichment analysis. There were both common pathways and pathways specifically enriched in the transcriptome, proteome, and secretome. At the transcript level, the most significantly upregulated pathways were related to ribosomes, DNA replication, and RNA processing (Fig. 2A and Table S8). The few downregulated pathways included cell adhesion and calcium binding (Fig. 2A and Table S8). In the proteome, downregulated pathways were mostly related to the cell cycle, RNA processing, and ribosome biogenesis, while upregulated pathways were enriched for nucleoside metabolism, immune responses, homeostasis, and regulation of actin filaments (Fig. 2B and Table S8). Downregulated pathways in the secretome were mainly related to ribosomes, RNA processing, and amide biosynthesis, and upregulated pathways were related to adhesion, proteoglycan metabolism, lysosomes, and synapse assembly (Fig. 2C and Table S8).
Interestingly, several of the pathways upregulated in the transcriptome were downregulated in the proteome and/or secretome (Fig. 2D): These pathways were related to RNA processing, DNA replication, and ribosome biogenesis. On the other hand, cell adhesion and calcium binding were downregulated in the transcriptome but upregulated in the secretome. Enriched pathways in the proteome and secretome showed the same direction of change.
Ribosome biogenesis factors are dysregulated in the transcriptome and proteome of CDR2L-knockout cells, but in opposite directions
CDR2L has previously been shown to localize to ribosomes and interact with ribosomal proteins [15]. As gene set analysis of CDR2L-knockout cells showed enrichment of ribosome biogenesis at the transcriptome and proteome level, we evaluated the expression of all genes with the GO term “ribosome biogenesis” significantly DE in at least one of the datasets (n = 206). In the transcriptome, the majority (80.4%) of DE mRNAs (n = 168) involved in ribosome biogenesis were upregulated, while in the proteome, the majority (84.2%) of DE proteins (n = 120) were downregulated. Eighty-one genes were DE in both transcriptome and proteome. Ribosome biogenesis is a complex process involving several highly orchestrated steps. To determine which steps were affected by the loss of CDR2L, we classified the genes according to a recent review of ribosome biogenesis factors [18]. In both the transcriptome and the proteome, most of the DE genes were associated with the small subunit processome, maturation of the 60S subunit, and ribosomal RNA (rRNA) processing (Fig. 3). That the same steps of ribosome biogenesis are affected in both transcriptome and proteome, albeit in opposite directions, strengthens the hypothesis that CDR2L is intimately involved in ribosomal function.
Dysregulated ribosome biogenesis interferes with the assembly of functional ribosomes and consequently inhibits protein translation. Therefore, we investigated the expression of individual ribosomal proteins and eukaryotic initiation factors (eIFs) that are essential for initiation of translation. In the transcriptome of CDR2L-knockout cells, genes encoding 39 ribosomal proteins were DE; of these 24 were downregulated. In the proteome and secretome, 37 and 53 ribosomal proteins were DE, respectively, of which all were downregulated. We also observed dysregulation of eIFs. A number were upregulated in the transcriptome, including EIF3B, EIF4A3, EIF5A, EIF5B, and EIF6, while EIF4G3, EIF4B, and EIF3K were downregulated. In the proteome, eIFs EIF1, EIF3L, EIF4B and EIF4G3 were downregulated. In the secretome, all DE eIFs were downregulated. These changes suggests that loss of CDR2L disrupts protein synthesis in the cells.
Knockout of CDR2L causes downregulation of key regulators of the cell cycle
Gene set enrichment analysis of CDR2L-knockout cells resulted in more enriched pathways for the proteome (n = 325) than for the transcriptome (n = 61) or secretome (n = 139). The large number of pathways in the proteome rendered interpretation difficult. Therefore, we created a term network to investigate the functional relationships between these pathways. Most of the identified functional groups were related to the cell cycle, including cell cycle phase transition, mitotic spindle organization, chromatid segregation, and nuclear division (Fig. 4A). Other enriched pathways were related to rRNA processing and nucleoside metabolism – processes that are tightly connected to the cell cycle. Pathways related to the cell cycle and rRNA processing were all downregulated, whereas pathways related to nucleoside metabolism were upregulated in cells lacking CDR2L.
In relation to this finding, analysis of individual genes showed that knockout of CDR2L caused altered expression of multiple important regulators of the cell cycle in the proteome and/or transcriptome, including genes encoding Aurora kinases AURKA and AURKB, checkpoint regulators (BUB1, BUB1B, CHEK1, and RB1), cyclins (CCNA2, CCNB1, and CCNB2), cyclin-dependent kinases (CDK2, CDK4, and CDK6), cell division cycle proteins (CDC7, CDC20, CDC23, and CDC27), and cyclin-dependent kinase inhibitors (CDKN1A, CDKN1B, and CDKN2A). Other dysregulated proteins included components of the minichromosome maintenance protein complex, which is essential for DNA replication (MCM2, MCM3, MCM4, MCM5, and MCM7), subunits of the condensin complex, which plays a central role in chromosome condensation (NCAPD2, NCAPD3, NCAPG, NCAPG2, and NCAPH), centromere proteins, which are involved in chromosome segregation (CENPE, CENPF, CENPM, and CENPU), and several transcription factors (E2F, FOS, JUN, and MYC). Similar to ribosome biogenesis factors, several of the cell cycle-related genes were upregulated in the transcriptome but downregulated in the proteome (Fig. 4B). Proapoptotic factors and key mediators of apoptosis like BAD, BAX, BID, and PYCARD were upregulated in the CDR2L-knockout proteome. Taken together, these data suggest that loss of CDR2L suppresses cell cycle progression in ovarian cancer, which may trigger apoptosis.
Utilizing transcription factor enrichment analysis to identify the likely contribution of transcription factors to the changes in the transcriptome following loss of CDR2L, we identified several transcription factors which regulate the cell cycle and proliferation (Table S9). These included MYC, ERF, PA2G4, and subunits of the AP-1 transcription factor (FOSL1, FOSL2, and JUN) which regulate AURKB, BAX and several eIFs, among others, and members and interaction partners of the E2F family of transcription factors (E2F1, E2F4, E2F7, and TFDP1) which regulate CDK4 and CDK6. FOSL1 showed a particularly large upregulation in the transcriptome (log2 fold change = 4.4).
CDR2L interacts with ribosome biogenesis factors and regulators of the cell cycle
In our earlier study, we used co-immunoprecipitation coupled with mass spectrometry to identify potential interaction partners of CDR2L in OVCAR-3 cells [15]. In light of our current findings, we re-examined our unpublished results and found that several of the dysregulated ribosome biogenesis factors were among the identified interaction partners of CDR2L including HEATR1, GTPBP4, PA2G4, PDCD11, SBDS, and UTP20, nucleolar proteins (NOP56 and NOP58), RNA helicases (DDX17, DDX21, DDX3X, and DHX29), and WD repeat-containing proteins (WDR18, WDR36, and WDR43). Several eIFs (e.g., EIF4G2, EIF2A, EIF3L, EIF6), cyclin-dependent kinases (CDK4 and CDK6), and components of the minichromosome maintenance protein complex (MCM2, MCM3, MCM4, MCM5, and MCM7) were also identified as potential interaction partners of CDR2L (Table S10). This strengthens the relationship between CDR2L and important regulators of ribosome biogenesis, protein synthesis and cell cycle progression.
Ribosome-associated genes are key hubs in the network of genes differentially expressed in both transcriptome and proteome of CDR2L-knockout cells
We identified 1958 genes that were altered in both the transcriptome and the proteome of the CDR2L-knockout cells (Fig. 5A). Although most of these genes showed a consistent pattern of regulation, that is they were upregulated or downregulated in both datasets, 31.2% of the genes were oppositely altered in the transcriptome versus the proteome (Fig. 5B). This proportion was twice as high in the CDR2L-knockout cells as in the other CDR-knockout cells (Figure S2).
The genes with altered expression in both transcriptome and proteome of CDR2L-knockout cells were used to create a protein-protein interaction network using the STRING database. It consisted of a total 1958 nodes and 2745 edges. The genes in the network were enriched for many essential cellular processes including organelle organization, cell cycle and cell division, ribonucleoprotein complex biogenesis, gene expression, and intracellular transport (Table S11). By ranking the genes using the maximal clique centrality algorithm (see methods), we identified key hub genes with high connectivity within the network. The top 10 hub genes were all ribosome-associated genes including small ribosomal subunit proteins (RPS7 and RPS9), components of the small subunit processome (UTP4, UTP11, UTP20, HEATR1, WDR46, NOP56, and KRR1), and a nucleolar complex protein (NOC4L), again highlighting the effect of CDR2L knockout on ribosomal pathways and their interaction partners both at transcript and protein levels.
Next, to identify clusters of tightly connected genes within the network, we performed Markov clustering analysis followed by enrichment analysis of the genes in the clusters. Three of the clusters showed enrichment for pathways related to ribosome biogenesis, RNA processing and translation: clusters 1, 3, and 4 (Fig. 6). The largest cluster (cluster 1) contained several components of the small subunit processome (UTP4, UTP11, UTP20, UTP25, NOC4L, NOP56, KRR1, and WDR46), RNA helicases (DDX47 and DHX37), translation initiation factor 3 subunits (EIF3I, EIF3G, EIF3C, and EIF3L), and small ribosomal subunit proteins (RPSA, RPS4X, RPS5, RPS9, RPS11, RPS16, and RPS23). Other clusters were distinctly enriched for mitochondrial pathways. Cluster 2 contained several subunits of complex 1 of the mitochondrial respiratory chain, which were downregulated in both transcriptome and proteome. Cluster 5 contained several mitochondrial ribosomal proteins, which were mostly upregulated in both datasets. Cluster 7 was enriched for protein folding. Taken together, these data indicate that loss of CDR2L results in dysregulation of protein synthesis and mitochondrial function.
CDR2L knockout suppressed proliferation of ovarian cancer cells
To study the effect of CDR2L knockout on proliferation rates in vitro, we performed live cell imaging of WT cells and CDR2L-knockout cells. The WT cells reached 90% confluence after approximately 40 h, whereas CDR2L-knockout cells did not reach the same level of confluence until after approximately 65 h (Fig. 7). This demonstrates that CDR2L knockout markedly impaired cell proliferation in OVCAR-3 cells.
Discussion
The pathological mechanisms underlying PCD are largely unknown. Tumour expression and genetic alterations of the CDR proteins, particularly CDR2L, is thought to trigger an immune response which targets both the cancer cells and the Purkinje neurons endogenously expressing the CDR proteins [19]. However, the mechanisms of immune activation and neural injury induced by these intracellular neural proteins remain unknown. This is in part due to our limited understanding of the biological properties of the CDR proteins. Here, we used transcriptomics and proteomics to evaluate the effects of CRISPR/Cas9-mediated knockout of CDR-encoding genes in ovarian cancer cells. As the CDR2L protein is likely to be the major target of the onconeural anti-Yo antibody in PCD [11], we performed an in-depth characterization of the CDR2L-knockout cells. Previously, we demonstrated that CDR2L colocalizes with ribosomes in ovarian cancer cells and Purkinje neurons [15]. Here, we show that loss of CDR2L causes dysregulation of ribosome biogenesis and reduces cell proliferation.
Changes in mRNA levels are generally assumed to be reflected by similar changes in the abundances of the corresponding proteins [20]. However, we found only a moderate correlation between protein abundances and mRNA levels, in line with previous studies [21, 22]. Therefore, integrating transcriptomic data with proteomic data provided us with a more comprehensive view of the cellular changes induced by knockout of CDR2L. Interestingly, we found that nearly one third of genes with expression altered at both transcriptomic and proteomic levels showed opposite trends of regulation, with the majority upregulated in the transcriptome and downregulated in the proteome. The proportion of DE genes with opposite regulation was twice as high in the CDR2L-knockout cells as in the other knockout cells, suggesting that CDR2L specifically affects the balance of mRNAs and proteins in the cell. Multiple processes contribute to this balance, including modulation of the translation rate by RNA-binding proteins, degradation of transcripts through binding of microRNAs, protein degradation through the ubiquitin-proteasome pathway or autophagy, temporal delay of protein synthesis, and availability of ribosomes in controlling the abundance of proteins [23]. Our work confirms that transcript abundance is not the only determinant of protein abundance, and highlights the importance of post-transcriptional, translational, and degradation regulation in controlling the abundance of proteins.
Early studies using serum-derived antibodies from PCD patients localized the Yo antigen to ribosomes and rough endoplasmic reticulum [13, 14], thereby providing the first indications that CDR2L, the major antigen of anti-Yo antibodies [11], was a regulator of protein synthesis. Using both Yo antibodies and commercial CDR2L antibodies, we confirmed the localization of CDR2L to ribosomes and identified the ribosomal protein RPS6 as a potential interaction partner of CDR2L [15]. In the present study, we found that knockout of CDR2L caused dysregulation of expression of multiple ribosome biogenesis factors, RNA processing factors, ribosomal proteins, and eukaryotic translation initiation factors. The majority of these were upregulated in the transcriptome and downregulated in the proteome and secretome. Downregulation of ribosome biogenesis factors, ribosomal proteins, and translation initiation factors should lead to reduced protein synthesis, with global effects favouring decreased protein abundances. In CDR2L-knockout cells, we identified similar numbers of increased and decreased mRNAs, but 1.5-fold more downregulated than upregulated proteins.
Upregulation of mRNAs involved in ribosome biogenesis, rRNA processing, and protein translation could be a consequence of attempts by cells to compensate for loss of ribosomal function and reduced protein synthesis. In support of this hypothesis, Myc, a proto-oncogenic transcription factor and an important regulator of ribosome biogenesis and cell proliferation [24], was upregulated in the transcriptome of CDR2L-knockout cells. Interestingly, previous studies have shown that CDR2 and c-Myc interact to synergistically regulate c-Myc-dependent transcription during passage through mitosis, and CDR2 was found to interact with and downregulate c-Myc function [25, 26]. As PCD sera blocked this interaction, it was suggested that increased c-Myc activity could contribute to the degeneration of Purkinje neurons seen in PCD. As these studies used a combination of antibodies from PCD patients and commercial anti-CDR2 antibodies, it is possible that some of these findings can be attributed to CDR2L. Further investigations into the CDR2L-Myc relationship are therefore warranted.
Several ribosome biogenesis factors, ribosomal proteins, and eukaryotic translation initiation factors were identified in the secretome of the OVCAR-3 cells, many of which were downregulated in the CDR2L-knockout cells. Although generally considered intracellular proteins, these can be found in extracellular vesicles (EVs) [27]. Ribosomal proteins are selectively incorporated into EVs and released from both normal and malignant cells [28]. The amounts and types of ribosomal proteins in EVs are altered under pathological conditions. For example, Dabbah et al. reported increased levels and repertoire of ribosomal proteins in EVs from mesenchymal stem cells in multiple myeloma compared with EVs from healthy controls [29]. These EVs were internalized by neighbouring cells and promoted their proliferation – a processes that was dependent on the ribosomal content of the EVs. The dysregulated ribosome biogenesis in CDR2L-knockout cells may therefore not only disrupt the activity of intracellular ribosomes but may also affect the levels and repertoire of ribosomal proteins released from cells in EVs, which in turn may affect the behaviour of neighbouring cells.
Ribosome biogenesis influences the cell cycle as it regulates cell size and growth. Disruption of ribosome biogenesis causes ribosomal stress, which results in p53-dependent cell cycle arrest and apoptosis [30]. In the present study, we showed that multiple pathways related to the cell cycle were downregulated in CDR2L-knockout cells. Important regulators of the cell cycle, such as cyclins (CCNA2, CCNB1, and CCNB2), cyclin-dependent kinases (CDK2, CDK4, and CDK6), and checkpoint regulators (BUB1, BUB1B, and CHEK1), were downregulated, whereas proapoptotic factors (BAD, BAX, and BID) were upregulated. CDK6, which together with CDK4 and cyclin D facilitate the progression of cells through the early G1 phase of the cell cycle [31], was highly downregulated (log2 fold change = -3.5). In line with downregulation of cell cycle proteins, knockout of CDR2L impaired cell proliferation in vitro. Similarly, reduction of levels of the ribosomal protein RPS6 in ovarian cancer cells caused downregulation of CDK2, CDK4, CDK6, cyclin E, and cyclin D1, thereby blocking the transition from G0/G1 to S phase and suppressing cell proliferation [32]. Previously, we showed that CDR2L interacts with RPS6 [15]. The observation that loss of either CDR2L or RPS6 causes similar effects on the cell cycle in ovarian cancer cells further supports an association between the two proteins. However, the effect of knockout of CDR2L on the cell cycle seems to be independent of changes in RPS6 expression as RPS6 was not differentially expressed following the knockout. Rather, knockout of CDR2L may alter the function of RPS6 or disrupt common downstream signalling events.
The current study provides an extensive bioinformatic analysis of the effect of CDR-knockout which can serve as a valuable resource for designing future experiments to determine the biological properties of the CDR proteins. For example, the effect on the cell cycle should be further explored using flow cytometry which would provide insight into which phases of the cell cycle is affected by the absence of the CDR2L protein. Rescue experiments in which CDR2L is re-expressed would provide further validation of the observed knockout phenotype. Comparative studies with all CDR-knockout cells would determine whether the effect on the cell cycle is unique to CDR2L, or whether it is a shared effect among the CDR proteins. Further, an orthotopic xenograft model could be used to explore the effects of CDR-knockout on the tumorigenic potential of ovarian cancer cells, including assessment of tumour growth rate and metastasis.
Conclusions
In conclusion, we found that knockout of CDR2L in an ovarian cancer line dysregulates genes involved in ribosome biogenesis, protein synthesis and cell cycle-related processes, ultimately impairing cell proliferation. Yo antibodies from PCD patients are internalized by Purkinje neurons where they bind to CDR2L [9, 11]. Whether binding of anti-Yo to CDR2L interferes with CDR2L and thereby disrupts ribosome homeostasis or protein synthesis in Purkinje neurons remains to be explored. Deficits in ribosome biogenesis in neurons have been shown to cause dendritic degeneration and loss of synaptic plasticity [33]. The data provided by this study serves as a valuable resource for exploring the biological properties of the CDR proteins.
Materials and methods
Generation of knockout OVCAR-3 cells
CDR1-knockout, CDR2-knockout, and CDR2L-knockout OVCAR-3 cell lines and WT (#YC-D019) cells were provided by Ubigene Biosciences Co., Ltd. The single guide RNAs (sgRNAs) were designed using the online CRISPR design tool (https://en.rc-crispr.com/). The pair of oligonucleotides corresponding to each sgRNA were annealed and ligated into the YKO-RP003 vector (Ubigene Biosciences Co., Ltd.). The YKO-RP003-[sgRNA] plasmids were transfected into cells with Lipofectamine 3000 (Thermo Fisher Scientific). At 24–48 h after the transfection, puromycin was added. After antibiotic selection, cells were diluted by limited dilution methods and inoculated into 96-well plate. Selection of single clones was performed after 2–4 weeks, and knockout was validated by polymerase chain reaction and Sanger sequencing.
Cell culture
OVCAR-3 cells were cultured in RPMI 1640 medium supplemented with 20% fetal bovine serum (Thermo Fisher Scientific) and 1% penicillin/streptomycin (Thermo Fisher Scientific) in a humidified atmosphere with 5% CO2 at 37 °C. Cells were used for experiments within 10 passages of thawing. Cells were routinely tested for mycoplasma (MycoAlert PLUS Mycoplasma Detection Kit).
Cell proliferation assay
WT and CDR2L-knockout OVCAR-3 cells were seeded at 10,000 cells/well in Incucyte Imagelock 96-well plates (Sartorius). Following cell seeding, plates were placed in an Incucyte live-cell analysis system (Sartorius) and imaged every 2 h for 96 h. The rate of cell proliferation was determined based on changes in confluency of each imaged well. Analysis was performed using the Incucyte S3 software (Sartorius). The experiment was repeated three times with similar results.
RNA sequencing
Total RNA was extracted from all four cell lines in biological triplicates using miRNeasy Tissue/Cells Advanced (Qiagen) according to the manufacturer’s protocol. Library preparation and sequencing was performed by BGI Genomics (https://www.bgi.com/global/). Briefly, RNA quality and concentration was assessed using the Agilent 2100 Bioanalyzer. mRNA was enriched using oligo (dT)-functionalized magnetic beads, fragmented, converted to double-stranded cDNA, and amplified using PCR. Samples were sequenced through combinatorial probe-anchor synthesis on the DNBSEQ sequencing system using paired-end 100-bp sequencing. Raw reads were filtered using SOAPnuke version 2.1.3. Reads containing the adaptor, reads with N content greater than 5%, and low-quality reads were removed.
RNA expression quantification
Salmon version 1.10.2 [34] was used to quantify the abundance at the transcript level against the Ensembl release 75 transcriptome. Transcript-level quantification was collapsed onto gene-level quantification using the tximport R package version 1.26.1 [35]. Genes with zero counts and zero variance, genes without gene symbol, and mitochondria-encoded genes were excluded. Lowly expressed genes were excluded by only retaining genes where the counts across samples in the 80th percentile were greater than 0.1 counts per million. Differential gene expression analysis was performed using the DESeq2 R package version 1.32.2 [36] with default parameters. Multiple hypothesis testing was performed with the default automatic filtering of DESeq2 followed by FDR calculation by the Benjamini-Hochberg procedure.
Liquid chromatography and mass spectrometry analysis
Whole cell lysate and conditioned cell media (CCM) was prepared from all four cell lines in biological triplicates and analyzed by mass spectrometry-based proteomics. Cells were lysed with RIPA lysis buffer (G-Biosciences) supplemented with protease inhibitors (Thermo Fisher Scientific) on ice for 15 min. Protein lysates were collected by centrifugation at 14,000 x g for 15 min at 4 °C. Protein concentration was measured with the BCA assay (Thermo Fisher Scientific). For CCM, cells cultured in T175 flasks (Corning) were incubated with media without any additives for 24 h. The CCM was collected, centrifuged at 3000 x g for 15 min and filtered through a 0.2-µm syringe filter (Thermo Fisher Scientific). The CCM was concentrated using Ultra-15 3-kD cut-off centrifugal filters (Millipore) and centrifugation at 3000 x g for 1 h.
The samples were mixed with paramagnetic beads (Sera-Mag Speed beads, GE Healthcare) in a 1:10 ratio. A 100% ethanol to 70% ethanol solution was added to the peptide-bead mixture, followed by agitation at 1000 rpm for 7 min at room temperature. The beads were washed twice in 80% ethanol to remove the lysis buffer, and 50 µl digestion buffer containing 100 mM AmBic, 1 mM CaCl2, and 0.2 µg/µl trypsin was added to each sample, followed by sonication for 30 s. Samples were incubated for 16 h at 37 °C in a thermomixer at 1,000 rpm. After digestion, the samples were centrifuged at 13,000 rpm at 24 °C for 3 min. Following magnetic separation, the supernatants were transferred to a fresh tube. NaCl (0.5 M) was added to the magnetic beads, the sample was sonicated for 30 s and centrifuged a second time at 13,000 rpm at 24 °C for 3 min. The supernatants of these tubes were added to the tubes with the previous supernatant. Trifluoroacetic acid (TFA, 0.1%, 200 µl) was added to each sample.
Oasis 96-well cartridges (Waters) were used for desalting. Cartridges were activated by adding 500 µl of 80% acetonitrile (ACN), 0.1% formic acid (FA) and centrifuged at 200 x g for 1 min. The cartridges were then washed three times with 0.1% TFA. After adding the samples, the cartridges were centrifuged at 100 x g for 3 min. The flow through was discarded, and cartridges were washed twice with 0.1% TFA followed by a short centrifugation. Samples were eluted with 100 µl of 80% ACN, 0.1% FA. The samples were freeze dried prior to tandem mass tag labelling.
Approximately 0.5 µg protein as tryptic peptides dissolved in 2% ACN and 0.5% FA, were injected into an Ultimate 3000 RSLC system (Thermo Fisher Scientific) connected online to Orbitrap Eclipse mass spectrometer (Thermo Fisher Scientific) equipped with EASY-spray nano-electrospray ion source (Thermo Fisher Scientific). For trapping and desalting process, the samples were loaded and desalted on a pre-column (Acclaim PepMap 100, 2 cm × 75 μm ID nanoViper column, packed with 3 μm C18 beads) at a flow rate of 5 µl/min for 5 min with 0.1% TFA. Peptides were separated during a biphasic ACN gradient from two nanoflow UPLC pumps (flow rate of 250 nl/min) on a 25-cm analytical column (PepMap RSLC, 50 cm × 75 μm ID EASY-Spray column, packed with 2-µm C18 beads). Solvent A and B were 0.1% FA (vol/vol) in distilled H2O and 100% ACN, respectively. The gradient composition was 5% B during trapping (5 min) followed by 5–7% B over 30 s, 8–22% B for the next 145 min, 22–28% B over 16 min, and 35–80% B over 15 min. Elution of very hydrophobic peptides and conditioning of the column were performed during a 15 min isocratic elution with 90% B and a 20 min isocratic elution with 5% B, respectively. The peptides eluted from the LC-column were ionized in the electrospray and analyzed by the Orbitrap Eclipse. The mass spectrometer was operated in the data-dependent-acquisition mode to automatically switch between full scan MS and MS/MS acquisition. Instrument control was through Tune 2.7.0 and Xcalibur 4.4.16.14.
Survey full-scan MS spectra (from m/z 375 to 1500) were acquired in the Orbitrap with resolution R = 120,000 at m/z 200 after accumulation to a target value of 4e5 in the C-trap, ion accumulation time was set as auto. FAIMS was enabled using two compensation voltages (CVs), -45 V and − 65 V respectively. During each CV, the mass spectrometer was operated in data-dependent-acquisition mode to automatically switch between full scan MS and MS/MS acquisition. The cycle time was maintained at 0.9 s/CV. The most intense peptides with charge states 2 to 6 were sequentially isolated to a target value (AGC) of 2e5 and maximum IT of 120 ms in the C-trap, and isolation width maintained at 0.7 m/z. Fragmentation was performed with a normalized collision energy of 30%, and fragments were detected in the Orbitrap at a resolution of 30,000 at m/z 200, with first mass fixed at m/z 110. The spray and ion-source parameters were as follows: ion spray voltage = 1900 V, no sheath and auxiliary gas flow, and capillary temperature of 275 °C. The raw files were analyzed using Proteome Discoverer 2.4 (Thermo Fisher Scientific). Spectra were matched against the Homo sapiens database obtained from UniProt. The final protein intensity files were exported as txt files and used for differential expression analysis.
Protein expression quantification
Protein intensity txt files were analyzed for differentially expressed proteins using R version 4.2.2. First, proteins labelled as “Contaminant” by Proteome Discoverer were removed from the analysis. Second, to be included for downstream analysis, proteins were required to have non-zero intensity in all replicates in at least one experimental group. The filtered proteomic dataset was composed of 7447 proteins. Missing values (n = 1830) were imputed using the missForest function from the missForest R package version 1.5 [37], which applies a random forest imputation algorithm shown to perform well with proteomics data [38]. Additionally, we removed proteins with a median abundance value less than the 5th percentile. To test for differential expression we used the lmFit and eBayes functions from the limma R package version 3.54.1 [39].
Pathway enrichment analysis
Gene set enrichment analysis was performed using the fgsea R package version 1.24.0 [40]. The specific parameters used are documented in the code for the analysis. Gene sets derived from the GO [41] and the KEGG [42] databases were downloaded from the Molecular Signatures Database [43]. Gene sets with an FDR adjusted p-value < 0.05 were considered significant. Over-representation analysis of differentially expressed genes and proteins was performed using the WebGestaltR R package version 0.4.4 [44], with GO and KEGG gene sets. Gene sets with an FDR adjusted p-value < 0.05 were considered significant.
Protein-protein interaction network
The protein-protein interaction network was generated based on the interactions derived from the STRING database using the Cytoscape STRING app version 2.0.1 [45] with a strict confidence score cutoff of 0.8 and network type set to “physical subnetwork”. The network was visualized using Cytoscape version 3.9.1 [46], and nodes in the network were coloured according to their log2 fold change using the Omics Visualizer app version 1.3.0 [47]. Key hub genes were identified using the maximal clique centrality algorithm in the CytoHubba app version 0.1 [48]. To cluster the network based on their interactions from STRING, we used the clusterMaker2 version 2.3.4 [49] and performed Markov clustering with default parameters. Enrichment analysis of the genes in the clusters was performed using the Cytoscape STRING app. Gene sets with an FDR adjusted p-value < 0.05 were considered significant.
Pathway term network
Enriched pathways were grouped into functionally organized networks using ClueGO version 2.5.10 [50] and Cluepedia version 1.5.10 [51]. The networks were visualized using Cytoscape.
Transcription factor enrichment analysis
Transcription factor enrichment analysis was performed by analyzing the differentially expressed mRNAs (FDR < 0.05) using the ChIP-X Enrichment Analysis version 3 (ChEA3) [52] with default parameters.
Data availability
All data generated or analysed during this study are included in this published article and its supplementary files. The source code and all necessary files to perform the analyses are available on GitHub under the GPL public license v3.0: https://github.com/ETSunhome.
References
Graus F, Vogrig A, Muñiz-Castrillo S, Antoine JG, Desestret V, Dubey D et al. Updated Diagnostic Criteria for Paraneoplastic Neurologic Syndromes. Neurology(R) neuroimmunology & neuroinflammation. 2021;8(4). Epub 2021/05/20. https://doi.org/10.1212/nxi.0000000000001014. PubMed PMID: 34006622; PubMed Central PMCID: PMCPMC8237398.
Shams’ili S, Grefkens J, de Leeuw B, van den Bent M, Hooijkaas H, van der Holt B, et al. Paraneoplastic cerebellar degeneration associated with antineuronal antibodies: analysis of 50 patients. Brain. 2003;126(Pt 6):1409–18. https://doi.org/10.1093/brain/awg133. Epub 2003/05/24.
Vogrig A, Gigli GL, Segatti S, Corazza E, Marini A, Bernardini A, et al. Epidemiology of paraneoplastic neurological syndromes: a population-based study. J Neurol. 2020;267(1):26–35. https://doi.org/10.1007/s00415-019-09544-1. Epub 2019/09/26.
Vedeler CA, Antoine JC, Giometto B, Graus F, Grisold W, Hart IK et al. Management of paraneoplastic neurological syndromes: report of an EFNS Task Force. European journal of neurology. 2006;13(7):682 – 90. Epub 2006/07/13. https://doi.org/10.1111/j.1468-1331.2006.01266.x. PubMed PMID: 16834698.
Furneaux HM, Dropcho EJ, Barbut D, Chen YT, Rosenblum MK, Old LJ, et al. Characterization of a cDNA encoding a 34-kDa Purkinje neuron protein recognized by sera from patients with paraneoplastic cerebellar degeneration. Proc Natl Acad Sci USA. 1989;86(8):2873–7. https://doi.org/10.1073/pnas.86.8.2873. Epub 1989/04/01.
Cunningham J, Graus F, Anderson N, Posner JB. Partial characterization of the Purkinje cell antigens in paraneoplastic cerebellar degeneration. Neurology. 1986;36(9):1163–8. https://doi.org/10.1212/wnl.36.9.1163. Epub 1986/09/01.
Fathallah-Shaykh H, Wolf S, Wong E, Posner JB, Furneaux HM. Cloning of a leucine-zipper protein recognized by the sera of patients with antibody-associated paraneoplastic cerebellar degeneration. Proc Natl Acad Sci USA. 1991;88(8):3451–4. https://doi.org/10.1073/pnas.88.8.3451. Epub 1991/04/15.
Eichler TW, Totland C, Haugen M, Qvale TH, Mazengia K, Storstein A, et al. CDR2L antibodies: a New Player in Paraneoplastic Cerebellar Degeneration. PLoS ONE. 2013;8(6):e66002. https://doi.org/10.1371/journal.pone.0066002. Epub 2013/07/05.
Schubert M, Panja D, Haugen M, Bramham CR, Vedeler CA. Paraneoplastic CDR2 and CDR2L antibodies affect Purkinje cell calcium homeostasis. Acta Neuropathol. 2014;128(6):835–52. https://doi.org/10.1007/s00401-014-1351-6. Epub 2014/10/25.
Greenlee JE, Clawson SA, Hill KE, Wood BL, Tsunoda I, Carlson NG. Purkinje cell death after uptake of anti-yo antibodies in cerebellar slice cultures. J Neuropathol Exp Neurol. 2010;69(10):997–1007. https://doi.org/10.1097/NEN.0b013e3181f0c82b. Epub 2010/09/15.
Kråkenes T, Herdlevaer I, Raspotnig M, Haugen M, Schubert M, Vedeler CA. CDR2L is the Major Yo Antibody Target in paraneoplastic cerebellar degeneration. Ann Neurol. 2019;86(2):316–21. https://doi.org/10.1002/ana.25511. Epub 2019/05/31.
Vialatte de Pémille C, Berzero G, Small M, Psimaras D, Giry M, Daniau M, et al. Transcriptomic immune profiling of ovarian cancers in paraneoplastic cerebellar degeneration associated with anti-yo antibodies. Br J Cancer. 2018;119(1):105–13. https://doi.org/10.1038/s41416-018-0125-7. Epub 2018/06/15.
Rodriguez M, Truh LI, O’Neill BP, Lennon VA. Autoimmune paraneoplastic cerebellar degeneration: ultrastructural localization of antibody-binding sites in Purkinje cells. Neurology. 1988;38(9):1380–6. Epub 1988/09/01. PubMed PMID: 3045692.
Hida C, Tsukamoto T, Awano H, Yamamoto T. Ultrastructural localization of anti-purkinje cell antibody-binding sites in paraneoplastic cerebellar degeneration. Arch Neurol. 1994;51(6):555–8. Epub 1994/06/01. PubMed PMID: 8198465.
Herdlevaer I, Kråkenes T, Schubert M, Vedeler CA. Localization of CDR2L and CDR2 in paraneoplastic cerebellar degeneration. Ann Clin Transl Neurol. 2020;7(11):2231–42. https://doi.org/10.1002/acn3.51212. Epub 2020/10/04.
Nanjappa V, Thomas JK, Marimuthu A, Muthusamy B, Radhakrishnan A, Sharma R et al. Plasma Proteome Database as a resource for proteomics research: 2014 update. Nucleic acids research. 2014;42(Database issue):D959-65. Epub 2013/12/07. https://doi.org/10.1093/nar/gkt1251. PubMed PMID: 24304897; PubMed Central PMCID: PMCPMC3965042.
Chitti SV, Gummadi S, Kang T, Shahi S, Marzan AL, Nedeva C, et al. Vesiclepedia 2024: an extracellular vesicles and extracellular particles repository. Nucleic Acids Res. 2023. https://doi.org/10.1093/nar/gkad1007. PubMed PMID: 37953359. Epub 2023/11/13.
Dörner K, Ruggeri C, Zemp I, Kutay U. Ribosome biogenesis factors-from names to functions. EMBO J. 2023;42(7):e112699. https://doi.org/10.15252/embj.2022112699. Epub 2023/02/11.
Galli J, Greenlee J, F1000Research. Paraneoplastic Diseases of the Central Nervous System. 2020;9. Epub 2020/03/21. https://doi.org/10.12688/f1000research.21309.1. PubMed PMID: 32194943; PubMed Central PMCID: PMCPMC7062203.
Zhang M, Matyunina LV, Walker LD, Chen W, Xiao H, Benigno BB, et al. Evidence for the importance of post-transcriptional regulatory changes in ovarian cancer progression and the contribution of miRNAs. Sci Rep. 2017;7(1):8171. https://doi.org/10.1038/s41598-017-08502-z.
Gry M, Rimini R, Strömberg S, Asplund A, Pontén F, Uhlén M, et al. Correlations between RNA and protein expression profiles in 23 human cell lines. BMC Genomics. 2009;10(1):365. https://doi.org/10.1186/1471-2164-10-365.
Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13(4):227–32. https://doi.org/10.1038/nrg3185.
Liu Y, Beyer A, Aebersold R. On the Dependency of Cellular protein levels on mRNA abundance. Cell. 2016;165(3):535–50. https://doi.org/10.1016/j.cell.2016.03.014.
van Riggelen J, Yetil A, Felsher DW. MYC as a regulator of ribosome biogenesis and protein synthesis. Nat Rev Cancer. 2010;10(4):301–9. https://doi.org/10.1038/nrc2819.
Okano HJ, Park WY, Corradi JP, Darnell RB. The cytoplasmic Purkinje onconeural antigen cdr2 down-regulates c-Myc function: implications for neuronal and tumor cell survival. Genes Dev. 1999;13(16):2087–97. https://doi.org/10.1101/gad.13.16.2087. Epub 1999/08/31.
O’Donovan KJ, Diedler J, Couture GC, Fak JJ, Darnell RB. The onconeural antigen cdr2 is a novel APC/C target that acts in mitosis to regulate c-myc target genes in mammalian tumor cells. PLoS ONE. 2010;5(4):e10045. https://doi.org/10.1371/journal.pone.0010045. Epub 2010/04/13.
Lischnig A, Bergqvist M, Ochiya T, Lässer C. Quantitative Proteomics Identifies Proteins Enriched in Large and Small Extracellular Vesicles. Molecular & cellular proteomics: MCP. 2022;21(9):100273. Epub 2022/08/03. https://doi.org/10.1016/j.mcpro.2022.100273. PubMed PMID: 35918030; PubMed Central PMCID: PMCPMC9486130.
Ochkasova A, Arbuzov G, Malygin A, Graifer D. Two edges in our knowledge on the functions of ribosomal proteins: the revealed contributions of their regions to translation mechanisms and the issues of their Extracellular Transport by Exosomes. Int J Mol Sci. 2023;24(14). https://doi.org/10.3390/ijms241411458. PubMed PMID: 37511213; PubMed Central PMCID: PMCPMC10380927. Epub 2023/07/29.
Dabbah M, Lishner M, Jarchowsky-Dolberg O, Tartakover-Matalon S, Brin YS, Pasmanik-Chor M, et al. Ribosomal proteins as distinct passengers of microvesicles: new semantics in myeloma and mesenchymal stem cells’ communication. Translational Res. 2021;236:117–32. https://doi.org/10.1016/j.trsl.2021.04.002.
Wang W, Nag S, Zhang X, Wang MH, Wang H, Zhou J, et al. Ribosomal proteins and human diseases: pathogenesis, molecular mechanisms, and therapeutic implications. Med Res Rev. 2015;35(2):225–85. https://doi.org/10.1002/med.21327. Epub 2014/08/29.
Nebenfuehr S, Kollmann K, Sexl V. The role of CDK6 in cancer. Int J Cancer. 2020;147(11):2988–95. https://doi.org/10.1002/ijc.33054. Epub 2020/05/15.
Yang X, Xu L, Yang YE, Xiong C, Yu J, Wang Y, et al. Knockdown of ribosomal protein S6 suppresses proliferation, migration, and invasion in epithelial ovarian cancer. J Ovarian Res. 2020;13(1):100. https://doi.org/10.1186/s13048-020-00707-7. Epub 2020/08/31.
Slomnicki LP, Pietrzak M, Vashishta A, Jones J, Lynch N, Elliot S, et al. Requirement of neuronal ribosome synthesis for growth and maintenance of the dendritic tree. J Biol Chem. 2016;291(11):5721–39. https://doi.org/10.1074/jbc.M115.682161. Epub 2016/01/14.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9. https://doi.org/10.1038/nmeth.4197. Epub 2017/03/07.
Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Research. 2015;4:1521. https://doi.org/10.12688/f1000research.7563.2. Epub 2016/03/01.
Love MI, Huber W, Anders S. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8. Epub 2014/12/18. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.
Stekhoven DJ, Bühlmann P. MissForest–non-parametric missing value imputation for mixed-type data. Bioinf (Oxford England). 2012;28(1):112–8. https://doi.org/10.1093/bioinformatics/btr597. Epub 2011/11/01.
Jin L, Bi Y, Hu C, Qu J, Shen S, Wang X, et al. A comparative study of evaluating missing value imputation methods in label-free proteomics. Sci Rep. 2021;11(1):1760. https://doi.org/10.1038/s41598-021-81279-4.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Nucleic Acids Res. 2015;43(7):e47. https://doi.org/10.1093/nar/gkv007. Epub 2015/01/22. limma powers differential expression analyses for RNA-sequencing and microarray studies.
Sergushichev AA. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation. bioRxiv. 2016;060012. https://doi.org/10.1101/060012.
Consortium TGO, Aleksander SA, Balhoff J, Carbon S, Cherry JM, Drabkin HJ, et al. The Gene Ontology knowledgebase in 2023. Genetics. 2023;224(1). https://doi.org/10.1093/genetics/iyad031.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92. https://doi.org/10.1093/nar/gkac963. PubMed PMID: 36300620; PubMed Central PMCID: PMCPMC9825424. Epub 2022/10/28.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Proc Natl Acad Sci USA. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102. Epub 2005/10/04. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.
Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B, WebGestalt. 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic acids research. 2019;47(W1):W199-w205. Epub 2019/05/23. https://doi.org/10.1093/nar/gkz401. PubMed PMID: 31114916; PubMed Central PMCID: PMCPMC6602449.
Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J Proteome Res. 2019;18(2):623–32. https://doi.org/10.1021/acs.jproteome.8b00702. Epub 2018/11/20.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303. Epub 2003/11/05.
Legeay M, Doncheva NT, Morris JH, Jensen LJ. Visualize omics data on networks with Omics Visualizer, a Cytoscape App. F1000Research. 2020;9:157. https://doi.org/10.12688/f1000research.22280.2. Epub 2020/07/15.
Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC systems biology. 2014;8 Suppl 4(Suppl 4):S11. Epub 2014/12/19. https://doi.org/10.1186/1752-0509-8-s4-s11. PubMed PMID: 25521941; PubMed Central PMCID: PMCPMC4290687.
Morris JH, Apeltsin L, Newman AM, Baumbach J, Wittkop T, Su G, et al. clusterMaker: a multi-algorithm clustering plugin for Cytoscape. BMC Bioinformatics. 2011;12:436. https://doi.org/10.1186/1471-2105-12-436. Epub 2011/11/11.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinf (Oxford England). 2009;25(8):1091–3. https://doi.org/10.1093/bioinformatics/btp101. Epub 2009/02/25.
Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinf (Oxford England). 2013;29(5):661–3. https://doi.org/10.1093/bioinformatics/btt019. Epub 2013/01/18.
Keenan AB, Torre D, Lachmann A, Leong AK, Wojciechowicz ML, Utti V, et al. ChEA3: transcription factor enrichment analysis by orthogonal omics integration. Nucleic Acids Res. 2019;47(W1):W212–24. https://doi.org/10.1093/nar/gkz446.
Acknowledgements
We thank Laurence Bindoff for valuable discussion and advice.
Funding
Open access funding provided by University of Bergen. This study was funded by unrestrictive grants from the legacy of Gerda Meyer Nyquist Gulbranson & Gerdt Meyer Nyquist. The founding body had no impact on the design, analysis, and interpretation of data, nor the writing of the manuscript.
Open access funding provided by University of Bergen.
Author information
Authors and Affiliations
Contributions
E.T.S., T.K., C.T. and C.A.V. conceptualized and designed the study. E.B. performed the mass spectrometry analysis. E.T.S., F.D. and Y.G. analysed and interpreted the omics data. I.H. analyzed the mass spectrometry data for interaction partners. E.T.S., Y.G., F.D. and C.A.V. were the major contributors in writing the manuscript. All authors revised and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
12885_2024_12399_MOESM1_ESM.xlsx
Supplementary Material 1. supplementary_table_S1.xlsx. Transcriptome dataset. Expression data provided as counts per million. Genes filtered for mitochondria-encoded genes and lowly expressed genes.
12885_2024_12399_MOESM2_ESM.xlsx
Supplementary Material 2. supplementary_table_S2.xlsx. Proteome dataset. Expression data provided as normalized abundance values. Proteins were filtered for contaminants and lowly expressed proteins.
12885_2024_12399_MOESM3_ESM.xlsx
Supplementary Material 3. supplementary_table_S3.xlsx. Secretome dataset. Expression data provided as normalized abundance values. Proteins were filtered for contaminants and lowly expressed proteins.
12885_2024_12399_MOESM4_ESM.xlsx
Supplementary Material 4. supplementary_table_S4.xlsx. Results of differential expression analysis of knockout CDR2L cells vs wildtype cells.
12885_2024_12399_MOESM5_ESM.xlsx
Supplementary Material 5. supplementary_table_S5.xlsx. Results of differential expression analysis of knockout CDR2 cells vs wildtype cells.
12885_2024_12399_MOESM6_ESM.xlsx
Supplementary Material 6. supplementary_table_S6.xlsx. Results of differential expression analysis of knockout CDR1 cells vs wildtype cells.
12885_2024_12399_MOESM7_ESM.xlsx
Supplementary Material 7. supplementary_table_S7.xlsx. Results of over-representation analysis of differentially expressed genes common across all knockout cell lines.
12885_2024_12399_MOESM8_ESM.xlsx
Supplementary Material 8. supplementary_table_S8.xlsx. Results of gene set enrichment analysis of knockout CDR2L cells vs wildtype cells.
12885_2024_12399_MOESM9_ESM.xlsx
Supplementary Material 9. supplementary_table_S9.xlsx. Results of transcription factor enrichment analysis of knockout CDR2L cells. Analysis was performed using ChEA3 and results are given as mean rank.
12885_2024_12399_MOESM10_ESM.xlsx
Supplementary Material 10. supplementary_table_S10.xlsx. Proteins identified as potential interaction partners of CDR2L using co-immunoprecipitation and mass spectrometry. Listed are the proteins along with the antibodies used for the co-immunoprecipitation.
12885_2024_12399_MOESM11_ESM.xlsx
Supplementary Material 11. supplementary_table_S11.xlsx. Results of enrichment analysis of co-differentially expressed genes from the STRING network.
12885_2024_12399_MOESM12_ESM.png
Supplementary Material 12. supplementary_figure_S1.png. Scatterplots with associated Spearman’s rank correlation of gene expression between datasets. (A) RNA expression in counts per million (CPM; x-axis) and proteome protein expression as log2 transformed abundance values (y-axis). (B) RNA expression in CPM (x-axis) and secretome protein expression as log2 transformed abundance values (y-axis). (C) Proteome protein expression as log2 transformed abundance values (x-axis) and secretome protein expression as log2 transformed abundance values (y-axis).
12885_2024_12399_MOESM13_ESM.png
Supplementary Material 13. supplementary_figure_S2.png. A-B, Scatter plots showing log2 fold change of all genes significantly differentially expressed (FDR < 0.05) in both transcriptome (x-axis) and proteome (y-axis) in knockout CDR1 cells (A) and knockout CDR2 cells (B). The colour represents the four groups of possible combinations of direction of change. Correlation between log2 fold changes in transcriptome and proteome was assessed using Spearman’s rank correlation.
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.
About this article
Cite this article
Solheim, E.T., Gerking, Y., Kråkenes, T. et al. Multi-omics profiling reveals dysregulated ribosome biogenesis and impaired cell proliferation following knockout of CDR2L. BMC Cancer 24, 645 (2024). https://doi.org/10.1186/s12885-024-12399-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12885-024-12399-z