Comparison of tumor-informed and tumor-naïve sequencing assays for ctDNA detection in breast cancer.

Analysis of circulating tumor DNA (ctDNA) to monitor cancer dynamics and detect minimal residual disease has been an area of increasing interest. Multiple methods have been proposed but few studies have compared the performance of different approaches. Here, we compare detection of ctDNA in serial plasma samples from patients with breast cancer using different tumor-informed and tumor-naïve assays designed to detect structural variants (SVs), single nucleotide variants (SNVs), and/or somatic copy-number aberrations, by multiplex PCR, hybrid capture, and different depths of whole-genome sequencing. Our results demonstrate that the ctDNA dynamics and allele fractions (AFs) were highly concordant when analyzing the same patient samples using different assays. Tumor-informed assays showed the highest sensitivity for detection of ctDNA at low concentrations. Hybrid capture sequencing targeting between 1,347 and 7,491 tumor-identified mutations at high depth was the most sensitive assay, detecting ctDNA down to an AF of 0.00024% (2.4 parts per million, ppm). Multiplex PCR targeting 21-47 tumor-identified SVs per patient detected ctDNA down to 0.00047% AF (4.7 ppm) and has potential as a clinical assay.

Metformin suppresses pro-inflammatory cytokines in vitreous of diabetes patients and human retinal vascular endothelium.

Metformin is a traditional anti-hyperglycemic medication that has recently been shown to benefit vascular complications of diabetes via an anti-inflammatory mechanism other than glycemic control. This study aims to test the hypothesis that metformin suppresses diabetic retinopathy (DR) associated intraocular inflammation. Human vitreous from control and proliferative diabetic retinopathy (PDR) patients with or without long-term metformin treatment (> 5 years) were collected for multiple inflammatory cytokines measurements with a cytokine array kit. The vast majority of the measurable cytokines in PDR vitreous has a lower level in metformin group than non-metformin group. Although the p values are not significant due to a relatively small sample size and large deviations, the 95% confidence interval (CI) for the mean difference between the two groups shows some difference in the true values should not be neglected. Using quantitative ELISA, soluble intercellular adhesion molecule -1 (ICAM-1) and monocyte chemoattractant protein -1 (MCP-1) presented with significantly lower concentrations in metformin group versus non-metformin group. Metformin group also has significantly less up-regulated cytokines and diminished positive correlations among the cytokines when compared to non-metformin group. Possible role of AMP-activated protein kinase (AMPK) and nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) in metformin's anti-inflammatory effects were studied in human retinal vascular endothelial cells (hRVECs) cultured in normal glucose (NG) and high glucose (HG) conditions. Metformin inhibited HG-induced ICAM-1, IL-8, and MCP-1 via AMPK activation, whereas pharmacological AMPK inhibition had no effect on its inhibition of NF-κB p65, sICAM-1, and tumor necrosis factor-α (TNF-α). Metformin-induced suppression of the inflammatory cytokines could also be mediated through its direct inhibition of NF-κB, independent of AMPK pathway. This is a proof-of-concept study that found metformin treatment was associated with reduced inflammatory responses in vitreous of diabetes patients and retinal vascular endothelial cells, supporting the rationale for using metformin to treat DR at an early stage.

Rearrangement processes and structural variations show evidence of selection in oesophageal adenocarcinomas.

Oesophageal adenocarcinoma (OAC) provides an ideal case study to characterize large-scale rearrangements. Using whole genome short-read sequencing of 383 cases, for which 214 had matched whole transcriptomes, we observed structural variations (SV) with a predominance of deletions, tandem duplications and inter-chromosome junctions that could be identified as LINE-1 mobile element (ME) insertions. Complex clusters of rearrangements resembling breakage-fusion-bridge cycles or extrachromosomal circular DNA accounted for 22% of complex SVs affecting known oncogenes. Counting SV events affecting known driver genes substantially increased the recurrence rates of these drivers. After excluding fragile sites, we identified 51 candidate new drivers in genomic regions disrupted by SVs, including ETV5, KAT6B and CLTC. RUNX1 was the most recurrently altered gene (24%), with many deletions inactivating the RUNT domain but preserved the reading frame, suggesting an altered protein product. These findings underscore the importance of identification of SV events in OAC with implications for targeted therapies.

NRG1 fusions in breast cancer.

