Anke Wienecke-Baldacchino

Interferon-γ-induced activation of Signal Transducer and Activator of Transcription 1 (STAT1) up-regulates the tumor suppressing microRNA-29 family in melanoma cells

BackgroundThe type-II-cytokine IFN-γ is a pivotal player in innate immune responses but also assumes functions in controlling tumor cell growth by orchestrating cellular responses against neoplastic cells. The role of IFN-γ in melanoma is not fully understood: it is a well-known growth inhibitor of melanoma cells in vitro. On the other hand, IFN-γ may also facilitate melanoma progression. While interferon-regulated genes encoding proteins have been intensively studied since decades, the contribution of miRNAs to effects mediated by interferons is an emerging area of research.We recently described a distinct and dynamic regulation of a whole panel of microRNAs (miRNAs) after IFN-γ-stimulation. The aim of this study was to analyze the transcriptional regulation of miR-29 family members in detail, identify potential interesting target genes and thus further elucidate a potential signaling pathway IFN-γ → Jak→ P-STAT1 → miR-29 → miR-29 target genes and its implication for melanoma growth.ResultsHere we show that IFN-γ induces STAT1-dependently a profound up-regulation of the miR-29 primary cluster pri-29a~b-1 in melanoma cell lines. Furthermore, expression levels of pri-29a~b-1 and mature miR-29a and miR-29b were elevated while the pri-29b-2~c cluster was almost undetectable. We observed an inverse correlation between miR-29a/b expression and the proliferation rate of various melanoma cell lines. This finding could be corroborated in cells transfected with either miR-29 mimics or inhibitors. The IFN-γ-induced G1-arrest of melanoma cells involves down-regulation of CDK6, which we proved to be a direct target of miR-29 in these cells. Compared to nevi and normal skin, and metastatic melanoma samples, miR-29a and miR-29b levels were found strikingly elevated in certain patient samples derived from primary melanoma.ConclusionsOur findings reveal that the miR-29a/b1 cluster is to be included in the group of IFN- and STAT-regulated genes. The up-regulated miR-29 family members may act as effectors of cytokine signalling in melanoma and other cancer cells as well as in the immune system.

Gene-pair expression signatures reveal lineage control

The distinct cell types of multicellular organisms arise due to constraints imposed by gene regulatory networks on the collective change of gene expression across the genome, creating self-stabilizing expression states, or attractors. We compiled a resource of curated human expression data comprising 166 cell types and 2,602 transcription regulating genes and developed a data driven method built around the concept of expression reversal defined at the level of gene pairs, such as those participating in toggle switch circuits. This approach allows us to organize the cell types into their ontogenetic lineage-relationships and to reflect regulatory relationships among genes that explain their ability to function as determinants of cell fate. We show that this method identifies genes belonging to regulatory circuits that control neuronal fate, pluripotency and blood cell differentiation, thus offering a novel large-scale perspective on lineage specification.The distinct cell types of multicellular organisms arise owing to constraints imposed by gene regulatory networks on the collective change of gene expression across the genome, creating self-stabilizing expression states, or attractors. We curated human expression data comprising 166 cell types and 2,602 transcription-regulating genes and developed a data-driven method for identifying putative determinants of cell fate built around the concept of expression reversal of gene pairs, such as those participating in toggle-switch circuits. This approach allows us to organize the cell types into their ontogenic lineage relationships. Our method identifies genes in regulatory circuits that control neuronal fate, pluripotency and blood cell differentiation, and it may be useful for prioritizing candidate factors for direct conversion of cell fate.

Dataset integration identifies transcriptional regulation of microRNA genes by PPARγ in differentiating mouse 3T3-L1 adipocytes

