Single-cell Genomic DNA Research Articles

Raman-Activated Cell Ejection for Validating the Reliability of the Raman Fingerprint Database of Foodborne Pathogens

Raman spectroscopy for rapid identification of foodborne pathogens based on phenotype has attracted increasing attention, and the reliability of the Raman fingerprint database through genotypic determination is crucial. In the research, the classification model of four foodborne pathogens was established based on t-distributed stochastic neighbor embedding (t-SNE) and support vector machine (SVM); the recognition accuracy was 97.04%. The target bacteria named by the model were ejected through Raman-activated cell ejection (RACE), and then single-cell genomic DNA was amplified for species analysis. The accuracy of correct matches between the predicted phenotype and the actual genotype of the target cells was at least 83.3%. Furthermore, all anticipant sequencing results brought into correspondence with the species were predicted through the model. In sum, the Raman fingerprint database based on Raman spectroscopy combined with machine learning was reliable and promising in the field of rapid detection of foodborne pathogens.

Abstract 2702: Analysis of low-pass sequencing data reveals extensive loss-of-heterozygosity in circulating multiple myeloma cells

Abstract Background: Clonal evolution and heterogeneity are among the factors that make treatment of Multiple Myeloma (MM) challenging, and disease monitoring requires invasive bone marrow biopsies. Tumor heterogeneity has previously been demonstrated at copy-number level, showing patterns of recurring Copy-Number Alterations (CNA). Conversely, Loss-of-Heterozygosity (LoH) profiling has been investigated less extensively. Presence of genomic patterns of LoH has been associated with Homologous Recombination Deficiency (HRD) in MM, suggesting a role for therapy with PARP inhibitors. Here we report the genome-wide LoH profiling in single-CMMCs isolated from enriched peripheral blood of four different MM patients using a non-invasive approach that combines CellSearch® and DEPArray™ technology. Methods: CMMCs were obtained from peripheral blood of four MM patients using CellSearch for enrichment and DEPArray NxT for isolation. CMMCs enrichment was obtained using a custom kit with anti-CD138 or anti-CD138/CD38 antibody-conjugated ferrofluids for positive enrichment and CD38-PE, CD19/CD45-APC immunofluorescent staining for detection. After single cell isolation, Ampli1™ Whole Genome Amplification (WGA) kit was used to amplify single-cell genomic DNA. Whole genome sequencing libraries were prepared from WGA products using Ampli1™ LowPass kit and low-pass sequencing was performed on HiSeq 2500 Illumina® platform. Genome-wide single-cell Low-Pass Copy Number Alteration (LPCNA) analysis was performed using the cloud-based bioinformatic suite MSBiosuite. Single-cell genome-wide LoH profiles were obtained by a custom algorithm based on the ratio between mono- and bi-allelic loci in non-overlapping segments of uniform copy-number. Results: Out of 315 cells analyzed from the four MM patients, 244 (77%) showed LoH events. Extensive LoH regions were detected in all patients, with 231 cells showing copy-neutral LoH events and 121 showing LoH events in copy-gain regions. Globally, 460 copy-neutral LoHs, and 292 LoHs in copy-gain regions were detected. Two highly conserved copy-neutral events (9p; 16q) were detected in one of the patients analyzed, suggesting they are early events. Interestingly, the first of the 2 events includes CDKN2A, a tumor suppressor which has been implicated in predisposition to MM. Notably, in a second patient, chromosome 11 is associated both to a copy-neutral LoH and to an LoH event in a copy-gain segment including the CCND1 gene, whose amplification in MM has been associated with multidrug resistance expression. Conclusions: Here we present a non-invasive workflow for the molecular characterization of single CMMCs isolated from peripheral blood. Specifically, we reported the presence of chromosomal or sub-chromosomal LoH events in all MM patients, suggesting a pervasive phenomenon in MM which deserves to be further explored as potential therapy target. Citation Format: Claudio Forcato, Alberto Ferrarini, Genny Buson, Paola Tononi, Marianna Garonzi, Valentina del Monaco, Andrea Raspadori, Steven Gross, Francesca Fontana, Gianni Medoro, Mark Connelly, Nicolò Manaresi. Analysis of low-pass sequencing data reveals extensive loss-of-heterozygosity in circulating multiple myeloma cells [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 2702.

Using SPAdes De Novo Assembler.

SPAdes-St. Petersburg genome Assembler-was originally developed for de novo assembly of genome sequencing data produced for cultivated microbial isolates and for single-cell genomic DNA sequencing. With time, the functionality of SPAdes was extended to enable assembly of IonTorrent data, as well as hybrid assembly from short and long reads (PacBio and Oxford Nanopore). In this article we present protocols for five different assembly pipelines that comprise the SPAdes package and that are used for assembly of metagenomes and transcriptomes as well as assembly of putative plasmids and biosynthetic gene clusters from whole-genome sequencing and metagenomic datasets. In addition, we present guidelines for understanding results with use cases for each pipeline, and several additional support protocols that help in using SPAdes properly. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Assembling isolate bacterial datasets Basic Protocol 2: Assembling metagenomic datasets Basic Protocol 3: Assembling sets of putative plasmids Basic Protocol 4: Assembling transcriptomes Basic Protocol 5: Assembling putative biosynthetic gene clusters Support Protocol 1: Installing SPAdes Support Protocol 2: Providing input via command line Support Protocol 3: Providing input data via YAML format Support Protocol 4: Restarting previous run Support Protocol 5: Determining strand-specificity of RNA-seq data.

Abstract 2911: Genome-wide copy number profiling of single circulating multiple myeloma cells (CMMCs) reveals intra-patient convergent copy-number alterations (CNAs)

Abstract Introduction: Multiple Myeloma (MM) evolution and heterogeneity is interesting for its potential translational relevance. Recent studies using bulk sequencing of Smoldering MM cells obtained from Bone Marrow (BM) report recurring CNA patterns. By analyzing single-CMMCs isolated from enriched peripheral blood, we show here, for the first time in MM, evidence of frequent convergent lesions, i.e. alterations developed independently across different evolutionary branches, including CNAs often found in MM and previously reported as common truncal alterations. Methods: Peripheral blood samples (4.0 ml) were obtained from n=3 patients with MM. CMMCs were enriched with CellSearch® AutoPrep® using a custom kit with anti-CD138 or anti-CD138/CD38 antibody-conjugated ferrofluids for positive enrichment and CD38-PE, CD19/CD45-APC immunofluorescent staining for detection. Cell enumeration was based on the co-localization of nuclear DAPI staining and CD38-PE on CellSearch CTAII®. Single CMMCs (CD38+/CD19- and CD45-/DAPI+) and White Blood Cells (WBCs: CD38-/CD19+ or CD45+/DAPI+) were then isolated with DEPArray NxT system. Single-cell genomic DNA was amplified using Ampli1™ Whole Genome Amplification (WGA) kit, Illumina®-compatible libraries were obtained using Ampli1™ LowPass kit and the process was automated on a Hamilton STARLet Liquid handler. Multiplexed, low-pass whole-genome sequencing was performed on HiSeq 2500 Illumina® platform. Genome-wide single-cell Low-Pass Copy Number Alteration (LPCNA) analysis was performed using the cloud-based bioinformatic suite MSBio Suite (Menarini Silicon Biosystems). Results: 186/215 (86%) single CMMC in total passed QC criteria and were included in the analysis. Single WBCs were also included as normal controls. Cumulatively, CNA profiles of single CMMCs showed patterns typical of MM, including 1q gain, 13 monosomy, sub-chromosomal gain or trisomy 3, 5, 7, 9, 11, 15, 21. Of these, intra-patient single-cell profiling surprisingly revealed -in all three patients- convergent lesions, i.e. alterations developed independently across different evolutionary branches, along with conserved (common truncal), and divergent alterations (found only in specific sub-clusters). In addition, we found evidence that 1q gain, 13q deletion and 6p gain were actually subclonal, in contrast with recent publications reporting them as truncal early-onset lesions. Conclusion: Single CMMCs CNA profiling reveals patterns of frequent convergent alterations developed independently through branched evolution, undetectable through bulk sequencing. Citation Format: Claudio Forcato, Andrea Raspadori, Alberto Ferrarini, Mario Terracciano, Valentina del Monaco, Marianna Garonzi, Carrie Morano, Steven Gross, Genny Buson, Chiara Bolognesi, Francesca Fontana, Gianni Medoro, Mark Connelly, Nicolò Manaresi. Genome-wide copy number profiling of single circulating multiple myeloma cells (CMMCs) reveals intra-patient convergent copy-number alterations (CNAs) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2911.

Hydrogel droplet single-cell processing: DNA purification, handling, release, and on-chip linearization.

The preparation and handling of mammalian single-cell genomic DNA is limited by the complexity bottleneck inherent to performing multi-step, multi-reagent operations in a microfluidic environment. We have developed a method for benchtop preparation of high-molecular weight, intact, single-cell genomes and demonstrate the extraction of long nucleic acid molecules in a microfluidic system. Lymphoblasts are encapsulated inside of alginate microparticles using a droplet microfluidics, and cells are lysed in bulk. The purified genomes are then delivered to and imaged on a dedicated microfluidic device. High-molecular weight DNA is protected from shear and retains its original cellular identity. Using this encapsulation protocol, we were able to extract individual nucleic acid strands on the millimeter scale inside of a microfluidic channel.

New library construction method for single-cell genomes

A central challenge in sequencing single-cell genomes is the accurate determination of point mutations, phasing of these mutations, and identifying copy number variations with few assumptions. Ideally, this is accomplished under as low sequencing coverage as possible. Here we report our attempt to meet these goals with a novel library construction and library amplification methodology. In our approach, single-cell genomic DNA is first fragmented with saturated transposition to make a primary library that uniformly covers the whole genome by short fragments. The library is then amplified by a carefully optimized PCR protocol in a uniform and synchronized fashion for next-generation sequencing. Each step of the protocol can be quantitatively characterized. Our shallow sequencing data show that the library is tightly distributed and is useful for the determination of copy number variations.

Down the Rabbit Hole of Single-Cell Genome Analysis

Chen etal. (2017) demonstrate whole-genome amplification of single-cell genomic DNA using linear nucleic acid amplification. This provides reliable single-nucleotide variation (SNV) detection across the single-cell genome, facilitating an understanding of cell-to-cell similarities and distinctions.

An evaluation of multiple annealing and looping based genome amplification using a synthetic bacterial community

The low biomass in environmental samples is a major challenge for microbial metagenomic studies. The amplification of a genomic DNA was frequently applied to meeting the minimum requirement of the DNA for a high-throughput next-generation-sequencing technology. Using a synthetic bacterial community, the amplification efficiency of the Multiple Annealing and Looping Based Amplification Cycles (MALBAC) kit that is originally developed to amplify the single-cell genomic DNA of mammalian organisms is examined. The DNA template of 10 pg in each reaction of the MALBAC amplification may generate enough DNA for Illumina sequencing. Using 10 pg and 100 pg templates for each reaction set, the MALBAC kit shows a stable and homogeneous amplification as indicated by the highly consistent coverage of the reads from the two amplified samples on the contigs assembled by the original unamplified sample. Although GenomePlex whole genome amplification kit allows one to generate enough DNA using 100 pg of template in each reaction, the minority of the mixed bacterial species is not linearly amplified. For both of the kits, the GC-rich regions of the genomic DNA are not efficiently amplified as suggested by the low coverage of the contigs with the high GC content. The high efficiency of the MALBAC kit is supported for the amplification of environmental microbial DNA samples, and the concerns on its application are also raised to bacterial species with the high GC content.

Abstract IA14: Single-cell sequencing in cancer genomics

Abstract Comprehensive analyses of cancer genomes promise to inform prognoses and precise cancer treatments. A major barrier, however, is inaccessibility of metastatic tissue. A potential solution is to characterize circulating tumor cells (CTCs), but this requires overcoming multiple hurdles. This talk will present a collaborative research program that has generated an integrated process to isolate, qualify, and sequence whole exomes of CTCs with high fidelity, using a census-based sequencing strategy. Power calculations suggest that mapping of >99.995% of the standard exome is possible in CTCs. We validated our process in two prostate cancer patients including one for whom we sequenced CTCs, a lymph node metastasis, and nine cores of the primary tumor. 51 of 73 CTC mutations (70%) were observed in matched tissue. Moreover, we identified 10 early-trunk and 56 metastatic-trunk mutations in the non-CTC tumor samples and found 90% and 73% of these, respectively, in CTC exomes. The talk will also present some technical considerations for enabling single-cell sequencing of tumor cells, including a description of statistical methods to quantitatively assess the amplification bias resulting from whole-genome amplification of single-cell genomic DNA. This study establishes a foundation for CTC genomics in the clinic, and initial applications to patient monitoring will be described. Citation Format: J. Christopher Love. Single-cell sequencing in cancer genomics. [abstract]. In: Proceedings of the AACR Special Conference on Translation of the Cancer Genome; Feb 7-9, 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 1):Abstract nr IA14.

