Onepot-Seq: capturing single-cell transcriptomes simultaneously in a continuous medium via transient localization of mRNA.

The development of single-cell RNA-seq has broadened the spectrum for biological research by providing a high-resolution analysis of cellular heterogeneity. However, the requirement for sophisticated devices for the compartmentalization of cells has limited its widespread applicability. Here, we develop Onepot-Seq, a device-free method, that harnesses the transient localization of mRNA after lysis to capture single-cell transcriptomes simultaneously in a continuous fluid medium. In mixed-species experiments, we obtained high-quality single-cell profiles. Further, cell type-specific poly(A)-conjugated antibodies allow Onepot-Seq to effectively capture target cells in complex populations. Chemical perturbations to cells can be profiled by Onepot-Seq at single-cell resolution. Onepot-Seq should allow routine transcriptional profiling at single-cell resolution, accelerating clinical and scientific discoveries in many fields of science.

This article describes a strategy for single-cell RNA sequencing that avoids the usual imitation required by first physically compartmentalizing single cells. Based on the slow rate of mRNA diffusion relative to mRNA capture on beads, the authors devised a method that profiles single-cell mRNA among a multiple cell population, within one reaction mixture, termed 'Onepot-Seq'. This approach uses a time-constrained, transient reaction chamber that forms around each cell, with only slow diffusion of the cell contents after gentle lysis. The sparse distribution of cells and beads allows the beads to capture the profile of a single cell with minimal intercellular mRNA contamination.

Personalised circulating tumour DNA assay with large-scale mutation coverage for sensitive minimal residual disease detection in colorectal cancer.

BACKGROUND: Postoperative minimal residual disease (MRD) detection using circulating-tumour DNA (ctDNA) requires a highly sensitive analysis platform. We have developed a tumour-informed, hybrid-capture ctDNA sequencing MRD assay. METHODS: Personalised target-capture panels for ctDNA detection were designed using individual variants identified in tumour whole-exome sequencing of each patient. MRD status was determined using ultra-high-depth sequencing data of plasma cell-free DNA. The MRD positivity and its association with clinical outcome were analysed in Stage II or III colorectal cancer (CRC). RESULTS: In 98 CRC patients, personalised panels for ctDNA sequencing were built from tumour data, including a median of 185 variants per patient. In silico simulation showed that increasing the number of target variants increases MRD detection sensitivity in low fractions (<0.01%). At postoperative 3-week, 21.4% of patients were positive for MRD by ctDNA. Postoperative positive MRD was strongly associated with poor disease-free survival (DFS) (adjusted hazard ratio 8.40, 95% confidence interval 3.49-20.2). Patients with a negative conversion of MRD after adjuvant therapy showed significantly better DFS (P < 0.001). CONCLUSION: Tumour-informed, hybrid-capture-based ctDNA assay monitoring a large number of patient-specific mutations is a sensitive strategy for MRD detection to predict recurrence in CRC.

Dynamic changes in longitudinal circulating tumour DNA profile during metastatic colorectal cancer treatment.

BACKGROUND: Circulating tumour DNA (ctDNA) has been spotlighted as an attractive biomarker because of its easy accessibility and real-time representation of tumour genetic profile. However, the clinical utility of longitudinal ctDNA monitoring has not been clearly defined. METHODS: Serial blood samples were obtained from metastatic colorectal cancer patients undergoing first-line chemotherapy. ctDNA was sequenced using a targeted next-generation sequencing platform which included 106 genes. Changes in ctDNA profile and treatment outcome were comprehensively analysed. RESULTS: A total of 272 samples from 62 patients were analysed. In all, 90.3% of patients had detectable ctDNA mutation before treatment. ctDNA clearance after chemotherapy was associated with longer progression-free survival which was independent of radiological response (adjusted hazard ratio 0.22, 95% confidence interval 0.10-0.46). Longitudinal monitoring was able to detect ctDNA progression which preceded radiological progressive disease (PD) in 58.1% (median 3.3 months). Diverse resistant mutations (34.9%) and gene amplification (7.0%) at the time of PD were discovered. For 16.3% of the PD patients, the newly identified mutations could be potential candidates of targeted therapy or clinical trial. CONCLUSION: ctDNA profile provided a more accurate landscape of tumour and dynamic changes compared to radiological evaluation. Longitudinal ctDNA monitoring may improve personalised treatment decision-making.

