Population Dynamics is a Cancer Driver.

Most tissues are continuously renovated through the division of stem cells and the death of old or damaged cells, which is known as cell turnover rate (CTOR). Despite being in steady state, tissues have different population dynamics and leading to diverse clonality levels. Here, we propose and test that cell population dynamics can be a cancer driver. We employed the evolutionary software esiCancer to show that CTOR, within a range comparable to what is observed in human tissues, can amplify the risk of a mutation due to ancestral selection (ANSEL). In a high CTOR tissue, a mutated ancestral cell is likely to be selected and persist over generations, which leads to a scenario of elevated ANSEL profile, characterized by few niches of large clones, which does not occur in low CTOR. We found that CTOR is significantly associated with the risk of developing cancer, even when correcting for mutation load, indicating that population dynamics per se is a cancer driver. This concept is central to understanding cancer risk and for the design of new therapeutic interventions that minimize the contribution of ANSEL in cancer growth.

Impact of the thyroid hormone T3 and its nuclear receptor TRα1 on colon cancer stem cell phenotypes and response to chemotherapies.

Colorectal cancers (CRCs) are highly heterogeneous and show a hierarchical organization, with cancer stem cells (CSCs) responsible for tumor development, maintenance, and drug resistance. Our previous studies showed the importance of thyroid hormone-dependent signaling on intestinal tumor development and progression through action on stem cells. These results have a translational value, given that the thyroid hormone nuclear receptor TRα1 is upregulated in human CRCs, including in the molecular subtypes associated with CSC features. We used an established spheroid model generated from the human colon adenocarcinoma cell line Caco2 to study the effects of T3 and TRα1 on spheroid formation, growth, and response to conventional chemotherapies. Our results show that T3 treatment and/or increased TRα1 expression in spheroids impaired the response to FOLFIRI and conferred a survival advantage. This was achieved by stimulating drug detoxification pathways and increasing ALDH1A1-expressing cells, including CSCs, within spheroids. These results suggest that clinical evaluation of the thyroid axis and assessing TRα1 levels in CRCs could help to select optimal therapeutic regimens for patients with CRC. Proposed mechanism of action of T3/TRα1 in colon cancer spheroids. In the control condition, TRα1 participates in maintaining homeostatic cell conditions. The presence of T3 in the culture medium activates TRα1 action on target genes, including the drug efflux pumps ABCG2 and ABCB1. In the case of chemotherapy FOLFIRI, the increased expression of ABC transcripts and proteins induced by T3 treatment is responsible for the augmented efflux of 5-FU and Irinotecan from the cancer cells. Taken together, these mechanisms contribute to the decreased efficacy of the chemotherapy and allow cells to escape the treatment. Created with BioRender.com .

SERBP1 interacts with PARP1 and is present in PARylation-dependent protein complexes regulating splicing, cell division, and ribosome biogenesis.

RNA binding proteins (RBPs) containing intrinsically disordered regions (IDRs) are present in diverse molecular complexes where they function as dynamic regulators. Their characteristics promote liquid-liquid phase separation (LLPS) and the formation of membraneless organelles such as stress granules and nucleoli. IDR-RBPs are particularly relevant in the nervous system and their dysfunction is associated with neurodegenerative diseases and brain tumor development. SERBP1 is a unique member of this group, being mostly disordered and lacking canonical RNA-binding domains. Using a proteomics approach followed by functional analysis, we defined SERBP1's interactome. We uncovered novel SERBP1 roles in splicing, cell division, and ribosomal biogenesis and showed its participation in pathological stress granules and Tau aggregates in Alzheimer's disease brains. SERBP1 preferentially interacts with other G-quadruplex (G4) binders, implicated in different stages of gene expression, suggesting that G4 binding is a critical component of SERBP1 function in different settings. Similarly, we identified important associations between SERBP1 and PARP1/polyADP-ribosylation (PARylation). SERBP1 interacts with PARP1 and its associated factors and influences PARylation. Moreover, protein complexes in which SERBP1 participates contain mostly PARylated proteins and PAR binders. Based on these results, we propose a feedback regulatory model in which SERBP1 influences PARP1 function and PARylation, while PARylation modulates SERBP1 functions and participation in regulatory complexes.

