SOX17 enables immune evasion of early colorectal adenomas and cancers.

A hallmark of cancer is the avoidance of immune destruction. This process has been primarily investigated in locally advanced or metastatic cancer1-3; however, much less is known about how pre-malignant or early invasive tumours evade immune detection. Here, to understand this process in early colorectal cancers (CRCs), we investigated how naive colon cancer organoids that were engineered in vitro to harbour Apc-null, KrasG12D and Trp53-null (AKP) mutations adapted to the in vivo native colonic environment. Comprehensive transcriptomic and chromatin analyses revealed that the endoderm-specifying transcription factor SOX17 became strongly upregulated in vivo. Notably, whereas SOX17 loss did not affect AKP organoid propagation in vitro, its loss markedly reduced the ability of AKP tumours to persist in vivo. The small fraction of SOX17-null tumours that grew displayed notable interferon-γ (IFNγ)-producing effector-like CD8+ T cell infiltrates in contrast to the immune-suppressive microenvironment in wild-type counterparts. Mechanistically, in both endogenous Apc-null pre-malignant adenomas and transplanted organoid-derived AKP CRCs, SOX17 suppresses the ability of tumour cells to sense and respond to IFNγ, preventing anti-tumour T cell responses. Finally, SOX17 engages a fetal intestinal programme that drives differentiation away from LGR5+ tumour cells to produce immune-evasive LGR5- tumour cells with lower expression of major histocompatibility complex class I (MHC-I). We propose that SOX17 is a transcription factor that is engaged during the early steps of colon cancer to orchestrate an immune-evasive programme that permits CRC initiation and progression.

Strand-resolved mutagenicity of DNA damage and repair.

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Convergence of coronary artery disease genes onto endothelial cell programs.

Linking variants from genome-wide association studies (GWAS) to underlying mechanisms of disease remains a challenge1-3. For some diseases, a successful strategy has been to look for cases in which multiple GWAS loci contain genes that act in the same biological pathway1-6. However, our knowledge of which genes act in which pathways is incomplete, particularly for cell-type-specific pathways or understudied genes. Here we introduce a method to connect GWAS variants to functions. This method links variants to genes using epigenomics data, links genes to pathways de novo using Perturb-seq and integrates these data to identify convergence of GWAS loci onto pathways. We apply this approach to study the role of endothelial cells in genetic risk for coronary artery disease (CAD), and discover 43 CAD GWAS signals that converge on the cerebral cavernous malformation (CCM) signalling pathway. Two regulators of this pathway, CCM2 and TLNRD1, are each linked to a CAD risk variant, regulate other CAD risk genes and affect atheroprotective processes in endothelial cells. These results suggest a model whereby CAD risk is driven in part by the convergence of causal genes onto a particular transcriptional pathway in endothelial cells. They highlight shared genes between common and rare vascular diseases (CAD and CCM), and identify TLNRD1 as a new, previously uncharacterized member of the CCM signalling pathway. This approach will be widely useful for linking variants to functions for other common polygenic diseases.

Structures, functions and adaptations of the human LINE-1 ORF2 protein.

The LINE-1 (L1) retrotransposon is an ancient genetic parasite that has written around one-third of the human genome through a 'copy and paste' mechanism catalysed by its multifunctional enzyme, open reading frame 2 protein (ORF2p)1. ORF2p reverse transcriptase (RT) and endonuclease activities have been implicated in the pathophysiology of cancer2,3, autoimmunity4,5 and ageing6,7, making ORF2p a potential therapeutic target. However, a lack of structural and mechanistic knowledge has hampered efforts to rationally exploit it. We report structures of the human ORF2p 'core' (residues 238-1061, including the RT domain) by X-ray crystallography and cryo-electron microscopy in several conformational states. Our analyses identified two previously undescribed folded domains, extensive contacts to RNA templates and associated adaptations that contribute to unique aspects of the L1 replication cycle. Computed integrative structural models of full-length ORF2p show a dynamic closed-ring conformation that appears to open during retrotransposition. We characterize ORF2p RT inhibition and reveal its underlying structural basis. Imaging and biochemistry show that non-canonical cytosolic ORF2p RT activity can produce RNA:DNA hybrids, activating innate immune signalling through cGAS/STING and resulting in interferon production6-8. In contrast to retroviral RTs, L1 RT is efficiently primed by short RNAs and hairpins, which probably explains cytosolic priming. Other biochemical activities including processivity, DNA-directed polymerization, non-templated base addition and template switching together allow us to propose a revised L1 insertion model. Finally, our evolutionary analysis demonstrates structural conservation between ORF2p and other RNA- and DNA-dependent polymerases. We therefore provide key mechanistic insights into L1 polymerization and insertion, shed light on the evolutionary history of L1 and enable rational drug development targeting L1.

