MicroRNAs play an important functional role in post-transcriptional gene regulation. One of the largest known microRNA clusters is located within the imprinted Dlk1/Gtl2 region on human chromosome 14 and mouse chromosome 12. This cluster contains more than 40 microRNA genes that are expressed only from the maternal chromosome in mouse.

A large subset of mammalian imprinted genes show extra-embryonic lineage (EXEL) specific imprinted expression that is restricted to placental trophectoderm lineages and to visceral yolk sac endoderm (ysE). Isolated ysE provides a homogenous in vivo model of a mid-gestation extra-embryonic tissue to examine the mechanism of EXEL-specific imprinted gene silencing, but an in vitro model of ysE to facilitate more rapid and cost-effective experiments is not available. Reports indicate that ES cells differentiated into cystic embryoid bodies (EBs) contain ysE, so here we investigate if cystic EBs model ysE imprinted expression. The imprinted expression pattern of cystic EBs is shown to resemble fetal liver and not ysE. To investigate the reason for this we characterized the methylome and transcriptome of cystic EBs in comparison to fetal liver and ysE, by whole genome bisulphite sequencing and RNA-seq. Cystic EBs show a fetal liver pattern of global hypermethylation and low expression of repeats, while ysE shows global hypomethylation and high expression of IAPEz retroviral repeats, as reported for placenta. Transcriptome analysis confirmed that cystic EBs are more similar to fetal liver than ysE and express markers of early embryonic endoderm. Genome-wide analysis shows that ysE shares epigenetic and repeat expression features with placenta. Contrary to previous reports, we show that cystic EBs do not contain ysE, but are more similar to the embryonic endoderm of fetal liver. This explains why cystic EBs reproduce the imprinted expression seen in the embryo but not that seen in the ysE.

CD8 coreceptor expression is dynamically regulated during thymocyte development and is tightly controlled by the activity of at least 5 different cis-regulatory elements. Despite the detailed characterization of the Cd8 loci, the regulation of the complex expression pattern of CD8 cannot be fully explained by the activity of the known Cd8 enhancers. In this study, we revisited the Cd8ab gene complex with bioinformatics and transgenic reporter gene expression approaches to search for additional Cd8 cis-regulatory elements. This led to the identification of an ECR (ECR-4), which in transgenic reporter gene expression assays, directed expression preferentially in CD44(hi)CD62L(+) CD8(+) T cells, including innate-like CD8(+) T cells. ECR-4, designated as Cd8 enhancer E8VI, was bound by Runx/CBFβ complexes and Bcl11b, indicating that E8VI is part of the cis-regulatory network that recruits transcription factors to the Cd8ab gene complex in CD8(+) T cells. Transgenic reporter expression was maintained in LCMV-specific CD8(+) T cells upon infection, although short-term, in vitro activation led to a down-regulation of E8VI activity. Finally, E8VI directed transgene expression also in CD8αα(+) DCs but not in CD8αα-expressing IELs. Taken together, we have identified a novel Cd8 enhancer that directs expression in CD44(hi)CD62L(+) CD8(+) T cells, including innate-like and antigen-specific effector/memory CD8(+) T cells and in CD8αα(+) DCs, and thus, our data provide further insight into the cis-regulatory networks that control CD8 expression.

Neurons are postmitotic and thus irreplaceable cells of the central nervous system (CNS). Accordingly, CNS inflammation with resulting neuronal damage can have devastating consequences. We investigated molecular mediators and structural consequences of CD8(+) T lymphocyte (CTL) attack on neurons in vivo. In a viral encephalitis model in mice, disease depended on CTL-derived interferon-γ (IFN-γ) and neuronal IFN-γ signaling. Downstream STAT1 phosphorylation and nuclear translocation in neurons were associated with dendrite and synapse loss (deafferentation). Analogous molecular and structural alterations were also found in human Rasmussen encephalitis, a CTL-mediated human autoimmune disorder of the CNS. Importantly, therapeutic intervention by IFN-γ blocking antibody prevented neuronal deafferentation and clinical disease without reducing CTL responses or CNS infiltration. These findings identify neuronal IFN-γ signaling as a novel target for neuroprotective interventions in CTL-mediated CNS disease.

