Mariann Micsinai

Genome-wide remodeling of the epigenetic landscape during myogenic differentiation

We have examined changes in the chromatin landscape during muscle differentiation by mapping the genome-wide location of ten key histone marks and transcription factors in mouse myoblasts and terminally differentiated myotubes, providing an exceptionally rich dataset that has enabled discovery of key epigenetic changes underlying myogenesis. Using this compendium, we focused on a well-known repressive mark, histone H3 lysine 27 trimethylation, and identified novel regulatory elements flanking the myogenin gene that function as a key differentiation-dependent switch during myogenesis. Next, we examined the role of Polycomb-mediated H3K27 methylation in gene repression by systematically ablating components of both PRC1 and PRC2 complexes. Surprisingly, we found mechanistic differences between transient and permanent repression of muscle differentiation and lineage commitment genes and observed that the loss of PRC1 and PRC2 components produced opposing differentiation defects. These phenotypes illustrate striking differences as compared to embryonic stem cell differentiation and suggest that PRC1 and PRC2 do not operate sequentially in muscle cells. Our studies of PRC1 occupancy also suggested a “fail-safe” mechanism, whereby PRC1/Bmi1 concentrates at genes specifying nonmuscle lineages, helping to retain H3K27me3 in the face of declining Ezh2-mediated methyltransferase activity in differentiated cells.

TrAp: a tree approach for fingerprinting subclonal tumor composition

Revealing the clonal composition of a single tumor is essential for identifying cell subpopulations with metastatic potential in primary tumors or with resistance to therapies in metastatic tumors. Sequencing technologies provide only an overview of the aggregate of numerous cells. Computational approaches to de-mix a collective signal composed of the aberrations of a mixed cell population of a tumor sample into its individual components are not available. We propose an evolutionary framework for deconvolving data from a single genome-wide experiment to infer the composition, abundance and evolutionary paths of the underlying cell subpopulations of a tumor. We have developed an algorithm (TrAp) for solving this mixture problem. In silico analyses show that TrAp correctly deconvolves mixed subpopulations when the number of subpopulations and the measurement errors are moderate. We demonstrate the applicability of the method using tumor karyotypes and somatic hypermutation data sets. We applied TrAp to Exome-Seq experiment of a renal cell carcinoma tumor sample and compared the mutational profile of the inferred subpopulations to the mutational profiles of single cells of the same tumor. Finally, we deconvolve sequencing data from eight acute myeloid leukemia patients and three distinct metastases of one melanoma patient to exhibit the evolutionary relationships of their subpopulations.

Picking ChIP-seq peak detectors for analyzing chromatin modification experiments

Numerous algorithms have been developed to analyze ChIP-Seq data. However, the complexity of analyzing diverse patterns of ChIP-Seq signals, especially for epigenetic marks, still calls for the development of new algorithms and objective comparisons of existing methods. We developed Qeseq, an algorithm to detect regions of increased ChIP read density relative to background. Qeseq employs critical novel elements, such as iterative recalibration and neighbor joining of reads to identify enriched regions of any length. To objectively assess its performance relative to other 14 ChIP-Seq peak finders, we designed a novel protocol based on Validation Discriminant Analysis (VDA) to optimally select validation sites and generated two validation datasets, which are the most comprehensive to date for algorithmic benchmarking of key epigenetic marks. In addition, we systematically explored a total of 315 diverse parameter configurations from these algorithms and found that typically optimal parameters in one dataset do not generalize to other datasets. Nevertheless, default parameters show the most stable performance, suggesting that they should be used. This study also provides a reproducible and generalizable methodology for unbiased comparative analysis of high-throughput sequencing tools that can facilitate future algorithmic development.

The mammalian Sin3 proteins are required for muscle development and sarcomere specification.

ABSTRACT The highly related mammalian Sin3A and Sin3B proteins provide a versatile platform for chromatin-modifying activities. Sin3-containing complexes play a role in gene repression through deacetylation of nucleosomes. Here, we explore a role for Sin3 in myogenesis by examining the phenotypes resulting from acute somatic deletion of both isoforms in vivo and from primary myotubes in vitro. Myotubes ablated for Sin3A alone, but not Sin3B, displayed gross defects in sarcomere structure that were considerably enhanced upon simultaneous ablation of both isoforms. Massively parallel sequencing of Sin3A- and Sin3B-bound genomic loci revealed a subset of target genes directly involved in sarcomere function that are positively regulated by Sin3A and Sin3B proteins. Both proteins were coordinately recruited to a substantial number of genes. Interestingly, depletion of Sin3B led to compensatory increases in Sin3A recruitment at certain target loci, but Sin3B was never found to compensate for Sin3A loss. Thus, our analyses describe a novel transcriptional role for Sin3A and Sin3B proteins associated with maintenance of differentiated muscle cells.

RUNX Transcription Factor-Mediated Association of Cd4 and Cd8 Enables Coordinate Gene Regulation

