Whole genome sequencing of cyanobacterium Nostoc sp. CCCryo 231-06 using microfluidic single cell technology

The Nostoc sp. strain CCCryo 231-06 is a cyanobacterial strain capable of surviving under extreme conditions and thus is of great interest for the astrobiology community. The knowledge of its complete genome sequence would serve as a guide for further studies. However, a major concern has been placed on the effects of contamination on the quality of sequencing data without a reference genome. Here, we report the use of microfluidic technology combined with single cell sequencing and de novo assembly to minimize the contamination and recover the complete genome of the Nostoc strain CCCryo 231-06 with high quality. 100% of the whole genome was recovered with all contaminants removed and a strongly supported phylogenetic tree. The data reported can be useful for comparative genomics for phylogenetic and taxonomic studies. The method used in this work can be applied to studies that require high-quality assemblies of genomes of unknown microorganisms.

•

This work uses a microfluidic platform for Nostoc single cell sequencing

•

This technology provides minimal contamination in single cell sequencing

•

Complete genome of the Nostoc strain CCCryo 231-06 was recovered with high quality

Direct comparison of circulating tumor DNA sequencing assays with targeted large gene panels

Next generation sequencing (NGS) assays with large targeted gene panels can comprehensively profile cancer somatic mutations in a tumor sample. Given the rapid adoption of such assays for circulating tumor DNA (ctDNA) analysis in clinical oncology, it is essential for the community to understand their analytical performance in liquid biopsy settings. Here, we directly compared five ctDNA NGS assays, most of which having a panel of 400 or more genes, with simulated samples harboring mutations relevant to solid tumors or myeloid malignancy. Our results indicate that the detection sensitivity and reproducibility of all five assays was 90% or higher when the mutations were at 0.5% or 1.0% allele frequency, and with optimal DNA input of 30 ng or 50 ng per vendor’s protocol. The performances decreased and varied dramatically, when mutations were at a 0.1% allele frequency and/or when a lower genomic input of 10 ng DNA was used. Interestingly, one of the assays repeatedly showed higher rate of false positivity than the others across two different sample sets. Multiple intrinsic technical factors pertaining to the NGS assays were further investigated. Notable differences among the assays were seen for depth of coverage and background noise, which profoundly impacted assay performance. The results derived from this study are highly informative and provide a framework to assess and select suitable assays for specific application in cancer monitoring and potential clinical use.

Improved Characterization of Complex β-Globin Gene Cluster Structural Variants Using Long-Read Sequencing

Complex insertion-deletion (indel) events in the globin genes manifest in widely variable clinical phenotypes. Many are incompletely characterized because of a historic lack of efficient methods. A more complete assessment enables improved prediction of clinical impact, which guides emerging therapeutic choices. Current methods have limited capacity for breakpoint assignment and accurate assessment of mutation extent, especially in cases containing duplications or multiple deletions and insertions. Technology, such as long-read sequencing, holds promise for significant impact in the characterization of indel events because of read lengths that span large regions, resulting in improved resolution. Four known complex β-globin gene cluster indel types were assessed using single-molecule, real-time sequencing technology and showed high correlation with previous reports, including the Caribbean locus control deletion (g.5,305,478_5,310,336del), a large β-gene duplication containing the Hb S mutation (g.4,640,335_5,290,171dup with g.5,248,232T>A, c.20A>T; variant allele fraction, 64%), and two nested variants (double deletions with intervening inversion): the Indian ^Gγ(^Aγδβ)⁰-thalassemia (g.5,246,804-5,254,275del, g.5,254,276_5,269,600inv, and g.5,269,601_5,270,442del) and the Turkish/Macedonian (δβ)⁰ thalassemia (g.5,235,064_5,236,652del, g.5,236,653_5,244,280inv, and g.5,244,281_5,255,766del). Our data confirm long-read sequencing as an efficient and accurate method to identify these clinically significant complex events. Limitations include high-complexity sample preparation requirements, which hinder routine use in clinical laboratories. Continued improvements in sample and data workflow processes are needed to accommodate volumes in a tertiary clinical laboratory.

Intratumor Epigenetic Heterogeneity-A Panel Gene Methylation Study in Thyroid Cancer

Background

Thyroid cancer (TC) is the most common endocrine malignancy, and the incidence is increasing very fast. Surgical resection and radioactive iodine ablation are major therapeutic methods, however, around 10% of differentiated thyroid cancer and all anaplastic thyroid carcinoma (ATC) are failed. Comprehensive understanding the molecular mechanisms may provide new therapeutic strategies for thyroid cancer. Even though genetic heterogeneity is rigorously studied in various cancers, epigenetic heterogeneity in human cancer remains unclear.

Methods

A total of 405 surgical resected thyroid cancer samples were employed (three spatially isolated specimens were obtained from different regions of the same tumor). Twenty-four genes were selected for methylation screening, and frequently methylated genes in thyroid cancer were used for further validation. Methylation specific PCR (MSP) approach was employed to detect the gene promoter region methylation.

Results

Five genes (AP2, CDH1, DACT2, HIN1, and RASSF1A) are found frequently methylated (>30%) in thyroid cancer. The five genes panel is used for further epigenetic heterogeneity analysis. AP2 methylation is associated with gender (P < 0.05), DACT2 methylation is associated with age, gender and tumor size (all P < 0.05), HIN1 methylation is associated to tumor size (P < 0.05) and extra-thyroidal extension (P < 0.01). RASSF1A methylation is associated with lymph node metastasis (P < 0.01). For heterogeneity analysis, AP2 methylation heterogeneity is associated with tumor size (P < 0.01), CDH1 methylation heterogeneity is associated with lymph node metastasis (P < 0.05), DACT2 methylation heterogeneity is associated with tumor size (P < 0.01), HIN1 methylation heterogeneity is associated with tumor size and extra-thyroidal extension (all P < 0.01). The multivariable analysis suggested that the risk of lymph node metastasis is 2.5 times in CDH1 heterogeneous methylation group (OR = 2.512, 95% CI 1.135, 5.557, P = 0.023). The risk of extra-thyroidal extension is almost 3 times in HIN1 heterogeneous methylation group (OR = 2.607, 95% CI 1.138, 5.971, P = 0.023).

Conclusion

Five of twenty-four genes were found frequently methylated in human thyroid cancer. Based on 5 genes panel analysis, epigenetic heterogeneity is an universal event. Epigenetic heterogeneity is associated with cancer development and progression.

Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study

