Epigenomic Blood-Based Early Detection of Pancreatic Cancer Employing Cell-Free DNA

BACKGROUND & AIMS

Early detection of pancreatic cancer (PaC) can drastically improve survival rates. Approximately 25% of subjects with PaC have type 2 diabetes diagnosed within 3 years prior to the PaC diagnosis, suggesting that subjects with type 2 diabetes are at high risk of occult PaC. We have developed an early-detection PaC test, based on changes in 5-hydroxymethylcytosine (5hmC) signals in cell-free DNA from plasma.

METHODS

Blood was collected from 132 subjects with PaC and 528 noncancer subjects to generate epigenomic and genomic feature sets yielding a predictive PaC signal algorithm. The algorithm was validated in a blinded cohort composed of 102 subjects with PaC, 2048 noncancer subjects, and 1524 subjects with non-PaCs.

RESULTS

5hmC differential profiling and additional genomic features enabled the development of a machine learning algorithm capable of distinguishing subjects with PaC from noncancer subjects with high specificity and sensitivity. The algorithm was validated with a sensitivity for early-stage (stage I/II) PaC of 68.3% (95% confidence interval [CI], 51.9%-81.9%) and an overall specificity of 96.9% (95% CI, 96.1%-97.7%).

CONCLUSIONS

The PaC detection test showed robust early-stage detection of PaC signal in the studied cohorts with varying type 2 diabetes status. This assay merits further clinical validation for the early detection of PaC in high-risk individuals.

Early detection of pancreatic cancer using 5-hydroxymethylation profiles in plasma-derived cell-free DNA

672

Background: Pancreatic cancer is one of the deadliest cancers, with approximately 15-20% of patients who present at diagnosis with a resectable disease. The major barrier to better outcomes is the lack of early-detection molecular tools to enable timely intervention. We have developed a test that enables the detection of pancreatic cancer from a simple blood draw. The test incorporates a novel, genome-wide sequencing-based epigenomics detection method that enriches for DNA loci that undergo active de-methylation. The measurement of 5-hydroxymethylcytosine (5hmC) provides a unique and stable biomarker for the early detection of cancer including pancreatic cancer. Methods: Whole-blood was obtained from a training cohort of 660 individuals (consisting of 132 pancreatic cancers (PaCa) and 528 non-cancers) and a validation cohort of 2,150 individuals (consisting of 102 PaCa and 2,048 non-cancers). Cell-free DNA (cfDNA) was isolated from plasma from which 5hmC and whole-genome libraries were generated and sequenced. Logistic regression algorithms were employed using 5hmC feature sets combined with physical characteristics of DNA fragments to optimally partition cancer from non-cancer samples. Results: Cross validation of the training model yielded an overall sensitivity of 65.9%,(95% CI, 57.2%–73.9%), early-stage (stage I-II) sensitivity of 57.1% (95% CI, 44%–69.5%) and a specificity of 98%. The model was further validated in a separate, non-overlapping set of blinded and independently processed samples and yielded an early-stage sensitivity of 68.3% (95% CI, 51.9%–81.9%) and a specificity of 96.9% (95% CI, 96.0%–97.6%). Conclusions: Our results demonstrate that plasma-derived cfDNA 5hmC profiles enable the accurate detection of early-stage PaCa, providing a valuable non-invasive tool especially for those individuals at high risk for the disease, including individuals with genetic predisposition and newly diagnosed type 2 diabetes. A larger clinical study (NODMED - NCT05188586) is ongoing and will provide clinical validation for the detection in individuals at high risk for this deadly disease.

Pancreatic cancer detection using 5-hydroxymethylation signatures in plasma-derived cell free DNA in high-risk patients with new onset diabetes

539