Summary T cell fate is associated with mutually exclusive expression of CD4 or CD8 in helper and cytotoxic T cells, respectively. How expression of one locus is temporally coordinated with repression of the other has been a long-standing enigma, though we know RUNX transcription factors activate the Cd8 locus, silence the Cd4 locus, and repress the Zbtb7b locus (encoding the transcription factor ThPOK), which is required for CD4 expression. Here we found that nuclear organization was altered by interplay among members of this transcription factor circuitry: RUNX binding mediated association of Cd4 and Cd8 whereas ThPOK binding kept the loci apart. Moreover, targeted deletions within Cd4 modulated CD8 expression and pericentromeric repositioning of Cd8. Communication between Cd4 and Cd8 thus appears to enable long-range epigenetic regulation to ensure that expression of one excludes the other in mature CD4 or CD8 single-positive (SP) cells.

The RAG2 C-terminus and ATM protect genome integrity by controlling antigen receptor gene cleavage

Tight control of antigen-receptor gene rearrangement is required to preserve genome integrity and prevent the occurrence of leukemia and lymphoma. Nonetheless, mistakes can happen, leading to the generation of aberrant rearrangements, such as Tcra/d-Igh inter-locus translocations that are a hallmark of ATM deficiency. Current evidence indicates that these translocations arise from the persistence of unrepaired breaks converging at different stages of thymocyte differentiation. Here we show that a defect in feedback control of RAG2 activity gives rise to bi-locus breaks and damage on Tcra/d and Igh in the same T cell at the same developmental stage, which provides a direct mechanism for generating these inter-locus rearrangements. Both the RAG2 C-terminus and ATM prevent bi-locus RAG-mediated cleavage through modulation of 3D conformation (higher order loops) and nuclear organization of the two loci. This limits the number of potential substrates for translocation and provides an important mechanism for protecting genome stability.

IL-7 Functionally Segregates the Pro-B Cell Stage by Regulating Transcription of Recombination Mediators across Cell Cycle

Ag receptor diversity involves the introduction of DNA double-stranded breaks during lymphocyte development. To ensure fidelity, cleavage is confined to the G0-G1 phase of the cell cycle. One established mechanism of regulation is through periodic degradation of the RAG2 recombinase protein. However, there are additional levels of protection. In this paper, we show that cyclical changes in the IL-7R signaling pathway functionally segregate pro-B cells according to cell cycle status. In consequence, the level of a downstream effector of IL-7 signaling, phospho-STAT5, is inversely correlated with cell cycle expression of Rag, a key gene involved in recombination. Higher levels of phopho-STAT5 in S-G2 correlate with decreased Rag expression and Rag relocalization to pericentromeric heterochromatin. These cyclical changes in transcription and locus repositioning are ablated upon transformation with v-Abl, which renders STAT5 constitutively active across the cell cycle. We propose that this activity of the IL-7R/STAT5 pathway plays a critical protective role in development, complementing regulation of RAG2 at the protein level, to ensure that recombination does not occur during replication. Our data, suggesting that pro-B cells are not a single homogeneous population, explain inconsistencies in the role of IL-7 signaling in regulating Igh recombination.

Combined Immunofluorescence and DNA FISH on 3D-preserved Interphase Nuclei to Study Changes in 3D Nuclear Organization

Fluorescent in situ hybridization using DNA probes on 3-dimensionally preserved nuclei followed by 3D confocal microscopy (3D DNA FISH) represents the most direct way to visualize the location of gene loci, chromosomal sub-regions or entire territories in individual cells. This type of analysis provides insight into the global architecture of the nucleus as well as the behavior of specific genomic loci and regions within the nuclear space. Immunofluorescence, on the other hand, permits the detection of nuclear proteins (modified histones, histone variants and modifiers, transcription machinery and factors, nuclear sub-compartments, etc). The major challenge in combining immunofluorescence and 3D DNA FISH is, on the one hand to preserve the epitope detected by the antibody as well as the 3D architecture of the nucleus, and on the other hand, to allow the penetration of the DNA probe to detect gene loci or chromosome territories (1-5). Here we provide a protocol that combines visualization of chromatin modifications with genomic loci in 3D preserved nuclei.

CONDEX: Copy number detection in exome sequences

Exome sequencing by hybrid capture facilitates obtaining cost-effective, comprehensive data on coding sequence variation from short reads. Standard analysis tools focus on detecting and characterizing variants of single or a few nucleotides while copy number variants (CNVs) that span multiple regions of exon baits have not yet been considered. Here, we present a Hidden Markov Model based method to identify CNVs from exome sequencing data. Using depth coverage and the heterozygosity of SNPs, we call CNVs with per-exon training data from other samples. The method has >90% accuracy in identification of deletions and insertions in simulations. Availability - http://code.google.com/p/condr/ (Java)

Response to Casellas et al.

Casellas and colleagues challenge the conclusions from our recent paper and claim that AID-dependent translocations occur independent of Igh proximity. In contrast to our study (Rocha et al., 2012), they focus solely on hotspots, which they claim represent the only ‘true’ AID mediated translocations (Hakim et al., 2012). However, their contention that other sites with translocation capture (TC) reads reflect AID independent rearrangements is not supported by any genome wide subtraction analysis examining signal enrichment in the IghI-SceI AID sufficient versus deficient sample. Nonetheless, they state that this is precisely the analysis that is missing from our study and that without this we could not distinguish between AID dependent and AID independent events. By implication, as per their suggestion we should find considerable overlap and comparable signal strength in the two data sets, which would account for the 90% of sites that they claim are AID independent translocations. Furthermore, non-overlapping regions or sites with signal enrichment in the AID sufficient data set should be restricted to AID induced hotspots. This is not what we observed.