Petri Pehkonen

PIAS proteins: pleiotropic interactors associated with SUMO.

The interactions and functions of protein inhibitors of activated STAT (PIAS) proteins are not restricted to the signal transducers and activators of transcription (STATs), but PIAS1, -2, -3 and -4 interact with and regulate a variety of distinct proteins, especially transcription factors. Although the majority of PIAS-interacting proteins are prone to modification by small ubiquitin-related modifier (SUMO) proteins and the PIAS proteins have the capacity to promote the modification as RING-type SUMO ligases, they do not function solely as SUMO E3 ligases. Instead, their effects are often independent of their Siz/PIAS (SP)-RING finger, but dependent on their capability to noncovalently interact with SUMOs or DNA through their SUMO-interacting motif and scaffold attachment factor-A/B, acinus and PIAS domain, respectively. Here, we present an overview of the cellular regulation by PIAS proteins and propose that many of their functions are due to their capability to mediate and facilitate SUMO-linked protein assemblies.

Gene expression in the brain and kidney of rainbow trout in response to handling stress

BackgroundMicroarray technologies are rapidly becoming available for new species including teleost fishes. We constructed a rainbow trout cDNA microarray targeted at the identification of genes which are differentially expressed in response to environmental stressors. This platform included clones from normalized and subtracted libraries and genes selected through functional annotation. Present study focused on time-course comparisons of stress responses in the brain and kidney and the identification of a set of genes which are diagnostic for stress response.ResultsFish were stressed with handling and samples were collected 1, 3 and 5 days after the first exposure. Gene expression profiles were analysed in terms of Gene Ontology categories. Stress affected different functional groups of genes in the tissues studied. Mitochondria, extracellular matrix and endopeptidases (especially collagenases) were the major targets in kidney. Stress response in brain was characterized with dramatic temporal alterations. Metal ion binding proteins, glycolytic enzymes and motor proteins were induced transiently, whereas expression of genes involved in stress and immune response, cell proliferation and growth, signal transduction and apoptosis, protein biosynthesis and folding changed in a reciprocal fashion. Despite dramatic difference between tissues and time-points, we were able to identify a group of 48 genes that showed strong correlation of expression profiles (Pearson r > |0.65|) in 35 microarray experiments being regulated by stress. We evaluated performance of the clone sets used for preparation of microarray. Overall, the number of differentially expressed genes was markedly higher in EST than in genes selected through Gene Ontology annotations, however 63% of stress-responsive genes were from this group.Conclusions1. Stress responses in fish brain and kidney are different in function and time-course. 2. Identification of stress-regulated genes provides the possibility for measuring stress responses in various conditions and further search for the functionally related genes.

Nuclear hormone 1α,25-dihydroxyvitamin D3 elicits a genome-wide shift in the locations of VDR chromatin occupancy

A global understanding of the actions of the nuclear hormone 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3) and its vitamin D receptor (VDR) requires a genome-wide analysis of VDR binding sites. In THP-1 human monocytic leukemia cells we identified by ChIP-seq 2340 VDR binding locations, of which 1171 and 520 occurred uniquely with and without 1α,25(OH)2D3 treatment, respectively, while 649 were common. De novo identified direct repeat spaced by 3 nucleotides (DR3)-type response elements (REs) were strongly associated with the ligand-responsiveness of VDR occupation. Only 20% of the VDR peaks diminishing most after ligand treatment have a DR3-type RE, in contrast to 90% for the most growing peaks. Ligand treatment revealed 638 1α,25(OH)2D3 target genes enriched in gene ontology categories associated with immunity and signaling. From the 408 upregulated genes, 72% showed VDR binding within 400 kb of their transcription start sites (TSSs), while this applied only for 43% of the 230 downregulated genes. The VDR loci showed considerable variation in gene regulatory scenarios ranging from a single VDR location near the target gene TSS to very complex clusters of multiple VDR locations and target genes. In conclusion, ligand binding shifts the locations of VDR occupation to DR3-type REs that surround its target genes and occur in a large variety of regulatory constellations.

