Benchmarking blood collection tubes and processing intervals for extracellular vesicle performance metrics.

The analysis of extracellular vesicles (EV) in blood samples is under intense investigation and holds the potential to deliver clinically meaningful biomarkers for health and disease. Technical variation must be minimized to confidently assess EV-associated biomarkers, but the impact of pre-analytics on EV characteristics in blood samples remains minimally explored. We present the results from the first large-scale EV Blood Benchmarking (EVBB) study in which we systematically compared 11 blood collection tubes (BCT; six preservation and five non-preservation) and three blood processing intervals (BPI; 1, 8 and 72 h) on defined performance metrics (n = 9). The EVBB study identifies a significant impact of multiple BCT and BPI on a diverse set of metrics reflecting blood sample quality, ex-vivo generation of blood-cell derived EV, EV recovery and EV-associated molecular signatures. The results assist the informed selection of the optimal BCT and BPI for EV analysis. The proposed metrics serve as a framework to guide future research on pre-analytics and further support methodological standardization of EV studies.

mRNA Capture Sequencing and RT-qPCR for the Detection of Pathognomonic, Novel, and Secondary Fusion Transcripts in FFPE Tissue: A Sarcoma Showcase.

We assess the performance of mRNA capture sequencing to identify fusion transcripts in FFPE tissue of different sarcoma types, followed by RT-qPCR confirmation. To validate our workflow, six positive control tumors with a specific chromosomal rearrangement were analyzed using the TruSight RNA Pan-Cancer Panel. Fusion transcript calling by FusionCatcher confirmed these aberrations and enabled the identification of both fusion gene partners and breakpoints. Next, whole-transcriptome TruSeq RNA Exome sequencing was applied to 17 fusion gene-negative alveolar rhabdomyosarcoma (ARMS) or undifferentiated round cell sarcoma (URCS) tumors, for whom fluorescence in situ hybridization (FISH) did not identify the classical pathognomonic rearrangements. For six patients, a pathognomonic fusion transcript was readily detected, i.e., PAX3-FOXO1 in two ARMS patients, and EWSR1-FLI1, EWSR1-ERG, or EWSR1-NFATC2 in four URCS patients. For the 11 remaining patients, 11 newly identified fusion transcripts were confirmed by RT-qPCR, including COPS3-TOM1L2, NCOA1-DTNB, WWTR1-LINC01986, PLAA-MOB3B, AP1B1-CHEK2, and BRD4-LEUTX fusion transcripts in ARMS patients. Additionally, recurrently detected secondary fusion transcripts in patients diagnosed with EWSR1-NFATC2-positive sarcoma were confirmed (COPS4-TBC1D9, PICALM-SYTL2, SMG6-VPS53, and UBE2F-ALS2). In conclusion, this study shows that mRNA capture sequencing enhances the detection rate of pathognomonic fusions and enables the identification of novel and secondary fusion transcripts in sarcomas.

Toward best practice in cancer mutation detection with whole-genome and whole-exome sequencing.

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

Establishing community reference samples, data and call sets for benchmarking cancer mutation detection using whole-genome sequencing.

The lack of samples for generating standardized DNA datasets for setting up a sequencing pipeline or benchmarking the performance of different algorithms limits the implementation and uptake of cancer genomics. Here, we describe reference call sets obtained from paired tumor-normal genomic DNA (gDNA) samples derived from a breast cancer cell line-which is highly heterogeneous, with an aneuploid genome, and enriched in somatic alterations-and a matched lymphoblastoid cell line. We partially validated both somatic mutations and germline variants in these call sets via whole-exome sequencing (WES) with different sequencing platforms and targeted sequencing with >2,000-fold coverage, spanning 82% of genomic regions with high confidence. Although the gDNA reference samples are not representative of primary cancer cells from a clinical sample, when setting up a sequencing pipeline, they not only minimize potential biases from technologies, assays and informatics but also provide a unique resource for benchmarking 'tumor-only' or 'matched tumor-normal' analyses.

