RAGE engagement by SARS-CoV-2 enables monocyte infection and underlies COVID-19 severity.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has fueled the COVID-19 pandemic with its enduring medical and socioeconomic challenges because of subsequent waves and long-term consequences of great concern. Here, we chart the molecular basis of COVID-19 pathogenesis by analyzing patients' immune responses at single-cell resolution across disease course and severity. This approach confirms cell subpopulation-specific dysregulation in COVID-19 across disease course and severity and identifies a severity-associated activation of the receptor for advanced glycation endproducts (RAGE) pathway in monocytes. In vitro THP1-based experiments indicate that monocytes bind the SARS-CoV-2 S1-receptor binding domain (RBD) via RAGE, pointing to RAGE-Spike interaction enabling monocyte infection. Thus, our results demonstrate that RAGE is a functional receptor of SARS-CoV-2 contributing to COVID-19 severity.

A perfusion-based three-dimensional cell culture system to model alveolar rhabdomyosarcoma pathological features.

Although a rare disease, rhabdomyosarcoma (RMS) is one of the most common cancers in children the more aggressive and metastatic subtype is the alveolar RMS (ARMS). Survival outcomes with metastatic disease remain dismal and the need for new models that recapitulate key pathological features, including cell-extracellular matrix (ECM) interactions, is warranted. Here, we report an organotypic model that captures cellular and molecular determinants of invasive ARMS. We cultured the ARMS cell line RH30 on a collagen sponge in a perfusion-based bioreactor (U-CUP), obtaining after 7 days a 3D construct with homogeneous cell distribution. Compared to static culture, perfusion flow induced higher cell proliferation rates (20% vs. 5%), enhanced secretion of active MMP-2, and upregulation of the Rho pathway, associated with cancer cell dissemination. Consistently, the ECM genes LAMA1 and LAMA2, the antiapoptotic gene HSP90, identified in patient databases as hallmarks of invasive ARMS, were higher under perfusion flow at mRNA and protein level. Our advanced ARMS organotypic model mimics (1) the interactions cells-ECM, (2) the cell growth maintenance, and (3) the expression of proteins that characterize tumor expansion and aggressiveness. In the future, the perfusion-based model could be used with primary patient-derived cell subtypes to create a personalized ARMS chemotherapy screening system.

3D geometry orchestrates the transcriptional landscape of metastatic neuroblastoma cells in a multicellular in vitro bone model.

A key challenge for the discovery of novel molecular targets and therapeutics against pediatric bone metastatic disease is the lack of bona fide in vitro cell models. Here, we show that a beta-tricalcium phosphate (ß-TCP) multicellular 3D in vitro bone microtissue model reconstitutes key phenotypic and transcriptional patterns of native metastatic tumor cells while promoting their stemness and proinvasive features. Comparing planar with interconnected channeled scaffolds, we identified geometry as a dominant orchestrator of proangiogenic traits in neuroblastoma tumor cells. On the other hand, the ß-TCP-determined gene signature was DNA replication related. Jointly, the geometry and chemical impact of ß-TCP revealed a prometastatic landscape of the engineered tumor microenvironment. The proposed 3D multicellular in vitro model of pediatric bone metastatic disease may advance further analysis of the molecular, genetic and metabolic bases of the disease and allow more efficient preclinical target validations.

Ambra1 deficiency impairs mitophagy in skeletal muscle.