CRISPR-Cap: multiplexed double-stranded DNA enrichment based on the CRISPR system.

Existing methods to enrich target regions of genomic DNA based on PCR, hybridization capture, or molecular inversion probes have various drawbacks, including long experiment times and low throughput and/or enrichment quality. We developed CRISPR-Cap, a simple and scalable CRISPR-based method to enrich target regions of dsDNA, requiring only two short experimental procedures that can be completed within two hours. We used CRISPR-Cap to enrich 10 target genes 355.7-fold on average from Escherichia coli genomic DNA with a maximum on-target ratio of 81% and high enrichment uniformity. We also used CRISPR-Cap to measure gene copy numbers and detect rare alleles with frequencies as low as 1%. Finally, we enriched coding sequence regions of 20 genes from the human genome. We envision that CRISPR-Cap can be used as an alternative to other widely used target-enrichment methods, which will broaden the scope of CRISPR applications to the field of target enrichment field.

Highly selective retrieval of accurate DNA utilizing a pool of in situ-replicated DNA from multiple next-generation sequencing platforms.

Scalable and cost-effective production of error-free DNA is critical to meet the increased demand for such DNA in the field of biological science. Methods based on 'Dial-out PCR' have enabled the high-throughput error-free DNA synthesis from a microarray-synthesized DNA pool by labeling with retrieval PCR tags, and retrieving error-free DNA of which the sequence is identified via next generation sequencing (NGS). However, most of the retrieved products contain byproducts due to background amplification of redundantly labeled DNAs. Here, we present a highly selective retrieval method of desired DNA from a pool of millions of DNA clones from NGS platforms. Our strategy is based on replicating entire sequence-verified DNA molecules from NGS plates to obtain population-controlled DNA pool. Using the NGS-replica pool, we could perform improved and selective retrieval of desired DNA from the replicated DNA pool compared to other dial-out PCR based methods. To evaluate the method, we tested this strategy by using 454, Illumina, and Ion Torrent platforms for producing NGS-replica pool. As a result, we observed a highly selective retrieval yield of over 95%. We anticipate that applications based on this method will enable the preparation of high-fidelity sequenced DNA from heterogeneous collections of DNA molecules.

High-throughput construction of multiple cas9 gene variants via assembly of high-depth tiled and sequence-verified oligonucleotides.

Selective retrieval of sequence-verified oligonucleotides (oligos) from next-generation sequencing (NGS) flow cells, termed megacloning, promises accurate and reliable gene synthesis. However, gene assembly requires a complete collection of overlapping sense and nonsense oligos, and megacloning does not typically guarantee the complete production of sequence-verified oligos. Therefore, missing oligos must be provided via repetitive rounds of megacloning, which introduces a bottleneck for scaled-up efforts at gene assembly. Here, we introduce the concept of high-depth tiled oligo design to successfully utilize megacloned oligos for gene synthesis. Using acquired oligos from a single round of the megacloning process, we assembled 72 of 81 target Cas9-coding gene variants. We further validated 62 of these cas9 constructs, and deposited the plasmids to Addgene for subsequent functional characterization by the scientific community. This study demonstrates the utility of using sequence-verified oligos for DNA assembly and provides a practical and reliable optimized method for high-throughput gene synthesis.

Association of pathway mutation with survival after recurrence in colorectal cancer patients treated with adjuvant fluoropyrimidine and oxaliplatin chemotherapy.

