Search

Objective: Temporal lobe epilepsy (TLE) is the most common form of refractory focal epilepsy, and the current clinical diagnosis is based on EEG, clinical neurological history and neuroimaging findings. Methods: So far, there are no blood-based molecular biomarkers of TLE to support clinical diagnosis, despite the pathogenic mechanisms underlying TLE involving defects in the regulation of gene expression. MicroRNAs (miRNAs) have emerged as important post-transcriptional regulators of gene expression. Results: Recent studies show the feasibility of detecting miRNAs in body fluids; circulating miRNAs have emerged as potential clinical biomarkers in epilepsy, although the TLE miRNA profile needs to be addressed. Here, we analysed the diagnostic potential of 8 circulating miRNAs in sera of 52 TLE patients and 40 age- and sex-matched donor controls by RT-qPCR analyses. Conclusion: We found that miR-34a-5p, -106b-5p, -130a-3p, -146a-5p, and -19a-3p are differently expressed in TLE compared to control subjects, suggesting a diagnostic role. Furthermore, we found that miR-34a-5p, -106b-5p, -146a-5p and miR-451a could become prognostic biomarkers, being differentially expressed between drug-resistant and drug-responsive TLE subjects. Therefore, serum miRNAs are diagnostic and drug-resistance predictive molecules of TLE.

Modern society grew rapidly over the last few decades and this led to an alarming increase in air pollutants and a worsening of the human health, especially in relation to the respiratory system. Indeed, chronic respiratory diseases were the third main cause of death in 2017, with over 3 million of deaths. Furthermore, the pollution has considerable consequences both for burden medical expenses and environmental. However, the mechanisms linking pollutants to the onset of these diseases remain unclear. Thus, in this study we addressed this problem through the United Kingdom BioBank database, analyzing 170 genome-wide association studies (103 related to respiratory diseases and 67 related to pollutants). We analyzed the genetic correlations and causal relationships of these traits, leveraging the summary statistics and bioinformatics packages such as Linkage Disequilibrium Score Regression and Latent Causal Variable. We obtained 158 significant genetic correlations and subsequently we analyzed them through the Latent Causal Variable analysis, obtaining 20 significant causal relationships. The most significant were between "Workplace full of chemicals or other fumes: Sometimes" and “Condition that has ever been diagnosed by a doctor: Asthma” and between “Workplace very dusty: Sometimes” and “Condition that has ever been diagnosed by a doctor: Emphysema or chronic bronchitis”. Finally, we identified single nucleotide polymorphisms independently associated with sveral pollutants to analyze the genes and pathways that could be involved in the onset of the aforementioned respiratory system disorders and that could be useful clinical target. This study highlighted how crucial are the air condition of the working environments and the type of transport used in the onset of respiratory-related morbidity. Based on that, we also suggested some interventions, in order to improve quality life and develop new and eco-friendly society and life style, such as improving indoor air cir

Sporadic Alzheimer’s disease (AD) is a complex neurological disorder characterized by many risk loci with potential associations with different traits and diseases. AD, characterized by a progressive loss of neuronal functions, manifests with different symptoms such as decline in memory, movement, coordination, and speech. The mechanisms underlying the onset of AD are not always fully understood, but involve a multiplicity of factors. Early diagnosis of AD plays a central role as it can offer the possibility of early treatment, which can slow disease progression. Currently, the methods of diagnosis are cognitive testing, neuroimaging, or cerebrospinal fluid analysis that can be time-consuming, expensive, invasive, and not always accurate. In the present study, we performed a genetic correlation analysis using genome-wide association statistics from a large study of AD and UK Biobank, to examine the association of AD with other human traits and disorders. In addition, since hippocampus, a part of cerebral cortex could play a central role in several traits that are associated with AD; we analyzed the gene expression profiles of hippocampus of AD patients applying 4 different artificial neural network models. We found 65 traits correlated with AD grouped into 9 clusters: medical conditions, fluid intelligence, education, anthropometric measures, employment status, activity, diet, lifestyle, and sexuality. The comparison of different 4 neural network models along with feature selection methods on 5 Alzheimer’s gene expression datasets showed that the simple basic neural network model obtains a better performance (66% of accuracy) than other more complex methods with dropout and weight regularization of the network.