Affinity proteomics reveals human host factors implicated in discrete stages of LINE-1 retrotransposition.

LINE-1s are active human DNA parasites that are agents of genome dynamics in evolution and disease. These streamlined elements require host factors to complete their life cycles, whereas hosts have developed mechanisms to combat retrotransposition's mutagenic effects. As such, endogenous L1 expression levels are extremely low, creating a roadblock for detailed interactomic analyses. Here, we describe a system to express and purify highly active L1 RNP complexes from human suspension cell culture and characterize the copurified proteome, identifying 37 high-confidence candidate interactors. These data sets include known interactors PABPC1 and MOV10 and, with in-cell imaging studies, suggest existence of at least three types of compositionally and functionally distinct L1 RNPs. Among the findings, UPF1, a key nonsense-mediated decay factor, and PCNA, the polymerase-delta-associated sliding DNA clamp, were identified and validated. PCNA interacts with ORF2p via a PIP box motif; mechanistic studies suggest that this occurs during or immediately after target-primed reverse transcription.

Signatures of TOP1 transcription-associated mutagenesis in cancer and germline.

The mutational landscape is shaped by many processes. Genic regions are vulnerable to mutation but are preferentially protected by transcription-coupled repair1. In microorganisms, transcription has been demonstrated to be mutagenic2,3; however, the impact of transcription-associated mutagenesis remains to be established in higher eukaryotes4. Here we show that ID4-a cancer insertion-deletion (indel) mutation signature of unknown aetiology5 characterized by short (2 to 5 base pair) deletions -is due to a transcription-associated mutagenesis process. We demonstrate that defective ribonucleotide excision repair in mammals is associated with the ID4 signature, with mutations occurring at a TNT sequence motif, implicating topoisomerase 1 (TOP1) activity at sites of genome-embedded ribonucleotides as a mechanistic basis. Such TOP1-mediated deletions occur somatically in cancer, and the ID-TOP1 signature is also found in physiological settings, contributing to genic de novo indel mutations in the germline. Thus, although topoisomerases protect against genome instability by relieving topological stress6, their activity may also be an important source of mutations in the human genome.

Identification of a localized nonsense-mediated decay pathway at the endoplasmic reticulum.

Nonsense-mediated decay (NMD) is a translation-dependent RNA quality control mechanism that occurs in the cytoplasm. However, it is unknown how NMD regulates the stability of RNAs translated at the endoplasmic reticulum (ER). Here, we identify a localized NMD pathway dedicated to ER-translated mRNAs. We previously identified NBAS, a component of the Syntaxin 18 complex involved in Golgi-to-ER trafficking, as a novel NMD factor. Furthermore, we show that NBAS fulfills an independent function in NMD. This ER-NMD pathway requires the interaction of NBAS with the core NMD factor UPF1, which is partially localized at the ER in the proximity of the translocon. NBAS and UPF1 coregulate the stability of ER-associated transcripts, in particular those associated with the cellular stress response. We propose a model where NBAS recruits UPF1 to the membrane of the ER and activates an ER-dedicated NMD pathway, thus providing an ER-protective function by ensuring quality control of ER-translated mRNAs.

Enzymatic removal of ribonucleotides from DNA is essential for mammalian genome integrity and development.

The presence of ribonucleotides in genomic DNA is undesirable given their increased susceptibility to hydrolysis. Ribonuclease (RNase) H enzymes that recognize and process such embedded ribonucleotides are present in all domains of life. However, in unicellular organisms such as budding yeast, they are not required for viability or even efficient cellular proliferation, while in humans, RNase H2 hypomorphic mutations cause the neuroinflammatory disorder Aicardi-Goutières syndrome. Here, we report that RNase H2 is an essential enzyme in mice, required for embryonic growth from gastrulation onward. RNase H2 null embryos accumulate large numbers of single (or di-) ribonucleotides embedded in their genomic DNA (>1,000,000 per cell), resulting in genome instability and a p53-dependent DNA-damage response. Our findings establish RNase H2 as a key mammalian genome surveillance enzyme required for ribonucleotide removal and demonstrate that ribonucleotides are the most commonly occurring endogenous nucleotide base lesion in replicating cells.