BACKGROUND: Although the prognostic biomarkers associated with colorectal cancer (CRC) survival are well known, there are limited data on the association between the molecular characteristics and survival after recurrence (SAR). The purpose of this study was to assess the association between pathway mutations and SAR. METHODS: Of the 516 patients with stage III or high risk stage II CRC patients treated with surgery and adjuvant chemotherapy, 87 who had distant recurrence were included in the present study. We analyzed the association between SAR and mutations of 40 genes included in the five critical pathways of CRC (WNT, P53, RTK-RAS, TGF-ß, and PI3K). RESULTS: Mutation of genes within the WNT, P53, RTK-RAS, TGF-ß, and PI3K pathways were shown in 69(79.3%), 60(69.0%), 57(65.5%), 21(24.1%), and 19(21.8%) patients, respectively. Patients with TGF-ß pathway mutation were younger and had higher incidence of mucinous adenocarcinoma (MAC) histology and microsatellite instability-high. TGF-ß pathway mutation (median SAR of 21.6 vs. 44.4 months, p = 0.021) and MAC (20.0 vs. 44.4 months, p = 0.003) were associated with poor SAR, and receiving curative resection after recurrence was associated with favorable SAR (Not reached vs. 23.6 months, p <  0.001). Mutations in genes within other critical pathways were not associated with SAR. When MAC was excluded as a covariate, multivariate analysis revealed TGF-ß pathway mutation and curative resection after distant recurrence as an independent prognostic factor for SAR. The impact of TGF-ß pathway mutations were predicted using the PROVEAN, SIFT, and PolyPhen-2. Among 25 mutations, 23(92.0%)-24(96.0%) mutations were predicted to be damaging mutation. CONCLUSIONS: Mutation in genes within TGF-ß pathway may have negative prognostic role for SAR in CRC. Other pathway mutations were not associated with SAR.

Association between mutations of critical pathway genes and survival outcomes according to the tumor location in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) develops through the alteration of several critical pathways. This study was aimed at evaluating the influence of critical pathways on survival outcomes for patients with CRC. METHODS: Targeted next-generation sequencing of 40 genes included in the 5 critical pathways of CRC (WNT, P53, RTK-RAS, phosphatidylinositol-4,5-bisphosphate 3-kinase [PI3K], and transforming growth factor ß [TGF-ß]) was performed for 516 patients with stage III or high-risk stage II CRC treated with surgery followed by adjuvant fluoropyrimidine and oxaliplatin chemotherapy. The associations between critical pathway mutations and relapse-free survival (RFS) and overall survival were analyzed. The associations were further analyzed according to the tumor location. RESULTS: The mutation rates for the WNT, P53, RTK-RAS, PI3K, and TGF-ß pathways were 84.5%, 69.0%, 60.7%, 30.0%, and 28.9%, respectively. A mutation in the PI3K pathway was associated with longer RFS (adjusted hazard ratio [HR], 0.59; 95% confidence interval [CI], 0.36-0.99), whereas a mutation in the RTK-RAS pathway was associated with shorter RFS (adjusted HR, 1.60; 95% CI, 1.01-2.52). Proximal tumors showed a higher mutation rate than distal tumors, and the mutation profile was different according to the tumor location. The mutation rates of Kirsten rat sarcoma viral oncogene homolog (KRAS), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit α (PIK3CA), and B-Raf proto-oncogene serine/threonine kinase (BRAF) were higher in proximal tumors, and the mutation rates of adenomatous polyposis coli (APC), tumor protein 53 (TP53), and neuroblastoma RAS viral oncogene homolog (NRAS) were higher in distal tumors. The better RFS with the PI3K pathway mutation was significant only for proximal tumors, and the worse RFS with the RTK-RAS pathway mutation was significant only for distal tumors. CONCLUSIONS: A PI3K pathway mutation was associated with better RFS for CRC patients treated with adjuvant chemotherapy, and an RTK-RAS pathway mutation was associated with worse RFS. The significance of the prognostic impact differed according to the tumor location. Cancer 2017;123:3513-23. © 2017 American Cancer Society.

microDuMIP: target-enrichment technique for microarray-based duplex molecular inversion probes.

Molecular inversion probe (MIP)-based capture is a scalable and effective target-enrichment technology that can use synthetic single-stranded oligonucleotides as probes. Unlike the straightforward use of synthetic oligonucleotides for low-throughput target capture, high-throughput MIP capture has required laborious protocols to generate thousands of single-stranded probes from DNA microarray because of multiple enzymatic steps, gel purifications and extensive PCR amplifications. Here, we developed a simple and efficient microarray-based MIP preparation protocol using only one enzyme with double-stranded probes and improved target capture yields by designing probes with overlapping targets and unique barcodes. To test our strategy, we produced 11 510 microarray-based duplex MIPs (microDuMIPs) and captured 3554 exons of 228 genes in a HapMap genomic DNA sample (NA12878). Under our protocol, capture performance and precision of calling were compatible to conventional MIP capture methods, yet overlapping targets and unique barcodes allowed us to precisely genotype with as little as 50 ng of input genomic DNA without library preparation. microDuMIP method is simpler and cheaper, allowing broader applications and accurate target sequencing with a scalable number of targets.