Uniform and accurate single-cell sequencing based on emulsion whole-genome amplification

Whole-genome amplification (WGA) for next-generation sequencing has seen wide applications in biology and medicine when characterization of the genome of a single cell is required. High uniformity and fidelity of WGA is needed to accurately determine genomic variations, such as copy number variations (CNVs) and single-nucleotide variations (SNVs). Prevailing WGA methods have been limited by fluctuation of the amplification yield along the genome, as well as false-positive and -negative errors for SNV identification. Here, we report emulsion WGA (eWGA) to overcome these problems. We divide single-cell genomic DNA into a large number (10(5)) of picoliter aqueous droplets in oil. Containing only a few DNA fragments, each droplet is led to reach saturation of DNA amplification before demulsification such that the differences in amplification gain among the fragments are minimized. We demonstrate the proof-of-principle of eWGA with multiple displacement amplification (MDA), a popular WGA method. This easy-to-operate approach enables simultaneous detection of CNVs and SNVs in an individual human cell, exhibiting significantly improved amplification evenness and accuracy.

Calibrating genomic and allelic coverage bias in single-cell sequencing.

Artifacts introduced in whole-genome amplification (WGA) make it difficult to derive accurate genomic information from single-cell genomes and require different analytical strategies from bulk genome analysis. Here, we describe statistical methods to quantitatively assess the amplification bias resulting from whole-genome amplification of single-cell genomic DNA. Analysis of single-cell DNA libraries generated by different technologies revealed universal features of the genome coverage bias predominantly generated at the amplicon level (1-10 kb). The magnitude of coverage bias can be accurately calibrated from low-pass sequencing (∼0.1 × ) to predict the depth-of-coverage yield of single-cell DNA libraries sequenced at arbitrary depths. We further provide a benchmark comparison of single-cell libraries generated by multi-strand displacement amplification (MDA) and multiple annealing and looping-based amplification cycles (MALBAC). Finally, we develop statistical models to calibrate allelic bias in single-cell whole-genome amplification and demonstrate a census-based strategy for efficient and accurate variant detection from low-input biopsy samples.

Single-Cell Forensic Short Tandem Repeat Typing within Microfluidic Droplets

A short tandem repeat (STR) typing method is developed for forensic identification of individual cells. In our strategy, monodisperse 1.5 nL agarose-in-oil droplets are produced with a high frequency using a microfluidic droplet generator. Statistically dilute single cells, along with primer-functionalized microbeads, are randomly compartmentalized in the droplets. Massively parallel single-cell droplet polymerase chain reaction (PCR) is performed to transfer replicas of desired STR targets from the single-cell genomic DNA onto the coencapsulated microbeads. These DNA-conjugated beads are subsequently harvested and reamplified under statistically dilute conditions for conventional capillary electrophoresis (CE) STR fragment size analysis. The 9-plex STR profiles of single cells from both pure and mixed populations of GM09947 and GM09948 human lymphoid cells show that all alleles are correctly called and allelic drop-in/drop-out is not observed. The cell mixture study exhibits a good linear relationship between the observed and input cell ratios in the range of 1:1 to 10:1. Additionally, the STR profile of GM09947 cells could be deduced even in the presence of a high concentration of cell-free contaminating 9948 genomic DNA. Our method will be valuable for the STR analysis of samples containing mixtures of cells/DNA from multiple contributors and for low-concentration samples.