TransFind--predicting transcriptional regulators for gene sets

Cell-state dependent regulation of PPAR γ signaling by ZBTB9 in adipocytes

Adipocytes play a critical role in metabolic homeostasis. Peroxisome proliferator-activated receptor- γ (PPAR γ ) is a nuclear hormone receptor that is a master regulator of adipocyte differentiation and function. ZBTB9 was predicted to interact with PPAR γ based on large-scale protein interaction experiments. In addition, GWAS studies in the type 2 diabetes (T2D) Knowledge Portal revealed associations between Z btb9 and both BMI and T2D risk. Here we show that ZBTB9 positively regulates PPAR γ activity in mature adipocytes. Surprisingly Z btb9 knockdown (KD) also increased adipogenesis in 3T3-L1 cells and human preadipocytes. E2F activity was increased and E2F downstream target genes were upregulated in Zbtb9 -KD preadipocytes. Accordingly, RB phosphorylation, which regulates E2F activity, was enhanced in Zbtb9 -KD preadipocytes. Critically, an E2F1 inhibitor blocked the effects of Zbtb9 deficiency on adipogenic gene expression and lipid accumulation. Collectively, these results demonstrate that Zbtb9 inhibits adipogenesis as a negative regulator of Pparg expression via altered RB-E2F1 signaling. Our findings reveal complex cell-state dependent roles of ZBTB9 in adipocytes, identifying a new molecule that regulates adipogenesis and adipocyte biology as both a positive and negative regulator of PPAR γ signaling depending on the cellular context, and thus may be important in the pathogenesis and treatment of obesity and T2D.

Interferon regulatory factor 7 mediates obesity-associated MCP-1 transcription

Hypertrophy, associated with adipocyte dysfunction, causes increased pro-inflammatory adipokine, and abnormal glucose and lipid metabolism, leading to insulin resistance and obesity-related-health problems. By combining DNA microarray and genomic data analyses to predict DNA binding motifs, we identified the transcription factor Interferon Regulatory Factor 7 (IRF7) as a possible regulator of genes related to adipocyte hypertrophy. To investigate the role of IRF7 in adipocytes, we examined gene expression patterns in 3T3-L1 cells infected with a retrovirus carrying the IRF7 gene and found that enforced IRF7 expression induced the expression of monocyte chemoattractant protein-1 (MCP-1), a key initial adipokine in the chronic inflammation of obesity. CRISPR/Cas9 mediated-suppression of IRF7 significantly reduced MCP-1 mRNA. Luciferase assays, chromatin immunoprecipitation PCR analysis and gel shift assay showed that IRF7 transactivates the MCP-1 gene by binding to its proximal Interferon Stimulation Response Element (ISRE), a putative IRF7 binding motif. IRF7 knockout mice exhibited lower expression of MCP-1 in epidydimal white adipose tissue under high-fat feeding conditions, suggesting the transcription factor is physiologically important for inducing MCP-1. Taken together, our results suggest that IRF7 transactivates MCP-1 mRNA in adipocytes, and it may be involved in the adipose tissue inflammation associated with obesity.

Transcriptome sequencing profiles of cervical cancer tissues and SiHa cells

High-risk human papillomavirus (HPV) is a causal factor for cervical cancer, of which HPV16 is the predominant genotype, but the detailed mechanism remains to be elucidated. In this study, we performed transcriptome sequencing in cervical cancer tissues with HPV16-positive and normal tissues with HPV16-negative, and SiHa cells with or without HPV16 E6/E7 knockdown, and identified 140 differential expressed genes (DEGs) in two data sets. We carried out a series of bioinformatic analyses to learn more about the 140 DEGs, and found that 140 DEGs were mostly enriched in cell cycle and DNA repair through Kyoto Encyclopedia of Genes and Genomes pathway enrichment, Gene Ontology annotation, and gene set enrichment analysis. A total of 20 genes including RMI1, MKI67, FANCB, KIF14, CENPI, RACGAP1, EXO1, KIF4A, FOXM1, C19orf57, PSRC1, NUSAP1, CIT, NDC80, MCM7, GINS2, MCM6, ORC1, TLX2, and UHRF1 were screened by co-expression analysis; of those, the expressions of 6 (CENPI, FANCB, KIF14, ORC1, RACGAP1, and RMI1) were verified by qRT-PCR. Further, we found that E2F family, NF-Y, AhR:Arnt, and KROX family may be involved in modulating DEGs by TransFind prediction. TF2DNA database and co-expression analysis suggested that 12 TFs (ZNF367, TLX2, DEPDC1B, E2F8, ZNF541, EGR2, ZMAT3, HES6, CEBPA, MYBL2, FOXM1, and RAD51) were upstream modulators of DEGs. Our findings may provide a new understanding for effects of HPV oncogenes in the maintenance of cancerous state at the transcriptional level.

Neurotrophic Factors and Their Receptors Are Altered by the Mere Partial IGF-1 Deficiency

