Hendrik G. Stunnenberg

Analysis of the Plasmodium falciparum proteome by high-accuracy mass spectrometry

The annotated genomes of organisms define a ‘blueprint’ of their possible gene products. Post-genome analyses attempt to confirm and modify the annotation and impose a sense of the spatial, temporal and developmental usage of genetic information by the organism. Here we describe a large-scale, high-accuracy (average deviation less than 0.02 Da at 1,000 Da) mass spectrometric proteome analysis of selected stages of the human malaria parasite Plasmodium falciparum. The analysis revealed 1,289 proteins of which 714 proteins were identified in asexual blood stages, 931 in gametocytes and 645 in gametes. The last two groups provide insights into the biology of the sexual stages of the parasite, and include conserved, stage-specific, secreted and membrane-associated proteins. A subset of these proteins contain domains that indicate a role in cell–cell interactions, and therefore can be evaluated as potential components of a malaria vaccine formulation. We also report a set of peptides with significant matches in the parasite genome but not in the protein set predicted by computational methods.

The Transcriptional and Epigenomic Foundations of Ground State Pluripotency

Summary Mouse embryonic stem (ES) cells grown in serum exhibit greater heterogeneity in morphology and expression of pluripotency factors than ES cells cultured in defined medium with inhibitors of two kinases (Mek and GSK3), a condition known as “2i” postulated to establish a naive ground state. We show that the transcriptome and epigenome profiles of serum- and 2i-grown ES cells are distinct. 2i-treated cells exhibit lower expression of lineage-affiliated genes, reduced prevalence at promoters of the repressive histone modification H3K27me3, and fewer bivalent domains, which are thought to mark genes poised for either up- or downregulation. Nonetheless, serum- and 2i-grown ES cells have similar differentiation potential. Precocious transcription of developmental genes in 2i is restrained by RNA polymerase II promoter-proximal pausing. These findings suggest that transcriptional potentiation and a permissive chromatin context characterize the ground state and that exit from it may not require a metastable intermediate or multilineage priming.

Genome-wide profiling of PPARγ:RXR and RNA polymerase II occupancy reveals temporal activation of distinct metabolic pathways and changes in RXR dimer composition during adipogenesis

The nuclear receptor peroxisome proliferator-activated receptor gamma (PPARgamma) is a key regulator of adipocyte differentiation in vivo and ex vivo and has been shown to control the expression of several adipocyte-specific genes. In this study, we used chromatin immunoprecipitation combined with deep sequencing to generate genome-wide maps of PPARgamma and retinoid X receptor (RXR)-binding sites, and RNA polymerase II (RNAPII) occupancy at very high resolution throughout adipocyte differentiation of 3T3-L1 cells. We identify >5000 high-confidence shared PPARgamma:RXR-binding sites in adipocytes and show that during early stages of differentiation, many of these are preoccupied by non-PPARgamma RXR-heterodimers. Different temporal and compositional patterns of occupancy are observed. In addition, we detect co-occupancy with members of the C/EBP family. Analysis of RNAPII occupancy uncovers distinct clusters of similarly regulated genes of different biological processes. PPARgamma:RXR binding is associated with the majority of induced genes, and sites are particularly abundant in the vicinity of genes involved in lipid and glucose metabolism. Our analyses represent the first genome-wide map of PPARgamma:RXR target sites and changes in RNAPII occupancy throughout adipocyte differentiation and indicate that a hitherto unrecognized high number of adipocyte genes of distinctly regulated pathways are directly activated by PPARgamma:RXR.

Quantitative comparison of genome-wide DNA methylation mapping technologies

DNA methylation plays a key role in regulating eukaryotic gene expression. Although mitotically heritable and stable over time, patterns of DNA methylation frequently change in response to cell differentiation, disease and environmental influences. Several methods have been developed to map DNA methylation on a genomic scale. Here, we benchmark four of these approaches by analyzing two human embryonic stem cell lines derived from genetically unrelated embryos and a matched pair of colon tumor and adjacent normal colon tissue obtained from the same donor. Our analysis reveals that methylated DNA immunoprecipitation sequencing (MeDIP-seq), methylated DNA capture by affinity purification (MethylCap-seq), reduced representation bisulfite sequencing (RRBS) and the Infinium HumanMethylation27 assay all produce accurate DNA methylation data. However, these methods differ in their ability to detect differentially methylated regions between pairs of samples. We highlight strengths and weaknesses of the four methods and give practical recommendations for the design of epigenomic case-control studies.DNA methylation plays a key role in regulating eukaryotic gene expression. Although mitotically heritable and stable over time, patterns of DNA methylation frequently change in response to cell differentiation, disease and environmental influences. Several methods have been developed to map DNA methylation on a genomic scale. Here, we benchmark four of these approaches by analyzing two human embryonic stem cell lines derived from genetically unrelated embryos and a matched pair of colon tumor and adjacent normal colon tissue obtained from the same donor. Our analysis reveals that methylated DNA immunoprecipitation sequencing (MeDIP-seq), methylated DNA capture by affinity purification (MethylCap-seq), reduced representation bisulfite sequencing (RRBS) and the Infinium HumanMethylation27 assay all produce accurate DNA methylation data. However, these methods differ in their ability to detect differentially methylated regions between pairs of samples. We highlight strengths and weaknesses of the four methods and give practical recommendations for the design of epigenomic case-control studies.