Identification of gene expression changes in transgenic C. elegans overexpressing human α-synuclein

Alpha-synuclein containing cellular inclusions are a hallmark of Parkinson Disease, Lewy Body Dementia, and Multiple System Atrophy. A genome wide expression screen was performed in C. elegans overexpressing both wild-type and A53T human alpha-synuclein. 433 genes were up- and 67 genes down-regulated by statistical and fold change (> or <2) criteria. Gene ontology (GO) categories within the regulated gene lists indicated over-representation of development and reproduction, mitochondria, catalytic activity, and histone groups. Seven genes (pdr-1, ubc-7, pas-5, pas-7, pbs-4, RPT2, PSMD9) with function in the ubiquitin-proteasome system and 35 mitochondrial function genes were up-regulated. Nine genes that form histones H1, H2B, and H4 were down-regulated. These results demonstrate the effects of alpha-synuclein on proteasome and mitochondrial complex gene expression and provide further support for the role of these complexes in mediating neurotoxicity. The results also indicate an effect on nuclear protein genes that suggests a potential new avenue for investigation.

Theme discovery from gene lists for identification and viewing of multiple functional groups.

BackgroundHigh throughput methods of the genome era produce vast amounts of data in the form of gene lists. These lists are large and difficult to interpret without advanced computational or bioinformatic tools. Most existing methods analyse a gene list as a single entity although it is comprised of multiple gene groups associated with separate biological functions. Therefore it is imperative to define and visualize gene groups with unique functionality within gene lists.ResultsIn order to analyse the functional heterogeneity within a gene list, we have developed a method that clusters genes to groups with homogenous functionalities. The method uses Non-negative Matrix Factorization (NMF) to create several clustering results with varying numbers of clusters. The obtained clustering results are combined into a simple graphical presentation showing the functional groups over-represented in the analyzed gene list. We demonstrate its performance on two data sets and show results that improve upon existing methods. The comparison also shows that our method creates a more simplified view that aids in discovery of biological themes within the list and discards less informative classes from the results.ConclusionThe presented method and associated software are useful for the identification and interpretation of biological functions associated with gene lists and are especially useful for the analysis of large lists.

Genome-wide landscape of liver X receptor chromatin binding and gene regulation in human macrophages

BackgroundThe liver X receptors (LXRs) are oxysterol sensing nuclear receptors with multiple effects on metabolism and immune cells. However, the complete genome-wide cistrome of LXR in cells of human origin has not yet been provided.ResultsWe performed ChIP-seq in phorbol myristate acetate-differentiated THP-1 cells (macrophage-type) after stimulation with the potent synthetic LXR ligand T0901317 (T09). Microarray gene expression analysis was performed in the same cellular model. We identified 1357 genome-wide LXR locations (FDR < 1%), of which 526 were observed after T09 treatment. De novo analysis of LXR binding sequences identified a DR4-type element as the major motif. On mRNA level T09 up-regulated 1258 genes and repressed 455 genes. Our results show that LXR actions are focused on 112 genomic regions that contain up to 11 T09 target genes per region under the control of highly stringent LXR binding sites with individual constellations for each region. We could confirm that LXR controls lipid metabolism and transport and observed a strong association with apoptosis-related functions.ConclusionsThis first report on genome-wide binding of LXR in a human cell line provides new insights into the transcriptional network of LXR and its target genes with their link to physiological processes, such as apoptosis.The gene expression microarray and sequence data have been submitted collectively to the NCBI Gene Expression Omnibus http://www.ncbi.nlm.nih.gov/geo under accession number GSE28319.

Dynamics of 1α,25-dihydroxyvitamin D3-dependent chromatin accessibility of early vitamin D receptor target genes