Charting Extracellular Transcriptomes in The Human Biofluid RNA Atlas.

Extracellular RNAs present in biofluids have emerged as potential biomarkers for disease. Where most studies focus on blood-derived fluids, other biofluids may be more informative. We present an atlas of messenger, circular, and small RNA transcriptomes of a comprehensive collection of 20 human biofluids. By means of synthetic spike-in controls, we compare RNA content across biofluids, revealing a 10,000-fold difference in concentration. The circular RNA fraction is increased in most biofluids compared to tissues. Each biofluid transcriptome is enriched for RNA molecules derived from specific tissues and cell types. Our atlas enables an informed selection of the most relevant biofluid to monitor particular diseases. To verify the biomarker potential in these biofluids, four validation cohorts representing a broad spectrum of diseases were profiled, revealing numerous differential RNAs between case and control subjects. Spike-normalized data are publicly available in the R2 web portal for further exploration.

Integrated RNA and metabolite profiling of urine liquid biopsies for prostate cancer biomarker discovery.

Sensitive and specific diagnostic and prognostic biomarkers for prostate cancer (PCa) are urgently needed. Urine samples are a non-invasive means to obtain abundant and readily accessible "liquid biopsies". Herein we used urine liquid biopsies to identify and characterize a novel group of urine-enriched RNAs and metabolites in patients with PCa and normal individuals with or without benign prostatic disease. Differentially expressed RNAs were identified in urine samples by deep sequencing and metabolites in urine were measured by mass spectrometry. mRNA and metabolite profiles were distinct in patients with benign and malignant disease. Integrated analysis of urinary gene expression and metabolite signatures unveiled an aberrant glutamate metabolism and tricarboxylic acid (TCA) cycle node in prostate cancer-derived cells. Functional validation supported a role for glutamate metabolism and glutamate oxaloacetate transaminase 1 (GOT1)-dependent redox balance in PCa, which could be exploited for novel biomarkers and therapies. In this study, we discovered cancer-specific changes in urinary RNAs and metabolites, paving the way for the development of sensitive and specific urinary PCa diagnostic biomarkers either alone or in combination. Our methodology was based on single void urine samples (i.e., without prostatic massage). The integrated analysis of metabolomic and transcriptomic data from these liquid biopsies revealed a glutamate metabolism and tricarboxylic acid cycle node that was specific to prostate-derived cancer cells and cancer-specific metabolic changes in urine.

Author Correction: Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

RNA sequencing of esophageal adenocarcinomas identifies novel fusion transcripts, including NPC1-MELK, arising from a complex chromosomal rearrangement.

BACKGROUND: Studies of chromosomal rearrangements and fusion transcripts have elucidated mechanisms of tumorigenesis and led to targeted cancer therapies. This study was aimed at identifying novel fusion transcripts in esophageal adenocarcinoma (EAC). METHODS: To identify new fusion transcripts associated with EAC, targeted RNA sequencing and polymerase chain reaction (PCR) verification were performed in 40 EACs and matched nonmalignant specimens from the same patients. Genomic PCR and Sanger sequencing were performed to find the breakpoint of fusion genes. RESULTS: Five novel in-frame fusion transcripts were identified and verified in 40 EACs and in a validation cohort of 15 additional EACs (55 patients in all): fibroblast growth factor receptor 2 (FGFR2)-GRB2-associated binding protein 2 (GAB2) in 2 of 55 or 3.6%, Niemann-Pick C1 (NPC1)-maternal embryonic leucine zipper kinase (MELK) in 2 of 55 or 3.6%, ubiquitin-specific peptidase 54 (USP54)-calcium/calmodulin dependent protein kinase II γ (CAMK2G) in 2 of 55 or 3.6%, megakaryoblastic leukemia (translocation) 1 (MKL1)-fibulin 1 (FBLN1) in 1 of 55 or 1.8%, and CCR4-NOT transcription complex subunit 2 (CNOT2)-chromosome 12 open reading frame 49 (C12orf49) in 1 of 55 or 1.8%. A genomic analysis indicated that NPC1-MELK arose from a complex interchromosomal translocation event involving chromosomes 18, 3, and 9 with 3 rearrangement points, and this was consistent with chromoplexy. CONCLUSIONS: These data indicate that fusion transcripts occur at a stable frequency in EAC. Furthermore, our results indicate that chromoplexy is an underlying mechanism that generates fusion transcripts in EAC. These and other fusion transcripts merit further study as diagnostic markers and potential therapeutic targets in EAC. Cancer 2017;123:3916-24. © 2017 American Cancer Society.

