Targetable NOTCH1 rearrangements in reninoma

Reninomas are exceedingly rare renin-secreting kidney tumours that derive from juxtaglomerular cells, specialised smooth muscle cells that reside at the vascular inlet of glomeruli. They are the central component of the juxtaglomerular apparatus which controls systemic blood pressure through the secretion of renin. We assess somatic changes in reninoma and find structural variants that generate canonical activating rearrangements of, NOTCH1 whilst removing its negative regulator, NRARP. Accordingly, in single reninoma nuclei we observe excessive renin and NOTCH1 signalling mRNAs, with a concomitant non-excess of NRARP expression. Re-analysis of previously published reninoma bulk transcriptomes further corroborates our observation of dysregulated Notch pathway signalling in reninoma. Our findings reveal NOTCH1 rearrangements in reninoma, therapeutically targetable through existing NOTCH1 inhibitors, and indicate that unscheduled Notch signalling may be a disease-defining feature of reninoma.

De novo detection of somatic mutations in high-throughput single-cell profiling data sets

Characterization of somatic mutations at single-cell resolution is essential to study cancer evolution, clonal mosaicism and cell plasticity. Here, we describe SComatic, an algorithm designed for the detection of somatic mutations in single-cell transcriptomic and ATAC-seq (assay for transposase-accessible chromatin sequence) data sets directly without requiring matched bulk or single-cell DNA sequencing data. SComatic distinguishes somatic mutations from polymorphisms, RNA-editing events and artefacts using filters and statistical tests parameterized on non-neoplastic samples. Using >2.6 million single cells from 688 single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) data sets spanning cancer and non-neoplastic samples, we show that SComatic detects mutations in single cells accurately, even in differentiated cells from polyclonal tissues that are not amenable to mutation detection using existing methods. Validated against matched genome sequencing and scRNA-seq data, SComatic achieves F1 scores between 0.6 and 0.7 across diverse data sets, in comparison to 0.2-0.4 for the second-best performing method. In summary, SComatic permits de novo mutational signature analysis, and the study of clonal heterogeneity and mutational burdens at single-cell resolution.

Mapping single-cell transcriptomes in the intra-tumoral and associated territories of kidney cancer

Tumor behavior is intricately dependent on the oncogenic properties of cancer cells and their multi-cellular interactions. To understand these dependencies within the wider microenvironment, we studied over 270,000 single-cell transcriptomes and 100 microdissected whole exomes from 12 patients with kidney tumors, prior to validation using spatial transcriptomics. Tissues were sampled from multiple regions of the tumor core, the tumor-normal interface, normal surrounding tissues, and peripheral blood. We find that the tissue-type location of CD8+ T cell clonotypes largely defines their exhaustion state with intra-tumoral spatial heterogeneity that is not well explained by somatic heterogeneity. De novo mutation calling from single-cell RNA-sequencing data allows us to broadly infer the clonality of stromal cells and lineage-trace myeloid cell development. We report six conserved meta-programs that distinguish tumor cell function, and find an epithelial-mesenchymal transition meta-program highly enriched at the tumor-normal interface that co-localizes with IL1B-expressing macrophages, offering a potential therapeutic target.

A Phase II study of neoadjuvant axitinib for reducing the extent of venous tumour thrombus in clear cell renal cell cancer with venous invasion (NAXIVA)

BACKGROUND

Surgery for renal cell carcinoma (RCC) with venous tumour thrombus (VTT) extension into the renal vein (RV) and/or inferior vena cava (IVC) has high peri-surgical morbidity/mortality. NAXIVA assessed the response of VTT to axitinib, a potent tyrosine kinase inhibitor.

METHODS

NAXIVA was a single-arm, multi-centre, Phase 2 study. In total, 20 patients with resectable clear cell RCC and VTT received upto 8 weeks of pre-surgical axitinib. The primary endpoint was percentage of evaluable patients with VTT improvement by Mayo level on MRI. Secondary endpoints were percentage change in surgical approach and VTT length, response rate (RECISTv1.1) and surgical morbidity.

RESULTS

In all, 35% (7/20) patients with VTT had a reduction in Mayo level with axitinib: 37.5% (6/16) with IVC VTT and 25% (1/4) with RV-only VTT. No patients had an increase in Mayo level. In total, 75% (15/20) of patients had a reduction in VTT length. Overall, 41.2% (7/17) of patients who underwent surgery had less invasive surgery than originally planned. Non-responders exhibited lower baseline microvessel density (CD31), higher Ki67 and exhausted or regulatory T-cell phenotype.

CONCLUSIONS

NAXIVA provides the first Level II evidence that axitinib downstages VTT in a significant proportion of patients leading to reduction in the extent of surgery.

CLINICAL TRIAL REGISTRATION

NCT03494816.

