A simple bead-based method for generating cost-effective co-barcoded sequence reads

doi:10.1038/protex.2018.116

Method Article

A simple bead-based method for generating cost-effective co-barcoded sequence reads

https://doi.org/10.1038/protex.2018.116

This work is licensed under a CC BY 4.0 License

This protocol has been posted on Protocol Exchange, an open repository of community-contributed protocols sponsored by Nature Portfolio. These protocols are posted directly on the Protocol Exchange by authors and are made freely available to the scientific community for use and comment.

Version 1

posted 11 Oct, 2018

You are reading this latest protocol version

Advances in genome sequencing continue to decrease cost and improve quality. However, the diploid information found in most genomes is still lost by the most common sequencing methods. There are several methods currently available for obtaining this haplotype information. However, most of these methods are either time consuming, require expensive equipment and reagents, or both. Here we describe single tube long fragment read (stLFR), a simple barcoded bead-based process capable of near perfect whole genome variant calling and haplotyping. Our method uses equipment available in nearly all molecular biology laboratories and costs approximately 30 dollars a sample. In addition, the data generated from this process is capable of detecting and phasing structural variations and scaffolding contigs. We expect in the future that this process will enable affordable diploid de novo assembly.

Genetics

Biotechnology

Molecular Biology

stLFR

co-barcoded reads

de novo assembly

haplotyping

NGS

There are currently many methods now available for generating haplotype information (Zhang et al. 2006; Ma et al. 2010; Fan et al. 2011; Kitzman et al. 2011; Suk et al. 2011; Duitama et al. 2012; Peters et al. 2012; Selvaraj et al. 2013; Amini et al. 2014; Kuleshov et al. 2014; Zheng et al. 2016; Zhang et al. 2017). The majority are based on co-barcoding short sequencing templates originating from long genomic DNA molecules (Peters et al. 2014), a process first described by Drmanac (Drmanac 2006). We first enabled co-barcoding with the development of long fragment read (LFR). While the method (McElwain et al. 2017) was quite useful and resulted in the complete phasing of hundreds of genomes (Peters et al. 2012; Schaaf et al. 2013; Peters et al. 2015; Ciotlos et al. 2016; Hellner et al. 2016; Mao et al. 2016; Gulbahce et al. ; Walker et al. 2017; Mao et al. 2018), including over 100 from the personal genome project (Mao et al. 2016), and the accurate sequencing from difficult samples where as few as five cells were available (Peters et al. 2012; Peters et al. 2015; Hellner et al. 2016; Gulbahce et al. 2017; Mao et al. 2018), it had several limitations. Primarily, each sample was split across a 384 well plate, so processing multiple samples at a time meant handling multiple plates. In addition, because the process relied on multiple displacement amplification, it had base calling bias in regions of high GC base content and sequencing coverage was more variable necessitating _{100X coverage to achieve high quality whole genome variant calling. Finally, because only 384 compartments were used, this resulted in}10% of a haploid genome per barcode and limited the phasing N50 length, made it more difficult to detect structural variations, and made de novo assembly impossible.

Here we introduce a method called single tube long fragment read (stLFR). This process was first proposed by us in a patent application (Drmanac 2013) and a similar, but limited implementation of this method was published by Zhang et al. (Zhang et al. 2017). Like the LFR process described above stLFR utilizes the concept of co-barcoding short reads belonging to a long DNA fragment. Unlike LFR, the stLFR process uses the surface of micron sized magnetic beads as compartments and as the name suggests this process take place in a single tube. In a typical stLFR experiment 10-50 million beads are used. Each bead contains hundreds of thousands of copies of the same barcode and there are a total of 3.6 billion different barcodes. This means in an stLFR sample there is practically no overlap between barcode sequences. Also, because only 1-10 ngs of genomic DNA is used, most beads capture less than 10 long DNA molecules. For de novo assembly applications it is possible to achieve over 80% of long DNA fragments having sequencing reads co-barcoded by a unique barcode (captured by a single bead) by using 1 ng of DNA and 50 million beads.

stLFR is a simple process that can be performed without any special equipment. The process begins by inserting a transposon with a capture sequence at a regular interval along the genomic DNA molecules. These transposons inserted DNA molecules are then added to barcode labeled 3 μM magnetic beads that contain a complementary sequence to the capture sequence on the transposon. After hybridization, the barcode is transferred to the transposon sequence through a ligation step. At this point two different methods can be used. The first and the simplest is to directly PCR amplify using primers to the transposon sequences followed by analysis on a next generation sequencing (NGS) platform. If more coverage of each original long genomic DNA molecule is desired an additional ligation step with another adapter sequence can be performed followed by PCR and sequencing. Both methods are described in detail in the following protocol. A full analysis of stLFR data was performed by Wang et al. (Wang et al. 2018).

1 Kb Plus DNA Ladder (ThermoFisher, cat. no. 10787018)

100Kd MWCO Biotech CE dialysis tubing (Spectrum Labs, cat. no. 131486)

384-well Armadillo PCR plate (ThermoFisher, cat. no. AB2384)

Agencourt® AMPure XP beads (Beckman Coulter, cat. no. A63882)

APE 1 (10,000 units/mL) (New England Biolabs, cat. no. M0282L)

ATP (100 mM) (Teknova, cat. no. A1210)

Barcoded bead construction oligonucleotides (IDT) (see note)

