Time-course gene expression experiments are useful tools for exploring biological processes. In this type of experiments, gene expression changes are monitored along time. Unfortunately, replication of time series is still costly and usually long time course do not have replicates. Many approaches have been proposed to deal with this data structure, but none of them in the field of pathway analysis. Pathway analyses have acquired great relevance for helping the interpretation of gene expression data. Several methods have been proposed to this aim: from the classical enrichment to the more complex topological analysis that gains power from the topology of the pathway. None of them were devised to identify temporal variations in time course data.

The production rate of gene expression data is nothing less than astounding. However, with the benefit of hindsight we can assert that, since we completely ignored the non-coding part of the transcriptome, we spent the last decade to study cell mechanisms having few data in our hands. In this scenario, microRNAs, which are key post-trascriptional regulators, deserve special attention. Given the state of knowledge about their biogenesis, mechanisms of action and the numerous experimentally validated target genes, miRNAs are also gradually appearing in the formal pathway representations such as KEGG and Reactome maps. However, the number of miRNAs annotated in pathway maps are very few and pathway analyses exploiting this new regulatory layer are still lacking. To fill these gaps, we present 'micrographite' a new pipeline to perform topological pathway analysis integrating gene and miRNA expression profiles. Here, micrographite is used to study and dissect the epithelial ovarian cancer gene and miRNA transcriptome defining and validating a new regulatory circuit related to ovarian cancer histotype specificity.

The existence of cancer stem cells (CSCs) within a tumor bulk has been demonstrated for many solid tumors including epithelial ovarian carcinoma (EOC). CSCs have been associated to tumor invasion, metastasis and development of chemoresistant recurrences. In this context, we aim to characterize EOC CSCs from the molecular point of view in order to identify potential biomarkers associated with chemoresistance.

RNA-Seq allows a theoretically unbiased analysis of both genome-wide transcription levels and mutation status of a tumor. Using this technique we sought to identify novel candidate therapeutic targets expressed in epithelial ovarian cancer (EOC).

During myogenesis, myoblasts fuse into multinucleated myotubes that acquire the contractile fibrils and accessory structures typical of striated skeletal muscle fibers. To support the high energy requirements of muscle contraction, myogenesis entails an increase in mitochondrial (mt) mass with stimulation of mtDNA synthesis and consumption of DNA precursors (dNTPs). Myotubes are quiescent cells and as such down-regulate dNTP production despite a high demand for dNTPs. Although myogenesis has been studied extensively, changes in dNTP metabolism have not been examined specifically. In differentiating cultures of C2C12 myoblasts and purified myotubes, we analyzed expression and activities of enzymes of dNTP biosynthesis, dNTP pools, and the expansion of mtDNA. Myotubes exibited pronounced post-mitotic modifications of dNTP synthesis with a particularly marked down-regulation of de novo thymidylate synthesis. Expression profiling revealed the same pattern of enzyme down-regulation in adult murine muscles. The mtDNA increased steadily after myoblast fusion, turning over rapidly, as revealed after treatment with ethidium bromide. We individually down-regulated p53R2 ribonucleotide reductase, thymidine kinase 2, and deoxyguanosine kinase by siRNA transfection to examine how a further reduction of these synthetic enzymes impacted myotube development. Silencing of p53R2 had little effect, but silencing of either mt kinase caused 50% mtDNA depletion and an unexpected decrease of all four dNTP pools independently of the kinase specificity. We suggest that during development of myotubes the shortage of even a single dNTP may affect all four pools through dysregulation of ribonucleotide reduction and/or dissipation of the non-limiting dNTPs during unproductive elongation of new DNA chains.

The application of high-throughput genomic tools in nutrition research is a widespread practice. However, it is becoming increasingly clear that the outcome of individual expression studies is insufficient for the comprehensive understanding of such a complex field. Currently, the availability of the large amounts of expression data in public repositories has opened up new challenges on microarray data analyses. We have focused on PPARalpha, a ligand-activated transcription factor functioning as fatty acid sensor controlling the gene expression regulation of a large set of genes in various metabolic organs such as liver, small intestine or heart. The function of PPARalpha is strictly connected to the function of its target genes and, although many of these have already been identified, major elements of its physiological function remain to be uncovered. To further investigate the function of PPARalpha, we have applied a cross-species meta-analysis approach to integrate sixteen microarray datasets studying high fat diet and PPARalpha signal perturbations in different organisms.