Clonal diversification and histogenesis of malignant germ cell tumours

Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells. They are uniquely able to recapitulate embryonic and extraembryonic tissues, which carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we study the relationship of histogenesis and clonal diversification in GCTs by analysing the genomes and transcriptomes of 547 microdissected histological units. We find no correlation between genomic and histological heterogeneity. However, we identify unifying features including the retention of fetal developmental transcripts across tissues, expression changes on chromosome 12p, and a conserved somatic evolutionary sequence of whole genome duplication followed by clonal diversification. While this pattern is preserved across all GCTs, the developmental timing of the duplication varies between prepubertal and postpubertal cases. In addition, tumours of younger children exhibit distinct substitution signatures which may lend themselves as potential biomarkers for risk stratification. Our findings portray the extensive diversification of GCT tissues and genetic subclones as randomly distributed, while identifying overarching transcriptional and genomic features.

Hyperpolarized 13C-Pyruvate Metabolism as a Surrogate for Tumor Grade and Poor Outcome in Renal Cell Carcinoma-A Proof of Principle Study

Differentiating aggressive clear cell renal cell carcinoma (ccRCC) from indolent lesions is challenging using conventional imaging. This work prospectively compared the metabolic imaging phenotype of renal tumors using carbon-13 MRI following injection of hyperpolarized [1-13C]pyruvate (HP-13C-MRI) and validated these findings with histopathology. Nine patients with treatment-naïve renal tumors (6 ccRCCs, 1 liposarcoma, 1 pheochromocytoma, 1 oncocytoma) underwent pre-operative HP-13C-MRI and conventional proton (1H) MRI. Multi-regional tissue samples were collected using patient-specific 3D-printed tumor molds for spatial registration between imaging and molecular analysis. The apparent exchange rate constant (kPL) between 13C-pyruvate and 13C-lactate was calculated. Immunohistochemistry for the pyruvate transporter (MCT1) from 44 multi-regional samples, as well as associations between MCT1 expression and outcome in the TCGA-KIRC dataset, were investigated. Increasing kPL in ccRCC was correlated with increasing overall tumor grade (ρ = 0.92, p = 0.009) and MCT1 expression (r = 0.89, p = 0.016), with similar results acquired from the multi-regional analysis. Conventional 1H-MRI parameters did not discriminate tumor grades. The correlation between MCT1 and ccRCC grade was confirmed within a TCGA dataset (p < 0.001), where MCT1 expression was a predictor of overall and disease-free survival. In conclusion, metabolic imaging using HP-13C-MRI differentiates tumor aggressiveness in ccRCC and correlates with the expression of MCT1, a predictor of survival. HP-13C-MRI may non-invasively characterize metabolic phenotypes within renal cancer.

The WIRE study a phase II, multi-arm, multi-centre, non-randomised window-of-opportunity clinical trial platform using a Bayesian adaptive design for proof-of-mechanism of novel treatment strategies in operable renal cell cancer - a study protocol

BACKGROUND

Window-of-opportunity trials, evaluating the engagement of drugs with their biological target in the time period between diagnosis and standard-of-care treatment, can help prioritise promising new systemic treatments for later-phase clinical trials. Renal cell carcinoma (RCC), the 7th commonest solid cancer in the UK, exhibits targets for multiple new systemic anti-cancer agents including DNA damage response inhibitors, agents targeting vascular pathways and immune checkpoint inhibitors. Here we present the trial protocol for the WIndow-of-opportunity clinical trial platform for evaluation of novel treatment strategies in REnal cell cancer (WIRE).

METHODS

WIRE is a Phase II, multi-arm, multi-centre, non-randomised, proof-of-mechanism (single and combination investigational medicinal product [IMP]), platform trial using a Bayesian adaptive design. The Bayesian adaptive design leverages outcome information from initial participants during pre-specified interim analyses to determine and minimise the number of participants required to demonstrate efficacy or futility. Patients with biopsy-proven, surgically resectable, cT1b+, cN0-1, cM0-1 clear cell RCC and no contraindications to the IMPs are eligible to participate. Participants undergo diagnostic staging CT and renal mass biopsy followed by treatment in one of the treatment arms for at least 14 days. Initially, the trial includes five treatment arms with cediranib, cediranib + olaparib, olaparib, durvalumab and durvalumab + olaparib. Participants undergo a multiparametric MRI before and after treatment. Vascularised and de-vascularised tissue is collected at surgery. A ≥ 30% increase in CD8+ T-cells on immunohistochemistry between the screening and nephrectomy is the primary endpoint for durvalumab-containing arms. Meanwhile, a reduction in tumour vascular permeability measured by Ktrans on dynamic contrast-enhanced MRI by ≥30% is the primary endpoint for other arms. Secondary outcomes include adverse events and tumour size change. Exploratory outcomes include biomarkers of drug mechanism and treatment effects in blood, urine, tissue and imaging.

