Accurate detection of subclonal variants in paired diagnosis-relapse acute myeloid leukemia samples by next generation Duplex Sequencing.

Mutations characterize diverse human cancers; there is a positive correlation between elevated mutation frequency and tumor progression. One exception is acute myeloid leukemia (AML), which has few clonal single nucleotide mutations. We used highly sensitive and accurate Duplex Sequencing (DS) to show now that AML, in addition, has an extensive repertoire of variants with low allele frequencies, < 1%, which is below the accurate detection limit of most other sequencing methodologies. The subclonal variants are unique to each individual and change in composition, frequency, and sequence context from diagnosis to relapse. Their functional significance is apparent by the observation that many are known variants and cluster within functionally important protein domains. Subclones provide a reservoir of variants that could expand and contribute to the development of drug resistance and relapse. In accord, we accurately identified subclonal variants in AML driver genes NRAS and RUNX1 at allele frequencies between 0.1% and 0.3% at diagnosis, which expanded to comprise a major fraction (14-53%) of the blast population at relapse. Early and accurate detection of subclonal variants with low allele frequency thus offers the opportunity for early intervention, prior to detection of clinical relapse, to improve disease outcome and enhance patient survival.

A high-resolution landscape of mutations in the BCL6 super-enhancer in normal human B cells.

The super-enhancers (SEs) of lineage-specific genes in B cells are off-target sites of somatic hypermutation. However, the inability to detect sufficient numbers of mutations in normal human B cells has precluded the generation of a high-resolution mutational landscape of SEs. Here we captured and sequenced 12 B cell SEs at single-nucleotide resolution from 10 healthy individuals across diverse ethnicities. We detected a total of approximately 9,000 subclonal mutations (allele frequencies <0.1%); of these, approximately 8,000 are present in the BCL6 SE alone. Within the BCL6 SE, we identified 3 regions of clustered mutations in which the mutation frequency is ∼7 × 10-4 Mutational spectra show a predominance of C > T/G > A and A > G/T > C substitutions, consistent with the activities of activation-induced-cytidine deaminase (AID) and the A-T mutator, DNA polymerase η, respectively, in mutagenesis in normal B cells. Analyses of mutational signatures further corroborate the participation of these factors in this process. Single base substitution signatures SBS85, SBS37, and SBS39 were found in the BCL6 SE. While SBS85 is a denoted signature of AID in lymphoid cells, the etiologies of SBS37 and SBS39 are unknown. Our analysis suggests the contribution of error-prone DNA polymerases to the latter signatures. The high-resolution mutation landscape has enabled accurate profiling of subclonal mutations in B cell SEs in normal individuals. By virtue of the fact that subclonal SE mutations are clonally expanded in B cell lymphomas, our studies also offer the potential for early detection of neoplastic alterations.

Detecting ultralow-frequency mutations by Duplex Sequencing.

Duplex Sequencing (DS) is a next-generation sequencing methodology capable of detecting a single mutation among >1 × 10(7) wild-type nucleotides, thereby enabling the study of heterogeneous populations and very-low-frequency genetic alterations. DS can be applied to any double-stranded DNA sample, but it is ideal for small genomic regions of <1 Mb in size. The method relies on the ligation of sequencing adapters harboring random yet complementary double-stranded nucleotide sequences to the sample DNA of interest. Individually labeled strands are then PCR-amplified, creating sequence 'families' that share a common tag sequence derived from the two original complementary strands. Mutations are scored only if the variant is present in the PCR families arising from both of the two DNA strands. Here we provide a detailed protocol for efficient DS adapter synthesis, library preparation and target enrichment, as well as an overview of the data analysis workflow. The protocol typically takes 1-3 d.

Sphingosine, a modulator of human translesion DNA polymerase activity.