Assessing the reproducibility, accuracy and utility of massively parallel DNA sequencing platforms remains an ongoing challenge. Here the Association of Biomolecular Resource Facilities (ABRF) Next-Generation Sequencing Study benchmarks the performance of a set of sequencing instruments (HiSeq/NovaSeq/paired-end 2 × 250-bp chemistry, Ion S5/Proton, PacBio circular consensus sequencing (CCS), Oxford Nanopore Technologies PromethION/MinION, BGISEQ-500/MGISEQ-2000 and GS111) on human and bacterial reference DNA samples. Among short-read instruments, HiSeq 4000 and X10 provided the most consistent, highest genome coverage, while BGI/MGISEQ provided the lowest sequencing error rates. The long-read instrument PacBio CCS had the highest reference-based mapping rate and lowest non-mapping rate. The two long-read platforms PacBio CCS and PromethION/MinION showed the best sequence mapping in repeat-rich areas and across homopolymers. NovaSeq 6000 using 2 × 250-bp read chemistry was the most robust instrument for capturing known insertion/deletion events. This study serves as a benchmark for current genomics technologies, as well as a resource to inform experimental design and next-generation sequencing variant calling.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10^-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

Integrating genomic features for non-invasive early lung cancer detection

Radiologic screening of high-risk adults reduces lung-cancer-related mortality1,2; however, a small minority of eligible individuals undergo such screening in the United States3,4. The availability of blood-based tests could increase screening uptake. Here we introduce improvements to cancer personalized profiling by deep sequencing (CAPP-Seq)5, a method for the analysis of circulating tumour DNA (ctDNA), to better facilitate screening applications. We show that, although levels are very low in early-stage lung cancers, ctDNA is present prior to treatment in most patients and its presence is strongly prognostic. We also find that the majority of somatic mutations in the cell-free DNA (cfDNA) of patients with lung cancer and of risk-matched controls reflect clonal haematopoiesis and are non-recurrent. Compared with tumour-derived mutations, clonal haematopoiesis mutations occur on longer cfDNA fragments and lack mutational signatures that are associated with tobacco smoking. Integrating these findings with other molecular features, we develop and prospectively validate a machine-learning method termed 'lung cancer likelihood in plasma' (Lung-CLiP), which can robustly discriminate early-stage lung cancer patients from risk-matched controls. This approach achieves performance similar to that of tumour-informed ctDNA detection and enables tuning of assay specificity in order to facilitate distinct clinical applications. Our findings establish the potential of cfDNA for lung cancer screening and highlight the importance of risk-matching cases and controls in cfDNA-based screening studies.

Constitutive Interferon Pathway Activation in Tumors as an Efficacy Determinant Following Oncolytic Virotherapy

Background

Attenuated measles virus (MV) strains are promising agents currently being tested against solid tumors or hematologic malignancies in ongoing phase I and II clinical trials; factors determining oncolytic virotherapy success remain poorly understood, however.

Methods

We performed RNA sequencing and gene set enrichment analysis to identify pathways differentially activated in MV-resistant (n = 3) and -permissive (n = 2) tumors derived from resected human glioblastoma (GBM) specimens and propagated as xenografts (PDX). Using a unique gene signature we identified, we generated a diagonal linear discriminant analysis (DLDA) classification algorithm to predict MV responders and nonresponders, which was validated in additional randomly selected GBM and ovarian cancer PDX and 10 GBM patients treated with MV in a phase I trial. GBM PDX lines were also treated with the US Food and Drug Administration-approved JAK inhibitor, ruxolitinib, for 48 hours prior to MV infection and virus production, STAT1/3 signaling and interferon stimulated gene expression was assessed. All statistical tests were two-sided.

Results

Constitutive interferon pathway activation, as reflected in the DLDA algorithm, was identified as the key determinant for MV replication, independent of virus receptor expression, in MV-permissive and -resistant GBM PDXs. Using these lines as the training data for the DLDA algorithm, we confirmed the accuracy of our algorithm in predicting MV response in randomly selected GBM PDX ovarian cancer PDXs. Using the DLDA prediction algorithm, we demonstrate that virus replication in patient tumors is inversely correlated with expression of this resistance gene signature (ρ = -0.717, P = .03). In vitro inhibition of the interferon response pathway with the JAK inhibitor ruxolitinib was able to overcome resistance and increase virus production (1000-fold, P = .03) in GBM PDX lines.

Conclusions

These findings document a key mechanism of tumor resistance to oncolytic MV therapy and describe for the first time the development of a prediction algorithm to preselect for oncolytic treatment or combinatorial strategies.

Bacterial Single Cell Whole Transcriptome Amplification in Microfluidic Platform Shows Putative Gene Expression Heterogeneity

Single cell RNA sequencing is a technology that provides the capability of analyzing the transcriptome of a single cell from a population. So far, single cell RNA sequencing has been focused mostly on human cells due to the larger starting amount of RNA template for subsequent amplification. One of the major challenges of applying single cell RNA sequencing to microbial cells is to amplify the femtograms of the RNA template to obtain sufficient material for downstream sequencing with minimal contamination. To achieve this goal, efforts have been focused on multiround RNA amplification, but would introduce additional contamination and bias. In this work, we for the first time coupled a microfluidic platform with multiple displacement amplification technology to perform single cell whole transcriptome amplification and sequencing of Porphyromonas somerae, a microbe of interest in endometrial cancer, as a proof-of-concept demonstration of using single cell RNA sequencing tool to unveil gene expression heterogeneity in single microbial cells. Our results show that the bacterial single-cell gene expression regulation is distinct across different cells, supporting widespread heterogeneity.

A prospective genome-wide study of prostate cancer metastases reveals association of wnt pathway activation and increased cell cycle proliferation with primary resistance to abiraterone acetate-prednisone

Background

Genomic aberrations have been identified in metastatic castration-resistant prostate cancer (mCRPC), but molecular predictors of resistance to abiraterone acetate/prednisone (AA/P) treatment are not known.

Patients and methods

In a prospective clinical trial, mCRPC patients underwent whole-exome sequencing (n = 82) and RNA sequencing (n = 75) of metastatic biopsies before initiating AA/P with the objective of identifying genomic alterations associated with resistance to AA/P. Primary resistance was determined at 12 weeks of treatment using criteria for progression that included serum prostate-specific antigen measurement, bone and computerized tomography imaging and symptom assessments. Acquired resistance was determined using the end point of time to treatment change (TTTC), defined as time from enrollment until change in treatment from progressive disease. Associations of genomic and transcriptomic alterations with primary resistance were determined using logistic regression, Fisher's exact test, single and multivariate analyses. Cox regression models were utilized for determining association of genomic and transcriptomic alterations with TTTC.

Results

