BLR: a flexible pipeline for haplotype analysis of multiple linked-read technologies.

Linked-read sequencing promises a one-method approach for genome-wide insights including single nucleotide variants (SNVs), structural variants, and haplotyping. We introduce Barcode Linked Reads (BLR), an open-source haplotyping pipeline capable of handling millions of barcodes and data from multiple linked-read technologies including DBS, 10× Genomics, TELL-seq and stLFR. Running BLR on DBS linked-reads yielded megabase-scale phasing with low (<0.2%) switch error rates. Of 13616 protein-coding genes phased in the GIAB benchmark set (v4.2.1), 98.6% matched the BLR phasing. In addition, large structural variants showed concordance with HPRC-HG002 reference assembly calls. Compared to diploid assembly with PacBio HiFi reads, BLR phasing was more continuous when considering switch errors. We further show that integrating long reads at low coverage (∼10×) can improve phasing contiguity and reduce switch errors in tandem repeats. When compared to Long Ranger on 10× Genomics data, BLR showed an increase in phase block N50 with low switch-error rates. For TELL-Seq and stLFR linked reads, BLR generated longer or similar phase block lengths and low switch error rates compared to results presented in the original publications. In conclusion, BLR provides a flexible workflow for comprehensive haplotype analysis of linked reads from multiple platforms.

Characterizing single extracellular vesicles by droplet barcode sequencing for protein analysis.

Small extracellular vesicles (sEVs) have in recent years evolved as a source of biomarkers for disease diagnosis and therapeutic follow up. sEV samples derived from multicellular organisms exhibit a high heterogeneous repertoire of vesicles which current methods based on ensemble measurements cannot capture. In this work we present droplet barcode sequencing for protein analysis (DBS-Pro) to profile surface proteins on individual sEVs, facilitating identification of sEV-subtypes within and between samples. The method allows for analysis of multiple proteins through use of DNA barcoded affinity reagents and sequencing as readout. High throughput single vesicle profiling is enabled through compartmentalization of individual sEVs in emulsion droplets followed by droplet barcoding through PCR. In this proof-of-concept study we demonstrate that DBS-Pro allows for analysis of single sEVs, with a mixing rate below 2%. A total of over 120,000 individual sEVs obtained from a NSCLC cell line and from malignant pleural effusion (MPE) fluid of NSCLC patients have been analyzed based on their surface proteins. We also show that the method enables single vesicle surface protein profiling and by extension characterization of sEV-subtypes, which is essential to identify the cellular origin of vesicles in heterogenous samples.

High throughput barcoding method for genome-scale phasing.

The future of human genomics is one that seeks to resolve the entirety of genetic variation through sequencing. The prospect of utilizing genomics for medical purposes require cost-efficient and accurate base calling, long-range haplotyping capability, and reliable calling of structural variants. Short-read sequencing has lead the development towards such a future but has struggled to meet the latter two of these needs. To address this limitation, we developed a technology that preserves the molecular origin of short sequencing reads, with an insignificant increase to sequencing costs. We demonstrate a novel library preparation method for high throughput barcoding of short reads where millions of random barcodes can be used to reconstruct megabase-scale phase blocks.

Fast and Efficient Fc-Specific Photoaffinity Labeling To Produce Antibody-DNA Conjugates.

Antibody-DNA conjugates are powerful tools for DNA-assisted protein analysis. Growing usage of these methods demands efficient production of high-quality conjugates. We developed an easy and fast synthesis route yielding covalent antibody-DNA conjugates with a defined conjugation site and low batch-to-batch variability. We utilize the Z domain from protein A, containing the unnatural amino acid 4-benzoylphenylalanine (BPA) for photoaffinity labeling of the antibodies' Fc region. Z(xBPA) domains are C-terminally modified with triple-glycine (G3)-modified DNA-oligonucleotides via enzymatic Sortase A coupling. We show reliable modification of the most commonly used IgG's. To prove our conjugates' functionality, we detected antibody-antigen binding events in an assay called Droplet Barcode Sequencing for Protein analysis (DBS-Pro). It confirms not only retained functionality for both conjugate parts but also the potential of using DBS-Pro for quantifying protein abundances. As intermediates are easily storable and our approach is modular, it offers a convenient strategy for screening various antibody-DNA conjugates using the same starting material.

