Multiplex accurate sensitive quantitation (MASQ) with application to minimal residual disease in acute myeloid leukemia.

Measuring minimal residual disease in cancer has applications for prognosis, monitoring treatment and detection of recurrence. Simple sequence-based methods to detect nucleotide substitution variants have error rates (about 10-3) that limit sensitive detection. We developed and characterized the performance of MASQ (multiplex accurate sensitive quantitation), a method with an error rate below 10-6. MASQ counts variant templates accurately in the presence of millions of host genomes by using tags to identify each template and demanding consensus over multiple reads. Since the MASQ protocol multiplexes 50 target loci, we can both integrate signal from multiple variants and capture subclonal response to treatment. Compared to existing methods for variant detection, MASQ achieves an excellent combination of sensitivity, specificity and yield. We tested MASQ in a pilot study in acute myeloid leukemia (AML) patients who entered complete remission. We detect leukemic variants in the blood and bone marrow samples of all five patients, after induction therapy, at levels ranging from 10-2 to nearly 10-6. We observe evidence of sub-clonal structure and find higher target variant frequencies in patients who go on to relapse, demonstrating the potential for MASQ to quantify residual disease in AML.

Utility of Single-Cell Genomics in Diagnostic Evaluation of Prostate Cancer.

A distinction between indolent and aggressive disease is a major challenge in diagnostics of prostate cancer. As genetic heterogeneity and complexity may influence clinical outcome, we have initiated studies on single tumor cell genomics. In this study, we demonstrate that sparse DNA sequencing of single-cell nuclei from prostate core biopsies is a rich source of quantitative parameters for evaluating neoplastic growth and aggressiveness. These include the presence of clonal populations, the phylogenetic structure of those populations, the degree of the complexity of copy-number changes in those populations, and measures of the proportion of cells with clonal copy-number signatures. The parameters all showed good correlation to the measure of prostatic malignancy, the Gleason score, derived from individual prostate biopsy tissue cores. Remarkably, a more accurate histopathologic measure of malignancy, the surgical Gleason score, agrees better with these genomic parameters of diagnostic biopsy than it does with the diagnostic Gleason score and related measures of diagnostic histopathology. This is highly relevant because primary treatment decisions are dependent upon the biopsy and not the surgical specimen. Thus, single-cell analysis has the potential to augment traditional core histopathology, improving both the objectivity and accuracy of risk assessment and inform treatment decisions.Significance: Genomic analysis of multiple individual cells harvested from prostate biopsies provides an indepth view of cell populations comprising a prostate neoplasm, yielding novel genomic measures with the potential to improve the accuracy of diagnosis and prognosis in prostate cancer. Cancer Res; 78(2); 348-58. ©2017 AACR.

Whole-genome single-cell copy number profiling from formalin-fixed paraffin-embedded samples.

A substantial proportion of tumors consist of genotypically distinct subpopulations of cancer cells. This intratumor genetic heterogeneity poses a substantial challenge for the implementation of precision medicine. Single-cell genomics constitutes a powerful approach to resolve complex mixtures of cancer cells by tracing cell lineages and discovering cryptic genetic variations that would otherwise be obscured in tumor bulk analyses. Because of the chemical alterations that result from formalin fixation, single-cell genomic approaches have largely remained limited to fresh or rapidly frozen specimens. Here we describe the development and validation of a robust and accurate methodology to perform whole-genome copy-number profiling of single nuclei obtained from formalin-fixed paraffin-embedded clinical tumor samples. We applied the single-cell sequencing approach described here to study the progression from in situ to invasive breast cancer, which revealed that ductal carcinomas in situ show intratumor genetic heterogeneity at diagnosis and that these lesions may progress to invasive breast cancer through a variety of evolutionary processes.

Single-cell copy number analysis of prostate cancer cells captured with geometrically enhanced differential immunocapture microdevices.

Limited access to tumor tissue makes repeated sampling and real-time tracking of cancer progression infeasible. Circulating tumor cells (CTCs) provide the capacity for real-time genetic characterization of a disseminating tumor cell population via a simple blood draw. However, there is no straightforward method to analyze broadscale genetic rearrangements in this heterogeneous cell population at the single cell level. We present a one-step controllable chemical extraction of whole nuclei from prostate cancer cells captured using geometrically enhanced differential immunocapture (GEDI) microdevices. We have successfully used copy number profile analysis to differentiate between two unique cancer cell line populations of metastatic origin (LNCaP and VCaP) and to analyze key mutations important in disease progression.

Rapid phenotypic and genomic change in response to therapeutic pressure in prostate cancer inferred by high content analysis of single circulating tumor cells.

