DNMT inhibitors reverse a specific signature of aberrant promoter DNA methylation and associated gene silencing in AML

<b>Background</b>. Myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) are neoplastic disorders of hematopoietic stem cells. DNA methyltransferase inhibitors (DNMTi), 5-azacytidine (AzaC) and 5-aza-2’-deoxycytidine (Decitabine), benefit some MDS/AML patients. However, the role of DNMTi-induced DNA hypomethylation in regulation of gene expression in AML is unclear.<p></p> <b>Results. </b> We compared the effects of AzaC on DNA methylation and gene expression using whole-genome single-nucleotide bisulfite-sequencing (WGBS) and RNA-sequencing in OCI-AML3 (AML3) cells. For data analysis, we used an approach recently developed for discovery of differential patterns of DNA methylation associated with changes in gene expression, that is tailored to single-nucleotide bisulfite-sequencing data (Washington University Interpolated Methylation Signatures (WIMSi)). By this approach, a subset of genes upregulated by AzaC was found to be characterized by AzaC-induced signature methylation loss flanking the transcription start site. These genes are enriched for genes increased in methylation and decreased in expression in AML3 cells compared to normal hematopoietic stem and progenitor cells. Moreover, these genes are preferentially upregulated by Decitabine in human primary AML blasts, and control cell proliferation, death and development. <p></p> <b>Conclusions.</b> Our WGBS and WIMSi data analysis approach has identified a set of genes whose is methylation and silencing in AML is reversed by DNMTi. These genes are good candidates for direct regulation by DNMTi, and their reactivation by DNMTi may contribute to therapeutic activity. This study also demonstrates the ability of WIMSi to reveal relationships between DNA methylation and gene expression, based on single-nucleotide bisulfite-sequencing and RNA-seq data.<p></p&gt

MicroRNA-23a promotes myelination in the central nervous system.

Demyelinating disorders including leukodystrophies are devastating conditions that are still in need of better understanding, and both oligodendrocyte differentiation and myelin synthesis pathways are potential avenues for developing treatment. Overexpression of lamin B1 leads to leukodystrophy characterized by demyelination of the central nervous system, and microRNA-23 (miR-23) was found to suppress lamin B1 and enhance oligodendrocyte differentiation in vitro. Here, we demonstrated that miR-23a-overexpressing mice have increased myelin thickness, providing in vivo evidence that miR-23a enhances both oligodendrocyte differentiation and myelin synthesis. Using this mouse model, we explored possible miR-23a targets and revealed that the phosphatase and tensin homologue/phosphatidylinositol trisphosphate kinase/Akt/mammalian target of rapamycin pathway is modulated by miR-23a. Additionally, a long noncoding RNA, 2700046G09Rik, was identified as a miR-23a target and modulates phosphatase and tensin homologue itself in a miR-23a-dependent manner. The data presented here imply a unique role for miR-23a in the coordination of proteins and noncoding RNAs in generating and maintaining healthy myelin

Cost effective, experimentally robust differential-expression analysis for human/mammalian, pathogen and dual-species transcriptomics.

As sequencing read length has increased, researchers have quickly adopted longer reads for their experiments. Here, we examine 14 pathogen or host-pathogen differential gene expression data sets to assess whether using longer reads is warranted. A variety of data sets was used to assess what genomic attributes might affect the outcome of differential gene expression analysis including: gene density, operons, gene length, number of introns/exons and intron length. No genome attribute was found to influence the data in principal components analysis, hierarchical clustering with bootstrap support, or regression analyses of pairwise comparisons that were undertaken on the same reads, looking at all combinations of paired and unpaired reads trimmed to 36, 54, 72 and 101 bp. Read pairing had the greatest effect when there was little variation in the samples from different conditions or in their replicates (e.g. little differential gene expression). But overall, 54 and 72 bp reads were typically most similar. Given differences in costs and mapping percentages, we recommend 54 bp reads for organisms with no or few introns and 72 bp reads for all others. In a third of the data sets, read pairing had absolutely no effect, despite paired reads having twice as much data. Therefore, single-end reads seem robust for differential-expression analyses, but in eukaryotes paired-end reads are likely desired to analyse splice variants and should be preferred for data sets that are acquired with the intent to be community resources that might be used in secondary data analyses