BACKGROUND: NRG1 gene fusions may be clinically actionable, since cancers carrying the fusion transcripts can be sensitive to tyrosine kinase inhibitors. The NRG1 gene encodes ligands for the HER2(ERBB2)-ERBB3 heterodimeric receptor tyrosine kinase, and the gene fusions are thought to lead to autocrine stimulation of the receptor. The NRG1 fusion expressed in the breast cancer cell line MDA-MB-175 serves as a model example of such fusions, showing the proposed autocrine loop and exceptional drug sensitivity. However, its structure has not been properly characterised, its oncogenic activity has not been fully explained, and there is limited data on such fusions in breast cancer. METHODS: We analysed genomic rearrangements and transcripts of NRG1 in MDA-MB-175 and a panel of 571 breast cancers. RESULTS: We found that the MDA-MB-175 fusion-originally reported as a DOC4(TENM4)-NRG1 fusion, lacking the cytoplasmic tail of NRG1-is in reality a double fusion, PPP6R3-TENM4-NRG1, producing multiple transcripts, some of which include the cytoplasmic tail. We hypothesise that many NRG1 fusions may be oncogenic not for lacking the cytoplasmic domain but because they do not encode NRG1's nuclear-localised form. The fusion in MDA-MB-175 is the result of a very complex genomic rearrangement, which we partially characterised, that creates additional expressed gene fusions, RSF1-TENM4, TPCN2-RSF1, and MRPL48-GAB2. We searched for NRG1 rearrangements in 571 breast cancers subjected to genome sequencing and transcriptome sequencing and found four cases (0.7%) with fusions, WRN-NRG1, FAM91A1-NRG1, ARHGEF39-NRG1, and ZNF704-NRG1, all splicing into NRG1 at the same exon as in MDA-MB-175. However, the WRN-NRG1 and ARHGEF39-NRG1 fusions were out of frame. We identified rearrangements of NRG1 in many more (8% of) cases that seemed more likely to inactivate than to create activating fusions, or whose outcome could not be predicted because they were complex, or both. This is not surprising because NRG1 can be pro-apoptotic and is inactivated in some breast cancers. CONCLUSIONS: Our results highlight the complexity of rearrangements of NRG1 in breast cancers and confirm that some do not activate but inactivate. Careful interpretation of NRG1 rearrangements will therefore be necessary for appropriate patient management.

Single-cell sequencing of genomic DNA resolves sub-clonal heterogeneity in a melanoma cell line.

We performed shallow single-cell sequencing of genomic DNA across 1475 cells from a cell-line, COLO829, to resolve overall complexity and clonality. This melanoma tumor-line has been previously characterized by multiple technologies and is a benchmark for evaluating somatic alterations. In some of these studies, COLO829 has shown conflicting and/or indeterminate copy number and, thus, single-cell sequencing provides a tool for gaining insight. Following shallow single-cell sequencing, we first identified at least four major sub-clones by discriminant analysis of principal components of single-cell copy number data. Based on clustering, break-point and loss of heterozygosity analysis of aggregated data from sub-clones, we identified distinct hallmark events that were validated within bulk sequencing and spectral karyotyping. In summary, COLO829 exhibits a classical Dutrillaux's monosomic/trisomic pattern of karyotype evolution with endoreduplication, where consistent sub-clones emerge from the loss/gain of abnormal chromosomes. Overall, our results demonstrate how shallow copy number profiling can uncover hidden biological insights.

Pan-cancer analysis of whole genomes identifies driver rearrangements promoted by LINE-1 retrotransposition.

About half of all cancers have somatic integrations of retrotransposons. Here, to characterize their role in oncogenesis, we analyzed the patterns and mechanisms of somatic retrotransposition in 2,954 cancer genomes from 38 histological cancer subtypes within the framework of the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. We identified 19,166 somatically acquired retrotransposition events, which affected 35% of samples and spanned a range of event types. Long interspersed nuclear element (LINE-1; L1 hereafter) insertions emerged as the first most frequent type of somatic structural variation in esophageal adenocarcinoma, and the second most frequent in head-and-neck and colorectal cancers. Aberrant L1 integrations can delete megabase-scale regions of a chromosome, which sometimes leads to the removal of tumor-suppressor genes, and can induce complex translocations and large-scale duplications. Somatic retrotranspositions can also initiate breakage-fusion-bridge cycles, leading to high-level amplification of oncogenes. These observations illuminate a relevant role of L1 retrotransposition in remodeling the cancer genome, with potential implications for the development of human tumors.

Fbxl17 is rearranged in breast cancer and loss of its activity leads to increased global O-GlcNAcylation.

