Heterogeneous resistance to quizartinib in acute myeloid leukemia revealed by single-cell analysis.

Genomic studies have revealed significant branching heterogeneity in cancer. Studies of resistance to tyrosine kinase inhibitor therapy have not fully reflected this heterogeneity because resistance in individual patients has been ascribed to largely mutually exclusive on-target or off-target mechanisms in which tumors either retain dependency on the target oncogene or subvert it through a parallel pathway. Using targeted sequencing from single cells and colonies from patient samples, we demonstrate tremendous clonal diversity in the majority of acute myeloid leukemia (AML) patients with activating FLT3 internal tandem duplication mutations at the time of acquired resistance to the FLT3 inhibitor quizartinib. These findings establish that clinical resistance to quizartinib is highly complex and reflects the underlying clonal heterogeneity of AML.

Validation of ITD mutations in FLT3 as a therapeutic target in human acute myeloid leukaemia.

Effective targeted cancer therapeutic development depends upon distinguishing disease-associated 'driver' mutations, which have causative roles in malignancy pathogenesis, from 'passenger' mutations, which are dispensable for cancer initiation and maintenance. Translational studies of clinically active targeted therapeutics can definitively discriminate driver from passenger lesions and provide valuable insights into human cancer biology. Activating internal tandem duplication (ITD) mutations in FLT3 (FLT3-ITD) are detected in approximately 20% of acute myeloid leukaemia (AML) patients and are associated with a poor prognosis. Abundant scientific and clinical evidence, including the lack of convincing clinical activity of early FLT3 inhibitors, suggests that FLT3-ITD probably represents a passenger lesion. Here we report point mutations at three residues within the kinase domain of FLT3-ITD that confer substantial in vitro resistance to AC220 (quizartinib), an active investigational inhibitor of FLT3, KIT, PDGFRA, PDGFRB and RET; evolution of AC220-resistant substitutions at two of these amino acid positions was observed in eight of eight FLT3-ITD-positive AML patients with acquired resistance to AC220. Our findings demonstrate that FLT3-ITD can represent a driver lesion and valid therapeutic target in human AML. AC220-resistant FLT3 kinase domain mutants represent high-value targets for future FLT3 inhibitor development efforts.

The diagnostic application of RNA sequencing in patients with thyroid cancer: an analysis of 851 variants and 133 fusions in 524 genes.

BACKGROUND: Thyroid carcinomas are known to harbor oncogenic driver mutations and advances in sequencing technology now allow the detection of these in fine needle aspiration biopsies (FNA). Recent work by The Cancer Genome Atlas (TCGA) Research Network has expanded the number of genetic alterations detected in papillary thyroid carcinomas (PTC). We sought to investigate the prevalence of these and other genetic alterations in diverse subtypes of thyroid nodules beyond PTC, including a variety of samples with benign histopathology. This is the first clinical evaluation of a large panel of TCGA-reported genomic alterations in thyroid FNAs. RESULTS: In FNAs, genetic alterations were detected in 19/44 malignant samples (43% sensitivity) and in 7/44 histopathology benign samples (84% specificity). Overall, after adding a cohort of tissue samples, 38/76 (50%) of histopathology malignant samples were found to harbor a genetic alteration, while 15/75 (20%) of benign samples were also mutated. The most frequently mutated malignant subtypes were medullary thyroid carcinoma (9/12, 75%) and PTC (14/30, 47%). Additionally, follicular adenoma, a benign subtype of thyroid neoplasm, was also found to harbor mutations (12/29, 41%). Frequently mutated genes in malignant samples included BRAF (20/76, 26%) and RAS (9/76, 12%). Of the TSHR variants detected, (6/7, 86%) were in benign nodules. In a direct comparison of the same FNA also tested by an RNA-based gene expression classifier (GEC), the sensitivity of genetic alterations alone was 42%, compared to the 91% sensitivity achieved by the GEC. The specificity based only on genetic alterations was 84%, compared to 77% specificity with the GEC. CONCLUSIONS: While the genomic landscape of all thyroid neoplasm subtypes will inevitably be elucidated, caution should be used in the early adoption of published mutations as the sole predictor of malignancy in thyroid. The largest set of such mutations known to date detects only a portion of thyroid carcinomas in preoperative FNAs in our cohort and thus is not sufficient to rule out cancer. Due to the finding that variants are also found in benign nodules, testing only GEC suspicious nodules may be helpful in avoiding false positives and altering the extent of treatment when selected mutations are found.

Direct detection of DNA methylation during single-molecule, real-time sequencing.

We describe the direct detection of DNA methylation, without bisulfite conversion, through single-molecule, real-time (SMRT) sequencing. In SMRT sequencing, DNA polymerases catalyze the incorporation of fluorescently labeled nucleotides into complementary nucleic acid strands. The arrival times and durations of the resulting fluorescence pulses yield information about polymerase kinetics and allow direct detection of modified nucleotides in the DNA template, including N6-methyladenine, 5-methylcytosine and 5-hydroxymethylcytosine. Measurement of polymerase kinetics is an intrinsic part of SMRT sequencing and does not adversely affect determination of primary DNA sequence. The various modifications affect polymerase kinetics differently, allowing discrimination between them. We used these kinetic signatures to identify adenine methylation in genomic samples and found that, in combination with circular consensus sequencing, they can enable single-molecule identification of epigenetic modifications with base-pair resolution. This method is amenable to long read lengths and will likely enable mapping of methylation patterns in even highly repetitive genomic regions.

A flexible and efficient template format for circular consensus sequencing and SNP detection.