The tapeworm interactome: inferring confidence scored protein-protein interactions from the proteome of Hymenolepis microstoma

BACKGROUND: Reference genome and transcriptome assemblies of helminths have reached a level of completion whereby secondary analyses that rely on accurate gene estimation or syntenic relationships can be now conducted with a high level of confidence. Recent public release of the v.3 assembly of the mouse bile-duct tapeworm, Hymenolepis microstoma, provides chromosome-level characterisation of the genome and a stabilised set of protein coding gene models underpinned by bioinformatic and empirical data. However, interactome data have not been produced. Conserved protein-protein interactions in other organisms, termed interologs, can be used to transfer interactions between species, allowing systems-level analysis in non-model organisms. RESULTS: Here, we describe a probabilistic, integrated network of interologs for the H. microstoma proteome, based on conserved protein interactions found in eukaryote model species. Almost a third of the 10,139 gene models in the v.3 assembly could be assigned interaction data and assessment of the resulting network indicates that topologically-important proteins are related to essential cellular pathways, and that the network clusters into biologically meaningful components. Moreover, network parameters are similar to those of single-species interaction networks that we constructed in the same way for S. cerevisiae, C. elegans and H. sapiens, demonstrating that information-rich, system-level analyses can be conducted even on species separated by a large phylogenetic distance from the major model organisms from which most protein interaction evidence is based. Using the interolog network, we then focused on sub-networks of interactions assigned to discrete suites of genes of interest, including signalling components and transcription factors, germline multipotency genes, and genes differentially-expressed between larval and adult worms. Results show not only an expected bias toward highly-conserved proteins, such as components of intracellular signal transduction, but in some cases predicted interactions with transcription factors that aid in identifying their target genes. CONCLUSIONS: With key helminth genomes now complete, systems-level analyses can provide an important predictive framework to guide basic and applied research on helminths and will become increasingly informative as new protein-protein interaction data accumulate

Single-cell epigenomic variability reveals functional cancer heterogeneity.

BackgroundCell-to-cell heterogeneity is a major driver of cancer evolution, progression, and emergence of drug resistance. Epigenomic variation at the single-cell level can rapidly create cancer heterogeneity but is difficult to detect and assess functionally.ResultsWe develop a strategy to bridge the gap between measurement and function in single-cell epigenomics. Using single-cell chromatin accessibility and RNA-seq data in K562 leukemic cells, we identify the cell surface marker CD24 as co-varying with chromatin accessibility changes linked to GATA transcription factors in single cells. Fluorescence-activated cell sorting of CD24 high versus low cells prospectively isolated GATA1 and GATA2 high versus low cells. GATA high versus low cells express differential gene regulatory networks, differential sensitivity to the drug imatinib mesylate, and differential self-renewal capacity. Lineage tracing experiments show that GATA/CD24hi cells have the capability to rapidly reconstitute the heterogeneity within the entire starting population, suggesting that GATA expression levels drive a phenotypically relevant source of epigenomic plasticity.ConclusionSingle-cell chromatin accessibility can guide prospective characterization of cancer heterogeneity. Epigenomic subpopulations in cancer impact drug sensitivity and the clonal dynamics of cancer evolution

A Scalable Approach for Discovering Conserved Active Subnetworks across Species

Overlaying differential changes in gene expression on protein interaction networks has proven to be a useful approach to interpreting the cell's dynamic response to a changing environment. Despite successes in finding active subnetworks in the context of a single species, the idea of overlaying lists of differentially expressed genes on networks has not yet been extended to support the analysis of multiple species' interaction networks. To address this problem, we designed a scalable, cross-species network search algorithm, neXus (Network - cross(X)-species - Search), that discovers conserved, active subnetworks based on parallel differential expression studies in multiple species. Our approach leverages functional linkage networks, which provide more comprehensive coverage of functional relationships than physical interaction networks by combining heterogeneous types of genomic data. We applied our cross-species approach to identify conserved modules that are differentially active in stem cells relative to differentiated cells based on parallel gene expression studies and functional linkage networks from mouse and human. We find hundreds of conserved active subnetworks enriched for stem cell-associated functions such as cell cycle, DNA repair, and chromatin modification processes. Using a variation of this approach, we also find a number of species-specific networks, which likely reflect mechanisms of stem cell function that have diverged between mouse and human. We assess the statistical significance of the subnetworks by comparing them with subnetworks discovered on random permutations of the differential expression data. We also describe several case examples that illustrate the utility of comparative analysis of active subnetworks