At 12 weeks, 32 patients in the cohort had progressed (nonresponders). Median study follow-up was 32.1 months by which time 58 patients had switched treatments due to progression. Median TTTC was 10.1 months (interquartile range: 4.4-24.1). Genes in the Wnt/β-catenin pathway were more frequently mutated and negative regulators of Wnt/β-catenin signaling were more frequently deleted or displayed reduced mRNA expression in nonresponders. Additionally, mRNA expression of cell cycle regulatory genes was increased in nonresponders. In multivariate models, increased cell cycle proliferation scores (≥ 50) were associated with shorter TTTC (hazard ratio = 2.11, 95% confidence interval: 1.17-3.80; P = 0.01).

Conclusions

Wnt/β-catenin pathway activation and increased cell cycle progression scores can serve as molecular markers for predicting resistance to AA/P therapy.

Predicting Treatment Response Based on RNA Expression in Large Datasets

PD-L1 expression levels derived from >16,000 samples guided the selection of tumor types likely to benefit from pembrolizuamb monotherapy in clinical trials. Although not fail-proof, FDA approvals for most of the prioritized indications speak to the power of RNA expression profiling and the value of large genomic datasets.See related article by Ayers et al., p. 1564.

Molecular signatures of multiple myeloma progression through single cell RNA-Seq

We used single cell RNA-Seq to examine molecular heterogeneity in multiple myeloma (MM) in 597 CD138 positive cells from bone marrow aspirates of 15 patients at different stages of disease progression. 790 genes were selected by coefficient of variation (CV) method and organized cells into four groups (L1–L4) using unsupervised clustering. Plasma cells from each patient clustered into at least two groups based on gene expression signature. The L1 group contained cells from all MGUS patients having the lowest expression of genes involved in the oxidative phosphorylation, Myc targets, and mTORC1 signaling pathways (p < 1.2 × 10⁻¹⁴). In contrast, the expression level of these pathway genes increased progressively and were the highest in L4 group containing only cells from MM patients with t(4;14) translocations. A 44 genes signature of consistently overexpressed genes among the four groups was associated with poorer overall survival in MM patients (APEX trial, p < 0.0001; HR, 1.83; 95% CI, 1.33–2.52), particularly those treated with bortezomib (p < 0.0001; HR, 2.00; 95% CI, 1.39–2.89). Our study, using single cell RNA-Seq, identified the most significantly affected molecular pathways during MM progression and provided a novel signature predictive of patient prognosis and treatment stratification.

Neoantigenic Potential of Complex Chromosomal Rearrangements in Mesothelioma

INTRODUCTION

Malignant pleural mesothelioma is a disease primarily associated with exposure to the carcinogen asbestos. Whereas other carcinogen-related tumors are associated with a high tumor mutation burden, mesothelioma is not. We sought to resolve this discrepancy.

METHODS

We used mate-pair (n = 22), RNA (n = 28), and T cell receptor sequencing along with in silico predictions and immunologic assays to understand how structural variants of chromosomes affect the transcriptome.

RESULTS

We observed that inter- or intrachromosomal rearrangements were present in every specimen and were frequently in a pattern of chromoanagenesis such as chromoplexy or chromothripsis. Transcription of rearrangement-related junctions was predicted to result in many potential neoantigens, some of which were proven to bind patient-specific major histocompatibility complex molecules and to expand intratumoral T cell clones. T cells responsive to these predicted neoantigens were also present in a patient's circulating T cell repertoire. Analysis of genomic array data from the mesothelioma cohort in The Cancer Genome Atlas suggested that multiple chromothriptic-like events negatively impact survival.

CONCLUSIONS

Our findings represent the discovery of potential neoantigen expression driven by structural chromosomal rearrangements. These results may have implications for the development of novel immunotherapeutic strategies and the selection of patients to receive immunotherapies.

Abstract A17: Establishment of a unique patient-derived tumor model positive for STRN-ALK fusion from a patient with stage IV lung adenocarcinoma

In lung adenocarcinoma, gene fusion involving the anaplastic lymphoma kinase (ALK) gene leads to overexpression of its tyrosine kinase domain, which serves as the oncogenic driver in lung tumors carrying such a molecular alteration. Similar fusions involving ROS1, RET and NTRK are also found in lung cancer and tumors with oncogenic fusions often respond well to targeted therapy. However, these tumors do develop resistance to the targeting drug and patients eventually relapse while on treatment. Although newer drugs have been developed, it remains a challenge to effectively control cancer recurrence after targeted therapy. In vivo tumor models carrying these fusions, particularly those that have progressed while on targeted treatment, will help to facilitate further drug evaluations prior to their administration to the patient as well as offer an opportunity to better understand the molecular mechanisms driving disease onset, progression, and drug resistance in tumors carrying the fusion of interest.

Here, we describe the first patient-derived tumor xenograft (PDX) carrying an STRN-ALK fusion. The patient was diagnosed with stage IV lung adenocarcinoma in 2011. She was initially treated with crizotinib and remained stable until 2015. At the time of disease progression, a biopsy was taken with patient consent under CT guidance from a chest wall mass using an 18-gauge needle and then implanted into two SCID-beige mice subcutaneously and into the fat pad of the kidney. Palpable tumor was observed in both mice at about 28 weeks post implantation and passaged at 36 weeks. They have been passaged two successive times and viably frozen. Pathologic analysis of the PDX showed that both tumors have the same mucinous morphology as observed in the original biopsy sample. RNA-seq and matepair analysis revealed that the tumors carry a rare fusion involving the first three exons of STRN and exon 20 to 3’ end of the ALK gene. The patient has been on ceritinib since January 2016 and remains stable as of last follow up.

Several unique features contributed to the success of this study. 1) The metastatic chest wall lesion used for the PDX developed was easily accessible on the chest wall using a relative large 18-gauge needle for CT-guided biopsies. 2) The slow growth of the tumor biopsy mirrors the clinical course of the disease observed in the patient. 3) The patient remains clinically stable and continues to respond to TKi. The availability of a viable PDX tumor model for this rare STRN-ALK fusion offers a unique opportunity to identify the most appropriate drug for the patient at the time of disease progression as well as the opportunity to examine the unique molecular features associated with this rare ALK fusion.

This work was supported in part by the National Foundation for Cancer Research Hillsberg Lung Cancer Translational Research Grant and by the Biomarker Discovery Program at the Mayo Clinic Center for Individualized Medicine.

Citation Format: HongZheng Res, XiaoNan Hou, Karlyn E. Pierson, Patric Eiken, Asha Nair, PingYang Yang, Joanne Eunhee Yi, Aaron S. Mansfield, Jin Jen. Establishment of a unique patient-derived tumor model positive for STRN-ALK fusion from a patient with stage IV lung adenocarcinoma [abstract]. In: Proceedings of the Fifth AACR-IASLC International Joint Conference: Lung Cancer Translational Science from the Bench to the Clinic; Jan 8-11, 2018; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2018;24(17_Suppl):Abstract nr A17.