Neurotrophic factors (NTFs) are a relevant group of secreted proteins that modulate growth, differentiation, repair, and survival of neurons, playing a role in the maintenance of the synaptic unions, dendrites, and axons and also being crucial for peripheral nervous system development and regulating plasticity in the adult central nervous system. On the other hand, insulin-like growth factor 1 (IGF-1) has been ascertained multiple beneficial actions in the brain: neuro-development, -protection, -genesis and plasticity. To further investigate the possible mechanisms underlying IGF-1 deficiency in the establishment of neurological disease, microarray and reverse transcription polymerase chain reaction gene expression analyses coupled with in silico processing were performed in an experimental model of partial IGF-1 deficiency. Results show that the mere IGF-1 deficiency seems to be responsible for an altered expression of genes coding for neurotrophic factors (particularly ciliary neurotrophic factor and mesencephalic astrocyte-derived neurotrophic factor), their receptors and signaling pathways (specially RET). The presented findings support that IGF-1 deficiency might be involved in the establishment and progression of neurodegenerative disorders.

Transcriptomic and proteomic analyses reveal new insights into the regulation of immune pathways during adenovirus type 2 infection

Background

Human adenovirus (Ad) infection leads to the changes of host cell gene expression and biosynthetic processes. Transcriptomics in adenovirus type 2 (Ad2)-infected lung fibroblasts (IMR-90) cells has previously been studied using RNA sequencing. However, this study included only two time points (12 and 24 hpi) using constrained 76 bp long sequencing reads. Therefore, a more detailed study of transcription at different phases of infection using an up-graded sequencing technique is recalled. Furthermore, the correlation between transcription and protein expression needs to be addressed.

Results

In total, 3556 unique cellular genes were identified as differentially expressed at the transcriptional level with more than 2-fold changes in Ad2-infected cells as compared to non-infected cells by using paired-end sequencing. Based on the kinetics of the gene expression changes at different times after infection, these RNAs fell into 20 clusters. Among them, cellular genes involved in immune response were highly up-regulated in the early phase before becoming down-regulated in the late phase. Comparison of differentially expressed genes at transcriptional and posttranscriptional levels revealed low correlation. Particularly genes involved in cellular immune pathways showed a negative correlation. Here, we highlight the genes which expose inconsistent expression profiles with an emphasis on key factors in cellular immune pathways including NFκB, JAK/STAT, caspases and MAVS. Different from their transcriptional profiles with up- and down-regulation in the early and late phase, respectively, these proteins were up-regulated in the early phase and were sustained in the late phase. A surprising finding was that the target genes of the sustained activators failed to show response.

Conclusion

There were features common to genes which play important roles in cellular immune pathways. Their expression was stimulated at both RNA and protein levels during the early phase. In the late phase however, their transcription was suppressed while protein levels remained stable. These results indicate that Ad2 and the host cell use different strategies to regulate cellular immune pathways. A control mechanism at the post-translational level must thus exist which is under the control of Ad2.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1375-5) contains supplementary material, which is available to authorized users.

From gene networks to drugs: systems pharmacology approaches for AUD

The alcohol research field has amassed an impressive number of gene expression datasets spanning key brain areas for addiction, species (humans as well as multiple animal models), and stages in the addiction cycle (binge/intoxication, withdrawal/negative effect, and preoccupation/anticipation). These data have improved our understanding of the molecular adaptations that eventually lead to dysregulation of brain function and the chronic, relapsing disorder of addiction. Identification of new medications to treat alcohol use disorder (AUD) will likely benefit from the integration of genetic, genomic, and behavioral information included in these important datasets. Systems pharmacology considers drug effects as the outcome of the complex network of interactions a drug has rather than a single drug-molecule interaction. Computational strategies based on this principle that integrate gene expression signatures of pharmaceuticals and disease states have shown promise for identifying treatments that ameliorate disease symptoms (called in silico gene mapping or connectivity mapping). In this review, we suggest that gene expression profiling for in silico mapping is critical to improve drug repurposing and discovery for AUD and other psychiatric illnesses. We highlight studies that successfully apply gene mapping computational approaches to identify or repurpose pharmaceutical treatments for psychiatric illnesses. Furthermore, we address important challenges that must be overcome to maximize the potential of these strategies to translate to the clinic and improve healthcare outcomes.

Network analysis reveals a causal role of mitochondrial gene activity in atherosclerotic lesion formation

BACKGROUND AND AIMS

Mitochondrial damage and augmented production of reactive oxygen species (ROS) may represent an intermediate step by which hypercholesterolemia exacerbates atherosclerotic lesion formation.

METHODS

To test this hypothesis, in mice with severe but genetically reversible hypercholesterolemia (i.e. the so called Reversa mouse model), we performed time-resolved analyses of mitochondrial transcriptome in the aortic arch employing a systems-level network approach.

RESULTS

During hypercholesterolemia, we observed a massive down-regulation (>28%) of mitochondrial genes, specifically at the time of rapid atherosclerotic lesion expansion and foam cell formation, i.e. between 30 and 40 weeks of age. Both phenomena - down-regulation of mitochondrial genes and lesion expansion - were largely reversible by genetically lowering plasma cholesterol (by >80%, from 427 to 54 ± 31 mg/L) at 30 weeks. Co-expression network analysis revealed that both mitochondrial signature genes were highly connected in two modules, negatively correlating with lesion size and supported as causal for coronary artery disease (CAD) in humans, as expression-associated single nucleotide polymorphisms (eSNPs) representing their genes overlapped markedly with established disease risk loci. Within these modules, we identified the transcription factor estrogen related receptor (ERR)-α and its co-factors PGC1-α and -β, i.e. two members of the peroxisome proliferator-activated receptor γ co-activator 1 family of transcription regulators, as key regulatory genes. Together, these factors are known as major orchestrators of mitochondrial biogenesis and antioxidant responses.