Leveraging genome-wide association statistics generated from a large study of amyotrophic lateral sclerosis (ALS; 29,612 cases and 122,656 controls) and UK Biobank (UKB; 4,024 phenotypes, up to 361,194 participants), we conducted a phenome-wide analysis of ALS genetic liability and identified 46 genetically correlated traits, such as fluid intelligence score (rg = − 0.21, p = 1.74 × 10–6), "spending time in pub or social club” (rg = 0.24, p = 2.77 × 10–6), non-work related walking (rg = − 0.25, p = 1.95 × 10–6), college education (rg = − 0.15, p = 7.08 × 10–5), “ever diagnosed with panic attacks (rg = 0.39, p = 4.24 × 10–5), and “self-reported other gastritis including duodenitis” (rg = 0.28, p = 1.4 × 10–3). To assess the putative directionality of these genetic correlations, we conducted a latent causal variable analysis, identifying significant genetic causality proportions (gĉp) linking ALS genetic liability to seven traits. While the genetic component of “self-reported other gastritis including duodenitis" showed a causal effect on ALS (gĉp = 0.50, p = 1.26 × 10–29), the genetic liability to ALS is potentially causal for multiple traits, also including an effect on "ever being diagnosed with panic attacks” (gĉp = 0.79, p = 5.011 × 10–15) and inverse effects on “other leisure/social group activities” (gĉp = 0.66, p = 1 × 10–4) and prospective memory result (gĉp = 0.35, p = 0.005). Our subsequent Mendelian randomization analysis indicated that some of these associations may be due to bidirectional effects. In conclusion, this phenome-wide investigation of ALS polygenic architecture highlights the widespread pleiotropy linking this disorder with several health domains.

Prostate cancer (PC) is a male common neoplasm and is the second leading cause of cancer death in American men. PC is traditionally diagnosed by the evaluation of prostate secreted antigen (PSA) in the blood. Due to the high levels of false positives, digital rectal examination and transrectal ultrasound guided biopsy are necessary in uncertain cases with elevated PSA levels. Nevertheless, the high mortality rate suggests that new PC biomarkers are urgently needed to help clinical diagnosis. In a previous study, we have identified a network of genes, altered in high Gleason Score (GS) PC (GS ≥ 7), being regulated by miR‐153. Until now, no publication has explained the mechanism of action of miR‐153 in PC. By in vitro studies, we found that the overexpression of miR‐ 153 in high GS cell lines is required to control cell proliferation, migration and invasion rates, targeting Kruppel‐like factor 5 (KLF5). Moreover, miR‐153 could be secreted by exosomes and microvesicles in the microenvironment and, once entered into the surrounding tissue, could influence cellular growth. Being upregulated in high GS human PC, miR‐153 could be proposed as a circulating biomarker for PC diagnosis.

Alzheimer’s disease (AD) is a neurodegenerative disease characterized by rapid brain cell degeneration affecting different areas of the brain. Hippocampus is one of the earliest involved brain regions in the disease. Modern technologies based on high-throughput data have identified transcriptional profiling of several neurological diseases, including AD, for a better comprehension of genetic mechanisms of the disease. In this study, we investigated differentially expressed genes (DEGs) from six Gene Expression Omnibus (GEO) datasets of hippocampus of AD patients. The identified DEGs were submitted to Weighted correlation network analysis (WGCNA) and ClueGo to explore genes with a higher degree centrality and to comprehend their biological role. Subsequently, MCODE was used to identify subnetworks of interconnected DEGs. Our study found 40 downregulated genes and 36 up-regulated genes as consensus DEGs. Analysis of the co-expression network revealed ACOT7, ATP8A2, CDC42, GAD1, GOT1, INA, NCALD, and WWTR1 to be genes with a higher degree centrality. ClueGO revealed the pathways that were mainly enriched, such as clathrin coat assembly, synaptic vesicle endocytosis, and DNA damage response signal transduction by p53 class mediator. In addition, we found a subnetwork of 12 interconnected genes (AMPH, CA10, CALY, NEFL, SNAP25, SNAP91, SNCB, STMN2, SV2B, SYN2, SYT1, and SYT13). Only CA10 and CALY are targets of known drugs while the others could be potential novel drug targets.

Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disease characterized by the neurodegeneration of motoneurons. About 10% of ALS is hereditary and involves mutation in 25 different genes, while 90% of the cases are sporadic forms of ALS (sALS). The diagnosis of ALS includes the detection of early symptoms and, as disease progresses, muscle twitching and then atrophy spreads from hands to other parts of the body. The disease causes high disability and has a high mortality rate; moreover, the therapeutic approaches for the pathology are not effective. miRNAs are small non-coding RNAs, whose activity has a major impact on the expression levels of coding mRNA. The literature identifies several miRNAs with diagnostic abilities on sALS, but a unique diagnostic profile is not defined. As miRNAs could be secreted, the identification of specific blood miRNAs with diagnostic ability for sALS could be helpful in the identification of the patients. In the view of personalized medicine, we performed a meta-analysis of the literature in order to select specific circulating miRNAs with diagnostic properties and, by bioinformatics approaches, we identified a panel of 10 miRNAs (miR-193b, miR-3911, miR-139-5p, miR-193b-1, miR-338-5p, miR-3911-1, miR-455-3p, miR-4687-5p, miR-4745-5p, and miR-4763-3p) able to classify sALS patients by blood analysis. Among them, the analysis of expression levels of the couple of blood miR-193b/miR-4745-5p could be translated in clinical practice for the diagnosis of sALS.

