Convergent Clonal Evolution of Signaling Gene Mutations Is a Hallmark of Myelodysplastic Syndrome Progression.

Progression from myelodysplastic syndromes (MDS) to secondary acute myeloid leukemia (AML) is associated with the acquisition and expansion of subclones. Our understanding of subclone evolution during progression, including the frequency and preferred order of gene mutation acquisition, remains incomplete. Sequencing of 43 paired MDS and secondary AML samples identified at least one signaling gene mutation in 44% of MDS and 60% of secondary AML samples, often below the level of standard sequencing detection. In addition, 19% of MDS and 47% of secondary AML patients harbored more than one signaling gene mutation, almost always in separate, coexisting subclones. Signaling gene mutations demonstrated diverse patterns of clonal evolution during disease progression, including acquisition, expansion, persistence, and loss of mutations, with multiple patterns often coexisting in the same patient. Multivariate analysis revealed that MDS patients who had a signaling gene mutation had a higher risk of AML progression, potentially providing a biomarker for progression. SIGNIFICANCE: Subclone expansion is a hallmark of progression from MDS to secondary AML. Subclonal signaling gene mutations are common at MDS (often at low levels), show complex and convergent patterns of clonal evolution, and are associated with future progression to secondary AML. See related article by Guess et al., p. 316 (33). See related commentary by Romine and van Galen, p. 270. This article is highlighted in the In This Issue feature, p. 265.

Kdm6a deficiency restricted to mouse hematopoietic cells causes an age- and sex-dependent myelodysplastic syndrome-like phenotype.