CONCLUSIONS

Using a network approach, we demonstrate how hypercholesterolemia could hamper mitochondrial activity during atherosclerosis progression and pinpoint potential therapeutic targets to counteract these processes.

Old age and the associated impairment of bones' adaptation to loading are associated with transcriptomic changes in cellular metabolism, cell-matrix interactions and the cell cycle

In old animals, bone's ability to adapt its mass and architecture to functional load-bearing requirements is diminished, resulting in bone loss characteristic of osteoporosis. Here we investigate transcriptomic changes associated with this impaired adaptive response. Young adult (19-week-old) and aged (19-month-old) female mice were subjected to unilateral axial tibial loading and their cortical shells harvested for microarray analysis between 1h and 24h following loading (36 mice per age group, 6 mice per loading group at 6 time points). In non-loaded aged bones, down-regulated genes are enriched for MAPK, Wnt and cell cycle components, including E2F1. E2F1 is the transcription factor most closely associated with genes down-regulated by ageing and is down-regulated at the protein level in osteocytes. Genes up-regulated in aged bone are enriched for carbohydrate metabolism, TNFα and TGFβ superfamily components. Loading stimulates rapid and sustained transcriptional responses in both age groups. However, genes related to proliferation are predominantly up-regulated in the young and down-regulated in the aged following loading, whereas those implicated in bioenergetics are down-regulated in the young and up-regulated in the aged. Networks of inter-related transcription factors regulated by E2F1 are loading-responsive in both age groups. Loading regulates genes involved in similar signalling cascades in both age groups, but these responses are more sustained in the young than aged. From this we conclude that cells in aged bone retain the capability to sense and transduce loading-related stimuli, but their ability to translate acute responses into functionally relevant outcomes is diminished.

Genome-wide Chromatin Profiling of Legionella pneumophila-Infected Human Macrophages Reveals Activation of the Probacterial Host Factor TNFAIP2

BACKGROUND

Legionella pneumophila is a causative agent of severe pneumonia. Infection leads to a broad host cell response, as evident, for example, on the transcriptional level. Chromatin modifications, which control gene expression, play a central role in the transcriptional response to L. pneumophila

METHODS

We infected human-blood-derived macrophages (BDMs) with L. pneumophila and used chromatin immunoprecipitation followed by sequencing to screen for gene promoters with the activating histone 4 acetylation mark.

RESULTS

We found the promoter of tumor necrosis factor α-induced protein 2 (TNFAIP2) to be acetylated at histone H4. This factor has not been characterized in the pathology of L. pneumophila TNFAIP2 messenger RNA and protein were upregulated in response to L. pneumophila infection of human-BDMs and human alveolar epithelial (A549) cells. We showed that L. pneumophila-induced TNFAIP2 expression is dependent on the NF-κB transcription factor. Importantly, knock down of TNFAIP2 led to reduced intracellular replication of L. pneumophila Corby in A549 cells.

CONCLUSIONS

Taken together, genome-wide chromatin analysis of L. pneumophila-infected macrophages demonstrated induction of TNFAIP2, a NF-κB-dependent factor relevant for bacterial replication.

Actin remodeling confers BRAF inhibitor resistance to melanoma cells through YAP/TAZ activation

The activation of transcriptional coactivators YAP and its paralog TAZ has been shown to promote resistance to anti-cancer therapies. YAP/TAZ activity is tightly coupled to actin cytoskeleton architecture. However, the influence of actin remodeling on cancer drug resistance remains largely unexplored. Here, we report a pivotal role of actin remodeling in YAP/TAZ-dependent BRAF inhibitor resistance in BRAF V600E mutant melanoma cells. Melanoma cells resistant to the BRAF inhibitor PLX4032 exhibit an increase in actin stress fiber formation, which appears to promote the nuclear accumulation of YAP/TAZ. Knockdown of YAP/TAZ reduces the viability of resistant melanoma cells, whereas overexpression of constitutively active YAP induces resistance. Moreover, inhibition of actin polymerization and actomyosin tension in melanoma cells suppresses both YAP/TAZ activation and PLX4032 resistance. Our siRNA library screening identifies actin dynamics regulator TESK1 as a novel vulnerable point of the YAP/TAZ-dependent resistance pathway. These results suggest that inhibition of actin remodeling is a potential strategy to suppress resistance in BRAF inhibitor therapies.

Global Analysis of Mouse Polyomavirus Infection Reveals Dynamic Regulation of Viral and Host Gene Expression and Promiscuous Viral RNA Editing

Mouse polyomavirus (MPyV) lytically infects mouse cells, transforms rat cells in culture, and is highly oncogenic in rodents. We have used deep sequencing to follow MPyV infection of mouse NIH3T6 cells at various times after infection and analyzed both the viral and cellular transcriptomes. Alignment of sequencing reads to the viral genome illustrated the transcriptional profile of the early-to-late switch with both early-strand and late-strand RNAs being transcribed at all time points. A number of novel insights into viral gene expression emerged from these studies, including the demonstration of widespread RNA editing of viral transcripts at late times in infection. By late times in infection, 359 host genes were seen to be significantly upregulated and 857 were downregulated. Gene ontology analysis indicated transcripts involved in translation, metabolism, RNA processing, DNA methylation, and protein turnover were upregulated while transcripts involved in extracellular adhesion, cytoskeleton, zinc finger binding, SH3 domain, and GTPase activation were downregulated. The levels of a number of long noncoding RNAs were also altered. The long noncoding RNA MALAT1, which is involved in splicing speckles and used as a marker in many late-stage cancers, was noticeably downregulated, while several other abundant noncoding RNAs were strongly upregulated. We discuss these results in light of what is currently known about the MPyV life cycle and its effects on host cell growth and metabolism.

