How tumor suppressor p53 selectively responds to specific signals, especially in normal cells, is poorly understood. We performed genome-wide profiling of p53 chromatin interactions and target gene expression in human embryonic stem cells (hESCs) in response to early differentiation, induced by retinoic acid, versus DNA damage, caused by adriamycin. Most p53-binding sites are unique to each state and define stimulus-specific p53 responses in hESCs. Differentiation-activated p53 targets include many developmental transcription factors and, in pluripotent hESCs, are bound by OCT4 and NANOG at chromatin enriched in both H3K27me3 and H3K4me3. Activation of these genes occurs with recruitment of p53 and H3K27me3-specific demethylases, UTX and JMJD3, to chromatin. In contrast, genes associated with cell migration and motility are bound by p53 specifically after DNA damage. Surveillance functions of p53 in cell death and cell cycle regulation are conserved in both DNA damage and differentiation. Comparative genomic analysis of p53-targets in mouse and human ESCs supports an inter-species divergence in p53 regulatory functions during evolution. Our findings expand the registry of p53-regulated genes to define p53-regulated opposition to pluripotency during early differentiation, a process highly distinct from stress-induced p53 response in hESCs.

Recent evidence suggests that lncRNAs play an integral regulatory role in numerous functions, including determination of cellular identity. We determined global expression (RNA-seq) and genome-wide profiles (ChIP-seq) of histone post-translational modifications and p53 binding in human embryonic stem cells (hESCs) undergoing differentiation to define a high-confidence set of 40 lncRNAs, which are p53 transcriptional targets. We focused on lncRNAs highly expressed in pluripotent hESCs and repressed by p53 during differentiation to identify lncPRESS1 as a p53-regulated transcript that maintains hESC pluripotency in concert with core pluripotency factors. RNA-seq of hESCs depleted of lncPRESS1 revealed that lncPRESS1 controls a gene network that promotes pluripotency. Further, we found that lncPRESS1 physically interacts with SIRT6 and prevents SIRT6 chromatin localization, which maintains high levels of histone H3K56 and H3K9 acetylation at promoters of pluripotency genes. In summary, we describe a p53-regulated, pluripotency-specific lncRNA that safeguards the hESC state by disrupting SIRT6 activity.

Recognition of modified histone species by distinct structural domains within 'reader' proteins plays a critical role in the regulation of gene expression. Readers that simultaneously recognize histones with multiple marks allow transduction of complex chromatin modification patterns into specific biological outcomes. Here we report that chromatin regulator tripartite motif-containing 24 (TRIM24) functions in humans as a reader of dual histone marks by means of tandem plant homeodomain (PHD) and bromodomain (Bromo) regions. The three-dimensional structure of the PHD-Bromo region of TRIM24 revealed a single functional unit for combinatorial recognition of unmodified H3K4 (that is, histone H3 unmodified at lysine 4, H3K4me0) and acetylated H3K23 (histone H3 acetylated at lysine 23, H3K23ac) within the same histone tail. TRIM24 binds chromatin and oestrogen receptor to activate oestrogen-dependent genes associated with cellular proliferation and tumour development. Aberrant expression of TRIM24 negatively correlates with survival of breast cancer patients. The PHD-Bromo of TRIM24 provides a structural rationale for chromatin activation through a non-canonical histone signature, establishing a new route by which chromatin readers may influence cancer pathogenesis.

