Search

Biotechnological advances in sequencing have led to an explosion of publicly available data via large international consortia such as The Cancer Genome Atlas (TCGA), The Encyclopedia of DNA Elements (ENCODE), and The NIH Roadmap Epigenomics Mapping Consortium (Roadmap). These projects have provided unprecedented opportunities to interrogate the epigenome of cultured cancer cell lines as well as normal and tumor tissues with high genomic resolution. The bioconductor project offers more than 1,000 open-source software and statistical packages to analyze high-throughput genomic data. However, most packages are designed for specific data types (e.g. expression, epigenetics, genomics) and there is no comprehensive tool that provides a complete integrative analysis harnessing the resources and data provided by all three public projects. A need to create an integration of these different analyses was recently proposed. In this workflow, we provide a series of biologically focused integrative downstream analyses of different molecular data. We describe how to download, process and prepare TCGA data and by harnessing several key bioconductor packages, we describe how to extract biologically meaningful genomic and epigenomic data and by using Roadmap and ENCODE data, we provide a workplan to identify candidate biologically relevant functional epigenomic elements associated with cancer. To illustrate our workflow, we analyzed two types of brain tumors: low-grade glioma (LGG) versus high-grade glioma (glioblastoma multiform or GBM). This workflow introduces the following Bioconductor packages: AnnotationHub, ChIPSeeker, ComplexHeatmap, pathview, ELMER, GAIA, MINET, RTCGAtoolbox, TCGAbiolinks., SCOPUS: ar.j, info:eu-repo/semantics/published

The natural species tropism of hepatitis C virus (HCV) is limited to humans and higher primates. Due to the limited availability of chimpanzees, the establishment of a small-animal model for the study of HCV infection has a high priority. The recent identification of the human cellular receptors that are essential for HCV entry together with the advances in mouse genome engineering provides the bases for creating an inbred mouse model of HCV infection and propagation. In our study we designed a gene targeting strategy to produce novel humanized mouse strains expressing human-specific HCV entry receptors. Human genes for HCV essential entry factors were introduced in the mouse genome by homologous recombination in murine embryonic stem cells. To achieve the native expression patterns of their murine orthologues, we adopted a knock-in approach replacing the murine genes with human coding sequences. Three novel humanized mouse strains were generated, expressing respectively human CD81, scavenger receptor class B member 1 and occludin. Molecular analyses of mice genomic DNA and messenger RNA confirmed the presence and germ line propagation of the human genes, and their transcription. The native regulated expression of human-specific HCV entry receptors can overcome the major block to HCV replication in murine cells. To verify the virus entry into the cells of our humanized mice, we successfully optimized isolation and culture of murine primary hepatocytes and HCV pseudotype particles infection assay. Genetic humanization of mice for species-specific HCV receptors could allow HCV glycoproteins mediated cell entry and could be an important tool for in vivo studying of virus-host interactions.

Despite producing a panoply of potential cancer-specific targets, the proteogenomic characterization of human tumors has yet to demonstrate value for precision cancer medicine. Integrative multi-omics using a machine-learning network identified master kinases responsible for effecting phenotypic hallmarks of functional glioblastoma subtypes. In subtype-matched patient-derived models, we validated PKCδ and DNA-PK as master kinases of glycolytic/plurimetabolic and proliferative/progenitor subtypes, respectively, and qualified the kinases as potent and actionable glioblastoma subtype-specific therapeutic targets. Glioblastoma subtypes were associated with clinical and radiomics features, orthogonally validated by proteomics, phospho-proteomics, metabolomics, lipidomics and acetylomics analyses, and recapitulated in pediatric glioma, breast and lung squamous cell carcinoma, including subtype specificity of PKCδ and DNA-PK activity. We developed a probabilistic classification tool that performs optimally with RNA from frozen and paraffin-embedded tissues, which can be used to evaluate the association of therapeutic response with glioblastoma subtypes and to inform patient selection in prospective clinical trials.