Mouse polyomavirus (MPyV) is a small 5.3kb circular double-stranded DNA virus capable of causing tumors in a variety of tissues in immunocompromised mice. It has been a subject of study for over 60 years, yielding insights into a number of processes including tumorigenesis, cell cycle signaling, and transformation. This study serves to provide a global view of the MPyV infection by utilizing Illumina sequencing to observe changes in total RNA from both the virus and the host cell as well as applying new methods to more directly confirm the extent of A-to-I editing of viral RNA by host ADAR enzymes. This allows for a simultaneous observation of both host and viral transcriptional changes that occur as a result of early gene expression and the viral switch from early to late genes that occurs coincident with the initiation of DNA replication.

Increased expression of interferon signaling genes in the bone marrow microenvironment of myelodysplastic syndromes

Introduction

The bone marrow (BM) microenvironment plays an important role in the pathogenesis of myelodysplastic syndromes (MDS) through a reciprocal interaction with resident BM hematopoietic cells. We investigated the differences between BM mesenchymal stromal cells (MSCs) in MDS and normal individuals and identified genes involved in such differences.

Materials and Methods

BM-derived MSCs from 7 MDS patients (3 RCMD, 3 RAEB-1, and 1 RAEB-2) and 7 controls were cultured. Global gene expression was analyzed using a microarray.

Result

We found 314 differentially expressed genes (DEGs) in RCMD vs. control, 68 in RAEB vs. control, and 51 in RAEB vs. RCMD. All comparisons were clearly separated from one another by hierarchical clustering. The overall similarity between differential expression signatures from the RCMD vs. control comparison and the RAEB vs. control comparison was highly significant (p = 0), which indicates a common transcriptomic response in these two MDS subtypes. RCMD and RAEB simultaneously showed an up-regulation of interferon alpha/beta signaling and the ISG15 antiviral mechanism, and a significant fraction of the RAEB vs. control DEGs were also putative targets of transcription factors IRF and ICSBP. Pathways that involved RNA polymerases I and III and mitochondrial transcription were down-regulated in RAEB compared to RCMD.

Conclusion

Gene expression in the MDS BM microenvironment was different from that in normal BM and exhibited altered expression according to disease progression. The present study provides genetic evidence that inflammation and immune dysregulation responses that involve the interferon signaling pathway in the BM microenvironment are associated with MDS pathogenesis, which suggests BM MSCs as a possible therapeutic target in MDS.

Transcriptomic portrait of human Mesenchymal Stromal/Stem Cells isolated from bone marrow and placenta

Background

Human Mesenchymal Stromal/Stem Cells (MSCs) are adult multipotent cells that behave in a highly plastic manner, inhabiting the stroma of several tissues. The potential utility of MSCs is nowadays strongly investigated in the field of regenerative medicine and cell therapy, although many questions about their molecular identity remain uncertain.

Results

MSC primary cultures from human bone marrow (BM) and placenta (PL) were derived and verified by their immunophenotype standard pattern and trilineage differentiation potential. Then, a broad characterization of the transcriptome of these MSCs was performed using RNA deep sequencing (RNA-Seq). Quantitative analysis of these data rendered an extensive expression footprint that includes 5,271 protein-coding genes. Flow cytometry assays of canonical MSC CD-markers were congruent with their expression levels detected by the RNA-Seq. Expression of other recently proposed MSC markers (CD146, Nestin and CD271) was tested in the placenta samples, finding only CD146 and Nestin. Functional analysis revealed enrichment in stem cell related genes and mesenchymal regulatory transcription factors (TFs). Analysis of TF binding sites (TFBSs) identified 11 meta-regulators, including factors KLF4 and MYC among them. Epigenetically, hypomethylated promoter patterns supported the active expression of the MSC TFs found. An interaction network of these TFs was built to show up their links and relations. Assessment of dissimilarities between cell origins (BM versus PL) disclosed two hundred differentially expressed genes enrolled in microenvironment processes related to the cellular niche, as regulation of bone formation and blood vessel morphogenesis for the case of BM-MSCs. By contrast genes overexpressed in PL-MSCs showed functional enrichment on mitosis, negative regulation of cell-death and embryonic morphogenesis that supported the higher growth rates observed in the cultures of these fetal cells and their closer links with development processes.

Conclusions

The results present a transcriptomic portrait of the human MSCs isolated from bone marrow and placenta. The data are released as a cell-specific resource, providing a comprehensive expression footprint of the MSCs useful to better understand their cellular and molecular biology and for further investigations on the isolation and biomedical use of these multipotent cells.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-910) contains supplementary material, which is available to authorized users.

NFkB disrupts tissue polarity in 3D by preventing integration of microenvironmental signals