Background: Pancreatic cancer (PaCa) is the third leading cause of cancer death in the United States despite a low incidence rate. It is often diagnosed when cancer has already metastasized to distant organs. Late diagnosis deprives patients of potentially curative treatments such as surgery and impacts survival rates. People with new onset diabetes (NOD) are at 6-8 fold increased risk for PaCa compared to the general population. Indeed, 0.85% of patients with NOD will be diagnosed with PaCa within 3 years. This population of PaCa patients with NOD constitute 25% of all new pancreatic cancer diagnoses. Surveillance of the NOD population for PaCa presents an opportunity to shift PaCa diagnosis to earlier stage. Methods: Whole blood was obtained from a cohort of 167 PaCa patients and 490 patients with cancers other than PaCa as well as 836 non-cancer controls with and without NOD. Plasma was processed to isolate cfDNA and 5hmC libraries were generated and sequenced. 5hmC data is used to generate models for PaCa detection using Bluestar Genomics’s technology platform. Results: To investigate whether PaCa can be detected in plasma, we interrogated plasma-derived cfDNA hydroxymethylation in PaCa patients and non-cancer controls. Models trained using 5hmC-based biomarkers from cfDNA consistently performed with a mean test sensitivity of 61.1% [95% confidence interval (CI): 35.7% to 82.7%] and a test specificity of 97.6% (CI: 93% to 99.5%) measured across 50 cross validation iterations within the training data set, which was composed of 48.3% early stage (Stages I & II) disease. The final model was trained using all of the training data, yielding 58.4% (CI: 47.5% to 68.8%) sensitivity at 98% (CI: 96.5% to 99.0%) specificity. This model was then tested on an independent test set with 22 PaCa patients (51.7% early stage, 15 of which was NOD) and 123 non-cancer control patients (53 of which were NOD) and yielded a classification performance of 59.1% (CI: 36.4% to 79.3%) sensitivity at 95.9% (CI: 90.8% to 98.7%) specificity. The model performance in the subset of patient cohort with NOD was 53.3% (CI: 26.6% to 78.7%) sensitivity at 94.3% (CI: 84.3% to 98.8%) specificity. Lastly, sensitivity observed on an independent validation set, composed of 56 PaCa and 117 ITTP samples, was 46.4% (CI: 33.0% to 60.2%) with 100% (CI: 96.8 to 100%) specificity. Conclusions: Our results demonstrate PaCa detection in plasma-derived cfDNA using 5hmC profiles. Overall, the model performed consistently between the training and independent validation datasets. A larger clinical study is under development to clinically validate the model described in this study with the goal of identifying occult PaCa within the NOD population in order to enable earlier detection and thus improve patient outcomes.

A community challenge to evaluate RNA-seq, fusion detection, and isoform quantification methods for cancer discovery

The accurate identification and quantitation of RNA isoforms present in the cancer transcriptome is key for analyses ranging from the inference of the impacts of somatic variants to pathway analysis to biomarker development and subtype discovery. The ICGC-TCGA DREAM Somatic Mutation Calling in RNA (SMC-RNA) challenge was a crowd-sourced effort to benchmark methods for RNA isoform quantification and fusion detection from bulk cancer RNA sequencing (RNA-seq) data. It concluded in 2018 with a comparison of 77 fusion detection entries and 65 isoform quantification entries on 51 synthetic tumors and 32 cell lines with spiked-in fusion constructs. We report the entries used to build this benchmark, the leaderboard results, and the experimental features associated with the accurate prediction of RNA species. This challenge required submissions to be in the form of containerized workflows, meaning each of the entries described is easily reusable through CWL and Docker containers at https://github.com/SMC-RNA-challenge. A record of this paper's transparent peer review process is included in the supplemental information.

eVIP2: Expression-based variant impact phenotyping to predict the function of gene variants