Cancer driver gene alterations influence cancer development, occurring in oncogenes, tumor suppressors, and dual role genes. Discovering dual role cancer genes is difficult because of their elusive context-dependent behavior. We define oncogenic mediators as genes controlling biological processes. With them, we classify cancer driver genes, unveiling their roles in cancer mechanisms. To this end, we present Moonlight, a tool that incorporates multiple -omics data to identify critical cancer driver genes. With Moonlight, we analyze 8000+ tumor samples from 18 cancer types, discovering 3310 oncogenic mediators, 151 having dual roles. By incorporating additional data (amplification, mutation, DNA methylation, chromatin accessibility), we reveal 1000+ cancer driver genes, corroborating known molecular mechanisms. Additionally, we confirm critical cancer driver genes by analysing cell-line datasets. We discover inactivation of tumor suppressors in intron regions and that tissue type and subtype indicate dual role status. These findings help explain tumor heterogeneity and could guide therapeutic decisions.

Amyotrophic lateral sclerosis (ALS) is a complex disease with a late onset and is characterized by the progressive loss of muscular and respiratory functions. Although recent studies have partially elucidated ALS’s mechanisms, many questions remain such as what the most important molecular pathways involved in ALS are and why there is such a large difference in ALS onset among different populations. In this study, we addressed this issue with a bioinformatics approach, using the United Kingdom Biobank (UKBB) and the European 1000 Genomes Project (1KG) in order to analyze the most ALS-representative single nucleotide polymorphisms (SNPs) that differ for minor allele frequency (MAF) between the United Kingdom population and some European populations including Finnish in Finland, Iberian population in Spain, and Tuscans in Italy. We found 84 SNPs associated with 46 genes that are involved in different pathways including: “Ca2+ activated K+ channels”, “cGMP effects”, ”Nitric oxide stimulates guanylate cyclase”, “Proton/oligopeptide co-transporters”, and “Signaling by MAPK mutants”. In addition, we revealed that 83% of the 84 SNPs can alter transcription factor-motives binding sites of 224 genes implicated in “Regulation of beta-cell development”, “Transcription-al regulation by RUNX3”, “Transcriptional regulation of pluripotent stem cells”, and “FOXO-mediated transcription of cell death genes”. In conclusion, the genes and pathways analyzed could explain the cause of the difference of ALS onset.

Background: There is extensive scientific evidence that radiation therapy (RT) is a crucial treatment, either alone or in combination with other treatment modalities, for many types of cancer, including breast cancer (BC). BC is a heterogeneous disease at both clinical and molecular levels, presenting distinct subtypes linked to the hormone receptor (HR) status and associated with different clinical outcomes. The aim of this study was to assess the molecular changes induced by high doses of ionizing radiation (IR) on immortalized and primary BC cell lines grouped according to Human epidermal growth factor receptor (HER2), estrogen, and progesterone receptors, to study how HR status influences the radiation response. Our genomic approach using in vitro and ex-vivo models (e.g., primary cells) is a necessary first step for a translational study to describe the common driven radio-resistance features associated with HR status. This information will eventually allow clinicians to prescribe more personalized total doses or associated targeted therapies for specific tumor subtypes, thus enhancing cancer radio-sensitivity. Methods: Nontumorigenic (MCF10A) and BC (MCF7 and MDA-MB-231) immortalized cell lines, as well as healthy (HMEC) and BC (BCpc7 and BCpcEMT) primary cultures, were divided into low grade, high grade, and healthy groups according to their HR status. At 24 h post-treatment, the gene expression profiles induced by two doses of IR treatment with 9 and 23 Gy were analyzed by cDNA microarray technology to select and compare the differential gene and pathway expressions among the experimental groups. Results: We present a descriptive report of the substantial alterations in gene expression levels and pathways after IR treatment in both immortalized and primary cell cultures. Overall, the IR-induced gene expression profiles and pathways appear to be cell-line dependent. The data suggest that some specific gene and pathway signatures seem to be linked to HR status