The microenvironment of cells controls their phenotype, and thereby the architecture of the emerging multicellular structure or tissue. We have reported more than a dozen microenvironmental factors whose signaling must be integrated in order to effect an organized, functional tissue morphology. However, the factors that prevent integration of signaling pathways that merge form and function are still largely unknown. We have identified nuclear factor kappa B (NFkB) as a transcriptional regulator that disrupts important microenvironmental cues necessary for tissue organization. We compared the gene expression of organized and disorganized epithelial cells of the HMT-3522 breast cancer progression series: the non-malignant S1 cells that form polarized spheres (‘acini’), the malignant T4-2 cells that form large tumor-like clusters, and the ‘phenotypically reverted’ T4-2 cells that polarize as a result of correction of the microenvironmental signaling. We identified 180 genes that display an increased expression in disorganized compared to polarized structures. Network, GSEA and transcription factor binding site analyses suggested that NFkB is a common activator for the 180 genes. NFkB was found to be activated in disorganized breast cancer cells, and inhibition of microenvironmental signaling via EGFR, beta1 integrin, MMPs, or their downstream signals suppressed its activation. The postulated role of NFkB was experimentally verified: Blocking the NFkB pathway with a specific chemical inhibitor or shRNA induced polarization and inhibited invasion of breast cancer cells in 3D cultures. These results may explain why NFkB holds promise as a target for therapeutic intervention: Its inhibition can reverse the oncogenic signaling involved in breast cancer progression and integrate the essential microenvironmental control of tissue architecture.

KRAS and YAP1 converge to regulate EMT and tumor survival

Cancer cells that express oncogenic alleles of RAS typically require sustained expression of the mutant allele for survival, but the molecular basis of this oncogene dependency remains incompletely understood. To identify genes that can functionally substitute for oncogenic RAS, we systematically expressed 15,294 open reading frames in a human KRAS-dependent colon cancer cell line engineered to express an inducible KRAS-specific shRNA. We found 147 genes that promoted survival upon KRAS suppression. In particular, the transcriptional coactivator YAP1 rescued cell viability in KRAS-dependent cells upon suppression of KRAS and was required for KRAS-induced cell transformation. Acquired resistance to Kras suppression in a Kras-driven murine lung cancer model also involved increased YAP1 signaling. KRAS and YAP1 converge on the transcription factor FOS and activate a transcriptional program involved in regulating the epithelial-mesenchymal transition (EMT). Together, these findings implicate transcriptional regulation of EMT by YAP1 as a significant component of oncogenic RAS signaling.

Yap1 activation enables bypass of oncogenic Kras addiction in pancreatic cancer

Activating mutations in KRAS are among the most frequent events in diverse human carcinomas and are particularly prominent in human pancreatic ductal adenocarcinoma (PDAC). An inducible Kras(G12D)-driven mouse model of PDAC has established a critical role for sustained Kras(G12D) expression in tumor maintenance, providing a model to determine the potential for and the underlying mechanisms of Kras(G12D)-independent PDAC recurrence. Here, we show that some tumors undergo spontaneous relapse and are devoid of Kras(G12D) expression and downstream canonical MAPK signaling and instead acquire amplification and overexpression of the transcriptional coactivator Yap1. Functional studies established the role of Yap1 and the transcriptional factor Tead2 in driving Kras(G12D)-independent tumor maintenance. The Yap1/Tead2 complex acts cooperatively with E2F transcription factors to activate a cell cycle and DNA replication program. Our studies, along with corroborating evidence from human PDAC models, portend a novel mechanism of escape from oncogenic Kras addiction in PDAC.

In silico discovery of novel transcription factors regulated by mTOR-pathway activities

The mammalian target of rapamycine (mTOR) pathway is a key regulator of cellular growth, development, and ageing, and unraveling its control is essential for understanding life and death of biological organisms. A motif-discovery workbench including nine tools was used to identify transcription factors involved in five basic (Insulin, MAPK, VEGF, Hypoxia, and mTOR core) activities of the mTOR pathway. Discovered transcription factors are classified as “process-specific” or “pathway-ubiquitous” with highlights toward their regulating/regulated activities within the mTOR pathway. Our transcription regulation results will facilitate further research on investigating the control mechanism in mTOR pathway.

Transcriptional regulation induced by cAMP elevation in mouse Schwann cells

In peripheral nerves, Schwann cell development is regulated by a variety of signals. Some of the aspects of Schwann cell differentiation can be reproduced in vitro in response to forskolin, an adenylyl cyclase activator elevating intracellular cAMP levels. Herein, the effect of forskolin treatment was investigated by a comprehensive genome-wide expression study on primary mouse Schwann cell cultures. Additional to myelin-related genes, many so far unconsidered genes were ascertained to be modulated by forskolin. One of the strongest differentially regulated gene transcripts was the transcription factor Olig1 (oligodendrocyte transcription factor 1), whose mRNA expression levels were reduced in treated Schwann cells. Olig1 protein was localized in myelinating and nonmyelinating Schwann cells within the sciatic nerve as well as in primary Schwann cells, proposing it as a novel transcription factor of the Schwann cell lineage. Data analysis further revealed that a number of differentially expressed genes in forskolin-treated Schwann cells were associated with the ECM (extracellular matrix), underlining its importance during Schwann cell differentiation in vitro. Comparison of samples derived from postnatal sciatic nerves and from both treated and untreated Schwann cell cultures showed considerable differences in gene expression between in vivo and in vitro, allowing us to separate Schwann cell autonomous from tissue-related changes. The whole data set of the cell culture microarray study is provided to offer an interactive search tool for genes of interest.