Timely characterization of a cancer's evolution is required to predict treatment efficacy and to detect resistance early. High content analysis of single Circulating Tumor Cells (CTCs) enables sequential characterization of genotypic, morphometric and protein expression alterations in real time over the course of cancer treatment. This concept was investigated in a patient with castrate-resistant prostate cancer progressing through both chemotherapy and targeted therapy. In this case study, we integrate across four timepoints 41 genome-wide copy number variation (CNV) profiles plus morphometric parameters and androgen receptor (AR) protein levels. Remarkably, little change was observed in response to standard chemotherapy, evidenced by the fact that a unique clone (A), exhibiting highly rearranged CNV profiles and AR+ phenotype was found circulating before and after treatment. However, clinical response and subsequent progression after targeted therapy was associated with the drastic depletion of clone A, followed by the sequential emergence of two distinct CTC sub-populations that differed in both AR genotype and expression phenotype. While AR- cells with flat or pseudo-diploid CNV profiles (clone B) were identified at the time of response, a new tumor lineage of AR+ cells (clone C) with CNV altered profiles was detected during relapse. We showed that clone C, despite phylogenetically related to clone A, possessed a unique set of somatic CNV alterations, including MYC amplification, an event linked to hormone escape. Interesting, we showed that both clones acquired AR gene amplification by deploying different evolutionary paths. Overall, these data demonstrate the timeframe of tumor evolution in response to therapy and provide a framework for the multi-scale analysis of fluid biopsies to quantify and monitor disease evolution in individual patients.

De novo gene disruptions in children on the autistic spectrum.

Exome sequencing of 343 families, each with a single child on the autism spectrum and at least one unaffected sibling, reveal de novo small indels and point substitutions, which come mostly from the paternal line in an age-dependent manner. We do not see significantly greater numbers of de novo missense mutations in affected versus unaffected children, but gene-disrupting mutations (nonsense, splice site, and frame shifts) are twice as frequent, 59 to 28. Based on this differential and the number of recurrent and total targets of gene disruption found in our and similar studies, we estimate between 350 and 400 autism susceptibility genes. Many of the disrupted genes in these studies are associated with the fragile X protein, FMRP, reinforcing links between autism and synaptic plasticity. We find FMRP-associated genes are under greater purifying selection than the remainder of genes and suggest they are especially dosage-sensitive targets of cognitive disorders.

Genome-wide copy number analysis of single cells.

Copy number variation (CNV) is increasingly recognized as an important contributor to phenotypic variation in health and disease. Most methods for determining CNV rely on admixtures of cells in which information regarding genetic heterogeneity is lost. Here we present a protocol that allows for the genome-wide copy number analysis of single nuclei isolated from mixed populations of cells. Single-nucleus sequencing (SNS), combines flow sorting of single nuclei on the basis of DNA content and whole-genome amplification (WGA); this is followed by next-generation sequencing to quantize genomic intervals in a genome-wide manner. Multiplexing of single cells is discussed. In addition, we outline informatic approaches that correct for biases inherent in the WGA procedure and allow for accurate determination of copy number profiles. All together, the protocol takes ∼3 d from flow cytometry to sequence-ready DNA libraries.

Tumour evolution inferred by single-cell sequencing.

Genomic analysis provides insights into the role of copy number variation in disease, but most methods are not designed to resolve mixed populations of cells. In tumours, where genetic heterogeneity is common, very important information may be lost that would be useful for reconstructing evolutionary history. Here we show that with flow-sorted nuclei, whole genome amplification and next generation sequencing we can accurately quantify genomic copy number within an individual nucleus. We apply single-nucleus sequencing to investigate tumour population structure and evolution in two human breast cancer cases. Analysis of 100 single cells from a polygenomic tumour revealed three distinct clonal subpopulations that probably represent sequential clonal expansions. Additional analysis of 100 single cells from a monogenomic primary tumour and its liver metastasis indicated that a single clonal expansion formed the primary tumour and seeded the metastasis. In both primary tumours, we also identified an unexpectedly abundant subpopulation of genetically diverse 'pseudodiploid' cells that do not travel to the metastatic site. In contrast to gradual models of tumour progression, our data indicate that tumours grow by punctuated clonal expansions with few persistent intermediates.

A RSC/nucleosome complex determines chromatin architecture and facilitates activator binding.

How is chromatin architecture established and what role does it play in transcription? We show that the yeast regulatory locus UASg bears, in addition to binding sites for the activator Gal4, sites bound by the RSC complex. RSC positions a nucleosome, evidently partially unwound, in a structure that facilitates Gal4 binding to its sites. The complex comprises a barrier that imposes characteristic features of chromatin architecture. In the absence of RSC, ordinary nucleosomes encroach over the UASg and compete with Gal4 for binding. Taken with our previous work, the results show that both prior to and following induction, specific DNA-binding proteins are the predominant determinants of chromatin architecture at the GAL1/10 genes. RSC/nucleosome complexes are also found scattered around the yeast genome. Higher eukaryotic RSC lacks the specific DNA-binding determinants found on yeast RSC, and evidently Gal4 works in those organisms despite whatever obstacle broadly positioned nucleosomes present.