Thermophoresis-Controlled Size-Dependent DNA Translocation through an Array of Nanopores.

Large arrays of nanopores can be used for high-throughput biomolecule translocation with applications toward size discrimination and sorting at the single-molecule level. In this paper, we propose to discriminate DNA length by the capture rate of the molecules to an array of relatively large nanopores (50-130 nm) by introducing a thermal gradient by laser illumination in front of the pores balancing the force from an external electric field. Nanopore arrays defined by photolithography were batch processed using standard silicon technology in combination with electrochemical etching. Parallel translocation of single, fluorophore-labeled dsDNA strands is recorded by imaging the array with a fast CMOS camera. The experimental data show that the capture rates of DNA molecules decrease with increasing DNA length due to the thermophoretic effect of the molecules. It is shown that the translocation can be completely turned off for the longer molecule using an appropriate bias, thus allowing a size discrimination of the DNA translocation through the nanopores. A derived analytical model correctly predicts the observed capture rate. Our results demonstrate that by combining a thermal and a potential gradient at the nanopores, such large nanopore arrays can potentially be used as a low-cost, high-throughput platform for molecule sensing and sorting.

Droplet Barcode Sequencing for targeted linked-read haplotyping of single DNA molecules.

Data produced with short-read sequencing technologies result in ambiguous haplotyping and a limited capacity to investigate the full repertoire of biologically relevant forms of genetic variation. The notion of haplotype-resolved sequencing data has recently gained traction to reduce this unwanted ambiguity and enable exploration of other forms of genetic variation; beyond studies of just nucleotide polymorphisms, such as compound heterozygosity and structural variations. Here we describe Droplet Barcode Sequencing, a novel approach for creating linked-read sequencing libraries by uniquely barcoding the information within single DNA molecules in emulsion droplets, without the aid of specialty reagents or microfluidic devices. Barcode generation and template amplification is performed simultaneously in a single enzymatic reaction, greatly simplifying the workflow and minimizing assay costs compared to alternative approaches. The method has been applied to phase multiple loci targeting all exons of the highly variable Human Leukocyte Antigen A (HLA-A) gene, with DNA from eight individuals present in the same assay. Barcode-based clustering of sequencing reads confirmed analysis of over 2000 independently assayed template molecules, with an average of 753 reads in support of called polymorphisms. Our results show unequivocal characterization of all alleles present, validated by correspondence against confirmed HLA database entries and haplotyping results from previous studies.

Transplanted Bone Marrow-Derived Cells Contribute to Human Adipogenesis.

Because human white adipocytes display a high turnover throughout adulthood, a continuous supply of precursor cells is required to maintain adipogenesis. Bone marrow (BM)-derived progenitor cells may contribute to mammalian adipogenesis; however, results in animal models are conflicting. Here we demonstrate in 65 subjects who underwent allogeneic BM or peripheral blood stem cell (PBSC) transplantation that, over the entire lifespan, BM/PBSC-derived progenitor cells contribute ∼10% to the subcutaneous adipocyte population. While this is independent of gender, age, and different transplantation-related parameters, body fat mass exerts a strong influence, with up to 2.5-fold increased donor cell contribution in obese individuals. Exome and whole-genome sequencing of single adipocytes suggests that BM/PBSC-derived progenitors contribute to adipose tissue via both differentiation and cell fusion. Thus, at least in the setting of transplantation, BM serves as a reservoir for adipocyte progenitors, particularly in obese subjects.