BACKGROUND: Maintaining healthy mitochondria is mandatory for muscle viability and function. An essential surveillance mechanism targeting defective and harmful mitochondria to degradation is the selective form of autophagy called mitophagy. Ambra1 is a multifaceted protein with well-known autophagic and mitophagic functions. However, the study of its role in adult tissues has been extremely limited due to the embryonic lethality caused by full-body Ambra1 deficiency. METHODS: To establish the role of Ambra1 as a positive regulator of mitophagy, we exploited in vivo overexpression of a mitochondria-targeted form of Ambra1 in skeletal muscle. To dissect the consequence of Ambra1 inactivation in skeletal muscle, we generated muscle-specific Ambra1 knockout (Ambra1fl/fl :Mlc1f-Cre) mice. Mitochondria-enriched fractions were obtained from muscles of fed and starved animals to investigate the dynamics of the mitophagic flux. RESULTS: Our data show that Ambra1 has a critical role in the mitophagic flux of adult murine skeletal muscle and that its genetic inactivation leads to mitochondria alterations and myofibre remodelling. Ambra1 overexpression in wild-type muscles is sufficient to enhance mitochondria clearance through the autophagy-lysosome system. Consistently with this, Ambra1-deficient muscles display an abnormal accumulation of the mitochondrial marker TOMM20 by +76% (n = 6-7; P < 0.05), a higher presence of myofibres with swollen mitochondria by +173% (n = 4; P < 0.05), and an alteration in the maintenance of the mitochondrial membrane potential and a 34% reduction in the mitochondrial respiratory complex I activity (n = 4; P < 0.05). Lack of Ambra1 in skeletal muscle leads to impaired mitophagic flux, without affecting the bulk autophagic process. This is due to a significantly decreased recruitment of DRP1 (n = 6-7 mice; P < 0.01) and Parkin (n = 6-7 mice; P < 0.05) to the mitochondrial compartment, when compared with controls. Ambra1-deficient muscles also show a marked dysregulation of the endolysosome compartment, as the incidence of myofibres with lysosomal accumulation is 20 times higher than wild-type muscles (n = 4; P < 0.05). Histologically, Ambra1-deficient muscles of both 3- and 6-month-old animals display a significant decrease of myofibre cross-sectional area and a 52% reduction in oxidative fibres (n = 6-7; P < 0.05), thus highlighting a role for Ambra1 in the proper structure and activity of skeletal muscle. CONCLUSIONS: Our study indicates that Ambra1 is critical for skeletal muscle mitophagy and for the proper maintenance of functional mitochondria.

The long non-coding RNA CDK6-AS1 overexpression impacts on acute myeloid leukemia differentiation and mitochondrial dynamics.

Patients with acute myeloid leukemia (AML) carrying high-risk genetic lesions or high residual disease levels after therapy are particularly exposed to the risk of relapse. Here, we identified the long non-coding RNA CDK6-AS1 able to cluster an AML subgroup with peculiar gene signatures linked to hematopoietic cell differentiation and mitochondrial dynamics. CDK6-AS1 silencing triggered hematopoietic commitment in healthy CD34+ cells, whereas in AML cells the pathological undifferentiated state was rescued. This latter phenomenon derived from RUNX1 transcriptional control, responsible for the stemness of hematopoietic precursors and for the block of differentiation in AML. By CDK6-AS1 silencing in vitro, AML mitochondrial mass decreased with augmented pharmacological sensitivity to mitochondria-targeting drugs. In vivo, the combination of tigecycline and cytarabine reduced leukemia progression in the AML-PDX model with high CDK6-AS1 levels, supporting the concept of a mitochondrial vulnerability. Together, these findings uncover CDK6-AS1 as crucial in myeloid differentiation and mitochondrial mass regulation.

Integrated CGH/WES Analyses Advance Understanding of Aggressive Neuroblastoma Evolution: A Case Study.

Neuroblastoma (NB) is the most common extra-cranial malignancy in preschool children. To portray the genetic landscape of an overly aggressive NB leading to a rapid clinical progression of the disease, tumor DNA collected pre- and post-treatment has been analyzed. Array comparative genomic hybridization (aCGH), whole-exome sequencing (WES), and pharmacogenetics approaches, respectively, have identified relevant copy number alterations (CNAs), single nucleotide variants (SNVs), and polymorphisms (SNPs) that were then combined into an integrated analysis. Spontaneously formed 3D tumoroids obtained from the recurrent mass have also been characterized. The results prove the power of combining CNAs, SNVs, and SNPs analyses to assess clonal evolution during the disease progression by evidencing multiple clones at disease onset and dynamic genomic alterations during therapy administration. The proposed molecular and cytogenetic integrated analysis empowers the disease follow-up and the prediction of tumor recurrence.

MCM2 and Carbonic Anhydrase 9 Are Novel Potential Targets for Neuroblastoma Pharmacological Treatment.