DISCUSSION

WIRE is the first trial using a window-of-opportunity design to demonstrate pharmacological activity of novel single and combination treatments in RCC in the pre-surgical space. It will provide rationale for prioritising promising treatments for later phase trials and support the development of new biomarkers of treatment effect with its extensive translational agenda.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT03741426 / EudraCT: 2018-003056-21 .

The mutational landscape of human somatic and germline cells

Over the course of an individual's lifetime, normal human cells accumulate mutations1. Here we compare the mutational landscape in 29 cell types from the soma and germline using multiple samples from the same individuals. Two ubiquitous mutational signatures, SBS1 and SBS5/40, accounted for the majority of acquired mutations in most cell types, but their absolute and relative contributions varied substantially. SBS18, which potentially reflects oxidative damage2, and several additional signatures attributed to exogenous and endogenous exposures contributed mutations to subsets of cell types. The rate of mutation was lowest in spermatogonia, the stem cells from which sperm are generated and from which most genetic variation in the human population is thought to originate. This was due to low rates of ubiquitous mutational processes and may be partially attributable to a low rate of cell division in basal spermatogonia. These results highlight similarities and differences in the maintenance of the germline and soma.

Extensive phylogenies of human development inferred from somatic mutations

Starting from the zygote, all cells in the human body continuously acquire mutations. Mutations shared between different cells imply a common progenitor and are thus naturally occurring markers for lineage tracing1,2. Here we reconstruct extensive phylogenies of normal tissues from three adult individuals using whole-genome sequencing of 511 laser capture microdissections. Reconstructed embryonic progenitors in the same generation of a phylogeny often contribute to different extents to the adult body. The degree of this asymmetry varies between individuals, with ratios between the two reconstructed daughter cells of the zygote ranging from 60:40 to 93:7. Asymmetries pervade subsequent generations and can differ between tissues in the same individual. The phylogenies resolve the spatial embryonic patterning of tissues, revealing contiguous patches of, on average, 301 crypts in the adult colonic epithelium derived from a most recent embryonic cell and also a spatial effect in brain development. Using data from ten additional men, we investigated the developmental split between soma and germline, with results suggesting an extraembryonic contribution to primordial germ cells. This research demonstrates that, despite reaching the same ultimate tissue patterns, early bottlenecks and lineage commitments lead to substantial variation in embryonic patterns both within and between individuals.

Single cell derived mRNA signals across human kidney tumors

Tumor cells may share some patterns of gene expression with their cell of origin, providing clues into the differentiation state and origin of cancer. Here, we study the differentiation state and cellular origin of 1300 childhood and adult kidney tumors. Using single cell mRNA reference maps of normal tissues, we quantify reference "cellular signals" in each tumor. Quantifying global differentiation, we find that childhood tumors exhibit fetal cellular signals, replacing the presumption of "fetalness" with a quantitative measure of immaturity. By contrast, in adult cancers our assessment refutes the suggestion of dedifferentiation towards a fetal state in most cases. We find an intimate connection between developmental mesenchymal populations and childhood renal tumors. We demonstrate the diagnostic potential of our approach with a case study of a cryptic renal tumor. Our findings provide a cellular definition of human renal tumors through an approach that is broadly applicable to human cancer.

Microvesicles released from pneumolysin-stimulated lung epithelial cells carry mitochondrial cargo and suppress neutrophil oxidative burst

Microvesicles (MVs) are cell-derived extracellular vesicles that have emerged as markers and mediators of acute lung injury (ALI). One of the most common pathogens in pneumonia-induced ALI is Streptococcus pneumoniae (Spn), but the role of MVs during Spn lung infection is largely unknown. In the first line of defense against Spn and its major virulence factor, pneumolysin (PLY), are the alveolar epithelial cells (AEC). In this study, we aim to characterize MVs shed from PLY-stimulated AEC and explore their contribution in mediating crosstalk with neutrophils. Using in vitro cell and ex vivo (human lung tissue) models, we demonstrated that Spn in a PLY-dependent manner stimulates AEC to release increased numbers of MVs. Spn infected mice also had higher levels of epithelial-derived MVs in their alveolar compartment compared to control. Furthermore, MVs released from PLY-stimulated AEC contain mitochondrial content and can be taken up by neutrophils. These MVs then suppress the ability of neutrophils to produce reactive oxygen species, a critical host-defense mechanism. Taken together, our results demonstrate that AEC in response to pneumococcal PLY release MVs that carry mitochondrial cargo and suggest that these MVs regulate innate immune responses during lung injury.

Lineage-Independent Tumors in Bilateral Neuroblastoma

