Pscan: finding over-represented transcription factor binding site motifs in sequences from co-regulated or co-expressed genes

CREB Regulates Distinct Adaptive Transcriptional Programs in Astrocytes and Neurons

The cyclic AMP response element binding protein (CREB) is a primary hub of activity-driven genetic programs in neurons controlling plasticity, neurogenesis and survival. By contrast, the gene networks coordinated by CREB in astrocytes are unknown despite the fact that the astrocytic CREB is also activity-driven and neuroprotective. Herein we identified the transcriptional programs regulated by CREB in astrocytes as compared to neurons using, as study materials, transcriptome databases of astrocyte exposed to well-known activators of CREB-dependent transcription as well as publicly available transcriptomes of neuronal cultures. Functional CREB signatures were extracted from the transcriptomes using Gene Ontology, adult-brain gene lists generated by Translating Ribosome Affinity Purification (TRAP) and CREB-target gene repositories. We found minimal overlap between CREB signatures in astrocytes and neurons. In astrocytes, the top triad of functions regulated by CREB consists of 'Gene expression', 'Mitochondria', and 'Signalling', while in neurons it is 'Neurotransmission', 'Signalling' and 'Gene expression', the latter two being represented by different genes from those in astrocytes. The newly generated databases will provide a tool to explore novel means whereby CREB impinges on brain functions requiring adaptive, long-lasting changes by coordinating transcriptional cascades in astrocytes.

SRF modulates seizure occurrence, activity induced gene transcription and hippocampal circuit reorganization in the mouse pilocarpine epilepsy model

A hallmark of temporal lobe epilepsy (TLE) is hippocampal neuronal demise and aberrant mossy fiber sprouting. In addition, unrestrained neuronal activity in TLE patients induces gene expression including immediate early genes (IEGs) such as Fos and Egr1.

We employed the mouse pilocarpine model to analyze the transcription factor (TF) serum response factor (SRF) in epileptogenesis, seizure induced histopathology and IEG induction. SRF is a neuronal activity regulated TF stimulating IEG expression as well as nerve fiber growth and guidance. Adult conditional SRF deficient mice (Srf^CaMKCreERT2) were more refractory to initial status epilepticus (SE) acquisition. Further, SRF deficient mice developed more spontaneous recurrent seizures (SRS). Genome-wide transcriptomic analysis uncovered a requirement of SRF for SE and SRS induced IEG induction (e.g. Fos, Egr1, Arc, Npas4, Btg2, Atf3). SRF was required for epilepsy associated neurodegeneration, mossy fiber sprouting and inflammation. We uncovered MAP kinase signaling as SRF target during epilepsy. Upon SRF ablation, seizure evoked induction of dual specific phosphatases (Dusp5 and Dusp6) was reduced. Lower expression of these negative ERK kinase regulators correlated with altered P-ERK levels in epileptic Srf mutant animals.

Overall, this study uncovered an SRF contribution to several processes of epileptogenesis in the pilocarpine model.

Single-cell profiling reveals heterogeneity and functional patterning of GPCR expression in the vascular system

G-protein-coupled receptor (GPCR) expression is extensively studied in bulk cDNA, but heterogeneity and functional patterning of GPCR expression in individual vascular cells is poorly understood. Here, we perform a microfluidic-based single-cell GPCR expression analysis in primary smooth muscle cells (SMC) and endothelial cells (EC). GPCR expression is highly heterogeneous in all cell types, which is confirmed in reporter mice, on the protein level and in human cells. Inflammatory activation in murine models of sepsis or atherosclerosis results in characteristic changes in the GPCR repertoire, and we identify functionally relevant subgroups of cells that are characterized by specific GPCR patterns. We further show that dedifferentiating SMC upregulate GPCRs such as Gpr39, Gprc5b, Gprc5c or Gpr124, and that selective targeting of Gprc5b modulates their differentiation state. Taken together, single-cell profiling identifies receptors expressed on pathologically relevant subpopulations and provides a basis for the development of new therapeutic strategies in vascular diseases.

GPCRs are key regulators of vascular functions. By analysing single-cell GPCRs expression in vascular smooth muscle and endothelial cells from healthy and diseased murine vessels, Kaur et al. show that GPCR expression is highly heterogeneous in all cell types and that disease causes GPCR repertoire changes depending on cell type and vascular localization.

DNA Methylation and Hydroxymethylation Profile of CD34+-Enriched Cell Products Intended for Autologous CD34+ Cell Transplantation

Epigenetic dysregulation has been shown to limit functional capacity of aging hematopoietic stem cells, which may contribute to impaired outcome of hematopoietic stem cell-based therapies. The aim of our study was to gain better insight into the epigenetic profile of CD34⁺-enriched cell products intended for autologous CD34⁺ cell transplantation in patients with cardiomyopathy. We found global DNA methylation content significantly higher in immunoselected CD34⁺ cells compared to leukocytes in leukapheresis products (2.33 ± 1.03% vs. 1.84 ± 0.86%, p = 0.04). Global DNA hydroxymethylation content did not differ between CD34⁺ cells and leukocytes (p = 0.30). By measuring methylation levels of 94 stem cell transcription factors on a ready-to-use array, we identified 15 factors in which average promoter methylation was significantly different between leukocytes and CD34⁺ cells. The difference was highest for HOXC12 (58.18 ± 6.47% vs. 13.34 ± 24.18%, p = 0.0009) and NR2F2 (51.65 ± 25.89% vs. 7.66 ± 21.43%, p = 0.0045) genes. Our findings suggest that global DNA methylation and hydroxymethylation patterns as well as target methylation profile of selected genes in CD34⁺-enriched cell products do not differ significantly compared to leukapheresis products and, thus, can tell us little about the functional capacity and regenerative properties of CD34⁺ cells. Future studies should examine other CD34⁺ cell graft characteristics, which may serve as prognostic tools for autologous CD34⁺ cell transplantation.

CONSTANS Imparts DNA Sequence Specificity to the Histone Fold NF-YB/NF-YC Dimer

Nuclear Factor Y (NF-Y) is a heterotrimeric transcription factor that binds CCAAT elements. The NF-Y trimer is composed of a Histone Fold Domain (HFD) dimer (NF-YB/NF-YC) and NF-YA, which confers DNA sequence specificity. NF-YA shares a conserved domain with the CONSTANS, CONSTANS-LIKE, TOC1 (CCT) proteins. We show that CONSTANS (CO/B-BOX PROTEIN1 BBX1), a master flowering regulator, forms a trimer with Arabidopsis thaliana NF-YB2/NF-YC3 to efficiently bind the CORE element of the FLOWERING LOCUS T promoter. We term this complex NF-CO. Using saturation mutagenesis, electrophoretic mobility shift assays, and RNA-sequencing profiling of co, nf-yb, and nf-yc mutants, we identify CCACA elements as the core NF-CO binding site. CO physically interacts with the same HFD surface required for NF-YA association, as determined by mutations in NF-YB2 and NF-YC9, and tested in vitro and in vivo. The co-7 mutation in the CCT domain, corresponding to an NF-YA arginine directly involved in CCAAT recognition, abolishes NF-CO binding to DNA. In summary, a unifying molecular mechanism of CO function relates it to the NF-YA paradigm, as part of a trimeric complex imparting sequence specificity to HFD/DNA interactions. It is likely that members of the large CCT family participate in similar complexes with At-NF-YB and At-NF-YC, broadening HFD combinatorial possibilities in terms of trimerization, DNA binding specificities, and transcriptional regulation.

Sequential EMT-MET induces neuronal conversion through Sox2

Direct neuronal conversion can be achieved with combinations of small-molecule compounds and growth factors. Here, by studying the first or induction phase of the neuronal conversion induced by defined 5C medium, we show that the Sox2-mediated switch from early epithelial-mesenchymal transition (EMT) to late mesenchymal-epithelial transition (MET) within a high proliferation context is essential and sufficient for the conversion from mouse embryonic fibroblasts (MEFs) to TuJ⁺ cells. At the early stage, insulin and basic fibroblast growth factor (bFGF)-induced cell proliferation, early EMT, the up-regulation of Stat3 and Sox2, and the subsequent activation of neuron projection. Up-regulated Sox2 then induced MET and directed cells towards a neuronal fate at the late stage. Inhibiting either stage of this sequential EMT-MET impaired the conversion. In addition, Sox2 could replace sequential EMT-MET to induce a similar conversion within a high proliferation context, and its functions were confirmed with other neuronal conversion protocols and MEFs reprogramming. Therefore, the critical roles of the sequential EMT-MET were implicated in direct cell fate conversion in addition to reprogramming, embryonic development and cancer progression.