While advancements in genome sequencing have identified millions of somatic mutations in cancer, their functional impact is poorly understood. We previously developed the expression-based variant impact phenotyping (eVIP) method to use gene expression data to characterize the function of gene variants. The eVIP method uses a decision tree-based algorithm to predict the functional impact of somatic variants by comparing gene expression signatures induced by introduction of wild-type (WT) versus mutant cDNAs in cell lines. The method distinguishes between variants that are gain-of-function, loss-of-function, change-of-function, or neutral. We present eVIP2, software that allows for pathway analysis (eVIP Pathways) and usage with RNA-seq data. To demonstrate the eVIP2 software and approach, we characterized two recurrent frameshift variants in RNF43, a negative regulator of Wnt signaling, frequently mutated in colorectal, gastric, and endometrial cancer. RNF43 WT, RNF43 R117fs, RNF43 G659fs, or GFP control cDNA were overexpressed in HEK293T cells. Analysis with eVIP2 predicted that the frameshift at position 117 was a loss-of-function mutation, as expected. The second frameshift at position 659 has been previously described as a passenger mutation that maintains the RNF43 WT function as a negative regulator of Wnt. Surprisingly, eVIP2 predicted G659fs to be a change-of-function mutation. Additional eVIP Pathways analysis of RNF43 G659fs predicted 10 pathways to be significantly altered, including TNF-α via NFκB signaling, KRAS signaling, and hypoxia, highlighting the benefit of a more comprehensive approach when determining the impact of gene variant function. To validate these predictions, we performed reporter assays and found that each pathway activated by expression of RNF43 G659fs, but not expression of RNF43 WT, was identified as impacted by eVIP2, supporting that RNF43 G659fs is a change-of-function mutation and its effect on the identified pathways. Pathway activation was further validated by Western blot analysis. Lastly, we show primary colon adenocarcinoma patient samples with R117fs and G659fs variants have transcriptional profiles similar to BRAF missense mutations with activated RAS/MAPK signaling, consistent with KRAS signaling pathways being GOF in both variants. The eVIP2 method is an important step towards overcoming the current challenge of variant interpretation in the implementation of precision medicine. eVIP2 is available at https://github.com/BrooksLabUCSC/eVIP2.

Pancreatic cancer detection using EpiDetect signatures in plasma-derived cell free DNA in high-risk patients with new onset diabetes

e16265

Background: Pancreatic cancer (PaCa) is the third leading cause of cancer death in the United States despite its low incidence rate, owing to a 5-year survival rate of 10%. It is often asymptomatic in early stage, resulting in the majority of diagnoses occurring when cancer has already metastasized to distant organs. Late diagnosis deprives patients of potentially curative treatments such as surgery and impacts survival rates. Diabetes can be an early symptom of PaCa. Indeed, 25% of PaCa patients had a preceding diabetes diagnosis. Among all people with new onset diabetes (NOD), 0.85% will be diagnosed with PaCa within 3 years, which represents 6-8 fold increased risk for PaCa compared to the general population. Surveillance of the NOD population for PaCa presents an opportunity to shift PaCa diagnosis to earlier stage by finding it sooner. Methods: Whole blood was obtained from a cohort of 117 PaCa patients as well as 800 non-cancer controls with and without NOD. Plasma was processed to isolate cfDNA and 5hmC and low pass whole genome libraries were generated and sequenced. The EpiDetect assay combines 5hmC and whole genome sequencing data and were generated using Bluestar Genomics’s technology platform. Results: To investigate whether PaCa can be detected in plasma, we interrogated plasma-derived cfDNA epigenomic and genomic signal from PaCa patients and non-cancer controls. We first trained stacked ensemble models on PaCa and non-cancer samples utilizing 5hmC, fragmentation and CNV-based biomarkers from cfDNA. These models performed stably with a median of 72.8% sensitivity and 90.1% specificity measured across 25 outer fold iterations using the training data set, which was composed of 50% early stage (Stages I & II) disease. The final binomial ensemble model was trained using all of the training data, yielding an area under the receiver operating characteristic curve (auROC) of 0.9, with 75% sensitivity and 89% specificity. This model was then tested on an independent validation data set from 33 PaCa patients (24 with diabetes, 15 of which was NOD) and 202 non-cancer control patients (76 with diabetes, 51 of which was NOD) and yielded a classification performance auROC of 0.9 with 67% sensitivity at 92% specificity. Lastly, model performance in the subset of patient cohort with NOD only had an auROC of 0.87 with 60% sensitivity at 88% specificity. Conclusions: Our results indicate that 5hmC profiles along with CNV and fragmentation patterns from cfDNA can be used to detect PaCa in plasma-derived cfDNA. Overall, model performance was stable and consistent between the training and independent validation datasets. A larger clinical study is under development to investigate the utility of the model described in this pilot study in identifying occult PaCa within the NOD population, with the aim of shifting diagnosis to early stage and potentially improving patient outcomes.