Childhood tumors that occur synchronously in different anatomical sites usually represent metastatic disease. However, such tumors can be independent neoplasms. We investigated whether cases of bilateral neuroblastoma represented independent tumors in two children with pathogenic germline mutations by genotyping somatic mutations shared between tumors and blood. Our results suggested that in both children, the lineages that had given rise to the tumors had segregated within the first cell divisions of the zygote, without being preceded by a common premalignant clone. In one patient, the tumors had parallel evolution, including distinct second hits in SMARCA4, a putative predisposition gene for neuroblastoma. These findings portray cases of bilateral neuroblastoma as having independent lesions mediated by a germline predisposition. (Funded by Children with Cancer UK and Wellcome.).

Extensive heterogeneity in somatic mutation and selection in the human bladder

The extent of somatic mutation and clonal selection in the human bladder remains unknown. We sequenced 2097 bladder microbiopsies from 20 individuals using targeted (n = 1914 microbiopsies), whole-exome (n = 655), and whole-genome (n = 88) sequencing. We found widespread positive selection in 17 genes. Chromatin remodeling genes were frequently mutated, whereas mutations were absent in several major bladder cancer genes. There was extensive interindividual variation in selection, with different driver genes dominating the clonal landscape across individuals. Mutational signatures were heterogeneous across clones and individuals, which suggests differential exposure to mutagens in the urine. Evidence of APOBEC mutagenesis was found in 22% of the microbiopsies. Sequencing multiple microbiopsies from five patients with bladder cancer enabled comparisons with cancer-free individuals and across histological features. This study reveals a rich landscape of mutational processes and selection in normal urothelium with large heterogeneity across clones and individuals.

Three-Dimensional Printed Molds for Image-Guided Surgical Biopsies: An Open Source Computational Platform

PURPOSE

Spatial heterogeneity of tumors is a major challenge in precision oncology. The relationship between molecular and imaging heterogeneity is still poorly understood because it relies on the accurate coregistration of medical images and tissue biopsies. Tumor molds can guide the localization of biopsies, but their creation is time consuming, technologically challenging, and difficult to interface with routine clinical practice. These hurdles have so far hindered the progress in the area of multiscale integration of tumor heterogeneity data.

METHODS

We have developed an open-source computational framework to automatically produce patient-specific 3-dimensional-printed molds that can be used in the clinical setting. Our approach achieves accurate coregistration of sampling location between tissue and imaging, and integrates seamlessly with clinical, imaging, and pathology workflows.

RESULTS

We applied our framework to patients with renal cancer undergoing radical nephrectomy. We created personalized molds for 6 patients, obtaining Dice similarity coefficients between imaging and tissue sections ranging from 0.86 to 0.96 for tumor regions and between 0.70 and 0.76 for healthy kidneys. The framework required minimal manual intervention, producing the final mold design in just minutes, while automatically taking into account clinical considerations such as a preference for specific cutting planes.

CONCLUSION

Our work provides a robust and automated interface between imaging and tissue samples, enabling the development of clinical studies to probe tumor heterogeneity on multiple spatial scales.

Planar Airy beam light-sheet for two-photon microscopy

We demonstrate the first planar Airy light-sheet microscope. Fluorescence light-sheet microscopy has become the method of choice to study large biological samples with cellular or sub-cellular resolution. The propagation-invariant Airy beam enables a ten-fold increase in field-of-view with single-photon excitation; however, the characteristic asymmetry of the light-sheet limits its potential for multi-photon excitation. Here we show how a planar light-sheet can be formed from the curved propagation-invariant Airy beam. The resulting symmetric light sheet excites two-photon fluorescence uniformly across an extended field-of-view without the need for deconvolution. We demonstrate the method for rapid two-photon imaging of large volumes of neuronal tissue.

Timing the initiation of multiple myeloma

The evolution and progression of multiple myeloma and its precursors over time is poorly understood. Here, we investigate the landscape and timing of mutational processes shaping multiple myeloma evolution in a large cohort of 89 whole genomes and 973 exomes. We identify eight processes, including a mutational signature caused by exposure to melphalan. Reconstructing the chronological activity of each mutational signature, we estimate that the initial transformation of a germinal center B-cell usually occurred during the first 2nd-3rd decades of life. We define four main patterns of activation-induced deaminase (AID) and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) mutagenesis over time, including a subset of patients with evidence of prolonged AID activity during the pre-malignant phase, indicating antigen-responsiveness and germinal center reentry. Our findings provide a framework to study the etiology of multiple myeloma and explore strategies for prevention and early detection.

Mechanisms generating cancer genome complexity from a single cell division error

