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.

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.

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.

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.

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.

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.

Flow cytometry for enrichment and titration in massively parallel DNA sequencing.

Massively parallel DNA sequencing is revolutionizing genomics research throughout the life sciences. However, the reagent costs and labor requirements in current sequencing protocols are still substantial, although improvements are continuously being made. Here, we demonstrate an effective alternative to existing sample titration protocols for the Roche/454 system using Fluorescence Activated Cell Sorting (FACS) technology to determine the optimal DNA-to-bead ratio prior to large-scale sequencing. Our method, which eliminates the need for the costly pilot sequencing of samples during titration is capable of rapidly providing accurate DNA-to-bead ratios that are not biased by the quantification and sedimentation steps included in current protocols. Moreover, we demonstrate that FACS sorting can be readily used to highly enrich fractions of beads carrying template DNA, with near total elimination of empty beads and no downstream sacrifice of DNA sequencing quality. Automated enrichment by FACS is a simple approach to obtain pure samples for bead-based sequencing systems, and offers an efficient, low-cost alternative to current enrichment protocols.

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.

Generations of sequencing technologies.

Advancements in the field of DNA sequencing are changing the scientific horizon and promising an era of personalized medicine for elevated human health. Although platforms are improving at the rate of Moore's Law, thereby reducing the sequencing costs by a factor of two or three each year, we find ourselves at a point in history where individual genomes are starting to appear but where the cost is still too high for routine sequencing of whole genomes. These needs will be met by miniaturized and parallelized platforms that allow a lower sample and template consumption thereby increasing speed and reducing costs. Current massively parallel, state-of-the-art systems are providing significantly improved throughput over Sanger systems and future single-molecule approaches will continue the exponential improvements in the field.

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.

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.

Expression profiling of signature gene sets with trinucleotide threading.

In recent years, studies have shown that expression profiling of carefully chosen intermediary gene sets, comprising approximately 10 to 100 genes, can convey the most relevant information compared to much more complex whole-genome studies. In this paper, we present a novel method suitable for expression profiling of moderate gene sets in a large number of samples. The assay implements the parallel amplification features of the trinucleotide threading technique (TnT), which encompasses linear transcript-based DNA thread formation in conjunction with exponential multiplexed thread amplification. The amplifications bestow the method with high sensitivity. The TnT procedure together with thread detection, relying on thread-specific primer extension followed by hybridization to universal tag arrays, allows for three distinction levels, thus offering high specificity. Additionally, the assay is easily automated and flexible. A gene set, comprising 18 protein epitope signature tags from the Swedish Human Protein Atlas program, was analyzed with the TnT-based approach and the data were compared with those generated by both real-time PCR and genome-wide cDNA arrays, with the highest correlation observed between TnT and real-time PCR. Taken together, expression profiling with trinucleotide threading represents a reliable approach for studies of intermediary gene sets.

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.

Allelotyping by massively parallel pyrosequencing of SNP-carrying trinucleotide threads.

Here we present an approach for allelotyping combining the multiplexing features of the trinucleotide threading (TnT) method with pooling of genomic DNA and massively parallel pyrosequencing, enabling reliable allele frequency estimation in large cohorts. The approach offers several benefits as compared to array-based methods and allows undertaking highly complex studies without compromising accuracy, while keeping the workload to a minimum. This proof-of-concept study involves formation of trinucleotide threads, targeting a total of 147 single-nucleotide polymorphisms (SNPs) related to obesity and cancer, for multiplex amplification and allele extraction from a pool of 462 genomes, followed by massively parallel pyrosequencing. Approximately 177k reads were approved, identified, and assigned to SNP-carrying threads rendering representative allele frequencies in the cohort.

Screening of common CYP1B1 mutations in Iranian POAG patients using a microarray-based PrASE protocol.