Genome-scale promoter engineering by coselection MAGE.

Multiplex automated genome engineering (MAGE) uses short oligonucleotides to scarlessly modify genomes; however, insertions >10 bases are still inefficient but can be improved substantially by selection of highly modified chromosomes. Here we describe 'coselection' MAGE (CoS-MAGE) to optimize biosynthesis of aromatic amino acid derivatives by combinatorially inserting multiple T7 promoters simultaneously into 12 genomic operons. Promoter libraries can be quickly generated to study gain-of-function epistatic interactions in gene networks.

Tumor evolution and intratumor heterogeneity of an epithelial ovarian cancer investigated using next-generation sequencing.

BACKGROUND: The extent to which metastatic tumors further evolve by accumulating additional mutations is unclear and has yet to be addressed extensively using next-generation sequencing of high-grade serous ovarian cancer. METHODS: Eleven spatially separated tumor samples from the primary tumor and associated metastatic sites and two normal samples were obtained from a Stage IIIC ovarian cancer patient during cytoreductive surgery prior to chemotherapy. Whole exome sequencing and copy number analysis were performed. Omental exomes were sequenced with a high depth of coverage to thoroughly explore the variants in metastatic lesions. Somatic mutations were further validated by ultra-deep targeted sequencing to sort out false positives and false negatives. Based on the somatic mutations and copy number variation profiles, a phylogenetic tree was generated to explore the evolutionary relationship among tumor samples. RESULTS: Only 6% of the somatic mutations were present in every sample of a given case with TP53 as the only known mutant gene consistently present in all samples. Two non-spatial clusters of primary tumors (cluster P1 and P2), and a cluster of metastatic regions (cluster M) were identified. The patterns of mutations indicate that cluster P1 and P2 diverged in the early phase of tumorigenesis, and that metastatic cluster M originated from the common ancestral clone of cluster P1 with few somatic mutations and copy number variations. CONCLUSIONS: Although a high level of intratumor heterogeneity was evident in high-grade serous ovarian cancer, our results suggest that transcoelomic metastasis arises with little accumulation of somatic mutations and copy number alterations in this patient.

'Shotgun DNA synthesis' for the high-throughput construction of large DNA molecules.

We developed a highly scalable 'shotgun' DNA synthesis technology by utilizing microchip oligonucleotides, shotgun assembly and next-generation sequencing technology. A pool of microchip oligonucleotides targeting a penicillin biosynthetic gene cluster were assembled into numerous random fragments, and tagged with 20 bp degenerate barcode primer pairs. An optimal set of error-free fragments were identified by high-throughput DNA sequencing, selectively amplified using the barcode sequences, and successfully assembled into the target gene cluster.

An open-source, 3D printed inkjet DNA synthesizer.

Synthetic oligonucleotides have become a fundamental tool in a wide range of biological fields, including synthetic biology, biosensing, and DNA storage. Reliable access to equipment for synthesizing high-density oligonucleotides in the laboratory ensures research security and the freedom of research expansion. In this study, we introduced the Open-Source Inkjet DNA Synthesizer (OpenIDS), an open-source inkjet-based microarray synthesizer that offers ease of construction, rapid deployment, and flexible scalability. Utilizing 3D printing, Arduino, and Raspberry Pi, this newly designed synthesizer achieved robust stability with an industrial inkjet printhead. OpenIDS maintains low production costs and is therefore suitable for self-fabrication and optimization in academic laboratories. Moreover, even non-experts can create and control the synthesizer with a high degree of freedom for structural modifications. Users can easily add printheads or alter the design of the microarray substrate according to their research needs. To validate its performance, we synthesized oligonucleotides on 144 spots on a 15 × 25-mm silicon wafer filled with controlled pore glass. The synthesized oligonucleotides were analyzed using urea polyacrylamide gel electrophoresis.

Cost and time-efficient construction of a 3'-end mRNA library from unpurified bulk RNA in a single tube.