Retro-miRs: novel and functional miRNAs originating from mRNA retrotransposition.

BACKGROUND: Reverse-transcribed gene copies (retrocopies) have emerged as major sources of evolutionary novelty. MicroRNAs (miRNAs) are small and highly conserved RNA molecules that serve as key post-transcriptional regulators of gene expression. The origin and subsequent evolution of miRNAs have been addressed but not fully elucidated. RESULTS: In this study, we performed a comprehensive investigation of miRNA origination through retroduplicated mRNA sequences (retro-miRs). We identified 17 retro-miRs that emerged from the mRNA retrocopies. Four of these retro-miRs had de novo origins within retrocopied sequences, while 13 retro-miRNAs were located within exon regions and duplicated along with their host mRNAs. We found that retro-miRs were primate-specific, including five retro-miRs conserved among all primates and two human-specific retro-miRs. All retro-miRs were expressed, with predicted and experimentally validated target genes except miR-10527. Notably, the target genes of retro-miRs are involved in key biological processes such as metabolic processes, cell signaling, and regulation of neurotransmitters in the central nervous system. Additionally, we found that these retro-miRs play a potential oncogenic role in cancer by targeting key cancer genes and are overexpressed in several cancer types, including liver hepatocellular carcinoma and stomach adenocarcinoma. CONCLUSIONS: Our findings demonstrated that mRNA retrotransposition is a key mechanism for the generation of novel miRNAs (retro-miRs) in primates. These retro-miRs are expressed, conserved, have target genes with important cellular functions, and play important roles in cancer.

Acute Myeloid Leukemia Expresses a Specific Group of Olfactory Receptors.

Acute myeloid leukemia (AML) is the most common form of acute leukemia in adults, with a 5-year overall survival rate of approximately 30%. Despite recent advances in therapeutic options, relapse remains the leading cause of death and poor survival outcomes. New drugs benefit specific small subgroups of patients with actionable therapeutic targets. Thus, finding new targets with greater applicability should be pursued. Olfactory receptors (ORs) are seven transmembrane G-protein coupled receptors preferentially expressed in sensory neurons with a critical role in recognizing odorant molecules. Recent studies have revealed ectopic expression and putative function of ORs in nonolfactory tissues and pathologies, including AML. Here, we investigated OR expression in 151 AML samples, 6400 samples of 15 other cancer types, and 11,200 samples of 51 types of healthy tissues. First, we identified 19 ORs with a distinct and major expression pattern in AML, which were experimentally validated by RT-PCR in an independent set of 13 AML samples, 13 healthy donors, and 8 leukemia cell lines. We also identified an OR signature with prognostic potential for AML patients. Finally, we found cancer-related genes coexpressed with the ORs in the AML samples. In summary, we conducted an extensive study to identify ORs that can be used as novel biomarkers for the diagnosis of AML and as potential drug targets.

The paralogues MAGOH and MAGOHB are oncogenic factors in high-grade gliomas and safeguard the splicing of cell division and cell cycle genes.