Abstract 5726: Rearrangement-related peptides with neoantigenic potential in malignant pleural mesothelioma

Malignant pleural mesothelioma is a disease primarily associated with exposure to the carcinogen asbestos. Consistent with this carcinogenic exposure, cytogenetic analyses have identified multiple recurrent structural chromosomal abnormalities in this malignancy, but more recent high-throughput sequencing evaluations of point mutations suggest that there is a low mutational burden in mesothelioma. Since tumor mutational burden has been correlated with responses to treatment with immune checkpoint inhibitors such as nivolumab, it was not consistent that patients with mesothelioma and low mutation burdens would have similar response rates in clinical trials with immune checkpoint inhibitors as patients with non-small cell lung cancer which is associated with a high mutation burden. In order to reconcile these differences, and given the potential for an improved understanding of the molecular pathogenesis of mesothelioma to improve therapeutic options, we used mate-pair sequencing (MPseq) and RNA sequencing (RNAseq) to understand how structural variants affect the transcriptome. MPseq differs from standard next generation sequencing approaches by tiling the whole genome with larger fragments (2-5kb) to reliably detect structural variants such as insertions, deletions and rearrangements. Amongst 22 mesothelioma specimens there were 1535 chromosomal rearrangements (median 41, range 3-298 per specimen), that resulted in junctions or novel fusions of non-coding DNA or genes. Six-hundred thirty-seven of these rearrangements (median 22, range 5-103 range per specimen) resulted in novel fusions of genes. Many of these inter- or intra-chromosomal rearrangements were consistent with a pattern of chromoanagesis such as chromoplexy or chromothripsis. Chromosomal rearrangements detected by MPseq were used to guide analysis of RNAseq data and revealed that these chromosomal junctions resulted in the expression of 179 novel amino acid sequences (median 5, 0-51 range per specimen). To determine whether transcription of chromosomal rearrangement-related junctions have neoantigenic potential, we used in silico tools to determine whether any of the expressed junctions contained peptides that could be presented by patient-specific HLA molecules. The top candidate rearrangement-related peptides with neoantigenic potential bound patient-specific HLA molecules nearly as well or as well as a positive control in competitive binding assays. Our findings represent the discovery of potential neoantigen expression driven by structural chromosomal rearrangements. These results may have implications for the development of novel therapeutic strategies, the selection of patients to receive immunotherapy, and blood-based treatment monitoring strategies.

Citation Format: Aaron S. Mansfield, Tobias Peikert, James B. Smadbeck, Julia B. Udell, Farhad Kosari, Stephen J. Murphy, Hongzheng Ren, Vishnu V. Serla, Janet L. Schaefer Klein, Giannoula Karagouga, Faye R. Harris, Carlos Sosa, Sarah H. Johnson, Wendy Nevala, Svetomir N. Markovic, Aaron O. Bungum, Eric S. Edell, Haidong Dong, John C. Cheville, Marie Christine Aubry, Jin Jen, George Vasmatzis. Rearrangement-related peptides with neoantigenic potential in malignant pleural mesothelioma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5726.

A Comprehensive Approach to Sequence-oriented IsomiR annotation (CASMIR): demonstration with IsomiR profiling in colorectal neoplasia

Background

MicroRNA (miRNA) profiling is an important step in studying biological associations and identifying marker candidates. miRNA exists in isoforms, called isomiRs, which may exhibit distinct properties. With conventional profiling methods, limitations in assay and analysis platforms may compromise isomiR interrogation.

Results

We introduce a comprehensive approach to sequence-oriented isomiR annotation (CASMIR) to allow unbiased identification of global isomiRs from small RNA sequencing data. In this approach, small RNA reads are maintained as independent sequences instead of being summarized under miRNA names. IsomiR features are identified through step-wise local alignment against canonical forms and precursor sequences. Through customizing the reference database, CASMIR is applicable to isomiR annotation across species. To demonstrate its application, we investigated isomiR profiles in normal and neoplastic human colorectal epithelia. We also ran miRDeep2, a popular miRNA analysis algorithm to validate isomiRs annotated by CASMIR. With CASMIR, specific and biologically relevant isomiR patterns could be identified. We note that specific isomiRs are often more abundant than their canonical forms. We identify isomiRs that are commonly up-regulated in both colorectal cancer and advanced adenoma, and illustrate advantages in targeting isomiRs as potential biomarkers over canonical forms.

Conclusions

Studying miRNAs at the isomiR level could reveal new insight into miRNA biology and inform assay design for specific isomiRs. CASMIR facilitates comprehensive annotation of isomiR features in small RNA sequencing data for isomiR profiling and differential expression analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4794-7) contains supplementary material, which is available to authorized users.

CX3CR1 identifies PD-1 therapy-responsive CD8+ T cells that withstand chemotherapy during cancer chemoimmunotherapy

Although immune checkpoint inhibitors have resulted in durable clinical benefits in a subset of patients with advanced cancer, some patients who did not respond to initial anti-PD-1 therapy have been found to benefit from the addition of salvage chemotherapy. However, the mechanism responsible for the successful chemoimmunotherapy is not completely understood. Here we show that a subset of circulating CD8+ T cells expressing the chemokine receptor CX3CR1 are able to withstand the toxicity of chemotherapy and are increased in patients with metastatic melanoma who responded to chemoimmunotherapy (paclitaxel and carboplatin plus PD-1 blockade). These CX3CR1+CD8+ T cells have effector memory phenotypes and the ability to efflux chemotherapy drugs via the ABCB1 transporter. In line with clinical observation, our preclinical models identified an optimal sequencing of chemoimmunotherapy that resulted in an increase of CX3CR1+CD8+ T cells. Taken together, we found a subset of PD-1 therapy-responsive CD8+ T cells that were capable of withstanding chemotherapy and executing tumor rejection with their unique abilities of drug efflux (ABCB1), cytolytic activity (granzyme B and perforin), and migration to and retention (CX3CR1 and CD11a) at tumor sites. Future strategies to monitor and increase the frequency of CX3CR1+CD8+ T cells may help to design effective chemoimmunotherapy to overcome cancer resistance to immune checkpoint blockade therapy.

Detection and localization of surgically resectable cancers with a multi-analyte blood test

Earlier detection is key to reducing cancer deaths. Here, we describe a blood test that can detect eight common cancer types through assessment of the levels of circulating proteins and mutations in cell-free DNA. We applied this test, called CancerSEEK, to 1005 patients with nonmetastatic, clinically detected cancers of the ovary, liver, stomach, pancreas, esophagus, colorectum, lung, or breast. CancerSEEK tests were positive in a median of 70% of the eight cancer types. The sensitivities ranged from 69 to 98% for the detection of five cancer types (ovary, liver, stomach, pancreas, and esophagus) for which there are no screening tests available for average-risk individuals. The specificity of CancerSEEK was greater than 99%: only 7 of 812 healthy controls scored positive. In addition, CancerSEEK localized the cancer to a small number of anatomic sites in a median of 83% of the patients.