Rare and common epilepsies converge on a shared gene regulatory network providing opportunities for novel antiepileptic drug discovery

Background The relationship between monogenic and polygenic forms of epilepsy is poorly understood, and the extent to which the genetic and acquired epilepsies share common pathways is unclear. Here, we use an integrated systems-level analysis of brain gene expression data to identify molecular networks disrupted in epilepsy. Results We identify a co-expression network of 320 genes (M30), which is significantly enriched for non-synonymous de novo mutations ascertained from patients with monogenic epilepsy, and for common variants associated with polygenic epilepsy. The genes in M30 network are expressed widely in the human brain under tight developmental control, and encode physically interacting proteins involved in synaptic processes. The most highly connected proteins within M30 network are preferentially disrupted by deleterious de novo mutations for monogenic epilepsy, in line with the centrality-lethality hypothesis. Analysis of M30 expression revealed consistent down-regulation in the epileptic brain in heterogeneous forms of epilepsy including human temporal lobe epilepsy, a mouse model of acquired temporal lobe epilepsy, and a mouse model of monogenic Dravet (SCN1A) disease. These results suggest functional disruption of M30 via gene mutation or altered expression as a convergent mechanism regulating susceptibility to epilepsy broadly. Using the large collection of drug-induced gene expression data from Connectivity Map, several drugs were predicted to preferentially restore the down-regulation of M30 in epilepsy toward health, most notably valproic acid, whose effect on M30 expression was replicated in neurons. Conclusions Taken together, our results suggest targeting the expression of M30 as a potential new therapeutic strategy in epilepsy

Lineage-specific interface proteins match up the cell cycle and differentiation in embryo stem cells.

The shortage of molecular information on cell cycle changes along embryonic stem cell (ESC) differentiation prompts an in silico approach, which may provide a novel way to identify candidate genes or mechanisms acting in coordinating the two programs. We analyzed germ layer specific gene expression changes during the cell cycle and ESC differentiation by combining four human cell cycle transcriptome profiles with thirteen in vitro human ESC differentiation studies. To detect cross-talk mechanisms we then integrated the transcriptome data that displayed differential regulation with protein interaction data. A new class of non-transcriptionally regulated genes was identified, encoding proteins which interact systematically with proteins corresponding to genes regulated during the cell cycle or cell differentiation, and which therefore can be seen as interface proteins coordinating the two programs. Functional analysis gathered insights in fate-specific candidates of interface functionalities. The non-transcriptionally regulated interface proteins were found to be highly regulated by post-translational ubiquitylation modification, which may synchronize the transition between cell proliferation and differentiation in ESCs

Long noncoding RNAs (lncRNAs) dynamics evidence immunomodulation during ISAV-Infected Atlantic salmon (Salmo salar)

Despite evidence for participation in the host response to infection, the roles of many long non-coding RNAs (lncRNAs) remain unknown. Therefore, the aims of this study were to identify lncRNAs in Atlantic salmon (Salmo salar) and evaluate their transcriptomic regulation during ISA virus (ISAV) infection, an Orthomyxoviridae virus associated with high mortalities in salmonid aquaculture. Using next-generation sequencing, whole-transcriptome analysis of theSalmo salarresponse to ISAV infection was performed, identifying 5,636 putative lncRNAs with a mean length of 695 base pairs. The transcriptional modulation evidenced a similar number of differentially expressed lncRNAs in the gills (3,294), head-kidney (3,275), and liver (3,325) over the course of the infection. Moreover, analysis of a subset of these lncRNAs showed the following: (i) Most were similarly regulated in response to ISA virus infection; (ii) The transcript subsets were uniquely modulated in each tissue (gills, liver, and head-kidney); and (iii) A subset of lncRNAs were upregulated for each tissue and time analysed, indicating potential markers for ISAV infection. These findings represent the first discovery of widespread differential expression of lncRNAs in response to virus infection in non-model species, suggesting that lncRNAs could be involved in regulating the host response during ISAV infection