Polyglutamine expansion in the androgen receptor (AR) causes spinobulbar muscular atrophy (SBMA). Skeletal muscle is a primary site of toxicity; however, the current understanding of the early pathological processes that occur and how they unfold during disease progression remains limited. Using transgenic and knock-in mice and patient-derived muscle biopsies, we show that SBMA mice in the presymptomatic stage develop a respiratory defect matching defective expression of genes involved in excitation-contraction coupling (ECC), altered contraction dynamics, and increased fatigue. These processes are followed by stimulus-dependent accumulation of calcium into mitochondria and structural disorganization of the muscle triads. Deregulation of expression of ECC genes is concomitant with sexual maturity and androgen raise in the serum. Consistent with the androgen-dependent nature of these alterations, surgical castration and AR silencing alleviate the early and late pathological processes. These observations show that ECC deregulation and defective mitochondrial respiration are early but reversible events followed by altered muscle force, calcium dyshomeostasis, and dismantling of triad structure.

Patients with colorectal liver metastasis (CLM) present with heterogenous clinical outcomes and improved classification is needed to ameliorate the therapeutic output. Macrophages (Mϕ) hold promise as prognostic classifiers and therapeutic targets. Here, stemming from a single-cell analysis of mononuclear phagocytes infiltrating human CLM, we identified two Mϕ markers associated with distinct populations with opposite clinical relevance. The invasive margin of CLM was enriched in pro-inflammatory monocyte-derived Mϕ (MoMϕ) expressing the monocytic marker SERPINB2, and a more differentiated population, tumor-associated Mϕ (TAM), expressing glycoprotein nonmetastatic melanoma protein B (GPNMB). SERPINB2+ MoMϕ had an early inflammatory profile, whereas GPNMB+ TAMs were enriched in pathways of matrix degradation, angiogenesis, and lipid metabolism and were found closer to the tumor margin, as confirmed by spatial transcriptomics on CLM specimens. In a cohort of patients, a high infiltration of SERPINB2+ cells independently associated with longer disease-free survival (DFS; P = 0.033), whereas a high density of GPNMB+ cells correlated with shorter DFS (P = 0.012) and overall survival (P = 0.002). Cell-cell interaction analysis defined opposing roles for MoMϕ and TAMs, suggesting that SERPINB2+ and GPNMB+ cells are discrete populations of Mϕ and may be exploited for further translation to an immune-based stratification tool. This study provides evidence of how multi-omics approaches can identify nonredundant, clinically relevant markers for further translation to immune-based patient stratification tools and therapeutic targets. GPNMB has been shown to set Mϕ in an immunosuppressive mode. Our high dimensional analyses provide further evidence that GPNMB is a negative prognostic indicator and a potential player in the protumor function of Mϕ populations.

We report here for the first time, a comprehensive characterization of biological and clinical features of early-stage triple negative Invasive Lobular Carcinomas(TN-ILCs) METHODS: We analyzed all consecutive patients with early-stage TN-ILC operated at two reference cancer-centers between 1994 and 2012. Primary objective was to assess the invasive disease-free survival(iDFS). Co-primary objective was to assess biological features of TN-ILCs, including molecular intrinsic subtypes based on PAM-50 assay, expression of androgen receptor (AR) and mutational status of ERBB2-gene. Additionally, DNA mutational status of an independent cohort of 45 TN-ILCs from three databases were analyzed, to confirm mutations in ERBB2-gene and to identify other recurrently mutated genes.

Physical activity improves the prognosis of cancer patients, partly by contrasting the associated muscle wasting (cachexia), through still unknown mechanisms. We asked whether aerobic exercise causes secretion by skeletal muscles of proteins (myokines) that may contrast cachexia. Media conditioned by peroxisome proliferator-activated receptor γ coactivator 1α (PGC1α)-expressing myotubes, reproducing some metabolic adaptations of aerobic exercise, as increased mitochondrial biogenesis and oxidative phosphorylation, restrained constitutively active Forkhead box-containing subfamily O3 (caFoxO3)-induced proteolysis. Microarray analysis identified amphiregulin (AREG), natriuretic peptide precursor B (NppB), musclin and fibroblast growth factor 18 (FGF18) as myokines highly induced by PGC1α. Notably, only musclin tended to be low in muscle of mice with a rare human renal carcinoma; it was reduced in plasma and in muscles of C26-bearing mice and in atrophying myotubes, where PGC1α expression is impaired. Therefore, we electroporated the Tibialis Anterior (TA) of C26-bearing mice with musclin or (its receptor) natriuretic peptide receptor 3 (Npr3)-encoding plasmids and found a preserved fiber area, as a result of restrained proteolysis. Musclin knockout (KO) mice lose more muscle tissue during growth of two distinct cachexia-causing tumors. Running protected C26-bearing mice from cachexia, not changing tumor growth, and rescued the C26-induced downregulation of musclin in muscles and plasma. Musclin expression did not change in overloaded plantaris of mice, recapitulating partially muscle adaptations to anaerobic exercise. Musclin might, therefore, be beneficial to cancer patients who cannot exercise and are at risk of cachexia and may help to explain how aerobic exercise alleviates cancer-induced muscle wasting.