Contraction of T cell richness in lung cancer brain metastases

Very little is known about how the adaptive immune system responds to clonal evolution and tumor heterogeneity in non-small cell lung cancer. We profiled the T-cell receptor β complementarity determining region 3 in 20 patients with fully resected non-small cell lung cancer primary lesions and paired brain metastases. We characterized the richness, abundance and overlap of T cell clones between pairs, in addition to the tumor mutation burden and predicted neoantigens. We found a significant contraction in the number of unique T cell clones in brain metastases compared to paired primary cancers. The vast majority of T cell clones were specific to a single lesion, and there was minimal overlap in T cell clones between paired lesions. Despite the contraction in the number of T cell clones, brain metastases had higher non-synonymous mutation burdens than primary lesions. Our results suggest that there is greater richness of T cell clones in primary lung cancers than their paired metastases despite the higher mutation burden observed in metastatic lesions. These results may have implications for immunotherapy.

Gene Expression Signatures Characterized by Longitudinal Stability and Interindividual Variability Delineate Baseline Phenotypic Groups with Distinct Responses to Immune Stimulation

Human immunity exhibits remarkable heterogeneity among individuals, which engenders variable responses to immune perturbations in human populations. Population studies reveal that, in addition to interindividual heterogeneity, systemic immune signatures display longitudinal stability within individuals, and these signatures may reliably dictate how given individuals respond to immune perturbations. We hypothesize that analyzing relationships among these signatures at the population level may uncover baseline immune phenotypes that correspond with response outcomes to immune stimuli. To test this, we quantified global gene expression in peripheral blood CD4⁺ cells from healthy individuals at baseline and following CD3/CD28 stimulation at two time points 1 mo apart. Systemic CD4⁺ cell baseline and poststimulation molecular immune response signatures (MIRS) were defined by identifying genes expressed at levels that were stable between time points within individuals and differential among individuals in each state. Iterative differential gene expression analyses between all possible phenotypic groupings of at least three individuals using the baseline and stimulated MIRS gene sets revealed shared baseline and response phenotypic groupings, indicating the baseline MIRS contained determinants of immune responsiveness. Furthermore, significant numbers of shared phenotype-defining sets of determinants were identified in baseline data across independent healthy cohorts. Combining the cohorts and repeating the analyses resulted in identification of over 6000 baseline immune phenotypic groups, implying that the MIRS concept may be useful in many immune perturbation contexts. These findings demonstrate that patterns in complex gene expression variability can be used to define immune phenotypes and discover determinants of immune responsiveness.

Trinucleotide Repeat Expansion in the Transcription Factor 4 (TCF4) Gene Leads to Widespread mRNA Splicing Changes in Fuchs' Endothelial Corneal Dystrophy

Purpose

To identify RNA missplicing events in human corneal endothelial tissue isolated from Fuchs' endothelial corneal dystrophy (FECD).

Methods

Total RNA was isolated and sequenced from corneal endothelial tissue obtained during keratoplasty from 12 patients with FECD and 4 patients undergoing keratoplasty or enucleation for other indications. The length of the trinucleotide repeat (TNR) CTG in the transcription factor 4 (TCF4) gene was determined using leukocyte-derived DNA analyzed by a combination of Southern blotting and Genescan analysis. Commercial statistical software was used to quantify expression of alternatively spliced genes. Validation of selected alternative splicing events was performed by using RT-PCR. Gene sets identified were analyzed for overrepresentation using Web-based analysis system.

Results

Corneal endothelial tissue from FECD patients containing a CTG TNR expansion sequence in the TCF4 gene revealed widespread changes in mRNA splicing, including a novel splicing event involving FGFR2. Differential splicing of NUMA1, PPFIBP1, MBNL1, and MBNL2 transcripts were identified in all FECD samples containing a TNR expansion. The differentially spliced genes were enriched for products that localize to the cell cortex and bind cytoskeletal and cell adhesion proteins.

Conclusions

Corneal endothelium from FECD patients harbors a unique signature of mis-splicing events due to CTG TNR expansion in the TCF4 gene, consistent with the hypothesis that RNA toxicity contributes to the pathogenesis of FECD. Changes to the endothelial barrier function, a known event in the development of FECD, was identified as a key biological process influenced by the missplicing events.

CK1δ/GSK3β/FBXW7α axis promotes degradation of the ZNF322A oncoprotein to suppress lung cancer progression

Overexpression of Cys2His2 zinc-finger 322A (ZNF322A) oncogenic transcription factor is associated with lung tumorigenesis. However, the mechanism of ZNF322A overexpression remains poorly understood. Here, we discover that protein stability of ZNF322A is regulated by coordinated phosphorylation and ubiquitination through the CK1δ/GSK3β/FBXW7α axis. CK1δ and GSK3β kinases sequentially phosphorylate ZNF322A at serine-396 and then serine-391. Moreover, the doubly phosphorylated ZNF322A protein creates a destruction motif for the ubiquitin ligase FBXW7α leading to ZNF322A protein destruction. Overexpression of FBXW7α induces ZNF322A protein degradation, thereby blocks ZNF322A transcription activity and suppresses ZNF322A-induced tumor growth and metastasis in vitro and in vivo. Clinically, overexpression of ZNF322A correlates with low FBXW7α or defective CK1δ/GSK3β-mediated phosphorylation in lung cancer patients. Multivariate Cox regression analysis indicates that patients with ZNF322A high/FBXW7 low expression profile can be used as an independent factor to predict the clinical outcome in lung cancer patients. Our results reveal a new mechanism of ZNF322A oncoprotein destruction regulated by the CK1δ/GSK3β/FBXW7α axis. Deregulation of this signaling axis results in ZNF322A overexpression and promotes cancer progression.

Meeting Report of the Fifth International Cancer Epigenetics Conference in Beijing, China, October 2016

Fifth International Cancer Epigenetics Conference, Beijing, China, 21-23 October 2016 This meeting reported many new findings in the field of cancer epigenetics, including basic science, translational and clinical studies. In this report, we summarize some of the main advancements and prospects in cancer epigenetics presented at this meeting.

Androgen Receptor Variant AR-V9 Is Coexpressed with AR-V7 in Prostate Cancer Metastases and Predicts Abiraterone Resistance