Phasing of single DNA molecules by massively parallel barcoding.

High-throughput sequencing platforms mainly produce short-read data, resulting in a loss of phasing information for many of the genetic variants analysed. For certain applications, it is vital to know which variant alleles are connected to each individual DNA molecule. Here we demonstrate a method for massively parallel barcoding and phasing of single DNA molecules. First, a primer library with millions of uniquely barcoded beads is generated. When compartmentalized with single DNA molecules, the beads can be used to amplify and tag any target sequences of interest, enabling coupling of the biological information from multiple loci. We apply the assay to bacterial 16S sequencing and up to 94% of the hypothesized phasing events are shown to originate from single molecules. The method enables use of widely available short-read-sequencing platforms to study long single molecules within a complex sample, without losing phase information.

Toward Rare Blood Cell Preservation for RNA Sequencing.

Cancer is driven by various events leading to cell differentiation and disease progression. Molecular tools are powerful approaches for describing how and why these events occur. With the growing field of next-generation DNA sequencing, there is an increasing need for high-quality nucleic acids derived from human cells and tissues-a prerequisite for successful cell profiling. Although advances in RNA preservation have been made, some of the largest biobanks still do not employ RNA blood preservation as standard because of limitations in low blood-input volume and RNA stability over the whole gene body. Therefore, we have developed a robust protocol for blood preservation and long-term storage while maintaining RNA integrity. Furthermore, we explored the possibility of using the protocol for preserving rare cell samples, such as circulating tumor cells. The results of our study confirmed that gene expression was not impacted by the preservation procedure (r(2) > 0.88) or by long-term storage (r(2) = 0.95), with RNA integrity number values averaging over 8. Similarly, cell surface antigens were still available for antibody selection (r(2) = 0.95). Lastly, data mining for fusion events showed that it was possible to detect rare tumor cells among a background of other cells present in blood irrespective of fixation. Thus, the developed protocol would be suitable for rare blood cell preservation followed by RNA sequencing analysis.

Parallel barcoding of antibodies for DNA-assisted proteomics.

DNA-assisted proteomics technologies enable ultra-sensitive measurements in multiplex format using DNA-barcoded affinity reagents. Although numerous antibodies are available, nowadays targeting nearly the complete human proteome, the majority is not accessible at the quantity, concentration, or purity recommended for most bio-conjugation protocols. Here, we introduce a magnetic bead-assisted DNA-barcoding approach, applicable for several antibodies in parallel, as well as reducing required reagents quantities up to a thousand-fold. The success of DNA-barcoding and retained functionality of antibodies were demonstrated in sandwich immunoassays and standard quantitative Immuno-PCR assays. Specific DNA-barcoding of antibodies for multiplex applications was presented on suspension bead arrays with read-out on a massively parallel sequencing platform in a procedure denoted Immuno-Sequencing. Conclusively, human plasma samples were analyzed to indicate the functionality of barcoded antibodies in intended proteomics applications.

Magnetic bead assisted labeling of antibodies at nanogram scale.

There are currently several initiatives that aim to produce binding reagents for proteome-wide analysis. To enable protein detection, visualization, and target quantification, covalent coupling of reporter molecules to antibodies is essential. However, current labeling protocols recommend considerable amount of antibodies, require antibody purity and are not designed for automation. Given that small amounts of antibodies are often sufficient for downstream analysis, we developed a labeling protocol that combines purification and modification of antibodies at submicrogram quantities. With the support of magnetic microspheres, automated labeling of antibodies in parallel using biotin or fluorescent dyes was achieved.

A novel method for rapid hybridization of DNA to a solid support.