Kdm6a/Utx, a gene on the X chromosome, encodes a histone H3K27me3 demethylase that has an orthologue on the Y chromosome (Uty) (Zheng et al. 2018). We previously identified inactivating mutations of Kdm6a in approximately 50% of mouse acute promyelocytic leukemia samples; however, somatic mutations of KDM6A are more rare in human AML samples, ranging in frequency from 2-15% in different series of patients, where their role in pathogenesis is not yet clear. In this study, we show that female Kdm6aflox/flox mice (with allele inactivation initiated by Vav1-Cre in hematopoietic stem and progenitor cells (HSPCs) have a sex-specific phenotype that emerges with aging, with features resembling a myelodysplastic syndrome (MDS). Female Kdm6a-knockout (KO) mice have an age-dependent expansion of their HSPCs with aberrant self-renewal, but they did not differentiate normally into downstream progeny. These mice became mildly anemic and thrombocytopenic, but did not develop overt leukemia, or die from these cytopenias. ChIP-seq and ATAC-seq studies showed only minor changes in H3K27me3, H3K27ac, H3K4me, H3K4me3 and chromatin accessibility between Kdm6a-WT and Kdm6a-KO mice. Utilizing scRNA-seq, Kdm6a loss was linked to the transcriptional repression of genes that mediate hematopoietic cell fate determination. These data demonstrate that Kdm6a plays an important role in normal hematopoiesis, and that its inactivation may contribute to AML pathogenesis.

Unexpected suppression of tumorigenesis by c-MYC via TFAP4-dependent restriction of stemness in B lymphocytes.

The proliferative burst of B lymphocytes is essential for antigen receptor repertoire diversification during the development and selective expansion of antigen-specific clones during immune responses. High proliferative activity inevitably promotes oncogenesis, the risk of which is further elevated in B lymphocytes by endogenous gene rearrangement and somatic mutations. However, B-cell-derived cancers are rare, perhaps owing to putative intrinsic tumor-suppressive mechanisms. We show that c-MYC facilitates B-cell proliferation as a protumorigenic driver and unexpectedly coengages counteracting tumor suppression through its downstream factor TFAP4. TFAP4 is mutated in human lymphoid malignancies, particularly in >10% of Burkitt lymphomas, and reduced TFAP4 expression was associated with poor survival of patients with MYC-high B-cell acute lymphoblastic leukemia. In mice, insufficient TFAP4 expression accelerated c-MYC-driven transformation of B cells. Mechanistically, c-MYC suppresses the stemness of developing B cells by inducing TFAP4 and restricting self-renewal of proliferating B cells. Thus, the pursuant transcription factor cascade functions as a tumor suppressor module that safeguards against the transformation of developing B cells.

Remethylation of Dnmt3a-/- hematopoietic cells is associated with partial correction of gene dysregulation and reduced myeloid skewing.

Mutations in the DNA methyltransferase 3A (DNMT3A) gene are the most common cause of age-related clonal hematopoiesis (ARCH) in older individuals, and are among the most common initiating events for acute myeloid leukemia (AML). The most frequent DNMT3A mutation in AML patients (R882H) encodes a dominant-negative protein that reduces methyltransferase activity by ∼80% in cells with heterozygous mutations, causing a focal, canonical DNA hypomethylation phenotype; this phenotype is partially recapitulated in murine Dnmt3a-/- bone marrow cells. To determine whether the hypomethylation phenotype of Dnmt3a-/- hematopoietic cells is reversible, we developed an inducible transgene to restore expression of DNMT3A in transplanted bone marrow cells from Dnmt3a-/- mice. Partial remethylation was detected within 1 wk, but near-complete remethylation required 6 mo. Remethylation was accurate, dynamic, and highly ordered, suggesting that differentially methylated regions have unique properties that may be relevant for their functions. Importantly, 22 wk of DNMT3A addback partially corrected dysregulated gene expression, and mitigated the expansion of myeloid cells. These data show that restoring DNMT3A expression can alter the epigenetic "state" created by loss of Dnmt3a activity; this genetic proof-of-concept experiment suggests that this approach could be relevant for patients with ARCH or AML caused by loss-of-function DNMT3A mutations.

Mutation Clearance after Transplantation for Myelodysplastic Syndrome.

BACKGROUND: Allogeneic hematopoietic stem-cell transplantation is the only curative treatment for patients with myelodysplastic syndrome (MDS). The molecular predictors of disease progression after transplantation are unclear. METHODS: We sequenced bone marrow and skin samples from 90 adults with MDS who underwent allogeneic hematopoietic stem-cell transplantation after a myeloablative or reduced-intensity conditioning regimen. We detected mutations before transplantation using enhanced exome sequencing, and we evaluated mutation clearance by using error-corrected sequencing to genotype mutations in bone marrow samples obtained 30 days after transplantation. In this exploratory study, we evaluated the association of a mutation detected after transplantation with disease progression and survival. RESULTS: Sequencing identified at least one validated somatic mutation before transplantation in 86 of 90 patients (96%); 32 of these patients (37%) had at least one mutation with a maximum variant allele frequency of at least 0.5% (equivalent to 1 heterozygous mutant cell in 100 cells) 30 days after transplantation. Patients with disease progression had mutations with a higher maximum variant allele frequency at 30 days than those who did not (median maximum variant allele frequency, 0.9% vs. 0%; P<0.001). The presence of at least one mutation with a variant allele frequency of at least 0.5% at day 30 was associated with a higher risk of progression (53.1% vs. 13.0%; conditioning regimen-adjusted hazard ratio, 3.86; 95% confidence interval [CI], 1.96 to 7.62; P<0.001) and a lower 1-year rate of progression-free survival than the absence of such a mutation (31.3% vs. 59.3%; conditioning regimen-adjusted hazard ratio for progression or death, 2.22; 95% CI, 1.32 to 3.73; P=0.005). The rate of progression-free survival was lower among patients who had received a reduced-intensity conditioning regimen and had at least one persistent mutation with a variant allele frequency of at least 0.5% at day 30 than among patients with other combinations of conditioning regimen and mutation status (P≤0.001). Multivariate analysis confirmed that patients who had a mutation with a variant allele frequency of at least 0.5% detected at day 30 had a higher risk of progression (hazard ratio, 4.48; 95% CI, 2.21 to 9.08; P<0.001) and a lower 1-year rate of progression-free survival than those who did not (hazard ratio for progression or death, 2.39; 95% CI, 1.40 to 4.09; P=0.002). CONCLUSIONS: The risk of disease progression was higher among patients with MDS in whom persistent disease-associated mutations were detected in the bone marrow 30 days after transplantation than among those in whom these mutations were not detected. (Funded by the Leukemia and Lymphoma Society and others.).

Subclones dominate at MDS progression following allogeneic hematopoietic cell transplant.

Allogeneic hematopoietic cell transplantation (alloHCT) is a potentially curative treatment for myelodysplastic syndromes (MDS), but patients who relapse after transplant have poor outcomes. In order to understand the contribution of tumor clonal evolution to disease progression,we applied exome and error-corrected targeted sequencing coupled with copy number analysis to comprehensively define changes in the clonal architecture of MDS in response to therapy using 51 serially acquired tumor samples from 9 patients who progressed after an alloHCT. We show that small subclones before alloHCT can drive progression after alloHCT. Notably, at least one subclone expanded or emerged at progression in all patients. Newly acquired structural variants (SVs) were present in an emergent/expanding subclone in 8 of 9 patients at progression, implicating the acquisition of SVs as important late subclonal progression events. In addition, pretransplant therapy with azacitidine likely influenced the mutation spectrum and evolution of emergent subclones after alloHCT. Although subclone evolution is common, founding clone mutations are always present at progression and could be detected in the bone marrow as early as 30 and/or 100 days after alloHCT in 6 of 8 (75%) patients, often prior to clinical progression. In conclusion, MDS progression after alloHCT is characterized by subclonal expansion and evolution, which can be influenced by pretransplant therapy.

Haploinsufficiency for DNA methyltransferase 3A predisposes hematopoietic cells to myeloid malignancies.

The gene that encodes de novo DNA methyltransferase 3A (DNMT3A) is frequently mutated in acute myeloid leukemia genomes. Point mutations at position R882 have been shown to cause a dominant negative loss of DNMT3A methylation activity, but 15% of DNMT3A mutations are predicted to produce truncated proteins that could either have dominant negative activities or cause loss of function and haploinsufficiency. Here, we demonstrate that 3 of these mutants produce truncated, inactive proteins that do not dimerize with WT DNMT3A, strongly supporting the haploinsufficiency hypothesis. We therefore evaluated hematopoiesis in mice heterozygous for a constitutive null Dnmt3a mutation. With no other manipulations, Dnmt3a+/- mice developed myeloid skewing over time, and their hematopoietic stem/progenitor cells exhibited a long-term competitive transplantation advantage. Dnmt3a+/- mice also spontaneously developed transplantable myeloid malignancies after a long latent period, and 3 of 12 tumors tested had cooperating mutations in the Ras/MAPK pathway. The residual Dnmt3a allele was neither mutated nor downregulated in these tumors. The bone marrow cells of Dnmt3a+/- mice had a subtle but statistically significant DNA hypomethylation phenotype that was not associated with gene dysregulation. These data demonstrate that haploinsufficiency for Dnmt3a alters hematopoiesis and predisposes mice (and probably humans) to myeloid malignancies by a mechanism that is not yet clear.

Association Between Mutation Clearance After Induction Therapy and Outcomes in Acute Myeloid Leukemia.

IMPORTANCE: Tests that predict outcomes for patients with acute myeloid leukemia (AML) are imprecise, especially for those with intermediate risk AML. OBJECTIVES: To determine whether genomic approaches can provide novel prognostic information for adult patients with de novo AML. DESIGN, SETTING, AND PARTICIPANTS: Whole-genome or exome sequencing was performed on samples obtained at disease presentation from 71 patients with AML (mean age, 50.8 years) treated with standard induction chemotherapy at a single site starting in March 2002, with follow-up through January 2015. In addition, deep digital sequencing was performed on paired diagnosis and remission samples from 50 patients (including 32 with intermediate-risk AML), approximately 30 days after successful induction therapy. Twenty-five of the 50 were from the cohort of 71 patients, and 25 were new, additional cases. EXPOSURES: Whole-genome or exome sequencing and targeted deep sequencing. Risk of identification based on genetic data. MAIN OUTCOMES AND MEASURES: Mutation patterns (including clearance of leukemia-associated variants after chemotherapy) and their association with event-free survival and overall survival. RESULTS: Analysis of comprehensive genomic data from the 71 patients did not improve outcome assessment over current standard-of-care metrics. In an analysis of 50 patients with both presentation and documented remission samples, 24 (48%) had persistent leukemia-associated mutations in at least 5% of bone marrow cells at remission. The 24 with persistent mutations had significantly reduced event-free and overall survival vs the 26 who cleared all mutations. Patients with intermediate cytogenetic risk profiles had similar findings. [table: see text]. CONCLUSIONS AND RELEVANCE: The detection of persistent leukemia-associated mutations in at least 5% of bone marrow cells in day 30 remission samples was associated with a significantly increased risk of relapse, and reduced overall survival. These data suggest that this genomic approach may improve risk stratification for patients with AML.

A comprehensive and high-resolution genome-wide response of p53 to stress.

Tumor suppressor p53 regulates transcription of stress-response genes. Many p53 targets remain undiscovered because of uncertainty as to where p53 binds in the genome and the fact that few genes reside near p53-bound recognition elements (REs). Using chromatin immunoprecipitation followed by exonuclease treatment (ChIP-exo), we associated p53 with 2,183 unsplit REs. REs were positionally constrained with other REs and other regulatory elements, which may reflect structurally organized p53 interactions. Surprisingly, stress resulted in increased occupancy of transcription factor IIB (TFIIB) and RNA polymerase (Pol) II near REs, which was reduced when p53 was present. A subset associated with antisense RNA near stress-response genes. The combination of high-confidence locations for p53/REs, TFIIB/Pol II, and their changes in response to stress allowed us to identify 151 high-confidence p53-regulated genes, substantially increasing the number of p53 targets. These genes composed a large portion of a predefined DNA-damage stress-response network. Thus, p53 plays a comprehensive role in regulating the stress-response network, including regulating noncoding transcription.

Unusual combinatorial involvement of poly-A/T tracts in organizing genes and chromatin in Dictyostelium.

Dictyostelium discoideum is an amoebozoa that exists in both a free-living unicellular and a multicellular form. It is situated in a deep branch in the evolutionary tree and is particularly noteworthy in having a very A/T-rich genome. Dictyostelium provides an ideal system to examine the extreme to which nucleotide bias may be employed in organizing promoters, genes, and nucleosomes across a genome. We find that Dictyostelium genes are demarcated precisely at their 5' ends by poly-T tracts and precisely at their 3' ends by poly-A tracts. These tracts are also associated with nucleosome-free regions and are embedded with precisely positioned TATA boxes. Homo- and heteropolymeric tracts of A and T demarcate nucleosome border regions. Together, these findings reveal the presence of a variety of functionally distinct polymeric A/T elements. Strikingly, Dictyostelium chromatin may be organized in di-nucleosome units but is otherwise organized as in animals. This includes a +1 nucleosome in a position that predicts the presence of a paused RNA polymerase II. Indeed, we find a strong phylogenetic relationship between the presence of the NELF pausing factor and positioning of the +1 nucleosome. Pausing and +1 nucleosome positioning may have coevolved in animals.

Phylogenetic Profiles Reveal Structural and Functional Determinants of Lipid-binding.

One of the major challenges in the genomic era is annotating structure/function to the vast quantities of sequence information now available. Indeed, most of the protein sequence database lacks comprehensive annotation, even when experimental evidence exists. Further, within structurally resolved and functionally annotated protein domains, additional functionalities contained in these domains are not apparent. To add further complication, small changes in the amino-acid sequence can lead to profound changes in both structure and function, underscoring the need for rapid and reliable methods to analyze these types of data. Phylogenetic profiles provide a quantitative method that can relate the structural and functional properties of proteins, as well as their evolutionary relationships. Using all of the structurally resolved Src-Homology-2 (SH2) domains, we demonstrate that knowledge-bases can be used to create single-amino acid phylogenetic profiles which reliably annotate lipid-binding. Indeed, these measures isolate the known phosphotyrosine and hydrophobic pockets as integral to lipid-binding function. In addition, we determined that the SH2 domain of Tec family kinases bind to lipids with varying affinity and specificity. Simulating mutations in Bruton's tyrosine kinase (BTK) that cause X-Linked Agammaglobulinemia (XLA) predict that these mutations alter lipid-binding, which we confirm experimentally. In light of these results, we propose that XLA-causing mutations in the SH3-SH2 domain of BTK alter lipid-binding, which could play a causative role in the XLA-phenotype. Overall, our study suggests that the number of lipid-binding proteins is drastically underestimated and, with further development, phylogenetic profiles can provide a method for rapidly increasing the functional annotation of protein sequences.

Phylogenetic profiles reveal structural/functional determinants of TRPC3 signal-sensing antennae.

Biochemical assessment of channel structure/function is incredibly challenging. Developing computational tools that provide these data would enable translational research, accelerating mechanistic experimentation for the bench scientist studying ion channels. Starting with the premise that protein sequence encodes information about structure, function and evolution (SF&E), we developed a unified framework for inferring SF&E from sequence information using a knowledge-based approach. The Gestalt Domain Detection Algorithm-Basic Local Alignment Tool (GDDA-BLAST) provides phylogenetic profiles that can model, ab initio, SF&E relationships of biological sequences at the whole protein, single domain and single-amino acid level.1,2 In our recent paper,4 we have applied GDDA-BLAST analysis to study canonical TRP (TRPC) channels1 and empirically validated predicted lipid-binding and trafficking activities contained within the TRPC3 TRP_2 domain of unknown function. Overall, our in silico, in vitro, and in vivo experiments support a model in which TRPC3 has signal-sensing antennae which are adorned with lipid-binding, trafficking and calmodulin regulatory domains. In this Addendum, we correlate our functional domain analysis with the cryo-EM structure of TRPC3.3 In addition, we synthesize recent studies with our new findings to provide a refined model on the mechanism(s) of TRPC3 activation/deactivation.

Phylogenetic profiles reveal evolutionary relationships within the "twilight zone" of sequence similarity.

Inferring evolutionary relationships among highly divergent protein sequences is a daunting task. In particular, when pairwise sequence alignments between protein sequences fall <25% identity, the phylogenetic relationships among sequences cannot be estimated with statistical certainty. Here, we show that phylogenetic profiles generated with the Gestalt Domain Detection Algorithm-Basic Local Alignment Tool (GDDA-BLAST) are capable of deriving, ab initio, phylogenetic relationships for highly divergent proteins in a quantifiable and robust manner. Notably, the results from our computational case study of the highly divergent family of retroelements accord with previous estimates of their evolutionary relationships. Taken together, these data demonstrate that GDDA-BLAST provides an independent and powerful measure of evolutionary relationships that does not rely on potentially subjective sequence alignment. We demonstrate that evolutionary relationships can be measured with phylogenetic profiles, and therefore propose that these measurements can provide key insights into relationships among distantly related and/or rapidly evolving proteins.