Genetic dissection of histidine biosynthesis in Arabidopsis.

The biosynthesis of histidine (His) in microorganisms, long studied through the isolation and characterization of auxotrophic mutants, has emerged as a paradigm for the regulation of metabolism and gene expression. Much less is known about His biosynthesis in flowering plants. One limiting factor has been the absence of large collections of informative auxotrophs. We describe here the results of a systematic screen for His auxotrophs of Arabidopsis (Arabidopsis thaliana). Ten insertion mutants disrupted in four different biosynthetic genes (HISN2, HISN3, HISN4, HISN6A) were identified through a combination of forward and reverse genetics and were shown to exhibit an embryo-defective phenotype that could be rescued by watering heterozygous plants with His. Male transmission of the mutant allele was in several cases reduced. Knockouts of two redundant genes (HISN1B and HISN5A) had no visible phenotype. Another mutant blocked in the final step of His biosynthesis (hisn8) and a double mutant altered in the redundant first step of the pathway (hisn1a hisn1b) exhibited a combination of gametophytic and embryonic lethality in heterozygotes. Homozygous mutant seedlings and callus tissue produced from rescued seeds appeared normal when grown in the presence of His but typically senesced after continued growth in the absence of His. These knockout mutants document the importance of His biosynthesis for plant growth and development, provide valuable insights into amino acid transport and source-sink relationships during seed development, and represent a significant addition to the limited collection of well-characterized auxotrophs in flowering plants.

Regulation of lysine catabolism in Arabidopsis through concertedly regulated synthesis of the two distinct gene products of the composite AtLKR/SDH locus.

Lysine catabolism in plants is initiated by a bifunctional LKR/SDH (lysine-ketoglutarate reductase/saccharopine dehydrogenase) enzyme encoded by a single LKR/SDH gene. Yet, the AtLKR/SDH gene of Arabidopsis also encodes a second gene product, namely a monofunctional SDH. To elucidate the regulation of lysine catabolism in Arabidopsis through these two gene products of the AtLKR/SDH gene, an analysis was carried out on the effects of the hormones, abscisic acid and jasmonate, as well as various metabolic and stress signals, including lysine itself, on their mRNA and protein levels. The response of the two gene products to the various treatments was only partially co-ordinated, but the levels of the monofunctional SDH mRNA and protein were always in excess over their bifunctional LKR/SDH counterparts. These results suggest that lysine catabolism is regulated primarily by the first enzyme LKR, while the excess level of SDH enables efficient flux of lysine catabolism following the LKR step. Analysis of transgenic plants expressing beta-glucoronidase fusion constructs with the AtLKR/SDH and monofunctional AtSDH promoters demonstrated that transcriptional regulation contributes to the modulation of expression of the bifunctional LKR/SDH and monofunctional SDH gene products in response to hormonal and metabolic signals. To test whether the enhanced expression of the LKR/SDH gene under various hormonal and metabolic signals is correlated with enhanced lysine catabolism, wild-type Arabidopsis and a knockout mutant lacking lysine catabolism were exposed to abscisic acid and sugar starvation. Free lysine accumulated to significantly higher levels in this knockout mutant than in the wild-type plants.

Synthesis of the Arabidopsis bifunctional lysine-ketoglutarate reductase/saccharopine dehydrogenase enzyme of lysine catabolism is concertedly regulated by metabolic and stress-associated signals.

In plants, excess cellular lysine (Lys) is catabolized into glutamic acid and acetyl-coenzyme A; yet, it is still not clear whether this pathway has other functions in addition to balancing Lys levels. To address this issue, we examined the effects of stress-related hormones, abscisic acid (ABA), and jasmonate, as well as various metabolic signals on the production of the mRNA and polypeptide of the bifunctional Lys-ketoglutarate reductase (LKR)/saccharopine dehydrogenase (SDH) enzyme, which contains the first two linked enzymes of Lys catabolism. The level of LKR/SDH was strongly enhanced by ABA, jasmonate, and sugar starvation, whereas excess sugars and nitrogen starvation reduced its level; thus this pathway appears to fulfill multiple functions in stress-related and carbon/nitrogen metabolism. Treatments with combination of hormones and/or metabolites, as well as use of ABA mutants in conjunction with the tester sugars mannose and 3-O-methyl-glucose further supported the idea that the hormonal and metabolic signals apparently operate through different signal transduction cascades. The stimulation of LKR/SDH protein expression by ABA is regulated by a signal transduction cascade that contains the ABI1-1 and ABI2-1 protein phosphatases. By contrast, the stimulation of LKR/SDH protein expression by sugar starvation is regulated by the hexokinase-signaling cascade in a similar manner to the repression of many photosynthetic genes by sugars. These findings suggest a metabolic and mechanistic link between Lys catabolism and photosynthesis-related metabolism in the regulation of carbon/nitrogen partitioning.