Genes encoding histone proteins are recurrently mutated in tumor samples, and these mutations may impact nucleosome stability, histone post-translational modification, or chromatin dynamics. The prevalence of histone mutations across diverse cancer types suggest that normal chromatin structure is a barrier to tumorigenesis. Oncohistone mutations disrupt chromatin structure and gene regulatory mechanisms, resulting in aberrant gene expression and the development of cancer phenotypes. Examples of oncohistones include the histone H3 K27M mutation found in pediatric brain cancers that blocks post-translational modification of the H3 N-terminal tail and the histone H2B E76K mutation found in some solid tumors that disrupts nucleosome stability. Oncohistones may comprise a limited fraction of the total histone pool yet cause global effects on chromatin structure and drive cancer phenotypes. Here, we survey histone mutations in cancer and review their function and role in tumorigenesis.
Introduction
Histone proteins play a fundamental role in protecting and regulating access to DNA during replication, repair, and transcription. The four core histones (H2A, H2B, H3, and H4) are well-conserved basic proteins ranging in size from 102 to 135 amino acids. H2A forms a dimer with H2B, while an H3 dimer interacts with H4 to form the H3/H4 tetramer [1]. Two H2A/H2B dimers bind with one H3/H4 tetramer to form the histone octamer that interacts with ∼147 bp of negatively charged DNA to form a nucleosome [1,2]. Core histone proteins have a positively charged N- or C-terminal tail that protrudes from the nucleosome and histone-fold domains that associate with 27–28 bp of DNA [2]. In addition, the linker histone H1 dynamically interacts with DNA and nucleosomes to condense and stabilize DNA into higher-order structures [3]. Nucleosome formation, histone eviction from DNA and chromatin remodeling, which shifts the distance between nucleosomes, allow for dynamic regulation of chromatin accessibility to gene expression mechanisms.
There are 118 histone genes: 36 H2A genes, 31 H2B genes, 23 H3 genes, 16 H4 genes, and 12 H1 genes, though some are believed to be pseudogenes not transcribed or translated [4,5]. Histone genes are primarily clustered at the HIST1 locus of chromosome 6 with the second-largest cluster on chromosome 1 (Figure 1A) [4,6]. Many histone genes encode proteins that diverge by a few amino acids, which may subtly adjust nucleosome stability or mobility and provide additional regulation for chromatin accessibility and gene expression [4,7].
Histones are frequently mutated in cancer and coincide with worse overall survival.
There are two classes of histones: canonical and variant [8]. Expression of canonical histones is orchestrated by promoter sequences which guide S-phase specific transcription to increase the histone supply necessary to complex with newly synthesized DNA [8,9]. Rather than being polyadenylated, canonical histone genes have a stem loop sequence and a histone downstream element (HDE) in the 3′UTR [10]. The stem loop and HDE are required for efficient translation and catabolism of the canonical histone RNA, restricting its expression to S phase [10]. Notably, canonical histones do not contain introns. Variant histones are expressed in a replication-independent manner throughout the cell cycle [11]. Variant histones may replace canonical histones in the nucleosome to confer distinct functional or structural features, impacting mechanisms such as gene expression and DNA repair [12]. Histones are subject to post-translation modifications that may affect histone stability, histone–DNA interaction, chromatin accessibility or binding of nucleosome interacting proteins to regulate replication, transcription, or repair. Consequently, patterns of histone post-translational modifications create chromatin states associated with either gene expression or gene repression [13–19]. Within the past decade, canonical and variant histone mutations have been revealed as a new class of epigenetic mutations in cancer. [20–26]. Here we catalog mutations of histone H2A, H2B, H3, H4, and H1 in cancer and summarize recent findings for how these mutations may contribute to tumorigenesis.
Histone mutation frequency in cancer
To determine the landscape of recurrent histone mutations in cancer, we surveyed 205 non-redundant cancer studies (65 489 patient samples) using the cBioPortal mirror hosted by the University of Florida Health Cancer Center, which allowed simultaneous search of all histone genes. After excluding 32 622 patient samples not profiled for every histone gene, we found genomic alterations (mutations, structural variants, copy number alterations) of H1 genes in 7% of patients, H2A genes in 11% of patients, H2B genes in 9% of patients, H3 genes in 10% of patients, and H4 genes in 6% of patients. Genomic alterations in any histone gene occurred with a high frequency in a wide range of cancers and were found in more than 50% of non-melanoma skin cancer patients, and up to 30% of patients with tobacco-related cancers, such as bladder urothelial carcinoma, head and neck cancer, esophagogastric adenocarcinoma, or non-small cell lung cancer (Figure 1B). Importantly, cancer patients carrying any histone genomic alteration had a significantly lower overall survival compared with cancer patients without a histone mutation (Figure 1C). In addition, patients carrying genomic alterations of either H2A, H2B, H3, or H4 had a significantly decreased overall survival across all cancers.
About 85% of recorded mutations in each histone were missense mutations (Figure 1D). To identify which histone amino acids were most often found mutated in cancer patients, we combined cBioportal data from multiple gene paralogs whose amino acid sequences aligned well. We found that most histone consensus amino acid residues had a missense mutation in at least one patient, and on average each residue was mutated in 4–12 patient samples. A Z-score was calculated at each amino acid residue of each histone by subtracting the number of mutations at a particular amino acid from the mean number of mutations across all amino acids and dividing the result by the standard deviation across all amino acids. Any amino acid with a Z-score >2 was considered a significant hotspot of recurrent mutation (Figures 2 and 4). Hotspot mutations were found in both histone tails and the histone-fold domains (Figures 2 and 4). Many histone missense mutations caused a charge swap (E to K mutations) or removed a positive charge (R to C mutations). Mutational hotspots in tumors suggest gain or change of function that promotes tumorigenesis. Indeed, some of these recurrent histone mutations have been previously studied in cancer such as histone H3.3 K27 and H3.3 G34 mutations [20,21,25–27]. Our prior analysis of histone mutations found only two histone mutant hotspots, H2A S1 and R35, had single nucleotide polymorphisms (SNP) present with similar frequency in the SNP database (dbSNP) and cBioPortal, indicating that most of the hotspot mutations are not simply passenger SNPs [25]. Also, some histone residues are rarely found mutated, which may be due to carcinogen or sequence specificity of DNA mutations or reflect that these histone mutations have deleterious effects on cell growth.
Survey of the most frequent core histone mutations in cancer patient samples.
Analysis of H2B E76 mutations in cancer.
Survey of histone H1 mutations in cancer.
H2A histone mutations in cancer
Histone H2A missense mutations constituted ∼20% of histone missense mutations in cancer (Figure 1D). The most common mutated H2A residue, E121, demonstrated decreased nucleosome sliding via ATP-utilizing chromatin remodeling and assembly factor (ACF) when mutated in vitro and impaired yeast growth (Table 1) [27]. The second most common mutated H2A residue, R29, is a site of post-translational modification by methylation and a point of contact with nucleosomal DNA [2,28]. Mutations at this residue demonstrated increased nucleosome sliding via ACF in vitro, decreased nucleosome thermal stability, increased histone H2A/H2B dimer exchange via histone chaperone nuclear assembly protein-1 (Nap1), and impaired yeast growth (Table 1) [27]. Interestingly, many H2A missense mutations occur in the histone-fold domain or the acidic patch on the histone (Figure 2A). The acidic patch consists of six H2A (E56, E61, E64, D90, E91, E92) and two H2B (E105, E113) amino acids on the nucleosome surface that provide a charged surface for protein binding [29]. Mutations occurring on acidic patch residues increased Nap1-mediated histone H2A/H2B dimer exchange (H2A E56), decreased thermal stability (H2A E56), increased (H2A E56) or decreased (H2A D90, E92) nucleosome sliding via ACF or ISWI, and impaired yeast growth (H2A E56, D90, E91, E92) [27,30,31] (Table 1). Recently, H2A E92K, which reverses the amino acid charge, was found to decrease binding of DNMT3A and the deacetylase SIRT6 to the nucleosome while increasing binding of bromodomain proteins BRD2/3/4 [32]. Mutations in the histone-fold domain of H2A were found to increase nucleosome sliding (H2A K74, K75) and impair yeast growth (H2A D72, K75) [27] (Table 1).
Histone amino acid . | Function/phenotype of mutant . | References . |
---|---|---|
H2A | ||
G4 | 1 | [27] |
R11 | Unknown | |
R29 | 1,2,3,4 | [27] |
G37 | Unknown | |
E56 | 1,2,3,4 | [27] |
D72 | 1 | [27] |
K74 | 3 | [27] |
K75 | 1,3 | [27] |
D90 | 1,3 and inhibited chromatin remodeling by ISWI | [27,30,31] |
E91 | 1 | [27] |
E92 | 1,3 and inhibited chromatin remodeling by ISWI | [27,30,31] |
E121 | 1,3 | [27] |
H2B | ||
E2 | 1 | [27] |
S36 | 1 | [27] |
G53 | 1,3 and decreased DNA interaction, increased migration, increased RNA Pol II passage, altered gene expression. | [27,34,35,37] |
D68 | 1,2,4, unstable histone octamer, and impaired H2A.Z-H2B incorporation into chromatin | [27,78,89] |
F70 | 2 | [27] |
E71 | 1,2,4,5 | [27] |
E76 | 1,2,4,5 and nucleosome instability, increased H2A–H2B dimer mobility, altered gene expression, increased proliferation and colony formation, increased chromatin accessibility, increased interaction with NAP1L1 and NAP1L2 chaperones. | [24–27,33] |
E113 | 1,5 and increased chromatin remodeling by ISWI | [27,31] |
H3 | ||
R26 | Unknown | |
K27 | 5 and global loss & focal gain of H3K27me3 modification, suppression of PRC2, altered proliferation, increased H3K27ac histone modification, increased DNA hypomethylation | [20,21,39–45,90] |
G34 | 5 and increased colony formation, increased infiltration, increased proliferation, DNA hypomethylation, decreased H3K36 methylation in cis, increased H3K27 methylation. | [20,39,59,60,61,62,66,91] |
K36 | 1,5 and increased colony formation, suppression of SETD2, global decrease in H3K36 di- and trimethylation, tumor formation in vivo, decreased sensitivity to apoptotic stimuli, impaired homologous recombination | [27,62,66–68] |
E50 | 1,2 | [27] |
E73 | 1 | [27] |
E97 | 1,2,4 and enhanced colony formation, nucleosome instability | [24,27,70] |
E105 | 1 and yeast sin mutant | [27,71] |
H4 | ||
G42 | 1,2,3 | [27] |
R45 | 1,3 and Increased resistance to UV damage, increased nucleotide excision repair, yeast sin mutant, prevents Mg2+ dependent chromatin condensation, loss of interaction with nucleosome DNA, decreased RNA Pol II pausing | [7,27,74–76,77,92] |
E63 | 1 | [27] |
D68 | 1,2,4 | [27] |
R92 | 1,6 | [27,73] |
Histone amino acid . | Function/phenotype of mutant . | References . |
---|---|---|
H2A | ||
G4 | 1 | [27] |
R11 | Unknown | |
R29 | 1,2,3,4 | [27] |
G37 | Unknown | |
E56 | 1,2,3,4 | [27] |
D72 | 1 | [27] |
K74 | 3 | [27] |
K75 | 1,3 | [27] |
D90 | 1,3 and inhibited chromatin remodeling by ISWI | [27,30,31] |
E91 | 1 | [27] |
E92 | 1,3 and inhibited chromatin remodeling by ISWI | [27,30,31] |
E121 | 1,3 | [27] |
H2B | ||
E2 | 1 | [27] |
S36 | 1 | [27] |
G53 | 1,3 and decreased DNA interaction, increased migration, increased RNA Pol II passage, altered gene expression. | [27,34,35,37] |
D68 | 1,2,4, unstable histone octamer, and impaired H2A.Z-H2B incorporation into chromatin | [27,78,89] |
F70 | 2 | [27] |
E71 | 1,2,4,5 | [27] |
E76 | 1,2,4,5 and nucleosome instability, increased H2A–H2B dimer mobility, altered gene expression, increased proliferation and colony formation, increased chromatin accessibility, increased interaction with NAP1L1 and NAP1L2 chaperones. | [24–27,33] |
E113 | 1,5 and increased chromatin remodeling by ISWI | [27,31] |
H3 | ||
R26 | Unknown | |
K27 | 5 and global loss & focal gain of H3K27me3 modification, suppression of PRC2, altered proliferation, increased H3K27ac histone modification, increased DNA hypomethylation | [20,21,39–45,90] |
G34 | 5 and increased colony formation, increased infiltration, increased proliferation, DNA hypomethylation, decreased H3K36 methylation in cis, increased H3K27 methylation. | [20,39,59,60,61,62,66,91] |
K36 | 1,5 and increased colony formation, suppression of SETD2, global decrease in H3K36 di- and trimethylation, tumor formation in vivo, decreased sensitivity to apoptotic stimuli, impaired homologous recombination | [27,62,66–68] |
E50 | 1,2 | [27] |
E73 | 1 | [27] |
E97 | 1,2,4 and enhanced colony formation, nucleosome instability | [24,27,70] |
E105 | 1 and yeast sin mutant | [27,71] |
H4 | ||
G42 | 1,2,3 | [27] |
R45 | 1,3 and Increased resistance to UV damage, increased nucleotide excision repair, yeast sin mutant, prevents Mg2+ dependent chromatin condensation, loss of interaction with nucleosome DNA, decreased RNA Pol II pausing | [7,27,74–76,77,92] |
E63 | 1 | [27] |
D68 | 1,2,4 | [27] |
R92 | 1,6 | [27,73] |
Phenotype key: .
Impaired yeast growth;
Increased Nap1-mediated histone dimer exchange;
Altered nucleosome sliding;
Decreased thermal stability;
Impaired differentiation of cells;
Hypersensitivity to DNA damaging agents.
H2B histone mutations in cancer
Histone H2B mutations account for ∼21% of all histone missense mutants in cancer samples (Figure 1D). Several H2B mutations have been identified for their roles in destabilizing the nucleosome or weakening histone–DNA interactions [24,25,27,33–35]. The most frequent H2B mutation occurs at residue 76, in which glutamate is often mutated to a lysine or, less commonly, to a glutamine (Figure 3A). In a pan-cancer analysis using data from cBioPortal, we found 41/54 patient samples with a glutamate to lysine mutation (H2B E76K), and 10/54 patient samples with a glutamate to glutamine mutation (H2B E76Q) (Figure 3A). While it is unknown why these missense mutations are most common, only a single nucleotide substitution would change glutamate to lysine (G/C > A/T) or glutamine (G/C > C/G). In our cBioPortal analysis, the most common single base substitutions for the H2B-E76K mutation were C > T or G > A transition mutations with about half carrying a C > T mutation, a signature cancer-related substitution [36]. H2B E76 mutations were distributed across 16 of 22 canonical H2B genes surveyed (Figure 3B). The variant allele frequency (VAF) of H2B E76 mutations in cancer was ∼20%, suggesting these are subclonal mutations contributing to tumor progression rather than initiating mutations (Figure 3C). H2B E76 mutations occur in numerous cancer types such as bladder (24%) and lung (24%) cancers (Figure 3D). We examined whether the top 20 genomic alterations found in bladder or lung cancers coincide with the H2B E76 mutation in patient samples with these cancers. In bladder cancer patient samples, H2B E76 mutations significantly co-occur with genomic alterations in RB1, KMT2D, ERBB2, and ERCC2, while in lung cancer patient samples the E76 mutation coincides with TP53, PIK3CA, KEAP1, ARID1A, KMT2C, and CDKN2A genomic alterations and these mutations may synergize to propel tumor progression (Figure 3E).
X-ray crystallography demonstrated that the H2B E76K mutation repels the H4 R92 residue, disrupting an important H2B/H4 interaction point and destabilizing the histone octamer [24–26]. Thus, H2B E76K is not readily able to form a histone octamer without DNA present [24,25,27]. H2A/H2B dimer exchange is increased in cells expressing H2B E76K and is dependent on the histone chaperone Nap1 (nuclear assembly protein-1) [24,25,27] (Table 1). Others have reported increased interaction between H2B E76K histones and chaperone proteins Nap1l1 and Nap1l2 that may contribute to increased H2A/H2B dimer exchange in H2B E76K expressing cells [27,33]. Unlike other histone mutants, H2B E76K does not affect global histone post-translational modification levels [33].
The H2B E76K oncohistone has been reported to promote chromatin accessibility at promoters and enhancers [25,33]. Multiple studies report that H2B E76K alters gene expression [25,27,33]. These changes in gene expression translate to phenotypic changes such as increased proliferation, increased colony formation and inhibition of differentiation that is observed in cells expressing H2B E76K [24,25,27,33]. Why only certain sets of genes are activated is not certain, but it has been reported that H2B E76K leads to enhanced transcriptional elongation rates, which might favor expression of subsets of genes.
The second most mutated histone H2B residue, H2B E113, is part of the nucleosome acidic patch. Previous studies have demonstrated that H2B E113 is most often mutated to lysine (H2B E113K), which causes increased nucleosome remodeling via ACF or ISWI [27,31]. Like other H2B mutants, H2B E113K has been reported to inhibit differentiation of murine mesenchymal cells and impair yeast growth [27] (Table 1).
Another well-characterized H2B oncohistone is H2B G53D originally found in pancreatic ductal adenocarcinoma (PDAC), glioblastoma multiforme, and lung cancer [35,37]. This residue lies close to the DNA entry/exit point on the nucleosome. Nucleosome pulling and holding assays suggest that the H2B G53D acquired negative charge may interfere with histone–DNA interactions [35]. Though DNA replication is not affected, H2B G53D expressing cells have increased RNA polymerase II passage through mutant nucleosomes [35]. Consistent with a weakened histone–DNA interaction, H2B G53D nucleosomes display increased nucleosome sliding via ACF in vitro [27] (Table 1). Expression of H2B G53D also causes gene expression changes, affecting pathways such as growth factor signaling, cell adhesion and cancer development [34] (Table 1). Based on analysis of CUT&RUN, gene expression, and PRO-seq data from H2B G53D expressing cells, H2B G53D increases active RNA-polymerase II and increased gene expression at a subset of genes [34]. One oncogene found in this subset, ANXA3, was demonstrated to have increased transcription due to H2B G53D, and ANXA3 expression was greater in cancer patients carrying the H2B G53D mutation, supporting the notion that expression of H2B G53D drives expression of the ANXA3 oncogene [34]. In addition to chromatin and gene expression changes, the H2B G53D mutation also promotes gap closure and transwell migration in a pancreatic cell line [35], but impairs yeast growth [27] (Table 1).
A common phenotype for multiple histone H2B mutations such as amino acid E2, S36, G53, D68, E71, E76, or E113, is impaired yeast growth [27] (Table 1). How this translates to human oncogenic phenotypes is not yet certain. In addition, mutations at H2B D68 destabilize the histone octamer, like the H2B E76K nucleosome [25,27]. Furthermore, histone H2A/H2B dimer exchange via Nap1 histone chaperone was also found to be increased when residues H2B D68, F70, or E71 were mutated in vitro [27]. Expression of H2B E71K inhibits differentiation of murine mesenchymal cells, as observed with many other histone mutants [27].
Histone H3 mutations in cancer
Almost 25% of histone missense mutations in the cBioPortal are in Histone H3 (Figure 1D). The first oncohistones identified were mutations in genes H3F3A and H3F3B that encode variant histone H3.3. Recurrent H3.3 mutations affecting two critical amino acids (K27 and G34) have been identified in one-third of pediatric glioblastoma (GBM) and nearly 80% of diffuse intrinsic pontine gliomas (DIPG). Our cBioPortal analysis indicates that the most frequently mutated histone H3 residue is K27, usually in H3.3 (Figure 2C). This amino acid on the N-terminal tail of H3 can be post-translationally modified by acetylation, methylation, ubiquitination, and ADP-ribosylation [38]. H3.3 K27 is commonly mutated to a methionine (H3 K27M), which can no longer be methylated or acetylated, leading to decreased global H3 K27me2/me3 and increased H3 K27 acetylation levels [39–44]. Polycomb repressive complex 2 (PRC2) methylates H3 K27 via its enhancer of zeste homolog 2 (EZH2) catalytic subunit. Initial reports suggested PRC2 became sequestered on the K27M mutant histone, leading to a global decrease in H3K27me3 [39,41,45]. Such experiments showed enhanced co-precipitation of H3 K27M with PRC2 and increased affinity of H3 K27M for PRC2 compared with WT H3 [39,41]. Additionally, H3 K27M was found to have decreased association with corepressor protein CYL, which normally connects PRC2 with neighboring nucleosomes to promote H3K27me3 spreading [46,47]. However, recent reports offer differing results, highlighting that the mechanism of action of H3 K27M may not fully be elucidated. In contrast with the sequestration model: (1) H3K27me3 marks are focally retained in chromatin, primarily on CpG islands and a subset of genes such as the tumor suppressor p16/INK4A [41,42,48–51]; (2) PRC2 has not consistently been shown to have greater affinity for H3 K27M over H3 WT [50,52,53]; (3) H3 K27M does not always co-localize with PRC2 on chromatin, in some cases PRC2 has been excluded [54–56]; (4) H3K27ac marks are enriched on H3 K27M nucleosomes, which further exclude PRC2 [55]; and (5) there is a persistent decrease in PRC2 activity following its transient interaction with H3 K27M [41,50]. Collectively these data indicate that H3K27M disrupts normal patterns of histone methylation, with large regions of H3K27me3 loss but other areas of focal gain which together lead to disruption of normal gene expression contributing to oncogenesis. Accordingly, H3 K27M increases proliferation of cells and the combination of H3 K27M, p53 knockdown, and oncogenic PDGFRA inhibited differentiation of astrocytes [38]. Deposition of K27M at active enhancers of neurodevelopment genes leads to H3K27ac loss, decreased chromatin accessibility and reduced expression of neurodevelopmental genes in neural stem cells [57]. In addition, because DNMT3A/B DNA methyltransferases associate with EZH2, K27M mutant tumors have a substantially different DNA methylation profile from H3.3 WT pediatric high-grade gliomas, with an overall reduction in global DNA methylation and enrichment of hypomethylated regions [41,58].
Another H3 oncohistone is at amino acid G34 on the N-terminal tail. In pediatric glioblastoma, H3.3 G34 is frequently mutated to arginine or valine (H3.3 G34R/V) [20]. In giant cell tumors of bone (GCTB) this residue is most often mutated to tryptophan (H3.3 G34W) [59]. Primary cells derived from GCTB biopsies displayed increased proliferation, invasion, and colony formation [60]. In addition, H3.1 G34 mutations are found in carcinomas [26]. Although not a target of post-translational modification itself, mutations at G34 disrupt histone H3 K36 and H3 K27 methylation in cis (on the same histone tail) [39,59,61,62]. Mutations changing glycine 34 to a bulkier residue have been reported to block histone tail binding by SETD2 and H3K36 methylation [63]. Despite the inhibitory effect of H3.3 G34W on SETD2, genome-wide studies indicate more areas of H3K36me3 gain than loss [64]. Furthermore, H3 G34W has been associated with global DNA hypomethylation as well as specific focal DNA hypermethylation that impairs differentiation in vitro [59]. Mass spectrometry has identified that H3.3 G34W binding partners include SWI/SNF components ARID2, SMARCA1, and HDAC2 [60]. Thus, G34W may recruit proteins from other sites in the genome to alter local chromatin architecture as well as alter the accumulation of H3.3 through aberrant interaction with histone chaperones. Furthermore, in GCTB cell lines H3.3 G34W made up only 3% of histone H3, yet H3K27me3 levels on non-mutant histones H3.3 and H3.1 decreased by 30% to 40%, suggesting an important trans effect of the mutant histone [64]. In the presence of G34W, H3K27me3 and PRC2 were significantly depleted from intergenic regions, and this may play a role in the observed chromatin and gene expression changes observed in H3.3 G34W mutant tumors [64]. In addition, a recent study revealed that H3.3 G34R/V/W germline mutations have distinct and substantial impacts on genome-wide distribution of epigenetic marks and interacting proteins that may account for the phenotypic diversity associated with these mutants [65]. Specifically, G34R greatly depleted H3K36me2/3 in cis but had no obvious effect on H3K27 methylation [65]. In contrast, G34V/W significantly increased H3K27me3 but had little effect on H3K36me2/3 [65]. Furthermore, the G34R-mediated redistribution of H3K36me2 altered the recruitment and redistribution of DNMT3A, which in turn changed the landscape of DNA methylation, leading to deregulated expression of complement genes and suppression of neuronal developmental genes [65].
About one year after the discoveries of the H3.3 K27 and H3.3 G34 oncohistones, the H3.3 K36M N-terminal tail mutation was discovered in chondroblastoma patient samples [66]. H3 K36 is a site of post-translational modification by acetylation, methylation or ubiquitylation [38]. H3 K36M causes a global decrease in H3 K36 di- and trimethylation levels [62,67,68] at least in part by trapping NSD enzymes and suppressing SETD2, the histone methyltransferases responsible for H3 K36 methylation [62,69]. Cells expressing H3.3 K36M displayed increased colony formation and resistance to apoptosis [68]. The H3 K36M mutation impaired adipocyte and osteocyte differentiation of murine mesenchymal cells, and mesenchymal progenitor cells expressing H3 K36M formed an undifferentiated sarcoma after injection into mice [27,67,68].
Though not as thoroughly studied as H3.3 K27, G34, and K36 mutations, H3.1 mutations at amino acids E97 and E105 were found mutated at a high frequency in our pan-cancer analysis (Figure 2C). Both H3 E97 and E105 residues are in the histone-fold domain, as opposed to the N-terminal tail, and are often found mutated to lysine from our cBioPortal results (Figure 2C). H3 E97 mutations increased Nap1-mediated histone H2A/H2B dimer exchange, destabilized the nucleosome, and decreased thermal stability [24,27]. Mutations at both H3 E97 and E105 impair yeast growth, and H3 E97 mutants also increase colony formation [24,27,70]. In yeast, the H3 E105K causes very little change in nucleosome structure and stability yet causes SWI/SNF independence (SIN) i.e. this mutation alleviates the transcriptional defects caused by inactivation of the yeast SWI/SNF complex [71].
Histone H4 mutations in cancer
The only histone to have a single variant form [72], histone H4 binds with histone H3 to form the crucial H3/H4 tetramer and histone H4 interacts with histone H2B to form the histone octamer [1]. The most frequent histone H4 missense mutations in cancer patient samples were located at amino acids R45 and D68. Other, less frequent histone H4 mutations at G42, E63, and R92 were found to impair yeast growth while mutation of H4 at R92 has been reported to increase sensitivity to DNA damaging agents [27,73] (Table 1).
H4 R45 is inserted into the minor groove of the nucleosomal DNA, interacting with DNA via van der Waals forces rather than site specific contacts [1,74]. Based on our cBioPortal search, histone H4 R45 is most often mutated to a glutamine or cysteine. Mutations to this residue has been shown to increase nucleosome sliding, decrease RNA polymerase II pausing, and alter interactions with nucleosomal DNA [27,74,75] (Table 1). In yeast, the H4 R45C mutation functions as a SIN mutant yet still impaired yeast growth [76]. In addition, H4 R45C increased nucleotide excision repair and resistance to UV damage in yeast [27,76]. The R45C mutation does not disrupt nucleosome positioning or the accessibility of nucleosomal DNA but rather has been reported to eliminate Mg2+-dependent, intramolecular folding of nucleosomal arrays, suggesting that SIN histone mutants may alleviate the need for SWI/SNF in vivo by disrupting higher-order chromatin folding [77].
The most frequently mutated amino acid of histone H4 is located at D68 (Figure 2D). This residue forms a stable electrostatic interaction between histones H2B E76 and H4 [24–26]. Histone H4 D68 is commonly mutated to an asparagine and these mutations have been reported to increase Nap1-mediated histone H2A/H2B dimer exchange in vitro [27]. In addition, H4 mutations at D68 impair binding of histone H2A.Z-H2B dimer to chromatin [78].
Histone H1 mutations in cancer
The linker histone H1 binds to DNA exiting the nucleosome core to catalyze condensation of nucleosomal arrays [79]. H1 is a larger, highly mobile histone that interacts with DNA in a more rapid fashion than core histones comprising the nucleosome [79]. Similar to other histones, histone H1 is expressed from multiple genes, though the sequences are less conserved [79]. In our cBioPortal pan-cancer analysis, we found missense mutations in 10 of 12 histone H1 genes (Figure 4A). Histone H1 missense mutations occurred in a broad range of cancers, including skin, prostrate and endometrial cancers as well as mature B-cell neoplasms (Figure 4B). H1 missense mutations were detected in cancer patient samples at nearly every residue, and the two most common ‘hotspot’ mutations were observed at A64/67 and A163/164 (Figure 4A). The two most common histone H1 mutations were overwhelmingly observed in mature B-cell neoplasms (Figure 4C), a cancer type previously reported to carry histone H1 mutations and deletions [22,80–85].
Histone H1 may have a tumor suppressor function. Loss of histone H1 was associated with chromatin decompaction as well as increased H3 K36 methylation and decreased H3 K27 methylation, causing B to A compartment shifts and transcriptional activation [81,86,87]. For instance, genes related to T-cell activation were expressed with H1 loss [81]. Also, H1 depletion by deletion of three H1 genes in pluripotent embryonic stem cells blocked differentiation due to unchecked expression of pluripotency genes [88]. Germinal center B cells may be especially sensitive to H1 dosage or mutation due to their rapid proliferation rate [81]. In Diffuse Large B-cell Lymphoma (DLBCL), samples carrying H1 mutations have been reported to display an increased mutational burden [84]. In an analysis of mutations driving follicular lymphoma initiation and progression, mutations affecting at least one histone H1 gene were observed in 28% of samples with most being missense mutations in the globular domains of H1, H1–2, and H1–4 genes [80]. When either WT or a S101F mutation of histone H1–2 was expressed in H1 knockout murine embryonic stem cells, the S101F mutant was not associated with chromatin to the same degree as wild-type histone, suggesting this is a loss of function mutation that reduces the affinity of H1 with chromatin [80]. In our pan-cancer analysis, histone H1–2 residues G102 and S103 on H1–2 are more often found mutated than S101 (Figure 4A), and these mutations may disrupt a similar region of H1 as S101 necessary for DNA interaction. When GFP-tagged S101F, G102A, and S103F were expressed in 3T3 cells more rapid fluorescence recovery after photobleaching was observed than in cells expressing GFP-tagged wild-type H1, and the dissociation constant of H1 S101F was much higher than wild-type H1 in a mononucleosome binding assay [81]. Thus, H1 globular domain mutations may impair the association of H1 with chromatin, while C-terminal mutants remain able to bind chromatin normally but disrupt compaction [81].
Summary
Histone mutations are pervasive in all cancer types. They may disrupt histone post-translational modifications, protein interactions with the nucleosome acid patch, nucleosome stability, remodeling, or chromatin compaction (Figure 5). Oncohistone mutations disrupt chromatin function to cause a state of gene expression favoring tumorigenesis. Many facets of histone mutations in cancer remain to be elucidated. For instance, histones are encoded by multiple genes with small sequence variations of unknown significance that may provide functional diversity. Various tissues have distinct histone gene expression patterns, but the importance of these patterns is unknown. Explicit histone variants may be required by specific cell types during development or other cellular processes. Histone mutations are rarely found on their own in cancer patient samples but instead often co-occur with well-known oncogenes. It isn't yet clear how histone mutations and destabilization of chromatin structure or function may impact or enable oncogene activity.
Histone mutations may alter chromatin structure and function by diverse mechanisms.
Perspectives
Histone mutations are frequently found in patient samples of all cancer types and associate with worse survival. Yet how these mutations may affect chromatin structure and regulation of gene expression is unknown.
Oncohistone mutations may disrupt nucleosome stability, histone post-translational modifications, or chromatin dynamics to undermine gene regulatory mechanisms, resulting in gene expression programs that favor oncogenesis.
As the functions of oncohistone mutations become better understood, we will gain further insights into how loss of nucleosomal mediated gene repression contributes to tumorigenesis and whether this may be countered therapeutically.
Competing Interests
The authors declare that there are no competing interests associated with the manuscript.
Funding
R01CA266078 and the UF Health Cancer Center, supported by state appropriations provided in Florida Statute § 381.915 and P30CA247796.
Open Access
Open access for this article was enabled by the participation of University of Florida in an all-inclusive Read & Publish agreement with Portland Press and the Biochemical Society under a transformative agreement with Individual.
Author Contributions
K.N.E.P., J.S. and R.L.B. conceived and wrote the manuscript as well as created the figures, J.D.L. edited the manuscript and contributed to the design and direction of the work.