Sequencing by Binding rivals error-corrected Sequencing by Synthesis technology for accurate detection and quantification of minor (<0.1%) subpopulation variants.

Detecting very minor (< 1%) subpopulations using next-generation sequencing is a critical need for multiple applications including detection of drug resistant pathogens and somatic variant detection in oncology. To enable these applications, wet lab enhancements and bioinformatic error correction methods have been developed for 'sequencing by synthesis' technology to reduce its inherent sequencing error rate. A recently available sequencing approach termed 'sequencing by binding' claims to have higher base calling accuracy data "out of the box." This paper evaluates the utility of using 'sequencing by binding' for the detection of ultra-rare subpopulations down to 0.001%.

Developmental validation of the MiSeq FGx Forensic Genomics System for Targeted Next Generation Sequencing in Forensic DNA Casework and Database Laboratories.

Human DNA profiling using PCR at polymorphic short tandem repeat (STR) loci followed by capillary electrophoresis (CE) size separation and length-based allele typing has been the standard in the forensic community for over 20 years. Over the last decade, Next-Generation Sequencing (NGS) matured rapidly, bringing modern advantages to forensic DNA analysis. The MiSeq FGx™ Forensic Genomics System, comprised of the ForenSeq™ DNA Signature Prep Kit, MiSeq FGx™ Reagent Kit, MiSeq FGx™ instrument and ForenSeq™ Universal Analysis Software, uses PCR to simultaneously amplify up to 231 forensic loci in a single multiplex reaction. Targeted loci include Amelogenin, 27 common, forensic autosomal STRs, 24 Y-STRs, 7 X-STRs and three classes of single nucleotide polymorphisms (SNPs). The ForenSeq™ kit includes two primer sets: Amelogenin, 58 STRs and 94 identity informative SNPs (iiSNPs) are amplified using DNA Primer Set A (DPMA; 153 loci); if a laboratory chooses to generate investigative leads using DNA Primer Set B, amplification is targeted to the 153 loci in DPMA plus 22 phenotypic informative (piSNPs) and 56 biogeographical ancestry SNPs (aiSNPs). High-resolution genotypes, including detection of intra-STR sequence variants, are semi-automatically generated with the ForenSeq™ software. This system was subjected to developmental validation studies according to the 2012 Revised SWGDAM Validation Guidelines. A two-step PCR first amplifies the target forensic STR and SNP loci (PCR1); unique, sample-specific indexed adapters or "barcodes" are attached in PCR2. Approximately 1736 ForenSeq™ reactions were analyzed. Studies include DNA substrate testing (cotton swabs, FTA cards, filter paper), species studies from a range of nonhuman organisms, DNA input sensitivity studies from 1ng down to 7.8pg, two-person human DNA mixture testing with three genotype combinations, stability analysis of partially degraded DNA, and effects of five commonly encountered PCR inhibitors. Calculations from ForenSeq™ STR and SNP repeatability and reproducibility studies (1ng template) indicate 100.0% accuracy of the MiSeq FGx™ System in allele calling relative to CE for STRs (1260 samples), and >99.1% accuracy relative to bead array typing for SNPs (1260 samples for iiSNPs, 310 samples for aiSNPs and piSNPs), with >99.0% and >97.8% precision, respectively. Call rates of >99.0% were observed for all STRs and SNPs amplified with both ForenSeq™ primer mixes. Limitations of the MiSeq FGx™ System are discussed. Results described here demonstrate that the MiSeq FGx™ System meets forensic DNA quality assurance guidelines with robust, reliable, and reproducible performance on samples of various quantities and qualities.

A somatic reference standard for cancer genome sequencing.

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

Open-access synthetic spike-in mRNA-seq data for cancer gene fusions.