Background: Modern high-throughput genomic technologies represent a comprehensive hallmark of molecular changes in pan-cancer studies. Although different cancer gene signatures have been revealed, the mechanism of tumourigenesis has yet to be completely understood. Pathways and networks are important tools to explain the role of genes in functional genomic studies. However, few methods consider the functional non-equal roles of genes in pathways and the complex gene-gene interactions in a network. Results: We present a novel method in pan-cancer analysis that identifies de-regulated genes with a functional role by integrating pathway and network data. A pan-cancer analysis of 7158 tumour/normal samples from 16 cancer types identified 895 genes with a central role in pathways and de-regulated in cancer. Comparing our approach with 15 current tools that identify cancer driver genes, we found that 35.6% of the 895 genes identified by our method have been found as cancer driver genes with at least 2/15 tools. Finally, we applied a machine learning algorithm on 16 independent GEO cancer datasets to validate the diagnostic role of cancer driver genes for each cancer. We obtained a list of the top-ten cancer driver genes for each cancer considered in this study. Conclusions: Our analysis 1) confirmed that there are several known cancer driver genes in common among different types of cancer, 2) highlighted that cancer driver genes are able to regulate crucial pathways

Purpose: We demonstrated that Hsa-miR-567 expression is significantly downregulated in poor prognosis breast cancer, compared to better prognosis breast cancer, having a role in the control of cell proliferation and migration by regulating KPNA4 gene. Methods and results: In this study, based on our previously published in silico results, we proved both in vitro (cell line studies) and ex vivo (clinical studies), that Hsa-miR-567 expression is significantly downregulated in breast cancer with poor prognosis when compared to breast cancer with better prognosis. More intriguingly, we demonstrated that the ectopic expression of Hsa-miR-567 in poor prognosis breast cancer cell line strongly inhibits in vitro cell proliferation and migration. Furthermore, we showed in vivo that breast cancer cells, stably expressing Hsa-miR-567, xenografted in mouse, reduce tumor growth ability. Consistently, we found that karyopherin 4 (KPNA4), predicted target gene of Hsa-miR-567 as identified by our in silico analysis, is upregulated in highly aggressive MDA-MB-231 breast cancer cell line and patient tissues with poor prognosis with respect to good prognosis. Conclusions: Our results suggest a potential role of Hsa-miR-567 as a novel prognostic biomarker for BC and as regulator of KPNA4, Purpose: We demonstrated that Hsa-miR-567 expression is significantly downregulated in poor prognosis breast cancer, compared to better prognosis breast cancer, having a role in the control of cell proliferation and migration by regulating KPNA4 gene. Methods and results: In this study, based on our previously published in silico results, we proved both in vitro (cell line studies) and ex vivo (clinical studies), that Hsa-miR-567 expression is significantly downregulated in breast cancer with poor prognosis when compared to breast cancer with better prognosis. More intriguingly, we demonstrated that the ectopic expression of Hsa-miR-567 in poor prognosis breast cancer cell line strongly inhibits in vitro cell proliferation and migration. Furthermore, we showed in vivo that breast cancer cells, stably expressing Hsa-miR-567, xenografted in mouse, reduce tumor growth ability. Consistently, we found that karyopherin 4 (KPNA4), predicted target gene of Hsa-miR-567 as identified by our in silico analysis, is upregulated in highly aggressive MDA-MB-231 breast cancer cell line and patient tissues with poor prognosis with respect to good prognosis. Conclusions: Our results suggest a potential role of Hsa-miR-567 as a novel prognostic biomarker for BC and as regulator of KPNA4.