To overcome the lack of effective pharmacological treatments for high-risk neuroblastoma (HR-NB), the development of novel in vitro and in vivo models that better recapitulate the disease is required. Here, we used an in vitro multiclonal cell model encompassing NB cell differentiation stages, to identify potential novel pharmacological targets. This model allowed us to identify, by low-density RT-PCR arrays, two gene sets, one over-expressed during NB cell differentiation, and the other up-regulated in more malignant cells. Challenging two HR-NB gene expression datasets, we found that these two gene sets are related to high and low survival, respectively. Using mouse NB cisplatin-treated xenografts, we identified two genes within the list associated to the malignant stage (MCM2 and carbonic anhydrase 9), whose expression is positively correlated with tumor growth. Thus, we tested their pharmacological targeting as potential therapeutic strategy. We measured mice survival and tumor growth rate after xenografts of human NB treated with cisplatin in the presence of MCM2/carbonic anhydrase 9 inhibitors (ciprofloxacin and acetazolamide). MCM2 or carbonic anhydrase 9 inhibition significantly increased cisplatin activity, supporting their possible testing for NB therapy.

A BAG's life: Every connection matters in cancer.

The members of the BCL-2 associated athanogene (BAG) family participate in the regulation of a variety of interrelated physiological processes, such as autophagy, apoptosis, and protein homeostasis. Under normal circumstances, the six BAG members described in mammals (BAG1-6) principally assist the 70 kDa heat-shock protein (HSP70) in protein folding; however, their role as oncogenes is becoming increasingly evident. Deregulation of the BAG multigene family has been associated with cell transformation, tumor recurrence, and drug resistance. In addition to BAG overexpression, BAG members are also involved in many oncogenic protein-protein interactions (PPIs). As such, either the inhibition of overloading BAGs or of specific BAG-client protein interactions could have paramount therapeutic value. In this review, we will examine the role of each BAG family member in different malignancies, focusing on their modular structure, which enables interaction with a variety of proteins to exert their pro-tumorigenic role. Lastly, critical remarks on the unmet needs for proposing effective BAG inhibitors will be pointed out.

Next-generation sequencing of PTEN mutations for monitoring minimal residual disease in T-cell acute lymphoblastic leukemia.

Minimal residual disease (MRD) analysis has become a powerful indicator to refine therapy in acute lymphoblastic leukemia (ALL). Here, we present an MRD detection based on the next-generation sequencing of PTEN exon 7 mutations (NGS-PTEN) in 30 pediatric T-cell ALL patients. By comparing the NGS-PTEN results with current quantitative PCR of antigen receptor gene rearrangements (qPCR-Ig/TR), an overall concordance of 80% was found between the two methods. However, the NGS-PTEN qualified a lower number of high-risk patients than qPCR-Ig/TR. These findings suggest that NGS-PTEN is a promising tool that could potentially be used to support current MRD methodologies for T-ALL patients.

Correction to: LIN28B increases neural crest cell migration and leads to transformation of trunk sympathoadrenal precursors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

LIN28B increases neural crest cell migration and leads to transformation of trunk sympathoadrenal precursors.

The RNA-binding protein LIN28B regulates developmental timing and determines stem cell identity by suppressing the let-7 family of microRNAs. Postembryonic reactivation of LIN28B impairs cell commitment to differentiation, prompting their transformation. In this study, we assessed the extent to which ectopic lin28b expression modulates the physiological behavior of neural crest cells (NCC) and governs their transformation in the trunk region of developing embryos. We provide evidence that the overexpression of lin28b inhibits sympathoadrenal cell differentiation and accelerates NCC migration in two vertebrate models, Xenopus leavis and Danio rerio. Our results highlight the relevance of ITGA5 and ITGA6 in the LIN28B-dependent regulation of the invasive motility of tumor cells. The results also establish that LIN28B overexpression supports neuroblastoma onset and the metastatic potential of malignant cells through let-7a-dependent and let-7a-independent mechanisms.

HIF-1α/Wnt signaling-dependent control of gene transcription regulates neuronal differentiation of glioblastoma stem cells.

