Non-random genetic alterations in the cyanobacterium Nostoc sp. exposed to space conditions.

Understanding the impact of long-term exposure of microorganisms to space is critical in understanding how these exposures impact the evolution and adaptation of microbial life under space conditions. In this work we subjected Nostoc sp. CCCryo 231-06, a cyanobacterium capable of living under many different ecological conditions, and also surviving in extreme ones, to a 23-month stay at the International Space Station (the Biology and Mars Experiment, BIOMEX, on the EXPOSE-R2 platform) and returned it to Earth for single-cell genome analysis. We used microfluidic technology and single cell sequencing to identify the changes that occurred in the whole genome of single Nostoc cells. The variant profile showed that biofilm and photosystem associated loci were the most altered, with an increased variant rate of synonymous base pair substitutions. The cause(s) of these non-random alterations and their implications to the evolutionary potential of single bacterial cells under long-term cosmic exposure warrants further investigation.

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.

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γδß)0-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 (δß)0 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.

LPCAT1-TERT fusions are uniquely recurrent in epithelioid trophoblastic tumors and positively regulate cell growth.

Gestational trophoblastic disease (GTD) is a heterogeneous group of lesions arising from placental tissue. Epithelioid trophoblastic tumor (ETT), derived from chorionic-type trophoblast, is the rarest form of GTD with only approximately 130 cases described in the literature. Due to its morphologic mimicry of epithelioid smooth muscle tumors and carcinoma, ETT can be misdiagnosed. To date, molecular characterization of ETTs is lacking. Furthermore, ETT is difficult to treat when disease spreads beyond the uterus. Here using RNA-Seq analysis in a cohort of ETTs and other gestational trophoblastic lesions we describe the discovery of LPCAT1-TERT fusion transcripts that occur in ETTs and coincide with underlying genomic deletions. Through cell-growth assays we demonstrate that LPCAT1-TERT fusion proteins can positively modulate cell proliferation and therefore may represent future treatment targets. Furthermore, we demonstrate that TERT upregulation appears to be a characteristic of ETTs, even in the absence of LPCAT1-TERT fusions, and that it appears linked to copy number gains of chromosome 5. No evidence of TERT upregulation was identified in other trophoblastic lesions tested, including placental site trophoblastic tumors and placental site nodules, which are thought to be the benign chorionic-type trophoblast counterpart to ETT. These findings indicate that LPCAT1-TERT fusions and copy-number driven TERT activation may represent novel markers for ETT, with the potential to improve the diagnosis, treatment, and outcome for women with this rare form of GTD.

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.

Chromosomal Junction Detection from Whole-Genome Sequencing on Formalin-Fixed, Paraffin-Embedded Tumors.

DNA junctions (DNAJs) frequently impact clinically relevant genes in tumors and are important for diagnostic and therapeutic purposes. Although routinely screened through fluorescence in situ hybridization assays, such testing only allows the interrogation of single-gene regions or known fusion partners. Comprehensive assessment of DNAJs present across the entire genome can only be determined from whole-genome sequencing. Structural variance analysis from whole-genome paired-end sequencing data is, however, frequently restricted to copy number changes without DNAJ detection. Through optimized whole-genome sequencing and specialized bioinformatics algorithms, complete structural variance analysis is reported, including DNAJs, from formalin-fixed DNA. Selective library assembly from larger fragments (>500 bp) and economical sequencing depths (300 to 400 million reads) provide representative genomic coverage profiles and increased allelic coverage to levels compatible with DNAJ calling (40× to 60×). Although consistently fragmented, more recently formalin-fixed, specimens (<2 years' storage) revealed consistent populations of larger DNA fragments. Optimized bioinformatics efficiently detected >90% of DNAJs in two prostate tumors (approximately 60% tumor) previously analyzed by mate-pair sequencing on fresh frozen tissue, with evidence of at least one spanning-read in 99% of DNAJs. Rigorous masking with data from unrelated formalin-fixed tissue progressively eliminated many false-positive DNAJs, without loss of true positives, resulting in low numbers of false-positive passing current filters. This methodology enables more comprehensive clinical genomics testing on formalin-fixed clinical 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.

High frequency electrical stimulation promotes expression of extracellular matrix proteins from human astrocytes.

Therapeutic benefits of deep brain stimulation (DBS), a neurosurgical treatment for certain movement disorders and other neurologic conditions, are well documented, but DBS mechanisms remain largely unexplained. DBS is thought to modulate pathological neural activity. However, although astrocytes, the most numerous cell type in the brain, play a significant role in neurotransmission, chemical homeostasis and synaptic plasticity, their role in DBS has not been fully examined. To investigate astrocytic function in DBS, we applied DBS-like high frequency electrical stimulation for 24 h to human astrocytes in vitro and analyzed single cell transcriptome mRNA profile. We found that DBS-like high frequency stimulation negatively impacts astrocyte metabolism and promotes the release of extracellular matrix (matricellular) proteins, including IGFBP3, GREM1, IGFBP5, THBS1, and PAPPA. Our results suggest that astrocytes are involved in the long-term modulation of extra cellular matrix environments and that they may influence persistent cell-to-cell interaction and help maintain neuromodulation over time.

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.

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-14). 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.

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.

Development and Verification of an RNA Sequencing (RNA-Seq) Assay for the Detection of Gene Fusions in Tumors.

We assessed the performance characteristics of an RNA sequencing (RNA-Seq) assay designed to detect gene fusions in 571 genes to help manage patients with cancer. Polyadenylated RNA was converted to cDNA, which was then used to prepare next-generation sequencing libraries that were sequenced on an Illumina HiSeq 2500 instrument and analyzed with an in-house developed bioinformatic pipeline. The assay identified 38 of 41 gene fusions detected by another method, such as fluorescence in situ hybridization or RT-PCR, for a sensitivity of 93%. No false-positive gene fusions were identified in 15 normal tissue specimens and 10 tumor specimens that were negative for fusions by RNA sequencing or Mate Pair NGS (100% specificity). The assay also identified 22 fusions in 17 tumor specimens that had not been detected by other methods. Eighteen of the 22 fusions had not previously been described. Good intra-assay and interassay reproducibility was observed with complete concordance for the presence or absence of gene fusions in replicates. The analytical sensitivity of the assay was tested by diluting RNA isolated from gene fusion-positive cases with fusion-negative RNA. Gene fusions were generally detectable down to 12.5% dilutions for most fusions and as little as 3% for some fusions. This assay can help identify fusions in patients with cancer; these patients may in turn benefit from both US Food and Drug Administration-approved and investigational targeted therapies.

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.