Purpose: Androgen receptor (AR) variant AR-V7 is a ligand-independent transcription factor that promotes prostate cancer resistance to AR-targeted therapies. Accordingly, efforts are under way to develop strategies for monitoring and inhibiting AR-V7 in castration-resistant prostate cancer (CRPC). The purpose of this study was to understand whether other AR variants may be coexpressed with AR-V7 and promote resistance to AR-targeted therapies.Experimental Design: We utilized complementary short- and long-read sequencing of intact AR mRNA isoforms to characterize AR expression in CRPC models. Coexpression of AR-V7 and AR-V9 mRNA in CRPC metastases and circulating tumor cells was assessed by RNA-seq and RT-PCR, respectively. Expression of AR-V9 protein in CRPC models was evaluated with polyclonal antisera. Multivariate analysis was performed to test whether AR variant mRNA expression in metastatic tissues was associated with a 12-week progression-free survival endpoint in a prospective clinical trial of 78 CRPC-stage patients initiating therapy with the androgen synthesis inhibitor, abiraterone acetate.Results: AR-V9 was frequently coexpressed with AR-V7. Both AR variant species were found to share a common 3' terminal cryptic exon, which rendered AR-V9 susceptible to experimental manipulations that were previously thought to target AR-V7 uniquely. AR-V9 promoted ligand-independent growth of prostate cancer cells. High AR-V9 mRNA expression in CRPC metastases was predictive of primary resistance to abiraterone acetate (HR = 4.0; 95% confidence interval, 1.31-12.2; P = 0.02).Conclusions: AR-V9 may be an important component of therapeutic resistance in CRPC. Clin Cancer Res; 23(16); 4704-15. ©2017 AACR.

Oncogenic MCT-1 activation promotes YY1-EGFR-MnSOD signaling and tumor progression

Tumor cells often produce high levels of reactive oxygen species (ROS) and display an increased ROS scavenging system. However, the molecular mechanism that balances antioxidative and oxidative stress in cancer cells is unclear. Here, we determined that oncogenic multiple copies in T-cell malignancy 1 (MCT-1) activity promotes the generation of intracellular ROS and mitochondrial superoxide. Overexpression of MCT-1 suppresses p53 accumulation but elevates the manganese-dependent superoxide dismutase (MnSOD) level via the YY1-EGFR signaling cascade, which protects cells against oxidative damage. Conversely, restricting ROS generation and/or targeting YY1 in lung cancer cells effectively inhibits the EGFR-MnSOD signaling pathway and cell invasiveness induced by MCT-1. Significantly, MCT-1 overexpression in lung cancer cells promotes tumor progression, necrosis and angiogenesis, and increases the number of tumor-promoting M2 macrophages and cancer-associated fibroblasts in the microenvironment. Clinical evidence further confirms that high expression of MCT-1 is associated with an increase in YY1, EGFR and MnSOD expression, accompanied by tumor recurrence, poor overall survival and EGFR mutation status in patients with lung cancers. Together, these data indicate that the MCT-1 oncogenic pathway is implicated in oxidative metabolism and lung carcinogenesis.

Oncoprotein ZNF322A transcriptionally deregulates alpha-adducin, cyclin D1 and p53 to promote tumor growth and metastasis in lung cancer

This corrects the article DOI: 10.1038/onc.2015.296.

Custom Gene Capture and Next-Generation Sequencing to Resolve Discordant ALK Status by FISH and IHC in Lung Adenocarcinoma

INTRODUCTION

We performed a genomic study in lung adenocarcinoma cases with discordant anaplastic lymphoma receptor tyrosine kinase gene (ALK) status by fluorescent in situ hybridization (FISH) and immunohistochemical (IHC) analysis.

METHODS

DNA from formalin-fixed paraffin-embedded tissues of 16 discordant (four FISH-positive/IHC-negative and 12 FISH-negative/IHC-positive) cases by Vysis ALK Break Apart FISH and ALK IHC testing (ALK1 clone) were subjected to whole gene capture and next-generation sequencing (NGS) of nine genes, including ALK, echinoderm microtubule associated protein like 4 gene (EML4), kinesin family member 5B gene (KIF5B), staphylococcal nuclease and tudor domain containing 1 gene (SND1), BRAF, ret proto-oncogene (RET), ezrin gene (EZR), ROS1, and telomerase reverse transcriptase (TERT). All discordant cases (except one FISH-negative/IHC-positive case without sufficient tissue) were analyzed by IHC with D5F3 antibody. In one case with fresh frozen tissue, whole transcriptome sequencing was also performed. Twenty-six concordant (16 FISH-positive/IHC-positive and 10 FISH-negative/IHC-negative) cases were included as controls.

RESULTS

In four ALK FISH-positive/IHC-negative cases, no EML4-ALK fusion gene was observed by NGS, but in one case using fresh frozen tissue, we identified EML4-baculoviral AIP repeat containing 6 gene (BIRC6) and AP2 associated kinase 1 gene (AAK1)-ALK fusion genes. Whole transcriptome sequencing revealed a highly expressed EML4-BIRC6 fusion transcript and a minimally expressed AAK1 transcript. Among the 12 FISH-negative/IHC-positive cases, no evidence of ALK gene rearrangement was detected by NGS. Eleven of 12 FISH-negative/IHC-positive cases detected by ALK1 clone were concordant by repeat ALK IHC with D5F3 antibody (i.e., FISH-negative/IHC-negative by D5F3 clone). Among the 16 ALK FISH-positive/IHC-positive positive controls, whole gene capture identified ALK gene fusion in 15 cases, including in one case with Huntington interacting protein 1 gene (HIP1)-ALK. No ALK fusion gene was observed in any of the 10 FISH-negative/IHC-negative cases. Other fusion genes involving ROS1, EZR, BRAF, and SND1 were also found.

CONCLUSIONS

ALK FISH results appeared to be false-positive in three of four FISH-positive/IHC-negative cases, whereas no false-negative ALK FISH case was identified among 12 ALK FISH-negative/IHC-positive cases by ALK1 clone, which was in keeping with the concordant FISH-negative/IHC-negative status by D5F3 clone. Our targeted whole gene capture approach using formalin-fixed paraffin embedded samples was effective for detecting rearrangements involving ALK and other actionable oncogenes.

Impact of RNA degradation on fusion detection by RNA-seq

Background

RNA-seq is a well-established method for studying the transcriptome. Popular methods for library preparation in RNA-seq such as Illumina TruSeq® RNA v2 kit use a poly-A pulldown strategy. Such methods can cause loss of coverage at the 5′ end of genes, impacting the ability to detect fusions when used on degraded samples. The goal of this study was to quantify the effects RNA degradation has on fusion detection when using poly-A selected mRNA and to identify the variables involved in this process.

Results