The exon junction complex (EJC) plays key roles throughout the lifespan of RNA and is particularly relevant in the nervous system. We investigated the roles of two EJC members, the paralogs MAGOH and MAGOHB, with respect to brain tumour development. High MAGOH/MAGOHB expression was observed in 14 tumour types; glioblastoma (GBM) showed the greatest difference compared to normal tissue. Increased MAGOH/MAGOHB expression was associated with poor prognosis in glioma patients, while knockdown of MAGOH/MAGOHB affected different cancer phenotypes. Reduced MAGOH/MAGOHB expression in GBM cells caused alterations in the splicing profile, including re-splicing and skipping of multiple exons. The binding profiles of EJC proteins indicated that exons affected by MAGOH/MAGOHB knockdown accumulated fewer complexes on average, providing a possible explanation for their sensitivity to MAGOH/MAGOHB knockdown. Transcripts (genes) showing alterations in the splicing profile are mainly implicated in cell division, cell cycle, splicing, and translation. We propose that high MAGOH/MAGOHB levels are required to safeguard the splicing of genes in high demand in scenarios requiring increased cell proliferation (brain development and GBM growth), ensuring efficient cell division, cell cycle regulation, and gene expression (splicing and translation). Since differentiated neuronal cells do not require increased MAGOH/MAGOHB expression, targeting these paralogs is a potential option for treating GBM.

ELF4 is a critical component of a miRNA-transcription factor network and is a bridge regulator of glioblastoma receptor signaling and lipid dynamics.

BACKGROUND: The loss of neurogenic tumor suppressor microRNAs miR-124, miR-128, and miR-137 is associated with glioblastoma's undifferentiated state. Most of their impact comes via the repression of a network of oncogenic transcription factors. We conducted a high-throughput functional siRNA screen in glioblastoma cells and identify E74 like ETS transcription factor 4 (ELF4) as the leading contributor to oncogenic phenotypes. METHODS: In vitro and in vivo assays were used to assess ELF4 impact on cancer phenotypes. We characterized ELF4's mechanism of action via genomic and lipidomic analyses. A MAPK reporter assay verified ELF4's impact on MAPK signaling, and qRT-PCR and western blotting were used to corroborate ELF4 regulatory role on most relevant target genes. RESULTS: ELF4 knockdown resulted in significant proliferation delay and apoptosis in GBM cells and long-term growth delay and morphological changes in glioma stem cells (GSCs). Transcriptomic analyses revealed that ELF4 controls two interlinked pathways: 1) Receptor tyrosine kinase signaling and 2) Lipid dynamics. ELF4 modulation directly affected receptor tyrosine kinase (RTK) signaling, as mitogen-activated protein kinase (MAPK) activity was dependent upon ELF4 levels. Furthermore, shotgun lipidomics revealed that ELF4 depletion disrupted several phospholipid classes, highlighting ELF4's importance in lipid homeostasis. CONCLUSIONS: We found that ELF4 is critical for the GBM cell identity by controlling genes of two dependent pathways: RTK signaling (SRC, PTK2B, and TNK2) and lipid dynamics (LRP1, APOE, ABCA7, PLA2G6, and PITPNM2). Our data suggest that targeting these two pathways simultaneously may be therapeutically beneficial to GBM patients.

Tissue-resident glial cells associate with tumoral vasculature and promote cancer progression.

Cancer cells are embedded within the tissue and interact dynamically with its components during cancer progression. Understanding the contribution of cellular components within the tumor microenvironment is crucial for the success of therapeutic applications. Here, we reveal the presence of perivascular GFAP+/Plp1+ cells within the tumor microenvironment. Using in vivo inducible Cre/loxP mediated systems, we demonstrated that these cells derive from tissue-resident Schwann cells. Genetic ablation of endogenous Schwann cells slowed down tumor growth and angiogenesis. Schwann cell-specific depletion also induced a boost in the immune surveillance by increasing tumor-infiltrating anti-tumor lymphocytes, while reducing immune-suppressor cells. In humans, a retrospective in silico analysis of tumor biopsies revealed that increased expression of Schwann cell-related genes within melanoma was associated with improved survival. Collectively, our study suggests that Schwann cells regulate tumor progression, indicating that manipulation of Schwann cells may provide a valuable tool to improve cancer patients' outcomes.

Interplay between Tumor Mutational Burden and Mutational Profile and Its Effect on Overall Survival: A Pilot Study of Metastatic Patients Treated with Immune Checkpoint Inhibitors.