Here we present a novel approach entitled Magnetic Forced Hybridization (MFH) that provides the means for efficient and direct hybridization of target nucleic acids to complementary probes immobilized on a glass surface in less than 15 seconds at ambient temperature. In addition, detection is carried out instantly since the beads become visible on the surface. The concept of MFH was tested for quality control of array manufacturing, and was combined with a multiplex competitive hybridization (MUCH) approach for typing of Human Papilloma Virus (HPV). Magnetic Forced Hybridization of bead-DNA constructs to a surface achieves a significant reduction in diagnostic testing time. In addition, readout of results by visual inspection of the unassisted eye eliminates the need for additional expensive instrumentation. The method uses the same set of beads throughout the whole process of manipulating and washing DNA constructs prior to detection, as in the actual detection step itself.

Targeted transcript profiling by sequencing.

In this work we present a targeted gene expression strategy employing trinucleotide threading (TnT) amplification and massive parallel sequencing. We have previously shown that TnT combined with array readout accurately monitors expression levels. However, with this detection strategy spurious products go undetected. Accordingly, we adapted the TnT protocol to massive parallel sequencing to acquire an unbiased view of the entire TnT-generated product population. In this manner we investigated the identity of undesired products, their extent at different oligonucleotide:RNA ratios and their effect on the expression levels. We demonstrate that TnT gene expression profiling with massive sequencing readout renders reliable expression data from as low as 3.5 ng of total RNA. Moreover, using 350 ng of total RNA results in only 0.7% to 1.1% undesired products. When lowering the amount of input material, the undesired product fraction increases but this does not influence the expression profiles.

Narrow genetic basis for the Australian dingo confirmed through analysis of paternal ancestry.

The dingo (Canis lupus dingo) is an iconic animal in the native culture of Australia, but archaeological and molecular records indicate a relatively recent history on the continent. Studies of mitochondrial DNA (mtDNA) imply that the current dingo population was founded by a small population of already tamed dogs from Southeast Asia. However, the maternal genetic data might give a unilateral picture, and the gene pool has yet to be screened for paternal ancestry. We sequenced 14,437 bp of the Y-chromosome (Y-chr) from two dingoes and one New Guinea Singing Dog (NGSD). This positioned dingo and NGSD within the domestic dog Y-chr phylogeny, and produced one haplotype not detected before. With this data, we characterized 47 male dingoes in 30 Y-chr single-nucleotide polymorphism sites using protease-mediated allele-specific extension technology. Only two haplotypes, H3 and H60, were found among the dingoes, at frequencies of 68.1 and 31.9 %, respectively, compared to 27 haplotypes previously established in the domestic dog. While H3 is common among Southeast Asian dogs, H60 was specifically found in dingoes and the NGSD, but was related to Southeast Asian dog Y-chr haplotypes. H3 and H60 were observed exclusively in the western and eastern parts of Australia, respectively, but had a common range in Southeast. Thus, the Y-chr diversity was very low, similar to previous observations for d-loop mtDNA. Overall genetic evidence suggests a very restricted introduction of the first dingoes into Australia, possibly from New Guinea. This study further confirms the dingo as an isolated feral dog.

Decoding a substantial set of samples in parallel by massive sequencing.

There has been a dramatic increase of throughput of sequenced bases in the last years but sequencing a multitude of samples in parallel has not yet developed equally. Here we present a novel strategy where the combination of two tags is used to link sequencing reads back to their origins from a pool of samples. By incorporating the tags in two steps sample-handling complexity is lowered by nearly 100 times compared to conventional indexing protocols. In addition, the method described here enables accurate identification and typing of thousands of samples in parallel. In this study the system was designed to test 4992 samples using only 122 tags. To prove the concept of the two-tagging method, the highly polymorphic 2(nd) exon of DLA-DRB1 in dogs and wolves was sequenced using the 454 GS FLX Titanium Chemistry. By requiring a minimum sequence depth of 20 reads per sample, 94% of the successfully amplified samples were genotyped. In addition, the method allowed digital detection of chimeric fragments. These results demonstrate that it is possible to sequence thousands of samples in parallel without complex pooling patterns or primer combinations. Furthermore, the method is highly scalable as only a limited number of additional tags leads to substantial increase of the sample size.