Using both artificially and naturally degraded samples, we found that there is a reduced ability to detect fusions as the distance of the breakpoint from the 3′ end of the gene increases. The median transcript coverage decreases exponentially as a function of the distance from the 3′ end and there is a linear relationship between the coverage decay rate and the RNA integrity number (RIN). Based on these findings we developed plots that show the probability of detecting a gene fusion (“sensitivity”) as a function of the distance of the fusion breakpoint from the 3′ end.

Conclusions

This study developed a strategy to assess the impact that RNA degradation has on the ability to detect gene fusions by RNA-seq.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3161-9) contains supplementary material, which is available to authorized users.

Chromoplectic TPM3-ALK rearrangement in a patient with inflammatory myofibroblastic tumor who responded to ceritinib after progression on crizotinib

Ceritinib resulted in a significant, durable response of a metastatic inflammatory myofibroblastic tumor (IMT) after failure of crizotinib. A chromoplectic TPM3–ALK rearrangement involving many known oncogenes was found in the residual IMT. Ceritinib may be useful for patients with IMT after failure of crizotinib, and chromoplexy may have a role in the oncogenesis or treatment resistance of IMTs.

Background

Inflammatory myofibroblastic tumors (IMTs) are rare sarcomas that can occur at any age. Surgical resection is the primary treatment for patients with localized disease; however, these tumors frequently recur. Less commonly, patients with IMTs develop or present with metastatic disease. There is no standard of care for these patients and traditional cytotoxic therapy is largely ineffective. Most IMTs are associated with oncogenic ALK, ROS1 or PDGFRβ fusions and may benefit from targeted therapy.

Patient and methods

We sought to understand the genomic abnormalities of a patient who presented for management of metastatic IMT after progression of disease on crizotinib and a significant and durable partial response to the more potent ALK inhibitor ceritinib.

Results

The residual IMT was resected based on the recommendations of a multidisciplinary tumor sarcoma tumor board and analyzed by whole-genome mate pair sequencing. Analysis of the residual, resected tumor identified a chromoplectic TPM3–ALK rearrangement that involved many other known oncogenes and was confirmed by rtPCR.

Conclusions

In our analysis of the treatment-resistant, residual IMT, we identified a complex pattern of genetic rearrangements consistent with chromoplexy. Although it is difficult to know for certain if these chromoplectic rearrangements preceded treatment, their presence suggests that chromoplexy has a role in the oncogenesis of IMTs. Furthermore, this patient's remarkable response suggests that ceritinib should be considered as an option after progression on crizotinib for patients with metastatic or unresectable IMT and ALK mutations.

Clinical Outcomes of TP53 Mutations in Cancers

High-throughput sequencing of cancer genomes is increasingly becoming an essential tool of clinical oncology that facilitates target identification and targeted therapy within the context of precision medicine. The cumulative profiles of somatic mutations in cancer yielded by comprehensive molecular studies also constitute a fingerprint of historical exposures to exogenous and endogenous mutagens, providing insight into cancer evolution and etiology. Mutational signatures that were first established by inspection of the TP53 gene somatic landscape have now been confirmed and expanded by comprehensive sequencing studies. Further, the degree of granularity achieved by deep sequencing allows detection of low-abundance mutations with clinical relevance. In tumors, they represent the emergence of small aggressive clones; in normal tissues, they signal a mutagenic exposure related to cancer risk; and, in blood, they may soon become effective surveillance tools for diagnostic purposes and for monitoring of cancer prognosis and recurrence.

Abstract B34: Integrated Approaches to Treating Lung Adenocarcinoma Resistant to Targeted Therapy

Mutations in the EGFR gene or transcript fusion involving EML4-ALK result in oncogenic activation of these oncogenes driving lung adenocarcinoma growth in tumors carrying such an alteration. For lung adenocarcinoma with these genetic changes, small-molecule tyrosine kinase inhibitors have been highly effective at reducing tumor burden and improve the outcome of the patients. Unfortunately, drug resistance eventually develops in these tumors and challenges remain in controlling lung cancer recurrence after targeted therapy.

To overcome recurrence, we initiated a prove-of-principle study to evaluate the feasibility of an integrated strategy utilizing genomic profiles of the tumor and patient-specific tumor xenografts derived (PDX) from biopsies for ex vivo evaluation of antitumor drugs to help guide personalized treatment of lung cancer in patients who have developed resistance to targeted therapy drugs. Our study has three main objectives. 1) Use tumor biopsies obtained at the time of recurrence for oncogene mutation and RNAseq analyses to identify molecular changes in oncogenes and gene pathways that can be potentially targeted. 2) Evaluate drug efficacy and optimize personalized therapy for each patient using PDX models. 3) Assess clinical feasibility and our experience using integrated approaches for lung cancer patients who failed targeted therapy.

This study was approved by the Mayo Clinic Institutional Review Board (IRB) and utilizes both clinical and research biopsies. A dedicated nurse study coordinator reviews clinical patient list and history on weekly basis and informs the study team of each potential candidate patient. A total of 30 patients having ALK positive lung cancers have been identified and followed up from nearly 500 potentially ALK positive cases between April 2014 to Sept. 2015. Most patients were never smokers but six were former smokers and two were current smokers. Age at diagnoses ranged from 27-78 years (median = 61). Seven patients (23%) were 45 years or younger at the time of diagnosis. A total of 22 patients are currently being treated with crizotinib while eight are on ceritinib or alectinib. For each biopsy, tumor tissues are obtained using a 20 gauge needle, transferred on ice in preserving media and implanted within one hour into 6-8 week old NOD/SCID mice. Many patients have been on treatment for 2-3 years with stable disease so they do not require biopsy. Using a similar strategy, we also obtained biopsies from patients with EGFR gene mutations in their tumors and have progressed while on targeted therapy.

A total of seven cases have been implanted. We will report our experiences in patient selection, clinical follow up, patient consent, PDX development, time from biopsy to tumor establishment, and the results of molecular analyses. Our study enabled us to gain new insights regarding the molecular changes associated with recurring tumors after they have failed targeted therapy as well as clinical experiences on how to utilize state-of-the art approaches and comprehensive genomic information to further improve cancer patient care upon disease progression.

Acknowledgements: This work is supported in part by the Hillsberg Award from the National Foundation for Cancer Research and by the Biomarker Discovery Program at the Mayo Clinic Center for Individualized Medicine.

Citation Format: Jin Jen, Aaron Mansfield, Patrick W. Eiken, Shawn Stoddard, Karlyn Pierson, Xiaonan Hou, Hong Zheng Ren, Julian Molina, Joanne Yi, Ping Yang. Integrated Approaches to Treating Lung Adenocarcinoma Resistant to Targeted Therapy. [abstract]. In: Proceedings of the AACR Special Conference: Patient-Derived Cancer Models: Present and Future Applications from Basic Science to the Clinic; Feb 11-14, 2016; New Orleans, LA. Philadelphia (PA): AACR; Clin Cancer Res 2016;22(16_Suppl):Abstract nr B34.

