Analytical Validation of an Early Detection Pancreatic Cancer Test Using 5-Hydroxymethylation Signatures.

Early detection of pancreatic cancer has been shown to improve patient survival rates. However, effective early detection tools to detect pancreatic cancer do not currently exist. The Avantect Pancreatic Cancer Test, leveraging the 5-hydroxymethylation [5-hydroxymethylcytosine (5hmC)] signatures in cell-free DNA, was developed and analytically validated to address this unmet need. We report a comprehensive analytical validation study encompassing precision, sample stability, limit of detection, interfering substance studies, and a comparison with an alternative method. The assay performance on an independent case-control patient cohort was previously reported with a sensitivity for early-stage (stage I/II) pancreatic cancer of 68.3% (95% CI, 51.9%-81.9%) and an overall specificity of 96.9% (95% CI, 96.1%-97.7%). Precision studies showed a cancer classification of 100% concordance in biological replicates. The sample stability studies revealed stable assay performance for up to 7 days after blood collection. The limit of detection studies revealed equal results between early- and late-stage cancer samples, emphasizing strong early-stage performance characteristics. Comparisons of concordance of the Avantect assay with the enzymatic methyl sequencing (EM-Seq) method, which measures both methylation (5-methylcytosine) and 5hmC, were >95% for all samples tested. The Avantect Pancreatic Cancer Test showed strong analytical validation in multiple validation studies required for laboratory-developed test accreditation. The comparison of 5hmC versus EM-Seq further validated the 5hmC approach as a robust and reproducible assay.

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.

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.

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.

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.

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.

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.

Free energy map for the co-oligomerization of formaldehyde and ammonia.

Density functional theory calculations, including Poisson-Boltzmann implicit solvent and free energy corrections, are applied to construct a free energy map of formaldehyde and ammonia co-oligomerization in aqueous solution at pH 7. The stepwise route to forming hexamethylenetetramine (HMTA), the one clearly identified major product in a complex mixture, involves a series of addition reactions of formaldehyde and ammonia coupled with dehydration and cyclization reactions at key steps in the pathway. The free energy map also allows us to propose the possible identity of some major peaks observed by mass spectroscopy in the reaction mixture being the result of stable species not along the pathway to HMTA, in particular those formed by intramolecular condensation of hydroxyl groups to form six-membered rings with ether linkages. Our study complements a baseline free energy map previously worked out for the self-oligomerization of formaldehyde in solution at pH 7 using the same computational protocol and published in this journal (J. Phys. Chem. A 2013, 117, 12658).

Thermodynamics and kinetics of imidazole formation from glyoxal, methylamine, and formaldehyde: a computational study.

Density functional theory calculations, including Poisson-Boltzmann implicit solvent and free energy corrections, are applied to study the mechanism of experimentally observed imidazole formation from the reaction of glyoxal and methylamine in solution. Our calculations suggest that a diimine species is an important intermediate in the reaction. Under acidic conditions, we find that the diimine acts as a nucleophile in attacking the carbonyl group of either formaldehyde or glyoxal to first generate an acyclic enol intermediate, which then goes on to close the ring and form imidazoles. Our results confirm that formaldehyde and, by extension, other small aldehydes are likely to be incorporated into imidazole ions in the presence of glyoxal and primary amines in clouds and aqueous aerosol. This is a new mechanism of aerosol formation by formaldehyde, the most abundant aldehyde in the atmosphere. The amount of aerosol formed will depend on the amounts of glyoxal and amines present.

Gas-phase water-mediated equilibrium between methylglyoxal and its geminal diol.

In aqueous solution, aldehydes, and to a lesser extent ketones, hydrate to form geminal diols. We investigate the hydration of methylglyoxal (MG) in the gas phase, a process not previously considered to occur in water-restricted environments. In this study, we spectroscopically identified methylglyoxal diol (MGD) and obtained the gas-phase partial pressures of MG and MGD. These results, in conjunction with the relative humidity, were used to obtain the equilibrium constant, K(P), for the water-mediated hydration of MG in the gas phase. The Gibbs free energy for this process, DeltaG(o), obtained as a result, suggests a larger than expected gas-phase diol concentration. This may have significant implications for understanding the role of organics in atmospheric chemistry.