Tumor-Educated Platelet RNA for the Detection and (Pseudo)progression Monitoring of Glioblastoma.

Tumor-educated platelets (TEPs) are potential biomarkers for cancer diagnostics. We employ TEP-derived RNA panels, determined by swarm intelligence, to detect and monitor glioblastoma. We assessed specificity by comparing the spliced RNA profile of TEPs from glioblastoma patients with multiple sclerosis and brain metastasis patients (validation series, n = 157; accuracy, 80%; AUC, 0.81 [95% CI, 0.74-0.89; p < 0.001]). Second, analysis of patients with glioblastoma versus asymptomatic healthy controls in an independent validation series (n = 347) provided a detection accuracy of 95% and AUC of 0.97 (95% CI, 0.95-0.99; p < 0.001). Finally, we developed the digitalSWARM algorithm to improve monitoring of glioblastoma progression and demonstrate that the TEP tumor scores of individual glioblastoma patients represent tumor behavior and could be used to distinguish false positive progression from true progression (validation series, n = 20; accuracy, 85%; AUC, 0.86 [95% CI, 0.70-1.00; p < 0.012]). In conclusion, TEPs have potential as a minimally invasive biosource for blood-based diagnostics and monitoring of glioblastoma patients.

Expression of let-7i and miR-192 is associated with resistance to cisplatin-based chemoradiotherapy in patients with larynx and hypopharynx cancer.

OBJECTIVES: The majority of patients with locally advanced larynx or hypopharynx squamous cell carcinoma are treated with organ-preserving chemoradiotherapy (CRT). Clinical outcome following CRT varies greatly. We hypothesized that tumor microRNA (miRNA) expression is predictive for outcome following CRT. METHODS: Next-generation sequencing (NGS) miRNA profiling was performed on 37 formalin-fixed paraffin-embedded (FFPE) tumor samples. Patients with a recurrence-free survival (RFS) of less than 2 years and patients with late/no recurrence within 2 years were compared by differential expression analysis. Tumor-specific miRNAs were selected based on normal mucosa miRNA expression data from The Cancer Genome Atlas database. A model was constructed to predict outcome using group-regularized penalized logistic ridge regression. Candidate miRNAs were validated by RT-qPCR in the initial sample set as well as in 46 additional samples. RESULTS: Thirteen miRNAs were differentially expressed (p < 0.05, FDR < 0.1) according to outcome group. Initial class prediction in the NGS cohort (n = 37) resulted in a model combining five miRNAs and disease stage, able to predict CRT outcome with an area under the curve (AUC) of 0.82. In the RT-qPCR cohort (n = 83), 25 patients (30%) experienced early recurrence (median RFS 8 months; median follow-up 42 months). Class prediction resulted in a model combining let-7i-5p, miR-192-5p and disease stage, able to discriminate patients with good versus poor clinical outcome (AUC:0.80). CONCLUSION: The combined miRNA expression and disease stage prediction model for CRT outcome is superior to using either factor alone. This study indicates NGS miRNA profiling using FFPE specimens is feasible, resulting in clinically relevant biomarkers.

Swarm Intelligence-Enhanced Detection of Non-Small-Cell Lung Cancer Using Tumor-Educated Platelets.

Blood-based liquid biopsies, including tumor-educated blood platelets (TEPs), have emerged as promising biomarker sources for non-invasive detection of cancer. Here we demonstrate that particle-swarm optimization (PSO)-enhanced algorithms enable efficient selection of RNA biomarker panels from platelet RNA-sequencing libraries (n = 779). This resulted in accurate TEP-based detection of early- and late-stage non-small-cell lung cancer (n = 518 late-stage validation cohort, accuracy, 88%; AUC, 0.94; 95% CI, 0.92-0.96; p < 0.001; n = 106 early-stage validation cohort, accuracy, 81%; AUC, 0.89; 95% CI, 0.83-0.95; p < 0.001), independent of age of the individuals, smoking habits, whole-blood storage time, and various inflammatory conditions. PSO enabled selection of gene panels to diagnose cancer from TEPs, suggesting that swarm intelligence may also benefit the optimization of diagnostics readout of other liquid biopsy biosources.

Apelin: A putative novel predictive biomarker for bevacizumab response in colorectal cancer.