Wounding induces dedifferentiation of epidermal Gata6+ cells and acquisition of stem cell properties

The epidermis is maintained by multiple stem cell populations whose progeny differentiate along diverse, and spatially distinct, lineages. Here we show that the transcription factor Gata6 controls the identity of the previously uncharacterized sebaceous duct (SD) lineage and identify the Gata6 downstream transcription factor network that specifies a lineage switch between sebocytes and SD cells. During wound healing differentiated Gata6⁺ cells migrate from the SD into the interfollicular epidermis and dedifferentiate, acquiring the ability to undergo long-term self-renewal and differentiate into a much wider range of epidermal lineages than in undamaged tissue. Our data not only demonstrate that the structural and functional complexity of the junctional zone is regulated by Gata6, but also reveal that dedifferentiation is a previously unrecognized property of post-mitotic, terminally differentiated cells that have lost contact with the basement membrane. This resolves the long-standing debate about the contribution of terminally differentiated cells to epidermal wound repair.

Selective inhibition of EZH2 by a small molecule inhibitor regulates microglial gene expression essential for inflammation

Multiple studies have documented that Enhancer of zeste homolog 2 (EZH2) could play a role in inflammation and a wide range of malignancies; however, the underlying mechanisms remain largely unaddressed. Microglial activation is a key process in the production and release of numerous pro-inflammatory mediators that play important roles in inflammation and neurodegeneration in the central nervous system (CNS). Therefore, our aim was to investigate whether inhibition of EZH2 with the selective small molecule inhibitor EPZ-6438 protects against neonatal microglial activation. First, in mouse primary microglial cells and a microglial cell line, we found that LPS can rapidly increase EZH2 mRNA level and we subsequently performed gene expression profiling and constructed networks in resting, EPZ-6438-treated, LPS-treated and LPS+EPZ-6438-treated primary microglial cells and a microglial cell line using transcriptome RNA sequencing and bioinformatics analyses. By examining the RNA sequencing, we identified EPZ-6438 target genes and co-regulated modules that were critical for inflammation. We also identified unexpected relationships between the inducible transcription factors (TFs), motif strength, and the transcription of key inflammatory mediators. Furthermore, we showed that EPZ-6438 controls important inflammatory gene targets by modulating interferon regulatory factor (IRF) 1, IRF8, and signal transducer and activator of transcription (STAT) 1 levels at their promoter sites. Our unprecedented findings demonstrate that pharmacological interventions built upon EZH2 inhibition by EPZ-6438 could be a useful therapeutic approach for the treatment of neuroinflammatory diseases associated with microglial activation.

FTY720 (fingolimod) regulates key target genes essential for inflammation in microglial cells as defined by high-resolution mRNA sequencing

Although microglial cells have an essential role in the host defense of the brain, the abnormal activation of microglia can lead to devastating outcomes, such as neuroinflammation and neurodegeneration. Emerging evidence indicates that FTY720 (fingolimod), an FDA-approved drug, has beneficial effects on brain cells in the central nervous system (CNS) and, more recently, immunosuppressive activities in microglia via modulation of the sphingosine 1 phosphate (S1P) 1 receptor. However, the exact molecular aspects of FTY720 contribution in microglia remain largely unaddressed. To understand the molecular mechanisms underlying the roles of FTY720 in microglia, we performed gene expression profiling in resting, FTY720, LPS and LPS + FTY720 challenged primary microglial (PM) cells isolated from 3-day-old ICR mice, and we identified FTY720 target genes and co-regulated modules that were critical in inflammation. By examining RNA sequencing and binding motif datasets from FTY720 suppressed LPS-induced inflammatory mediators, we also identified unexpected relationships between the inducible transcription factors (TFs), motif strength, and the transcription of key inflammatory mediators. Furthermore, we showed that FTY720 controls important inflammatory genes targets by modulating STAT1 and IRF8 levels at their promoter site. Our unprecedented findings demonstrate that FTY720 could be a useful therapeutic application for neuroinflammatory diseases associated with microglia activation, as well as provide a rich resource and framework for future analyses of FTY720 effects on microglia interaction.

SWIM: a computational tool to unveiling crucial nodes in complex biological networks

SWItchMiner (SWIM) is a wizard-like software implementation of a procedure, previously described, able to extract information contained in complex networks. Specifically, SWIM allows unearthing the existence of a new class of hubs, called “fight-club hubs”, characterized by a marked negative correlation with their first nearest neighbors. Among them, a special subset of genes, called “switch genes”, appears to be characterized by an unusual pattern of intra- and inter-module connections that confers them a crucial topological role, interestingly mirrored by the evidence of their clinic-biological relevance. Here, we applied SWIM to a large panel of cancer datasets from The Cancer Genome Atlas, in order to highlight switch genes that could be critically associated with the drastic changes in the physiological state of cells or tissues induced by the cancer development. We discovered that switch genes are found in all cancers we studied and they encompass protein coding genes and non-coding RNAs, recovering many known key cancer players but also many new potential biomarkers not yet characterized in cancer context. Furthermore, SWIM is amenable to detect switch genes in different organisms and cell conditions, with the potential to uncover important players in biologically relevant scenarios, including but not limited to human cancer.

Abnormal epigenetic changes during differentiation of human skeletal muscle stem cells from obese subjects

BACKGROUND

Human skeletal muscle stem cells are important for muscle regeneration. However, the combined genome-wide DNA methylation and expression changes taking place during adult myogenesis have not been described in detail and novel myogenic factors may be discovered. Additionally, obesity is associated with low relative muscle mass and diminished metabolism. Epigenetic alterations taking place during myogenesis might contribute to these defects.

METHODS

We used Infinium HumanMethylation450 BeadChip Kit (Illumina) and HumanHT-12 Expression BeadChip (Illumina) to analyze genome-wide DNA methylation and transcription before versus after differentiation of primary human myoblasts from 14 non-obese and 14 obese individuals. Functional follow-up experiments were performed using siRNA mediated gene silencing in primary human myoblasts and a transgenic mouse model.

RESULTS

We observed genome-wide changes in DNA methylation and expression patterns during differentiation of primary human muscle stem cells (myoblasts). We identified epigenetic and transcriptional changes of myogenic transcription factors (MYOD1, MYOG, MYF5, MYF6, PAX7, MEF2A, MEF2C, and MEF2D), cell cycle regulators, metabolic enzymes and genes previously not linked to myogenesis, including IL32, metallothioneins, and pregnancy-specific beta-1-glycoproteins. Functional studies demonstrated IL-32 as a novel target that regulates human myogenesis, insulin sensitivity and ATP levels in muscle cells. Furthermore, IL32 transgenic mice had reduced insulin response and muscle weight. Remarkably, approximately 3.7 times more methylation changes (147,161 versus 39,572) were observed during differentiation of myoblasts from obese versus non-obese subjects. In accordance, DNMT1 expression increased during myogenesis only in obese subjects. Interestingly, numerous genes implicated in metabolic diseases and epigenetic regulation showed differential methylation and expression during differentiation only in obese subjects.

CONCLUSIONS

Our study identifies IL-32 as a novel myogenic regulator, provides a comprehensive map of the dynamic epigenome during differentiation of human muscle stem cells and reveals abnormal epigenetic changes in obesity.

Major intrinsic proteins repertoire of Morus notabilis and their expression profiles in different species

Leaf moisture content in Morus is a significant trait regulating the yield of silk production. Studies have shown that fresh leaves or leaves with high water content are preferably eaten by silk worm. Water and certain other molecules transport in plants is known to be regulated by aquaporins or Major Intrinsic Proteins (MIPs). Members of the MIP gene family have also been implicated in plant development and stress responsiveness. To understand how members of MIP gene family are regulated and evolved, we carried out an extensive analysis of the gene family. We identified a total of 36 non redundant MIPs in Morus notabilis genome, belonging to five subfamilies PIPs, TIPs, NIPs, XIPs and SIPs) have been identified. We performed a Gene ontology (GO) term enrichment analysis and looked at distribution of cis elements in their 2K upstream regulatory region to reveal their putative roles in various stresses and developmental aspects. Expression analysis in developmental stages revealed their tissue preferential expression pattern in diverse vegetative and reproductive tissues. Comparison of expression profiles in the leaves of three species including Morus notabilis, Morus serrata and Morus laevigata led to identification of differential expression in these species. In all, this study elaborates a basic insight into the structure, function and evolutionary analysis of MIP gene family in Morus which is hitherto unavailable. Our analysis will provide a ready reference to the mulberry research community involved in the Morus improvement program.