BACKGROUND: Oncogenic fusion genes underlie the mechanism of several common cancers. Next-generation sequencing based RNA-seq analyses have revealed an increasing number of recurrent fusions in a variety of cancers. However, absence of a publicly available gene-fusion focused RNA-seq data impedes comparative assessment and collaborative development of novel gene fusions detection algorithms. We have generated nine synthetic poly-adenylated RNA transcripts that correspond to previously reported oncogenic gene fusions. These synthetic RNAs were spiked at known molarity over a wide range into total RNA prior to construction of next-generation sequencing mRNA libraries to generate RNA-seq data. RESULTS: Leveraging a priori knowledge about replicates and molarity of each synthetic fusion transcript, we demonstrate utility of this dataset to compare multiple gene fusion algorithms' detection ability. In general, more fusions are detected at higher molarity, indicating that our constructs performed as expected. However, systematic detection differences are observed based on molarity or algorithm-specific characteristics. Fusion-sequence specific detection differences indicate that for applications where specific sequences are being investigated, additional constructs may be added to provide quantitative data that is specific for the sequence of interest. CONCLUSIONS: To our knowledge, this is the first publicly available synthetic RNA-seq data that specifically leverages known cancer gene-fusions. The proposed method of designing multiple gene-fusion constructs over a wide range of molarity allows granular performance analyses of multiple fusion-detection algorithms. The community can leverage and augment this publicly available data to further collaborative development of analytical tools and performance assessment frameworks for gene fusions from next-generation sequencing data.

A highly scalable peptide-based assay system for proteomics.

We report a scalable and cost-effective technology for generating and screening high-complexity customizable peptide sets. The peptides are made as peptide-cDNA fusions by in vitro transcription/translation from pools of DNA templates generated by microarray-based synthesis. This approach enables large custom sets of peptides to be designed in silico, manufactured cost-effectively in parallel, and assayed efficiently in a multiplexed fashion. The utility of our peptide-cDNA fusion pools was demonstrated in two activity-based assays designed to discover protease and kinase substrates. In the protease assay, cleaved peptide substrates were separated from uncleaved and identified by digital sequencing of their cognate cDNAs. We screened the 3,011 amino acid HCV proteome for susceptibility to cleavage by the HCV NS3/4A protease and identified all 3 known trans cleavage sites with high specificity. In the kinase assay, peptide substrates phosphorylated by tyrosine kinases were captured and identified by sequencing of their cDNAs. We screened a pool of 3,243 peptides against Abl kinase and showed that phosphorylation events detected were specific and consistent with the known substrate preferences of Abl kinase. Our approach is scalable and adaptable to other protein-based assays.

Detection of low prevalence somatic mutations in solid tumors with ultra-deep targeted sequencing.

Ultra-deep targeted sequencing (UDT-Seq) can identify subclonal somatic mutations in tumor samples. Early assays' limited breadth and depth restrict their clinical utility. Here, we target 71 kb of mutational hotspots in 42 cancer genes. We present novel methods enhancing both laboratory workflow and mutation detection. We evaluate UDT-Seq true sensitivity and specificity (> 94% and > 99%, respectively) for low prevalence mutations in a mixing experiment and demonstrate its utility using six tumor samples. With an improved performance when run on the Illumina Miseq, the UDT-Seq assay is well suited for clinical applications to guide therapy and study clonal selection in heterogeneous samples.

HIV-1 Rev protein assembles on viral RNA one molecule at a time.

Oligomerization of the HIV-1 protein Rev on the Rev Response Element (RRE) regulates nuclear export of genomic viral RNA and partially spliced viral mRNAs encoding for structural proteins. Single-molecule fluorescence spectroscopy has been used to dissect the multistep assembly pathway of this essential ribonucleoprotein, revealing dynamic intermediates and the mechanism of assembly. Assembly is initiated by binding of Rev to a high-affinity site in stem-loop IIB of the RRE and proceeds rapidly by addition of single Rev monomers, facilitated by cooperative Rev-Rev interactions on the RRE. Dwell-time analysis of fluorescence trajectories recorded during individual Rev-RRE assembly reactions has revealed the microscopic rate constants for several of the Rev monomer binding and dissociation steps. The high-affinity binding of multiple Rev monomers to the RRE is achieved on a much faster timescale than reported in previous bulk kinetic studies of Rev-RRE association, indicating that oligomerization is an early step in complex assembly.