The major drawbacks of RNA sequencing (RNA-seq), a remarkably accurate transcriptome profiling method, is its high cost and poor scalability. Here, we report a highly scalable and cost-effective method for transcriptomics profiling called Bulk transcriptOme profiling of cell Lysate in a single poT (BOLT-seq), which is performed using unpurified bulk 3'-end mRNA in crude cell lysates. During BOLT-seq, RNA/DNA hybrids are directly subjected to tagmentation, and second-strand cDNA synthesis and RNA purification are omitted, allowing libraries to be constructed in 2 h of hands-on time. BOLT-seq was successfully used to cluster small molecule drugs based on their mechanisms of action and intended targets. BOLT-seq competes effectively with alternative library construction and transcriptome profiling methods.

Analytical and Clinical Validation of a Highly Sensitive NGS-Based ctDNA Assay with Real-World Concordance in NSCLC.

Purpose: There have been needs to improve the sensitivity of liquid biopsy. This report aims to report the analytical and clinical validation of a next generation sequencing (NGS)-based circulating tumor DNA (ctDNA) assay. Materials and Methods: Analytical validation was conducted in vitro by evaluating the limit of detection (LOD), precision, and specificity for various genomic aberrations. The real-world performance in non-small cell lung cancer (NSCLC) was assessed by comparing the results of AlphaLiquid®100 to the tissue-based results. Results: The LODs with 30 ng input DNA were 0.11%, 0.11%, 0.06%, 0.21%, and 2.13 copies for detecting SNVs, insertions, deletions, fusions, and copy number alterations (CNA), respectively. Quantitatively, SNV/INDELs, fusions, and CNAs showed a good correlation (R2=0.91, 0.40, and 0.65; y=0.95, 1.06, and 1.19) to the manufacturer's values, and per-base specificities for all types of variants were near 100%. In real-world NSCLC (n=122), key actionable mutations in NSCLC were detected in 60.7% (74/122) with the ctDNA assay. Comparative analysis against the NGS-based tissue results for all key mutations showed positive percent agreement (PPA) of 85.3%. For individual genes, the PPA was as high as 95.7% for EGFR mutations and 83.3% for ALK translocations. AlphaLiquid 100 detected drug-sensitive EGFR mutation at a variant allele frequency as low as 0.02% and also identified an EGFR mutation in a case where tissue sample missed. Blood samples collected post-targeted therapies revealed additional acquired mutations. Conclusion: The AlphaLiquid®100 ctDNA assay demonstrates robust analytical validity, offering clinically important information for NSCLC patients.

Total chemical synthesis of fully functional Photoactive Yellow Protein.

The Photoactive Yellow Protein (PYP) is a structural prototype for the PAS superfamily of proteins, which includes hundreds of receptor and regulatory proteins from all three kingdoms of life. PYP itself is a small globular protein that undergoes a photocycle involving a series of conformational changes in response to light excitation of its p-coumaric acid chromophore, making it an excellent model system to study the molecular basis of signaling in the PAS super family. To enable novel chemical approaches to elucidating the structural changes that accompany signaling in PYP, we have chemically synthesized the 125 amino acid residue protein molecule using a combination of Boc chemistry solid phase peptide synthesis and native chemical ligation. Synthetic PYP exhibits the wildtype photocycle, as determined in photobleaching studies. Planned future studies include incorporation of site-specific isotopic labels into specific secondary structural elements to determine which structural elements are involved in signaling state formation using difference FTIR spectroscopy.

EBS-seq: enrichment-based method for accurate analysis of 5-hydroxymethylcytosine at single-base resolution.

BACKGROUND: A growing body of research has emphasized 5-hydroxymethylcytosine (5hmC) as an important epigenetic mark. High-resolution methods to detect 5hmC require high sequencing depth and are therefore expensive. Many studies have used enrichment-based methods to detect 5hmC; however, conventional enrichment-based methods have limited resolution. To overcome these limitations, we developed EBS-seq, a cost-efficient method for 5hmC detection with single-base resolution that combines the advantages of high-resolution methods and enrichment-based methods. RESULTS: EBS-seq uses selective labeling of 5hmC, deamination of cytosine and 5-methylcytosine, pull-down of labeled 5hmC, and C-to-T conversion during DNA amplification. Using this method, we profiled 5hmC in HEK293T cells and two colorectal cancer samples. Compared with conventional enrichment-based 5hmC detection, EBS-seq improved 5hmC signals by localizing them at single-base resolution. Furthermore, EBS-seq was able to determine 5hmC levels in CpG-dense regions where distortion of signals can occur, such as CpG islands and CpG shores. Comparing EBS-seq and conventional high-resolution 5hmC detection by ACE-seq, we showed that EBS-seq is more effective at finding 5hmC sites. Using EBS-seq, we found strong associations between gene expression and gene-body 5hmC content in both HEK293T cells and colorectal cancer samples. CONCLUSIONS: EBS-seq is a reliable and cost-efficient method for 5hmC detection because it simultaneously enriches 5hmC-containing DNA fragments and localizes 5hmC signals at single-base resolution. This method is a promising choice for 5hmC detection in challenging clinical samples with low 5hmC levels, such as cancer tissues.