Prostate cancer (PC) includes several phenotypes, from indolent to highly aggressive cancer. Actual diagnostic and prognostic tools have several limitations, and there is a need for new biomarkers to stratify patients and assign them optimal therapies by taking into account potential genetic and epigenetic differences. MicroRNAs (miRNAs) are small sequences of non-coding RNA regulating specific genes involved in the onset and development of PC. Stable miRNAs have been found in biofluids, such as serum and plasma; thus, the measurement of PC-associated miRNAs is emerging as a non-invasive tool for PC detection and monitoring. In this study, we conduct an in-depth literature review focusing on miRNAs that may contribute to the diagnosis and prognosis of PC. The role of miRNAs as a potential theranostic tool in PC is discussed. Using a meta-analysis approach, we found a group of 29 miRNAs with diagnostic properties and a group of seven miRNAs with prognostic properties, which were found already expressed in both biofluids and PC tissues. We tested the two miRNA groups on The Cancer Genome Atlas dataset of PC tissue samples with a machine-learning approach. Our results suggest that these 29 miRNAs should be considered as potential panel of biomarkers for the diagnosis of PC, both as in vivo non-invasive test and ex vivo confirmation test.

The Cancer Genome Atlas (TCGA) research network has made public a large collection of clinical and molecular phenotypes of more than 10 000 tumor patients across 33 different tumor types. Using this cohort, TCGA has published over 20 marker papers detailing the genomic and epigenomic alterations associated with these tumor types. Although many important discoveries have been made by TCGA's research network, opportunities still exist to implement novel methods, thereby elucidating new biological pathways and diagnostic markers. However, mining the TCGA data presents several bioinformatics challenges, such as data retrieval and integration with clinical data and other molecular data types (e.g. RNA and DNA methylation). We developed an R/Bioconductor package called TCGAbiolinks to address these challenges and offer bioinformatics solutions by using a guided workflow to allow users to query, download and perform integrative analyses of TCGA data. We combined methods from computer science and statistics into the pipeline and incorporated methodologies developed in previous TCGA marker studies and in our own group. Using four different TCGA tumor types (Kidney, Brain, Breast and Colon) as examples, we provide case studies to illustrate examples of reproducibility, integrative analysis and utilization of different Bioconductor packages to advance and accelerate novel discoveries.

Dysregulation of microRNAs (miRNAs) is involved in the initiation and progression of several human cancers, including breast cancer (BC), as strong evidence has been found that miRNAs can act as oncogenes or tumor suppressor genes. This review presents the state of the art on the role of miRNAs in the diagnosis, prognosis, and therapy of BC. Based on the results obtained in the last decade, some miRNAs are emerging as biomarkers of BC for diagnosis (i.e., miR-9, miR-10b, and miR-17-5p), prognosis (i.e., miR-148a and miR-335), and prediction of therapeutic outcomes (i.e., miR-30c, miR-187, and miR-339-5p) and have important roles in the control of BC hallmark functions such as invasion, metastasis, proliferation, resting death, apoptosis, and genomic instability. Other miRNAs are of interest as new, easily accessible, affordable, non-invasive tools for the personalized management of patients with BC because they are circulating in body fluids (e.g., miR-155 and miR-210). In particular, circulating multiple miRNA profiles are showing better diagnostic and prognostic performance as well as better sensitivity than individual miRNAs in BC. New miRNA-based drugs are also promising therapy for BC (e.g., miR-9, miR-21, miR34a, miR145, and miR150), and other miRNAs are showing a fundamental role in modulation of the response to other non-miRNA treatments, being able to increase their efficacy (e.g., miR-21, miR34a, miR195, miR200c, and miR203 in combination with chemotherapy).

In this work an integrated approach was used to identify functional miRNAs regulating gene pathway cross-talk in breast cancer (BC). We first integrated gene expression profiles and biological pathway information to explore the underlying associations between genes differently expressed among normal and BC samples and pathways enriched from these genes. For each pair of pathways, a score was derived from the distribution of gene expression levels by quantifying their pathway cross-talk. Random forest classification allowed the identification of pairs of pathways with high cross-talk. We assessed miRNAs regulating the identified gene pathways by a mutual information analysis. A Fisher test was applied to demonstrate their significance in the regulated pathways. Our results suggest interesting networks of pathways that could be key regulatory of target genes in BC, including stem cell pluripotency, coagulation, and hypoxia pathways and miRNAs that control these networks could be potential biomarkers for diagnostic, prognostic, and therapeutic development in BC. This work shows that standard methods of predicting normal and tumor classes such as differentially expressed miRNAs or transcription factors could lose intrinsic features; instead our approach revealed the responsible molecules of the disease.