Integrative radiation systems biology

Maximisation of the ratio of normal tissue preservation and tumour cell reduction is the main concept of radiotherapy alone or combined with chemo-, immuno- or biologically targeted therapy. The foremost parameter influencing this ratio is radiation sensitivity and its modulation towards a more efficient killing of tumour cells and a better preservation of normal tissue at the same time is the overall aim of modern therapy schemas. Nevertheless, this requires a deep understanding of the molecular mechanisms of radiation sensitivity in order to identify its key players as potential therapeutic targets. Moreover, the success of conventional approaches that tried to statistically associate altered radiation sensitivity with any molecular phenotype such as gene expression proofed to be somewhat limited since the number of clinically used targets is rather sparse. However, currently a paradigm shift is taking place from pure frequentistic association analysis to the rather holistic systems biology approach that seeks to mathematically model the system to be investigated and to allow the prediction of an altered phenotype as the function of one single or a signature of biomarkers. Integrative systems biology also considers the data from different molecular levels such as the genome, transcriptome or proteome in order to partially or fully comprehend the causal chain of molecular mechanisms. An example for the application of this concept currently carried out at the Clinical Cooperation Group "Personalized Radiotherapy in Head and Neck Cancer" of the Helmholtz-Zentrum München and the LMU Munich is described. This review article strives for providing a compact overview on the state of the art of systems biology, its actual challenges, potential applications, chances and limitations in radiation oncology research working towards improved personalised therapy concepts using this relatively new methodology.

Gene expression profile and functional analysis of Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is a complex neurodegenerative disorder with largely unknown genetic mechanisms. Identifying altered neuronal gene expression in brain regions differentially affected by AD may provide the diagnostic or therapeutic targets of AD.

METHODS

The gene expression profile of AD was analyzed with bioinformatics. Function analysis was performed with Database for Annotation, Visualization and Integrated Discovery (DAVID), and TransFind was used to predict the possible transcriptional regulators in AD. Finally, connectivity map (cMap) database was used to explore small molecules targeted for AD.

RESULTS

The AD gene signatures associated with 6 different brain regions were identified. Functional analysis revealed that biological processes involved with metabolism, protein ubiquitination, and vasculature development were found dysregulated, and synaptic signaling pathways were found perturbated in AD. The WT1 was identified as an important transcriptional regulator in AD, and cMap database predicted that small molecules, such as histone deacetylase (HDAC) inhibitor, may be candidate drugs in the treatment of AD.

CONCLUSION

According to our in silico analysis, Wilms' tumor suppressor may play regulatory roles in AD development and progress. The HDAC inhibitor could possibly be used to treat AD.

Genome-wide profiling of the activity-dependent hippocampal transcriptome

Activity-dependent gene expression is central for sculpting neuronal connectivity in the brain. Despite the importance for synaptic plasticity, a comprehensive analysis of the temporal changes in the transcriptomic response to neuronal activity is lacking. In a genome wide survey we identified genes that were induced at 1, 4, 8, or 24 hours following neuronal activity in the hippocampus. According to their distinct expression kinetics we assigned these genes to five clusters, each containing approximately 200 genes. Using in situ hybridizations the regulated expression of 24 genes was validated. Apart from known activity-dependent genes our study reveals a large number of unknown induced genes with distinct expression kinetics. Among these we identified several genes with complex temporal expression patterns. Furthermore, our study provides examples for activity-induced exon switching in the coding region of genes and activity-induced alternative splicing of the 3'-UTR. One example is Zwint. In contrast to the constitutively expressed variant, the induced Zwint transcript harbors multiple regulatory elements in the 3'-UTR. Taken together, our study provides a comprehensive analysis of the transcriptomic response to neuronal activity and sheds new light on expression kinetics and alternative splicing events.

LASAGNA: a novel algorithm for transcription factor binding site alignment

BACKGROUND

Scientists routinely scan DNA sequences for transcription factor (TF) binding sites (TFBSs). Most of the available tools rely on position-specific scoring matrices (PSSMs) constructed from aligned binding sites. Because of the resolutions of assays used to obtain TFBSs, databases such as TRANSFAC, ORegAnno and PAZAR store unaligned variable-length DNA segments containing binding sites of a TF. These DNA segments need to be aligned to build a PSSM. While the TRANSFAC database provides scoring matrices for TFs, nearly 78% of the TFs in the public release do not have matrices available. As work on TFBS alignment algorithms has been limited, it is highly desirable to have an alignment algorithm tailored to TFBSs.

RESULTS

We designed a novel algorithm named LASAGNA, which is aware of the lengths of input TFBSs and utilizes position dependence. Results on 189 TFs of 5 species in the TRANSFAC database showed that our method significantly outperformed ClustalW2 and MEME. We further compared a PSSM method dependent on LASAGNA to an alignment-free TFBS search method. Results on 89 TFs whose binding sites can be located in genomes showed that our method is significantly more precise at fixed recall rates. Finally, we described LASAGNA-ChIP, a more sophisticated version for ChIP (Chromatin immunoprecipitation) experiments. Under the one-per-sequence model, it showed comparable performance with MEME in discovering motifs in ChIP-seq peak sequences.