The chromosome breakage-fusion-bridge (BFB) cycle is a mutational process that produces gene amplification and genome instability. Signatures of BFB cycles can be observed in cancer genomes alongside chromothripsis, another catastrophic mutational phenomenon. We explain this association by elucidating a mutational cascade that is triggered by a single cell division error-chromosome bridge formation-that rapidly increases genomic complexity. We show that actomyosin forces are required for initial bridge breakage. Chromothripsis accumulates, beginning with aberrant interphase replication of bridge DNA. A subsequent burst of DNA replication in the next mitosis generates extensive DNA damage. During this second cell division, broken bridge chromosomes frequently missegregate and form micronuclei, promoting additional chromothripsis. We propose that iterations of this mutational cascade generate the continuing evolution and subclonal heterogeneity characteristic of many human cancers.

The mutational landscape of normal human endometrial epithelium

All normal somatic cells are thought to acquire mutations, but understanding of the rates, patterns, causes and consequences of somatic mutations in normal cells is limited. The uterine endometrium adopts multiple physiological states over a lifetime and is lined by a gland-forming epithelium1,2. Here, using whole-genome sequencing, we show that normal human endometrial glands are clonal cell populations with total mutation burdens that increase at about 29 base substitutions per year and that are many-fold lower than those of endometrial cancers. Normal endometrial glands frequently carry 'driver' mutations in cancer genes, the burden of which increases with age and decreases with parity. Cell clones with drivers often originate during the first decades of life and subsequently progressively colonize the epithelial lining of the endometrium. Our results show that mutational landscapes differ markedly between normal tissues-perhaps shaped by differences in their structure and physiology-and indicate that the procession of neoplastic change that leads to endometrial cancer is initiated early in life.

Comprehensive characterization of cell-free tumor DNA in plasma and urine of patients with renal tumors

BACKGROUND

Cell-free tumor-derived DNA (ctDNA) allows non-invasive monitoring of cancers, but its utility in renal cell cancer (RCC) has not been established.

METHODS

Here, a combination of untargeted and targeted sequencing methods, applied to two independent cohorts of patients (n = 91) with various renal tumor subtypes, were used to determine ctDNA content in plasma and urine.

RESULTS

Our data revealed lower plasma ctDNA levels in RCC relative to other cancers of similar size and stage, with untargeted detection in 27.5% of patients from both cohorts. A sensitive personalized approach, applied to plasma and urine from select patients (n = 22) improved detection to ~ 50%, including in patients with early-stage disease and even benign lesions. Detection in plasma, but not urine, was more frequent amongst patients with larger tumors and in those patients with venous tumor thrombus. With data from one extensively characterized patient, we observed that plasma and, for the first time, urine ctDNA may better represent tumor heterogeneity than a single tissue biopsy. Furthermore, in a subset of patients (n = 16), longitudinal sampling revealed that ctDNA can track disease course and may pre-empt radiological identification of minimal residual disease or disease progression on systemic therapy. Additional datasets will be required to validate these findings.

CONCLUSIONS

These data highlight RCC as a ctDNA-low malignancy. The biological reasons for this are yet to be determined. Nonetheless, our findings indicate potential clinical utility in the management of patients with renal tumors, provided improvement in isolation and detection approaches.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

The evolutionary history of 2,658 cancers

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Embryonal precursors of Wilms tumor

Adult cancers often arise from premalignant clonal expansions. Whether the same is true of childhood tumors has been unclear. To investigate whether Wilms tumor (nephroblastoma; a childhood kidney cancer) develops from a premalignant background, we examined the phylogenetic relationship between tumors and corresponding normal tissues. In 14 of 23 cases studied (61%), we found premalignant clonal expansions in morphologically normal kidney tissues that preceded tumor development. These clonal expansions were defined by somatic mutations shared between tumor and normal tissues but absent from blood cells. We also found hypermethylation of the H19 locus, a known driver of Wilms tumor development, in 58% of the expansions. Phylogenetic analyses of bilateral tumors indicated that clonal expansions can evolve before the divergence of left and right kidney primordia. These findings reveal embryonal precursors from which unilateral and multifocal cancers develop.

Spatiotemporal immune zonation of the human kidney

Tissue-resident immune cells are important for organ homeostasis and defense. The epithelium may contribute to these functions directly or by cross-talk with immune cells. We used single-cell RNA sequencing to resolve the spatiotemporal immune topology of the human kidney. We reveal anatomically defined expression patterns of immune genes within the epithelial compartment, with antimicrobial peptide transcripts evident in pelvic epithelium in the mature, but not fetal, kidney. A network of tissue-resident myeloid and lymphoid immune cells was evident in both fetal and mature kidney, with postnatal acquisition of transcriptional programs that promote infection-defense capabilities. Epithelial-immune cross-talk orchestrated localization of antibacterial macrophages and neutrophils to the regions of the kidney most susceptible to infection. Overall, our study provides a global overview of how the immune landscape of the human kidney is zonated to counter the dominant immunological challenge.