Peroxisome proliferator-activated receptor γ (PPARγ) is a key transcription factor in mammalian adipogenesis. Genome-wide approaches have identified thousands of PPARγ binding sites in mouse adipocytes and PPARγ upregulates hundreds of protein-coding genes during adipogenesis. However, no microRNA (miRNA) genes have been identified as primary PPARγ-targets. By integration of four separate datasets of genome-wide PPARγ binding sites in 3T3-L1 adipocytes we identified 98 miRNA clusters with PPARγ binding within 50 kb from miRNA transcription start sites. Nineteen mature miRNAs were upregulated ≥2-fold during adipogenesis and for six of these miRNA loci the PPARγ binding sites were confirmed by at least three datasets. The upregulation of five miRNA genes miR-103-1 (host gene Pank3), miR-148b (Copz1), miR-182/96/183, miR-205 and miR-378 (Ppargc1b) followed that of Pparg. The PPARγ-dependence of four of these miRNA loci was demonstrated by PPARγ knock-down and the loci of miR-103-1 (Pank3), miR-205 and miR-378 (Ppargc1b) were also responsive to the PPARγ ligand rosiglitazone. Finally, chromatin immunoprecipitation analysis validated in silico predicted PPARγ binding sites at all three loci and H3K27 acetylation was analyzed to confirm the activity of these enhancers. In conclusion, we identified 22 putative PPARγ target miRNA genes, showed the PPARγ dependence of four of these genes and demonstrated three as direct PPARγ target genes in mouse adipogenesis.

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.

Dynamic regulation of microRNA expression following Interferon-γ-induced gene transcription

MicroRNAs are major players in post-transcriptional gene regulation. Even small changes in miRNA levels may have profound consequences for the expression levels of target genes. Hence, miRNAs themselves need to be tightly, albeit dynamically, regulated. Here, we investigated the dynamic behavior of miRNAs over a wide time range following stimulation of melanoma cells with interferon-γ (IFN-γ), which activates the transcription factor STAT1. By applying several bioinformatic and statistical software tools for visualization and identification of differentially expressed miRNAs derived from time-series microarray experiments, 8.9% of 1105 miRNAs appeared to be directly or indirectly regulated by STAT1. Focusing on distinct dynamic expression patterns, we found that the majority of robust miRNA expression changes occurred in the intermediate time range (24–48 h). Three miRNAs (miR-27a, miR-30a, miR-34a) had a delayed regulation occurring at 72 h while none showed significant expression changes at early time points between 30 min and 6 h. Expression patterns of individual miRNAs were altered gradually over time or abruptly increased or decreased between two time points. Furthermore, we observed coordinated dynamic transcription of most miRNA clusters while few were found to be regulated independently of their genetic cluster. Most interestingly, several “star” or passenger strand sequences were specifically regulated over time while their “guide” strands were not.

The Effects of Sequence Variation on Genome-wide NRF2 Binding—New Target Genes and Regulatory SNPs

Transcription factor binding specificity is crucial for proper target gene regulation. Motif discovery algorithms identify the main features of the binding patterns, but the accuracy on the lower affinity sites is often poor. Nuclear factor E2-related factor 2 (NRF2) is a ubiquitous redox-activated transcription factor having a key protective role against endogenous and exogenous oxidant and electrophile stress. Herein, we decipher the effects of sequence variation on the DNA binding sequence of NRF2, in order to identify both genome-wide binding sites for NRF2 and disease-associated regulatory SNPs (rSNPs) with drastic effects on NRF2 binding. Interactions between NRF2 and DNA were studied using molecular modelling, and NRF2 chromatin immunoprecipitation-sequence datasets together with protein binding microarray measurements were utilized to study binding sequence variation in detail. The binding model thus generated was used to identify genome-wide binding sites for NRF2, and genomic binding sites with rSNPs that have strong effects on NRF2 binding and reside on active regulatory elements in human cells. As a proof of concept, miR-126–3p and -5p were identified as NRF2 target microRNAs, and a rSNP (rs113067944) residing on NRF2 target gene (Ferritin, light polypeptide, FTL) promoter was experimentally verified to decrease NRF2 binding and result in decreased transcriptional activity.

Bioinformatics as a driver, not a passenger, of translational biomedical research: Perspectives from the 6th Benelux bioinformatics conference

The 6th Benelux Bioinformatics Conference (BBC11) held in Luxembourg on 12 and 13 December 2011 attracted around 200 participants, including internationally-renowned guest speakers and more than 100 peer-reviewed submissions from 3 continents. Researchers from the public and private sectors convened at BBC11 to discuss advances and challenges in a wide spectrum of application areas. A key theme of the conference was the contribution of bioinformatics to enable and accelerate translational and clinical research. The BBC11 stressed the need for stronger collaborating efforts across disciplines and institutions. The demonstration of the clinical relevance of systems approaches and of next-generation sequencing-based measurement technologies are among the existing opportunities for increasing impact in translational research. Translational bioinformatics will benefit from research models that strike a balance between the importance of protecting intellectual property and the need to openly access scientific and technological advances. The full conference proceedings are freely available at http://www.bbc11.lu.