Water reabsorption along the collecting duct is dependent on the function of aquaporin 2 (AQP2). Currently, information on microRNA (miRNA)-mediated, post-transcriptional regulation of AQP2, which may influence water reabsorption, is limited. In mice, ablation of the Dicer enzyme (crucial for miRNA maturation) in AQP2-expressing cells induces nephrogenic diabetes insipidus (NDI) with dysregulation of the miRNA profile. A major finding is the identification of miRNAs associated with NDI through mediating epigenetic control of AQP2. This study offers novel targets for AQP2 regulation and potential treatment for governing renal water reabsorption.Background MicroRNAs (miRNAs), formed by cleavage of pre-microRNA by the endoribonuclease Dicer, are critical modulators of cell function by post-transcriptionally regulating gene expression.Methods Selective ablation of Dicer in AQP2-expressing cells (DicerAQP2Cre+ mice) was used to investigate the role of miRNAs in the kidney collecting duct of mice.Results The mice had severe polyuria and nephrogenic diabetes insipidus, potentially due to greatly reduced AQP2 and AQP4 levels. Although epithelial sodium channel levels were decreased in cortex and increased in inner medulla, amiloride-sensitive sodium reabsorption was equivalent in DicerAQP2Cre+ mice and controls. Small-RNA sequencing and proteomic analysis revealed 31 and 178 significantly regulated miRNAs and proteins, respectively. Integrated bioinformatic analysis of the miRNAome and proteome suggested alterations in the epigenetic machinery and various transcription factors regulating AQP2 expression in DicerAQP2Cre+ mice. The expression profile and function of three miRNAs (miR-7688-5p, miR-8114, and miR-409-3p) whose predicted targets were involved in epigenetic control (Phf2, Kdm5c, and Kdm4a) or transcriptional regulation (GATA3, GATA2, and ELF3) of AQP2 were validated. Luciferase assays could not demonstrate direct interaction of AQP2 or the three potential transcription factors with miR-7688-5p, miR-8114, and miR-409textendash3p. However, transfection of respective miRNA mimics reduced AQP2 expression. Chromatin immunoprecipitation assays demonstrated decreased Phf2 and significantly increased Kdm5c interactions at the Aqp2 gene promoter in DicerAQP2Cre+ mice, resulting in decreased RNA Pol II association.Conclusions Novel evidence indicates miRNA-mediated epigenetic regulation of AQP2 expression.

Pancreatic organogenesis is a multistep process that requires the cooperation of several signaling pathways. In this context, the role of pancreatic mesenchyme is important to define the epithelium development, nevertheless, the precise space&ndash, temporal signaling activation still needs to be clarified. This study reports a dissection of the pancreatic embryogenesis, highlighting the molecular network surrounding the epithelium&ndash, mesenchyme interaction. To investigate this crosstalk, pancreatic epithelium and surrounding mesenchyme, at embryonic day 10.5, were collected through laser capture microdissection (LCM) and characterized based on their global gene expression. We performed a bioinformatic analysis to hypothesize crosstalk interactions, validating the most promising genes and verifying the precise localization of their expression in the compartments, by RNA in situ hybridization (ISH). Our analyses pointed out also the c-Met gene, a very well-known factor involved in stimulating motility, morphogenesis, and organ regeneration. We also highlighted the potential crosstalk between Versican (Vcan) and Syndecan4 (Sdc4) since these genes are involved in pancreatic tissue repair, strengthening the concept that the same signaling pathways required during pancreatic embryogenesis are also involved in tissue repair. This finding leads to novel strategies for obtaining functional pancreatic stem cells for cell replacement therapies.

The progressive and physiological decline in ovarian function depends on the rate of follicular loss by atresia, contributing to the reduction in ovarian reserve. Genetics and environmental factors play important roles in ovarian senescence and in the onset of ovarian dysfunctions such as diminished ovarian reserve. A better understanding of the mechanisms underlying ovarian aging and their regulation by genetic and environmental factors is needed to evaluate ovarian reserve and to predict fertility potential by identification of more accurate and less invasive markers. We report transcriptomic data (i) implicating novel (e.g. EIF2 signalling) and well-known pathways (e.g. TGFβ signalling), and (ii) defining a unique set of non-coding RNA (ncRNA), both associated with ovarian function. The latter includes miRNAs (e.g. Mir143 and Mir145), snoRNAs (e.g. Snord16a and Snora34), and one lncRNA (Gas5), which are differentially expressed in middle-aged ovaries (12 months) vs young-aged (3 months) from CD1 mice. Experimental analysis confirms that ovary lifespan varies across genetic backgrounds in mice and, genetics influences the response to environmental perturbations such as diet. Moreover, the identified ncRNAs were verified in a model of reproductive dysfunction promoted by the environmental toxicant ethylenthiourea. We also report the increase of miRNA143 and miRNA145 in follicular fluid of women with diminished ovarian reserve. Their levels inversely correlate with the hormonal profile and with the number of the oocytes recruited upon hormonal stimulation. Overall, we report a transcriptomic signature for ovarian dysfunction in vivo that provides a valuable resource for translational research in human reproductive aging.