Complete genome sequence of a KI polyomavirus isolated from an otherwise healthy child with severe lower respiratory tract infection.

Unbiased, deep sequencing of a nasal specimen from an otherwise healthy 13-month-old boy hospitalized in intensive care revealed high gene expression and the complete genome of a novel isolate of KI polyomavirus (KIPyV). Further investigation detected minimal gene expression of additional viruses, suggesting that KIPyV was potentially the causal agent. Analysis of the complete genome of isolate NMKI001 revealed it is different from all previously reported genomes and contains two amino acid differences as compared to the closest virus isolate, Stockholm 380 (EF127908). J. Med. Virol. 89:926-930, 2017. © 2016 Wiley Periodicals, Inc.

Cloche is a bHLH-PAS transcription factor that drives haemato-vascular specification.

Vascular and haematopoietic cells organize into specialized tissues during early embryogenesis to supply essential nutrients to all organs and thus play critical roles in development and disease. At the top of the haemato-vascular specification cascade lies cloche, a gene that when mutated in zebrafish leads to the striking phenotype of loss of most endothelial and haematopoietic cells and a significant increase in cardiomyocyte numbers. Although this mutant has been analysed extensively to investigate mesoderm diversification and differentiation and continues to be broadly used as a unique avascular model, the isolation of the cloche gene has been challenging due to its telomeric location. Here we used a deletion allele of cloche to identify several new cloche candidate genes within this genomic region, and systematically genome-edited each candidate. Through this comprehensive interrogation, we succeeded in isolating the cloche gene and discovered that it encodes a PAS-domain-containing bHLH transcription factor, and that it is expressed in a highly specific spatiotemporal pattern starting during late gastrulation. Gain-of-function experiments show that it can potently induce endothelial gene expression. Epistasis experiments reveal that it functions upstream of etv2 and tal1, the earliest expressed endothelial and haematopoietic transcription factor genes identified to date. A mammalian cloche orthologue can also rescue blood vessel formation in zebrafish cloche mutants, indicating a highly conserved role in vertebrate vasculogenesis and haematopoiesis. The identification of this master regulator of endothelial and haematopoietic fate enhances our understanding of early mesoderm diversification and may lead to improved protocols for the generation of endothelial and haematopoietic cells in vivo and in vitro.

Complete Genome Sequence of a Novel Human WU Polyomavirus Isolate Associated with Acute Respiratory Infection.

We report here the complete genome sequence of a WU polyomavirus (WUPyV) isolate, NM040708, collected from a patient with an acute respiratory infection in New Mexico. The double-stranded DNA (dsDNA) genome of NM040708 is 5,229 bp in length and differs from the WUPyV reference with accession no. NC_009539 by 6 nucleotides and 2 amino acids.

LIN28A Modulates Splicing and Gene Expression Programs in Breast Cancer Cells.