Bevacizumab (bvz) is currently employed as an anti-angiogenic therapy across several cancer indications. Bvz response heterogeneity has been well documented, with only 10-15% of colorectal cancer (CRC) patients benefitting in general. For other patients, clinical efficacy is limited and side effects are significant. This reinforces the need for a robust predictive biomarker of response. To identify such a biomarker, we performed a DNA microarray-based transcriptional profiling screen with primary endothelial cells (ECs) isolated from normal and tumour colon tissues. Thirteen separate populations of tumour-associated ECs and 10 of normal ECs were isolated using fluorescence-activated cell sorting. We hypothesised that VEGF-induced genes were overexpressed in tumour ECs; these genes could relate to bvz response and serve as potential predictive biomarkers. Transcriptional profiling revealed a total of 2,610 differentially expressed genes when tumour and normal ECs were compared. To explore their relation to bvz response, the mRNA expression levels of top-ranked genes were examined using quantitative PCR in 30 independent tumour tissues from CRC patients that received bvz in the adjuvant setting. These analyses revealed that the expression of MMP12 and APLN mRNA was significantly higher in bvz non-responders compared to responders. At the protein level, high APLN expression was correlated with poor progression-free survival in bvz-treated patients. Thus, high APLN expression may represent a novel predictive biomarker for bvz unresponsiveness.

RNA-Seq of Tumor-Educated Platelets Enables Blood-Based Pan-Cancer, Multiclass, and Molecular Pathway Cancer Diagnostics.

Tumor-educated blood platelets (TEPs) are implicated as central players in the systemic and local responses to tumor growth, thereby altering their RNA profile. We determined the diagnostic potential of TEPs by mRNA sequencing of 283 platelet samples. We distinguished 228 patients with localized and metastasized tumors from 55 healthy individuals with 96% accuracy. Across six different tumor types, the location of the primary tumor was correctly identified with 71% accuracy. Also, MET or HER2-positive, and mutant KRAS, EGFR, or PIK3CA tumors were accurately distinguished using surrogate TEP mRNA profiles. Our results indicate that blood platelets provide a valuable platform for pan-cancer, multiclass cancer, and companion diagnostics, possibly enabling clinical advances in blood-based "liquid biopsies".

DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges, including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with ∼0.1× genome coverage. We improve on previous methods by first implementing a combined correction for sequence mappability and GC content, and second, by applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions was previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1000 samples, most of which were obtained from the fixed tissue archives of more than 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches and better copy number data than high-resolution microarrays at a substantially lower cost.

CD98 marks a subpopulation of head and neck squamous cell carcinoma cells with stem cell properties.

Patients with advanced head and neck squamous cell carcinomas (HNSCCs) are often treated with concomitant chemotherapy and radiotherapy, but only 50% is cured. A possible explanation for treatment failure is therapy resistance of the cancer stem cells (CSCs). The application of compounds specifically targeting these CSCs, in addition to routinely used therapeutics, would likely improve clinical outcome. We demonstrate that the previously described monoclonal antibody K984 recognizes the CD98 cell surface protein, which is specifically expressed by cells forming the squamous basal cell layer, the region where the squamous stem cells reside. Moreover, CD98 is highly resistant to the proteolytic enzymes required for CSC enrichment procedures. We show that CD98(high) cells, in contrast to CD98(low) cells, are able to generate tumors in immunodeficient mice, indicating that CD98(high) cells have stem cell characteristics. Furthermore, the CD98(high) subpopulation expresses high levels of cell cycle control and DNA repair genes, while the CD98(low) fraction shows expression patterns that represent the more differentiated cells forming the bulk of the tumor. CD98 is a promising CSC enrichment marker in HNSCC. Our data support the CSC concept in head and neck cancer and the potential relevance of these cells for treatment outcome.

Gene expression profile in the diaphragm following contractile inactivity during thoracic surgery.

