High- or low-dose preoperative ipilimumab plus nivolumab in stage III urothelial cancer: the phase 1B NABUCCO trial.

Cohort 1 of the phase 1B NABUCCO trial showed high pathological complete response (pCR) rates with preoperative ipilimumab plus nivolumab in stage III urothelial cancer (UC). In cohort 2, the aim was dose adjustment to optimize responses. Additionally, we report secondary endpoints, including efficacy and tolerability, in cohort 2 and the association of presurgical absence of circulating tumor DNA (ctDNA) in urine and plasma with clinical outcome in both cohorts. Thirty patients received two cycles of either ipilimumab 3 mg kg-1 plus nivolumab 1 mg kg-1 (cohort 2A) or ipilimumab 1 mg kg-1 plus nivolumab 3 mg kg-1 (cohort 2B), both followed by nivolumab 3 mg kg-1. We observed a pCR in six (43%) patients in cohort 2A and a pCR in one (7%) patient in cohort 2B. Absence of urinary ctDNA correlated with pCR in the bladder (ypT0Nx) but not with progression-free survival (PFS). Absence of plasma ctDNA correlated with pCR (odds ratio: 45.0; 95% confidence interval (CI): 4.9-416.5) and PFS (hazard ratio: 10.4; 95% CI: 2.9-37.5). Our data suggest that high-dose ipilimumab plus nivolumab is required in stage III UC and that absence of ctDNA in plasma can predict PFS. ClinicalTrials.gov registration: NCT03387761 .

Substitution mutational signatures in whole-genome-sequenced cancers in the UK population.

Whole-genome sequencing (WGS) permits comprehensive cancer genome analyses, revealing mutational signatures, imprints of DNA damage and repair processes that have arisen in each patient's cancer. We performed mutational signature analyses on 12,222 WGS tumor-normal matched pairs, from patients recruited via the UK National Health Service. We contrasted our results to two independent cancer WGS datasets, the International Cancer Genome Consortium (ICGC) and Hartwig Foundation, involving 18,640 WGS cancers in total. Our analyses add 40 single and 18 double substitution signatures to the current mutational signature tally. Critically, we show for each organ, that cancers have a limited number of 'common' signatures and a long tail of 'rare' signatures. We provide a practical solution for utilizing this concept of common versus rare signatures in future analyses.

Substantial somatic genomic variation and selection for BCOR mutations in human induced pluripotent stem cells.

We explored human induced pluripotent stem cells (hiPSCs) derived from different tissues to gain insights into genomic integrity at single-nucleotide resolution. We used genome sequencing data from two large hiPSC repositories involving 696 hiPSCs and daughter subclones. We find ultraviolet light (UV)-related damage in ~72% of skin fibroblast-derived hiPSCs (F-hiPSCs), occasionally resulting in substantial mutagenesis (up to 15 mutations per megabase). We demonstrate remarkable genomic heterogeneity between independent F-hiPSC clones derived during the same round of reprogramming due to oligoclonal fibroblast populations. In contrast, blood-derived hiPSCs (B-hiPSCs) had fewer mutations and no UV damage but a high prevalence of acquired BCOR mutations (26.9% of lines). We reveal strong selection pressure for BCOR mutations in F-hiPSCs and B-hiPSCs and provide evidence that they arise in vitro. Directed differentiation of hiPSCs and RNA sequencing showed that BCOR mutations have functional consequences. Our work strongly suggests that detailed nucleotide-resolution characterization is essential before using hiPSCs.

A systematic CRISPR screen defines mutational mechanisms underpinning signatures caused by replication errors and endogenous DNA damage.