In cancer, many genes are mutated by genome rearrangement, but our understanding of the functional consequences of this remains rudimentary. Here we report the F-box protein encoded by FBXL17 is disrupted in the region of the gene that encodes its substrate-binding leucine rich repeat (LRR) domain. Truncating Fbxl17 LRRs impaired its association with the other SCF holoenzyme subunits Skp1, Cul1 and Rbx1, and decreased ubiquitination activity. Loss of the LRRs also differentially affected Fbxl17 binding to its targets. Thus, genomic rearrangements in FBXL17 are likely to disrupt SCFFbxl17-regulated networks in cancer cells. To investigate the functional effect of these rearrangements, we performed a yeast two-hybrid screen to identify Fbxl17-interacting proteins. Among the 37 binding partners Uap1, an enzyme involved in O-GlcNAcylation of proteins was identified most frequently. We demonstrate that Fbxl17 binds to UAP1 directly and inhibits its phosphorylation, which we propose regulates UAP1 activity. Knockdown of Fbxl17 expression elevated O-GlcNAcylation in breast cancer cells, arguing for a functional role for Fbxl17 in this metabolic pathway.

The route and rate of thiamethoxam soil degradation in laboratory and outdoor incubated tests, and field studies following seed treatments or spray application.

BACKGROUND: The route and rate of degradation of thiamethoxam in the laboratory and field was investigated. The effect of dark incubation versus light/dark cycles, seed treatment versus spray, and watering-in for spray application was explored in side-by-side trials. RESULTS: Geometric mean DT50 values were 75.4 days in OECD307 studies, and 18.3 (spray) and 16.5 (seed treatment) days in the field. In laboratory soil core studies DT50 values were 24.9 to 43.5 days, with the lowest value from the light/dark incubated soil core. Mean clothianidin formation was 19.7% applied thiamethoxam [mol/mol] in OECD307 studies and 17.5 (spray) and 3.4% (seed) in field trials. CONCLUSION: Soil DT50 values decreased with increasingly realistic tests (laboratory OECD307 to soil cores to soil cores with a light/dark cycle to field trials). The majority of the differences were associated with the soil treatment in OECD307 studies which destroys soil structure and retards the degradation rate; and from the impact on soil pore water movement in light/dark conditions. Degradation rates in the field were comparable between spray application and seed treatments. Maximum clothianidin concentrations were four-fold lower for seed treatments than for spray application in field studies. © 2018 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Complex multi-enhancer contacts captured by genome architecture mapping.

The organization of the genome in the nucleus and the interactions of genes with their regulatory elements are key features of transcriptional control and their disruption can cause disease. Here we report a genome-wide method, genome architecture mapping (GAM), for measuring chromatin contacts and other features of three-dimensional chromatin topology on the basis of sequencing DNA from a large collection of thin nuclear sections. We apply GAM to mouse embryonic stem cells and identify enrichment for specific interactions between active genes and enhancers across very large genomic distances using a mathematical model termed SLICE (statistical inference of co-segregation). GAM also reveals an abundance of three-way contacts across the genome, especially between regions that are highly transcribed or contain super-enhancers, providing a level of insight into genome architecture that, owing to the technical limitations of current technologies, has previously remained unattainable. Furthermore, GAM highlights a role for gene-expression-specific contacts in organizing the genome in mammalian nuclei.

Whole-genome sequencing of nine esophageal adenocarcinoma cell lines.

Esophageal adenocarcinoma (EAC) is highly mutated and molecularly heterogeneous. The number of cell lines available for study is limited and their genome has been only partially characterized. The availability of an accurate annotation of their mutational landscape is crucial for accurate experimental design and correct interpretation of genotype-phenotype findings. We performed high coverage, paired end whole genome sequencing on eight EAC cell lines-ESO26, ESO51, FLO-1, JH-EsoAd1, OACM5.1 C, OACP4 C, OE33, SK-GT-4-all verified against original patient material, and one esophageal high grade dysplasia cell line, CP-D. We have made available the aligned sequence data and report single nucleotide variants (SNVs), small insertions and deletions (indels), and copy number alterations, identified by comparison with the human reference genome and known single nucleotide polymorphisms (SNPs). We compare these putative mutations to mutations found in primary tissue EAC samples, to inform the use of these cell lines as a model of EAC.

Mutational signatures in esophageal adenocarcinoma define etiologically distinct subgroups with therapeutic relevance.