Decoding transcriptional states in cancer

Gene regulatory networks determine cellular identity. In cancer, aberrations of gene networks are caused by driver mutations that often affect transcription factors and chromatin modifiers. Nevertheless, gene transcription in cancer follows the same cis-regulatory rules as normal cells, and cancer cells have served as convenient model systems to study transcriptional regulation. Tumours often show regulatory heterogeneity, with subpopulations of cells in different transcriptional states, which has important therapeutic implications. Here, we review recent experimental and computational techniques to reverse engineer cancer gene networks using transcriptome and epigenome data. New algorithms, data integration strategies, and increasing amounts of single cell genomics data provide exciting opportunities to model dynamic regulatory states at unprecedented resolution.

Bdf1 Bromodomains Are Essential for Meiosis and the Expression of Meiotic-Specific Genes

Bromodomain and Extra-terminal motif (BET) proteins play a central role in transcription regulation and chromatin signalling pathways. They are present in unicellular eukaryotes and in this study, the role of the BET protein Bdf1 has been explored in Saccharomyces cerevisiae. Mutation of Bdf1 bromodomains revealed defects on both the formation of spores and the meiotic progression, blocking cells at the exit from prophase, before the first meiotic division. This phenotype is associated with a massive deregulation of the transcription of meiotic genes and Bdf1 bromodomains are required for appropriate expression of the key meiotic transcription factor NDT80 and almost all the Ndt80-inducible genes, including APC complex components. Bdf1 notably accumulates on the promoter of Ndt80 and its recruitment is dependent on Bdf1 bromodomains. In addition, the ectopic expression of NDT80 during meiosis partially bypasses this dependency. Finally, purification of Bdf1 partners identified two independent complexes with Bdf2 or the SWR complex, neither of which was required to complete sporulation. Taken together, our results unveil a new role for Bdf1 –working independently from its predominant protein partners Bdf2 and the SWR1 complex–as a regulator of meiosis-specific genes.

Chromatin modifying proteins play a central role in transcription regulation and chromatin signalling. In this study we investigated the functional role of the bromodomains of the chromatin protein Bdf1 during yeast gametogenesis. Our results show that the bromodomains of Bdf1 are essential for meiotic progression and the formation of mature spores. Bdf1 bromodomains are required for the expression of key meiotic genes and the master regulator NDT80. Forced expression of NDT80 can partially rescue the formation of spores when Bdf1 bromodomains are mutated. The results presented here indicate that Bdf1 forms two exclusive complexes, with Bdf2 or with the SWR complex. However, none of these complexes are required for sporulation progression. To conclude, our findings suggest that Bdf1 is a new regulator of the meiotic transcription program and of the expression of the master regulator NDT80.

DNA methylation patterns associated with oxidative stress in an ageing population

Background

Oxidative stress has been related to type 2 diabetes (T2D) and cardiovascular disease (CVD), the leading global cause of death. Contributions of environmental factors such as oxidative stress on complex traits and disease may be partly mediated through changes in epigenetic marks (e.g. DNA methylation). Studies relating differential methylation with intermediate phenotypes and disease endpoints may be useful in identifying additional candidate genes and mechanisms involved in disease.

Methods

To investigate the role of epigenetic variation in oxidative stress marker levels and subsequent development of CVD and T2D, we performed analyses of genome-wide DNA methylation in blood, ten markers of oxidative stress (total glutathione [TGSH], reduced glutathione [GSH], oxidised glutathione [GSSG], GSSG to GSH ratio, homocysteine [HCY], oxidised low-density lipoprotein (oxLDL), antibodies against oxLDL [OLAB], conjugated dienes [CD], baseline conjugated dienes [BCD]-LDL and total antioxidant capacity [TAOC]) and incident disease in up to 966 age-matched individuals.

Results

In total, we found 66 cytosine-guanine (CpG) sites associated with one or more oxidative stress markers (false discovery rate [FDR] <0.05). These sites were enriched in regulatory regions of the genome. Genes annotated to CpG sites showed enrichment in annotation clusters relating to phospho-metabolism and proteins with pleckstrin domains. We investigated the contribution of oxidative stress-associated CpGs to development of cardiometabolic disease. Methylation variation at CpGs in the 3'-UTR of HIST1H4D (cg08170869; histone cluster 1, H4d) and in the body of DVL1 (cg03465880; dishevelled-1) were associated with incident T2D events during 10 years of follow-up (all permutation p-values <0.01), indicating a role of epigenetic regulation in oxidative stress processes leading to development or progression of diabetes. Methylation QTL (meQTL) analysis showed significant associations with genetic sequence variants in cis at 28 (42%) of oxidative stress phenotype-associated sites (FDR < 0.05). Integrating cis-meQTLs with genotype-phenotype associations indicated that genetic effects on oxidative stress phenotype at one locus (cg07547695; BCL2L11) may be mediated through DNA methylation.

Conclusions

In conclusion, we report novel associations of DNA methylation with oxidative stress, some of which also show evidence of a relation with T2D incidence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-016-0235-0) contains supplementary material, which is available to authorized users.

Actin is an evolutionarily-conserved damage-associated molecular pattern that signals tissue injury in Drosophila melanogaster

Damage-associated molecular patterns (DAMPs) are molecules released by dead cells that trigger sterile inflammation and, in vertebrates, adaptive immunity. Actin is a DAMP detected in mammals by the receptor, DNGR-1, expressed by dendritic cells (DCs). DNGR-1 is phosphorylated by Src-family kinases and recruits the tyrosine kinase Syk to promote DC cross-presentation of dead cell-associated antigens. Here we report that actin is also a DAMP in invertebrates that lack DCs and adaptive immunity. Administration of actin to Drosophila melanogaster triggers a response characterised by selective induction of STAT target genes in the fat body through the cytokine Upd3 and its JAK/STAT-coupled receptor, Domeless. Notably, this response requires signalling via Shark, the Drosophila orthologue of Syk, and Src42A, a Drosophila Src-family kinase, and is dependent on Nox activity. Thus, extracellular actin detection via a Src-family kinase-dependent cascade is an ancient means of detecting cell injury that precedes the evolution of adaptive immunity.

Distinct Cell Stress Responses Induced by ATP Restriction in Quiescent Human Fibroblasts

Quiescence is the prevailing state of many cell types under homeostatic conditions. Yet, surprisingly little is known about how quiescent cells respond to energetic and metabolic challenges. To better understand compensatory responses of quiescent cells to metabolic stress, we established, in human primary dermal fibroblasts, an experimental ‘energy restriction’ model. Quiescence was achieved by short-term culture in serum-deprived media and ATP supply restricted using a combination of glucose transport inhibitors and mitochondrial uncouplers. In aggregate, these measures led to markedly reduced intracellular ATP levels while not compromising cell viability over the observation period of 48 h. Analysis of the transcription factor (TF) landscape induced by this treatment revealed alterations in several signal transduction nodes beyond the expected biosynthetic adaptations. These included increased abundance of NF-κB regulated TFs and altered TF subsets regulated by Akt and p53. The observed changes in gene regulation and corresponding alterations in key signaling nodes are likely to contribute to cell survival at intracellular ATP concentrations substantially below those achieved by growth factor deprivation alone. This experimental model provides a benchmark for the investigation of cell survival pathways and related molecular targets that are associated with restricted energy supply associated with biological aging and metabolic diseases.

The Arabidopsis transcription factor ABIG1 relays ABA signaled growth inhibition and drought induced senescence

Drought inhibits plant growth and can also induce premature senescence. Here we identify a transcription factor, ABA INSENSITIVE GROWTH 1 (ABIG1) required for abscisic acid (ABA) mediated growth inhibition, but not for stomatal closure. ABIG1 mRNA levels are increased both in response to drought and in response to ABA treatment. When treated with ABA, abig1 mutants remain greener and produce more leaves than comparable wild-type plants. When challenged with drought, abig1 mutants have fewer yellow, senesced leaves than wild-type. Induction of ABIG1 transcription mimics ABA treatment and regulates a set of genes implicated in stress responses. We propose a model in which drought acts through ABA to increase ABIG1 transcription which in turn restricts new shoot growth and promotes leaf senescence. The results have implications for plant breeding: the existence of a mutant that is both ABA resistant and drought resistant points to new strategies for isolating drought resistant genetic varieties.