Adult-Onset Hemiplegic Migraine with Cortical Enhancement and Oedema

We present genetically identical twin patients who experienced late-onset migraine with visual and somatosensory auras and later developed hemiplegic migraines associated with severe cortical oedema and enhancement. Both positron emission tomography and electroencephalography showed an increase in activity contralateral to the hemiplegic side. Brain biopsy during the attack showed reactive astrogliosis and microgliosis. Mutations in CACNA1A, ATP1A2, SLC1A3 and NOTCH3 were ruled out by sequencing. This report shows the clinical and genetic evaluation of a severe form of familial hemiplegic migraine as well as the evolution of the imaging changes.

Measure transcript integrity using RNA-seq data

Background

Stored biological samples with pathology information and medical records are invaluable resources for translational medical research. However, RNAs extracted from the archived clinical tissues are often substantially degraded. RNA degradation distorts the RNA-seq read coverage in a gene-specific manner, and has profound influences on whole-genome gene expression profiling.

Result

We developed the transcript integrity number (TIN) to measure RNA degradation. When applied to 3 independent RNA-seq datasets, we demonstrated TIN is a reliable and sensitive measure of the RNA degradation at both transcript and sample level. Through comparing 10 prostate cancer clinical samples with lower RNA integrity to 10 samples with higher RNA quality, we demonstrated that calibrating gene expression counts with TIN scores could effectively neutralize RNA degradation effects by reducing false positives and recovering biologically meaningful pathways. When further evaluating the performance of TIN correction using spike-in transcripts in RNA-seq data generated from the Sequencing Quality Control consortium, we found TIN adjustment had better control of false positives and false negatives (sensitivity = 0.89, specificity = 0.91, accuracy = 0.90), as compared to gene expression analysis results without TIN correction (sensitivity = 0.98, specificity = 0.50, accuracy = 0.86).

Conclusion

TIN is a reliable measurement of RNA integrity and a valuable approach used to neutralize in vitro RNA degradation effect and improve differential gene expression analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0922-z) contains supplementary material, which is available to authorized users.

Circulating tumor cells are associated with poor overall survival in patients with cholangiocarcinoma

Circulating tumor cells (CTCs) in blood are associated with poor survival of patients with breast, prostate, or colon cancer. We hypothesized that CTCs are associated with poor survival of patients with cholangiocarcinoma (CCA). Eighty-eight patients with CCA were prospectively enrolled at Mayo Clinic Rochester between June 2010 and September 2014. The CellSearch system by Veridex was used for detection of CTCs in peripheral blood. Associations between CTC, patient and tumor characteristics, and survival were examined using the Cox's proportional hazards model. Fifteen patients (17%) were positive for CTC ≥2 and 8 patients (9%) for CTC ≥5. CTCs were associated with tumor extent. CTC ≥2 (hazard ratio [HR]: 2.5; 95% confidence interval [CI]: 1.1-5.4; P = 0.02) and CTC ≥5 (HR, 4.1; 95% CI: 1.4-10.8; P = 0.01) were both independent predictors of survival. In subgroup analyses, CTC ≥2 (HR, 8.2; 95% CI: 1.8-57.5; P < 0.01) and CTC ≥5 (HR, 7.7; 95% CI: 1.4-42.9; P = 0.02) were both associated with shorter survival among patients with metastasis. There was a trend toward association of CTC ≥5 with shorter survival in patients with nonmetastatic CCA (HR, 4.3; 95% CI: 1.0-13.8; P = 0.06). CTC ≥2 (HR, 10.5; 95% CI: 2.2-40.1; P < 0.01) and CTC ≥5 (HR, 10.2; 95% CI: 1.5-42.3; P = 0.02) were both associated with shorter survival among patients with perihilar/distal CCA. CTC ≥5 was associated with shorter survival of patients with intrahepatic CCA (HR, 4.2; 95% CI: 1.1-14.1; P = 0.04).

CONCLUSION

CTCs were associated with more-aggressive tumor characteristics and independently associated with survival in patients with CCA. Assessment of CTCs may be useful for identifying CCA patients at risk of early mortality.

Parkin Regulates Mitosis and Genomic Stability through Cdc20/Cdh1

Mutations in the E3 ubiquitin ligase Parkin have been linked to familial Parkinson's disease. Parkin has also been implicated in mitosis through mechanisms that are unclear. Here we show that Parkin interacts with anaphase promoting complex/cyclosome (APC/C) coactivators Cdc20 and Cdh1 to mediate the degradation of several key mitotic regulators independent of APC/C. We demonstrate that ordered progression through mitosis is orchestrated by two distinct E3 ligases through the shared use of Cdc20 and Cdh1. Furthermore, Parkin is phosphorylated and activated by polo-like kinase 1 (Plk1) during mitosis. Parkin deficiency results in overexpression of its substrates, mitotic defects, genomic instability, and tumorigenesis. These results suggest that the Parkin-Cdc20/Cdh1 complex is an important regulator of mitosis.

Expression Analysis of p16, c-Myc, and mSin3A in Non-small Cell Lung Cancer by Computer Aided Scoring and Analysis (CASA)

BACKGROUND

Immunohistochemical analysis (IHC) of tissue microarray (TMA) slides enables large sets of tissue samples to be analyzed simultaneously on a single slide. However, manual evaluation of small cores on a TMA slide is time consuming and error prone.

METHODS

We describe a computer aided scoring and analysis (CASA) method to allow facile and reliable scoring of IHC staining using TMA containing 300 non-small cell lung cancer (NSCLC) cases. In the two previous published papers utilizing our TMA slides of lung cancer we examined 18 proteins involved in the chromatin machinery. We developed our study using more proteins of the chromatin complex and several transcription factors that facilitate the chromatin machinery. Then, a total of 78 antibodies were evaluated by CASA to derive a normalized intensity value that correlated with the overall staining status of the targeting protein. The intensity values for TMA cores were then examined for association to clinical variables and predictive significance individually and with other factors. RESULTs: Using our TMA, the intensity of several protein pairs were significantly correlated with an increased risk of death in NSCLC. These included c-Myc with p16, mSin3A with p16 and c-Myc with mSinA. Predictive values of these pairs remained significant when evaluated based on standard IHC scores.

CONCLUSIONS

Our results demonstrate the usefulness of CASA as a valuable tool for systematic assessment of TMA slides to identify potential predictive biomarkers using a large set of primary human tissues.