Two-photon absorption in three-dimensional chromophores based on [2.2]-paracyclophane.

A series of alpha,omega-bis donor substituted oligophenylenevinylene dimers held together by the [2.2]paracyclophane core were synthesized to probe how the number of repeat units and through-space delocalization influence two-photon absorption cross sections. Specifically, the paracyclophane molecules are tetra(4,7,12,15)-(4'-dihexylaminostyryl)[2.2]paracyclophane (3R(D)), tetra(4,7,12,15)-(4' '-(4'-dihexylaminostyryl)styryl)[2.2]paracyclophane (5R(D)), and tetra(4,7,12,15)-(4' "-(4' '-(4'-dihexylaminostyryl)styryl)styryl)[2.2]paracyclophane (7R(D)). The compounds bis(1,4)-(4'-dihexylaminostyryl)benzene (3R) and bis(1,4)-(4' '-(4'-dihexylaminostyryl)styryl)benzene (5R) were also synthesized to reveal the properties of the "monomeric" counterparts. The two-photon absorption cross sections were determined by the two-photon induced fluorescence method using both femtosecond and nanosecond pulsed lasers as excitation sources. While there is a red shift in the linear absorption spectra when going from the "monomer" chromophore to the paracyclophane "dimer" (i.e., 3R --> 3R(D), 5R --> 5R(D)), there is no shift in the two-photon absorption maxima. A theoretical treatment of these trends and the dependence of transition dipole moments on molecular structure rely on calculations that interfaced time-dependent density functional theory (TDDFT) techniques with the collective electronic oscillator (CEO) program. These theoretical and experimental results indicate that intermolecular interactions can strongly affect B(u) states but weakly perturb A(g) states, due to the small dipole-dipole coupling between A(g) states on the chromophores in the dimer.

Metal-ion sensing fluorophores with large two-photon absorption cross sections: aza-crown ether substituted donor-acceptor-donor distyrylbenzenes.

Chromophores based on a donor-acceptor-donor structure possessing a large two-photon absorption cross section and one or two mono-aza-15-crown-5 ether moieties, which can bind metal cations, have been synthesized. The influence of Mg(2+) binding on their one- and two-photon spectroscopic properties has been investigated. Upon binding, the two-photon action cross sections at 810 nm decrease by a factor of up to 50 at high Mg(2+) concentrations and this results in a large contrast in the two-photon excited fluorescence signal between the bound and unbound forms, for excitation in the range of 730 to 860 nm. Experimental and computational results indicate that there is a significant reduction of the electron donating strength of the aza-crown nitrogen atom(s) upon metal ion binding and that this leads to a blue shift in the position as well as a reduction in the strength of the lowest-energy two-photon absorption band. The molecules reported here can serve as models for the design of improved two-photon excitable metal-ion sensing fluorophores.

Real time differentiation of G-protein coupled receptor (GPCR) agonist and antagonist by two photon fluorescence laser microscopy.

Receptor-based signaling mechanisms are the primary source of cellular regulation. The superfamily of G protein-coupled receptors (GPCR) is the largest and most ubiquitous of the receptor-mediated processes. Desensitization of G-protein-coupled receptors is a fundamental mechanism regulating the cellular response to agonists. We have recently studied the agonist and antagonist of the human melanocortin receptors (hMC1, hMC3, hMC4, and hMC5 receptors), the human delta opioid receptor, and the human gluacagon receptor with the help of synthetic fluorescent labeled ligands and fluorescent protein-labeled beta-arrestin-receptors that shed new insight on cellular signaling and rapid screening of drugs in real time. It was demonstrated that stimulation of these receptors by the cognate agonist triggers the rapid internalization of ligand-receptor complexes, while the interaction of the receptor with antagonists does not follow this pathway. Furthermore, receptor internalization is dependent upon beta-arrestin, which has been shown to be responsible for the rapid desensitization of cAMP-signaling processes.