Genome-wide 5hmC profiles to enable cancer detection and tissue of origin classification in breast, colorectal, lung, ovarian, and pancreatic cancers

3044

Background: Epigenomics assays have recently become popular tools for identification of molecular biomarkers, both in tissue and in plasma. In particular 5-hydroxymethyl-cytosine (5hmC) method, has been shown to enable the epigenomic regulation of gene expression and subsequent gene activity, with different patterns, across several tumor and normal tissues types. In this study we show that 5hmC profiles enable discrete classification of tumor and normal tissue for breast, colorectal, lung ovary and pancreas. Such classification was also recapitulated in cfDNA from patient with breast, colorectal, lung, ovarian and pancreatic cancers. Methods: DNA was isolated from 176 fresh frozen tissues from breast, colorectal, lung, ovary and pancreas (44 per tumor per tissue type and up to 11 tumor tissues for each stage (I-IV)) and up to 10 normal tissues per tissue type. cfDNA was isolated from plasma from 783 non-cancer individuals and 569 cancer patients. Plasma-isolated cfDNA and tumor genomic DNA, were enriched for the 5hmC fraction using chemical labelling, sequenced, and aligned to a reference genome to construct features sets of 5hmC patterns. Results: 5hmC multinomial logistic regression analysis was employed across tumor and normal tissues and identified a set of specific and discrete tumor and normal tissue gene-based features. This indicates that we can classify samples regardless of source, with a high degree of accuracy, based on tissue of origin and also distinguish between normal and tumor status.Next, we employed a stacked ensemble machine learning algorithm combining multiple logistic regression models across diverse feature sets to the cfDNA dataset composed of 783 non cancers and 569 cancers comprising 67 breast, 118 colorectal, 210 Lung, 71 ovarian and 100 pancreatic cancers. We identified a genomic signature that enable the classification of non-cancer versus cancers with an outer fold cross validation sensitivity of 49% (CI 45%-53%) at 99% specificity. Further, individual cancer outer fold cross validation sensitivity at 99% specificity, was measured as follows: breast 30% (CI 119% -42%); colorectal 41% (CI 32%-50%); lung 49% (CI 42%-56%); ovarian 72% (CI 60-82%); pancreatic 56% (CI 46%-66%). Conclusions: This study demonstrates that 5hmC profiles can distinguish cancer and normal tissues based on their origin. Further, 5hmC changes in cfDNA enables detection of the several cancer types: breast, colorectal, lung, ovarian and pancreatic cancers. Our technology provides a non-invasive tool for cancer detection with low risk sample collection enabling improved compliance than current screening methods. Among other utilities, we believe our technology could be applied to asymptomatic high-risk individuals thus enabling enrichment for those subjects that most need a diagnostic imaging follow up.

Whole-genome characterization of lung adenocarcinomas lacking the RTK/RAS/RAF pathway

RTK/RAS/RAF pathway alterations (RPAs) are a hallmark of lung adenocarcinoma (LUAD). In this study, we use whole-genome sequencing (WGS) of 85 cases found to be RPA(-) by previous studies from The Cancer Genome Atlas (TCGA) to characterize the minority of LUADs lacking apparent alterations in this pathway. We show that WGS analysis uncovers RPA(+) in 28 (33%) of the 85 samples. Among the remaining 57 cases, we observe focal deletions targeting the promoter or transcription start site of STK11 (n = 7) or KEAP1 (n = 3), and promoter mutations associated with the increased expression of ILF2 (n = 6). We also identify complex structural variations associated with high-level copy number amplifications. Moreover, an enrichment of focal deletions is found in TP53 mutant cases. Our results indicate that RPA(-) cases demonstrate tumor suppressor deletions and genome instability, but lack unique or recurrent genetic lesions compensating for the lack of RPAs. Larger WGS studies of RPA(-) cases are required to understand this important LUAD subset.

Abstract 783: Epigenomic detection of multiple cancers in plasma derived cell free DNA