Mutational signatures are imprints of pathophysiological processes arising through tumorigenesis. We generated isogenic CRISPR-Cas9 knockouts (Δ) of 43 genes in human induced pluripotent stem cells, cultured them in the absence of added DNA damage, and performed whole-genome sequencing of 173 subclones. ΔOGG1, ΔUNG, ΔEXO1, ΔRNF168, ΔMLH1, ΔMSH2, ΔMSH6, ΔPMS1, and ΔPMS2 produced marked mutational signatures indicative of being critical mitigators of endogenous DNA modifications. Detailed analyses revealed mutational mechanistic insights, including how 8-oxo-dG elimination is sequence-context-specific while uracil clearance is sequence-context-independent. Mismatch repair (MMR) deficiency signatures are engendered by oxidative damage (C>A transversions), differential misincorporation by replicative polymerases (T>C and C>T transitions), and we propose a 'reverse template slippage' model for T>A transversions. ΔMLH1, ΔMSH6, and ΔMSH2 signatures were similar to each other but distinct from ΔPMS2. Finally, we developed a classifier, MMRDetect, where application to 7,695 WGS cancers showed enhanced detection of MMR-deficient tumors, with implications for responsiveness to immunotherapies.

Salt dependent mesoscopic model for RNA at multiple strand concentrations.

Mesoscopic models can be used for the description of the thermodynamic properties of RNA duplexes. With the use of experimental melting temperatures, its parametrization can provide important insights into its hydrogen bonds and stacking interactions as has been done for high sodium concentrations. However, the RNA parametrization for lower salt concentrations is still missing due to the limited amount of published melting temperature data. While the Peyrard-Bishop (PB) parametrization was found to be largely independent of strand concentrations, it requires that all temperatures are provided at the same strand concentrations. Here we adapted the PB model to handle multiple strand concentrations and in this way we were able to make use of an experimental set of temperatures to model the hydrogen bond and stacking interactions at low and intermediate sodium concentrations. For the parametrizations we make a distinction between terminal and internal base pairs, and the resulting potentials were qualitatively similar as we obtained previously for DNA. The main difference from DNA parameters, was the Morse potentials at low sodium concentrations for terminal r(AU) which is stronger than d(AT), suggesting higher hydrogen bond strength.

Spectrum of mutational signatures in T-cell lymphoma reveals a key role for UV radiation in cutaneous T-cell lymphoma.

T-cell non-Hodgkin's lymphomas develop following transformation of tissue resident T-cells. We performed a meta-analysis of whole exome sequencing data from 403 patients with eight subtypes of T-cell non-Hodgkin's lymphoma to identify mutational signatures and associated recurrent gene mutations. Signature 1, indicative of age-related deamination, was prevalent across all T-cell lymphomas, reflecting the derivation of these malignancies from memory T-cells. Adult T-cell leukemia-lymphoma was specifically associated with signature 17, which was found to correlate with the IRF4 K59R mutation that is exclusive to Adult T-cell leukemia-lymphoma. Signature 7, implicating UV exposure was uniquely identified in cutaneous T-cell lymphoma (CTCL), contributing 52% of the mutational burden in mycosis fungoides and 23% in Sezary syndrome. Importantly this UV signature was observed in CD4 + T-cells isolated from the blood of Sezary syndrome patients suggesting extensive re-circulation of these T-cells through skin and blood. Analysis of non-Hodgkin's T-cell lymphoma cases submitted to the national 100,000 WGS project confirmed that signature 7 was only identified in CTCL strongly implicating UV radiation in the pathogenesis of cutaneous T-cell lymphoma.

A practical framework and online tool for mutational signature analyses show inter-tissue variation and driver dependencies.