Betaine (5 M) (Sigma-Aldrich, cat. no. B0300-5VL)

BSA (20mg/mL) (New England Biolabs, cat. no. B9000S)

Common adapter oligonucleotides (IDT)

DMF (~100%) (Sigma-Aldrich, cat. no. D4551-250ML)

DMSO (100%) (Sigma-Aldrich, cat. no. D9170-5VL)

dNTPs (25 mM) (ThermoFisher, cat. no. R1121)

Dialysis tubing (1,000 kD MWCO) (Spectrum Laboratories, Inc., cat. no. 131486)

DTT(Sigma-Aldrich, cat. no. 11583786001)

Dynabeads™ M-280 Streptavidin (ThermoFisher, cat. no. 60210)

EDTA (0.5 M, pH 8.0) (Sigma-Aldrich, cat. no. 03690-100ML)

ET SSB (500 μg/ml) (New England Biolabs, cat. no. M2401S)

Exonuclease I (20,000 units/mL) (New England Biolabs, cat. no. M0293L)

Exonuclease III (100,000 units/mL) (New England Biolabs, cat. no. M0206L)

Formamide (100%, 250 mL) (Sigma-Aldrich, cat. no. 47671-250ML-F)

Glycerol (100%) (Sigma-Aldrich, cat. no. G5516-100ML)

KCl (Sigma-Aldrich, cat. no. P9333-1KG)

KH2PO4 (Sigma-Aldrich, cat. no. 795488-1KG)

KOH (Sigma-Aldrich, cat. no. P5958-1KG)

MgCl2 (1 M) (Sigma-Aldrich, cat. no. 63069-500ML)

MgSO4 (1 M) (Sigma-Aldrich, cat. no. M3409-100ML)

MicroAmp Clear Adhesive Film (ThermoFisher, cat. no. 4306311)

NaCl (5M) (ThermoFisher, cat. no. AM9760G)

Na2HPO4 (Sigma-Aldrich, cat. no. S7907-1KG)

NaOH (10M) (Sigma-Aldrich, cat. no. 72068-100ML)

NEB2 Buffer (10X) (New England Biolabs, cat. no. B7002S)

PEG-8000 (50%) (Rigaku, cat. no. 1008063)

Pfu Turbo Cx Hotstart DNA Polymerase (Agilent, cat. no. 600414)

Proteinase K, recombinant, PCR grade solution (14-22 mg/mL) (Roche, cat. no. 03115844001)

RiboRuler Low Range RNA Ladder (Thermofisher, cat. no. SM1831)

RNase-IT ribonuclease cocktail (Agilent, cat. no. 400720)

SDS (10%) (ThermoFisher, cat. no. 15553027)

Sucrose (Sigma-Aldrich, cat. no. S7903-1KG)

T4 DNA ligase (2x106 units/mL) (New England Biolabs, cat. no. M0202M)

TA Buffer (10X) (Teknova, cat. no. T0379)

TAPS-NaOH (1 M, pH 8.5) (Boston BioProducts, cat. no. BB-2375)

TBE (10X) (ThermoFisher, cat. no. 15581028)

TE Buffer (10X) (Fisher Scientific, cat. no. BP24771)

Tn5 enzyme (see note)

Transposon oligonucleotides (IDT)

Tris-HCl (1 M, pH 7.5) (ThermoFisher, cat. no. 15567027)

Tris-HCl (2 M, pH 7.8) (Amresco, cat. no. J837-500ML)

Triton™ X-100 (10%) (Sigma-Aldrich, cat. no. 93443-100ML)

TWEEN® 20 (10%) (Roche, cat. no. 11332465001)

UDG (5,000 units/mL) (New England Biolabs, cat. no. M0280L)

Reagent Setup

3’ Branch Ligation Buffer (3X)

6 mL of 50% PEG-8000

0.75 mL 2 M Tris-HCl (pH 7.8)

0.3 mL 1 M MgCl2

0.3 mL 0.1 M ATP

15 μL 1 M DTT

75 μL 20 mg/mL BSA

2.560 mL of sterile dH2O

Store at -20 °C for 1 year.

Annealing Buffer (3X)

3 mL of 1 M Tris-HCl, pH 7.5

6 mL of 5 M NaCl

91 mL of sterile dH2O

Store at room temperature for 1 year.

Buffer D (10X)

224 mg of KOH

50 μL of 0.5 M EDTA

2.45 mL of sterile dH2O

Make aliquots and store at -20 °C for 1 month.

Coupling Buffer (1X)

5 mL of 1X TE

5 mL of 100% glycerol

Store at -20 °C for 1 year.

Digestion Buffer (1X, pH 8.0)

1.75 g of Na2HPO4

0.2 g of KCl

0.2 g of KH2PO4

27.4 mL 5 M NaCl

20 mL of 0.5 M EDTA (pH 8.0)

800 mL of sterile dH2O

Adjust the pH to 8.0 with 1 M NaOH.

Add sterile dH2O to a final volume of 1 liter.

Filter sterilize.

Store at room temperature for 1 year.

High-Salt Bead Binding Buffer (1X)

5 mL of 1 M Tris-HCl (pH7.5)

6 mL of 5 M NaCl

20 μL of 0.5 M EDTA

88.98 mL of sterile dH2O

Store at room temperature for 1 year.

High-Salt Wash Buffer (1X)