A novel template design for single-molecule sequencing is introduced, a structure we refer to as a SMRTbell template. This structure consists of a double-stranded portion, containing the insert of interest, and a single-stranded hairpin loop on either end, which provides a site for primer binding. Structurally, this format resembles a linear double-stranded molecule, and yet it is topologically circular. When placed into a single-molecule sequencing reaction, the SMRTbell template format enables a consensus sequence to be obtained from multiple passes on a single molecule. Furthermore, this consensus sequence is obtained from both the sense and antisense strands of the insert region. In this article, we present a universal method for constructing these templates, as well as an application of their use. We demonstrate the generation of high-quality consensus accuracy from single molecules, as well as the use of SMRTbell templates in the identification of rare sequence variants.

Quantitative and comprehensive decomposition of the ion atmosphere around nucleic acids.

The ion atmosphere around nucleic acids critically affects biological and physical processes such as chromosome packing, RNA folding, and molecular recognition. However, the dynamic nature of the ion atmosphere renders it difficult to characterize. The basic thermodynamic description of this atmosphere, a full accounting of the type and number of associated ions, has remained elusive. Here we provide the first complete accounting of the ion atmosphere, using buffer equilibration and atomic emission spectroscopy (BE-AES) to accurately quantitate the cation association and anion depletion. We have examined the influence of ion size and charge on ion occupancy around simple, well-defined DNA molecules. The relative affinity of monovalent and divalent cations correlates inversely with their size. Divalent cations associate preferentially over monovalent cations; e.g., with Na+ in 4-fold excess of Mg2+ (20 vs 5 mM), the ion atmosphere nevertheless has 3-fold more Mg2+ than Na+. Further, the dicationic polyamine putrescine2+ does not compete effectively for association relative to divalent metal ions, presumably because of its lower charge density. These and other BE-AES results can be used to evaluate and guide the improvement of electrostatic treatments. As a first step, we compare the BE-AES results to predictions from the widely used nonlinear Poisson Boltzmann (NLPB) theory and assess the applicability and precision of this theory. In the future, BE-AES in conjunction with improved theoretical models, can be applied to complex binding and folding equilibria of nucleic acids and their complexes, to parse the electrostatic contribution from the overall thermodynamics of important biological processes.

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding.

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.

Low specificity of metal ion binding in the metal ion core of a folded RNA.

The structure and activity of nucleic acids depend on their interactions with metal ions. Fundamental to these interactions is the degree of specificity observed between the metal ions and nucleic acids, and a complete description of nucleic acid folding requires that we understand the nature of the interactions with metal ions, including specificity. The prior demonstration that high concentrations of monovalent cations prevent nonspecific association of divalent ions with nucleic acids provides a novel and powerful means to examine site-specific metal ion binding isolated from complicating effects of the ion atmosphere. Using these high monovalent cation solution conditions we have monitored the affinity of a series of divalent metal ions for two site-specific metal ion binding sites in the P4-P6 domain of the Tetrahymena group I intron ribozyme. The metal ion core of this highly structured RNA binds two divalent metal ions under these conditions. Despite multiple metal ion-RNA interactions observed in the X-ray crystallographic structure of P4-P6 RNA at the metal ion binding sites, these sites exhibit low specificity among Mn(2+), Mg(2+), Ca(2+), Ni(2+), and Zn(2+). Nevertheless, the largest divalent metal ions tested, Sr(2+) and Ba(2+), were excluded from binding, exhibiting affinities at least two orders of magnitude weaker than observed for the other metal ions. Thus, a picture emerges of two metal ion binding sites, each with a high tolerance for metal ions with different properties but also with limits to accommodation. These limits presumably arise from steric or electrostatic features of the metal ion binding sites.

Measuring the folding transition time of single RNA molecules.

We describe a new, time-apertured photon correlation method for resolving the transition time between two states of RNA in folding--i.e., the time of the transition between states rather than the time spent in each state. Single molecule fluorescence resonance energy transfer and fluorescence correlation spectroscopy are used to obtain these measurements. Individual RNA molecules are labeled with fluorophores such as Cy3 and Cy5. Those molecules are then immobilized on a surface and observed for many seconds during which time the molecules spontaneously switch between two conformational states with different levels of flourescence resonance energy transfer efficiency. Single photons are counted from each fluorophore and cross correlated in a small window around a transition. The average of over 1000 cross correlations can be fit to a polynomial, which can determine transition times as short as the average photon emission interval. We applied the method to the P4-P6 domain of the Tetrahymena group I self-splicing intron to yield the folding transition time of 240 micros. The unfolding time is found to be too short to measure with this method.

Determining the Mg2+ stoichiometry for folding an RNA metal ion core.

The folding and catalytic function of RNA molecules depend on their interactions with divalent metal ions, such as magnesium. As with every molecular process, the most basic knowledge required for understanding the close relationship of an RNA with its metal ions is the stoichiometry of the interaction. Unfortunately, inventories of the numbers of divalent ions associated with unfolded and folded RNA states have been unattainable. A common approach has been to interpret Hill coefficients fit to folding equilibria as the number of metal ions bound upon folding. However, this approach is vitiated by the presence of diffusely associated divalent ions in a dynamic ion atmosphere and by the likelihood of multiple transitions along a folding pathway. We demonstrate that the use of molar concentrations of background monovalent salt can alleviate these complications. These simplifying solution conditions allow a precise determination of the stoichiometry of the magnesium ions involved in folding the metal ion core of the P4-P6 domain of the Tetrahymena group I ribozyme. Hill analysis of hydroxyl radical footprinting data suggests that the P4-P6 RNA core folds cooperatively upon the association of two metal ions. This unexpectedly small stoichiometry is strongly supported by counting magnesium ions associated with the P4-P6 RNA via fluorescence titration and atomic emission spectroscopy. By pinpointing the metal ion stoichiometry, these measurements provide a critical but previously missing step in the thermodynamic dissection of the coupling between metal ion binding and RNA folding.