High-grade serous ovarian carcinoma (HGSOC) is generally associated with a very dismal prognosis. Nevertheless, patients with similar clinicopathological characteristics can have markedly different clinical outcomes. Our aim was the identification of novel molecular determinants influencing survival.

Background: Dimethyl fumarate (DMF) is a disease-modifying drug for relapsing-remitting multiple sclerosis. Among others, DMF impedes immune activation by shifting the balance between inflammatory and regulatory cell types and by inducing apoptosis-triggered lymphopenia. Although the decrease in lymphocyte count is an early effect of the drug in several patients, the long-term impact on lymphocyte subsets is largely unknown. Methods: We performed a 2-years observational study on total lymphocyte count and subsets thereof by flow cytometry of peripheral blood of 38 multiple sclerosis patients in treatment with DMF. Data were collected at the beginning and after 3, 6, 12, and 24 months of therapy. Results: Total lymphocyte count decreased in relation to time of exposure to DMF. Mean absolute B cell count decreased by 34.1% (p < 0.001) within the first 3 months of therapy and then remained stable over time. Mean absolute CD3+ T cells count decrement reached 47.5% after 12 months of treatment (p < 0.001). NK cells count showed a heterogeneous trend, increasing by 85.9% (p < 0.001) after 2 years of treatment. CD4+ T cells and CD8+ T cells substantially decreased, with a significant increase of CD4+/CD8+ ratio during the first year of therapy. Conclusions: NK cells showed a heterogeneous behavior during DMF treatment with a significant increase over time. Since NK cells may also have a regulatory effect on immune system modulation, their increase during DMF treatment might play a role in the efficacy and safety of the drug.

Skeletal muscle is composed of different myofiber types that preferentially use glucose or lipids for ATP production. How fuel preference is regulated in these post-mitotic cells is largely unknown, making this issue a key question in the fields of muscle and whole-body metabolism. Here, we show that microRNAs (miRNAs) play a role in defining myofiber metabolic profiles. mRNA and miRNA signatures of all myofiber types obtained at the single-cell level unveiled fiber-specific regulatory networks and identified two master miRNAs that coordinately control myofiber fuel preference and mitochondrial morphology. Our work provides a complete and integrated mouse myofiber type-specific catalog of gene and miRNA expression and establishes miR-27a-3p and miR-142-3p as regulators of lipid use in skeletal muscle.

Dysregulation of miRNAs expression plays a critical role in the pathogenesis of genetic, multifactorial disorders and in human cancers. We exploited sequence, genomic and expression information to investigate two main aspects of post-transcriptional regulation in miRNA biogenesis, namely strand selection regulation and expression relationships between intragenic miRNAs and host genes. We considered miRNAs expression profiles, measured in five sizeable microarray datasets, including samples from different normal cell types and tissues, as well as different tumours and disease states. First, the study of expression profiles of "sister" miRNA pairs (miRNA/miRNA*, 5' and 3' strands of the same hairpin precursor) showed that the strand selection is highly regulated since it shows tissue-/cell-/condition-specific modulation. We used information about the direction and the strength of the strand selection bias to perform an unsupervised cluster analysis for the sample classification evidencing that is able to distinguish among different tissues, and sometimes between normal and malignant cells. Then, considering a minimum expression threshold, in few miRNA pairs only one mature miRNA is always present in all considered cell types, whereas the majority of pairs were concurrently expressed in some cell types and alternatively in others. In a significant fraction of concurrently expressed pairs, the major and the minor forms found at comparable levels may contribute to post-transcriptional gene silencing, possibly in a coordinate way. In the second part of the study, the behaved tendency to co-expression of intragenic miRNAs and their "host" mRNA genes was confuted by expression profiles examination, suggesting that the expression profile of a given host gene can hardly be a good estimator of co-transcribed miRNA(s) for post-transcriptional regulatory networks inference. Our results point out the regulatory importance of post-transcriptional phases of miRNAs biogenesis, reinforcing the role of such layer of miRNA biogenesis in miRNA-based regulation of cell activities.