The signaling cascade of the transcription factor vitamin D receptor (VDR) is triggered by its specific ligand 1α,25-dihydroxyvitamin D3 (1α,25(OH)2D3). In this study we demonstrate that in THP-1 human monocytic leukemia cells 87.4% of the 1034 most prominent genome-wide VDR binding sites co-localize with loci of open chromatin. At 165 of them 1α,25(OH)2D3 strongly increases chromatin accessibility and has at further 217 sites weaker effects. Interestingly, VDR binding sites in 1α,25(OH)2D3-responsive chromatin regions are far more often composed of direct repeats with 3 intervening nucleotides (DR3s) than those in ligand insensitive regions. DR3-containing VDR sites are enriched in the neighborhood of genes that are involved in controling cellular growth, while non-DR3 VDR binding is often found close to genes related to immunity. At the example of six early VDR target genes we show that the slope of their 1α,25(OH)2D3-induced transcription correlates with the basal chromatin accessibility of their major VDR binding regions. However, the chromatin loci controlling these genes are indistinguishable in their VDR association kinetics. Taken together, ligand responsive chromatin loci represent dynamically regulated contact points of VDR with the genome, from where it controls early 1α,25(OH)2D3 target genes.

POXO: a web-enabled tool series to discover transcription factor binding sites

We present POXO, a comprehensive tool series to discover transcription factor binding sites from co-expressed genes (). POXO manages tasks such as functional evaluation and grouping of genes, sequence retrieval, pattern discovery and pattern verification. It also allows users to tailor analytical pipelines from these tools, with single mouse clicks. One typical pipeline of POXO begins by examining the biological functions that a set of co-expressed genes are involved in. In this examination, the functional coherence of the gene set is evaluated and representative functions are associated with the gene set. This examination can also be used to group genes into functionally similar subsets, if several biological processes are affected in the experiment. The next step in the pipeline is then to discover over-represented nucleotide patterns from the upstream sequences of the selected gene sets. This enables to investigate the possibility that the genes are co-regulated by common cis-elements. If over-represented patterns are found, similar ones can then be clustered together and be verified. The performance of POXO is demonstrated by analysing expression data from pathogen treated Arabidopsis thaliana. In this example, POXO detected activated gene sets and suggested transcription factors responsible for their regulation.

Overexpression of SUMO perturbs the growth and development of Caenorhabditis elegans

Small ubiquitin-related modifiers (SUMOs) are important regulator proteins. Caenorhabditis elegans contains a single SUMO ortholog, SMO-1, necessary for the reproduction of C. elegans. In this study, we constructed transgenic C. elegans strains expressing human SUMO-1 under the control of pan-neuronal (aex-3) or pan-muscular (myo-4) promoter and SUMO-2 under the control of myo-4 promoter. Interestingly, muscular overexpression of SUMO-1 or -2 resulted in morphological changes of the posterior part of the nematode. Movement, reproduction and aging of C. elegans were perturbed by the overexpression of SUMO-1 or -2. Genome-wide expression analyses revealed that several genes encoding components of SUMOylation pathway and ubiquitin-proteasome system were upregulated in SUMO-overexpressing nematodes. Since muscular overexpression of SMO-1 also brought up reproductive and mobility perturbations, our results imply that the phenotypes were largely due to an excess of SUMO, suggesting that a tight control of SUMO levels is important for the normal development of multicellular organisms.

The unique characteristics of the placental transcriptome and the hormonal metabolism enzymes in placenta

The placenta regulates the passage of both endogenous and xenobiotic compounds to the fetus during pregnancy. A small number of xenobiotic metabolizing CYP enzymes are constantly expressed in placenta, those include CYP19A1 which mainly converts androgens to estrogens, and CYP1A1 whose substrates include steroid hormones and xenobiotics. We performed an analysis of co-expression of xenobiotic metabolizing enzyme-coding genes and transcription factors in placenta and in 86 other tissues to discern the unique characteristics of the placental transcriptome. Transcription factors (TFs) driving the expression of proteins involved in phase I and II xenobiotic metabolism in the liver were not expressed in the placenta, nor were the bulk of xenobiotic metabolizing hepatic CYP enzymes. In contrast, TFs whose co-expression correlated with CYP1A1, i.e. AHR, PPARG, and CEBPB were highly expressed in placenta. The placenta is a hormonal tissue, and one needs to maintain the tissue-specific focus by removing the hepatic spectacles.