Background: Development of human cancer can proceed through the accumulation of different genetic changes affecting the structure and function of the genome. Combined analyses of molecular data at multiple levels, such as DNA copy-number alteration, mRNA and miRNA expression, can clarify biological functions and pathways deregulated in cancer. The integrative methods that are used to investigate these data involve different fields, including biology, bioinformatics, and statistics. Results: These methodologies are presented in this review, and their implementation in breast cancer is discussed with a focus on integration strategies. We report current applications, recent studies and interesting results leading to the identification of candidate biomarkers for diagnosis, prognosis, and therapy in breast cancer by using both individual and combined analyses. Conclusion: This review presents a state of art of the role of different technologies in breast cancer based on the integration of genetics and epigenetics, and shares some issues related to the new opportunities and challenges offered by the application of such integrative approaches.

Like other cancer diseases, prostate cancer (PC) is caused by the accumulation of genetic alterations in the cells that drives malignant growth. These alterations are revealed by gene profiling and copy number alteration (CNA) analysis. Moreover, recent evidence suggests that also microRNAs have an important role in PC development. Despite efforts to profile PC, the alterations (gene, CNA, and miRNA) and biological processes that correlate with disease development and progression remain partially elusive. Many gene signatures proposed as diagnostic or prognostic tools in cancer poorly overlap. The identification of co-expressed genes, that are functionally related, can identify a core network of genes associated with PC with a better reproducibility. By combining different approaches, including the integration of mRNA expression profiles, CNAs, and miRNA expression levels, we identified a gene signature of four genes overlapping with other published gene signatures and able to distinguish, in silico, high Gleason-scored PC from normal human tissue, which was further enriched to 19 genes by gene co-expression analysis. From the analysis of miRNAs possibly regulating this network, we found that hsa-miR-153 was highly connected to the genes in the network. Our results identify a four-gene signature with diagnostic and prognostic value in PC and suggest an interesting gene network that could play a key regulatory role in PC development and progression. Furthermore, hsa-miR-153, controlling this network, could be a potential biomarker for theranostics in high Gleason-scored PC.

Copy–number alterations (CNAs) represent an important component of genetic variations and play a significant role in many human diseases. Such alterations are related to certain types of cancers, including those of the pancreas, colon, and breast, among others. CNAs have been used as biomarkers for cancer prognosis in multiple studies, but few works report on the relation of CNAs with the disease progression. In this paper, we provide cases where the inference on the disease progression improves when exploiting CNA information. To this aim, a specific dissimilarity-based representation of patients is given. The employed framework outperforms a typical approach where patients are represented through a set of available attribute values. Three datasets were employed to validate the results of our analysis.

Gene Regulatory Networks (GRNs) control many biological systems, but how such network coordination is shaped is still unknown. GRNs can be subdivided into basic connections that describe how the network members interact e.g., co-expression, physical interaction, co-localization, genetic influence, pathways, and shared protein domains. The important regulatory mechanisms of these networks involve miRNAs. We developed an R/Bioconductor package, namely SpidermiR, which offers an easy access to both GRNs and miRNAs to the end user, and integrates this information with differentially expressed genes obtained from The Cancer Genome Atlas. Specifically, SpidermiR allows the users to: (i) query and download GRNs and miRNAs from validated and predicted repositories; (ii) integrate miRNAs with GRNs in order to obtain miRNA–gene–gene and miRNA–protein–protein interactions, and to analyze miRNA GRNs in order to identify miRNA–gene communities; and (iii) graphically visualize the results of the analyses. These analyses can be performed through a single interface and without the need for any downloads. The full data sets are then rapidly integrated and processed locally.

Cancer heterogeneity represents a major hurdle in the development of effective theranostic strategies, as it prevents to devise unique and maximally efficient diagnostic, prognostic and therapeutic procedures even for patients affected by the same tumor type. Computational techniques can nowadays leverage the huge and ever increasing amount of (epi)genomic data to tackle this problem, therefore providing new and valuable instruments for decision support to biologists and pathologists, in the broad sphere of precision medicine. In this context, we here introduce a novel cancer subtype classifier from gene expression data and we apply it to two different Breast Cancer datasets, from TCGA and GEO repositories. The classifier is based on Support Vector Machines and relies on the information about the relevant pathways involved in breast cancer development to reduce the huge variable space. Among the main results, we show that the classifier accuracy is preserved at excellent values even when the variable space is reduced by a 20-fold, hence providing a precious tool for cancer patient profiling even in case of limited experimental resources.