DNA lesion bypass and the stochastic dynamics of transcription-coupled repair.

DNA base damage is a major source of oncogenic mutations and disruption to gene expression. The stalling of RNA polymerase II (RNAP) at sites of DNA damage and the subsequent triggering of repair processes have major roles in shaping the genome-wide distribution of mutations, clearing barriers to transcription, and minimizing the production of miscoded gene products. Despite its importance for genetic integrity, key mechanistic features of this transcription-coupled repair (TCR) process are controversial or unknown. Here, we exploited a well-powered in vivo mammalian model system to explore the mechanistic properties and parameters of TCR for alkylation damage at fine spatial resolution and with discrimination of the damaged DNA strand. For rigorous interpretation, a generalizable mathematical model of DNA damage and TCR was developed. Fitting experimental data to the model and simulation revealed that RNA polymerases frequently bypass lesions without triggering repair, indicating that small alkylation adducts are unlikely to be an efficient barrier to gene expression. Following a burst of damage, the efficiency of transcription-coupled repair gradually decays through gene bodies with implications for the occurrence and accurate inference of driver mutations in cancer. The reinitation of transcription from the repair site is not a general feature of transcription-coupled repair, and the observed data is consistent with reinitiation never taking place. Collectively, these results reveal how the directional but stochastic activity of TCR shapes the distribution of mutations following DNA damage.

Discrete measurements of RNA polymerase and reverse transcriptase fidelity reveal evolutionary tuning.

Direct methods for determining the fidelity of DNA polymerases are robust, with relatively little sample manipulation before sequencing. In contrast, methods for measuring RNA polymerase and reverse transcriptase fidelities are complicated by additional preparation steps that introduce ambiguity and error. Here, we describe a sequencing method, termed Roll-Seq, for simultaneously determining the individual fidelities of RNA polymerases and reverse transcriptases (RT) using Pacific Biosciences Single Molecule Real-Time sequencing. By employing reverse transcriptases with high rolling-circle activity, Roll-Seq generates long concatemeric cDNA from a circular RNA template. To discern the origin of a mutation, errors are recorded and determined to occur within a single concatemer (reverse transcriptase error) or all concatemers (RNA polymerase error) over the cDNA strand. We used Roll-Seq to measure the fidelities of T7 RNA polymerases, a Group II intron-encoded RT (Induro), and two LINE RTs (Fasciolopsis buski R2-RT and human LINE-1). Substitution rates for Induro and R2-RT are the same for cDNA and second strand synthesis while LINE-1 has 2.5-fold lower fidelity when performing second strand synthesis. Deletion and insertion rates increase for all RTs during second strand synthesis. In addition, we find that a structured RNA template impacts fidelity for both RNA polymerase and RT. The accuracy and precision of Roll-Seq enable this method to be applied as a complementary analysis to structural and mechanistic characterization of RNA polymerases and reverse transcriptases or as a screening method for RNAP and RT fidelity.

Pervasive lesion segregation shapes cancer genome evolution.

Cancers arise through the acquisition of oncogenic mutations and grow by clonal expansion1,2. Here we reveal that most mutagenic DNA lesions are not resolved into a mutated DNA base pair within a single cell cycle. Instead, DNA lesions segregate, unrepaired, into daughter cells for multiple cell generations, resulting in the chromosome-scale phasing of subsequent mutations. We characterize this process in mutagen-induced mouse liver tumours and show that DNA replication across persisting lesions can produce multiple alternative alleles in successive cell divisions, thereby generating both multiallelic and combinatorial genetic diversity. The phasing of lesions enables accurate measurement of strand-biased repair processes, quantification of oncogenic selection and fine mapping of sister-chromatid-exchange events. Finally, we demonstrate that lesion segregation is a unifying property of exogenous mutagens, including UV light and chemotherapy agents in human cells and tumours, which has profound implications for the evolution and adaptation of cancer genomes.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome.

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Comparative transcriptomics of primary cells in vertebrates.