Despite its economical, cultural, and biological importance, there has not been a large scale sequencing project to date for Camelus dromedarius. With the goal of sequencing complete DNA of the organism, we first established and sequenced camel EST libraries, generating 70,272 reads. Following trimming, chimera check, repeat masking, cluster and assembly, we obtained 23,602 putative gene sequences, out of which over 4,500 potentially novel or fast evolving gene sequences do not carry any homology to other available genomes. Functional annotation of sequences with similarities in nucleotide and protein databases has been obtained using Gene Ontology classification. Comparison to available full length cDNA sequences and Open Reading Frame (ORF) analysis of camel sequences that exhibit homology to known genes show more than 80% of the contigs with an ORF>300 bp and approximately 40% hits extending to the start codons of full length cDNAs suggesting successful characterization of camel genes. Similarity analyses are done separately for different organisms including human, mouse, bovine, and rat. Accompanying web portal, CAGBASE (http://camel.kacst.edu.sa/), hosts a relational database containing annotated EST sequences and analysis tools with possibility to add sequences from public domain. We anticipate our results to provide a home base for genomic studies of camel and other comparative studies enabling a starting point for whole genome sequencing of the organism.

Activation-induced cytidine deaminase, AICDA or AID, is a driver of somatic hypermutation and class-switch recombination in immunoglobulins. In addition, this deaminase belonging to the APOBEC family may have off-target effects genome-wide, but its effects at pan-cancer level are not well elucidated. Here, we used different pan-cancer datasets, totaling more than 50,000 samples analyzed by whole-genome, whole-exome, or targeted sequencing. AID mutations are present at pan-cancer level with higher frequency in hematological cancers and higher presence at transcriptionally active TAD domains. AID synergizes initial hotspot mutations by a second composite mutation. AID mutational load was found to be independently associated with a favorable outcome in immune-checkpoint inhibitors (ICI) treated patients across cancers after analyzing 2000 samples. Finally, we found that AID-related neoepitopes, resulting from mutations at more frequent hotspots if compared to other mutational signatures, enhance CXCL13/CCR5 expression, immunogenicity, and T-cell exhaustion, which may increase ICI sensitivity.

DNA cytosine methylation in mammals modulates gene expression and chromatin accessibility. It also impacts mutation rates, via spontaneous oxidative deamination of 5-methylcytosine (5mC) to thymine. In most cases the resulting T:G mismatches are repaired, following T excision by one of the thymine DNA glycosylases, TDG or MBD4. We found that C-to-T mutations are enriched in the binding sites of CCAAT/enhancer binding proteins (CEBP). Within a CEBP site, the presence of a T:G mismatch increased CEBPβ binding affinity by a factor of >60 relative to the normal C:G base pair. This enhanced binding to a mismatch inhibits its repair by both TDG and MBD4 in vitro. Furthermore, repair of the deamination product of unmethylated cytosine, which yields a U:G DNA mismatch that is normally repaired via uracil DNA glycosylase, is also inhibited by CEBPβ binding. Passage of a replication fork over either a T:G or U:G mismatch, before repair can occur, results in a C-to-T mutation in one of the daughter duplexes. Our study thus provides a plausible mechanism for accumulation of C-to-T human somatic mutations.

Chondrosarcomas are the most common primary bone tumor in adults. Isocitrate dehydrogenase 1 (IDH1) and IDH2 mutations are prevalent. We aimed to assess the clinico-genomic properties of IDH mutant versus IDH wild-type (WT) chondrosarcomas as well as alterations in other genes.

Chromatin is folded into successive layers to organize linear DNA. Genes within the same topologically associating domains (TADs) demonstrate similar expression and histone-modification profiles, and boundaries separating different domains have important roles in reinforcing the stability of these features. Indeed, domain disruptions in human cancers can lead to misregulation of gene expression. However, the frequency of domain disruptions in human cancers remains unclear. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), which aggregated whole-genome sequencing data from 2,658 cancers across 38 tumor types, we analyzed 288,457 somatic structural variations (SVs) to understand the distributions and effects of SVs across TADs. Notably, SVs can lead to the fusion of discrete TADs, and complex rearrangements markedly change chromatin folding maps in the cancer genomes. Notably, only 14% of the boundary deletions resulted in a change in expression in nearby genes of more than twofold.

Higher-order chromatin structure is important for the regulation of genes by distal regulatory sequences1,2. Structural variants (SVs) that alter three-dimensional (3D) genome organization can lead to enhancer-promoter rewiring and human disease, particularly in the context of cancer3. However, only a small minority of SVs are associated with altered gene expression4,5, and it remains unclear why certain SVs lead to changes in distal gene expression and others do not. To address these questions, we used a combination of genomic profiling and genome engineering to identify sites of recurrent changes in 3D genome structure in cancer and determine the effects of specific rearrangements on oncogene activation. By analysing Hi-C data from 92 cancer cell lines and patient samples, we identified loci affected by recurrent alterations to 3D genome structure, including oncogenes such as MYC, TERT and CCND1. By using CRISPR-Cas9 genome engineering to generate de novo SVs, we show that oncogene activity can be predicted by using 'activity-by-contact' models that consider partner region chromatin contacts and enhancer activity. However, activity-by-contact models are only predictive of specific subsets of genes in the genome, suggesting that different classes of genes engage in distinct modes of regulation by distal regulatory elements. These results indicate that SVs that alter 3D genome organization are widespread in cancer genomes and begin to illustrate predictive rules for the consequences of SVs on oncogene activation.

Somatic mutations in driver genes may ultimately lead to the development of cancer. Understanding how somatic mutations accumulate in cancer genomes and the underlying factors that generate somatic mutations is therefore crucial for developing novel therapeutic strategies. To understand the interplay between spatial genome organization and specific mutational processes, we studied 3,000 tumor-normal-pair whole-genome datasets from 42 different human cancer types. Our analyses reveal that the change in somatic mutational load in cancer genomes is co-localized with topologically-associating-domain boundaries. Domain boundaries constitute a better proxy to track mutational load change than replication timing measurements. We show that different mutational processes lead to distinct somatic mutation distributions where certain processes generate mutations in active domains, and others generate mutations in inactive domains. Overall, the interplay between three-dimensional genome organization and active mutational processes has a substantial influence on the large-scale mutation-rate variations observed in human cancers.