CONCLUSIONS

We conclude that the LASAGNA algorithm is simple and effective in aligning variable-length binding sites. It has been integrated into a user-friendly webtool for TFBS search and visualization called LASAGNA-Search. The tool currently stores precomputed PSSM models for 189 TFs and 133 TFs built from TFBSs in the TRANSFAC Public database (release 7.0) and the ORegAnno database (08Nov10 dump), respectively. The webtool is available at http://biogrid.engr.uconn.edu/lasagna_search/.

Phylogenetic footprinting: a boost for microbial regulatory genomics

Phylogenetic footprinting is a method for the discovery of regulatory elements in a set of homologous regulatory regions, usually collected from multiple species. It does so by identifying the best conserved motifs in those homologous regions. There are two popular sets of methods-alignment-based and motif-based, which are generally employed for phylogenetic methods. However, serious efforts have lacked to develop a tool exclusively for phylogenetic footprinting, based on either of these methods. Nevertheless, a number of software and tools exist that can be applied for prediction of phylogenetic footprinting with variable degree of success. The output from these tools may get affected by a number of factors associated with current state of knowledge, techniques and other resources available. We here present a critical apprehension of various phylogenetic approaches with reference to prokaryotes outlining the available resources and also discussing various factors affecting footprinting in order to make a clear idea about the proper use of this approach on prokaryotes.

Genome-wide profiling of methylation identifies novel targets with aberrant hypermethylation and reduced expression in low-risk myelodysplastic syndromes

Gene expression profiling signatures may be used to classify the subtypes of Myelodysplastic syndrome (MDS) patients. However, there are few reports on the global methylation status in MDS. The integration of genome-wide epigenetic regulatory marks with gene expression levels would provide additional information regarding the biological differences between MDS and healthy controls. Gene expression and methylation status were measured using high-density microarrays. A total of 552 differentially methylated CpG loci were identified as being present in low-risk MDS; hypermethylated genes were more frequent than hypomethylated genes. In addition, mRNA expression profiling identified 1005 genes that significantly differed between low-risk MDS and the control group. Integrative analysis of the epigenetic and expression profiles revealed that 66.7% of the hypermethylated genes were underexpressed in low-risk MDS cases. Gene network analysis revealed molecular mechanisms associated with the low-risk MDS group, including altered apoptosis pathways. The two key apoptotic genes BCL2 and ETS1 were identified as silenced genes. In addition, the immune response and micro RNA biogenesis were affected by the hypermethylation and underexpression of IL27RA and DICER1. Our integrative analysis revealed that aberrant epigenetic regulation is a hallmark of low-risk MDS patients and could have a central role in these diseases.

Sustained liver regeneration after portal vein embolization --a human molecular pilot study

BACKGROUND

Portal vein embolization is a treatment option to achieve a sufficient future remnant liver volume for patients with central liver tumours requiring an extended resection with an extensive parenchymal loss. However, molecular mechanisms of this intervention are up to now poorly understood. The objective of this prospective pilot study was the characterization of molecular events leading to late hypertrophy of the non-embolized liver tissue in the human liver.

METHODS

Liver tissue of ten patients was collected before and intraoperatively more than one month after embolization. Investigation of molecular features was performed by pangenomic chips, polymerase chain reaction, immunostaining of proliferation marker Ki-67 and immunofluorescence measurements.

RESULTS

Significantly elevated genes hint towards angiogenesis and signalling by insulin-like growth factor and associated binding proteins. Increased transcript levels of activator protein 1 complex members like c-jun were reflecting potential molecular events of liver growth after embolization. Immunofluorescence data confirmed a predominant upregulation of β-catenin and c-jun (p<0.1) supported by Ki-67 (p<0.05) in the non-embolized liver. In silico analysis of transcriptomic dysplasia and hepatocellular carcinoma data showed divergent signatures compared to embolization.

CONCLUSIONS

Our findings indicate a sustained regeneration after portal vein embolization reflected in hyperplasia and angiogenesis in the human liver and provide novel molecular mechanisms of interlobe crosstalk.

Inferring the regulatory network behind a gene expression experiment

Transcription factors (TFs) and miRNAs are the most important dynamic regulators in the control of gene expression in multicellular organisms. These regulatory elements play crucial roles in development, cell cycling and cell signaling, and they have also been associated with many diseases. The Regulatory Network Analysis Tool (RENATO) web server makes the exploration of regulatory networks easy, enabling a better understanding of functional modularity and network integrity under specific perturbations. RENATO is suitable for the analysis of the result of expression profiling experiments. The program analyses lists of genes and search for the regulators compatible with its activation or deactivation. Tests of single enrichment or gene set enrichment allow the selection of the subset of TFs or miRNAs significantly involved in the regulation of the query genes. RENATO also offers an interactive advanced graphical interface that allows exploring the regulatory network found.RENATO is available at: http://renato.bioinfo.cipf.es/.

LOESS correction for length variation in gene set-based genomic sequence analysis

MOTIVATION

Sequence analysis algorithms are often applied to sets of DNA, RNA or protein sequences to identify common or distinguishing features. Controlling for sequence length variation is critical to properly score sequence features and identify true biological signals rather than length-dependent artifacts.

RESULTS