Neurofibromatosis type 1 (NF1) is a common tumor predisposition syndrome in which glioma is one of the prevalent tumors. Gliomagenesis in NF1 results in a heterogeneous spectrum of low- to high-grade neoplasms occurring during the entire lifespan of patients. The pattern of genetic and epigenetic alterations of glioma that develops in NF1 patients and the similarities with sporadic glioma remain unknown. Here, we present the molecular landscape of low- and high-grade gliomas in patients affected by NF1 (NF1-glioma). We found that the predisposing germline mutation of the NF1 gene was frequently converted to homozygosity and the somatic mutational load of NF1-glioma was influenced by age and grade. High-grade tumors harbored genetic alterations of TP53 and CDKN2A, frequent mutations of ATRX associated with Alternative Lengthening of Telomere, and were enriched in genetic alterations of transcription/chromatin regulation and PI3 kinase pathways. Low-grade tumors exhibited fewer mutations that were over-represented in genes of the MAP kinase pathway. Approximately 50% of low-grade NF1-gliomas displayed an immune signature, T lymphocyte infiltrates, and increased neo-antigen load. DNA methylation assigned NF1-glioma to LGm6, a poorly defined Isocitrate Dehydrogenase 1 wild-type subgroup enriched with ATRX mutations. Thus, the profiling of NF1-glioma defined a distinct landscape that recapitulates a subset of sporadic tumors. An integrated analysis of glioma samples from patients with neurofibromatosis 1 annotates their mutational, epigenetic, transcriptional, and immunological features and uncovers similitudes with a subset of sporadic gliomas.

In vitro Omics analysis (i.e. transcriptome) is suggested to predict in vivo toxicity and adverse effects in humans, although the causal link between high-throughput data and effects in vivo is not easily established. Indeed, the chemical-organism interaction can involve processes, such as adaptation, not established in cell cultures. Starting from this consideration we investigate the transcriptomic response of immortalized thyrocytes to ethylenthiourea and chlorpyrifos. In vitro data revealed specific and common genes/mechanisms of toxicity, controlling the proliferation/survival of the thyrocytes and unrelated hematopoietic cell lineages. These results were phenotypically confirmed in vivo by the reduction of circulating T4 hormone and the development of pancytopenia after long exposure. Our data imply that in vitro toxicogenomics is a powerful tool in predicting adverse effects in vivo, experimentally confirming the vision described as Tox21c (Toxicity Testing in the 21st century) although not fully recapitulating the biocomplexity of a living animal.

Epidemiological and experimental data highlighted the thyroid-disrupting activity of bisphenol A (BPA). Although pivotal to identify the mechanisms of toxicity, direct low-dose BPA effects on thyrocytes have not been assessed. Here, we report the results of microarray experiments revealing that the transcriptome reacts dynamically to low-dose BPA exposure, adapting the changes in gene expression to the exposure duration. The response involves many genes, enriching specific pathways and biological functions mainly cell death/proliferation or DNA repair. Their expression is only slightly altered but, since they enrich specific pathways, this results in major effects as shown here for transcripts involved in the DNA repair pathway. Indeed, even though no phenotypic changes are induced by the treatment, we show that the exposure to BPA impairs the cell response to further stressors. We experimentally verify that prolonged exposure to low doses of BPA results in a delayed response to UV-C-induced DNA damage, due to impairment of p21-Tp53 axis, with the BPA-treated cells more prone to cell death and DNA damage accumulation. The present findings shed light on a possible mechanism by which BPA, not able to directly cause genetic damage at environmental dose, may exert an indirect genotoxic activity.

Background Comparison of toxicogenomic data facilitates the identification of deregulated gene patterns and maximizes health risk prediction in human. Results Here, we performed phenotypic anchoring on the effects of acute exposure to low-grade polluted groundwater using mouse and zebrafish. Also, we evaluated two windows of chronic exposure in mouse, starting in utero and at the end of lactation. Bioinformatic analysis of livers microarray data showed that the number of deregulated biofunctions and pathways is higher after acute exposure, compared to the chronic one. It also revealed specific profiles of altered gene expression in all treatments, pointing to stress response/mitochondrial pathways as major players of environmental toxicity. Of note, dysfunction of steroid hormones was also predicted by bioinformatic analysis and verified in both models by traditional approaches, serum estrogens measurement and vitellogenin mRNA determination in mice and zebrafish, respectively. Conclusions In our report, phenotypic anchoring in two vertebrate model organisms highlights the toxicity of low-grade pollution, with varying susceptibility based on exposure window. The overlay of zebrafish and mice deregulated pathways, more than single genes, is useful in risk identification from chemicals implicated in the observed effects. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-1067) contains supplementary material, which is available to authorized users.