Genomic landscape and chronological reconstruction of driver events in multiple myeloma

The multiple myeloma (MM) genome is heterogeneous and evolves through preclinical and post-diagnosis phases. Here we report a catalog and hierarchy of driver lesions using sequences from 67 MM genomes serially collected from 30 patients together with public exome datasets. Bayesian clustering defines at least 7 genomic subgroups with distinct sets of co-operating events. Focusing on whole genome sequencing data, complex structural events emerge as major drivers, including chromothripsis and a novel replication-based mechanism of templated insertions, which typically occur early. Hyperdiploidy also occurs early, with individual trisomies often acquired in different chronological windows during evolution, and with a preferred order of acquisition. Conversely, positively selected point mutations, whole genome duplication and chromoplexy events occur in later disease phases. Thus, initiating driver events, drawn from a limited repertoire of structural and numerical chromosomal changes, shape preferred trajectories of evolution that are biologically relevant but heterogeneous across patients.

Interactive molecular dynamics in virtual reality from quantum chemistry to drug binding: An open-source multi-person framework

As molecular scientists have made progress in their ability to engineer nanoscale molecular structure, we face new challenges in our ability to engineer molecular dynamics (MD) and flexibility. Dynamics at the molecular scale differs from the familiar mechanics of everyday objects because it involves a complicated, highly correlated, and three-dimensional many-body dynamical choreography which is often nonintuitive even for highly trained researchers. We recently described how interactive molecular dynamics in virtual reality (iMD-VR) can help to meet this challenge, enabling researchers to manipulate real-time MD simulations of flexible structures in 3D. In this article, we outline various efforts to extend immersive technologies to the molecular sciences, and we introduce "Narupa," a flexible, open-source, multiperson iMD-VR software framework which enables groups of researchers to simultaneously cohabit real-time simulation environments to interactively visualize and manipulate the dynamics of molecular structures with atomic-level precision. We outline several application domains where iMD-VR is facilitating research, communication, and creative approaches within the molecular sciences, including training machines to learn potential energy functions, biomolecular conformational sampling, protein-ligand binding, reaction discovery using "on-the-fly" quantum chemistry, and transport dynamics in materials. We touch on iMD-VR's various cognitive and perceptual affordances and outline how these provide research insight for molecular systems. By synergistically combining human spatial reasoning and design insight with computational automation, technologies such as iMD-VR have the potential to improve our ability to understand, engineer, and communicate microscopic dynamical behavior, offering the potential to usher in a new paradigm for engineering molecules and nano-architectures.

Pilot study of head conformation changes over time in the Cavalier King Charles spaniel breed

Modern interpretation of head conformation in the Cavalier King Charles spaniel (CKCS) has favoured a smaller, more exaggerated, brachycephalic type than originally described in the 1929 breed standard. Recent research studies identified brachycephaly and reduced hind cranium as two conformational (dysmorphic) features that increase risk for symptomatic Chiari-like malformation and secondary syringomyelia (SM). A prospective pilot study investigated the hypothesis that dysmorphic head features could be assessed visually and correlated with risk of SM. Thirteen CKCS, selected from anonymised photographic evidence, were physically appraised by authorised Kennel Club judges using a head shape checklist. These subjective evaluations were then matched with objective measurements of the cranium (cephalic index and rostrocaudal doming) and their subsequent MRI. A positive correlation (P=0.039) between the judges’ checklist score and rostrocaudal doming (hindskull ratio) and a positive correlation between the cephalic index and hindskull ratio (P=0.042) were identified. Five CKCS had no SM and their status tallied with 62 per cent of the judges’ evaluation. Although the ability of adjudicators to identify differences in head conformation varied, there was sufficient association between the dysmorphic parameters and the risk of SM to cause concern and propose a larger study in CKCS breed.

Genomics and clinical correlates of renal cell carcinoma

PURPOSE

Clear cell, papillary cell, and chromophobe renal cell carcinomas (RCCs) have now been well characterised thanks to large collaborative projects such as The Cancer Genome Atlas (TCGA). Not only has knowledge of the genomic landscape helped inform the development of new drugs, it also promises to fine tune prognostication.

METHODS

A literature review was performed summarising the current knowledge on the genetic basis of RCC.

RESULTS

The Von Hippel-Lindau (VHL) tumour suppressor gene undergoes bi-allelic knockout in the vast majority of clear cell RCCs. The next most prevalent aberrations include a cohort of chromatin-modifying genes with diverse roles including PBRM1, SETD2, BAP1, and KMD5C. The most common non-clear cell renal cancers have also undergone genomic profiling and are characterised by distinct genomic landscapes. Many recurrent mutations have prognostic value and show promise in aiding decisions regarding treatment stratification. Intra-tumour heterogeneity appears to hamper the clinical applicability of sparsely sampled tumours. Ways to abrogate heterogeneity will be required to optimise the genomic classification of tumours.

