Functional diversity and cooperativity between subclonal populations of pediatric glioblastoma and diffuse intrinsic pontine glioma cells

The failure to develop effective therapies for pediatric glioblastoma (pGBM) and diffuse intrinsic pontine glioma (DIPG) is in part due to their intrinsic heterogeneity. We aimed to quantitatively assess the extent to which this was present in these tumors through subclonal genomic analyses and to determine whether distinct tumor subpopulations may interact to promote tumorigenesis by generating subclonal patient-derived models in vitro and in vivo. Analysis of 142 sequenced tumors revealed multiple tumor subclones, spatially and temporally coexisting in a stable manner as observed by multiple sampling strategies. We isolated genotypically and phenotypically distinct subpopulations that we propose cooperate to enhance tumorigenicity and resistance to therapy. Inactivating mutations in the H4K20 histone methyltransferase KMT5B (SUV420H1), present in <1% of cells, abrogate DNA repair and confer increased invasion and migration on neighboring cells, in vitro and in vivo, through chemokine signaling and modulation of integrins. These data indicate that even rare tumor subpopulations may exert profound effects on tumorigenesis as a whole and may represent a new avenue for therapeutic development. Unraveling the mechanisms of subclonal diversity and communication in pGBM and DIPG will be an important step toward overcoming barriers to effective treatments.

Control Materials and Digital PCR Methods for Evaluation of Circulating Cell-Free DNA Extractions from Plasma

Cell-free DNA is an accessible source of genetic material found naturally in plasma that could be used in many diagnostic applications. Translation of cfDNA analysis methods from research laboratories into the clinic would benefit from controls for monitoring the efficiency of patient sample purification and for quality control of the whole workflow from extraction through to analysis. Here we describe two types of control materials that can be "spiked" into plasma samples to monitor and evaluate different aspects of the workflow. The first control material is an internal control that enables evaluation of extraction efficiency, fragment size bias, and sample inhibition. The second control material serves as a parallel quality control material for measurement of specific genetic targets such as tumor mutations.

Comprehensive molecular characterisation of epilepsy-associated glioneuronal tumours

Glioneuronal tumours are an important cause of treatment-resistant epilepsy. Subtypes of tumour are often poorly discriminated by histological features and may be difficult to diagnose due to a lack of robust diagnostic tools. This is illustrated by marked variability in the reported frequencies across different epilepsy surgical series. To address this, we used DNA methylation arrays and RNA sequencing to assay the methylation and expression profiles within a large cohort of glioneuronal tumours. By adopting a class discovery approach, we were able to identify two distinct groups of glioneuronal tumour, which only partially corresponded to the existing histological classification. Furthermore, by additional molecular analyses, we were able to identify pathogenic mutations in BRAF and FGFR1, specific to each group, in a high proportion of cases. Finally, by interrogating our expression data, we were able to show that each molecular group possessed expression phenotypes suggesting different cellular differentiation: astrocytic in one group and oligodendroglial in the second. Informed by this, we were able to identify CCND1, CSPG4, and PDGFRA as immunohistochemical targets which could distinguish between molecular groups. Our data suggest that the current histological classification of glioneuronal tumours does not adequately represent their underlying biology. Instead, we show that there are two molecular groups within glioneuronal tumours. The first of these displays astrocytic differentiation and is driven by BRAF mutations, while the second displays oligodendroglial differentiation and is driven by FGFR1 mutations.

Digital PCR analysis of circulating tumor DNA: a biomarker for chondrosarcoma diagnosis, prognostication, and residual disease detection