LIN28 is an evolutionarily conserved RNA-binding protein with critical functions in developmental timing and cancer. However, the molecular mechanisms underlying LIN28's oncogenic properties are yet to be described. RNA-protein immunoprecipitation coupled with genome-wide sequencing (RIP-Seq) analysis revealed significant LIN28 binding within 843 mRNAs in breast cancer cells. Many of the LIN28-bound mRNAs are implicated in the regulation of RNA and cell metabolism. We identify heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1), a protein with multiple roles in mRNA metabolism, as a LIN28-interacting partner. Subsequently, we used a custom computational method to identify differentially spliced gene isoforms in LIN28 and hnRNP A1 small interfering RNA (siRNA)-treated cells. The results reveal that these proteins regulate alternative splicing and steady-state mRNA expression of genes implicated in aspects of breast cancer biology. Notably, cells lacking LIN28 undergo significant isoform switching of the ENAH gene, resulting in a decrease in the expression of the ENAH exon 11a isoform. The expression of ENAH isoform 11a has been shown to be elevated in breast cancers that express HER2. Intriguingly, analysis of publicly available array data from the Cancer Genome Atlas (TCGA) reveals that LIN28 expression in the HER2 subtype is significantly different from that in other breast cancer subtypes. Collectively, our data suggest that LIN28 may regulate splicing and gene expression programs that drive breast cancer subtype phenotypes.

Exome-level comparison of primary well-differentiated neuroendocrine tumors and their cell lines.

Neuroendocrine cancer cell lines are used to investigate therapeutic targets in neuroendocrine tumors (NET) and have been instrumental in the design of clinical trials targeting the PI3K/AKT/mTOR pathways, VEGF inhibitors, and somatostatin analogues. It remains unknown, however, whether the genomic makeup of NET cell lines reflect that of primary NET since comprehensive unbiased genome sequencing has not been performed on the cell lines. Four bronchopulmonary NET (BP-NET)-NCI-H720, NCI-H727, NCI-H835, and UMC11-and two pancreatic neuroendocrine tumors (panNET)-BON-1 and QGP1-were cultured. DNA was isolated, and exome sequencing was done. GATK and EXCAVATOR were used for bioinformatic analysis. We detected a total of 1,764 nonsynonymous single nucleotide variants at a rate of 8 per Mb in BP-NET and 4.3 per Mb in panNET cell lines, including 52 mutated COSMIC cancer genes in these cell lines, such as TP53, BRCA1, RB1, TSC2, NOTCH1, EP300, GNAS, KDR, STK11, and APC but not ATRX, DAXX, nor MEN1. Our data suggest that mutation rate, the pattern of copy number variations, and the mutational spectra in the BP-NET cell lines are more similar to the changes observed in small cell lung cancer than those found in primary BP-NET. Likewise, mutation rate and pattern including the absence of mutations in ATRX/DAXX, MEN1, and YY1 in the panNET cell lines BON1 and QGP1 suggest that these cell lines do not have the genetic signatures of a primary panNET. These results suggest that results from experiments with BP-NET and panNET cell lines need to be interpreted with caution.

Sequencing of Charcot-Marie-Tooth disease genes in a toxic polyneuropathy.

OBJECTIVE: Mutations in Charcot-Marie-Tooth disease (CMT) genes are the cause of rare familial forms of polyneuropathy. Whether allelic variability in CMT genes is also associated with common forms of polyneuropathy-considered "acquired" in medical parlance-is unknown. Chemotherapy-induced peripheral neuropathy (CIPN) occurs commonly in cancer patients and is individually unpredictable. We used CIPN as a clinical model to investigate the association of non-CMT polyneuropathy with CMT genes. METHODS: A total of 269 neurologically asymptomatic cancer patients were enrolled in the clinical trial Alliance N08C1 to receive the neurotoxic drug paclitaxel, while undergoing prospective assessments for polyneuropathy. Forty-nine CMT genes were analyzed by targeted massively parallel sequencing of genomic DNA from patient blood. RESULTS: A total of 119 (of 269) patients were identified from the 2 ends of the polyneuropathy phenotype distribution: patients that were most and least susceptible to paclitaxel polyneuropathy. The CMT gene PRX was found to be deleteriously mutated in patients who were susceptible to CIPN but not in controls (p = 8 × 10(-3)). Genetic variation in another CMT gene, ARHGEF10, was highly significantly associated with CIPN (p = 5 × 10(-4)). Three nonsynonymous recurrent single nucleotide variants contributed to the ARHGEF10 signal: rs9657362, rs2294039, and rs17683288. Of these, rs9657362 had the strongest effect (odds ratio = 4.8, p = 4 × 10(-4)). INTERPRETATION: The results reveal an association of CMT gene allelic variability with susceptibility to CIPN. The findings raise the possibility that other acquired polyneuropathies may also be codetermined by genetic etiological factors, of which some may be related to genes already known to cause the phenotypically related Mendelian disorders of CMT.