CONCLUSION

The somatic mutational landscape of the more common renal cancers is well known. Correlation with outcome needs to be more comprehensively furnished, particularly for small renal masses, rarer non-clear cell renal cancers, and for all tumours undergoing targeted therapy.

Single-cell transcriptomes from human kidneys reveal the cellular identity of renal tumors

Messenger RNA encodes cellular function and phenotype. In the context of human cancer, it defines the identities of malignant cells and the diversity of tumor tissue. We studied 72,501 single-cell transcriptomes of human renal tumors and normal tissue from fetal, pediatric, and adult kidneys. We matched childhood Wilms tumor with specific fetal cell types, thus providing evidence for the hypothesis that Wilms tumor cells are aberrant fetal cells. In adult renal cell carcinoma, we identified a canonical cancer transcriptome that matched a little-known subtype of proximal convoluted tubular cell. Analyses of the tumor composition defined cancer-associated normal cells and delineated a complex vascular endothelial growth factor (VEGF) signaling circuit. Our findings reveal the precise cellular identities and compositions of human kidney tumors.

Abstract 218: The evolutionary history of 2,658 cancers

Cancer develops through a continuous process of somatic evolution. Whole genome sequencing provides a snapshot of the tumor genome at the point of sampling, however, the data can contain information that permits the reconstruction of a tumor's evolutionary past.

Here, we apply such life history analyses on an unprecedented scale, to a set of 2,658 tumors spanning 39 cancer types. We estimated the timing of large chromosomal gains during tumor evolution, by comparing the rates of doubled to non-doubled point mutations within gained regions. Although we find that such events typically occur in the second half of clonal evolution, we also observe distinctive and early chromosomal gains in some cancer types, such as gains of chromosomes 7, 19 and 20 in glioblastoma, and isochromosome 17q in medulloblastoma. By integrating these results with the qualitative timing of individual driver mutations, we obtained an overall ranking, from early to late, of frequent somatic events per cancer type, which both identified novel patterns of tumor evolution, and incorporated additional detail into known models, such as the progression of APC-KRAS-TP53 in colorectal cancer proposed by Vogelstein and Fearon.

To estimate how mutational processes acting on the tumor genome change over time, we classified mutations in each sample according to three broad time periods (early clonal, late clonal, and subclonal), and quantified the activity of mutational signatures in each period. Most mutational processes appear to remain remarkably constant, however, certain signatures show clear and consistent changes during clonal evolution. Particularly, mutational signatures associated with exposure to carcinogens, such as smoking and UV light, tend to decrease over time. In contrast, signatures associated with defective endogenous processes, such as APOBEC mutagenesis and defective double strand break repair, show an increase between early and late phases of tumor evolution.

Making use of clock-like mutational signatures, we converted mutational time estimates for large events, such as whole genome duplication (WGD), and the emergence of the most recent common ancestor (MRCA), into real time estimates, which allowed us to combine our analyses into overall timelines of cancer evolution, per tumor type. For example, the typical timeline of ovarian adenocarcinoma development shows that early tumor evolution is characterized by mutations in TP53, and widespread genome instability, with WGD events taking place on average 8 years prior to diagnosis. In later stages of evolution, signatures of defective repair processes increase, and the MRCA emerges on average 1 year before diagnosis.

Taken together, these data reveal the common and divergent evolutionary trajectories available to a cancer, which might be crucial in understanding specific tumor biology, and in providing new opportunities for early detection and cancer prevention.

Citation Format: Clemency Jolly, Moritz Gerstung, Ignaty Leshchiner, Stefan C. Dentro, Santiago Gonzalez, Thomas J. Mitchell, Yulia Rubanova, Pavana Anur, Daniel Rosebrock, Kaixian Yu, Maxime Tarabichi, Amit Deshwar, Jeff Wintersinger, Kortine Kleinheinz, Ignacio Vásquez-García, Kerstin Haase, Subhajit Sengupta, Geoff Macintyre, Salem Malikic, Nilgun Donmez, Dimitri G. Livitz, Mark Cmero, Jonas Demeulemeester, Steve Schumacher, Yu Fan, Xiaotong Yao, Juhee Lee, Matthias Schlesner, Paul C. Boutros, David D. Bowtell, Hongtu Zhu, Gad Getz, Marcin Imielinski, Rameen Beroukhim, S Cenk Sahinalp, Yuan Ji, Martin Peifer, Florian Markowetz, Ville Mustonen, Ke Juan, Wenyi Wang, Quaid D. Morris, Paul T. Spellman, David C. Wedge, Peter Van Loo, PCAWG Evolution and Heterogeneity Working Group. The evolutionary history of 2,658 cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 218.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