Purpose: Solid tumors harboring tumor mutational burden (TMB) ≥10 mutations per megabase (mut/Mb) received agnostic approval for pembrolizumab. This work aims to analyze the somatic mutational profile's influence on the outcomes of patients with TMB-high tumors treated with immune checkpoint inhibitors (ICIs). Methods: This post-hoc analysis evaluated clinical and molecular features of patients with solid tumors treated with ICIs that could be either monoclonal antibody directed against programmed cell death protein-1 or monoclonal antibody directed against programmed cell death ligand 1 (anti-PD-1/anti-PD-L1), monoclonal antibody directed against cytotoxic T lymphocyte-associated antigen (anti-CTLA-4) or a combined treatment regimen including one anti-PD-1/anti-PD-L1 and one anti-CTLA-4 (ICIs combination). We performed OS analysis for TMB thresholds of ≥10, ≥20, and <10 mut/Mb. We assessed OS according to the mutational profile for a TMB ≥ 10 mut/Mb cutoff. For genes correlated with OS at the univariate assessment, we conducted a Cox multivariate analysis adjusted by median TMB, sex, age, microsatellite instability (MSI), and histology. Results: A total of 1661 patients were investigated; 488 with a TMB ≥10 mut/Mb (29.4%). The median OS was 42 months for TMB ≥10 or 20 mut/Mb, and 15 months for TMB <10 mut/Mb (p < 0.005). Among TMB ≥10 mut/Mb patients, mutations in E2F3 or STK11 correlated with worse OS, and mutations in NTRK3, PTPRD, RNF43, TENT5C, TET1, or ZFHX3 with better OS. These associations were confirmed with univariate and multivariate analyses (p < 0.05). Melanoma histology and TMB above the median endowed patients with better OS (p < 0.05), while MSI status, age, and gender did not have a statistically significant effect on OS. Conclusion: Combining TMB and mutation profiles in key cancer genes can better qualify patients for ICI treatment and predict their OS.

Regulation of the THRA gene, encoding the thyroid hormone nuclear receptor TRα1, in intestinal lesions.

The THRA gene, encoding the thyroid hormone nuclear receptor TRα1, is expressed in an increasing gradient at the bottom of intestinal crypts, overlapping with high Wnt and Notch activities. Importantly, THRA is upregulated in colorectal cancers, particularly in the high-Wnt molecular subtype. The basis of this specific and/or altered expression pattern has remained unknown. To define the mechanisms controlling THRA transcription and TRα1 expression, we used multiple in vitro and ex vivo approaches. Promoter analysis demonstrated that transcription factors important for crypt homeostasis and altered in colorectal cancers, such as transcription factor 7-like 2 (TCF7L2; Wnt pathway), recombining binding protein suppressor of hairless (RBPJ; Notch pathway), and homeobox protein CDX2 (epithelial cell identity), modulate THRA activity. Specifically, although TCF7L2 and CDX2 stimulated THRA, RBPJ induced its repression. In-depth analysis of the Wnt-dependent increase showed direct regulation of the THRA promoter in cells and of TRα1 expression in murine enteroids. Given our previous results on the control of the Wnt pathway by TRα1, our new results unveil a complex regulatory loop and synergy between these endocrine and epithelial-cell-intrinsic signals. Our work describes, for the first time, the regulation of the THRA gene in specific cell and tumor contexts.

Chemogenetic modulation of sensory neurons reveals their regulating role in melanoma progression.

Sensory neurons have recently emerged as components of the tumor microenvironment. Nevertheless, whether sensory neuronal activity is important for tumor progression remains unknown. Here we used Designer Receptors Exclusively Activated by a Designer Drug (DREADD) technology to inhibit or activate sensory neurons' firing within the melanoma tumor. Melanoma growth and angiogenesis were accelerated following inhibition of sensory neurons' activity and were reduced following overstimulation of these neurons. Sensory neuron-specific overactivation also induced a boost in the immune surveillance by increasing tumor-infiltrating anti-tumor lymphocytes, while reducing immune-suppressor cells. In humans, a retrospective in silico analysis of melanoma biopsies revealed that increased expression of sensory neurons-related genes within melanoma was associated with improved survival. These findings suggest that sensory innervations regulate melanoma progression, indicating that manipulation of sensory neurons' activity may provide a valuable tool to improve melanoma patients' outcomes.