Conventional chondrosarcoma is the most common primary bone tumor in adults. Prognosis corresponds with tumor grade but remains variable, especially for individuals with grade (G) II disease. There are currently no biomarkers available for monitoring or prognostication of chondrosarcoma. Circulating tumor DNA (ctDNA) has recently emerged as a promising biomarker for a broad range of tumor types. To date, little has been done to study the presence of ctDNA and its potential utility in the management of sarcomas, including chondrosarcoma. In this study, we have assessed ctDNA levels in a cohort of 71 patients, 32 with sarcoma, including 29 individuals with central chondrosarcoma (CS) and 39 with locally aggressive and benign bone and soft tissue tumors, using digital PCR. In patients with CS, ctDNA was detected in pretreatment samples in 14/29 patients, which showed clear correlation with tumor grade as demonstrated by the detection of ctDNA in all patients with GIII and dedifferentiated disease (n = 6) and in 8/17 patients with GII disease, but never associated with GI CS. Notably detection of ctDNA preoperatively in GII disease was associated with a poor outcome. A total of 14 patients with CS had ctDNA levels assessed at multiple time points and in most patients there was a clear reduction following surgical removal. This research lays the foundation for larger studies to assess the utility of ctDNA for chondrosarcoma diagnosis, prognostication, early detection of residual disease and monitoring disease progression.

MAPK pathway control of stem cell proliferation and differentiation in the embryonic pituitary provides insights into the pathogenesis of papillary craniopharyngioma

Despite the importance of the RAS-RAF-MAPK pathway in normal physiology and disease of numerous organs, its role during pituitary development and tumourigenesis remains largely unknown. Here, we show that the over-activation of the MAPK pathway, through conditional expression of the gain-of-function alleles BrafV600E and KrasG12D in the developing mouse pituitary, results in severe hyperplasia and abnormal morphogenesis of the gland by the end of gestation. Cell-lineage commitment and terminal differentiation are disrupted, leading to a significant reduction in numbers of most of the hormone-producing cells before birth, with the exception of corticotrophs. Of note, Sox2+ stem cells and clonogenic potential are drastically increased in the mutant pituitaries. Finally, we reveal that papillary craniopharyngioma (PCP), a benign human pituitary tumour harbouring BRAF p.V600E also contains Sox2+ cells with sustained proliferative capacity and disrupted pituitary differentiation. Together, our data demonstrate a crucial function of the MAPK pathway in controlling the balance between proliferation and differentiation of Sox2+ cells and suggest that persistent proliferative capacity of Sox2+ cells may underlie the pathogenesis of PCP.

The use of digital PCR to improve the application of quantitative molecular diagnostic methods for tuberculosis

BACKGROUND

Real-time PCR (qPCR) based methods, such as the Xpert MTB/RIF, are increasingly being used to diagnose tuberculosis (TB). While qualitative methods are adequate for diagnosis, the therapeutic monitoring of TB patients requires quantitative methods currently performed using smear microscopy. The potential use of quantitative molecular measurements for therapeutic monitoring has been investigated but findings have been variable and inconclusive. The lack of an adequate reference method and reference materials is a barrier to understanding the source of such disagreement. Digital PCR (dPCR) offers the potential for an accurate method for quantification of specific DNA sequences in reference materials which can be used to evaluate quantitative molecular methods for TB treatment monitoring.

METHODS

To assess a novel approach for the development of quality assurance materials we used dPCR to quantify specific DNA sequences in a range of prototype reference materials and evaluated accuracy between different laboratories and instruments. The materials were then also used to evaluate the quantitative performance of qPCR and Xpert MTB/RIF in eight clinical testing laboratories.

RESULTS

dPCR was found to provide results in good agreement with the other methods tested and to be highly reproducible between laboratories without calibration even when using different instruments. When the reference materials were analysed with qPCR and Xpert MTB/RIF by clinical laboratories, all laboratories were able to correctly rank the reference materials according to concentration, however there was a marked difference in the measured magnitude.

CONCLUSIONS

TB is a disease where the quantification of the pathogen could lead to better patient management and qPCR methods offer the potential to rapidly perform such analysis. However, our findings suggest that when precisely characterised materials are used to evaluate qPCR methods, the measurement result variation is too high to determine whether molecular quantification of Mycobacterium tuberculosis would provide a clinically useful readout. The methods described in this study provide a means by which the technical performance of quantitative molecular methods can be evaluated independently of clinical variability to improve accuracy of measurement results. These will assist in ultimately increasing the likelihood that such approaches could be used to improve patient management of TB.