Translesion (TLS) DNA polymerases are specialized, error-prone enzymes that synthesize DNA across bulky, replication-stalling DNA adducts. In so doing, they facilitate the progression of DNA synthesis and promote cell proliferation. To potentiate the effect of cancer chemotherapeutic regimens, we sought to identify inhibitors of TLS DNA polymerases. We screened five libraries of ∼ 3000 small molecules, including one comprising ∼ 600 nucleoside analogs, for their effect on primer extension activity of DNA polymerase η (Pol η). We serendipitously identified sphingosine, a lipid-signaling molecule that robustly stimulates the activity of Pol η by ∼ 100-fold at low micromolar concentrations but inhibits it at higher concentrations. This effect is specific to the Y-family DNA polymerases, Pols η, κ, and ι. The addition of a single phosphate group on sphingosine completely abrogates this effect. Likewise, the inclusion of other sphingolipids, including ceramide and sphingomyelin to extension reactions does not elicit this response. Sphingosine increases the rate of correct and incorrect nucleotide incorporation while having no effect on polymerase processivity. Endogenous Pol η activity is modulated similarly as the recombinant enzyme. Importantly, sphingosine-treated cells exhibit increased lesion bypass activity, and sphingosine tethered to membrane lipids mimics the effects of free sphingosine. Our studies have uncovered sphingosine as a modulator of TLS DNA polymerase activity; this property of sphingosine may be associated with its known role as a signaling molecule in regulating cell proliferation in response to cellular stress.

A rapid assay for measuring nucleotide excision repair by oligonucleotide retrieval.

Nucleotide excision repair (NER) excises bulky DNA lesions induced by mutagens and carcinogens. The repair process includes recognition of DNA damage, excision of a short patch of nucleotides containing the damaged base, re-synthesis of a new DNA strand and ligation of the nicks to restore the sequence integrity. Mutation or aberrant transcription of NER genes reduces repair efficiency and results in the accumulation of mutations that is associated with the development of cancer. Here we present a rapid, sensitive and quantitative assay to measure NER activity in human cells, which we term the Oligonucleotide Retrieval Assay (ORA). We used oligonucleotide constructs containing the UV-damaged adduct, cyclobutane pyrimidine dimer (CPD), to transfect human cells, and retrieved the oligonucleotides for quantification of the repaired, CPD-free DNA by real-time quantitative PCR. We demonstrate that ORA can quantify the extent of NER in diverse cell types, including immortalized, primary and stem-like cells.

The Werner syndrome exonuclease facilitates DNA degradation and high fidelity DNA polymerization by human DNA polymerase δ.

DNA Polymerase δ (Pol δ) and the Werner syndrome protein, WRN, are involved in maintaining cellular genomic stability. Pol δ synthesizes the lagging strand during replication of genomic DNA and also functions in the synthesis steps of DNA repair and recombination. WRN is a member of the RecQ helicase family, loss of which results in the premature aging and cancer-prone disorder, Werner syndrome. Both Pol δ and WRN encode 3' → 5' DNA exonuclease activities. Pol δ exonuclease removes 3'-terminal mismatched nucleotides incorporated during replication to ensure high fidelity DNA synthesis. WRN exonuclease degrades DNA containing alternate secondary structures to prevent formation and enable resolution of stalled replication forks. We now observe that similarly to WRN, Pol δ degrades alternate DNA structures including bubbles, four-way junctions, and D-loops. Moreover, WRN and Pol δ form a complex with enhanced ability to hydrolyze these structures. We also present evidence that WRN can proofread for Pol δ; WRN excises 3'-terminal mismatches to enable primer extension by Pol δ. Consistent with our in vitro observations, we show that WRN contributes to the maintenance of DNA synthesis fidelity in vivo. Cells expressing limiting amounts (∼10% of normal) of WRN have elevated mutation frequencies compared with wild-type cells. Together, our data highlight the importance of WRN exonuclease activity and its cooperativity with Pol δ in preserving genome stability, which is compromised by the loss of WRN in Werner syndrome.

A novel ING2 isoform, ING2b, synergizes with ING2a to prevent cell cycle arrest and apoptosis.

We identified a novel inhibitor of growth family member 2 (ING2) isoform, ING2b, which shares exon 2 with ING2a, but lacks the N-terminal p53 binding region. Contrary to ING2a, ING2b's promoter has no p53 binding sites. Consistently, activation of p53 led to suppression of ING2a, leaving ING2b unaffected. Through isoform-specific targeting, we showed that ING2a knockdown suppressed cell growth only in the presence of p53, ING2b knockdown had no effect on cell growth, and knockdown of both induced cell cycle arrest and apoptosis independently of p53. ING2a and ING2b have compensatory roles that protect cells from cell cycle arrest and apoptosis and may be involved in development of chemotherapeutic resistance.

Inhibitor of growth 4 suppresses cell spreading and cell migration by interacting with a novel binding partner, liprin alpha1.