Reboot: a straightforward approach to identify genes and splicing isoforms associated with cancer patient prognosis.

Nowadays, the massive amount of data generated by modern sequencing technologies provides an unprecedented opportunity to find genes associated with cancer patient prognosis, connecting basic and translational research. However, treating high dimensionality of gene expression data and integrating it with clinical variables are major challenges to perform these analyses. Here, we present Reboot, an integrative approach to find and validate genes and transcripts (splicing isoforms) associated with cancer patient prognosis from high dimensional expression datasets. Reboot innovates by using a multivariate strategy with penalized Cox regression (LASSO method) combined with a bootstrap approach, in addition to statistical tests and plots to support the findings. Applying Reboot on data from 154 glioblastoma patients, we identified a three-gene signature (IKBIP, OSMR, PODNL1) whose increased derived risk score was significantly associated with worse patients' prognosis. Similarly, Reboot was able to find a seven-splicing isoforms signature related to worse overall survival in 177 pancreatic adenocarcinoma patients with elevated risk scores after uni- and multivariate analyses. In summary, Reboot is an efficient, intuitive and straightforward way of finding genes or splicing isoforms signatures relevant to patient prognosis, which can democratize this kind of analysis and shed light on still under-investigated cancer-related genes and splicing isoforms.

Deciphering the Role of Intestinal Crypt Cell Populations in Resistance to Chemotherapy.

Intestinal crypts are composed of heterogeneous and highly plastic cell populations. Lgr5high-stem cells (SC) are responsible for homeostatic renewal, but other cells can revert to an SC-like phenotype to maintain epithelial integrity. Despite their distinct roles in orchestrating homeostasis, both populations have been designated as the putative "cell-of-origin" of colorectal cancer. However, their respective involvement in the emergence of drug-resistant cancer SCs (CSC), responsible for tumor relapse and associated with poor outcome of colorectal cancer, remains elusive. In this context, the intestinal SC/progenitor-marker Musashi1 (MSI1) is interesting as it plays important functions in intestinal homeostasis and is frequently overexpressed in human colorectal cancer. Therefore, our aims were: (i) to study the impact of chemotherapy on Lgr5-expressing and MSI1-expressing cell populations, (ii) to explore the effect of increased MSI1 levels in response to treatment, and (iii) to evaluate the relevance in human colorectal cancer. Engineered mouse models treated with the therapeutic agent 5-fluorouracil showed that upon increased MSI1 levels, Lgr5high SCs remain sensitive while Lgr5low progenitors reprogram to a drug-resistant phenotype. This resulted in the expansion of an MSI1-expressing cell subpopulation with improved resistance to DNA damage and increased detoxification, typical properties of dormant-CSCs that can reactivate after chemotherapy. Analysis in patients with colorectal cancer revealed a correlation between MSI1 levels and tumor grading, CSC phenotype, and chemoresistance. Altogether, these results shed new light on the biology and plasticity of normal crypt and cancer cell populations and also open new perspectives to target MSI1 to improve chemotherapy outcome. SIGNIFICANCE: This study unveils paradoxical roles for MSI1, underlining its importance in facilitating intestinal regeneration upon injury but also unraveling its new function in drug-resistant colorectal cancer stem cells.

Synergism of Proneurogenic miRNAs Provides a More Effective Strategy to Target Glioma Stem Cells.