HIF-1α has been suggested to interplay with Wnt signaling components in order to activate a neuronal differentiation process in both normal brain and glioblastoma (GBM). Based on these data, we explored the molecular mechanisms underlying the observed capability of GBM cells to acquire a neuronal phenotype upon Wnt signaling stimulation and how the microenvironment, particularly hypoxia, modulates this process. Methods: here, the employment of ChIP-seq techniques together with co-immunoprecipitation approaches allowed to reconstruct the molecular interactions responsible for activating specific pro-differentiating transcriptional programs in GBM cells. Moreover, gene silencing/over-expression approaches coupled with the functional analysis of cell phenotype were applied to confirm ChIP-driven hypotheses. Finally, we combined the use of publicly available gene expression datasets with protein expression data by immunohistochemistry to test the clinical relevance of obtained results. Results: our data clearly suggest that HIF-1α is recruited by the ß-catenin/TCF1 complex to foster neuronal differentiation gene transcription in hypoxic GBM cells. Conversely, at higher oxygen levels, the increased expression of TCF4 exerts a transcriptional inhibitory function on the same genomic regions, thus counteracting differentiation. Moreover, we demonstrate the existence of a positive correlation between the expression levels of HIF-1α, TCF1 and neuronal phenotype in GBM tumors, accompanied by the over-expression of several Wnt signaling components, finally affecting patient prognosis. Conclusion: we unveiled a peculiar mechanism by which TCF1 and HIF-1α can induce a reminiscent neuronal differentiation of hypoxic GBM cells, which is hampered, in normoxia, by high levels of TCF4, thus not only de facto controlling the balance between differentiation and stemness, but also impacting on intra-tumoral heterogeneity and eventually patient outcome.

Epigenetic heterogeneity affects the risk of relapse in children with t(8;21)RUNX1-RUNX1T1-rearranged AML.

The somatic translocation t(8;21)(q22;q22)/RUNX1-RUNX1T1 is one of the most frequent rearrangements found in children with standard-risk acute myeloid leukemia (AML). Despite the favorable prognostic role of this aberration, we recently observed a higher than expected frequency of relapse. Here, we employed an integrated high-throughput approach aimed at identifying new biological features predicting relapse among 34 t(8;21)-rearranged patients. We found that the DNA methylation status of patients who suffered from relapse was peculiarly different from that of children maintaining complete remission. The epigenetic signature, made up of 337 differentially methylated regions, was then integrated with gene and protein expression profiles, leading to a network, where cell-to-cell adhesion and cell-motility pathways were found to be aberrantly activated in relapsed patients. We identified most of these factors as RUNX1-RUNX1T1 targets, with Ras Homolog Family Member (RHOB) overexpression being the core of this network. We documented how RHOB re-organized the actin cytoskeleton through its downstream ROCK-LIMK-COFILIN axis: this increases blast adhesion by stress fiber formation, and reduces mitochondrial apoptotic cell death after chemotherapy treatment. Altogether, our data show an epigenetic heterogeneity within t(8;21)-rearranged AML patients at diagnosis able to influence the program of the chimeric transcript, promoting blast re-emergence and progression to relapse.

miR-194 as predictive biomarker of responsiveness to neoadjuvant chemoradiotherapy in patients with locally advanced rectal adenocarcinoma.

AIMS: Curative surgery remains the primary form of treatment for locally advanced rectal cancer (LARC). Recent data support the use of preoperative chemoradiotherapy (pCRT) to improve the prognosis of LARC with a significant reduction of local relapse and an increase of overall survival. Unfortunately, only 20% of the patients with LARC present complete pathological response after pCRT, whereas in 20%-40%, the response is poor or absent. METHODS: We investigated the expression level of miR-194 in n=38 patients with LARC using our public microRNA (miRNA) expression dataset. miR-194 expression was further validated by real-time quantitative PCR (qRT-PCR) and in situ hybridisation (ISH). Protein-protein interaction network and pathway enrichment analysis were performed on miR-194 targets. RESULTS AND DISCUSSION: Using biopsy samples collected at diagnosis, mir-194 was significantly upregulated in patients responding to treatment (p value=0.016). The data was confirmed with qRT-PCR (p value=0.0587) and ISH (p value=0.026). Protein-protein interaction network and pathway enrichment analysis reveal a possible mechanism of susceptibility to pCRT involving Wnt pathway via its downstream mediator TRAF6. Finally, we interrogated the Comparative Toxicogenomics Database database in order to identify those chemical compounds able to mimic the biological effects of miR-194 as new possible therapeutic option in LARC treatment. The present study combining miRNA expression profiling with integrative computational biology identified miR-194 as predictive biomarker of response to pCRT. Using known and predicted drug mechanism of action, we then identified possible chemical compounds for further in vitro validation.

TP-0903 inhibits neuroblastoma cell growth and enhances the sensitivity to conventional chemotherapy.