Systematic efforts to sequence the cancer genome have identified large numbers of mutations and copy number alterations in human cancers. However, elucidating the functional consequences of these variants, and their interactions to drive or maintain oncogenic states, remains a challenge in cancer research. We developed REVEALER, a computational method that identifies combinations of mutually exclusive genomic alterations correlated with functional phenotypes, such as the activation or gene dependency of oncogenic pathways or sensitivity to a drug treatment. We used REVEALER to uncover complementary genomic alterations associated with the transcriptional activation of β-catenin and NRF2, MEK-inhibitor sensitivity, and KRAS dependency. REVEALER successfully identified both known and new associations, demonstrating the power of combining functional profiles with extensive characterization of genomic alterations in cancer genomes.

Epithelial ovarian cancer (EOC) is a spectrum of different diseases, which makes their treatment a challenge. Forkhead box M1 (FOXM1) is an oncogene aberrantly expressed in many solid cancers including serous EOC, but its role in non-serous EOCs remains undefined. We examined FOXM1 expression and its correlation to prognosis across the three major EOC subtypes, and its role in tumorigenesis and chemo-resistance in vitro.

Many models and analysis of signaling pathways have been proposed. However, neither of them takes into account that a biological pathway is not a fixed system, but instead it depends on the organism, tissue and cell type as well as on physiological, pathological and experimental conditions.

MAGIA (miRNA and genes integrated analysis) is a novel web tool for the integrative analysis of target predictions, miRNA and gene expression data. MAGIA is divided into two parts: the query section allows the user to retrieve and browse updated miRNA target predictions computed with a number of different algorithms (PITA, miRanda and Target Scan) and Boolean combinations thereof. The analysis section comprises a multistep procedure for (i) direct integration through different functional measures (parametric and non-parametric correlation indexes, a variational Bayesian model, mutual information and a meta-analysis approach based on P-value combination) of mRNA and miRNA expression data, (ii) construction of bipartite regulatory network of the best miRNA and mRNA putative interactions and (iii) retrieval of information available in several public databases of genes, miRNAs and diseases and via scientific literature text-mining. MAGIA is freely available for Academic users at http://gencomp.bio.unipd.it/magia.

Topological gene-set analysis has emerged as a powerful means for omic data interpretation. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. Here, we propose a new method, called SourceSet, able to distinguish between the primary and the secondary dysregulation within a Gaussian graphical model context. The proposed method compares gene expression profiles in the control and in the perturbed condition and detects the differences in both the mean and the covariance parameters with a series of likelihood ratio tests. The resulting evidence is used to infer the primary and the secondary set, i.e. the genes responsible for the primary dysregulation, and the genes affected by the perturbation through network propagation. The proposed method demonstrates high specificity and sensitivity in different simulated scenarios and on several real biological case studies. In order to fit into the more traditional pathway analysis framework, SourceSet R package also extends the analysis from a single to multiple pathways and provides several graphical outputs, including Cytoscape visualization to browse the results.

We recently discovered that forebrain activation of the IL-1 receptor/Toll-like receptor (IL-1R1/TLR4) innate immunity signal plays a pivotal role in neuronal hyperexcitability underlying seizures in rodents. Since this pathway is activated in neurons and glia in human epileptogenic foci, it represents a potential target for developing drugs interfering with the mechanisms of epileptogenesis that lead to spontaneous seizures. The lack of such drugs represents a major unmet clinical need. We tested therefore novel therapies inhibiting the IL-1R1/TLR4 signaling in an established murine model of acquired epilepsy. We used an epigenetic approach by injecting a synthetic mimic of micro(mi)RNA-146a that impairs IL1R1/TLR4 signal transduction, or we blocked receptor activation with antiinflammatory drugs. Both interventions when transiently applied to mice after epilepsy onset, prevented disease progression and dramatically reduced chronic seizure recurrence, while the anticonvulsant drug carbamazepine was ineffective. We conclude that IL-1R1/TLR4 is a novel potential therapeutic target for attaining disease-modifications in patients with diagnosed epilepsy.