EpiMINE, a computational program for mining epigenomic data

Background

In epigenetic research, both the increasing ease of high-throughput sequencing and a greater interest in genome-wide studies have resulted in an exponential flooding of epigenetic-related data in public domain. This creates an opportunity for exploring data outside the limits of any specific query-centred study. Such data have to undergo standard primary analyses that are accessible with multiple well-stabilized programs. Further downstream analyses, such as genome-wide comparative, correlative and quantitative analyses, are critical in deciphering key biological features. However, these analyses are only accessible for computational researchers and completely lack platforms capable of handling, analysing and linking multiple interdisciplinary datasets with efficient analytical methods.

Results

Here, we present EpiMINE, a program for mining epigenomic data. It is a user-friendly, stand-alone computational program designed to support multiple datasets, for performing genome-wide correlative and quantitative analysis of ChIP-seq and RNA-seq data. Using data available from the ENCODE project, we illustrated several features of EpiMINE through different biological scenarios to show how easy some known observations can be verified. These results highlight how these approaches can be helpful in identifying novel biological features.

Conclusions

EpiMINE performs different kinds of genome-wide quantitative and correlative analyses, using ChIP-seq- and RNA-seq-related datasets. Its framework enables it to be used by both experimental and computational researchers. EpiMINE can be downloaded from https://sourceforge.net/projects/epimine/.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0095-z) contains supplementary material, which is available to authorized users.

Sense-antisense gene-pairs in breast cancer and associated pathological pathways

More than 30% of human protein-coding genes form hereditary complex genome architectures composed of sense-antisense (SA) gene pairs (SAGPs) transcribing their RNAs from both strands of a given locus. Such architectures represent important novel components of genome complexity contributing to gene expression deregulation in cancer cells. Therefore, the architectures might be involved in cancer pathways and, in turn, be used for novel drug targets discovery. However, the global roles of SAGPs in cancer pathways has not been studied. Here we investigated SAGPs associated with breast cancer (BC)-related pathways using systems biology, prognostic survival and experimental methods. Gene expression analysis identified 73 BC-relevant SAGPs that are highly correlated in BC. Survival modelling and metadata analysis of the 1161 BC patients allowed us to develop a novel patient prognostic grouping method selecting the 12 survival-significant SAGPs. The qRT-PCR-validated 12-SAGP prognostic signature reproducibly stratified BC patients into low- and high-risk prognostic subgroups. The 1381 SAGP-defined differentially expressed genes common across three studied cohorts were identified. The functional enrichment analysis of these genes revealed the GABPA gene network, including BC-relevant SAGPs, specific gene sets involved in cell cycle, spliceosomal and proteasomal pathways. The co-regulatory function of GABPA in BC cells was supported using siRNA knockdown studies. Thus, we demonstrated SAGPs as the synergistically functional genome architectures interconnected with cancer-related pathways and associated with BC patient clinical outcomes. Taken together, SAGPs represent an important component of genome complexity which can be used to identify novel aspects of coordinated pathological gene networks in cancers.

Transcriptomic analysis of mouse liver reveals a potential hepato-enteric pathogenic mechanism in acute Toxoplasma gondii infection

Background

Toxoplasma gondii is a worldwide spread pathogen which can infect all tissues of its host. The transcriptomic responses of infected brain and spleen have been reported. However, our knowledge of the global transcriptomic change in infected liver is limited. Additionally, T. gondii infection represents a highly dynamic process involving complex biological responses of the host at many levels. Herein, we describe such processes at a global level by discovering gene expression changes in mouse livers after acute infection with T. gondii ToxoDB#9 strain.

Results

Global transcriptomic analysis identified 2,758 differentially expressed transcripts in infected liver, of which 1,356 were significantly downregulated and 1,402 upregulated. GO and KEGG database analyses showed that host immune responses were upregulated, while the metabolic-related processes/pathways were downregulated, especially xenobiotic metabolism, fatty acid metabolism, energy metabolism, and bile biosynthesis and secretion. The metabolism of more than 800 chemical compounds including anti-Toxoplasma prescribed medicines were predicted to be modulated during acute T. gondii infection due to the downregulation of enzymes involved in xenobiotic metabolism.

Conclusions

To the best of our knowledge, this is the first global transcriptomic analysis of mouse liver infected by T. gondii. The present data indicate that during the early stage of liver infection, T. gondii can induce changes in liver xenobiotic metabolism, upregulating inflammatory response and downregulating hepatocellular PPAR signaling pathway, altering host bile biosynthesis and secretion pathway; these changes could enhance host intestinal dysbacteriosis and thus contribute to the pathological changes of both liver and intestine of infected mice. These findings describe the biological changes in infected liver, providing a potential mechanistic pathway that links hepatic and intestinal pathologies to T. gondii infection.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1716-x) contains supplementary material, which is available to authorized users.

Time-course gene profiling and networks in demethylated retinoblastoma cell line

Retinoblastoma, a very aggressive cancer of the developing retina, initiatiates by the biallelic loss of RB1 gene, and progresses very quickly following RB1 inactivation. While its genome is stable, multiple pathways are deregulated, also epigenetically. After reviewing the main findings in relation with recently validated markers, we propose an integrative bioinformatics approach to include in the previous group new markers obtained from the analysis of a single cell line subject to epigenetic treatment. In particular, differentially expressed genes are identified from time course microarray experiments on the WERI-RB1 cell line treated with 5-Aza-2′-deoxycytidine (decitabine; DAC). By inducing demethylation of CpG island in promoter genes that are involved in biological processes, for instance apoptosis, we performed the following main integrative analysis steps: i) Gene expression profiling at 48h, 72h and 96h after DAC treatment; ii) Time differential gene co-expression networks and iii) Context-driven marker association (transcriptional factor regulated protein networks, master regulatory paths). The observed DAC-driven temporal profiles and regulatory connectivity patterns are obtained by the application of computational tools, with support from curated literature. It is worth emphasizing the capacity of networks to reconcile multi-type evidences, thus generating testable hypotheses made available by systems scale predictive inference power. Despite our small experimental setting, we propose through such integrations valuable impacts of epigenetic treatment in terms of gene expression measurements, and then validate evidenced apoptotic effects.

Lysine-specific histone demethylase 1 inhibition promotes reprogramming by facilitating the expression of exogenous transcriptional factors and metabolic switch

Lysine-specific histone demethylase 1 (LSD1) regulates histone methylation and influences the epigenetic state of cells during the generation of induced pluripotent stem cells (iPSCs). Here we reported that LSD1 inhibition via shRNA or specific inhibitor, tranylcypromine, promoted reprogramming at early stage via two mechanisms. At early stage of reprogramming, LSD1 inhibition increased the retrovirus-mediated exogenous expression of Oct4, Klf4, and Sox2 by blocking related H3K4 demethylation. Since LSD1 inhibition still promoted reprogramming even when iPSCs were induced with small-molecule compounds in a virus-free system, additional mechanisms should be involved. When RNA-seq was used for analysis, it was found that LSD1 inhibition reversed some gene expression changes induced by OKS, which subsequently promoted reprogramming. For example, by partially rescuing the decreased expression of Hif1α, LSD1 inhibition reversed the up-regulation of genes in oxidative phosphorylation pathway and the down-regulation of genes in glycolysis pathway. Such effects facilitated the metabolic switch from oxidative phosphorylation to glycolysis and subsequently promoted iPSCs induction. In addition, LSD1 inhibition also promoted the conversion from pre-iPSCs to iPSCs by facilitating the similar metabolic switch. Therefore, LSD1 inhibition promotes reprogramming by facilitating the expression of exogenous transcriptional factors and metabolic switch.

Coordinate Regulation of Metabolite Glycosylation and Stress Hormone Biosynthesis by TT8 in Arabidopsis