Tumor suppressor microRNAs (miRNAs) have been explored as agents to target cancer stem cells. Most strategies use a single miRNA mimic and present many disadvantages, such as the amount of reagent required and the diluted effect on target genes. miRNAs work in a cooperative fashion to regulate distinct biological processes and pathways. Therefore, we propose that miRNA combinations could provide more efficient ways to target cancer stem cells. We have previously shown that miR-124, miR-128, and miR-137 function synergistically to regulate neurogenesis. We used a combination of these three miRNAs to treat glioma stem cells and showed that this treatment was much more effective than single miRNAs in disrupting cell proliferation and survival and promoting differentiation and response to radiation. Transcriptomic analyses indicated that transcription regulation, angiogenesis, metabolism, and neuronal differentiation are among the main biological processes affected by transfection of this miRNA combination. In conclusion, we demonstrated the value of using combinations of neurogenic miRNAs to disrupt cancer phenotypes and glioma stem cell growth. The synergistic effect of these three miRNA amplified the repression of oncogenic factors and the effect on cancer relevant pathways. Future therapeutic approaches would benefit from utilizing miRNA combinations, especially when targeting cancer-initiating cell populations.

The RNA-Binding Protein Musashi1 Regulates a Network of Cell Cycle Genes in Group 4 Medulloblastoma.

Medulloblastoma is the most common malignant brain tumor in children. Treatment with surgery, irradiation, and chemotherapy has improved survival in recent years, but patients are frequently left with devastating neurocognitive and other sequelae. Patients in molecular subgroups 3 and 4 still experience a high mortality rate. To identify new pathways contributing to medulloblastoma development and create new routes for therapy, we have been studying oncogenic RNA-binding proteins. We defined Musashi1 (Msi1) as one of the main drivers of medulloblastoma development. The high expression of Msi1 is prevalent in Group 4 and correlates with poor prognosis while its knockdown disrupted cancer-relevant phenotypes. Genomic analyses (RNA-seq and RIP-seq) indicated that cell cycle and division are the main biological categories regulated by Msi1 in Group 4 medulloblastoma. The most prominent Msi1 targets include CDK2, CDK6, CCND1, CDKN2A, and CCNA1. The inhibition of Msi1 with luteolin affected the growth of CHLA-01 and CHLA-01R Group 4 medulloblastoma cells and a synergistic effect was observed when luteolin and the mitosis inhibitor, vincristine, were combined. These findings indicate that a combined therapeutic strategy (Msi1 + cell cycle/division inhibitors) could work as an alternative to treat Group 4 medulloblastoma.

The RNA-binding protein SERBP1 functions as a novel oncogenic factor in glioblastoma by bridging cancer metabolism and epigenetic regulation.

BACKGROUND: RNA-binding proteins (RBPs) function as master regulators of gene expression. Alterations in RBP expression and function are often observed in cancer and influence critical pathways implicated in tumor initiation and growth. Identification and characterization of oncogenic RBPs and their regulatory networks provide new opportunities for targeted therapy. RESULTS: We identify the RNA-binding protein SERBP1 as a novel regulator of glioblastoma (GBM) development. High SERBP1 expression is prevalent in GBMs and correlates with poor patient survival and poor response to chemo- and radiotherapy. SERBP1 knockdown causes delay in tumor growth and impacts cancer-relevant phenotypes in GBM and glioma stem cell lines. RNAcompete identifies a GC-rich region as SERBP1-binding motif; subsequent genomic and functional analyses establish SERBP1 regulation role in metabolic routes preferentially used by cancer cells. An important consequence of these functions is SERBP1 impact on methionine production. SERBP1 knockdown decreases methionine levels causing a subsequent reduction in histone methylation as shown for H3K27me3 and upregulation of genes associated with neurogenesis, neuronal differentiation, and function. Further analysis demonstrates that several of these genes are downregulated in GBM, potentially through epigenetic silencing as indicated by the presence of H3K27me3 sites. CONCLUSIONS: SERBP1 is the first example of an RNA-binding protein functioning as a central regulator of cancer metabolism and indirect modulator of epigenetic regulation in GBM. By bridging these two processes, SERBP1 enhances glioma stem cell phenotypes and contributes to GBM poorly differentiated state.