Gene expression profiles in homologous tissues have been observed to be different between species, which may be due to differences between species in the gene expression program in each cell type, but may also reflect differences in cell type composition of each tissue in different species. Here, we compare expression profiles in matching primary cells in human, mouse, rat, dog, and chicken using Cap Analysis Gene Expression (CAGE) and short RNA (sRNA) sequencing data from FANTOM5. While we find that expression profiles of orthologous genes in different species are highly correlated across cell types, in each cell type many genes were differentially expressed between species. Expression of genes with products involved in transcription, RNA processing, and transcriptional regulation was more likely to be conserved, while expression of genes encoding proteins involved in intercellular communication was more likely to have diverged during evolution. Conservation of expression correlated positively with the evolutionary age of genes, suggesting that divergence in expression levels of genes critical for cell function was restricted during evolution. Motif activity analysis showed that both promoters and enhancers are activated by the same transcription factors in different species. An analysis of expression levels of mature miRNAs and of primary miRNAs identified by CAGE revealed that evolutionary old miRNAs are more likely to have conserved expression patterns than young miRNAs. We conclude that key aspects of the regulatory network are conserved, while differential expression of genes involved in cell-to-cell communication may contribute greatly to phenotypic differences between species.

Molecular Basis of Extramural Vascular Invasion (EMVI) in Colorectal Carcinoma.

BACKGROUND: Extramural vascular invasion (EMVI) is a known poor prognostic factor in colorectal carcinoma; however, its molecular basis has not been defined. This study aimed to assess the expression of molecular markers in EMVI positive colorectal carcinoma to understand their tumor microenvironment. METHODS: Immunohistochemistry was performed on tissue microarrays of surgically resected colorectal cancer specimens for immunological markers, and BRAFV600E mutation (and on the tissue blocks for mismatch repair proteins). Automated quantification was used for CD8, LAG3, FOXP3, PU1, and CD163, and manual quantification was used for PDL1, HLA I markers (beta-2 microglobulin, HC10), and HLA II. The Wilcoxon rank-sum test was used to compare EMVI positive and negative tumors. A logistic regression model was fitted to assess the predictive effect of biomarkers on EMVI. RESULTS: There were 340 EMVI positive and 678 EMVI negative chemo naïve tumors. PDL1 was barely expressed on tumor cells (median 0) in the entire cohort. We found a significantly lower expression of CD8, LAG3, FOXP3, PU1 cells, PDL1 positive macrophages, and beta-2 microglobulin on tumor cells in the EMVI positive subset (p ≤ 0.001). There was no association of BRAFV600E or deficient mismatch repair proteins (dMMR) with EMVI. PU1 (OR 0.8, 0.7-0.9) and low PDL1 (OR 1.6, 1.1-2.3) independently predicted EMVI on multivariate logistic regression among all biomarkers examined. CONCLUSION: There is a generalized blunting of immune response in EMVI positive colorectal carcinoma, which may contribute to a worse prognosis. Tumor-associated macrophages seem to play the most significant role in determining EMVI.

Contribution of clonal hematopoiesis to adult-onset hemophagocytic lymphohistiocytosis.

Adult-onset hemophagocytic lymphohistiocytosis (HLH) is a rare, life-threatening disease of immune hyperactivation. Unlike pediatric HLH, adult HLH is rarely driven by germline genetic variants. Although numerous precipitating etiologies have been identified, the reason that HLH occurs in only a subset of individuals and how other factors contribute to the disease remains unknown. We hypothesized that clonal hematopoiesis (CH), a state in which somatic mutations in blood cells cause an expanded population of mutant hematopoietic cells and drive an aberrant inflammatory state, could contribute to adult-onset HLH. In a highly annotated cohort of older adults with HLH we found that CH was more prevalent than in control cohorts. Using the adult-onset HLH mouse model in which repeated treatments of the TLR9 agonist, ODN1826, was delivered to the mouse, we observed that macrophages carrying mutations in Tet2, one of the most commonly mutated genes in CH, have an enhanced inflammatory response to TLR9 agonism. Finally, mice carrying Tet2 mutations in the hematopoietic compartment (a common model for CH) displayed an exaggerated response to TLR9 agonism, including worse splenomegaly and anemia. Our data suggest that CH is more common in individuals with adult-onset HLH and can contribute to the pathophysiology of this disease.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions.

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Thoracic nuclear protein in testis (NUT) carcinoma: expanded pathological spectrum with expression of thyroid transcription factor-1 and neuroendocrine markers.