Inhibitor of growth 4 (ING4) is a candidate tumor suppressor that plays a major role in gene regulation, cell cycle control, apoptosis, and angiogenesis. ING4 expression is down-regulated in glioblastoma cells and head and neck squamous cell carcinoma. Here, we identified liprin alpha1/PPFIA1, a cytoplasmic protein necessary for focal adhesion formation and axon guidance, as a novel interacting protein with ING4. ING4 and liprin alpha1 colocalized at lamellipodia in the vicinity of vinculin. Overexpressed ING4 suppressed cell spreading and cell migration. In contrast, overexpressed liprin alpha1 enhanced cell spreading and cell migration. Knockdown of endogenous ING4 with RNA interference induced cell motility, whereas knockdown of endogenous liprin alpha1 suppressed cell motility. ING4 also suppressed cell motility that was enhanced by liprin alpha1. However, ING4 did not further suppress cell motility when liprin alpha1 was suppressed with RNA interference, suggesting a functional and mechanistic interdependence between these proteins. In addition to its nuclear functions, cytoplasmic ING4 interacts with liprin alpha1 to regulate cell migration and, with its known antiangiogenic function, may prevent invasion and metastasis.

Novel splice variants of ING4 and their possible roles in the regulation of cell growth and motility.

The ING4 gene is a candidate tumor suppressor gene that functions in cell proliferation, contact inhibition, and angiogenesis. We identified three novel splice variants of ING4 with differing activities in controlling cell proliferation, cell spreading, and cell migration. ING4_v1 (the longest splice variant), originally identified as ING4, encodes an intact nuclear localization signal (NLS), whereas the other three splice variants (ING4_v2, ING4_v3, and ING4_v4) lack the full NLS, resulting in increased cytoplasmic localization of these proteins. We found that one of the three ING4 variants, ING4_v2, is expressed at the same level as the original ING4 (ING4_v1), suggesting that ING4 variants may have significant biological functions. Growth suppressive effects of the variants that have a partial NLS (ING4_v2 and ING4_v4) were attenuated by a weaker effect of the variants on p21(WAF1) promoter activation. ING4_v4 lost cell spreading and migration suppressive effects; on the other hand, ING4_v2 retained a cell migration suppressive effect but lost a cell spreading suppressive effect. Therefore, ING4_v2, which localized primarily into cytoplasm, might have an important role in the regulation of cell migration. We also found that ING4_v4 played dominant-negative roles in the induction of p21(WAF1) promoter activation and in the suppression of cell motility by ING4_v1. In addition, ING4 variants had different binding affinities to two cytoplasmic proteins, protein-tyrosine phosphatase, receptor type, f polypeptide (PTPRF), interacting protein (liprin), alpha1, and G3BP2a. Understanding the functions of the four splice variants may aid in defining their roles in human carcinogenesis.

ING2 regulates the onset of replicative senescence by induction of p300-dependent p53 acetylation.

ING2 is a candidate tumor suppressor gene that can activate p53 by enhancing its acetylation. Here, we demonstrate that ING2 is also involved in p53-mediated replicative senescence. ING2 protein expression increased in late-passage human primary cells, and it colocalizes with serine 15-phosphorylated p53. ING2 and p53 also complexed with the histone acetyltransferase p300. ING2 enhanced the interaction between p53 and p300 and acted as a cofactor for p300-mediated p53 acetylation. The level of ING2 expression directly modulated the onset of replicative senescence. While overexpression of ING2 induced senescence in young fibroblasts in a p53-dependent manner, expression of ING2 small interfering RNA delayed the onset of senescence. Hence, ING2 can act as a cofactor of p300 for p53 acetylation and thereby plays a positive regulatory role during p53-mediated replicative senescence.

Tumor suppressor p53 represses transcription of RECQ4 helicase.

RECQ4 is a member of the RecQ helicase family, which has been implicated in the regulation of DNA replication, recombination and repair. p53 modulates the functions of RecQ helicases including BLM and WRN. In this study, we demonstrate that p53 can regulate the transcription of RECQ4. Using nontransformed, immortalized normal human fibroblasts, we show that p53-dependent downregulation of RECQ4 expression occurred in G1-arrested cells, both in the absence or presence of exogenous DNA damage. Wild-type p53 (but not the tumor-derived mutant forms) repressed RECQ4 promoter activity. The camptothecin or etoposide-dependent p53-mediated repression was attenuated by trichostatin A (TSA), an inhibitor of histone deacetylases (HDACs). Repression of the RECQ4 promoter was accompanied with an increased accumulation of HDAC1, and the loss of SP1 and p53 binding to the promoter. The simultaneous formation of a camptothecin-dependent p53-SP1 complex indicated its occurrence outside of the RECQ4 promoter. These data suggest that p53-mediated repression of RECQ4 transcription during DNA damage results from the modulation of the promoter occupancy of transcription activators and repressors.