Gene-specific FACS sorting method for target selection in high-throughput amplicon sequencing.

BACKGROUND: In addition to shotgun sequencing, next generation sequencing has been shown to be suitable for deep sequencing of many specific PCR-amplified target genes in parallel. However, unspecific product formation is a common problem in amplicon sequencing since these fragments are difficult to fully remove by gel purification, and their presence inevitably reduces the number of mappable sequence reads that can be obtained in each sequencing run. RESULTS: We have used a novel flow cytometric sorting approach to specifically enrich Roche/454 DNA Capture beads carrying target DNA sequences on their surface, and reject beads carrying unspecific sequences. This procedure gives a nearly three-fold increase in the fraction of informative sequences obtained. Presented results also show that there are no significant differences in the distribution or presence of different genotypes between a FACS-enriched sample and a standard-enriched control sample. CONCLUSIONS: Target-specific FACS enrichment prior to Roche/454 sequencing provides a quick, inexpensive way of increasing the amount of high quality data obtained in a single sequencing run, without introducing any sequence bias.

PCR amplification and genetic analysis in a microwell cell culturing chip.

We have previously described a microwell chip designed for high throughput, long-term single-cell culturing and clonal analysis in individual wells providing a controlled way of studying high numbers of individual adherent or non-adherent cells. Here we present a method for the genetic analysis of cells cultured on-chip by PCR and minisequencing, demonstrated using two human adherent cell lines: one wild type and one with a single-base mutation in the p53 gene. Five wild type or mutated cells were seeded per well (in a defined set of wells, each holding 500 nL of culture medium) in a 672-microwell chip. The cell chip was incubated overnight, or cultured for up to five days, depending on the desired colony size, after which the cells were lysed and subjected to PCR directly in the wells. PCR products were detected, in the wells, using a biotinylated primer and a fluorescently labelled primer, allowing the products to be captured on streptavidin-coated magnetic beads and detected by a fluorescence microscope. In addition, to enable genetic analysis by minisequencing, the double-stranded PCR products were denatured and the immobilized strands were kept in the wells by applying a magnetic field from the bottom of the wells while the wells were washed, a minisequencing reaction mixture was added, and after incubation in appropriate conditions the expected genotypes were detected in the investigated microwells, simultaneously, by an array scanner. We anticipate that the technique could be used in mutation frequency screening, providing the ability to correlate cells' proliferative heterogeneity to their genetic heterogeneity, in hundreds of samples simultaneously. The presented method of single-cell culture and DNA amplification thus offers a potentially powerful alternative to single-cell PCR, with advantageous robustness and sensitivity.

Analysis of short tandem repeats by parallel DNA threading.

The majority of studies employing short tandem repeats (STRs) require investigation of several of these genetic markers. As such, we demonstrate the feasibility of the trinucleotide threading (TnT) approach for scalable analysis of STRs. The TnT method represents a parallel amplification alternative that addresses the obstacles associated with multiplex PCR. In this study, analysis of the STR fragments was performed with capillary gel electrophoresis; however, it should be possible to combine our approach with the massive 454 sequencing platform to considerably increase the number of targeted STRs.

Visual DNA as a diagnostic tool.

We report on the incorporation of the Visual DNA concept in a genotyping assay as a simple and straightforward detection tool. The principle of trapping streptavidin-coated superparamagnetic beads of micrometer size for visualization of genetic variances is used for PrASE-based detection of a panel of mutations in the severe and common genetic disorder of cystic fibrosis. The method allows a final investigation of genotypes by the naked eye and the output is easily documented using a regular hand-held device with an integrated digital camera. A number of samples were run through the assay, showing rapid and accurate detection using superparamagnetic beads and an off-the-shelf neodymium magnet. The assay emphasizes the power of Visual DNA and demonstrates the potential value of the method in future point-of-care tests.