RnBeads is a software tool for large-scale analysis and interpretation of DNA methylation data, providing a user-friendly analysis workflow that yields detailed hypertext reports (http://rnbeads.mpi-inf.mpg.de/). Supported assays include whole-genome bisulfite sequencing, reduced representation bisulfite sequencing, Infinium microarrays and any other protocol that produces high-resolution DNA methylation data. Notable applications of RnBeads include the analysis of epigenome-wide association studies and epigenetic biomarker discovery in cancer cohorts.

How functionally diverse populations of pathogen-specific killer T cells are generated during an immune response remains unclear. Here, we propose that fine-tuning of CD8αβ co-receptor levels via histone acetylation plays a role in lineage fate. We show that lysine acetyltransferase 6A (KAT6A) is responsible for maintaining permissive Cd8 gene transcription and enabling robust effector responses during infection. KAT6A-deficient CD8(+) T cells downregulated surface CD8 co-receptor expression during clonal expansion, a finding linked to reduced Cd8α transcripts and histone-H3 lysine 9 acetylation of the Cd8 locus. Loss of CD8 expression in KAT6A-deficient T cells correlated with reduced TCR signaling intensity and accelerated contraction of the effector-like memory compartment, whereas the long-lived memory compartment appeared unaffected, a result phenocopied by the removal of the Cd8 E8I enhancer element. These findings suggest a direct role of CD8αβ co-receptor expression and histone acetylation in shaping functional diversity within the cytotoxic T cell pool.

H1 linker histones facilitate higher-order chromatin folding and are essential for mammalian development. To achieve high-resolution mapping of H1 variants H1d and H1c in embryonic stem cells (ESCs), we have established a knock-in system and shown that the N-terminally tagged H1 proteins are functionally interchangeable to their endogenous counterparts in vivo. H1d and H1c are depleted from GC- and gene-rich regions and active promoters, inversely correlated with H3K4me3, but positively correlated with H3K9me3 and associated with characteristic sequence features. Surprisingly, both H1d and H1c are significantly enriched at major satellites, which display increased nucleosome spacing compared with bulk chromatin. While also depleted at active promoters and enriched at major satellites, overexpressed H1(0) displays differential binding patterns in specific repetitive sequences compared with H1d and H1c. Depletion of H1c, H1d, and H1e causes pericentric chromocenter clustering and de-repression of major satellites. These results integrate the localization of an understudied type of chromatin proteins, namely the H1 variants, into the epigenome map of mouse ESCs, and we identify significant changes at pericentric heterochromatin upon depletion of this epigenetic mark.

In contrast to its close homolog CDK4, the cell cycle kinase CDK6 is expressed at high levels in lymphoid malignancies. In a model for p185BCR-ABL+ B-acute lymphoid leukemia, we show that CDK6 is part of a transcription complex that induces the expression of the tumor suppressor p16INK4a and the pro-angiogenic factor VEGF-A. This function is independent of CDK6's kinase activity. High CDK6 expression thus suppresses proliferation by upregulating p16INK4a, providing an internal safeguard. However, in the absence of p16INK4a, CDK6 can exert its full tumor-promoting function by enhancing proliferation and stimulating angiogenesis. The finding that CDK6 connects cell-cycle progression to angiogenesis confirms CDK6's central role in hematopoietic malignancies and could underlie the selection pressure to upregulate CDK6 and silence p16INK4a.

Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for DNA methylation analysis in non-model organisms, we developed an integrated approach for studying DNA methylation differences independent of a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, carp, sea bass, and zebrafish). Bioinformatically, we developed the RefFreeDMA software to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions between samples or groups of individuals (http://RefFreeDMA.computational-epigenetics.org). The identified regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell-type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome.

Chronic lymphocytic leukaemia (CLL) is characterized by substantial clinical heterogeneity, despite relatively few genetic alterations. To provide a basis for studying epigenome deregulation in CLL, here we present genome-wide chromatin accessibility maps for 88 CLL samples from 55 patients measured by the ATAC-seq assay. We also performed ChIPmentation and RNA-seq profiling for ten representative samples. Based on the resulting data set, we devised and applied a bioinformatic method that links chromatin profiles to clinical annotations. Our analysis identified sample-specific variation on top of a shared core of CLL regulatory regions. IGHV mutation status-which distinguishes the two major subtypes of CLL-was accurately predicted by the chromatin profiles and gene regulatory networks inferred for IGHV-mutated versus IGHV-unmutated samples identified characteristic differences between these two disease subtypes. In summary, we discovered widespread heterogeneity in the chromatin landscape of CLL, established a community resource for studying epigenome deregulation in leukaemia and demonstrated the feasibility of large-scale chromatin accessibility mapping in cancer cohorts and clinical research.

Antiviral responses must be tightly regulated to defend rapidly against infection while minimizing inflammatory damage. Type 1 interferons (IFN-I) are crucial mediators of antiviral responses and their transcription is regulated by a variety of transcription factors; principal among these is the family of interferon regulatory factors (IRFs). The IRF gene regulatory networks are complex and contain multiple feedback loops. The tools of systems biology are well suited to elucidate the complex interactions that give rise to precise coordination of the interferon response. Here we have used an unbiased systems approach to predict that a member of the forkhead family of transcription factors, FOXO3, is a negative regulator of a subset of antiviral genes. This prediction was validated using macrophages isolated from Foxo3-null mice. Genome-wide location analysis combined with gene deletion studies identified the Irf7 gene as a critical target of FOXO3. FOXO3 was identified as a negative regulator of Irf7 transcription and we have further demonstrated that FOXO3, IRF7 and IFN-I form a coherent feed-forward regulatory circuit. Our data suggest that the FOXO3-IRF7 regulatory circuit represents a novel mechanism for establishing the requisite set points in the interferon pathway that balances the beneficial effects and deleterious sequelae of the antiviral response.

Arenaviruses such as Lassa fever virus cause significant mortality in endemic areas and represent potential bioterrorist weapons. The occurrence of arenaviral hemorrhagic fevers is largely confined to Third World countries with a limited medical infrastructure, and therefore live-attenuated vaccines have long been sought as a method of choice for prevention. Yet their rational design and engineering have been thwarted by technical limitations. In addition, viral genes had not been identified that are needed to cause disease but can be deleted or substituted to generate live-attenuated vaccine strains. Lymphocytic choriomeningitis virus, the prototype arenavirus, induces cell-mediated immunity against Lassa fever virus, but its safety for humans is unclear and untested. Using this virus model, we have developed the necessary methodology to efficiently modify arenavirus genomes and have exploited these techniques to identify an arenaviral Achilles' heel suitable for targeting in vaccine design. Reverse genetic exchange of the viral glycoprotein for foreign glycoproteins created attenuated vaccine strains that remained viable although unable to cause disease in infected mice. This phenotype remained stable even after extensive propagation in immunodeficient hosts. Nevertheless, the engineered viruses induced T cell-mediated immunity protecting against overwhelming systemic infection and severe liver disease upon wild-type virus challenge. Protection was established within 3 to 7 d after immunization and lasted for approximately 300 d. The identification of an arenaviral Achilles' heel demonstrates that the reverse genetic engineering of live-attenuated arenavirus vaccines is feasible. Moreover, our findings offer lymphocytic choriomeningitis virus or other arenaviruses expressing foreign glycoproteins as promising live-attenuated arenavirus vaccine candidates.

Enhancers are genetic elements that regulate spatiotemporal gene expression. Enhancer function requires transcription factor (TF) binding and correlates with histone modifications. However, the extent to which TF binding and histone modifications functionally define active enhancers remains unclear. Here, we combine chromatin immunoprecipitation with a massively parallel reporter assay (ChIP-STARR-seq) to identify functional enhancers in human embryonic stem cells (ESCs) genome-wide in a quantitative unbiased manner. Although active enhancers associate with TFs, only a minority of regions marked by NANOG, OCT4, H3K27ac, and H3K4me1 function as enhancers, with activity markedly changing under naive versus primed culture conditions. We identify an enhancer set associated with functions extending to non-ESC-specific processes. Moreover, although transposable elements associate with putative enhancers, only some exhibit activity. Similarly, within super-enhancers, large tracts are non-functional, with activity restricted to small sub-domains. This catalog of validated enhancers provides a valuable resource for further functional dissection of the regulatory genome.

The marsupial Tasmanian devil (Sarcophilus harrisii) faces extinction due to transmissible devil facial tumor disease (DFTD). To unveil the molecular underpinnings of this transmissible cancer, we combined pharmacological screens with an integrated systems-biology characterization. Sensitivity to inhibitors of ERBB tyrosine kinases correlated with their overexpression. Proteomic and DNA methylation analyses revealed tumor-specific signatures linked to the evolutionary conserved oncogenic STAT3. ERBB inhibition blocked phosphorylation of STAT3 and arrested cancer cells. Pharmacological blockade of ERBB or STAT3 prevented tumor growth in xenograft models and restored MHC class I expression. This link between the hyperactive ERBB-STAT3 axis and major histocompatibility complex class I-mediated tumor immunosurveillance provides mechanistic insights into horizontal transmissibility and puts forward a dual chemo-immunotherapeutic strategy to save Tasmanian devils from DFTD. VIDEO ABSTRACT.

Lipid metabolism is increasingly being appreciated to affect immunoregulation, inflammation and pathology. In this study we found that mice infected with lymphocytic choriomeningitis virus (LCMV) exhibit global perturbations of circulating serum lipids. Mice lacking the lipid-sensing surface receptor triggering receptor expressed on myeloid cells 2 (Trem2 -/-) were protected from LCMV-induced hepatitis and showed improved virus control despite comparable virus-specific T cell responses. Non-hematopoietic expression of TREM2 was found to be responsible for aggravated hepatitis, indicating a novel role for TREM2 in the non-myeloid compartment. These results suggest a link between virus-perturbed lipids and TREM2 that modulates liver pathogenesis upon viral infection. Targeted interventions of this immunoregulatory axis may ameliorate tissue pathology in hepatitis.

Tumor associated inflammation predicts response to immune checkpoint blockade in human melanoma. Current theories on regulation of inflammation center on anti-tumor T cell responses. Here we show that tumor associated B cells are vital to melanoma associated inflammation. Human B cells express pro- and anti-inflammatory factors and differentiate into plasmablast-like cells when exposed to autologous melanoma secretomes in vitro. This plasmablast-like phenotype can be reconciled in human melanomas where plasmablast-like cells also express T cell-recruiting chemokines CCL3, CCL4, CCL5. Depletion of B cells in melanoma patients by anti-CD20 immunotherapy decreases tumor associated inflammation and CD8+ T cell numbers. Plasmablast-like cells also increase PD-1+ T cell activation through anti-PD-1 blockade in vitro and their frequency in pretherapy melanomas predicts response and survival to immune checkpoint blockade. Tumor associated B cells therefore orchestrate and sustain melanoma inflammation and may represent a predictor for survival and response to immune checkpoint blockade therapy.

Interactions between commensal microbes and the immune system are tightly regulated and maintain intestinal homeostasis, but little is known about these interactions in humans. We investigated responses of human CD4+ T cells to the intestinal microbiota. We measured the abundance of T cells in circulation and intestinal tissues that respond to intestinal microbes and determined their clonal diversity. We also assessed their functional phenotypes and effects on intestinal resident cell populations, and studied alterations in microbe-reactive T cells in patients with chronic intestinal inflammation.

Maintenance of genome integrity via repair of DNA damage is a key biological process required to suppress diseases, including Fanconi anemia (FA). We generated loss-of-function human haploid cells for FA complementation group C (FANCC), a gene encoding a component of the FA core complex, and used genome-wide CRISPR libraries as well as insertional mutagenesis to identify synthetic viable (genetic suppressor) interactions for FA. Here we show that loss of the BLM helicase complex suppresses FANCC phenotypes and we confirm this interaction in cells deficient for FA complementation group I and D2 (FANCI and FANCD2) that function as part of the FA I-D2 complex, indicating that this interaction is not limited to the FA core complex, hence demonstrating that systematic genome-wide screening approaches can be used to reveal genetic viable interactions for DNA repair defects.

Glioblastoma might have widespread effects on the neural organization and cognitive function, and even focal lesions may be associated with distributed functional alterations. However, functional changes do not necessarily follow obvious anatomical patterns and the current understanding of this interrelation is limited. In this study, we used resting-state functional magnetic resonance imaging to evaluate changes in global functional connectivity patterns in 15 patients with glioblastoma. For six patients we followed longitudinal trajectories of their functional connectome and structural tumour evolution using bi-monthly follow-up scans throughout treatment and disease progression. In all patients, unilateral tumour lesions were associated with inter-hemispherically symmetric network alterations, and functional proximity of tumour location was stronger linked to distributed network deterioration than anatomical distance. In the longitudinal subcohort of six patients, we observed patterns of network alterations with initial transient deterioration followed by recovery at first follow-up, and local network deterioration to precede structural tumour recurrence by two months. In summary, the impact of focal glioblastoma lesions on the functional connectome is global and linked to functional proximity rather than anatomical distance to tumour regions. Our findings further suggest a relevance for functional network trajectories as a possible means supporting early detection of tumour recurrence.

B cell and plasma cell fates are controlled by different transcriptional networks, as exemplified by the mutually exclusive expression and cross-antagonism of the B cell identity factor Pax5 and the plasma cell regulator Blimp1. It has been postulated that repression of Pax5 by Blimp1 is essential for plasma cell development. Here, we challenged this hypothesis by analyzing the IghPax5/+ mouse, which expressed a Pax5 minigene from the immunoglobulin heavy-chain locus. Despite high Pax5 expression, plasma cells efficiently developed in young IghPax5/+ mice at steady state and upon immunization, while their number moderately declined in older mice. Although Pax5 significantly deregulated the plasma cell expression program, key plasma cell regulators were normally expressed in IghPax5/+ plasma cells. While IgM and IgA secretion by IghPax5/+ plasma cells was normal, IgG secretion was modestly decreased. Hence, Pax5 repression is not essential for robust plasma cell development and antibody secretion, although it is required for optimal IgG production and accumulation of long-lived plasma cells.