The N-terminal domain of the large subunit of human replication protein A binds to Werner syndrome protein and stimulates helicase activity.

Werner syndrome (WS) is a recessive inherited human disease characterized by the early onset of aging. The gene mutated in WS encodes a DNA helicase that unwinds the double helical structure of DNA in the 3'-->5' direction as well as a 3'-->5' exonuclease. Our previous studies indicated that the activity of Werner syndrome helicase (WRN) could be stimulated by human replication protein A (hRPA), a heterotrimeric single-stranded DNA binding protein. We now localize the interaction between WRN and hRPA by measuring the stimulation of helicase activity and the binding of WRN by hRPA and its derivatives. The large subunit of hRPA (hRPA70) stimulates WRN helicase to the same extent as the hRPA heterotrimer, whereas the dimer of the two smaller subunits (hRPA 32.14) does not stimulate. By examining hRPA70 mutants with progressive deletions from either the C- or N-terminus, we found that the domain responsible for stimulation lies in the N-terminal half of the protein. By using enzyme-linked immunosorbent assay (ELISA) to examine physical interaction between WRN and the same deletion mutants, we found that the WRN-binding motif is located within amino acids 100-300 and overlaps with the single-stranded DNA binding domain (amino acids 150-450). We suggest that hRPA, by engaging in both protein-protein and protein-DNA interactions, facilitates unwinding events catalyzed by WRN helicase during DNA synthetic processes. These data should help further elucidation of the molecular mechanisms of genetic instability and premature aging phenotypes manifested by WS.

Mutations in the alpha8 loop of human APE1 alter binding and cleavage of DNA containing an abasic site.

Recent crystallographic studies reveal loops in human AP endonuclease 1 (APE1) that interact with the major and minor grooves of DNA containing apurinic/apyrimidinic (AP) sites. These loops are postulated to stabilize the DNA helix and the flipped out AP residue. The loop alpha8 interacts with the major groove on the 3' side of the AP site. To determine the essentiality of the amino acids that constitute the alpha8 loop, we created a mutant library containing random nucleotides at codons 222-229 that, in wild-type APE1, specify the sequence NPKGNKKN. Upon expression of the library (2 x 10(6) different clones) in Escherichia coli and multiple rounds of selection with the alkylating agent methyl-methane sulfonate (MMS), we obtained approximately 2 x 10(5) active mutants that complemented the MMS sensitivity of AP endonuclease-deficient E. coli. DNA sequencing showed that active mutants tolerated amino acid substitutions at all eight randomized positions. Basic and uncharged polar amino acids together comprised the majority of substitutions, reflecting the positively charged, polar character of the wild-type loop. Asn-222, Asn-226, and Asn-229 exhibited the least mutability, consistent with x-ray data showing that each asparagine contacts a DNA phosphate. Substitutions at residues 226-229, located nearer to the AP site, that reduced basicity or hydrogen bonding potential, increased Km 2- to 6-fold and decreased AP site binding; substitutions at residues 222-225 exhibited lesser effects. This initial mutational analysis of the alpha8 loop supports and extends the conclusion of crystallographic studies that the loop is important for binding of AP.DNA and AP site incision.

The processing of Holliday junctions by BLM and WRN helicases is regulated by p53.

BLM, WRN, and p53 are involved in the homologous DNA recombination pathway. The DNA structure-specific helicases, BLM and WRN, unwind Holliday junctions (HJ), an activity that could suppress inappropriate homologous recombination during DNA replication. Here, we show that purified, recombinant p53 binds to BLM and WRN helicases and attenuates their ability to unwind synthetic HJ in vitro. The p53 248W mutant reduces abilities of both to bind HJ and inhibit helicase activities, whereas the p53 273H mutant loses these abilities. Moreover, full-length p53 and a C-terminal polypeptide (residues 373-383) inhibit the BLM and WRN helicase activities, but phosphorylation at Ser(376) or Ser(378) completely abolishes this inhibition. Following blockage of DNA replication, Ser(15) phospho-p53, BLM, and RAD51 colocalize in nuclear foci at sites likely to contain DNA replication intermediates in cells. Our results are consistent with a novel mechanism for p53-mediated regulation of DNA recombinational repair that involves p53 post-translational modifications and functional protein-protein interactions with BLM and WRN DNA helicases.