Introduction. Recent work revealed the development of marked muscle fiber weakness in the diaphragm, but not in the non-respiratory latissimus dorsi, during thoracic surgery. To disentangle the molecular processes that underlie the development of diaphragm muscle fiber weakness during thoracic surgery, we studied changes in the gene expression profile. Methods. Serial biopsies from the diaphragm and the latissimus dorsi muscle were obtained from four patients during thoracotomy for resection of a tumor in the right lung. Biopsies were taken as soon as the diaphragm had been exposed (t0) and again after two hours (t2). Gobal differences in gene expression in diaphragm biopsies were assessed by microarray analysis. Results. 346 differentially expressed gene transcripts were found in the diaphragm at t2 vs. t0. Pathway analysis revealed that genes associated with inflammation (83 genes; p<0.0001) and cell death (118 genes, p<0.0001) pathways were significantly overexpressed at t2. Of the 346 differentially expressed genes in the diaphragm at t2, 258 were also differential in the latissimus dorsi muscle, with the direction of change being identical for all differentially expressed genes. In addition, latissimus dorsi showed exclusive upregula-ton of negative regulators of cell death. Conclusions. Two hours of thoracic surgery result in rapid and profound changes in expression of inflammatory response and apoptotic genes in the diaphragm. The apoptotic response was stronger in the diaphragm than in the latissiums dorsi. These findings suggest that the development of selective diaphragm muscle fiber weakness in these patients might be related to an exaggerated apoptotic response.

Candidate driver genes in focal chromosomal aberrations of stage II colon cancer.

Chromosomal instable colorectal cancer is marked by specific large chromosomal copy number aberrations. Recently, focal aberrations of 3 Mb or smaller have been identified as a common phenomenon in cancer. Inherent to their limited size, these aberrations harbour one or few genes. The aim of this study was to identify recurrent focal chromosomal aberrations and their candidate driver genes in a well-defined series of stage II colon cancers and assess their potential clinical relevance. High-resolution DNA copy number profiles were obtained from 38 formalin-fixed, paraffin-embedded colon cancer samples with matched normal mucosa as a reference using array comparative genomic hybridization. In total, 81 focal chromosomal aberrations were identified that harboured 177 genes. Statistical validation of focal aberrations and identification of candidate driver genes were performed by enrichment analysis and mapping copy number and mutation data of colorectal, breast, and pancreatic cancer and glioblastomas to loci of focal aberrations in stage II colon cancer. This analysis demonstrated a significant overlap with previously identified focal amplifications in colorectal cancer, but not with cancers from other sites. In contrast, focal deletions seemed less tumour type-specific since they also showed significant overlap with focal deletions of other sites. Focal deletions detected were significantly enriched for cancer genes and genes frequently mutated in colorectal cancer. The mRNA expression of these genes was significantly correlated with DNA copy number status, supporting the relevance of focal aberrations. Loss of 5q34 and gain of 13q22.1 were identified as independent prognostic factors of survival in this series of patients. In conclusion, focal chromosomal copy number aberrations in stage II colon cancer are enriched in cancer genes that contribute to and drive the process of colorectal cancer development. DNA copy number status of these genes correlates with mRNA expression and some are associated with clinical outcome.

Co-chaperones are limiting in a depleted chaperone network.

To probe the limiting nodes in the chaperoning network which maintains cellular proteostasis, we expressed a dominant negative mutant of heat shock factor 1 (dnHSF1), the regulator of the cytoplasmic proteotoxic stress response. Microarray analysis of non-stressed dnHSF1 cells showed a two- or more fold decrease in the transcript level of 10 genes, amongst which are the (co-)chaperone genes HSP90AA1, HSPA6, DNAJB1 and HSPB1. Glucocorticoid signaling, which requires the Hsp70 and the Hsp90 folding machines, was severely impaired by dnHSF1, but fully rescued by expression of DNAJA1 or DNAJB1, and partially by ST13. Expression of DNAJB6, DNAJB8, HSPA1A, HSPB1, HSPB8, or STIP1 had no effect while HSP90AA1 even inhibited. PTGES3 (p23) inhibited only in control cells. Our results suggest that the DNAJ co-chaperones in particular become limiting in a depleted chaperoning network. Our results also suggest a difference between the transcriptomes of cells lacking HSF1 and cells expressing dnHSF1.

Gene expression profiling in autoantibody-positive patients with arthralgia predicts development of arthritis.