Neuroblastoma (NB) is an embryonal tumor with low cure rate for patients classified as high-risk. This class of NB tumors shows a very complex genomic background and requires aggressive treatment strategies. In this work we evaluated the efficacy of the novel multi-kinase inhibitor TP-0903 in impairing NB cells' growth, proliferation and motility. In vitro studies were performed using cell lines with different molecular background, and in vivo studies were done using the zebrafish experimental model. Our results confirmed a strong cytotoxicity of TP-0903 already at the sub-micro molar concentrations. The observed cytotoxicity of TP-0903 was irreversible and the resulting apoptosis was caspase dependent. In addition, TP-0903 impaired colony formation and neurosphere creation. Depending on the molecular background of the selected NB cell lines, TP-0903 influenced either their capacity to migrate, to complete their cell cycle or both. Likewise, TP-0903 reduced NB cells intravasation in vitro and in vivo. Importantly, TP-0903 showed remarkable pharmacological efficacy not only as a mono-treatment, but also in combination with conventional chemotherapy drugs (ATRA, cisplatin, and VP16) in different types of NB cells. In conclusion, the multi-kinase activity of TP-0903 allowed the impairment of several biological processes required for expansion of NB cells, making them more vulnerable to the conventional chemotherapeutics. Altogether, our results support the eligibility of TP-0903 for further (pre)clinical assessments in NB.

Transcription instability in high-risk neuroblastoma is associated with a global perturbation of chromatin domains.

Chromosome instability has a pivotal role among the hallmarks of cancer, but its transcriptional counterpart is rarely considered a relevant factor in cell destabilization. To examine transcription instability (TIN), we first devised a metric we named TIN index and used it to evaluate TIN on a dataset containing more than 500 neuroblastoma samples. We found that metastatic tumors from high-risk (HR) patients are characterized by significantly different TIN index values compared to low/intermediate-risk patients. Our results indicate that the TIN index is a good predictor of neuroblastoma patient's outcome, and a related TIN index gene signature (TIN-signature) is also able to predict the neuroblastoma patient's outcome with high confidence. Interestingly, we find that TIN-signature genes have a strong positional association with superenhancers in neuroblastoma tumors. Finally, we show that TIN is linked to chromatin structural domains and interferes with their integrity in HR neuroblastoma patients. This novel approach to gene expression analysis broadens the perspective of genome instability investigations to include functional aspects.

Serum miR-125b is a non-invasive predictive biomarker of the pre-operative chemoradiotherapy responsiveness in patients with rectal adenocarcinoma.

BACKGROUND: Therapeutic management of Locally Advanced Rectal Cancer (LARC) involves pre-operative chemoradiotherapy (pCRT) followed by surgery. However, after pCRT the complete pathological response is approximately 20%, whereas in 20 to 40% of patients the response is poor or absent. METHODS: Cancer biopsy specimens (n= 38) and serum samples (n= 34) obtained before pCRT from 38 LARC patients were included in the study. Patients were classified in responders (R, tumor regression grade [TRG] 1-2; n= 16) and non-responders (NR, TRG 3-5; n= 22) according to the pathological response observed upon surgery. We performed miRNA microarrays analysis on biopsy specimens, and validated the selected candidates both by qRT-PCR (tissue and serum) and by in situ hybridization (tissue, miR-125b) analyses. RESULTS: Eleven miRNAs were significantly different between R and NR (miR-154, miR-409-3p, miR-127-3p, miR-214*, miR-299-5p and miR-125b overexpressed in NR; miR-33a, miR-30e, miR-338-3p, miR-200a and miR-378 decreased). In particular, miR-125b resulted to be the best candidate to discriminate the two groups (AUC of 0.9026; 95% CI, 0.7618-1.043). Additionally, miR-125b serum levels were significantly overexpressed in NR patients compared to R (p-value=0.0087), with an excellent discriminating power (AUC of 0.782; 95% CI, 0.6123-0.9518). CONCLUSIONS: The obtained results further support the clinical impact of miRNA analysis. High miR-125b expression in tissue and serum were associated with a poor treatment response in LARC patients, therefore miR-125b could be considered as a possible novel non-invasive biomarker of response in LARC treatment.

Exome and deep sequencing of clinically aggressive neuroblastoma reveal somatic mutations that affect key pathways involved in cancer progression.