Functionally Convergent B Cell Receptor Sequences in Transgenic Rats Expressing a Human B Cell Repertoire in Response to Tetanus Toxoid and Measles Antigens

The identification and tracking of antigen-specific immunoglobulin (Ig) sequences within total Ig repertoires is central to high-throughput sequencing (HTS) studies of infections or vaccinations. In this context, public Ig sequences shared by different individuals exposed to the same antigen could be valuable markers for tracing back infections, measuring vaccine immunogenicity, and perhaps ultimately allow the reconstruction of the immunological history of an individual. Here, we immunized groups of transgenic rats expressing human Ig against tetanus toxoid (TT), Modified Vaccinia virus Ankara (MVA), measles virus hemagglutinin and fusion proteins expressed on MVA, and the environmental carcinogen benzo[a]pyrene, coupled to TT. We showed that these antigens impose a selective pressure causing the Ig heavy chain (IgH) repertoires of the rats to converge toward the expression of antibodies with highly similar IgH CDR3 amino acid sequences. We present a computational approach, similar to differential gene expression analysis, that selects for clusters of CDR3s with 80% similarity, significantly overrepresented within the different groups of immunized rats. These IgH clusters represent antigen-induced IgH signatures exhibiting stereotypic amino acid patterns including previously described TT- and measles-specific IgH sequences. Our data suggest that with the presented methodology, transgenic Ig rats can be utilized as a model to identify antigen-induced, human IgH signatures to a variety of different antigens.

Extraction and Analysis of RNA Isolated from Pure Bacteria-Derived Outer Membrane Vesicles

Outer membrane vesicles (OMVs) are released by commensal as well as pathogenic Gram-negative bacteria. These vesicles contain numerous bacterial components, such as proteins, peptidoglycans, lipopolysaccharides, DNA, and RNA. To examine if OMV-associated RNA molecules are bacterial degradation products and/or are functionally active, it is necessary to extract RNA from pure OMVs for subsequent analysis. Therefore, we describe here an isolation method of ultrapure OMVs and the subsequent extraction of RNA and basic steps of RNA-Seq analysis. Bacterial culture, extracellular supernatant concentration, OMV purification, and the subsequent RNA extraction out of OMVs are described. Specific pitfalls within the protocol and RNA contamination sources are highlighted.

High-throughput sequencing of murine immunoglobulin heavy chain repertoires using single side unique molecular identifiers on an Ion Torrent PGM

With the advent of high-throughput sequencing (HTS), profiling immunoglobulin (IG) repertoires has become an essential part of immunological research. Advances in sequencing technology enable the IonTorrent Personal Genome Machine (PGM) to cover the full-length of IG mRNA transcripts. Nucleotide insertions and deletions (indels) are the dominant errors of the PGM sequencing platform and can critically influence IG repertoire assessments. Here, we present a PGM-tailored IG repertoire sequencing approach combining error correction through unique molecular identifier (UID) barcoding and indel detection through ImMunoGeneTics (IMGT), the most commonly used sequence alignment database for IG sequences. Using artificially falsified sequences for benchmarking, we found that IMGTs underlying algorithms efficiently detect 98% of the introduced indels. Undetected indels are either located at the end of the sequences or produce masked frameshifts with an insertion and deletion in close proximity. The complementary determining regions 3 (CDR3s) are returned correct for up to 3 insertions or 3 deletions through conservative culling. We further show, that our PGM-tailored unique molecular identifiers result in highly accurate HTS data if combined with the presented processing strategy. In this regard, considering sequences with at least two copies from datasets with UID families of minimum 3 reads result in correct sequences with over 99% confidence. Finally, we show that the protocol can readily be used to generate homogenous datasets for bulk sequencing of murine bone marrow samples. Taken together, this approach will help to establish benchtop-scale sequencing of IG heavy chain transcripts in the field of IG repertoire research.