5 mL of 1 M Tris-HCl, pH 7.5

10 mL of 5 M NaCl

20 μL of 0.5 M EDTA

0.5 mL of 10% TWEEN® 20

84.48 mL of sterile dH2O

Store at room temperature for 1 year.

Hybridization Buffer (1X)

50 mL of 1 M Tris-HCl, pH 7.5

100 mL of 1 M MgCl2

5 mL of 10% TWEEN® 20

845 mL of Water

Store at room temperature for 1 year

Ligation Buffer (10X)

25 mL of 50% PEG-8000

12.5 mL of 2 M Tris-HCl (pH 7.8)

5 mL of 100 mM ATP

5 mL of 1 M MgCl2

2.5 mL of sterile dH2O

Store at -20 °C for 1 year.

Ligation Buffer, No MgCl2 (10X)

25 mL of 50% PEG-8000

12.5 mL of 2 M Tris-HCl (pH 7.8)

5 mL of 100 mM ATP

5 mL of 1 M DTT

2.5 mL of sterile dH2O

Store at -20 °C for 1 year.

Low-Salt Wash Buffer (1X)

5 mL of 1 M Tris-HCl, pH 7.5

3 mL of 5 M NaCl

0.5 mL of 10% TWEEN® 20

91.5 mL of sterile dH2O

Store at room temperature for 1 year.

Lysis Buffer (1X, pH 8.3)

0.22 g of KCl

120 g of sucrose

13 mL of 1 M Tris-HCl (pH 7.5)

2 mL of 0.5 M EDTA (pH 8.0)

28 mL of 5 M NaCl

10 mL of Triton® X-100

800 mL of sterile dH2O

Adjust the pH to 8.3

Add sterile dH2O to a final volume of 1 liter

Filter sterilize

Store at 4 °C for 1 year.

SSB binding buffer (1X)

0.5 mL of 1 M Tris-HCl (pH7.5)

0.2 mL of 1 M MgCl2

9.3 mL of sterile dH2O

Store at -20 °C for 1 year.

Transposase Buffer (5X)

0.5 mL of 1 M TAPS-NaOH (pH 8.5)

0.25 mL of 1 M MgCl2

5 mL of 100% DMF

4.25 mL of sterile dH2O

Store at -20 °C for 1 year.

PfuCx mix (2X)

2 mL of 10X PfuCx buffer (included with enzyme)

0.5 mL of 100% DMSO

2 mL of 5 M Betaine

60 μL of 1 M MgSO4

240 μL of 25 mM dNTPs

5.2 mL of sterile dH2O

Notes

Tn5 enzyme preparation

Tn5 was prepared following the protocol outlined in Picelli et al. (Picelli et al. 2014). Tn5 enzyme is also commercially available from Lucigen (cat. no. TNP92110) but requires the acceptance of Illumina’s license acknowledgment.

Barcoded bead construction oligos

Oligonucleotide sequences are available in Supplementary Materials. All barcoded oligonucleotides were synthesized at 100 nmol scale in 384 well format with standard desalting and delivered at a concentration of 200 μM in 1X TE (pH 8.0) by Integrated DNA Technologies (Coralville, IA). There are a total of 1,536 unique barcode oligos for each barcode set and there are 3 barcode sets. This enables up to ~3.6 billion different barcode combinations. This may be more than necessary for some applications and less barcode combinations can be achieved by ordering less plates of oligonucleotides. This particular design does require that at least one barcode oligo from each set is used to create the proper final sequence, however, slight modifications of the 6 base overlapping sequences between barcode sets can be made to remove an entire barcode set.

2.4 L tall polystyrene container (Click Clack, cat. no. 659030) or equivalent

DynaMag™-2 Magnet (ThermoFisher, cat. no. 12321D)

Easy 50 EasySep™ Magnet (Stem Cell Technologies, cat. no. 18002) or equivalent

Lab oven capable of holding Tube Revolver/Rotator

Magnetic plate stirrer

Medium sized magnetic stir bar

Standard lab vortexer

Tetrad PCR thermocycler (Bio-Rad, cat. no. PTC0240) or equivalent capable of 100 μL per well reaction volumes

Tube Revolver/Rotator (Thermo Fisher, cat. no. 88881001) or equivalent

High Molecular Weight DNA isolation from cells

This method is based on the RecoverEase™ DNA isolation kit protocol (Agent Technologies 2015), but is performed using much larger volumes so as to reduce the viscosity of the resulting solution.

