Snezana Drmanac

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.

cDNA screening by array hybridization.

Publisher Summary Sequencing by hybridization (SBH) applied to dense arrays of clones achieves an unprecedented throughput of analyzed samples that can be performed inexpensively. The adaptation of this technology to an industrial setting has already produced 1 million partial cDNA sequences per month and a database of more than 4 million individual cDNA clones in the facility. This chapter discusses the present protocols and the potential application of this technology to the fields of gene discovery and gene function. SBH technology applied on an array of clones provides a progressive approach from expression screening to complete sequencing of complex libraries. The same approach is applicable for mutation detection and complete sequencing of many genes amplified by polymerase chain reaction (PCR) from thousands of individual samples.

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.

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

BackgroundWe 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.ResultsStarting 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.ConclusionsWhile there is apparently no ideal platform for all challenges of miRNome analyses, cPAS shows high technical reproducibility and supplements the hitherto available platforms.

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

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.

Requirements in screening cDNA libraries for new genes and solutions offered by SBH technology

Under different assumptions about the total number of genes, the number of housekeeping and tissue-specific genes, and the difference in the number of mRNAs per cell for functional and nonfunctional genes, significantly different results can be expected from screening random cDNA clones. We have developed gene expression models as a guide for interpretation of experimental results. For statistical, biological, and technical reasons, the search for 100,000 plus genes and discrimination between nonfunctional, housekeeping, and tissue-specific genes requires the analysis of up to 10 million clones from 20 to 50 tissues. Oligonucleotide hybridization of dense clone blots is an inexpensive and fast way to screen such large clone sets. Our preliminary results on control clones and thousands of cDNA clones from an infant brain library demonstrate the feasibility of the method.

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.

Single tube bead-based DNA co-barcoding for cost effective and accurate sequencing, haplotyping, and assembly

Obtaining accurate sequences from long DNA molecules is very important for genome assembly and other applications. Here we describe single tube long fragment read (stLFR), a technology that enables this a low cost. It is based on adding the same barcode sequence to sub-fragments of the original long DNA molecule (DNA co-barcoding). 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 non-redundant co-barcoding with 50 million barcodes per sample. Using stLFR, we demonstrate efficient unique co-barcoding of over 8 million 20-300 kb genomic DNA fragments. Analysis of the genome of the human genome NA12878 with stLFR demonstrated high quality variant calling and phasing into contigs 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.Single tube long fragment read (stLFR) technology enables efficient WGS, haplotyping, and contig scaffolding. It is based on adding the same barcode sequence to sub-fragments of the original DNA molecule (DNA co-barcoding). To achieve this, stLFR uses the surface of microbeads to create millions of miniaturized compartments in a single tube. Using a combinatorial process over 1.8 billion unique barcode sequences were generated on beads, enabling practically non-redundant co-barcoding in reactions with 50 million barcodes. Using stLFR we demonstrate efficient unique co-barcoding of over 8 million 20300 kb genomic DNA fragments with near perfect variant calling and phasing of the genome of NA12878 into contigs up to N50 23.4 Mb. stLFR represents a low-cost single library solution that can enable long sequence data.

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