GNAS mutations are not detected in parosteal and low-grade central osteosarcomas

Parosteal osteosarcoma, low-grade central osteosarcoma, and fibrous dysplasia share similar histological features that may pose a diagnostic challenge. The detection of GNAS mutations in primary bone tumors has been useful in clinical practice for diagnosing fibrous dysplasia. However, the recent report of GNAS mutations being detected in a significant proportion of parosteal osteosarcoma challenges the specificity of this mutation. As the number of cases reported in this study was small we set out to determine if these results could be reproduced. We studied 97 formalin-fixed paraffin-embedded low-grade osteosarcomas from 90 patients including 62 parosteal osteosarcomas, of which MDM2 amplification was detected in 79%, 11 periosteal osteosarcomas and 24 low-grade central osteosarcoma samples. The mutational status of GNAS was analyzed in codons p.R201, p.Q227, and other less common GNAS alterations by bidirectional Sanger sequencing and/or next generation sequencing using the Life Technologies Ion Torrent platform. GNAS mutations were not detected in any of the low-grade osteosarcomas from which informative DNA was extracted. Our findings therefore support prior observations that GNAS mutations are highly specific for fibrous dysplasia and occur rarely, if ever, in parosteal and other low-grade osteosarcomas.

Glioblastoma adaptation traced through decline of an IDH1 clonal driver and macro-evolution of a double-minute chromosome

BACKGROUND

Glioblastoma (GBM) is the most common malignant brain cancer occurring in adults, and is associated with dismal outcome and few therapeutic options. GBM has been shown to predominantly disrupt three core pathways through somatic aberrations, rendering it ideal for precision medicine approaches.

METHODS

We describe a 35-year-old female patient with recurrent GBM following surgical removal of the primary tumour, adjuvant treatment with temozolomide and a 3-year disease-free period. Rapid whole-genome sequencing (WGS) of three separate tumour regions at recurrence was carried out and interpreted relative to WGS of two regions of the primary tumour.

RESULTS

We found extensive mutational and copy-number heterogeneity within the primary tumour. We identified a TP53 mutation and two focal amplifications involving PDGFRA, KIT and CDK4, on chromosomes 4 and 12. A clonal IDH1 R132H mutation in the primary, a known GBM driver event, was detectable at only very low frequency in the recurrent tumour. After sub-clonal diversification, evidence was found for a whole-genome doubling event and a translocation between the amplified regions of PDGFRA, KIT and CDK4, encoded within a double-minute chromosome also incorporating miR26a-2. The WGS analysis uncovered progressive evolution of the double-minute chromosome converging on the KIT/PDGFRA/PI3K/mTOR axis, superseding the IDH1 mutation in dominance in a mutually exclusive manner at recurrence, consequently the patient was treated with imatinib. Despite rapid sequencing and cancer genome-guided therapy against amplified oncogenes, the disease progressed, and the patient died shortly after.

CONCLUSION

This case sheds light on the dynamic evolution of a GBM tumour, defining the origins of the lethal sub-clone, the macro-evolutionary genomic events dominating the disease at recurrence and the loss of a clonal driver. Even in the era of rapid WGS analysis, cases such as this illustrate the significant hurdles for precision medicine success.

Towards standardisation of cell-free DNA measurement in plasma: controls for extraction efficiency, fragment size bias and quantification