The spectrum of somatic mutation of the most aggressive forms of neuroblastoma is not completely determined. We sought to identify potential cancer drivers in clinically aggressive neuroblastoma.Whole exome sequencing was conducted on 17 germline and tumor DNA samples from high-risk patients with adverse events within 36 months from diagnosis (HR-Event3) to identify somatic mutations and deep targeted sequencing of 134 genes selected from the initial screening in additional 48 germline and tumor pairs (62.5% HR-Event3 and high-risk patients), 17 HR-Event3 tumors and 17 human-derived neuroblastoma cell lines.We revealed 22 significantly mutated genes, many of which implicated in cancer progression. Fifteen genes (68.2%) were highly expressed in neuroblastoma supporting their involvement in the disease. CHD9, a cancer driver gene, was the most significantly altered (4.0% of cases) after ALK.Other genes (PTK2, NAV3, NAV1, FZD1 and ATRX), expressed in neuroblastoma and involved in cell invasion and migration were mutated at frequency ranged from 4% to 2%.Focal adhesion and regulation of actin cytoskeleton pathways, were frequently disrupted (14.1% of cases) thus suggesting potential novel therapeutic strategies to prevent disease progression.Notably BARD1, CHEK2 and AXIN2 were enriched in rare, potentially pathogenic, germline variants.In summary, whole exome and deep targeted sequencing identified novel cancer genes of clinically aggressive neuroblastoma. Our analyses show pathway-level implications of infrequently mutated genes in leading neuroblastoma progression.

Combating autophagy is a strategy to increase cytotoxic effects of novel ALK inhibitor entrectinib in neuroblastoma cells.

Neuroblastoma (NB) is a threatening childhood malignancy. Its prognosis is affected by several morphological, and biological characteristics, including the constitutive expression of ALK tyrosine kinase. In this study we examined the therapeutic potential of a novel ALK inhibitor, entrectinib, in obliterating NB tumor cells. Entrectinib showed the growth-inhibitory effects on NB cells with a 50% inhibitory concentration range of 0.03-5 µM. In the ALK-dependent cells, entrectinib mediated G1-arrest, which was associated with modified expression of multiple cell-cycle regulators. Down-regulation of Ki-67, and attenuated phosphorylation of ERK1/2, and STAT3, correlated with observed antiproliferative capacity of entrectinib. Initial cytostatic activity of entrectinib was followed by concentration-dependent apoptotic cell death, and Caspase-3 activation. However, we delineated a reduced sensitivity of ALK mutated NB cells to entrectinib, and demonstrated strong activation of autophagy in SH-SY5YF1174L NB cell line. Abrogation of autophagy by chloroquine increased significantly the toxicity of entrectinib, as confirmed by enhanced death rate, and PARP protein cleavage in SH-SY5YF1174L cells. In aggregate, our data show that entrectinib inhibits proliferation, and induces G1-arrest, and apoptosis in NB cells. We propose entrectinib for further consideration in treatment of NB, and recommend pharmacological inhibition of autophagy to be explored for a combined therapeutic approach in NB patients that might develop resistance to entrectinib.

An integrative approach for the identification of prognostic and predictive biomarkers in rectal cancer.

INTRODUCTION: Colorectal cancer is the third most common cancer in the world, a small fraction of which is represented by locally advanced rectal cancer (LARC). If not medically contraindicated, preoperative chemoradiotherapy, represent the standard of care for LARC patients. Unfortunately, patients shows a wide range of response rates in which approximately 20% has a complete pathological response, whereas in 20 to 40% the response is poor or absent. RESULTS: The following specific gene signature, able to discriminate responders' patients from non-responders, were founded: AKR1C3, CXCL11, CXCL10, IDO1, CXCL9, MMP12 and HLA-DRA. These genes are mainly involved in immune system pathways and interact with drugs traditionally used in the adjuvant treatment of rectal cancer. DISCUSSION: The present study suggests that new ideas for therapy could be found not only limited to studying genes differentially expressed between the two groups of patients but deepening the mechanisms, associated to response, in which they are involved. METHODS: Gene expression studies performed by: Agostini et al., Rimkus et al. and Kim et al. have been merged through a meta-analysis of the raw data. Gene expression data-sets have been processed using A-MADMAN. Common differentially expressed gene (DEG) were identified through SAM analysis. To further characterize the identified DEG we deeply investigated its biological role using an integrative computational biology approach.