PURPOSE: The gene coding cytochrome P4501B1 (CYP1B1) has been shown to be a major cause of primary congenital glaucoma in the Iranian population. More recently it was shown to also be important in juvenile-onset open angle glaucoma (JOAG). We aimed to further investigate the role of CYP1B1 in a larger cohort of primary open angle glaucoma (POAG) patients which included late-onset patients. We also aimed to set up a microarray based protocol for mutation screening with an intent of using the protocol in a future population level screening program. METHODS: Sixty three POAG patients, nine affected family members, and thirty three previously genotyped primary congenital glaucoma (PCG) patients were included in the study. Clinical examination included slit lamp biomicroscopy, IOP measurement, gonioscopic evaluation, fundus examination, and measurement of perimetry. G61E, R368H, R390H, and R469W were screened by a protocol that included multiplexed allele specific amplification in the presence of a protease (PrASE), use of sequence tagged primers, and hybridization to generic arrays on microarray slides. The entire coding sequences of CYP1B1 and myocilin (MYOC) genes were sequenced in all individuals assessed by the microarray assay to carry a mutation. Intragenic single nucleotide polymorphism (SNP) haplotpes were determined for mutated alleles. RESULTS: Genotypes assessed by the array-based PrASE methodology were in 100% concordance with sequencing results. Seven mutation carrying POAG patients (11.1%) were identified, and their distribution was quite skewed between the juvenile-onset individuals (5/21) as compared to late-onset cases (2/42). Four of the seven mutation carrying Iranian patients harbored two mutated alleles. CYP1B1 mutated alleles in Iranian PCG and POAG patients shared common haplotypes. MYOC mutations were not observed in any of the patients. CONCLUSIONS: The PrASE approach allowed reliable simultaneous genotyping of many individuals. It can be an appropriate tool for screening common mutations in large sample sizes. The results suggest that CYP1B1 is implicated in POAG among Iranians, notably in the juvenile-onset form. Contrary to POAG patients studied in other populations, many mutation harboring Iranian patients carry two mutated alleles. We propose an explanation for this observation.

Tri-nucleotide threading for parallel amplification of minute amounts of genomic DNA.

Efforts to correlate genetic variations with phenotypic differences are intensifying due to the availability of high-density maps of single nucleotide polymorphisms (SNPs) and the development of high throughput scoring methods. These recent advances have led to an increased interest for improved multiplex preparations of genetic material to facilitate such whole genome analyses. Here we propose a strategy for the parallel amplification of polymorphic loci based on a reduced set of nucleotides. The technique denoted Tri-nucleotide Threading (TnT), allows SNPs to be amplified via controlled linear amplification followed by complete removal of the target material and subsequent amplification with a pair of universal primers. A dedicated software tool was developed for this purpose and variable positions in genes associated with different forms of cancer were analyzed using sub-nanogram amounts of starting material. The amplified fragments were then successfully scored using a microarray-based PrASE technique. The results of this study, in which 75 SNPs were analyzed, show that the TnT technique circumvents potential problems associated with multiplex amplification of SNPs from minute amounts of material. The technique is specific, sensitive and can be readily adapted to equipment and genotyping techniques used in other research laboratories without requiring changes to the preferred typing method.

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.

Competitive enzymatic reaction to control allele-specific extensions.

Here, we present a novel method for SNP genotyping based on protease-mediated allele-specific primer extension (PrASE), where the two allele-specific extension primers only differ in their 3'-positions. As reported previously [Ahmadian,A., Gharizadeh,B., O'Meara,D., Odeberg,J. and Lundeberg,J. (2001), Nucleic Acids Res., 29, e121], the kinetics of perfectly matched primer extension is faster than mismatched primer extension. In this study, we have utilized this difference in kinetics by adding protease, a protein-degrading enzyme, to discriminate between the extension reactions. The competition between the polymerase activity and the enzymatic degradation yields extension of the perfectly matched primer, while the slower extension of mismatched primer is eliminated. To allow multiplex and simultaneous detection of the investigated single nucleotide polymorphisms (SNPs), each extension primer was given a unique signature tag sequence on its 5' end, complementary to a tag on a generic array. A multiplex nested PCR with 13 SNPs was performed in a total of 36 individuals and their alleles were scored. To demonstrate the improvements in scoring SNPs by PrASE, we also genotyped the individuals without inclusion of protease in the extension. We conclude that the developed assay is highly allele-specific, with excellent multiplex SNP capabilities.

Comparison of PrASE and Pyrosequencing for SNP Genotyping.

BACKGROUND: There is an imperative need for SNP genotyping technologies that are cost-effective per sample with retained high accuracy, throughput and flexibility. We have developed a microarray-based technique and compared it to Pyrosequencing. In the protease-mediated allele-specific extension (PrASE), the protease constrains the elongation reaction and thus prevents incorrect nucleotide incorporation to mismatched 3'-termini primers. RESULTS: The assay is automated for 48 genotyping reactions in parallel followed by a tag-microarray detection system. A script automatically visualizes the results in cluster diagrams and assigns the genotypes. Ten polymorphic positions suggested as prothrombotic genetic variations were analyzed with Pyrosequencing and PrASE technologies in 442 samples and 99.8 % concordance was achieved. In addition to accuracy, the robustness and reproducibility of the technique has been investigated. CONCLUSION: The results of this study strongly indicate that the PrASE technology can offer significant improvements in terms of accuracy and robustness and thereof increased number of typeable SNPs.