AIMS: Nuclear protein in testis (NUT) carcinoma, an aggressive tumour driven by NUTM1 rearrangements, often involves the lung/mediastinum and shows squamous differentiation. We encountered an index patient with a thoracic NUT carcinoma diagnosed by molecular testing, showing extensive pleural involvement and diffuse thyroid transcription factor-1 (TTF-1) expression, initially suggestive of lung adenocarcinoma with pseudomesotheliomatous growth. We thus gathered an institutional series of thoracic NUT carcinomas to examine their pathological spectrum. METHODS AND RESULTS: We searched for thoracic NUT carcinomas in our surgical pathology files and in 2289 consecutive patients with primary thoracic tumours investigated with RNA-based assays. We performed NUT immunohistochemistry on 425 additional lung adenocarcinomas. Collectively, we identified six patients (five men and one woman; age 31-80 years; four never-smokers) with thoracic NUT carcinomas confirmed by molecular testing (including five with positive NUT immunohistochemistry). They died at 2.3-12.9 months (median, 2.8 months) after presentation. Two patients were diagnosed by histopathological assessment, and the remaining four (including the index patient) were diagnosed by molecular testing. Analysis of the index case revealed expression of multiple neuroendocrine markers and TTF-1; no ultrastructural evidence of neuroendocrine differentiation was noted. No additional NUT-positive cases were found by immunohistochemical screening. CONCLUSIONS: Although NUT carcinoma classically shows squamous differentiation, it can rarely express TTF-1 (even diffusely) and/or multiple neuroendocrine markers. This immunophenotypic spectrum may lead to diagnostic confusion with pulmonary adenocarcinoma, neuroendocrine tumour, and others. To circumvent this pitfall, NUT immunohistochemistry and/or NUTM1 molecular testing should be considered in primitive-appearing tumours, regardless of their immunophenotypic features.

Lagging-strand replication shapes the mutational landscape of the genome.

The origin of mutations is central to understanding evolution and of key relevance to health. Variation occurs non-randomly across the genome, and mechanisms for this remain to be defined. Here we report that the 5' ends of Okazaki fragments have significantly increased levels of nucleotide substitution, indicating a replicative origin for such mutations. Using a novel method, emRiboSeq, we map the genome-wide contribution of polymerases, and show that despite Okazaki fragment processing, DNA synthesized by error-prone polymerase-α (Pol-α) is retained in vivo, comprising approximately 1.5% of the mature genome. We propose that DNA-binding proteins that rapidly re-associate post-replication act as partial barriers to Pol-δ-mediated displacement of Pol-α-synthesized DNA, resulting in incorporation of such Pol-α tracts and increased mutation rates at specific sites. We observe a mutational cost to chromatin and regulatory protein binding, resulting in mutation hotspots at regulatory elements, with signatures of this process detectable in both yeast and humans.

Timing But Not Patterns of Recurrence Is Different Between Node-negative and Node-positive Resected Pancreatic Cancer.

OBJECTIVE: Our aim was to evaluate recurrence patterns of surgically resected PDAC patients with negative (pN0) or positive (pN1) lymph nodes. SUMMARY BACKGROUND DATA: Pancreatic ductal adenocarcinoma (PDAC) is predicted to become the second leading cause of cancer death by 2030. This is mostly due to early local and distant metastasis, even after surgical resection. Knowledge about patterns of recurrence in different patient populations could offer new therapeutic avenues. METHODS: Clinicopathologic data were collected for 546 patients who underwent resection of their PDAC between 2005 and 2016 from 2 tertiary university centers. Patients were divided into an upfront resection group (n = 394) and a neoadjuvant group (n = 152). RESULTS: Tumor recurrence was significantly less common in pN0 patients as compared with pN1 patients, (upfront surgery: 55% vs. 77%, P < 0.001 and 64% vs. 78%, P = 0.040 in the neoadjuvant group). In addition, time to recurrence was significantly longer in pN0 versus pN1 patients in the upfront resected patients (median 16 mo pN0 vs. 10 mo pN1 P < 0.001), and the neoadjuvant group (pN0 21 mo vs. 11 mo pN1, P < 0.001). Of the patients who recurred, 62% presented with distant metastases (63% of pN0 and 62% of pN1, P = 0.553), 24% with local disease (27% of pN0 and 23% of pN1, P = 0.672) and 14% with synchronous local and distant disease (10% of pN0 and 15% of pN1, P = 0.292). Similarly, there was no difference in recurrence patterns between pN0 and pN1 in the neoadjuvant group, in which 68% recurred with distant metastases (76% of pN0 and 64% of pN1, P = 0.326) and 18% recurred with local disease (pN0: 22% and pN1: 15%, P = 0.435). CONCLUSION: Time to recurrence was significantly longer for pN0 patients. However, patterns of recurrence for pN0 vs. pN1 patients were identical. Lymph node status was predictive of time to recurrence, but not location of recurrence.