Efficacy of Olaparib in Treatment-Refractory, Metastatic Breast Cancer with Uncommon Somatic BRCA Mutations Detected in Circulating Tumor DNA.

Poly(ADP-ribose) polymerase inhibitors have been shown dramatic responses in patients with BRCAness. However, clinical studies have been limited to breast cancer patients with germline mutations. Here, we describe a patient with metastatic breast cancer who had a rare BRCA1 somatic mutation (BRCA1 c.4336G>T (p.E1446*)) detected by cell-free DNA analysis after failing standard therapies. This tier III variant of unknown significance was predicted to be a pathogenic variant in our assessment, leading us to consider off-label treatment with olaparib. The patient responded well to olaparib for several months, with a decrease in allele frequency of this BRCA1 somatic mutation in cell-free DNA. Olaparib resistance subsequently developed with an increase in the allele frequency and new BRCA1 reversion mutations. To our knowledge, this is the first report confirming BRCA1 c.4336G>T (p.E1446*) as a mutation sensitive to olaparib in breast cancer and describing the dynamic changes in the associated mutations using liquid biopsy.

Accurate Detection of Urothelial Bladder Cancer Using Targeted Deep Sequencing of Urine DNA.

Patients with hematuria are commonly given an invasive cystoscopy test to detect bladder cancer (BC). To avoid the risks associated with cystoscopy, several urine-based methods for BC detection have been developed, the most prominent of which is the deep sequencing of urine DNA. However, the current methods for urine-based BC detection have significant levels of false-positive signals. In this study, we report on uAL100, a method to precisely detect BC tumor DNA in the urine without tumor samples. Using urine samples from 43 patients with BC and 21 healthy donors, uAL100 detected BC with 83.7% sensitivity and 100% specificity. The mutations identified in the urine DNA by uAL100 for BC detection were highly associated with BC tumorigenesis and progression. We suggest that uAL100 has improved accuracy compared to other urine-based methods for early BC detection and can reduce unnecessary cystoscopy tests for patients with hematuria.

Globally shared TCR repertoires within the tumor-infiltrating lymphocytes of patients with metastatic gynecologic cancer.

Gynecologic cancer, including ovarian cancer and endometrial cancer, is characterized by morphological and molecular heterogeneity. Germline and somatic testing are available for patients to screen for pathogenic variants in genes such as BRCA1/2. Tissue expression levels of immunogenomic markers such as PD-L1 are also being used in clinical research. The basic therapeutic approach to gynecologic cancer combines surgery with chemotherapy. Immunotherapy, while not yet a mainstream treatment for gynecologic cancers, is advancing, with Dostarlimab recently receiving approval as a treatment for endometrial cancer. The goal remains to harness stimulated immune cells in the bloodstream to eradicate multiple metastases, a feat currently deemed challenging in a typical clinical setting. For the discovery of novel immunotherapy-based tumor targets, tumor-infiltrating lymphocytes (TILs) give a key insight on tumor-related immune activities by providing T cell receptor (TCR) sequences. Understanding the TCR repertoires of TILs in metastatic tissues and the circulation is important from an immunotherapy standpoint, as a subset of T cells in the blood have the potential to help kill tumor cells. To explore the relationship between distant tissue biopsy regions and blood circulation, we investigated the TCR beta chain (TCRß) in bulk tumor and matched blood samples from 39 patients with gynecologic cancer. We found that the TCR clones of TILs at different tumor sites were globally shared within patients and had high overlap with the TCR clones in peripheral blood.