CoolMPS: evaluation of antibody labeling based massively parallel non-coding RNA sequencing.

Results of massive parallel sequencing-by-synthesis vary depending on the sequencing approach. CoolMPS™ is a new sequencing chemistry that incorporates bases by labeled antibodies. To evaluate the performance, we sequenced 240 human non-coding RNA samples (dementia patients and controls) with and without CoolMPS. The Q30 value as indicator of the per base sequencing quality increased from 91.8 to 94%. The higher quality was reached across the whole read length. Likewise, the percentage of reads mapping to the human genome increased from 84.9 to 86.2%. For both technologies, we computed similar distributions between different RNA classes (miRNA, piRNA, tRNA, snoRNA and yRNA) and within the classes. While standard sequencing-by-synthesis allowed to recover more annotated miRNAs, CoolMPS yielded more novel miRNAs. The correlation between the two methods was 0.97. Evaluating the diagnostic performance, we observed lower minimal P-values for CoolMPS (adjusted P-value of 0.0006 versus 0.0004) and larger effect sizes (Cohen's d of 0.878 versus 0.9). Validating 19 miRNAs resulted in a correlation of 0.852 between CoolMPS and reverse transcriptase-quantitative polymerase chain reaction. Comparison to data generated with Illumina technology confirmed a known shift in the overall RNA composition. With CoolMPS we evaluated a novel sequencing-by-synthesis technology showing high performance for the analysis of non-coding RNAs.

CoolMPS for robust sequencing of single-nuclear RNAs captured by droplet-based method.

Massively-parallel single-cell and single-nucleus RNA sequencing (scRNA-seq, snRNA-seq) requires extensive sequencing to achieve proper per-cell coverage, making sequencing resources and availability of sequencers critical factors for conducting deep transcriptional profiling. CoolMPS is a novel sequencing-by-synthesis approach that relies on nucleotide labeling by re-usable antibodies, but whether it is applicable to snRNA-seq has not been tested. Here, we use a low-cost and off-the-shelf protocol to chemically convert libraries generated with the widely-used Chromium 10X technology to be sequenceable with CoolMPS technology. To assess the quality and performance of converted libraries sequenced with CoolMPS, we generated a snRNA-seq dataset from the hippocampus of young and old mice. Native libraries were sequenced on an Illumina Novaseq and libraries that were converted to be compatible with CoolMPS were sequenced on a DNBSEQ-400RS. CoolMPS-derived data faithfully replicated key characteristics of the native library dataset, including correct estimation of ambient RNA-contamination, detection of captured cells, cell clustering results, spatial marker gene expression, inter- and intra-replicate differences and gene expression changes during aging. In conclusion, our results show that CoolMPS provides a viable alternative to standard sequencing of RNA from droplet-based libraries.

Efficient and unique cobarcoding of second-generation sequencing reads from long DNA molecules enabling cost-effective and accurate sequencing, haplotyping, and de novo assembly.

Here, we describe single-tube long fragment read (stLFR), a technology that enables sequencing of data from long DNA molecules using economical second-generation sequencing technology. It is based on adding the same barcode sequence to subfragments of the original long DNA molecule (DNA cobarcoding). To achieve this efficiently, stLFR uses the surface of microbeads to create millions of miniaturized barcoding reactions in a single tube. Using a combinatorial process, up to 3.6 billion unique barcode sequences were generated on beads, enabling practically nonredundant cobarcoding with 50 million barcodes per sample. Using stLFR, we demonstrate efficient unique cobarcoding of more than 8 million 20- to 300-kb genomic DNA fragments. Analysis of the human genome NA12878 with stLFR demonstrated high-quality variant calling and phase block lengths up to N50 34 Mb. We also demonstrate detection of complex structural variants and complete diploid de novo assembly of NA12878. These analyses were all performed using single stLFR libraries, and their construction did not significantly add to the time or cost of whole-genome sequencing (WGS) library preparation. stLFR represents an easily automatable solution that enables high-quality sequencing, phasing, SV detection, scaffolding, cost-effective diploid de novo genome assembly, and other long DNA sequencing applications.

Reliable multiplex sequencing with rare index mis-assignment on DNB-based NGS platform.

BACKGROUND: Massively-parallel-sequencing, coupled with sample multiplexing, has made genetic tests broadly affordable. However, intractable index mis-assignments (commonly exceeds 1%) were repeatedly reported on some widely used sequencing platforms. RESULTS: Here, we investigated this quality issue on BGI sequencers using three library preparation methods: whole genome sequencing (WGS) with PCR, PCR-free WGS, and two-step targeted PCR. BGI's sequencers utilize a unique DNA nanoball (DNB) technology which uses rolling circle replication for DNA-nanoball preparation; this linear amplification is PCR free and can avoid error accumulation. We demonstrated that single index mis-assignment from free indexed oligos occurs at a rate of one in 36 million reads, suggesting virtually no index hopping during DNB creation and arraying. Furthermore, the DNB-based NGS libraries have achieved an unprecedentedly low sample-to-sample mis-assignment rate of 0.0001 to 0.0004% under recommended procedures. CONCLUSIONS: Single indexing with DNB technology provides a simple but effective method for sensitive genetic assays with large sample numbers.

Clinical and genetic analysis of a rare syndrome associated with neoteny.

PurposeWe describe a novel syndrome in seven female patients with extreme developmental delay and neoteny.MethodsAll patients in this study were female, aged 4 to 23 years, were well below the fifth percentile in height and weight, had failed to develop sexually, and lacked the use of language. Karyotype and array chromosome genomic hybridization analysis failed to identify large-scale structural variations. To further understand the underlying cause of disease in these patients, whole-genome sequencing was performed.ResultsIn five patients, coding de novo mutations (DNMs) were found in five different genes. These genes fell into similar functional categories of transcription regulation and chromatin modification. Comparison to a control population suggested that individuals with neotenic complex syndrome (NCS)-a name that we propose herein-could have an excess of rare inherited variants in genes associated with developmental delay and autism, although the difference was not significant.ConclusionWe describe an extreme form of developmental delay, with the defining characteristic of neoteny. In most patients we identified coding DNMs in a set of genes intolerant of haploinsufficiency; however, it is not clear whether these contributed to NCS. Rare inherited variants may also be associated with NCS, but more samples need to be analyzed to achieve statistical significance.

Isolation and whole genome sequencing of fetal cells from maternal blood towards the ultimate non-invasive prenatal testing.

OBJECTIVE: The purpose of this study were to develop a methodology of isolating fetal cells from maternal blood and use deep sequence demonstrating the promise for complete and accurate genetic screening compared to other non-invasive prenatal testing. METHODS: Here in this study, we developed a double negative selection (DNS) procedure to unbiasedly enrich fetal cells. After validated by short tandem repeat (STR), the isolated circulating fetal cells (CFCs) were subjected to deep whole genome sequencing analysis. RESULTS: Our DNS protocol significantly increasing the purity of the mimic fetal cells from 1 in 1 million nucleated cells in whole blood to 1:8 to 1:30 (12.5%-3.33%) after 2 steps of enrichment. Isolated single fetal cell obtained a coverage rate (86.8%) and allelic dropout rate (24.90%) comparative to the reported results of human cell line. Several disease-associated variants were identified in the whole genome sequencing data of isolated CFCs and further confirmed in the sequencing data of unamplified gDNA. CONCLUSION: In conclusion, the robustness of DNS and STR to collect CFCs from peripheral maternal blood for the first time coupled with deep sequencing technique demonstrates the possibility of comprehensive non-invasive prenatal testing of genetic disorders using isolated CFCs.

Quantitative Whole Genome Sequencing of Circulating Tumor Cells Enables Personalized Combination Therapy of Metastatic Cancer.

Much effort has been dedicated to developing circulating tumor cells (CTC) as a noninvasive cancer biopsy, but with limited success as yet. In this study, we combine a method for isolation of highly pure CTCs using immunomagnetic enrichment/fluorescence-activated cell sorting with advanced whole genome sequencing (WGS), based on long fragment read technology, to illustrate the utility of an accurate, comprehensive, phased, and quantitative genomic analysis platform for CTCs. Whole genomes of 34 CTCs from a patient with metastatic breast cancer were analyzed as 3,072 barcoded subgenomic compartments of long DNA. WGS resulted in a read coverage of 23× per cell and an ensemble call rate of >95%. These barcoded reads enabled accurate detection of somatic mutations present in as few as 12% of CTCs. We found in CTCs a total of 2,766 somatic single-nucleotide variants and 543 indels and multi-base substitutions, 23 of which altered amino acid sequences. Another 16,961 somatic single nucleotide variant and 8,408 indels and multi-base substitutions, 77 of which were nonsynonymous, were detected with varying degrees of prevalence across the 34 CTCs. On the basis of our whole genome data of mutations found in all CTCs, we identified driver mutations and the tissue of origin of these cells, suggesting personalized combination therapies beyond the scope of most gene panels. Taken together, our results show how advanced WGS of CTCs can lead to high-resolution analyses of cancers that can reliably guide personalized therapy. Cancer Res; 77(16); 4530-41. ©2017 AACR.

cPAS-based sequencing on the BGISEQ-500 to explore small non-coding RNAs.

BACKGROUND: We present the first sequencing data using the combinatorial probe-anchor synthesis (cPAS)-based BGISEQ-500 sequencer. Applying cPAS, we investigated the repertoire of human small non-coding RNAs and compared it to other techniques. RESULTS: Starting with repeated measurements of different specimens including solid tissues (brain and heart) and blood, we generated a median of 30.1 million reads per sample. 24.1 million mapped to the human genome and 23.3 million to the miRBase. Among six technical replicates of brain samples, we observed a median correlation of 0.98. Comparing BGISEQ-500 to HiSeq, we calculated a correlation of 0.75. The comparability to microarrays was similar for both BGISEQ-500 and HiSeq with the first one showing a correlation of 0.58 and the latter one correlation of 0.6. As for a potential bias in the detected expression distribution in blood cells, 98.6% of HiSeq reads versus 93.1% of BGISEQ-500 reads match to the 10 miRNAs with highest read count. After using miRDeep2 and employing stringent selection criteria for predicting new miRNAs, we detected 74 high-likely candidates in the cPAS sequencing reads prevalent in solid tissues and 36 candidates prevalent in blood. CONCLUSIONS: While there is apparently no ideal platform for all challenges of miRNome analyses, cPAS shows high technical reproducibility and supplements the hitherto available platforms.

De novo assembly of a haplotype-resolved human genome.

The human genome is diploid, and knowledge of the variants on each chromosome is important for the interpretation of genomic information. Here we report the assembly of a haplotype-resolved diploid genome without using a reference genome. Our pipeline relies on fosmid pooling together with whole-genome shotgun strategies, based solely on next-generation sequencing and hierarchical assembly methods. We applied our sequencing method to the genome of an Asian individual and generated a 5.15-Gb assembled genome with a haplotype N50 of 484 kb. Our analysis identified previously undetected indels and 7.49 Mb of novel coding sequences that could not be aligned to the human reference genome, which include at least six predicted genes. This haplotype-resolved genome represents the most complete de novo human genome assembly to date. Application of our approach to identify individual haplotype differences should aid in translating genotypes to phenotypes for the development of personalized medicine.

Combinatorial sequencing-by-hybridization: analysis of the NF1 gene.

Neurofibromatosis type 1 (NF1), one of the most common autosomal dominant disorders, is caused by mutations in the NF1 gene. A variety of methods are currently used in clinical settings to define disease-causing mutations. We describe microarray-based combinatorial sequencing-by-hybridization (cSBH), which overcomes some disadvantages associated with other techniques. Sequence readout of 2 kb was achieved on a single slide, with detection of base substitutions, insertions and small deletions. In addition, cSBH analysis of the entire NF1 gene demonstrates reproducibility, efficiency and reduced time; therefore, representing an alternative to extensive DNA sequence characterization.

Multiplexed SNP genotyping using nanobarcode particle technology.

Single-nucleotide polymorphisms (SNP) are the most common form of sequence variation in the human genome. Large-scale studies demand high-throughput SNP genotyping platforms. Here we demonstrate the potential of encoded nanowires for use in a particles-based universal array for high-throughput SNP genotyping. The particles are encoded sub-micron metallic nanorods manufactured by electroplating inert metals such as gold and silver into templates and releasing the resulting striped nanoparticles. The power of this technology is that the particles are intrinsically encoded by virtue of the different reflectivity of adjacent metal stripes, enabling the generation of many thousands of unique encoded substrates. Using SNP found within the cytochrome P450 gene family, and a universal short oligonucleotide ligation strategy, we have demonstrated the simultaneous genotyping of 15 SNP; a format requiring discrimination of 30 encoded nanowires (one per allele). To demonstrate applicability to real-world applications, 160 genotypes were determined from multiplex PCR products from 20 genomic DNA samples.

Identification of APC gene mutations in colorectal cancer using universal microarray-based combinatorial sequencing-by-hybridization.

Familial adenomatous polyposis (FAP) is an autosomal dominant inherited form of colorectal cancer, caused mostly by mutations in the APC gene. Due to the wide variety of mutations found and the large size of the APC gene, several methods of mutation detection are used, which can be time consuming and costly. Here we demonstrate a new method of mutation detection in the APC gene using an array-based approach termed combinatorial sequencing-by-hybridization (cSBH). In cSBH, a universal probe set is attached to a support and a second one is in solution. Two-probe ligation occurs when a DNA strand from the target PCR product consecutively anneals to both unlabeled array-bound and solution-phase dye-labeled probe, creating all target-complementary long-labeled probes attached to the surface. A standard array reader scores fluorescent signals at each array position. Cell lines and patient DNA with known APC gene mutations were analyzed using a cSBH-based HyChip trade mark product. Results show that this universal hexamer (6-mer) chip can successfully detect a range of mutations. Results are very robust for a continuous readout of 3.6 kb from a PCR target, with 99.97% accuracy on a single HyChip trade mark slide. cSBH is a fast, cost-efficient method for first stage mutation screening in the APC or any other gene.

Selective DNA amplification from complex genomes using universal double-sided adapters.

There is a rapidly developing need for new technologies to amplify millions of different targets from genomic DNA for high throughput genotyping and population gene-sequencing from diverse species. Here we describe a novel approach for the specific selection and amplification of genomic DNA fragments of interest that eliminates the need for costly and time consuming synthesis and testing of potentially millions of amplicon-specific primers. This technique relies upon Type IIs restriction enzyme digestion of genomic DNA and ligation of the fragments to double-sided adapters to form closed-circular DNA molecules. The novel use of double-sided adapters, assembled through the combinatorial use of two small universal sets of oligonucleotide building blocks, provides greater selection capacity by utilizing both sides of the adapter in a sequence-specific ligation event. As demonstrated, formation of circular structures results in protection of the desired molecules from nuclease treatment and enables a level of selectivity high enough to isolate single, or multiple, pre-defined fragments from the human genome when digested at over five million sites. Priming sites incorporated into the adapter allows the utilization of a common pair of primers for the amplification of any adapter-captured DNA fragment of interest.

Sequencing by hybridization (SBH): advantages, achievements, and opportunities.

Efficient DNA sequencing of the genomes of individual species and organisms is a critical task for the advancement of biological sciences, medicine and agriculture. Advances in modern sequencing methods are needed to meet the challenge of sequencing such megabase to gigabase quantities of DNA. Two possible strategies for DNA sequencing exist: direct methods, in which each base position in the DNA chain is determined individually (e.g., gel sequencing or pyrosequencing), and indirect methods, in which the DNA sequence is assembled based on experimental determination of oligonucleotide content of the DNA chain. One promising indirect method is sequencing by hybridization (SBH), in which sets of oligonucleotides are hybridized under conditions that allow detection of complementary sequences in the target nucleic acid. The unprecedented sequence search parallelism of the SBH method has allowed development of high-throughput, low-cost, miniaturized sequencing processes on arrays of DNA samples or probes. Newly developed SBH methods use DNA ligation to combine relatively small sets of short probes to score potentially tens of millions of longer oligonucleotide sequences in a target DNA. Such combinatorial approaches allow analysis of DNA samples of up to several kilobases (several times longer than allowed by current direct methods) for a variety of DNA sequence analysis applications, including de novo sequencing, resequencing, mutation/SNP discovery and genotyping, and expression monitoring. Future advances in biochemistry and implementation of detection methods that allow single-molecule sensitivity may provide the necessary miniaturization, specificity, and multiplexing efficiency to allow routine whole genome analysis in a single solution-based hybridization experiment.