Secondary metabolites play a key role in coordinating ecology and defense strategies of plants. Diversity of these metabolites arises by conjugation of core structures with diverse chemical moieties, such as sugars in glycosylation. Active pools of phytohormones, including those involved in plant stress response, are also regulated by glycosylation. While much is known about the enzymes involved in glycosylation, we know little about their regulation or coordination with other processes. We characterized the flavonoid pathway transcription factor TRANSPARENT TESTA8 (TT8) in Arabidopsis (Arabidopsis thaliana) using an integrative omics strategy. This approach provides a systems-level understanding of the cellular machinery that is used to generate metabolite diversity by glycosylation. Metabolomics analysis of TT8 loss-of-function and inducible overexpression lines showed that TT8 coordinates glycosylation of not only flavonoids, but also nucleotides, thus implicating TT8 in regulating pools of activated nucleotide sugars. Transcriptome and promoter network analyses revealed that the TT8 regulome included sugar transporters, proteins involved in sugar binding and sequestration, and a number of carbohydrate-active enzymes. Importantly, TT8 affects stress response, along with brassinosteroid and jasmonic acid biosynthesis, by directly binding to the promoters of key genes of these processes. This combined effect on metabolite glycosylation and stress hormones by TT8 inducible overexpression led to significant increase in tolerance toward multiple abiotic and biotic stresses. Conversely, loss of TT8 leads to increased sensitivity to these stresses. Thus, the transcription factor TT8 is an integrator of secondary metabolism and stress response. These findings provide novel approaches to improve broad-spectrum stress tolerance.

Transcriptome sequencing wide functional analysis of human mesenchymal stem cells in response to TLR4 ligand

Due to their multipotentiality and immunomodulation, human mesenchymal stem cells (hMSCs) are widely studied for the treatment of degenerative and inflammatory diseases. Transplantation of hMSCs to damaged tissue is a promising approach for tissue regeneration. However, the physiological mechanisms and regulatory processes of MSC trafficking to injured tissue are largely unexplored. Here, we evaluated the gene expression profile and migratory potential of hMSCs upon stimulation with the TLR4 ligand lipopolysaccharide (LPS). Using RNA sequencing, we identified unique induction patterns of interferon stimulated genes, cytokines and chemokines involved in chemotaxis and homing. The −950 to +50 bp regions of many of these LPS-responsive genes were enriched with putative binding motifs for the transcription factors (TFs) interferon regulatory factor (IRF1) and nuclear factor kappa B (NF-κB1, REL), which were also induced by LPS along with other TFs. Chromatin immunoprecipitation assays showed that IRF1 bound within their target genes promoter region. In addition, IRF1 attenuation significantly down-regulated interferon stimulated genes as well as key cytokines. Furthermore, using pharmacological inhibitors, we showed that the NF-κB and phosphatidylinositol 3-kinase (PI3K) pathways regulate the migratory and cytokines/chemokines response to LPS. These unprecedented data suggest that IRF1 and NF-κB orchestrate the TLR4-primed immunomodulatory response of hMSCs and that this response also involves the PI3K pathway.

Transcriptome sequencing reveals that LPS-triggered transcriptional responses in established microglia BV2 cell lines are poorly representative of primary microglia

Background

Microglia are resident myeloid cells in the CNS that are activated by infection, neuronal injury, and inflammation. Established BV2 microglial cell lines have been the primary in vitro models used to study neuroinflammation for more than a decade because they reduce the requirement of continuously maintaining cell preparations and animal experimentation models. However, doubt has recently been raised regarding the value of BV2 cell lines as a model system.

Methods

We used triplicate RNA sequencing (RNA-seq) to investigate the molecular signature of primary and BV2 microglial cell lines using two transcriptomic techniques: global transcriptomic biological triplicate RNA-seq and quantitative real-time PCR. We analyzed differentially expressed genes (DEGs) to identify transcription factor (TF) motifs (−950 to +50 bp of the 5′ upstream promoters) and epigenetic mechanisms.

Results

Sequencing assessment and quality evaluation revealed that primary microglia have a distinct transcriptomic signature and express a unique cluster of transcripts in response to lipopolysaccharide. This microglial signature was not observed in BV2 microglial cell lines. Importantly, we observed that previously unidentified TFs (i.e., IRF2, IRF5, IRF8, STAT1, STAT2, and STAT5A) and the epigenetic regulators KDM1A, NSD3, and SETDB2 were significantly and selectively expressed in primary microglia (PM). Although transcriptomic alterations known to occur in BV2 microglial cell lines were identified in PM, we also observed several novel transcriptomic alterations in PM that are not frequently observed in BV2 microglial cell lines.

Conclusions

Collectively, these unprecedented findings demonstrate that established BV2 microglial cell lines are probably a poor representation of PM, and we establish a resource for future studies of neuroinflammation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12974-016-0644-1) contains supplementary material, which is available to authorized users.

Primary cancer cell culture: mammary-optimized vs conditional reprogramming

The impact of different culture conditions on biology of primary cancer cells is not always addressed. Here, conditional reprogramming (CRC) was compared with mammary-optimized EpiCult-B (EpiC) for primary mammary epithelial cell isolation and propagation, allograft generation, and genome-wide transcriptional consequences using cancer and non-cancer mammary tissue from mice with different dosages of Brca1 and p53 Selective comparison to DMEM was included. Primary cultures were established with all three media, but CRC was most efficient for initial isolation (P<0.05). Allograft development was faster using cells grown in EpiC compared with CRC (P<0.05). Transcriptome comparison of paired CRC and EpiC cultures revealed 1700 differentially expressed genes by passage 20. CRC promoted Trp53 gene family upregulation and increased expression of epithelial differentiation genes, whereas EpiC elevated expression of epithelial-mesenchymal transition genes. Differences did not persist in allografts where both methods yielded allografts with relatively similar transcriptomes. Restricting passage (<7) reduced numbers of differentially expressed genes below 50. In conclusion, CRC was most efficient for initial cell isolation but EpiC was quicker for allograft generation. The extensive culture-specific gene expression patterns that emerged with longer passage could be limited by reducing passage number when both culture transcriptomes were equally similar to that of the primary tissue. Defining impact of culture condition and passage on the transcriptome of primary cells could assist experimental design and interpretation. For example, differences that appear with passage and culture condition are potentially exploitable for comparative studies targeting specific biological networks in different transcriptional environments.

ADAR1 deletion induces NFκB and interferon signaling dependent liver inflammation and fibrosis

Adenosine deaminase acting on RNA (ADAR) 1 binds and edits double-stranded (ds) RNA secondary structures found mainly within untranslated regions of many transcripts. In the current research, our aim was to study the role of ADAR1 in liver homeostasis. As previous studies show a conserved immunoregulatory function for ADAR1 in mammalians, we focused on its role in preventing chronic hepatic inflammation and the associated activation of hepatic stellate cells to produce extracellular matrix and promote fibrosis. We show that hepatocytes specific ADAR1 knock out (KO) mice display massive liver damage with multifocal inflammation and fibrogenesis. The bioinformatics analysis of the microarray gene-expression datasets of ADAR1 KO livers reveled a type-I interferons signature and an enrichment for immune response genes compared to control littermate livers. Furthermore, we found that in vitro silencing of ADAR1 expression in HepG2 cells leads to enhanced transcription of NFκB target genes, foremost of the pro-inflammatory cytokines IL6 and IL8. We also discovered immune cell-independent paracrine signaling among ADAR1-depleted HepG2 cells and hepatic stellate cells, leading to the activation of the latter cell type to adopt a profibrogenic phenotype. This paracrine communication dependent mainly on the production and secretion of the cytokine IL6 induced by ADAR1 silencing in hepatocytes. Thus, our findings shed a new light on the vital regulatory role of ADAR1 in hepatic immune homeostasis, chiefly its inhibitory function on the crosstalk between the NFκB and type-I interferons signaling cascades, restraining the development of liver inflammation and fibrosis.

Inhibition of the hexosamine biosynthetic pathway promotes castration-resistant prostate cancer

The precise molecular alterations driving castration-resistant prostate cancer (CRPC) are not clearly understood. Using a novel network-based integrative approach, here, we show distinct alterations in the hexosamine biosynthetic pathway (HBP) to be critical for CRPC. Expression of HBP enzyme glucosamine-phosphate N-acetyltransferase 1 (GNPNAT1) is found to be significantly decreased in CRPC compared with localized prostate cancer (PCa). Genetic loss-of-function of GNPNAT1 in CRPC-like cells increases proliferation and aggressiveness, in vitro and in vivo. This is mediated by either activation of the PI3K-AKT pathway in cells expressing full-length androgen receptor (AR) or by specific protein 1 (SP1)-regulated expression of carbohydrate response element-binding protein (ChREBP) in cells containing AR-V7 variant. Strikingly, addition of the HBP metabolite UDP-N-acetylglucosamine (UDP-GlcNAc) to CRPC-like cells significantly decreases cell proliferation, both in-vitro and in animal studies, while also demonstrates additive efficacy when combined with enzalutamide in-vitro. These observations demonstrate the therapeutic value of targeting HBP in CRPC.