Background: Our feasibility study employed a novel genomic detection methodology that enriches 5-hydroxymethylcytosine (5hmC) loci in cell free DNA (cfDNA) from the plasma of cancer patients using click chemistry coupled with sequencing and machine learning based classification methods. These classification methods were developed to detect the presence of disease in the plasma of cancer and control subjects. Cancer and control patient cfDNA cohorts were accrued from multiple sites consisting of 48 breast, 55 lung, 32 prostate and 2 pancreatic datasets consisting of 41 and 53 cancer subjects (Set 1 and 2). In addition, a control cohort of 260 subjects (non-cancer) was employed to match cancer patient demographics (age, sex and smoking status) in a case-control study design.

Methods: Machine learning methods, applied to each cancer case cohort individually, with a balancing non-cancer cohort, were able to classify cancer and control samples. Measures of predictive performance using 5-fold cross validation coupled with out-of-fold Area Under the Receiver Operating Characteristic Curve (AUROC) measures were employed. Gene sets selected as part of biomarker discovery were further analyzed for disease relevance using pathway analysis tools (GSEA, mSigDB).

Results: 260 controls and 229 cancers from four disease types (breast, lung, pancreas and prostate) were analyzed; more than 60% of cancer patients had early stage disease (I or II). Predictive performance, employing AUROC measures, was established for breast (0.89), lung (0.84), pancreas (set 1 - 0.95 and 2 - 0.93) and prostate (0.83). The genes defining each of these predictive models were enriched for pathways relevant to disease specific etiology, notably in the control of gene regulation in these same pathways. The breast cancer cohort consisted primarily of stage I and II patients including tumors < 2 cm and these samples exhibited a higher prediction probability score. The prostate cancer cohort consisted of both indolent and aggressive disease sample and prediction performance was equally high for both (AUROC for indolent vs aggressive was 0.81 and 0.77, respectively).

Conclusions: These findings suggest that 5hmC changes in cfDNA enable non-invasive detection of early stage breast, pancreatic, prostate, and lung cancers. Furthermore, 5hmC profiling in cfDNA may enable the prediction of clinically relevant features such as tumor size in breast adenocarcinoma or indolent disease in prostate cancer. Finally, this study identified a suite of 5hmC biomarkers that may be further validated in larger, and more diverse, patient cohorts.