Pellet up to 1 × 107 dispersed nucleated cells in a 15 or 50 mL conical tube (500xg for 5 min). Remove supernatant. Add 500 μL of lysis buffer to the cell pellet and vortex sample briefly for 3–5 seconds on medium speed and place the conical tube in refrigerator for ~10 minutes, swirling occasionally.
Prepare proteinase K solution by combining 250 μL of 10% SDS, 250 μL of Proteinase K, and 4 mL of 1X TE. Place on 50 °C heat block and warm briefly (~5 minutes).
Prepare the digestion solution by combining 20 μL of RNase-It ribonuclease cocktail with 4 mL of digestion buffer.
Add ~4 mL of the prepared digestion solution to the lysed cells and buffer from step one and gently rock the conical tube.
Place the conical tube in a 50 °C heat block after 5 minutes add 4.5 mL of the warmed proteinase K solution to the free-floating pellet. Swirl the conical tube gently to mix.
Recap the tube and incubate in a 50 °C heat block for 2 hours, swirling the tube gently every 30 minutes.
Cut approximately 13 cm of dialysis tubing (it has a capacity of approximately 1 mL/cm). Allow to equilibrate in 0.5X TE for 30 minutes. Seal one end with a dialysis clip.
Pour at least 1 L of 0.5X TE buffer into a dialysis reservoir.
Carefully pour viscous genomic DNA from the conical tube into the open end of dialysis tubing. Seal open end of dialysis tubing with dialysis clip. Attach float to one clip. Place dialysis tubing with float into dialysis reservoir.
Dialyze the genomic DNA at room temperature for 24 to 48 hours while stirring the buffer gently with a magnetic stir bar. Replace the TE buffer once during the dialysis period to maximize the purity of the recovered DNA.
Upon completion of dialysis, remove the dialysis tubing from TE buffer, remove the float and clip from the top of the dialysis tubing and gentle pour into a 15 mL conical tube. DNA can be used immediately without shearing.

Barcoded beads

Barcoded beads are constructed using a split and pool strategy with 3 sets of double-stranded barcode DNA molecules. Full length adaptors are constructed through successive ligations (Figure 1 and Supplementary Figure 1). Barcode oligonucleotides are supplied in 384-well plates (see Reagents Notes). Common adapter oligo nucleotides are supplied in tubes. Depending on what sequencing technology is being used, it may be necessary to alter the PCR primer sequence within the common adapter oligo.

Mix 10 μL of complementary oligonucleotide from each well of the source 384-well plates in 384-well PCR plates with 10 μL 3X Annealing Buffer. Mix 30 μL of common adapter oligonucleotides in one well of an 8-wel PCR strip tube.
Incubate at 70°C for 3 minutes followed by a slow ramp of 0.1°C/s to 20°C on a PCR thermocycler. Hybridized barcode oligonucleotides have a final concentration of 66 μM.
Mix 4.725 mL (157.5 μmol) of hybridized Bead Linker containing a 5’ dual-biotin with 3.225 mL of Ligation Buffer (10X), 460.8 μL (921,600 units) of T4 DNA Ligase, and 9.67 mL dH2O to a total volume of 18.081 mL.
Dispense 11.2 μL of the ligation mixture into each well of four new 384-well PCR plates. Then add 8.8 μL (580 pmol) from each well of the hybridized first barcode plates to each well containing the bead linker ligation mixture. Seal with MicroAmp Clear Adhesive Film, vortex, centrifuge, and incubate at room temperature for 1 hour.
Collect 100 billion (143 mL) of M-280 streptavidin coated magnetic beads by transferring 50 mL beads into empty 50 mL centrifuge tube. Place the 50 mL tube with beads in the Easy 50 EasySep™ Magnet for 5 minutes to collect beads to the side of the tube. Carefully remove supernatant by pipette. Transfer a second 50 mL beads into tube on magnet. Let sit for 5 minutes on magnet and carefully remove supernatant. Transfer a final 43 mL of beads to the 50 mL tube. Let sit for 5 minutes on magnet and carefully remove supernatant. Wash beads twice with Low-Salt Wash Buffer then resuspend well in 8 mL of High-Salt Bead Binding Buffer.
Dispense 5 μL of beads in High-Salt Bead Binding Buffer to each well of the plates containing ligation product. Vortex the beads source tube occasionally during dispensing to keep the beads well-suspended.
Seal plates with MicroAmp Clear Adhesive Film, vortex, and place onto the tube revolver for incubation at room temperature for 1 hour on “oscillating” mode.
Centrifuge plates at 300 x g for 5 seconds to remove beads from seal, but not allow a pellet to form. Remove seal and add 2.8 μL 0.1% SDS to each well. Seal plates again with MicroAmp Clear Adhesive Film, vortex briefly and incubate at room temperature for 10 minutes.
Vortex and then centrifuge plates at 300 x g for 5 seconds to remove beads from the plate seal. Remove the seal from each plate and invert plates onto a collection tray. Centrifuge at 500 x g for 2 minutes. Using a 10 mL serological pipette, collect beads into one new 50 mL tube.
Collect beads to the side of the tube on the Easy 50 EasySep™ Magnet for 5 minutes. Discard supernatant. Wash once with 10 mL of High Salt Wash Buffer and then twice with Low-Salt Wash Buffer. Resuspend beads 8 mL 1XLigation Buffer.
Dispense 5ul of beads into each well of four new 384-well PCR plates. Vortex the beads source tube occasionally during dispensing to keep the beads well-suspended.
To ligate the second set of barcodes, make a mixture containing 3.225 mL Ligation Buffer (10X), 460.8 μL (921,600 units) of T4 DNA Ligase, and 6.33 mL dH2O to a total volume of 10.02 mL. Dispense 6.2 μL of the second ligation mixture to each well of the four 384-well PCR plates containing beads. Next add 8.8 μL (580pmol) from each well of the hybridized second barcode plates to the corresponding wells of the 384-well PCR plates containing the bead and ligation mixture.
Repeat steps 18-22.
To ligate the third set of barcodes, make a ligation mixture containing 3.225 mL Ligation Buffer (10X), 460.8 μL (921,600 units) of T4 DNA Ligase, and 6.33 mL dH2O to a total volume of 10.02 mL. Dispense 6.2 μL of the third ligation mixture to each well of the four 384-well PCR plates containing beads. Next add 8.8 μL (580pmol) from each well of the hybridized third barcode plates to the corresponding wells of the 384-well PCR plates containing the bead and ligation mixture.
Repeat steps 18-22. The beads can now be stored at 4 °C for up to one year. In the current form the beads are almost completely double stranded and not yet in the correct form to be used for stLFR.
Count the beads with hemocytometer and take out 5 million beads for the QC step. Place the tube with beads onto the DynaMag™-2 Magnet for 5 minutes. Discard supernatant. Add 5 μL 100% formamide, 4 μL dH2O and 1 μL 10X loading buffer. Incubate it at 95 °C for 3 minutes on a PCR thermocycler. Immediately place on ice for 2 minutes. Place the tube with beads onto the DynaMag™-2 Magnet for 5 minutes. Collect supernatant, load on a 15% TBU gel, and run at 200 V for 40 minutes to check the oligo length and amount (Figure 2). Alternatively, beads can be examined using a flow cytometer by hybridizing a fluorescently labeled oligo to the 3’ end of the bead adapter sequence. We typically see about 25% of the total streptavidin bound sites have a full length constructed adapter sequence.