Esophageal adenocarcinoma (EAC) has a poor outcome, and targeted therapy trials have thus far been disappointing owing to a lack of robust stratification methods. Whole-genome sequencing (WGS) analysis of 129 cases demonstrated that this is a heterogeneous cancer dominated by copy number alterations with frequent large-scale rearrangements. Co-amplification of receptor tyrosine kinases (RTKs) and/or downstream mitogenic activation is almost ubiquitous; thus tailored combination RTK inhibitor (RTKi) therapy might be required, as we demonstrate in vitro. However, mutational signatures showed three distinct molecular subtypes with potential therapeutic relevance, which we verified in an independent cohort (n = 87): (i) enrichment for BRCA signature with prevalent defects in the homologous recombination pathway; (ii) dominant T>G mutational pattern associated with a high mutational load and neoantigen burden; and (iii) C>A/T mutational pattern with evidence of an aging imprint. These subtypes could be ascertained using a clinically applicable sequencing strategy (low coverage) as a basis for therapy selection.

Mobile element insertions are frequent in oesophageal adenocarcinomas and can mislead paired-end sequencing analysis.

BACKGROUND: Mobile elements are active in the human genome, both in the germline and cancers, where they can mutate driver genes. RESULTS: While analysing whole genome paired-end sequencing of oesophageal adenocarcinomas to find genomic rearrangements, we identified three ways in which new mobile element insertions appear in the data, resembling translocation or insertion junctions: inserts where unique sequence has been transduced by an L1 (Long interspersed element 1) mobile element; novel inserts that are confidently, but often incorrectly, mapped by alignment software to L1s or polyA tracts in the reference sequence; and a combination of these two ways, where different sequences within one insert are mapped to different loci. We identified nine unique sequences that were transduced by neighbouring L1s, both L1s in the reference genome and L1s not present in the reference. Many of the resulting inserts were small fragments that include little or no recognisable mobile element sequence. We found 6 loci in the reference genome to which sequence reads from inserts were frequently mapped, probably erroneously, by alignment software: these were either L1 sequence or particularly long polyA runs. Inserts identified from such apparent rearrangement junctions averaged 16 inserts/tumour, range 0-153 insertions in 43 tumours. However, many inserts would not be detected by mapping the sequences to the reference genome, because they do not include sufficient mappable sequence. To estimate total somatic inserts we searched for polyA sequences that were not present in the matched normal or other normals from the same tumour batch, and were not associated with known polymorphisms. Samples of these candidate inserts were verified by sequencing across them or manual inspection of surrounding reads: at least 85 % were somatic and resembled L1-mediated events, most including L1Hs sequence. Approximately 100 such inserts were detected per tumour on average (range zero to approximately 700). CONCLUSIONS: Somatic mobile elements insertions are abundant in these tumours, with over 75 % of cases having a number of novel inserts detected. The inserts create a variety of problems for the interpretation of paired-end sequencing data.

Suppression of choroidal neovascularization through inhibition of APE1/Ref-1 redox activity.

PURPOSE: The redox function of APE1/Ref-1 is a key regulator in pathological angiogenesis, such as retinal neovascularization and tumor growth. In this study, we examined whether inhibition of APE1/Ref-1 redox function by a small molecule inhibitor E3330 suppresses experimental choroidal neovascularization (CNV) in vitro and in vivo. METHODS: Primate choroid endothelial cells (CECs) received treatment of 0 to 100 µM E3330 alone or cotreatment of E3330 and 500 µg/mL anti-VEGF antibody bevacizumab. Choroid endothelial cell angiogenic function was examined by cell proliferation, migration, and tube formation assays. The effects of E3330 on NF-κB and STAT3 signaling pathways were determined by reporter gene assay, Western blot, and ELISA. Laser-induced CNV mouse model was used to test the effects of E3330 in vivo. Potential toxicity of E3330 was evaluated by TUNEL assay. RESULTS: The E3330 of 25 to 100 µM dose-dependently suppressed CEC proliferation, migration, and tube formation, in the absence of noticeable cell toxicity. Lower doses of E3330 (10-20 µM) reduced the transcriptional activity of NF-κB and STAT3 without affecting protein phosphorylation of both molecules. At the same time, E3330 downregulated MCP-1 production in CECs. The antiangiogenic effect of E3330 was comparable and additive to bevacizumab. The E3330 effectively attenuated the progression of laser-induced CNV in mice after a single intravitreal injection. CONCLUSIONS: The APE1/Ref-1 redox function regulates multiple transcription factors and inflammatory molecules, and is essential for CEC angiogenesis. Specific inhibition of APE1/Ref-1 redox function with E3330 may represent a promising novel treatment for wet AMD.

Ordering of mutations in preinvasive disease stages of esophageal carcinogenesis.