OBJECTIVE: To identify molecular features associated with the development of rheumatoid arthritis (RA), to understand the pathophysiology of preclinical development of RA, and to assign predictive biomarkers. METHODS: The study group comprised 109 anti-citrullinated protein antibody (ACPA)- and/or rheumatoid factor-positive patients with arthralgia who did not have arthritis but were at risk of RA, and 25 patients with RA. The gene expression profiles of blood samples obtained from these patients were determined by DNA microarray analysis and quantitative polymerase chain reaction. RESULTS: In 20 of the 109 patients with arthralgia who were at risk of RA, arthritis developed after a median of 7 months. Gene expression profiling of blood cells revealed heterogeneity among the at-risk patients, based on differential expression of immune-related genes. This report is the first to describe gene signatures relevant to the development of arthritis. Signatures significantly associated with arthritis development were involved in interferon (IFN)-mediated immunity, hematopoiesis, and chemokine/cytokine activity. Logistic regression analysis revealed that the odds ratio (OR) for developing arthritis within 12 months was 21.0 (95% confidence interval [95% CI] 2.8-156.1 [P = 0.003]) for the subgroup characterized by increased expression of genes involved in IFN-mediated immunity and/or cytokine/chemokine-activity. Genes involved in B cell immunology were associated with protection against progression to arthritis (OR 0.38, 95% CI 0.21-0.70 [P = 0.002]). These processes were reminiscent of those in patients with RA, implying that the preclinical phase of disease is associated with features of established disease. CONCLUSION: The results of this study indicate that IFN-mediated immunity, hematopoiesis, and cell trafficking specify processes relevant to the progression of arthritis independent of ACPA positivity. These findings strongly suggest that certain gene signatures have value for predicting the progression to arthritis, which will pave the way to preventive medicine.

Regulation of IFN response gene activity during infliximab treatment in rheumatoid arthritis is associated with clinical response to treatment.

INTRODUCTION: Cross-regulation between TNF and type I IFN has been postulated to play an important role in autoimmune diseases. Therefore, we determined the effect of TNF blockade in rheumatoid arthritis (RA) on the type I IFN response gene activity in relation to clinical response. METHODS: Peripheral blood from 33 RA patients was collected in PAXgene tubes before and after the start of infliximab treatment. In a first group of 15 patients the baseline expression of type I IFN-regulated genes was determined using cDNA microarrays and compared to levels one month after treatment. The remaining 18 patients were studied as an independent group for validation using quantitative polymerase chain reaction (qPCR). RESULTS: Gene expression analysis revealed that anti-TNF antibody treatment induced a significant increase in type I IFN response gene activity in a subset of RA patients, whereas expression levels remained similar or were slightly decreased in others. The findings appear clinically relevant since patients with an increased IFN response gene activity after anti-TNF therapy had a poor clinical outcome. This association was confirmed and extended for an IFN response gene set consisting of OAS1, LGALS3BP, Mx2, OAS2 and SERPING1 in five EULAR good and five EULAR poor responders, by qPCR. CONCLUSIONS: Regulation of IFN response gene activity upon TNF blockade in RA is not as consistent as previously described, but varies between patients. The differential changes in IFN response gene activity appear relevant to the clinical outcome of TNF blockade in RA.

Gene expression profiling in rheumatology.

In the last decade, the analysis of gene expression in tissues and cells has evolved from the analysis of a selected set of genes to an efficient high throughput whole-genome screening approach of potentially all genes expressed. Development of sophisticated methodologies such as microarray technology allows an open-ended survey to identify comprehensively the fraction of genes that are differentially expressed between samples and that define the samples' unique biology. By a global analysis of the genes that are expressed in cells and tissues of an individual under different conditions and during disease, we can build up "gene expression profiles (signatures)" which characterize the dynamic functioning of the genome under pathophysiological conditions. This strategy also provides the means to subdivide patients that suffer from a complex heterogeneous disease into more homogeneous subgroups. Such discovery-based research identifies biological processes that may include new genes with unknown function or genes not previously known to be involved in this process. The latter category may hold surprises that sometimes urge us to redirect our thinking. We have used microarrays to disclose the heterogeneity of rheumatoid arthritis (RA) patients at the level of gene expression of the affected synovial tissues. Analysis of the expression profiles of synovial tissues from different patients with RA revealed considerable variability, resulting in the identification of at least two molecularly distinct forms of RA tissues. One is characterized by genes that indicate an active inflammatory infiltrate with high immunoglobulin production, whereas the other type shows little immune activation and instead shows a higher stromal cell activity. These results confirm the heterogeneous nature of RA and suggest the existence of distinct pathogenic mechanisms that contribute to RA. The differences in expression profiles provide opportunities to stratify patients for intervention therapies based on molecular criteria.