Circulating cell-free DNA (cfDNA) is becoming an important clinical analyte for prenatal testing, cancer diagnosis and cancer monitoring. The extraction stage is critical in ensuring clinical sensitivity of analytical methods measuring minority nucleic acid fractions, such as foetal-derived sequences in predominantly maternal cfDNA. Consequently, quality controls are required for measurement of extraction efficiency, fragment size bias and yield for validation of cfDNA methods. We evaluated the utility of an external DNA spike for monitoring these parameters in a study comparing three specific cfDNA extraction methods [QIAamp® circulating nucleic acid (CNA) kit, NucleoSpin® Plasma XS (NS) kit and FitAmp™ plasma/serum DNA isolation (FA) kit] with the commonly used QIAamp DNA blood mini (DBM) kit. We found that the extraction efficiencies of the kits ranked in the order CNA kit > DBM kit > NS kit > FA kit, and the CNA and NS kits gave a better representation of smaller DNA fragments in the extract than the DBM kit. We investigated means of improved reporting of cfDNA yield by comparing quantitative PCR measurements of seven different reference gene assays in plasma samples and validating these with digital PCR. We noted that the cfDNA quantities based on measurement of some target genes (e.g. TERT) were, on average, more than twofold higher than those of other assays (e.g. ERV3). We conclude that analysis and averaging of multiple reference genes using a GeNorm approach gives a more reliable estimate of total cfDNA quantity.

Figure

Lineage tracing reveals multipotent stem cells maintain human adenomas and the pattern of clonal expansion in tumor evolution

The genetic and morphological development of colorectal cancer is a paradigm for tumorigenesis. However, the dynamics of clonal evolution underpinning carcinogenesis remain poorly understood. Here we identify multipotential stem cells within human colorectal adenomas and use methylation patterns of nonexpressed genes to characterize clonal evolution. Numerous individual crypts from six colonic adenomas and a hyperplastic polyp were microdissected and characterized for genetic lesions. Clones deficient in cytochrome c oxidase (CCO(-)) were identified by histochemical staining followed by mtDNA sequencing. Topographical maps of clone locations were constructed using a combination of these data. Multilineage differentiation within clones was demonstrated by immunofluorescence. Methylation patterns of adenomatous crypts were determined by clonal bisulphite sequencing; methylation pattern diversity was compared with a mathematical model to infer to clonal dynamics. Individual adenomatous crypts were clonal for mtDNA mutations and contained both mucin-secreting and neuroendocrine cells, demonstrating that the crypt contained a multipotent stem cell. The intracrypt methylation pattern was consistent with the crypts containing multiple competing stem cells. Adenomas were epigenetically diverse populations, suggesting that they were relatively mitotically old populations. Intratumor clones typically showed less diversity in methylation pattern than the tumor as a whole. Mathematical modeling suggested that recent clonal sweeps encompassing the whole adenoma had not occurred. Adenomatous crypts within human tumors contain actively dividing stem cells. Adenomas appeared to be relatively mitotically old populations, pocketed with occasional newly generated subclones that were the result of recent rapid clonal expansion. Relative stasis and occasional rapid subclone growth may characterize colorectal tumorigenesis.

Modelling the evolution of genetic instability during tumour progression

The role of genetic instability in driving carcinogenesis remains controversial. Genetic instability should accelerate carcinogenesis by increasing the rate of advantageous driver mutations; however, genetic instability can also potentially retard tumour growth by increasing the rate of deleterious mutation. As such, it is unclear whether genetically unstable clones would tend to be more selectively advantageous than their genetically stable counterparts within a growing tumour. Here, we show the circumstances where genetic instability evolves during tumour progression towards cancer. We employ a Wright-Fisher type model that describes the evolution of tumour subclones. Clones can acquire both advantageous and deleterious mutations, and mutator mutations that increase a cell's intrinsic mutation rate. Within the model, cancers evolve with a mutator phenotype when driver mutations bestow only moderate increases in fitness: very strong or weak selection for driver mutations suppresses the evolution of a mutator phenotype. Genetic instability occurs secondarily to selectively advantageous driver mutations. Deleterious mutations have relatively little effect on the evolution of genetic instability unless selection for additional driver mutations is very weak or if deleterious mutations are very common. Our model provides a framework for studying the evolution of genetic instability in tumour progression. Our analysis highlights the central role of selection in shaping patterns of mutation in carcinogenesis.