The molecular alterations driving anti-androgen resistance in prostate cancer are unclear. Here, the authors show, using a network-based approach, that inhibition of the hexosamine biosynthetic pathway is necessary to develop resistance and that increasing the activity of the pathway enhances the anti-androgen response.

TET2 Regulates Mast Cell Differentiation and Proliferation through Catalytic and Non-catalytic Activities

Dioxygenases of the TET family impact genome functions by converting 5-methylcytosine (5mC) in DNA to 5-hydroxymethylcytosine (5hmC). Here, we identified TET2 as a crucial regulator of mast cell differentiation and proliferation. In the absence of TET2, mast cells showed disrupted gene expression and altered genome-wide 5hmC deposition, especially at enhancers and in the proximity of downregulated genes. Impaired differentiation of Tet2-ablated cells could be relieved or further exacerbated by modulating the activity of other TET family members, and mechanistically it could be linked to the dysregulated expression of C/EBP family transcription factors. Conversely, the marked increase in proliferation induced by the loss of TET2 could be rescued exclusively by re-expression of wild-type or catalytically inactive TET2. Our data indicate that, in the absence of TET2, mast cell differentiation is under the control of compensatory mechanisms mediated by other TET family members, while proliferation is strictly dependent on TET2 expression.

Repeated FcεRI triggering reveals modified mast cell function related to chronic allergic responses in tissue

BACKGROUND

Activation of mast cells through FcεRI plays an important role in acute allergic reactions. However, little is known about the function of mast cells in patients with chronic allergic inflammation or the effect of repeated FcεRI triggering occurring in such responses.

OBJECTIVE

We aimed to identify changes in mast cell function after repeated FcεRI triggering and to correlate these changes to chronic allergic responses in tissue.

METHODS

Human cord blood-derived mast cells were treated for 2 weeks with anti-IgE. The function of naive or treated mast cells was analyzed by means of RNA sequencing, quantitative RT-PCR, flow cytometry, and functional assays. Protein secretion was measured with ELISAs and multiplex assays.

RESULTS

We observed several changes in mast cell function after repeated anti-IgE triggering. Although the acute response was dampened, we identified 289 genes significantly upregulated after repeated anti-IgE. Most of these genes (84%) were not upregulated after a single anti-IgE stimulus, indicating a significantly different response mode characterized by increased antigen presentation, response to bacteria, and chemotaxis. Changes in mast cell function were related to changes in expression of the transcription factors RXRA and BATF and others. Importantly, we found a substantial overlap between genes upregulated after repeated anti-IgE triggering and genes upregulated in tissue from patients with chronic allergy, in particular those of patients with chronic rhinosinusitis.

CONCLUSION

Our study provides evidence for intrinsic modulation of mast cell function on repeated FcεRI-mediated activation. The overlap with gene expression in tissues is suggestive of a direct link between repeated IgE-mediated activation of mast cells and chronic allergy.

Proteomic Profiling of Mouse Liver following Acute Toxoplasma gondii Infection

Toxoplasma gondii remains a global public health problem. However, its pathophysiology is still not-completely understood particularly the impact of infection on host liver metabolism. We performed iTRAQ-based proteomic analysis to evaluate early liver protein responses in BALB/c mice following infection with T. gondii PYS strain (genotype ToxoDB#9) infection. Our data revealed modification of protein expression in key metabolic pathways, as indicated by the upregulation of immune response and downregulation of mitochondrial respiratory chain, and the metabolism of fatty acids, lipids and xenobiotics. T. gondii seems to hijack host PPAR signaling pathway to downregulate the metabolism of fatty acids, lipids and energy in the liver. The metabolism of over 400 substances was affected by the downregulation of genes involved in xenobiotic metabolism. The top 10 transcription factors used by upregulated genes were Stat2, Stat1, Irf2, Irf1, Sp2, Egr1, Stat3, Klf4, Elf1 and Gabpa, while the top 10 transcription factors of downregulated genes were Hnf4A, Ewsr1, Fli1, Hnf4g, Nr2f1, Pparg, Rxra, Hnf1A, Foxa1 and Foxo1. These findings indicate global reprogramming of the metabolism of the mouse liver after acute T. gondii infection. Functional characterization of the altered proteins may enhance understanding of the host responses to T. gondii infection and lead to the identification of new therapeutic targets.

RNA Sequencing Reveals the Alteration of the Expression of Novel Genes in Ethanol-Treated Embryoid Bodies

Fetal alcohol spectrum disorder is a collective term representing fetal abnormalities associated with maternal alcohol consumption. Prenatal alcohol exposure and related anomalies are well characterized, but the molecular mechanism behind this phenomenon is not well characterized. In this present study, our aim is to profile important genes that regulate cellular development during fetal development. Human embryonic carcinoma cells (NCCIT) are cultured to form embryoid bodies and then treated in the presence and absence of ethanol (50 mM). We employed RNA sequencing to profile differentially expressed genes in the ethanol-treated embryoid bodies from NCCIT vs. EB, NCCIT vs. EB+EtOH and EB vs. EB+EtOH data sets. A total of 632, 205 and 517 differentially expressed genes were identified from NCCIT vs. EB, NCCIT vs. EB+EtOH and EB vs. EB+EtOH, respectively. Functional annotation using bioinformatics tools reveal significant enrichment of differential cellular development and developmental disorders. Furthermore, a group of 42, 15 and 35 transcription factor-encoding genes are screened from all of the differentially expressed genes obtained from NCCIT vs. EB, NCCIT vs. EB+EtOH and EB vs. EB+EtOH, respectively. We validated relative gene expression levels of several transcription factors from these lists by quantitative real-time PCR. We hope that our study substantially contributes to the understanding of the molecular mechanism underlying the pathology of alcohol-mediated anomalies and ease further research.

A Stringent Systems Approach Uncovers Gene-Specific Mechanisms Regulating Inflammation

Much has been learned about transcriptional cascades and networks from large-scale systems analyses of high-throughput datasets. However, analysis methods that optimize statistical power through simultaneous evaluation of thousands of ChIP-seq peaks or differentially expressed genes possess substantial limitations in their ability to uncover mechanistic principles of transcriptional control. By examining nascent transcript RNA-seq, ChIP-seq, and binding motif datasets from lipid A-stimulated macrophages with increased attention to the quantitative distribution of signals, we identified unexpected relationships between the in vivo binding properties of inducible transcription factors, motif strength, and transcription. Furthermore, rather than emphasizing common features of large clusters of co-regulated genes, our results highlight the extent to which unique mechanisms regulate individual genes with key biological functions. Our findings demonstrate the mechanistic value of stringent interrogation of well-defined sets of genes as a complement to broader systems analyses of transcriptional cascades and networks.

Analysis of Genomic Sequence Motifs for Deciphering Transcription Factor Binding and Transcriptional Regulation in Eukaryotic Cells

Eukaryotic genomes contain a variety of structured patterns: repetitive elements, binding sites of DNA and RNA associated proteins, splice sites, and so on. Often, these structured patterns can be formalized as motifs and described using a proper mathematical model such as position weight matrix and IUPAC consensus. Two key tasks are typically carried out for motifs in the context of the analysis of genomic sequences. These are: identification in a set of DNA regions of over-represented motifs from a particular motif database, and de novo discovery of over-represented motifs. Here we describe existing methodology to perform these two tasks for motifs characterizing transcription factor binding. When applied to the output of ChIP-seq and ChIP-exo experiments, or to promoter regions of co-modulated genes, motif analysis techniques allow for the prediction of transcription factor binding events and enable identification of transcriptional regulators and co-regulators. The usefulness of motif analysis is further exemplified in this review by how motif discovery improves peak calling in ChIP-seq and ChIP-exo experiments and, when coupled with information on gene expression, allows insights into physical mechanisms of transcriptional modulation.

Targeted activation of CREB in reactive astrocytes is neuroprotective in focal acute cortical injury

The clinical challenge in acute injury as in traumatic brain injury (TBI) is to halt the delayed neuronal loss that occurs hours and days after the insult. Here we report that the activation of CREB-dependent transcription in reactive astrocytes prevents secondary injury in cerebral cortex after experimental TBI. The study was performed in a novel bitransgenic mouse in which a constitutively active CREB, VP16-CREB, was targeted to astrocytes with the Tet-Off system. Using histochemistry, qPCR, and gene profiling we found less neuronal death and damage, reduced macrophage infiltration, preserved mitochondria, and rescued expression of genes related to mitochondrial metabolism in bitransgenic mice as compared to wild type littermates. Finally, with meta-analyses using publicly available databases we identified a core set of VP16-CREB candidate target genes that may account for the neuroprotective effect. Enhancing CREB activity in astrocytes thus emerges as a novel avenue in acute brain post-injury therapeutics.

Aging-associated DNA methylation changes in middle-aged individuals: the Young Finns study

Background

Chronological aging-associated changes in the human DNA methylome have been studied by multiple epigenome-wide association studies (EWASs). Certain CpG sites have been identified as aging-associated in multiple studies, and the majority of the sites identified in various studies show common features regarding location and direction of the methylation change. However, as a whole, the sets of aging-associated CpGs identified in different studies, even with similar tissues and age ranges, show only limited overlap. In this study, we further explore and characterize CpG sites that show close relationship between their DNA methylation level and chronological age during adulthood and which bear the relationship regardless of blood cell type heterogeneity.

Results

In this study, with a multivariable regression model adjusted for cell type heterogeneity, we identified 1202 aging-associated CpG sites (a-CpGs, FDR < 5 %), in whole blood in a population with an especially narrow age range (40 - 49 years). Repeatedly reported a-CpGs located in genes ELOVL2, FHL2, PENK and KLF14 were also identified. Regions with aging-associated hypermethylation were enriched regarding several gene ontology (GO) terms (especially in the cluster of developmental processes), whereas hypomethylated sites showed no enrichment. The genes with higher numbers of a-CpG hits were more often hypermethylated with advancing age. The comparison analysis revealed that of the 1202 a-CpGs identified in the present study, 987 were identified as differentially methylated also between nonagenarians and young adults in a previous study (The Vitality 90+ study), and importantly, the directions of changes were identical in the previous and in the present study.

Conclusions

Here we report that aging-associated DNA methylation features can be identified in a middle-aged population with an age range of only 9 years. A great majority of these sites have been previously reported as aging-associated in a population aged 19 to 90 years. Aging is associated with different types of changes in DNA methylation, clock-like as well as random. We speculate that the a-CpGs identified here in a population with a narrow age-range represent clock-like changes, as they showed concordant methylation behavior in population spanning whole adulthood as well.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2421-z) contains supplementary material, which is available to authorized users.