From single-cell to cell-pool transcriptomes: stochasticity in gene expression and RNA splicing.

Single-cell RNA-seq mammalian transcriptome studies are at an early stage in uncovering cell-to-cell variation in gene expression, transcript processing and editing, and regulatory module activity. Despite great progress recently, substantial challenges remain, including discriminating biological variation from technical noise. Here we apply the SMART-seq single-cell RNA-seq protocol to study the reference lymphoblastoid cell line GM12878. By using spike-in quantification standards, we estimate the absolute number of RNA molecules per cell for each gene and find significant variation in total mRNA content: between 50,000 and 300,000 transcripts per cell. We directly measure technical stochasticity by a pool/split design and find that there are significant differences in expression between individual cells, over and above technical variation. Specific gene coexpression modules were preferentially expressed in subsets of individual cells, including one enriched for mRNA processing and splicing factors. We assess cell-to-cell variation in alternative splicing and allelic bias and report evidence of significant differences in splice site usage that exceed splice variation in the pool/split comparison. Finally, we show that transcriptomes from small pools of 30-100 cells approach the information content and reproducibility of contemporary RNA-seq from large amounts of input material. Together, our results define an experimental and computational path forward for analyzing gene expression in rare cell types and cell states.

mRNA-Seq of single prostate cancer circulating tumor cells reveals recapitulation of gene expression and pathways found in prostate cancer.

Circulating tumor cells (CTC) mediate metastatic spread of many solid tumors and enumeration of CTCs is currently used as a prognostic indicator of survival in metastatic prostate cancer patients. Some evidence suggests that it is possible to derive additional information about tumors from expression analysis of CTCs, but the technical difficulty of isolating and analyzing individual CTCs has limited progress in this area. To assess the ability of a new generation of MagSweeper to isolate intact CTCs for downstream analysis, we performed mRNA-Seq on single CTCs isolated from the blood of patients with metastatic prostate cancer and on single prostate cancer cell line LNCaP cells spiked into the blood of healthy donors. We found that the MagSweeper effectively isolated CTCs with a capture efficiency that matched the CellSearch platform. However, unlike CellSearch, the MagSweeper facilitates isolation of individual live CTCs without contaminating leukocytes. Importantly, mRNA-Seq analysis showed that the MagSweeper isolation process did not have a discernible impact on the transcriptional profile of single LNCaPs isolated from spiked human blood, suggesting that any perturbations caused by the MagSweeper process on the transcriptional signature of isolated cells are modest. Although the RNA from patient CTCs showed signs of significant degradation, consistent with reports of short half-lives and apoptosis amongst CTCs, transcriptional signatures of prostate tissue and of cancer were readily detectable with single CTC mRNA-Seq. These results demonstrate that the MagSweeper provides access to intact CTCs and that these CTCs can potentially supply clinically relevant information.

Highly sensitive profiling of CD44+/CD24- breast cancer stem cells by combining global mRNA amplification and next generation sequencing: evidence for a hyperactive PI3K pathway.

We performed next generation sequencing- and microarray-based gene expression profiling of CD44(+)/CD24(-)/CD45(-) breast CSCs (cancer stem cells) isolated from primary ERα-positive breast cancer. By combining semi-automated dissociation of human tumor tissue, magnetic cell sorting and cDNA amplification less than 500 CSCs were required for transcriptome analyses. Besides overexpressing genes involved in maintenance of stemness, the CSCs showed higher levels of genes that drive the PI3K pathway, including EGFR, HB-EGF, PDGFRA/B, PDGF, MET, PIK3CA, PIK3R1 and PIK3R2. This suggests that, in CSCs of ERα-positive breast cancer, the PI3K pathway which is involved in endocrine resistance is hyperactive.

Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.