Background: An important challenge in cancer biology is to understand the complex aspects of the disease. It is increasingly evident that genes are not isolated from each other and the comprehension of how different genes are related to each other could explain biological mechanisms causing diseases. Biological pathways are important tools to reveal gene interaction and reduce the large number of genes to be studied by partitioning it into smaller paths. Furthermore, recent scientific evidence has proven that a combination of pathways, instead than a single element of the pathway or a single pathway, could be responsible for pathological changes in a cell. Results: In this paper we develop a new method that can reveal miRNAs able to regulate, in a coordinated way, networks of gene pathways. We applied the method to subtypes of breast cancer. The basic idea is the identification of pathways significantly enriched with differentially expressed genes among the different breast cancer subtypes and normal tissue. Looking at the pairs of pathways that were found to be functionally related, we created a network of dependent pathways and we focused on identifying miRNAs that could act as miRNA drivers in a coordinated regulation process. Conclusions: Our approach enables miRNAs identification that could have an important role in the development of breast cancer.

Understanding the dynamical evolution of cancer, with the final goal of developing effective techniques for diagnosis, prediction and treatment is one of the main challenges of modern biosciences. In this paper we approach the temporal ordering reconstruction problem, which refers to the temporal sorting of a collection of static biological data. The solution of this problem may help in better understanding the key principles and properties of the disease progression. By using a previously proposed technique for extracting temporal progressions from cross-sectional cancer gene expression data, we develop a novel methodology to be applied to static cross-sectional copy number alterations, and we test it on patients diagnosed with colorectal cancer at different stages. To capture distinct aspects of this complex phenomenon, we define several measures of chromosomal alterations and filters targeting significant portions of chromosomes. Results obtained with various measures and filters highlight the best setting for the problem, the most relevant chromosomal alterations and emphasise the influence that copy number alterations hitting key genes may have on the development of the disease.

Bioinformatics traditionally deals with computational approaches to the analysis of big data from high-throughput technologies as genomics, proteomics, and sequencing. Bioinformatics analysis allows extraction of new information from big data that might help to better assess the biological details at a molecular and cellular level. The wide-scale and high-dimensionality of Bioinformatics data has led to an increasing need of high performance computing and repository. In this chapter, the authors demonstrate the advantages of cloud computing in Bioinformatics research for high-throughput technologies.

Defining the aggressiveness and growth rate of a malignant cell population is a key step in the clinical approach to treating tumor disease. The correct grading of breast cancer (BC) is a fundamental part in determining the appropriate treatment. Biological variables can make it difficult to elucidate the mechanisms underlying BC development. To identify potential markers that can be used for BC classification, we analyzed mRNAs expression profiles, gene copy numbers, microRNAs expression and their association with tumor grade in BC microarray-derived datasets. From mRNA expression results, we found that grade 2 BC is most likely a mixture of grade 1 and grade 3 that have been misclassified, being described by the gene signature of either grade 1 or grade 3. We assessed the potential of the new approach of integrating mRNA expression profile, copy number alterations, and microRNA expression levels to select a limited number of genomic BC biomarkers. The combination of mRNA profile analysis and copy number data with microRNA expression levels led to the identification of two gene signatures of 42 and 4 altered genes (FOXM1, KPNA4, H2AFV and DDX19A) respectively, the latter obtained through a meta-analytical procedure. The 42-based gene signature identifies 4 classes of up- or down-regulated microRNAs (17 microRNAs) and of their 17 target mRNA, and the 4-based genes signature identified 4 microRNAs (Hsa-miR-320d, Hsa-miR-139-5p, Hsa-miR-567 and Hsa-let-7c). These results are discussed from a biological point of view with respect to pathological features of BC. Our identified mRNAs and microRNAs were validated as prognostic factors of BC disease progression, and could potentially facilitate the implementation of assays for laboratory validation, due to their reduced number.

