Esther Rheinbay

Genomewide Analysis of PRC1 and PRC2 Occupancy Identifies Two Classes of Bivalent Domains

In embryonic stem (ES) cells, bivalent chromatin domains with overlapping repressive (H3 lysine 27 tri-methylation) and activating (H3 lysine 4 tri-methylation) histone modifications mark the promoters of more than 2,000 genes. To gain insight into the structure and function of bivalent domains, we mapped key histone modifications and subunits of Polycomb-repressive complexes 1 and 2 (PRC1 and PRC2) genomewide in human and mouse ES cells by chromatin immunoprecipitation, followed by ultra high-throughput sequencing. We find that bivalent domains can be segregated into two classes—the first occupied by both PRC2 and PRC1 (PRC1-positive) and the second specifically bound by PRC2 (PRC2-only). PRC1-positive bivalent domains appear functionally distinct as they more efficiently retain lysine 27 tri-methylation upon differentiation, show stringent conservation of chromatin state, and associate with an overwhelming number of developmental regulator gene promoters. We also used computational genomics to search for sequence determinants of Polycomb binding. This analysis revealed that the genomewide locations of PRC2 and PRC1 can be largely predicted from the locations, sizes, and underlying motif contents of CpG islands. We propose that large CpG islands depleted of activating motifs confer epigenetic memory by recruiting the full repertoire of Polycomb complexes in pluripotent cells.

Evolution of pathogenicity and sexual reproduction in eight Candida genomes.

Candida species are the most common cause of opportunistic fungal infection worldwide. Here we report the genome sequences of six Candida species and compare these and related pathogens and non-pathogens. There are significant expansions of cell wall, secreted and transporter gene families in pathogenic species, suggesting adaptations associated with virulence. Large genomic tracts are homozygous in three diploid species, possibly resulting from recent recombination events. Surprisingly, key components of the mating and meiosis pathways are missing from several species. These include major differences at the mating-type loci (MTL); Lodderomyces elongisporus lacks MTL, and components of the a1/α2 cell identity determinant were lost in other species, raising questions about how mating and cell types are controlled. Analysis of the CUG leucine-to-serine genetic-code change reveals that 99% of ancestral CUG codons were erased and new ones arose elsewhere. Lastly, we revise the Candida albicans gene catalogue, identifying many new genes.

Phylogenetic and structural analysis of centromeric DNA and kinetochore proteins

BackgroundKinetochores are large multi-protein structures that assemble on centromeric DNA (CEN DNA) and mediate the binding of chromosomes to microtubules. Comprising 125 base-pairs of CEN DNA and 70 or more protein components, Saccharomyces cerevisiae kinetochores are among the best understood. In contrast, most fungal, plant and animal cells assemble kinetochores on CENs that are longer and more complex, raising the question of whether kinetochore architecture has been conserved through evolution, despite considerable divergence in CEN sequence.ResultsUsing computational approaches, ranging from sequence similarity searches to hidden Markov model-based modeling, we show that organisms with CENs resembling those in S. cerevisiae (point CENs) are very closely related and that all contain a set of 11 kinetochore proteins not found in organisms with complex CENs. Conversely, organisms with complex CENs (regional CENs) contain proteins seemingly absent from point-CEN organisms. However, at least three quarters of known kinetochore proteins are present in all fungi regardless of CEN organization. At least six of these proteins have previously unidentified human orthologs. When fungi and metazoa are compared, almost all have kinetochores constructed around Spc105 and three conserved multi-protein linker complexes (MIND, COMA, and the NDC80 complex).ConclusionOur data suggest that critical structural features of kinetochores have been well conserved from yeast to man. Surprisingly, phylogenetic analysis reveals that human kinetochore proteins are as similar in sequence to their yeast counterparts as to presumptive Drosophila melanogaster or Caenorhabditis elegans orthologs. This finding is consistent with evidence that kinetochore proteins have evolved very rapidly relative to components of other complex cellular structures.

Imp2 controls oxidative phosphorylation and is crucial for preserving glioblastoma cancer stem cells

Growth of numerous cancer types is believed to be driven by a subpopulation of poorly differentiated cells, often referred to as cancer stem cells (CSCs), that have the capacity for self-renewal, tumor initiation, and generation of nontumorigenic progeny. Despite their potentially key role in tumor establishment and maintenance, the energy requirements of these cells and the mechanisms that regulate their energy production are unknown. Here, we show that the oncofetal insulin-like growth factor 2 mRNA-binding protein 2 (IMP2, IGF2BP2) regulates oxidative phosphorylation (OXPHOS) in primary glioblastoma (GBM) sphere cultures (gliomaspheres), an established in vitro model for CSC expansion. We demonstrate that IMP2 binds several mRNAs that encode mitochondrial respiratory chain complex subunits and that it interacts with complex I (NADH:ubiquinone oxidoreductase) proteins. Depletion of IMP2 in gliomaspheres decreases their oxygen consumption rate and both complex I and complex IV activity that results in impaired clonogenicity in vitro and tumorigenicity in vivo. Importantly, inhibition of OXPHOS but not of glycolysis abolishes GBM cell clonogenicity. Our observations suggest that gliomaspheres depend on OXPHOS for their energy production and survival and that IMP2 expression provides a key mechanism to ensure OXPHOS maintenance by delivering respiratory chain subunit-encoding mRNAs to mitochondria and contributing to complex I and complex IV assembly.

Complementary genomic approaches highlight the PI3K/mTOR pathway as a common vulnerability in osteosarcoma