Dissection of transcriptional and cis-regulatory control of differentiation in human pancreatic cancer

The histological grade of carcinomas describes the ability of tumor cells to organize in differentiated epithelial structures and has prognostic and therapeutic impact. Here, we show that differential usage of the genomic repertoire of transcriptional enhancers leads to grade-specific gene expression programs in human pancreatic ductal adenocarcinoma (PDAC). By integrating gene expression profiling, epigenomic footprinting, and loss-of-function experiments in PDAC cell lines of different grade, we identified the repertoires of enhancers specific to high- and low-grade PDACs and the cognate set of transcription factors acting to maintain their activity. Among the candidate regulators of PDAC differentiation, KLF5 was selectively expressed in pre-neoplastic lesions and low-grade primary PDACs and cell lines, where it maintained the acetylation of grade-specific enhancers, the expression of epithelial genes such as keratins and mucins, and the ability to organize glandular epithelia in xenografts. The identification of the transcription factors controlling differentiation in PDACs will help clarify the molecular bases of its heterogeneity and progression.

Identification of key regulators and their controlling mechanism in a combinatorial apoptosis network: a systems biology approach

NFKB1, SP1 and hsa-let-7a, were identified as key regulators of apoptosis, by network theory through probability of signal propagation, hub-removal and motif analysis.

iRegulon and i-cisTarget: Reconstructing Regulatory Networks Using Motif and Track Enrichment

Gene expression profiling is often used to identify genes that are co-expressed in a biological process or disease. Downstream analyses of co-expressed gene sets using bioinformatics methods can reveal candidate transcription factors (TF) that co-regulate these genes, based on the presence of shared TF binding sites. Drawing gene regulatory networks that connect TFs to their predicted target genes can uncover gene modules that implement a particular function. Here, we describe several protocols to analyze any set of co-expressed genes using iRegulon and i-cisTarget. These tools perform regulatory sequence analysis (motif discovery) and integrate and mine large collections of existing regulatory data, such as ChIP-Seq, DHS-seq, and FAIRE-seq (track discovery). While iRegulon focuses on sets of co-expressed genes, i-cisTarget also analyses genomic regions as input. The following protocols describe how to install and use these tools, how to interpret the obtained results, and will thus help to create meaningful regulatory networks.

Nuclear Factor κB1/RelA Mediates Inflammation in Human Lung Epithelial Cells at Atmospheric Oxygen Levels

Oxygen levels range from 2% to 9% in vivo. Atmospheric O2 levels (21%) are known to induce cell proliferation defects and cellular senescence in primary cell cultures. However, the mechanistic basis of the deleterious effects of higher O2 levels is not fully understood. On the other hand, immortalized cells including cancer cell lines, which evade cellular senescence are normally cultured at 21% O2 and the effects of higher O2 on these cells are understudied. Here, we addressed this problem by culturing immortalized human bronchial epithelial (BEAS-2B) cells at ambient atmospheric, 21% O2 and lower, 10% O2. Our results show increased inflammatory response at 21% O2 but not at 10% O2. We found higher RelA binding at the NF-κB1/RelA target gene promoters as well as upregulation of several pro-inflammatory cytokines in cells cultured at 21% O2. RelA knockdown prevented the upregulation of the pro-inflammatory cytokines at 21% O2, suggesting NF-κB1/RelA as a major mediator of inflammatory response in cells cultured at 21% O2. Interestingly, unlike the 21% O2 cultured cells, exposure of 10% O2 cultured cells to H2O2 did not elicit inflammatory response, suggesting increased ability to tolerate oxidative stress in cells cultured at lower O2 levels.

Fine-tuned SRF activity controls asymmetrical neuronal outgrowth: implications for cortical migration, neural tissue lamination and circuit assembly

The stimulus-regulated transcription factor Serum Response Factor (SRF) plays an important role in diverse neurodevelopmental processes related to structural plasticity and motile functions, although its precise mechanism of action has not yet been established. To further define the role of SRF in neural development and distinguish between cell-autonomous and non cell-autonomous effects, we bidirectionally manipulated SRF activity through gene transduction assays that allow the visualization of individual neurons and their comparison with neighboring control cells. In vitro assays showed that SRF promotes survival and filopodia formation and is required for normal asymmetric neurite outgrowth, indicating that its activation favors dendrite enlargement versus branching. In turn, in vivo experiments demonstrated that SRF-dependent regulation of neuronal morphology has important consequences in the developing cortex and retina, affecting neuronal migration, dendritic and axonal arborization and cell positioning in these laminated tissues. Overall, our results show that the controlled and timely activation of SRF is essential for the coordinated growth of neuronal processes, suggesting that this event regulates the switch between neuronal growth and branching during developmental processes.

Loss of function of myosin chaperones triggers Hsf1-mediated transcriptional response in skeletal muscle cells

Background

Mutations in myosin chaperones Unc45b and Hsp90aa1.1 as well as in the Unc45b-binding protein Smyd1b impair formation of myofibrils in skeletal muscle and lead to the accumulation of misfolded myosin. The concomitant transcriptional response involves up-regulation of the three genes encoding these proteins, as well as genes involved in muscle development. The transcriptional up-regulation of unc45b, hsp90aa1.1 and smyd1b is specific to zebrafish mutants with myosin folding defects, and is not triggered in other zebrafish myopathy models.

Results

By dissecting the promoter of unc45b, we identify a Heat shock factor 1 (Hsf1) binding element as a mediator of unc45b up-regulation in myofibers lacking myosin folding proteins. Loss-of-function of Hsf1 abolishes unc45b up-regulation in mutants with defects in myosin folding.