Background: Copy number alterations (CNAs) represent an important component of genetic variations. Such alterations are related with certain type of cancer including those of the pancreas, colon, and breast, among others. CNAs have been used as biomarkers for cancer prognosis in multiple studies, but few works report on the relation of CNAs with the disease progression. Moreover, most studies do not consider the following two important issues. (I) The identification of CNAs in genes which are responsible for expression regulation is fundamental in order to define genetic events leading to malignant transformation and progression. (II) Most real domains are best described by structured data where instances of multiple types are related to each other in complex ways. Results: Our main interest is to check whether the colorectal cancer (CRC) progression inference benefits when considering both (I) the expression levels of genes with CNAs, and (II) relationships (i.e. dissimilarities) between patients due to expression level differences of the altered genes. We first evaluate the accuracy performance of a state-of-the-art inference method (support vector machine) when subjects are represented only through sets of available attribute values (i.e. gene expression level). Then we check whether the inference accuracy improves, when explicitly exploiting the information mentioned above. Our results suggest that the CRC progression inference improves when the combined data (i.e. CNA and expression level) and the considered dissimilarity measures are applied. Conclusions: Through our approach, classification is intuitively appealing and can be conveniently obtained in the resulting dissimilarity spaces. Different public datasets from Gene Expression Omnibus (GEO) were used to validate the results.

Specific genome copy number alterations, such as deletions and amplifications are an important factor in tumor development and progression, and are also associated with changes in gene expression. By combining analyses of gene expression and genome copy number we identified genes as candidate biomarkers of BC which were validated as prognostic factors of the disease progression. These results suggest that the proposed combined approach may become a valuable method for BC prognosis.

Tamoxifen is currently used for the treatment of breast cancer. Response to tamoxifen in metastatic conditions is a primary issue in cancer development. We used a cohort of breast cancer patients, treated or not with tamoxifen, and combined these data with the gene signature of metastatic samples in order to investigate the genetic mechanism of metastasis development, in search of a possible therapeutic effect of tamoxifen in metastatic conditions,. The analysis revealed a group of 21 genes common both to the set of up regulated genes in metastatic BC patients and to the set of down regulated genes in tamoxifen treated patients. These genes could be used as biomarkers for tamoxifen-sensitivity in order to optimize BC treatment

The development and progression of Colorectal Cancer (CRC) is a multi-step process leading to the accumulation of chromosomal instability (CIN) that occurs over the lifetime of a tumor (Shen et al, 2007; Vogelstein et al, 1988; Fearon et al, 1990). CINs include DNA copy number alterations (CNAs), i.e., regions of aberrantly increased or decreased DNA. For this reason, it is a challenge to identify both the regions that have CNAs and the genes whose expression could be deregulated because of gain or loss of their copies. In this paper we focus on the role of copy number alteration in assessing prognosis of patients with CRC. Specifically, we show that the inference of the CRC progression benefits from exploiting CNA information when a particular relational representation of patients is given. The employed framework outperforms standard approaches where patients are represented through a set of available attributes. The data set for this analysis was provided by IRCCS Istituto Nazionale dei Tumori Milano (INT) and deposited at NCBI Gene Expression Omnibus (GSE16125).

Cancer is a very complex disease and understanding its dynamics and evolution is one of the challenges of modern biosciences. As most available data on cancer is static, extracting dynamic information about its progression from “static” biological data would have a major significance. We are approaching the Temporal Ordering Reconstruction (TOR) problem, that is the sorting of a collection of multi-dimensional biological data to reflect an accurate temporal progression of the target disease. The most general form of the TOR problem has been studied from many points of view. Firstly, the TOR problem, as defined above has been tackled mostly in two works, which use gene expression data as the “raw’” data in the samples (Gupta and Bar-Joseph, 2008; Magwene et al., 2003). Secondly, another series of works start by analyzing comparative genomic hybridization data to build a plausible tree of possible gene mutation events and continue towards a use of Bayesian models to assess pathways variations in a disease (Desper et al., 1999; Pathare et al., 2009; Gerstung et al., 2011; Beerenwinkel et al., 2005). Our work is more focused on a specific approach to the TOR problem, previously proposed by Gupta and Bar-Joseph (Gupta and Bar-Joseph, 2008), which has been shown to work for gene expression data and we develop a methodology which enables us to apply this technique on a Copy Number Alterations (CNAs) data set. We also aim to provide a building block in an analysis pipeline that can be used to look at temporal reconstruction problems that assume an already (partially) ordered dataset (Ramakrishnan, 2010; Antoniotti, 2010).