Significance We present, to our knowledge, the first comprehensive next-generation sequencing of osteosarcoma in combination with a functional genomic screen in a genetically defined mouse model of osteosarcoma. Our data provide a strong rationale for targeting the phosphatidylinositol 3-kinase/mammalian target of rapamycin pathway in osteosarcoma and a foundation for rational clinical trial design. These findings present an immediate clinical opportunity because multiple inhibitors of this pathway are currently in clinical trials. Osteosarcoma is the most common primary bone tumor, yet there have been no substantial advances in treatment or survival in three decades. We examined 59 tumor/normal pairs by whole-exome, whole-genome, and RNA-sequencing. Only the TP53 gene was mutated at significant frequency across all samples. The mean nonsilent somatic mutation rate was 1.2 mutations per megabase, and there was a median of 230 somatic rearrangements per tumor. Complex chains of rearrangements and localized hypermutation were detected in almost all cases. Given the intertumor heterogeneity, the extent of genomic instability, and the difficulty in acquiring a large sample size in a rare tumor, we used several methods to identify genomic events contributing to osteosarcoma survival. Pathway analysis, a heuristic analytic algorithm, a comparative oncology approach, and an shRNA screen converged on the phosphatidylinositol 3-kinase/mammalian target of rapamycin (PI3K/mTOR) pathway as a central vulnerability for therapeutic exploitation in osteosarcoma. Osteosarcoma cell lines are responsive to pharmacologic and genetic inhibition of the PI3K/mTOR pathway both in vitro and in vivo.

H2A.Z landscapes and dual modifications in pluripotent and multipotent stem cells underlie complex genome regulatory functions

BackgroundThe histone variant H2A.Z has been implicated in nucleosome exchange, transcriptional activation and Polycomb repression. However, the relationships among these seemingly disparate functions remain obscure.ResultsWe mapped H2A.Z genome-wide in mammalian ES cells and neural progenitors. H2A.Z is deposited promiscuously at promoters and enhancers, and correlates strongly with H3K4 methylation. Accordingly, H2A.Z is present at poised promoters with bivalent chromatin and at active promoters with H3K4 methylation, but is absent from stably repressed promoters that are specifically enriched for H3K27 trimethylation. We also characterized post-translational modification states of H2A.Z, including a novel species dually-modified by ubiquitination and acetylation that is enriched at bivalent chromatin.ConclusionsOur findings associate H2A.Z with functionally distinct genomic elements, and suggest that post-translational modifications may reconcile its contrasting locations and roles.

EWS-FLI1 Utilizes Divergent Chromatin Remodeling Mechanisms to Directly Activate or Repress Enhancer Elements in Ewing Sarcoma

The aberrant transcription factor EWS-FLI1 drives Ewing sarcoma, but its molecular function is not completely understood. We find that EWS-FLI1 reprograms gene regulatory circuits in Ewing sarcoma by directly inducing or repressing enhancers. At GGAA repeat elements, which lack evolutionary conservation and regulatory potential in other cell types, EWS-FLI1 multimers induce chromatin opening and create de novo enhancers that physically interact with target promoters. Conversely, EWS-FLI1 inactivates conserved enhancers containing canonical ETS motifs by displacing wild-type ETS transcription factors. These divergent chromatin-remodeling patterns repress tumor suppressors and mesenchymal lineage regulators while activating oncogenes and potential therapeutic targets, such as the kinase VRK1. Our findings demonstrate how EWS-FLI1 establishes an oncogenic regulatory program governing both tumor survival and differentiation.

Recurrent and functional regulatory mutations in breast cancer

Genomic analysis of tumours has led to the identification of hundreds of cancer genes on the basis of the presence of mutations in protein-coding regions. By contrast, much less is known about cancer-causing mutations in non-coding regions. Here we perform deep sequencing in 360 primary breast cancers and develop computational methods to identify significantly mutated promoters. Clear signals are found in the promoters of three genes. FOXA1, a known driver of hormone-receptor positive breast cancer, harbours a mutational hotspot in its promoter leading to overexpression through increased E2F binding. RMRP and NEAT1, two non-coding RNA genes, carry mutations that affect protein binding to their promoters and alter expression levels. Our study shows that promoter regions harbour recurrent mutations in cancer with functional consequences and that the mutations occur at similar frequencies as in coding regions. Power analyses indicate that more such regions remain to be discovered through deep sequencing of adequately sized cohorts of patients.

A Tell-Tail Sign of Chromatin: Histone Mutations Drive Pediatric Glioblastoma

Recent genomic analyses of pediatric glioblastoma, a poorly understood tumor with dismal outcome, have identified mutations in histone H3 variants that affect critical amino acids in the tail. The findings extend discoveries of chromatin regulator inactivation and gain-of-function mutations by documenting alteration of a modifiable histone residue in human cancer.

Analysis of somatic microsatellite indels identifies driver events in human tumors

Microsatellites (MSs) are tracts of variable-length repeats of short DNA motifs that exhibit high rates of mutation in the form of insertions or deletions (indels) of the repeated motif. Despite their prevalence, the contribution of somatic MS indels to cancer has been largely unexplored, owing to difficulties in detecting them in short-read sequencing data. Here we present two tools: MSMuTect, for accurate detection of somatic MS indels, and MSMutSig, for identification of genes containing MS indels at a higher frequency than expected by chance. Applying MSMuTect to whole-exome data from 6,747 human tumors representing 20 tumor types, we identified >1,000 previously undescribed MS indels in cancer genes. Additionally, we demonstrate that the number and pattern of MS indels can accurately distinguish microsatellite-stable tumors from tumors with microsatellite instability, thus potentially improving classification of clinically relevant subgroups. Finally, we identified seven MS indel driver hotspots: four in known cancer genes (ACVR2A, RNF43, JAK1, and MSH3) and three in genes not previously implicated as cancer drivers (ESRP1, PRDM2, and DOCK3).