Genome-wide transcriptome analyses are routinely used to monitor tissue-, disease- and cell type­specific gene expression, but it has been technically challenging to generate expression profiles from single cells. Here we describe a robust mRNA-Seq protocol (Smart-Seq) that is applicable down to single cell levels. Compared with existing methods, Smart-Seq has improved read coverage across transcripts, which enhances detailed analyses of alternative transcript isoforms and identification of single-nucleotide polymorphisms. We determined the sensitivity and quantitative accuracy of Smart-Seq for single-cell transcriptomics by evaluating it on total RNA dilution series. We found that although gene expression estimates from single cells have increased noise, hundreds of differentially expressed genes could be identified using few cells per cell type. Applying Smart-Seq to circulating tumor cells from melanomas, we identified distinct gene expression patterns, including candidate biomarkers for melanoma circulating tumor cells. Our protocol will be useful for addressing fundamental biological problems requiring genome-wide transcriptome profiling in rare cells.

Base-pair resolution DNA methylation sequencing reveals profoundly divergent epigenetic landscapes in acute myeloid leukemia.

We have developed an enhanced form of reduced representation bisulfite sequencing with extended genomic coverage, which resulted in greater capture of DNA methylation information of regions lying outside of traditional CpG islands. Applying this method to primary human bone marrow specimens from patients with Acute Myelogeneous Leukemia (AML), we demonstrated that genetically distinct AML subtypes display diametrically opposed DNA methylation patterns. As compared to normal controls, we observed widespread hypermethylation in IDH mutant AMLs, preferentially targeting promoter regions and CpG islands neighboring the transcription start sites of genes. In contrast, AMLs harboring translocations affecting the MLL gene displayed extensive loss of methylation of an almost mutually exclusive set of CpGs, which instead affected introns and distal intergenic CpG islands and shores. When analyzed in conjunction with gene expression profiles, it became apparent that these specific patterns of DNA methylation result in differing roles in gene expression regulation. However, despite this subtype-specific DNA methylation patterning, a much smaller set of CpG sites are consistently affected in both AML subtypes. Most CpG sites in this common core of aberrantly methylated CpGs were hypermethylated in both AML subtypes. Therefore, aberrant DNA methylation patterns in AML do not occur in a stereotypical manner but rather are highly specific and associated with specific driving genetic lesions.

Detection of redundant fusion transcripts as biomarkers or disease-specific therapeutic targets in breast cancer.

Fusion genes and fusion gene products are widely employed as biomarkers and therapeutic targets in hematopoietic cancers, but their applications have yet to be appreciated in solid tumors. Here, we report the use of SnowShoes-FTD, a powerful new analytic pipeline that can identify fusion transcripts and assess their redundancy and tumor subtype-specific distribution in primary tumors. In a study of primary breast tumors, SnowShoes-FTD was used to analyze paired-end mRNA-Seq data from a panel of estrogen receptor (ER)(+), HER2(+), and triple-negative primary breast tumors, identifying tumor-specific fusion transcripts by comparison with mRNA-Seq data from nontransformed human mammary epithelial cell cultures plus the Illumina Body Map data from normal tissues. We found that every primary breast tumor that was analyzed expressed one or more fusion transcripts. Of the 131 tumor-specific fusion transcripts identified, 86 were "private" (restricted to a single tumor) and 45 were "redundant" (distributed among multiple tumors). Among the redundant fusion transcripts, 7 were unique to ER(+) tumors and 8 were unique to triple-negative tumors. In contrast, none of the redundant fusion transcripts were unique to HER2(+) tumors. Both private and redundant fusion transcripts were widely expressed in primary breast tumors, with many mapping to genomic loci implicated in breast carcinogenesis and/or risk. Our finding that some fusion transcripts are tumor subtype-specific suggests that these entities may be critical determinants in the etiology of breast cancer subtypes, useful as biomarkers for tumor stratification, or exploitable as cancer-specific therapeutic targets.