Conclusions

Taken together, our data show that skeletal muscle cells respond to defective myosin chaperones with a complex gene program and suggest that this response is mediated by Hsf1 activation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0825-8) contains supplementary material, which is available to authorized users.

The Regulatory Factor ZFHX3 Modifies Circadian Function in SCN via an AT Motif-Driven Axis

We identified a dominant missense mutation in the SCN transcription factor Zfhx3, termed short circuit (Zfhx3(Sci)), which accelerates circadian locomotor rhythms in mice. ZFHX3 regulates transcription via direct interaction with predicted AT motifs in target genes. The mutant protein has a decreased ability to activate consensus AT motifs in vitro. Using RNA sequencing, we found minimal effects on core clock genes in Zfhx3(Sci/+) SCN, whereas the expression of neuropeptides critical for SCN intercellular signaling was significantly disturbed. Moreover, mutant ZFHX3 had a decreased ability to activate AT motifs in the promoters of these neuropeptide genes. Lentiviral transduction of SCN slices showed that the ZFHX3-mediated activation of AT motifs is circadian, with decreased amplitude and robustness of these oscillations in Zfhx3(Sci/+) SCN slices. In conclusion, by cloning Zfhx3(Sci), we have uncovered a circadian transcriptional axis that determines the period and robustness of behavioral and SCN molecular rhythms.

Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

BACKGROUND

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

METHODS

We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls).

RESULTS

Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network.

CONCLUSION

We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development.

IMPACT

Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization.

Transcriptome sequencing of microglial cells stimulated with TLR3 and TLR4 ligands

Background

Resident macrophages in the CNS microglia become activated and produce proinflammatory molecules upon encountering bacteria or viruses. TLRs are a phylogenetically conserved diverse family of sensors that drive innate immune responses following interactions with PAMPs. TLR3 and TLR4 recognize viral dsRNA Poly (I:C) and bacterial endotoxin LPS, respectively. Importantly, these receptors differ in their downstream adaptor molecules. Thus far, only a few studies have investigated the effects of TLR3 and TLR4 in macrophages. However, a genome-wide search for the effects of these TLRs has not been performed in microglia using RNA-seq. Gene expression patterns were determined for the BV-2 microglial cell line when stimulated with viral dsRNA Poly (I:C) or bacterial endotoxin LPS to identify novel transcribed genes, as well as investigate how differences in downstream signaling could influence gene expression in innate immunity.

Results

Sequencing assessment and quality evaluation revealed that common and unique patterns of proinflammatory genes were significantly up-regulated in response to TLR3 and TLR4 stimulation. However, the IFN/viral response gene showed a stronger response to TLR3 stimulation than to TLR4 stimulation. Unexpectedly, TLR3 and TLR4 stimulation did not activate IFN-ß and IRF3 in BV-2 microglia. Most importantly, we observed that previously unidentified transcription factors (TFs) (i.e., IRF1, IRF7, and IRF9) and the epigenetic regulators KDM4A and DNMT3L were significantly up-regulated in both TLR3- and TLR4-stimulated microglia. We also identified 29 previously unidentified genes that are important in immune regulation. In addition, we confirmed the expressions of key inflammatory genes as well as pro-inflammatory mediators in the supernatants were significantly induced in TLR3-and TLR4-stimulated primary microglial cells. Moreover, transcriptional start sites (TSSs) and isoforms, as well as differential promoter usage, revealed a complex pattern of transcriptional and post-transcriptional gene regulation upon infection with TLR3 and TLR4. Furthermore, TF motif analysis (-950 to +50 bp of the 5′ upstream promoters) revealed that the DNA sequences for NF-κB, IRF1, and STAT1 were significantly enriched in TLR3- and TLR4-stimulated microglia.

Conclusions

These unprecedented findings not only permit a comparison of TLR3-and TLR4-stimulated genes but also identify new genes that have not been previously implicated in innate immunity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1728-5) contains supplementary material, which is available to authorized users.

Age-Related Gene Expression Differences in Monocytes from Human Neonates, Young Adults, and Older Adults

A variety of age-related differences in the innate and adaptive immune systems have been proposed to contribute to the increased susceptibility to infection of human neonates and older adults. The emergence of RNA sequencing (RNA-seq) provides an opportunity to obtain an unbiased, comprehensive, and quantitative view of gene expression differences in defined cell types from different age groups. An examination of ex vivo human monocyte responses to lipopolysaccharide stimulation or Listeria monocytogenes infection by RNA-seq revealed extensive similarities between neonates, young adults, and older adults, with an unexpectedly small number of genes exhibiting statistically significant age-dependent differences. By examining the differentially induced genes in the context of transcription factor binding motifs and RNA-seq data sets from mutant mouse strains, a previously described deficiency in interferon response factor-3 activity could be implicated in most of the differences between newborns and young adults. Contrary to these observations, older adults exhibited elevated expression of inflammatory genes at baseline, yet the responses following stimulation correlated more closely with those observed in younger adults. Notably, major differences in the expression of constitutively expressed genes were not observed, suggesting that the age-related differences are driven by environmental influences rather than cell-autonomous differences in monocyte development.

A novel function of nuclear nonmuscle myosin regulatory light chain in promotion of xanthine oxidase transcription after myocardial ischemia/reperfusion

Nuclear myosin regulates gene transcription and this novel function might be modulated through phosphorylation of the myosin regulatory light chain (p-MLC20). Nonmuscle MLC20 (nmMLC20) is also present in the nuclei of cardiomyocytes and a potential nmMLC20 binding sequence has been identified in the promoter of the xanthine oxidase (XO) gene. Thus, we investigated its function in the regulation of XO transcription after myocardial ischemia/reperfusion (IR). In a rat model of myocardial IR and a cardiomyocyte model of hypoxia/reoxygenation (HR) injury, the cardiac or cell injury, myosin light chain kinase (MLCK) content, XO expression and activity, XO-derived products, and level of nuclear p-nmMLC20 were detected. Coimmunoprecipitation (co-IP), chromatin immunoprecipitation, DNA pull-down, and luciferase reporter gene assays were used to decipher the molecular mechanisms through which nmMLC20 promotes XO expression. IR or HR treatment dramatically elevated nuclear p-nmMLC20 level, accompanied by increased XO expression, activity, and products (H2O2 and uric acid), as well as the IR or HR injury; these effects were ameliorated by inhibition of MLCK or knockdown of nmMLC20. Our findings from these experiments demonstrated that nuclear p-nmMLC20 binds to the consensus sequence GTCGCC in the XO gene promoter, interacts with RNA polymerase II and transcription factor IIB to form a transcription preinitiation complex, and hence activates XO gene transcription. These results suggest that nuclear p-nmMLC20 plays an important role in IR/HR injury by transcriptionally upregulating XO gene expression to increase oxidative stress in myocardium. Our findings demonstrate nuclear nmMLC20 as a potential new therapeutic target to combat cardiac IR injury.

Gene and MicroRNA Expression Responses to Exercise; Relationship with Insulin Sensitivity

Background

Healthy individuals on the lower end of the insulin sensitivity spectrum also have a reduced gene expression response to exercise for specific genes. The goal of this study was to determine the relationship between insulin sensitivity and exercise-induced gene expression in an unbiased, global manner.

Methods and Findings

Euglycemic clamps were used to measure insulin sensitivity and muscle biopsies were done at rest and 30 minutes after a single acute exercise bout in 14 healthy participants. Changes in mRNA expression were assessed using microarrays, and miRNA analysis was performed in a subset of 6 of the participants using sequencing techniques. Following exercise, 215 mRNAs were changed at the probe level (Bonferroni-corrected P<0.00000115). Pathway and Gene Ontology analysis showed enrichment in MAP kinase signaling, transcriptional regulation and DNA binding. Changes in several transcription factor mRNAs were correlated with insulin sensitivity, including MYC, r=0.71; SNF1LK, r=0.69; and ATF3, r= 0.61 (5 corrected for false discovery rate). Enrichment in the 5’-UTRs of exercise-responsive genes suggested regulation by common transcription factors, especially EGR1. miRNA species of interest that changed after exercise included miR-378, which is located in an intron of the PPARGC1B gene.

Conclusions

These results indicate that transcription factor gene expression responses to exercise depend highly on insulin sensitivity in healthy people. The overall pattern suggests a coordinated cycle by which exercise and insulin sensitivity regulate gene expression in muscle.