Proneural and mesenchymal glioma stem cells display major differences in splicing and lncRNA profiles.

Therapy resistance and recurrence in high-grade gliomas are driven by their populations of glioma stem cells (GSCs). Thus, detailed molecular characterization of GSCs is needed to develop more effective therapies. We conducted a study to identify differences in the splicing profile and expression of long non-coding RNAs in proneural and mesenchymal GSC cell lines. Genes related to cell cycle, DNA repair, cilium assembly, and splicing showed the most differences between GSC subgroups. We also identified genes distinctly associated with survival among patients of mesenchymal or proneural subgroups. We determined that multiple long non-coding RNAs with increased expression in mesenchymal GSCs are associated with poor survival of glioblastoma patients. In summary, our study established critical differences between proneural and mesenchymal GSCs in splicing profiles and expression of long non-coding RNA. These splicing isoforms and lncRNA signatures may contribute to the uniqueness of GSC subgroups, thus contributing to cancer phenotypes and explaining differences in therapeutic responses.

Antagonism between the RNA-binding protein Musashi1 and miR-137 and its potential impact on neurogenesis and glioblastoma development.

RNA-binding proteins (RBPs) and miRNAs are critical gene expression regulators that interact with one another in cooperative and antagonistic fashions. We identified Musashi1 (Msi1) and miR-137 as regulators of a molecular switch between self-renewal and differentiation. Msi1 and miR-137 have opposite expression patterns and functions, and Msi1 is repressed by miR-137. Msi1 is a stem-cell protein implicated in self-renewal while miR-137 functions as a proneuronal differentiation miRNA. In gliomas, miR-137 functions as a tumor suppressor while Msi1 is a prooncogenic factor. We suggest that the balance between Msi1 and miR-137 is a key determinant in cell fate decisions and disruption of this balance could contribute to neurodegenerative diseases and glioma development. Genomic analyses revealed that Msi1 and miR-137 share 141 target genes associated with differentiation, development, and morphogenesis. Initial results pointed out that these two regulators have an opposite impact on the expression of their target genes. Therefore, we propose an antagonistic model in which this network of shared targets could be either repressed by miR-137 or activated by Msi1, leading to different outcomes (self-renewal, proliferation, tumorigenesis).

A Methodology for the Development of RESTful Semantic Web Services for Gene Expression Analysis.

Gene expression studies are generally performed through multi-step analysis processes, which require the integrated use of a number of analysis tools. In order to facilitate tool/data integration, an increasing number of analysis tools have been developed as or adapted to semantic web services. In recent years, some approaches have been defined for the development and semantic annotation of web services created from legacy software tools, but these approaches still present many limitations. In addition, to the best of our knowledge, no suitable approach has been defined for the functional genomics domain. Therefore, this paper aims at defining an integrated methodology for the implementation of RESTful semantic web services created from gene expression analysis tools and the semantic annotation of such services. We have applied our methodology to the development of a number of services to support the analysis of different types of gene expression data, including microarray and RNASeq. All developed services are publicly available in the Gene Expression Analysis Services (GEAS) Repository at http://dcm.ffclrp.usp.br/lssb/geas. Additionally, we have used a number of the developed services to create different integrated analysis scenarios to reproduce parts of two gene expression studies documented in the literature. The first study involves the analysis of one-color microarray data obtained from multiple sclerosis patients and healthy donors. The second study comprises the analysis of RNA-Seq data obtained from melanoma cells to investigate the role of the remodeller BRG1 in the proliferation and morphology of these cells. Our methodology provides concrete guidelines and technical details in order to facilitate the systematic development of semantic web services. Moreover, it encourages the development and reuse of these services for the creation of semantically integrated solutions for gene expression analysis.