Mutational signatures are patterns of mutations that arise during tumorigenesis. We present an enhanced, practical framework for mutational signature analyses. Applying these methods on 3,107 whole genome sequenced (WGS) primary cancers of 21 organs reveals known signatures and nine previously undescribed rearrangement signatures. We highlight inter-organ variability of signatures and present a way of visualizing that diversity, reinforcing our findings in an independent analysis of 3,096 WGS metastatic cancers. Signatures with a high level of genomic instability are dependent on TP53 dysregulation. We illustrate how uncertainty in mutational signature identification and assignment to samples affects tumor classification, reinforcing that using multiple orthogonal mutational signature data is not only beneficial, it is essential for accurate tumor stratification. Finally, we present a reference web-based tool for cancer and experimentally-generated mutational signatures, called Signal (https://signal.mutationalsignatures.com), that also supports performing mutational signature analyses.

Whole-genome sequencing of triple-negative breast cancers in a population-based clinical study.

Whole-genome sequencing (WGS) brings comprehensive insights to cancer genome interpretation. To explore the clinical value of WGS, we sequenced 254 triple-negative breast cancers (TNBCs) for which associated treatment and outcome data were collected between 2010 and 2015 via the population-based Sweden Cancerome Analysis Network-Breast (SCAN-B) project (ClinicalTrials.gov ID:NCT02306096). Applying the HRDetect mutational-signature-based algorithm to classify tumors, 59% were predicted to have homologous-recombination-repair deficiency (HRDetect-high): 67% explained by germline/somatic mutations of BRCA1/BRCA2, BRCA1 promoter hypermethylation, RAD51C hypermethylation or biallelic loss of PALB2. A novel mechanism of BRCA1 abrogation was discovered via germline SINE-VNTR-Alu retrotransposition. HRDetect provided independent prognostic information, with HRDetect-high patients having better outcome on adjuvant chemotherapy for invasive disease-free survival (hazard ratio (HR) = 0.42; 95% confidence interval (CI) = 0.2-0.87) and distant relapse-free interval (HR = 0.31, CI = 0.13-0.76) compared to HRDetect-low, regardless of whether a genetic/epigenetic cause was identified. HRDetect-intermediate, some possessing potentially targetable biological abnormalities, had the poorest outcomes. HRDetect-low cancers also had inadequate outcomes: ~4.7% were mismatch-repair-deficient (another targetable defect, not typically sought) and they were enriched for (but not restricted to) PIK3CA/AKT1 pathway abnormalities. New treatment options need to be considered for now-discernible HRDetect-intermediate and HRDetect-low categories. This population-based study advocates for WGS of TNBC to better inform trial stratification and improve clinical decision-making.

The Genomic and Immune Landscapes of Lethal Metastatic Breast Cancer.

The detailed molecular characterization of lethal cancers is a prerequisite to understanding resistance to therapy and escape from cancer immunoediting. We performed extensive multi-platform profiling of multi-regional metastases in autopsies from 10 patients with therapy-resistant breast cancer. The integrated genomic and immune landscapes show that metastases propagate and evolve as communities of clones, reveal their predicted neo-antigen landscapes, and show that they can accumulate HLA loss of heterozygosity (LOH). The data further identify variable tumor microenvironments and reveal, through analyses of T cell receptor repertoires, that adaptive immune responses appear to co-evolve with the metastatic genomes. These findings reveal in fine detail the landscapes of lethal metastatic breast cancer.

Effect of 3'UTR RET Variants on RET mRNA Secondary Structure and Disease Presentation in Medullary Thyroid Carcinoma.

BACKGROUND: The RET S836S variant has been associated with early onset and increased risk for metastatic disease in medullary thyroid carcinoma (MTC). However, the mechanism by which this variant modulates MTC pathogenesis is still open to discuss. Of interest, strong linkage disequilibrium (LD) between RET S836S and 3'UTR variants has been reported in Hirschsprung's disease patients. OBJECTIVE: To evaluate the frequency of the RET 3'UTR variants (rs76759170 and rs3026785) in MTC patients and to determine whether these variants are in LD with S836S polymorphism. METHODS: Our sample comprised 152 patients with sporadic MTC. The RET S836S and 3'UTR (rs76759170 and rs3026785) variants were genotyped using Custom TaqMan Genotyping Assays. Haplotypes were inferred using the phase 2.1 program. RET mRNA structure was assessed by Vienna Package. RESULTS: The mean age of MTC diagnosis was 48.5±15.5 years and 57.9% were women. The minor allele frequencies of RET polymorphisms were as follows: S836S, 5.6%; rs76759170, 5.6%; rs3026785, 6.2%. We observed a strong LD among S836S and 3'UTR variants (|D'| = -1, r2 = 1 and |D'| = -1, r2 = 0,967). Patients harboring the S836S/3'UTR variants presented a higher percentage of lymph node and distant metastasis (P = 0.013 and P<0.001, respectively). Accordingly, RNA folding analyses demonstrated different RNA secondary structure predictions for WT(TCCGT), S836S(TTCGT) or 3'UTR(GTCAC) haplotypes. The S836S/3'UTR haplotype presented a greater number of double helices sections and lower levels of minimal free energy when compared to the wild-type haplotype, suggesting that these variants provides the most thermodynamically stable mRNA structure, which may have functional consequences on the rate of mRNA degradation. CONCLUSION: The RET S836S polymorphism is in LD with 3'UTR variants. In silico analysis indicate that the 3'UTR variants may affect the secondary structure of RET mRNA, suggesting that these variants might play a role in posttranscriptional control of the RET transcripts.

Evaluating Hydrogen Bonds and Base Stacking of Single, Tandem and Terminal GU Mismatches in RNA with a Mesoscopic Model.

Guanine-Uracil (GU) mismatches are crucial to the stability of the RNA double helix and need to be considered in RNA folding algorithms for numerous biotechnological applications. Yet despite its importance, many aspects of GU base pairs are still poorly understood. There is also a lack of parametrization which prevents it to be considered in mesoscopic models. Here, we adapted the mesoscopic Peyrard-Bishop model to deal with context-dependent hydrogen bonds of GU mismatches and calculated the model parameters related to hydrogen bonding and base stacking from available experimental melting temperatures. The context-dependence causes a proliferation of parameters which made the problem computationally very demanding. We were able to overcome this problem by systematically regrouping the parameters during the minimization procedure. Our results not only provide the much needed parametrization but also answer several questions about the general properties of GU base pairs, as they can be associated straightforwardly to hydrogen bonding and base stacking. In particular, we found a very small Morse potential for tandem 5'-GU-3', which confirms a single hydrogen bond for this configuration, answering a long-standing question over conflicting experimental findings. Terminal GU base pairs are known to increase the duplex stability, but it is not clear why. Our results suggest that the increased terminal stability is mostly due to stronger hydrogen bonding.

DNA terminal base pairs have weaker hydrogen bonds especially for AT under low salt concentration.

DNA base pairs are known to open more easily at the helix terminal, a process usually called end fraying, the details of which are still poorly understood. Here, we present a mesoscopic model calculation based on available experimental data where we consider separately the terminal base pairs of a DNA duplex. Our results show an important reduction of hydrogen bond strength for terminal cytosine-guanine (CG) base pairs which is uniform over the whole range of salt concentrations, while for AT base pairs, we obtain a nearly 1/3 reduction but only at low salt concentrations. At higher salt concentrations, terminal adenine-thymine (AT) pair has almost the same hydrogen bond strength than interior bases. The calculated terminal stacking interaction parameters display some peculiarly contrasting behavior. While there is mostly no perceptible difference to internal stacking, for some cases, we observe an unusually strong dependence with salt concentration which does not appear follow any pattern or trend.

Effects of airway surface liquid height on the kinetics of extracellular nucleotides in airway epithelia.

Experimental techniques aimed at measuring the concentration of signaling molecules in the airway surface liquid (ASL) often require an unrealistically large ASL volume to facilitate sampling. This experimental limitation, prompted by the difficulty of pipetting liquid from a very shallow layer (~15 µm), leads to dilution and the under-prediction of physiologic concentrations of signaling molecules that are vital to the regulation of mucociliary clearance. Here, we use a computational model to describe the effect of liquid height on the kinetics of extracellular nucleotides in the airway surface liquid coating respiratory epithelia. The model consists of a reaction-diffusion equation with boundary conditions that represent the enzymatic reactions occurring on the epithelial surface. The simulations reproduce successfully the kinetics of extracellular ATP following hypotonic challenge for ASL volumes ranging from 25 µl to 500 µl in a 12-mm diameter cell culture. The model reveals that [ATP] and [ADO] reach 1200 nM and 2200 nM at the epithelial surface, respectively, while their volumetric averages remain less than 200 nM at all times in experiments with a large ASL volume (500 µl). These findings imply that activation of P2Y2 and A2B receptors is robust after hypotonic challenge, in contrast to what could be concluded based on experimental measurements of volumetric concentrations in large ASL volumes. Finally, given the central role that ATP and ADO play in regulating mucociliary clearance, we investigated which enzymes, when inhibited, provide the greatest increase in ATP and ADO concentrations. Our findings suggest that inhibition of NTPDase1/highTNAP would cause the greatest increase in [ATP] after hypotonic challenge, while inhibition of the transporter CNT3 would provide the greatest increase in [ADO].