Pneumolysin activates neutrophil extracellular trap formation

The primary objective of the current study was to investigate the potential of the pneumococcal toxin, pneumolysin (Ply), to activate neutrophil extracellular trap (NET) formation in vitro. Isolated human blood neutrophils were exposed to recombinant Ply (5-20 ng ml(-1) ) for 30-90 min at 37°C and NET formation measured using the following procedures to detect extracellular DNA: (i) flow cytometry using Vybrant® DyeCycle™ Ruby; (ii) spectrofluorimetry using the fluorophore, Sytox(®) Orange (5 μM); and (iii) NanoDrop(®) technology. These procedures were complemented by fluorescence microscopy using 4', 6-diamino-2-phenylindole (DAPI) (nuclear stain) in combination with anti-citrullinated histone monoclonal antibodies to visualize nets. Exposure of neutrophils to Ply resulted in relatively rapid (detected within 30-60 min), statistically significant (P < 0·05) dose- and time-related increases in the release of cellular DNA impregnated with both citrullinated histone and myeloperoxidase. Microscopy revealed that NETosis appeared to be restricted to a subpopulation of neutrophils, the numbers of NET-forming cells in the control and Ply-treated systems (10 and 20 ng ml(-1) ) were 4·3 (4·2), 14.3 (9·9) and 16·5 (7·5), respectively (n = 4, P < 0·0001 for comparison of the control with both Ply-treated systems). Ply-induced NETosis occurred in the setting of retention of cell viability, and apparent lack of involvement of reactive oxygen species and Toll-like receptor 4. In conclusion, Ply induces vital NETosis in human neutrophils, a process which may either contribute to host defence or worsen disease severity, depending on the intensity of the inflammatory response during pneumococcal infection.

Quantitative high dynamic range beam profiling for fluorescence microscopy

Modern developmental biology relies on optically sectioning fluorescence microscope techniques to produce non-destructive in vivo images of developing specimens at high resolution in three dimensions. As optimal performance of these techniques is reliant on the three-dimensional (3D) intensity profile of the illumination employed, the ability to directly record and analyze these profiles is of great use to the fluorescence microscopist or instrument builder. Though excitation beam profiles can be measured indirectly using a sample of fluorescent beads and recording the emission along the microscope detection path, we demonstrate an alternative approach where a miniature camera sensor is used directly within the illumination beam. Measurements taken using our approach are solely concerned with the illumination optics as the detection optics are not involved. We present a miniature beam profiling device and high dynamic range flux reconstruction algorithm that together are capable of accurately reproducing quantitative 3D flux maps over a large focal volume. Performance of this beam profiling system is verified within an optical test bench and demonstrated for fluorescence microscopy by profiling the low NA illumination beam of a single plane illumination microscope. The generality and success of this approach showcases a widely flexible beam amplitude diagnostic tool for use within the life sciences.

Syringomyelia: determining risk and protective factors in the conformation of the Cavalier King Charles Spaniel dog

BACKGROUND

Syringomyelia (SM) is a painful neurological condition, prevalent in brachycephalic toy breeds including the Cavalier King Charles Spaniel (CKCS). In these breeds, SM is typically secondary to Chiari-like Malformation (CM). There has been much debate in the scientific and veterinary communities to what extent head shape is indicative of either pathology, especially as certain craniosynostosis syndromes in humans (highly associated with CM) have characteristic facial and cranial morphologies. Elucidating a risk morphology would allow for selection away from these traits and proffer further breeding guidelines for the condition. Dogs were measured in multiple countries by means of a standardised bony landmark measuring protocol and photo analysis by blinded, trained researchers.

RESULTS

The results found two significant risk factors in the conformation of the CKCS: extent of brachycephaly and distribution of cranium. The study identified a greater amount of cranium distributed caudally (relative to the amount distributed rostrally) to be significantly protective against syrinx development at the levels of three years of age, five years of age and when comparing a sample of SM clear individuals over the age of five to those affected younger than three years of age. A decreased cephalic index (decreasing brachycephaly) was significantly protective at the latter level. Cephalic index and caudal cranium distribution exhibited a negative, linear relationship. Cephalic index demonstrated a positive linear relationship with the amount of doming of the head.

CONCLUSIONS

This study proposes a risk phenotype of brachycephaly with resulting rostrocaudal doming that is more rostrally distributed and hence sloping caudally. The results of this study may allow for selection against risk aspects of conformation in the CKCS in combination with the British Veterinary Association/Kennel Club CM/SM scheme to enable reduction in CM/SM incidence. Further research comparing this external risk phenotype to the internal presentation upon MRI would determine how these features are indicative of syrinx development. Utilising breeds in which CM free individuals are more available may allow for validation of this risk phenotype for CM or determine alternatives.