Citation Format: Anna Bergamaschi, Jeremy Ku, Yuhong Ning, Francois Collin, Chris Ellison, Tierney Phillips, Erin McCarthy, Wendy Wang, Michael Antoine, David Haan, Aaron Scott, Paul Lloyd, Gulfem Guler, Alan Ashworth, Stephen Quake, Samuel Levy. Epigenomic detection of multiple cancers in plasma derived cell free DNA [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 783.

Monitoring immunotherapy response in non-small cell lung cancer patients using 5-hydroxymethylcytosine signatures in circulating cell free DNA

e21505

Background: Liquid biopsies are gaining prominence for not only cancer diagnosis but also patient monitoring. Mutational signals derived from cell-free DNA (cfDNA) show promise to assess response to cancer treatment, including immunotherapy. However, reliance of these methods on mutational data from tissue biopsies limit their applicability when a tumor biopsy is unavailable, or when mutational landscape of tumor changes under the selective pressures of cancer drug treatment. Epigenomic approaches have the potential to address these shortcomings. Methods: Blood draws were obtained from a cohort of non-small cell lung cancer (NSCLC) patients (n = 19) who went on to anti-PD1 treatment prior to therapy start and while on therapy. cfDNA was isolated from plasma and was subsequently processed to generate 5hmC genome-wide profiles. Results: We analyzed cfDNA from NSCLC patients undergoing anti-PD1 therapy to investigate whether immunotherapy response can be detected from plasma. Using a predictive model trained on lung cancer and non-cancer samples, we were able to detect changes in prediction scores in patient treated with immunotherapy that were consistent with RECIST. Patients with progressive disease (n = 3), determined by RECIST, had prediction scores that increased while they received treatment. On the other hand, majority of the patients that exhibited partial response to treatment (n = 12) had predictive scores that decreased with treatment, again consistent with RECIST. Furthermore, score changes consistent with RECIST was observed one cycle prior to the RECIST timepoint in all except one patient, where an extra blood draw after baseline was available (n = 7). Annotation of the regions that account for differential scoring identified enhancer, 5’UTR and promoter regions. Comparison of partial responders to patients with progressive disease revealed genes involved in metastasis, oncogenes and tumor suppressors that change in opposing directions between these patient groups, consistent with the underlying biology. Conclusions: Our results suggest that 5hmC profiles from cfDNA can be used to determine immunotherapy response in non-small cell lung cancer patients. Compared with mutation based liquid biopsy methods to assess response, epigenomics-based methods have the advantage of being agnostic to starting tumor mutations, and not relying on a mutational analysis from tumor biopsy. Future work will help determine applicability of this method to other kinds of therapies and cancer types.

Effect of detection of epigenomic changes in plasma-derived cell-free DNA on multicancer classification

1539

Background: Epigenomic changes in DNA methylation patterns are more precisely delineated by active demethylation events as marked by 5-hydroxymethylation (5hmC) of cytosine residues. 5hmC appears to be dynamically modulated in tumor tissues and can be employed as a cancer biomarker. Strategies which interrogate 5hmC genome-wide patterns in a liquid biopsy context may provide efficient and precise technology for early cancer screening and detection. In this study we identified genome-wide 5hmC changes in plasma based circulating free DNA (cfDNA) from breast, colorectal, lung, pancreatic and prostate cancer patients versus non-cancer individuals. Methods: cfDNA was isolated from plasma, enriched for the 5hmC fraction using novel click-chemistry protocol for labelling followed by sequencing and alignment to a reference genome to construct features sets of 5hmC patterns. Regularized classification models were constructed to classify cancer samples apart from non-cancer. Results: > 500 non-cancer individuals and > 500 cancer patients across five cancer types (breast, colorectal, lung, pancreas and prostate) were included in this study. About 60% of the cancer samples were early stage disease (I or II). The ability to classify non-cancer versus cancer patients was evaluated by 5-fold cross validation of our trained prediction models. Our models were able to classify all breast cancer with a test auROC of 0.86 while prediction model classification for ER negative samples had an auROC of 0.92. Colorectal performance auROC was 0.9; lung auROC = 0.92, pancreatic auROC = 0.97 and prostate auROC = 0.91. Overall sensitivity range, when allowing 2% false positive, was between 85% and 52%. Further using 5hmC signal in blood we were able to identify several signaling pathways specifically relevant to the biology of the cancers investigated. Conclusions: These findings further demonstrate that 5hmC changes in cfDNA enable non-invasive detection of breast, colorectal, lung pancreatic, and prostate cancers. Further, 5hmC signals enabled the identification of a suite of cancer signaling pathways differentially enriched in cancers versus non-cancers. These data suggest that dynamic changes in tumor cell methylation, detectable through 5-hydroxymethylation, are contained in the circulating blood and signal active disease biology.

Pathway and network analysis of more than 2500 whole cancer genomes

The catalog of cancer driver mutations in protein-coding genes has greatly expanded in the past decade. However, non-coding cancer driver mutations are less well-characterized and only a handful of recurrent non-coding mutations, most notably TERT promoter mutations, have been reported. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancer across 38 tumor types, we perform multi-faceted pathway and network analyses of non-coding mutations across 2583 whole cancer genomes from 27 tumor types compiled by the ICGC/TCGA PCAWG project that was motivated by the success of pathway and network analyses in prioritizing rare mutations in protein-coding genes. While few non-coding genomic elements are recurrently mutated in this cohort, we identify 93 genes harboring non-coding mutations that cluster into several modules of interacting proteins. Among these are promoter mutations associated with reduced mRNA expression in TP53, TLE4, and TCF4. We find that biological processes had variable proportions of coding and non-coding mutations, with chromatin remodeling and proliferation pathways altered primarily by coding mutations, while developmental pathways, including Wnt and Notch, altered by both coding and non-coding mutations. RNA splicing is primarily altered by non-coding mutations in this cohort, and samples containing non-coding mutations in well-known RNA splicing factors exhibit similar gene expression signatures as samples with coding mutations in these genes. These analyses contribute a new repertoire of possible cancer genes and mechanisms that are altered by non-coding mutations and offer insights into additional cancer vulnerabilities that can be investigated for potential therapeutic treatments.

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.

Using Transcriptional Signatures to Find Cancer Drivers with LURE

Cancer genome projects have produced multidimensional datasets on thousands of samples. Yet, depending on the tumor type, 5-50% of samples have no known driving event. We introduce a semi-supervised method called Learning UnRealized Events (LURE) that uses a progressive label learning framework and minimum spanning analysis to predict cancer drivers based on their altered samples sharing a gene expression signature with the samples of a known event. We demonstrate the utility of the method on the TCGA Pan-Cancer Atlas dataset for which it produced a high-confidence result relating 59 new connections to 18 known mutation events including alterations in the same gene, family, and pathway. We give examples of predicted drivers involved in TP53, telomere maintenance, and MAPK/RTK signaling pathways. LURE identifies connections between genes with no known prior relationship, some of which may offer clues for targeting specific forms of cancer. Code and Supplemental Material are available on the LURE website: https://sysbiowiki.soe.ucsc.edu/lure.

Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling

The American alligator, Alligator mississippiensis, like all crocodilians, has temperature-dependent sex determination, in which the sex of an embryo is determined by the incubation temperature of the egg during a critical period of development. The lack of genetic differences between male and female alligators leaves open the question of how the genes responsible for sex determination and differentiation are regulated. Insight into this question comes from the fact that exposing an embryo incubated at male-producing temperature to estrogen causes it to develop ovaries. Because estrogen response elements are known to regulate genes over long distances, a contiguous genome assembly is crucial for predicting and understanding their impact. We present an improved assembly of the American alligator genome, scaffolded with in vitro proximity ligation (Chicago) data. We use this assembly to scaffold two other crocodilian genomes based on synteny. We perform RNA sequencing of tissues from American alligator embryos to find genes that are differentially expressed between embryos incubated at male- versus female-producing temperature. Finally, we use the improved contiguity of our assembly along with the current model of CTCF-mediated chromatin looping to predict regions of the genome likely to contain estrogen-responsive genes. We find that these regions are significantly enriched for genes with female-biased expression in developing gonads after the critical period during which sex is determined by incubation temperature. We thus conclude that estrogen signaling is a major driver of female-biased gene expression in the post-temperature sensitive period gonads.

[Bronchospasm, apnea and heart arrest following dipyridamole perfusion scintigraphy]

The occurrence of fatal respiratory insufficiency following dipyridamole-thallium imaging is described. The patient, a 67-year-old man, had a history of chronic obstructive lung disease. Since patients with a history of chronic obstructive lung disease have an increased risk of developing bronchospasm after dipyridamole infusion, it is advised to be cautious in performing dipyridamole-thallium imaging in these patients. Dobutamine may be an acceptable alternative to dipyridamole in these patients.

Cholesterol embolism as a complication of left heart catheterisation. Report of seven cases

Cholesterol embolism after left heart catheterisation by the femoral approach was diagnosed in seven men (mean age 59.6 years) out of a total of 4587 catheterisations. Diabetes was present in four patients, systemic hypertension in three, and signs of extensive atherosclerosis in six; five patients were taking anticoagulant drugs. Acute pain in the legs or abdomen occurred in six patients, two of whom had abdominal angina; renal failure was present in six patients, cutaneous manifestations in five, and a cholesterol embolus was seen in the retina in one. Six out of six patients had an appreciable increase in the erythrocyte sedimentation rate and five out of five had eosinophilia within a week of catheterisation. Renal failure was progressive in five patients, one of whom required haemodialysis. Three patients required amputation of the toes because of gangrene. Four patients died within four and a half months of catheterisation from causes not directly related to cholesterol embolism. At necropsy cholesterol emboli were found in all four patients. Cholesterol embolism is a rare but serious complication of left heart catheterisation.