Epigenetic programming of monocyte-to-macrophage differentiation and trained innate immunity

Introduction Monocytes circulate in the bloodstream for up to 3 to 5 days. Concomitantly, immunological imprinting of either tolerance (immunosuppression) or trained immunity (innate immune memory) determines the functional fate of monocytes and monocyte-derived macrophages, as observed after infection or vaccination. The epigenome, DNase I accessibility, and transcriptome were characterized in purified human circulating monocytes, in vitro differentiated naïve, tolerized (immunosuppression), and trained macrophages (innate immune memory). This allowed the identification of pathways functionally implicated in innate immune memory. This epigenetic signature of human monocyte-to-macrophage differentiation and monocyte training generates hypotheses to understand and manipulate medically relevant immune conditions. Methods Purified circulating monocytes from healthy volunteers were differentiated under the homeostatic macrophage colony-stimulating factor concentrations present in human serum. During the first 24 hours, trained immunity was induced by β-glucan (BG) priming, and postsepsis immunoparalysis was mimicked by exposure to lipopolysaccharide (LPS), generating endotoxin-induced tolerance. Epigenomic profiling of the histone marks H3K4me1, H3K4me3, and H3K27ac, DNase I accessibility, and RNA sequencing were performed at both the start of the experiment (ex vivo monocytes) and at the end of the 6 days of in vitro culture (macrophages). Results Compared with monocytes (Mo), naïve macrophages (Mf ) display a remodeled metabolic enzyme repertoire and attenuated innate inflammatory pathways, most likely necessary to generate functional tissue macrophages. Epigenetic profiling uncovered about 8000 dynamic regions associated with about 11,000 DNase I hypersensitive sites. Changes in histone acetylation identified most dynamic events. Furthermore, these regions of differential histone marks displayed some degree of DNase I accessibility that was already present in monocytes. H3K4me1 mark increased in parallel with de novo H3K27ac deposition at distal regulatory regions; H3K4me1 mark remained even after the loss of H3K27ac, marking decommissioned regulatory elements. β-glucan priming specifically induced about 3000 distal regulatory elements, whereas LPS tolerization induced H3K27ac at about 500 distal regulatory regions. At the transcriptional level, we identified coregulated gene modules during monocyte-to-macrophage differentiation, as well as discordant modules between trained and tolerized cells. These indicate that training likely involves an increased expression of modules expressed in naïve macrophages, including genes that code for metabolic enzymes. On the other hand, endotoxin tolerance involves gene modules that are more active in monocytes than in naïve macrophages. About 12% of known human transcription factors display variation in expression during macrophage differentiation, training, and tolerance. We also observed transcription factor motifs in DNase I hypersensitive sites at condition-specific dynamic epigenomic regions, implying that specific transcription factors are required for trained and tolerized macrophage epigenetic and transcriptional programs. Finally, our analyses and functional validation indicate that the inhibition of cyclic adenosine monophosphate generation blocked trained immunity in vitro and during an in vivo model of lethal Candida albicans infection, abolishing the protective effects of trained immunity. Discussion We documented the importance of epigenetic regulation of the immunological pathways underlying monocyte-to-macrophage differentiation and trained immunity. These dynamic epigenetic elements may inform on potential pharmacological targets that modulate innate immunity. Altogether, we uncovered the epigenetic and transcriptional programs of monocyte differentiation to macrophages that distinguish tolerant and trained macrophage phenotypes, providing a resource to further understand and manipulate immune-mediated responses. A BLUEPRINT of immune cell development To determine the epigenetic mechanisms that direct blood cells to develop into the many components of our immune system, the BLUEPRINT consortium examined the regulation of DNA and RNA transcription to dissect the molecular traits that govern blood cell differentiation. By inducing immune responses, Saeed et al. document the epigenetic changes in the genome that underlie immune cell differentiation. Cheng et al. demonstrate that trained monocytes are highly dependent on the breakdown of sugars in the presence of oxygen, which allows cells to produce the energy needed to mount an immune response. Chen et al. examine RNA transcripts and find that specific cell lineages use RNA transcripts of different length and composition (isoforms) to form proteins. Together, the studies reveal how epigenetic effects can drive the development of blood cells involved in the immune system. Science, this issue 10.1126/science.1251086, 10.1126/science.1250684, 10.1126/science.1251033 Genome-wide approaches analyze human monocyte differentiation in vitro into functional macrophages. Monocyte differentiation into macrophages represents a cornerstone process for host defense. Concomitantly, immunological imprinting of either tolerance or trained immunity determines the functional fate of macrophages and susceptibility to secondary infections. We characterized the transcriptomes and epigenomes in four primary cell types: monocytes and in vitro–differentiated naïve, tolerized, and trained macrophages. Inflammatory and metabolic pathways were modulated in macrophages, including decreased inflammasome activation, and we identified pathways functionally implicated in trained immunity. β-glucan training elicits an exclusive epigenetic signature, revealing a complex network of enhancers and promoters. Analysis of transcription factor motifs in deoxyribonuclease I hypersensitive sites at cell-type–specific epigenetic loci unveiled differentiation and treatment-specific repertoires. Altogether, we provide a resource to understand the epigenetic changes that underlie innate immunity in humans.