Bead preparation for stLFR

To prepare beads for stLFR, they must first be denatured to single stranded DNA and then rehybridized with the bridge oligo.

Pipette 500 million constructed barcoded beads from step 26. of the previous section into a standard 1.5 mL microcentrifuge tube.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant.
Add 1 mL of a 1X dilution of Buffer D. Vortex briefly and incubate for 2 minutes at room temperature.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant.
Repeat steps 30 and 31 one more time.
Wash once in 1X Annealing buffer. Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant.
Mix 36 μL of 100 μM Bridge Oligo, 333.33 μL of Annealing Buffer (3X), and 630.67 μL of dH2O for a final volume of 1 mL. Add to mixture to beads. Vortex briefly.
Incubate at 60 °C for 5 minutes, and room temperature for 50 minutes.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and resuspend in 500 μL of Low-Salt Wash Buffer. These beads are now ready for stLFR and can be stored for 3 months at 4 °C.

Two transposons stLFR protocol

The protocol utilizes two transposons to create hybridization sequences and PCR primer sites along the length of genomic DNA molecules. This is the most simplified and quickest stLFR method, but has potentially 50% less coverage per long DNA fragment than the 3’ branch ligation protocol. It may be necessary to alter some of the transposon sequence after the mosaic region for compatibility with sequencing technologies other than BGISEQ-500. Check the sequencing primers being used prior to ordering these oligonucleotides. Information on all of the oligonucleotide sequences is available in the supplementary materials.

Hybridize the capture transposon oligos by combining 10 μL of Transposon1T (100 μM), 10 μL of TransposonB (100 μM), 10 μL of Annealing Buffer (3X) in the first well of an 8 well PCR strip tube and the non-captured transposon oligos by combining 10 μL of Transposon1T (100 μM), 10 μL of TransposonB (100 μM), 10 μL of Annealing Buffer (3X) in the second well of the same PCR strip tube.
Incubate at 70 °C for 3 minutes followed by a slow ramp of 0.1 °C/s to 20 °C on a PCR thermocycler. Combine the two transposons into the third well of the PCR strip tube.
Couple the Tn5 enzyme to the transposon mix by combining 9.6 μL of mix transposons with 23.53 μL of Tn5 (13.6 pmol/μL), and 46.87 μL of Coupling Buffer (1X).
Incubate at 30 °C for 1 hour. Use immediately or store at -20 °C for up to 1 month. For optimal performance and consistency between experiments we suggest making aliquots prior to storage.
Incorporate transposons into long genomic DNA by combining 12 μL of Transposase Buffer (5X), 0.5 μL of coupled transposon from step 40, and 40 ng of DNA in a total volume of 60 μL in one well of an 8-well strip tube. Note: this amount of DNA and the amount of coupled transposon can be adjusted in this step. It will be necessary to titrate the amount of Tn5 enzyme used as there can be variability between batches. Also, starting with less DNA is possible, but for the purposes of titration it is useful to use 40 ng so that some of the material can be run on an agarose gel to determine the efficiency of transposon incorporation (see later steps).
Incubate at 55 °C for 10 minutes.
Transfer 40 μL of transposon incorporated material to one well of a new 8-well strip tube. Add 4 μL of 1% SDS and incubate at room temperature for 10 minutes.
Load the material from step 43 on a 0.5X TBE 1% agarose gel and run at 150 V for 40 minutes. The transposed DNA should run between 200 to 1,500 bp on the gel (Figure 3). We typically want to see the brightest part of the DNA smear around 600 bp, this might be different based on which sequencing technology is chosen. We typically load controls that are put through the same steps but lack the transposon, the Tn5 enzyme, or genomic DNA. If the size of transposon integrated products looks correct proceed to step 45. If not, repeat the steps above but adjust the concentration of the coupling product until the smear is the desired size.
Dilute 1.5 μL of the remaining product of step 42, with 248.5 μL of 1x Hybridization buffer.
Transfer 50 μL of beads (50 million) from step 36 to a 1.5 mL microcentrifuge tube. Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and resuspend in 250 μL of Hybridization Buffer (1X).
Heat DNA and beads separately at 60 °C for 30 seconds.
Add 250 μL of diluted DNA to the 250 μL of beads, mix gently by flicking the bottom of the tube with a finger, and continue incubating at 60 °C for 10 minutes. Lightly mix the tube every few minutes with your finger.
Place on tube revolver for incubation in oven at 45 °C for 50 minutes on “oscillating” mode.
Make ligation mix by combining 100 μL of Ligation Buffer, No MgCl2 (10X), 2 μL of T4 DNA ligase (2x106 units/mL), and 398 μL of dH2O. Remove tube from rotator and add ligation mix for a total volume of 1 mL.
Incubate on tube revolver for 1 hour on “oscillating” mode at room temperature.
Add 110 μL of 1% SDS to tube and incubate for 10 minutes at room temperature.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and wash once with 500 μL of Low-Salt Wash Buffer and once with 500 μL of NEB2 buffer (1X).
Make capture oligo digestion mix by combining 10 μL of NEB2 buffer (10X), 2 μL of UDG (5,000 U/mL), 3 μL of APE1 (10,000 U/mL), 2 μL of Exonuclease 1 (20,000 units/mL), and 83 μL of dH2O. Remove wash buffer and add digestion mix to beads.
Vortex lightly to resuspend beads and incubate at 37 °C for 30 minutes.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and wash once with 500 μL of Low-Salt Wash Buffer and once with 500 μL of PfuCx Buffer (1X).
Prepare PCR master mix by adding 150 μL of PCR mix (2X), 4 μL of PCR Primer 1 (100 μM), 4 μL of PCR primer 2 (100 μM), 6 μL of PfuCx enzyme, and 136 μL of dH2O. Preheat the PCR mastermix at 95 °C for 3 minutes. Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove wash buffer and add PCR master mix to beads.
Vortex lightly to resuspend beads and cycle PCR reaction with the following conditions:

Step 1 72 °C 10 minutes

Step 2 95 °C 10 seconds

Step 3 58 °C 30 seconds

Step 4 72 °C 2 minutes

Step 6 Repeat steps 2-5 10-12 times

PCR should result in ~500 ng of DNA, run 20 ng of product on a 0.5X TBE 1% agarose gel for 40 minutes at 150 V. The material should be a smear with a peak around 500 bp (Figure 4).
Purify PCR product with 300 μL of Agencourt XP beads following the manufacturer’s protocol. This purified product is now ready to enter the sequencing process.

Single transposon 3’ branch ligation stLFR protocol

This protocol is based on the single transposon insertion and novel adapter ligation methods in a DNA gap (ref) and can enable higher coverage per fragment, which may be important for some sequencing strategies such as de novo assembly. This strategy is slightly more expensive due to additional reagents. It also takes 2.5 hours longer.

Hybridize the capture transposon oligos by combining 10 μL of Transposon1T (100 μM), 10 μL of TransposonB (100 μM), 10 μL of Annealing Buffer (3X) in the first well of an 8 well PCR strip tube and the 3’ branch ligation adapter by combing 10 μL of BranchT (100 μM), 10 μL of BranchB (100 μM), 10 μL of Annealing Buffer (3X) in the second well of the same PCR strip tube.
Incubate at 70 °C for 3 minutes followed by a slow ramp of 0.1 °C/s to 20 °C on a PCR thermocycler.
Couple the Tn5 enzyme to the transposon by combining 9.6 μL of hybridized capture transposon in step 61 with 23.53 μL of Tn5 (13.6 pmol/μL), and 46.87 μL of Coupling Buffer (1X).
Incubate at 30 °C for 1 hour. Use immediately or store at -20 °C for up to 1 month.
Follow steps 41-51.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and wash once with 500 μL of Low-Salt Wash Buffer.
Make adapter oligo digestion mix by combining 10 μL of TA Buffer (10X), 4.5 μL of Exonuclease I (20,000 U/mL), 1 μL of Exonuclease III (100,000 U/mL), and 74.5 μL of dH2O. Remove wash buffer and add digestion mix to beads.
Vortex lightly to resuspend beads and incubate on the tube revolver for 10 minutes at 37 °C on “oscillating” mode.
Add 11 μL of 1% SDS and incubate for 10 minutes at room temperature.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and wash once with 500 μL of Low-Salt Wash Buffer and once with 500 μL of NEB2 buffer (1X).
Make capture oligo digestion mix by combining 10 μL of NEB2 buffer (10X), 2 μL of UDG (5,000 U/mL), 3 μL of APE1 (10,000 U/mL), and 85 μL of dH2O. Remove wash buffer and add digestion mix to beads.
Lightly vortex to resuspend beads and incubate at 37 °C for 30 minutes.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and wash once with 500 μL of High-Salt Wash Buffer and once with 500 μL of Low-Salt Wash Buffer (1X).
Prepare SSB binding mix by combining 20 μL of SSB binding buffer (1X) and 4 μL of ET SSB (500 μg/mL). Remove wash buffer and add SSB binding mix to beads.
Vortex lightly to resuspend beads and incubate on the tube revolver for 30 minutes at 37 °C on “oscillating” mode.
Prepare 3’ branch ligation mix by combining 33.4 μL of 3’ Branch Ligation Buffer (3X), 18 μL of the 3’ branch ligation adapter (16.7 μM) prepared in step 61, 2 μL of T4 DNA ligase (2x106 units/mL), and 46.6 μL of dH2O. Remove wash buffer and add ligation mix to beads.
Vortex lightly to resuspend beads and incubate on the tube revolver for 2 hours at 25 °C on “oscillating” mode.
Place on DynaMag™-2 Magnet for 2 minutes to collect beads onto the side of the tube. Remove supernatant and wash once with 500 μL of High-Salt Wash Buffer and once with 500 μL of PCR buffer (1X).
Prepare PCR master mix by adding 150 μL of 2X PCR buffer, 4 μL of PCR Primer 1 (100 μM), 4 μL of PCR primer 2 (100 μM), 6 μL of PCR enzyme, and 136 μL of dH2O. Remove wash buffer and add PCR master mix to beads.
Vortex lightly to resuspend beads and cycle PCR reaction with the following conditions:

Step 1 95 C 3 minutes

Step 2 95 C 10 seconds

Step 3 58 C 30 seconds

Step 4 72 C 2 minutes

Step 5 Repeat steps 2-4 10-12 times

Follow steps 59-60 above.

Analyzing stLFR data

The starting point for this process is a FASTQ file. This is a standard format for read data that is generated by most sequencing technologies. The software we use to deconvolute the barcode information takes the FASTQ file and expects 42 bases of the barcode and common adapter sequence to be appended to the end of the first read. It matches the barcode read data to the expected 1536 sequences at each barcode position. The barcoding strategy used by stLFR enables error correction of barcodes that have a single base mismatch. The final output from our software is a FASTQ file with the barcode information appended to the end of the read ID with the format #Barcode1ID_Barcode2ID_Barcode3ID, where BarcodeID is a number from 0-1536. Zero for a barcode ID means it did not match any of the expected barcode sequences. The software can be downloaded here https://github.com/stLFR/stLFR_read_demux. We recommend using BWA-mem (Li and Durbin 2009) for mapping, GATK (McKenna et al. 2010) for variant calling, and HapCUT2 (Edge et al. 2017) for phasing. We also recommend mapping to Hg19 with decoy sequences.

Agent Technologies I. 2015. RecoverEase DNA Isolation Kit. Revision C.0(Revision C.0).

Amini S, Pushkarev D, Christiansen L, Kostem E, Royce T, Turk C, Pignatelli N, Adey A, Kitzman JO, Vijayan K et al. 2014. Haplotype-resolved whole-genome sequencing by contiguity-preserving transposition and combinatorial indexing. Nat Genet 46(12): 1343-1349.

Ciotlos S, Mao Q, Zhang RY, Li Z, Chin R, Gulbahce N, Liu SJ, Drmanac R, Peters BA. 2016. Whole genome sequence analysis of BT-474 using complete Genomics' standard and long fragment read technologies. Gigascience 5: 8.

Drmanac R 2006. Nucleic Acid Analysis by Random Mixtures of Non-Overlapping Fragments. Vol WO 2006/138284.

Drmanac R, Peters, B.A., Alexeev, A. 2013. Multiple tagging of individual DNA fragments. Vol WO 2014/145820 A2 (ed. USPTO).

Duitama J, McEwen GK, Huebsch T, Palczewski S, Schulz S, Verstrepen K, Suk EK, Hoehe MR. 2012. Fosmid-based whole genome haplotyping of a HapMap trio child: evaluation of Single Individual Haplotyping techniques. Nucleic Acids Res 40(5): 2041-2053.

Edge P, Bafna V, Bansal V. 2017. HapCUT2: robust and accurate haplotype assembly for diverse sequencing technologies. Genome Res 27(5): 801-812.

Fan HC, Wang J, Potanina A, Quake SR. 2011. Whole-genome molecular haplotyping of single cells. Nat Biotechnol 29(1): 51-57.

Gulbahce N, Magbanua MJM, Chin R, Agarwal MR, Luo X, Liu J, Hayden DM, Mao Q, Ciotlos S, Li Z et al. 2017. Quantitative Whole Genome Sequencing of Circulating Tumor Cells Enables Personalized Combination Therapy of Metastatic Cancer. Cancer Res 77(16): 4530-4541.

Hellner K, Miranda F, Fotso Chedom D, Herrero-Gonzalez S, Hayden DM, Tearle R, Artibani M, KaramiNejadRanjbar M, Williams R, Gaitskell K et al. 2016. Premalignant SOX2 overexpression in the fallopian tubes of ovarian cancer patients: Discovery and validation studies. EBioMedicine 10: 137-149.

Kitzman JO, Mackenzie AP, Adey A, Hiatt JB, Patwardhan RP, Sudmant PH, Ng SB, Alkan C, Qiu R, Eichler EE et al. 2011. Haplotype-resolved genome sequencing of a Gujarati Indian individual. Nat Biotechnol 29(1): 59-63.

Kuleshov V, Xie D, Chen R, Pushkarev D, Ma Z, Blauwkamp T, Kertesz M, Snyder M. 2014. Whole-genome haplotyping using long reads and statistical methods. Nat Biotechnol 32(3): 261-266.

Li H, Durbin R. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25(14): 1754-1760.

Ma L, Xiao Y, Huang H, Wang Q, Rao W, Feng Y, Zhang K, Song Q. 2010. Direct determination of molecular haplotypes by chromosome microdissection. Nat Methods 7(4): 299-301.