Cancer genome sequencing studies have identified numerous driver genes, but the relative timing of mutations in carcinogenesis remains unclear. The gradual progression from premalignant Barrett's esophagus to esophageal adenocarcinoma (EAC) provides an ideal model to study the ordering of somatic mutations. We identified recurrently mutated genes and assessed clonal structure using whole-genome sequencing and amplicon resequencing of 112 EACs. We next screened a cohort of 109 biopsies from 2 key transition points in the development of malignancy: benign metaplastic never-dysplastic Barrett's esophagus (NDBE; n=66) and high-grade dysplasia (HGD; n=43). Unexpectedly, the majority of recurrently mutated genes in EAC were also mutated in NDBE. Only TP53 and SMAD4 mutations occurred in a stage-specific manner, confined to HGD and EAC, respectively. Finally, we applied this knowledge to identify high-risk Barrett's esophagus in a new non-endoscopic test. In conclusion, mutations in EAC driver genes generally occur exceptionally early in disease development with profound implications for diagnostic and therapeutic strategies.

The relative timing of mutations in a breast cancer genome.

Many tumors have highly rearranged genomes, but a major unknown is the relative importance and timing of genome rearrangements compared to sequence-level mutation. Chromosome instability might arise early, be a late event contributing little to cancer development, or happen as a single catastrophic event. Another unknown is which of the point mutations and rearrangements are selected. To address these questions we show, using the breast cancer cell line HCC1187 as a model, that we can reconstruct the likely history of a breast cancer genome. We assembled probably the most complete map to date of a cancer genome, by combining molecular cytogenetic analysis with sequence data. In particular, we assigned most sequence-level mutations to individual chromosomes by sequencing of flow sorted chromosomes. The parent of origin of each chromosome was assigned from SNP arrays. We were then able to classify most of the mutations as earlier or later according to whether they occurred before or after a landmark event in the evolution of the genome, endoreduplication (duplication of its entire genome). Genome rearrangements and sequence-level mutations were fairly evenly divided earlier and later, suggesting that genetic instability was relatively constant throughout the life of this tumor, and chromosome instability was not a late event. Mutations that caused chromosome instability would be in the earlier set. Strikingly, the great majority of inactivating mutations and in-frame gene fusions happened earlier. The non-random timing of some of the mutations may be evidence that they were selected.

Structural analysis of the genome of breast cancer cell line ZR-75-30 identifies twelve expressed fusion genes.

BACKGROUND: It has recently emerged that common epithelial cancers such as breast cancers have fusion genes like those in leukaemias. In a representative breast cancer cell line, ZR-75-30, we searched for fusion genes, by analysing genome rearrangements. RESULTS: We first analysed rearrangements of the ZR-75-30 genome, to around 10kb resolution, by molecular cytogenetic approaches, combining array painting and array CGH. We then compared this map with genomic junctions determined by paired-end sequencing. Most of the breakpoints found by array painting and array CGH were identified in the paired end sequencing-55% of the unamplified breakpoints and 97% of the amplified breakpoints (as these are represented by more sequence reads). From this analysis we identified 9 expressed fusion genes: APPBP2-PHF20L1, BCAS3-HOXB9, COL14A1-SKAP1, TAOK1-PCGF2, TIAM1-NRIP1, TIMM23-ARHGAP32, TRPS1-LASP1, USP32-CCDC49 and ZMYM4-OPRD1. We also determined the genomic junctions of a further three expressed fusion genes that had been described by others, BCAS3-ERBB2, DDX5-DEPDC6/DEPTOR and PLEC1-ENPP2. Of this total of 12 expressed fusion genes, 9 were in the coamplification. Due to the sensitivity of the technologies used, we estimate these 12 fusion genes to be around two-thirds of the true total. Many of the fusions seem likely to be driver mutations. For example, PHF20L1, BCAS3, TAOK1, PCGF2, and TRPS1 are fused in other breast cancers. HOXB9 and PHF20L1 are members of gene families that are fused in other neoplasms. Several of the other genes are relevant to cancer-in addition to ERBB2, SKAP1 is an adaptor for Src, DEPTOR regulates the mTOR pathway and NRIP1 is an estrogen-receptor coregulator. CONCLUSIONS: This is the first structural analysis of a breast cancer genome that combines classical molecular cytogenetic approaches with sequencing. Paired-end sequencing was able to detect almost all breakpoints, where there was adequate read depth. It supports the view that gene breakage and gene fusion are important classes of mutation in breast cancer, with a typical breast cancer expressing many fusion genes.