Bacille Calmette-Guérin induces NOD2-dependent nonspecific protection from reinfection via epigenetic reprogramming of monocytes

Adaptive features of innate immunity, recently described as “trained immunity,” have been documented in plants, invertebrate animals, and mice, but not yet in humans. Here we show that bacille Calmette-Guérin (BCG) vaccination in healthy volunteers led not only to a four- to sevenfold increase in the production of IFN-γ, but also to a twofold enhanced release of monocyte-derived cytokines, such as TNF and IL-1β, in response to unrelated bacterial and fungal pathogens. The enhanced function of circulating monocytes persisted for at least 3 mo after vaccination and was accompanied by increased expression of activation markers such as CD11b and Toll-like receptor 4. These training effects were induced through the NOD2 receptor and mediated by increased histone 3 lysine 4 trimethylation. In experimental studies, BCG vaccination induced T- and B-lymphocyte–independent protection of severe combined immunodeficiency SCID mice from disseminated candidiasis (100% survival in BCG-vaccinated mice vs. 30% in control mice). In conclusion, BCG induces trained immunity and nonspecific protection from infections through epigenetic reprogramming of innate immune cells.

Genome-wide analysis of ETS-family DNA-binding in vitro and in vivo

Members of the large ETS family of transcription factors (TFs) have highly similar DNA‐binding domains (DBDs)—yet they have diverse functions and activities in physiology and oncogenesis. Some differences in DNA‐binding preferences within this family have been described, but they have not been analysed systematically, and their contributions to targeting remain largely uncharacterized. We report here the DNA‐binding profiles for all human and mouse ETS factors, which we generated using two different methods: a high‐throughput microwell‐based TF DNA‐binding specificity assay, and protein‐binding microarrays (PBMs). Both approaches reveal that the ETS‐binding profiles cluster into four distinct classes, and that all ETS factors linked to cancer, ERG, ETV1, ETV4 and FLI1, fall into just one of these classes. We identify amino‐acid residues that are critical for the differences in specificity between all the classes, and confirm the specificities in vivo using chromatin immunoprecipitation followed by sequencing (ChIP‐seq) for a member of each class. The results indicate that even relatively small differences in in vitro binding specificity of a TF contribute to site selectivity in vivo.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

We analysed whole genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. 93 protein-coding cancer genes carried likely driver mutations. Some non-coding regions exhibited high mutation frequencies but most have distinctive structural features probably causing elevated mutation rates and do not harbour driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed 12 base substitution and six rearrangement signatures. Three rearrangement signatures, characterised by tandem duplications or deletions, appear associated with defective homologous recombination based DNA repair: one with deficient BRCA1 function; another with deficient BRCA1 or BRCA2 function; the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operative, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

ChIP-Seq of ERalpha and RNA polymerase II defines genes differentially responding to ligands.

We used ChIP‐Seq to map ERα‐binding sites and to profile changes in RNA polymerase II (RNAPII) occupancy in MCF‐7 cells in response to estradiol (E2), tamoxifen or fulvestrant. We identify 10 205 high confidence ERα‐binding sites in response to E2 of which 68% contain an estrogen response element (ERE) and only 7% contain a FOXA1 motif. Remarkably, 596 genes change significantly in RNAPII occupancy (59% up and 41% down) already after 1 h of E2 exposure. Although promoter proximal enrichment of RNAPII (PPEP) occurs frequently in MCF‐7 cells (17%), it is only observed on a minority of E2‐regulated genes (4%). Tamoxifen and fulvestrant partially reduce ERα DNA binding and prevent RNAPII loading on the promoter and coding body on E2‐upregulated genes. Both ligands act differently on E2‐downregulated genes: tamoxifen acts as an agonist thus downregulating these genes, whereas fulvestrant antagonizes E2‐induced repression and often increases RNAPII occupancy. Furthermore, our data identify genes preferentially regulated by tamoxifen but not by E2 or fulvestrant. Thus (partial) antagonist loaded ERα acts mechanistically different on E2‐activated and E2‐repressed genes.

Affinity purification of histidine-tagged proteins.

Expression of recombinant proteins is a standard technique in molecular biology and a wide variety of prokaryotic as well as eukaryotic expression systems are currently in use. A limiting step is often the purification of the expressed recombinant protein, particularly if mammalian expression systems that yield low expression levels are employed. Here, we discuss the advantages and restrictions of tagging recombinant proteins with histidines and purifying them by Ni2+-NTA chromatography.