Mao Q, Chin R, Xie W, Deng Y, Zhang W, Xu H, Yu Zhang R, Shi Q, Peters EE, Gulbahce N et al. 2018. Advanced Whole-Genome Sequencing and Analysis of Fetal Genomes from Amniotic Fluid. Clinical chemistry.

Mao Q, Ciotlos S, Zhang RY, Ball MP, Chin R, Carnevali P, Barua N, Nguyen S, Agarwal MR, Clegg T et al. 2016. The whole genome sequences and experimentally phased haplotypes of over 100 personal genomes. Gigascience 5(1): 1-9.

McElwain MA, Zhang RY, Drmanac R, Peters BA. 2017. Long Fragment Read (LFR) Technology: Cost-Effective, High-Quality Genome-Wide Molecular Haplotyping. Methods Mol Biol 1551: 191-205.

McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M et al. 2010. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 20(9): 1297-1303.

Peters BA, Kermani BG, Alferov O, Agarwal MR, McElwain MA, Gulbahce N, Hayden DM, Tang YT, Zhang RY, Tearle R et al. 2015. Detection and phasing of single base de novo mutations in biopsies from human in vitro fertilized embryos by advanced whole-genome sequencing. Genome Res 25(3): 426-434.

Peters BA, Kermani BG, Sparks AB, Alferov O, Hong P, Alexeev A, Jiang Y, Dahl F, Tang YT, Haas J et al. 2012. Accurate whole-genome sequencing and haplotyping from 10 to 20 human cells. Nature 487(7406): 190-195.

Peters BA, Liu J, Drmanac R. 2014. Co-barcoded sequence reads from long DNA fragments: a cost-effective solution for "perfect genome" sequencing. Frontiers in genetics 5: 466.

Picelli S, Bjorklund AK, Reinius B, Sagasser S, Winberg G, Sandberg R. 2014. Tn5 transposase and tagmentation procedures for massively scaled sequencing projects. Genome Res 24(12): 2033-2040.

Schaaf CP, Gonzalez-Garay ML, Xia F, Potocki L, Gripp KW, Zhang B, Peters BA, McElwain MA, Drmanac R, Beaudet AL et al. 2013. Truncating mutations of MAGEL2 cause Prader-Willi phenotypes and autism. Nat Genet 45(11): 1405-1408.

Selvaraj S, J RD, Bansal V, Ren B. 2013. Whole-genome haplotype reconstruction using proximity-ligation and shotgun sequencing. Nat Biotechnol 31(12): 1111-1118.

Suk EK, McEwen GK, Duitama J, Nowick K, Schulz S, Palczewski S, Schreiber S, Holloway DT, McLaughlin S, Peckham H et al. 2011. A comprehensively molecular haplotype-resolved genome of a European individual. Genome Res 21(10): 1672–1685.

Walker RF, Ciotlos S, Mao Q, Chin R, Drmanac S, Barua N, Agarwal MR, Zhang RY, Li Z, Wu MKY et al. 2017. Clinical and genetic analysis of a rare syndrome associated with neoteny. Genetics In Medicine.

Wang O, Chin R, Cheng X, Wu M, Mao Q, Tang J, Sun Y, Lam H, Chen D, Zhou Y et al. 2018. Single tube bead-based DNA co-barcoding for cost effective and accurate sequencing, haplotyping, and assembly. bioRxiv.

Zhang F, Christiansen L, Thomas J, Pokholok D, Jackson R, Morrell N, Zhao Y, Wiley M, Welch E, Jaeger E et al. 2017. Haplotype phasing of whole human genomes using bead-based barcode partitioning in a single tube. Nat Biotechnol 35(9): 852-857.

Zhang K, Zhu J, Shendure J, Porreca GJ, Aach JD, Mitra RD, Church GM. 2006. Long-range polony haplotyping of individual human chromosome molecules. Nat Genet 38(3): 382-387.

Zheng GX, Lau BT, Schnall-Levin M, Jarosz M, Bell JM, Hindson CM, Kyriazopoulou-Panagiotopoulou S, Masquelier DA, Merrill L, Terry JM et al. 2016. Haplotyping germline and cancer genomes with high-throughput linked-read sequencing. Nat Biotechnol.

We would like to acknowledge the ongoing contributions and support of all Complete Genomics and BGI-Shenzhen employees, in particular the many highly skilled individuals that work in the libraries, reagents, and sequencing groups that make it possible to generate high quality whole genome data. This work was supported in part by the Shenzhen Peacock Plan (NO.KQTD20150330171505310). B.A.P. is a recipient of and this work was partially supported by the Research Fund for International Young Scientists, National Natural Science Foundation of China (31550110216).

Employees of BGI and Complete Genomics have stock holdings in BGI.

supplement0.pdf
Supplementary Figure 1 Barcode sequence assembly Three ligations are necessary to generate ~3.6 billion different barcodes. The expected sequence at each step of the barcode assembly is displayed.
supplement0.xlsx
Supplementary Materials Supplementary Materials

PDF

Version 1

posted 11 Oct, 2018

You are reading this latest protocol version

A simple bead-based method for generating cost-effective co-barcoded sequence reads

Status:

Version 1

Abstract

Figures

Introduction

Reagents

Equipment

Procedure

References

Acknowledgements

Competing Interests

Supplementary Files

Associated Publications

Status:

Version 1

Privacy Policy

Terms of Service