Several cis-regulatory module discovery algorithms exhibit a substantial dependence between DNA sequence score and sequence length. Our newly developed LOESS method is flexible in capturing diverse score-length relationships and is more effective in correcting DNA sequence scores for length-dependent artifacts, compared with four other approaches. Application of this method to genes co-expressed during Drosophila melanogaster embryonic mesoderm development or neural development scored by the Lever motif analysis algorithm resulted in successful recovery of their biologically validated cis-regulatory codes. The LOESS length-correction method is broadly applicable, and may be useful not only for more accurate inference of cis-regulatory codes, but also for detection of other types of patterns in biological sequences.

AVAILABILITY

Source code and compiled code are available from http://thebrain.bwh.harvard.edu/LM_LOESS/

The transcriptome of the adenovirus infected cell

Alternations of cellular gene expression following an adenovirus type 2 infection of human primary cells were studied by using superior sensitive cDNA sequencing. In total, 3791 cellular genes were identified as differentially expressed more than 2-fold. Genes involved in DNA replication, RNA transcription and cell cycle regulation were very abundant among the up-regulated genes. On the other hand, genes involved in various signaling pathways including TGF-β, Rho, G-protein, Map kinase, STAT and NF-κB stood out among the down-regulated genes. Binding sites for E2F, ATF/CREB and AP2 were prevalent in the up-regulated genes, whereas binding sites for SRF and NF-κB were dominant among the down-regulated genes. It is evident that the adenovirus has gained a control of the host cell cycle, growth, immune response and apoptosis at 24 h after infection. However, efforts from host cell to block the cell cycle progression and activate an antiviral response were also observed.

Validation of skeletal muscle cis-regulatory module predictions reveals nucleotide composition bias in functional enhancers

We performed a genome-wide scan for muscle-specific cis-regulatory modules (CRMs) using three computational prediction programs. Based on the predictions, 339 candidate CRMs were tested in cell culture with NIH3T3 fibroblasts and C2C12 myoblasts for capacity to direct selective reporter gene expression to differentiated C2C12 myotubes. A subset of 19 CRMs validated as functional in the assay. The rate of predictive success reveals striking limitations of computational regulatory sequence analysis methods for CRM discovery. Motif-based methods performed no better than predictions based only on sequence conservation. Analysis of the properties of the functional sequences relative to inactive sequences identifies nucleotide sequence composition can be an important characteristic to incorporate in future methods for improved predictive specificity. Muscle-related TFBSs predicted within the functional sequences display greater sequence conservation than non-TFBS flanking regions. Comparison with recent MyoD and histone modification ChIP-Seq data supports the validity of the functional regions.

For efficient identification of genomic sequences responsible for regulating gene expression, a number of computer programs have been developed for automatic annotation of these regulatory regions. We searched for potential regulatory regions responsible for controlling the expression of skeletal muscle-specific genes using these programs, and validated the predictions in a popular cell culture model for muscle. We were able to identify 19 previously uncharacterized regulatory regions for muscle genes. The accuracy of the predictions made by these programs leaves much to be desired, leading us to conclude that other signals in addition to the sequence information will be required to achieve sufficient predictive power for genome annotation. Genomic regions with confirmed regulatory function were compared against non-functional sequences, revealing sequence conservation, composition and chromatin modification properties as important signals in determining regulatory region functionality.

Deregulation upon DNA damage revealed by joint analysis of context-specific perturbation data

Background

Deregulation between two different cell populations manifests itself in changing gene expression patterns and changing regulatory interactions. Accumulating knowledge about biological networks creates an opportunity to study these changes in their cellular context.

Results

We analyze re-wiring of regulatory networks based on cell population-specific perturbation data and knowledge about signaling pathways and their target genes. We quantify deregulation by merging regulatory signal from the two cell populations into one score. This joint approach, called JODA, proves advantageous over separate analysis of the cell populations and analysis without incorporation of knowledge. JODA is implemented and freely available in a Bioconductor package 'joda'.

Conclusions

Using JODA, we show wide-spread re-wiring of gene regulatory networks upon neocarzinostatin-induced DNA damage in Human cells. We recover 645 deregulated genes in thirteen functional clusters performing the rich program of response to damage. We find that the clusters contain many previously characterized neocarzinostatin target genes. We investigate connectivity between those genes, explaining their cooperation in performing the common functions. We review genes with the most extreme deregulation scores, reporting their involvement in response to DNA damage. Finally, we investigate the indirect impact of the ATM pathway on the deregulated genes, and build a hypothetical hierarchy of direct regulation. These results prove that JODA is a step forward to a systems level, mechanistic understanding of changes in gene regulation between different cell populations.

Microarray data analysis and mining tools

Microarrays are one of the latest breakthroughs in experimental molecular biology that allow monitoring the expression levels of tens of thousands of genes simultaneously. Arrays have been applied to studies in gene expression, genome mapping, SNP discrimination, transcription factor activity, toxicity, pathogen identification and many other applications. In this paper we concentrate on discussing various bioinformatics tools used for microarray data mining tasks with its underlying algorithms, web resources and relevant reference. We emphasize this paper mainly for digital biologists to get an aware about the plethora of tools and programs available for microarray data analysis. First, we report the common data mining applications such as selecting differentially expressed genes, clustering, and classification. Next, we focused on gene expression based knowledge discovery studies such as transcription factor binding site analysis, pathway analysis, protein- protein interaction network analysis and gene enrichment analysis.