Single-tube bead-based DNA co-barcoding for sequencing, haplotype typing, and assembly

Abstract

The present invention describes methods and compositions for preparing a nucleic acid sequencing library, the methods comprising: (a) transposing an insert sequence into a first fragment of a target nucleic acid, wherein the insert sequence comprises a hybridization sequence, and wherein the transposition produces a nick in the first fragment; (b) combining (i) the first fragment of the target nucleic acid from (a), (ii) a splint oligonucleotide, and (iii) a population of beads into a single mixture, wherein each bead comprises a capture oligonucleotide immobilized thereon, and (c) ligating the capture oligonucleotide of a single bead to the inserted hybridization sequence of a single first fragment.

Description

[0001] Citations of Previous Applications

[0002] This application claims the benefit of priority to U.S. Provisional Patent Application No. 62 / 668,757, filed May 8, 2018; U.S. Provisional Patent Application No. 62 / 672,501, filed May 16, 2018; and U.S. Provisional Patent Application No. 62 / 687,159, filed June 19, 2018. The above priority applications are incorporated herein by reference in their entirety for all purposes. Background Art

[0003] To date, the vast majority of individual whole genome sequences lack information about the order of single- to multi-base variants that are delivered as contiguous blocks on homologous chromosomes. A number of technologies have recently been developed to achieve this. Most of these technologies are based on co-barcoding approaches (13), that is, the same barcode is added to subfragments of a single long genomic DNA molecule. After sequencing, the barcode information can be used to determine which reads came from the original long DNA molecule. This approach was originally described by Drmanac (14) and implemented as a 384-well plate assay by Peters et al. (6). However, these methods are technically challenging to implement, expensive, have low data quality, do not provide unique co-barcodes, or some combination of these four. In practice, most of these methods require the generation of a separate whole genome sequence by standard methods to improve variant calling. This has limited the use of these methods, as cost and ease of use are the main factors affecting which technology to use for WGS. Summary of the invention

[0004] Figures and Tables

[0005] Figures 1(A) to 1(D). Overview of stLFR. Figure 1(A): The first step of stLFR involves inserting a hybridization sequence approximately every 200-1000 base pairs into a long genomic DNA molecule. This is accomplished using a transposon. The transposon-integrated DNA is then mixed with beads, each containing The co-barcoded adapter sequence contains a unique barcode shared by all adapters on the bead, a common PCR primer site, and a common capture sequence complementary to that on the integrated transposon. After the genomic DNA is captured to the beads, the transposon is ligated to the barcoded adapter. There are a few additional library processing steps, and then the co-barcoded subfragments are sequenced on a BGISEQ-500 or equivalent sequencer. Figure 1 (B): Mapping of read data by barcodes results in read clustering analysis within a 10 to 350 kb region of the genome. The total coverage and barcode coverage of the 4 barcodes of the 1 ng stLFR-1 library within a small region on Chr11 are shown. Most barcodes are associated with only one read cluster in the genome. Figure 1 (C): The number of original long DNA fragments per barcode is plotted for the 1 ng libraries stLFR-1 and stLFR-2 (orange) and the 10 ng stLFR libraries stLFR-3 and stLFR-4. More than 80% of the fragments in the 1ng stLFR library were co-barcoded by a single unique barcode. Figure 1(D): The fraction of non-overlapping sequence reads covering each original long DNA fragment and captured subfragments (orange) are plotted for the 1ng stLFR-1 library. See also Fig.14 .

[0006] Figures 2(A) to 2(D). SV Detection. The previously reported deletion in NA12878 was also discovered using stLFR data. A heat map of barcode sharing for each deletion can be seen in Figure 10. Figure 2(A) plots the barcode sharing of each deletion on chromosome 8 using the previously described Jaccard index (12). Heatmap of shared barcodes within a 2 kb window of heterozygous deletion regions. Regions of high overlap are indicated in dark red. Those without overlap are indicated in beige. Arrows indicate how regions that are spatially distant from each other on chromosome 8 increase overlap, marking the location of the deletion. Figure 2(B) Co-barcoded reads are separated by haplotype and plotted by unique barcode on the y-axis and by chromosome 8 position on the x-axis. Heterozygous deletions are found in a single haplotype. Figure 2(C) also plots heatmaps for overlapping barcodes between chromosomes 5 and 12 for a patient cell line with a known translocation (26); and Figure 2(D) plots a heatmap for GM20759, a cell line with a known transversion on chromosome 2 (27).

[0007] Figure 3(A) and 3(B): Coverage distribution plots. The coverage of stLFR-2 (A) and standard (B) libraries sequenced on BGISEQ500 are plotted. The coverage of both samples was reduced to 30X. The Poisson distribution of the 30X genome is plotted in blue.

[0008] Figures 4(A) and 4(B): FP overlap between libraries. (A) The FPs of each stLFR library, the BGISEQ-500 standard library, and the PCR-free library sequenced by Illumina (library "HiSeq2500-TruSeq_PCR-Free_DNA_2x251_NA12878" downloaded from basespace) were plotted in a Venn diagram. 2078 FPs were shared between the four stLFR libraries. (B) The overlap of the stLFR library FPs and the Chromium library FPs showed that 1194 FPs were shared between the two different technologies, both of which used DNA isolated from GM12878 instead of the GIAB reference material for NA12878. 884 FPs were unique to the stLFR libraries.

[0009] Figure 5(A) to 5(D): stLFR-1 variant metrics. Read depth and barcode depth were analyzed for reference and variant alleles for all true positive variants, false positive variants, and shared false positive variants (green). Read depth for reference alleles (A) and alternative alleles (B), and barcode counts for reference alleles (C) and alternative alleles (D) are plotted. In general, shared false positives look more like true positives, indicating that there is some filtering criteria that distinguishes these variants from unshared false positives.

[0010] Figure 6(A) to 6(D): stLFR-3 variant metrics. Read depth and barcode depth were analyzed for reference and variant alleles for all true positive variants, false positive variants, and shared false positive variants (green). Read depth for reference alleles (A) and alternative alleles (B), and barcode counts for reference alleles (C) and alternative alleles (D) are plotted. In general, shared false positives look more like true positives, suggesting that there is some filtering criteria that distinguishes these variants from unshared false positives.

[0011] Figures 7(A) and 7(B): Distribution of shared false positive variants. The genomic distances separating the 2078 shared FP variants were summed in consecutive bins of 100 bp (dark blue), 1000 bp (orange), 10000 bp, 100000 bp, and 1000000 bp. Five sets of 2078 randomly selected variants from the stLFR-1 library are also plotted. For each sample, the total number of positions or the total number of variants is plotted. Only variants in a bin or bins where 2 or more variants were found were summed. (A) Before filtering, 219 shared FPs appear to be tightly clustered and are most likely the result of mapping errors. The remaining 1859 variants appear to have a distribution similar to that of a random set of variants. (B) After filtering, 1738 shared FPs remain, but only 72 are tightly clustered.

[0012] Figures 8(A) to 8(T): NA12878 deletion detection using barcode sharing heatmaps. Deletions in the stLFR-1 library were detected using the following data at the following positions: chr3:65189000-65213999 using 230Gb (A) or 100Gb (B) of read data, chr4:116167000-116176999 using 230Gb (C) or 100Gb (D) of read data, chr4:187094000-187097999 using 230Gb (E) or 100Gb (F) of read data, chr7:110182000-110187999 using 230Gb (G) or 100Gb (H) of read data, chr7:110182000-110187999 using 230Gb (I) or 100Gb (J) of read data At chr16:62545000-62549999, 230Gb(K) or 100Gb(L) using read data at chr1:189704509-189783359, 230Gb(M) or 100Gb(N) using read data at chr3:162512134-162569235, 230Gb(O) or 100Gb(P) using read data at chr5:104432113-104467893, at chr6:78967194-79001807, and 230Gb(S) or 100Gb(T) using read data at chr8:39232074-39309652.

[0013] Figures 9(A) to 9(L): Translocation and inversion detection using stLFR. Patient cell lines and cell line GM20759, which carry a translocation between chromosomes 5 and chromosome 12 and an inversion on chromosome 2, respectively, were analyzed with stLFR. For each library, the total sequence coverage was downsampled to investigate the detection capability at lower coverage. The translocation between chromosomes 12 and chromosome 5 was easily detected at total sequence coverage of 40Gb (A), 20Gb (B), 10Gb (C), and even 5Gb (D). The inversion of GM20759 was also easily detected at total sequence coverage of 46Gb (E), 20Gb (F), 10Gb (G), and 5Gb (H). In addition, we investigated these regions in the GM12878 cell line, which is not known to contain any of these SVs. No translocations between chromosomes 5 and 12 were evident from either 1 ng of the stLFR library with 230 Gb coverage (I) or 10 ng of the stLFR library with 126 Gb coverage (J). No transversions were found in the stLFR-1 library (K) or the stLFR-4 library (L).

[0014] Figures 10(A) to 10(C): Alignment dot plots of NA12878 scaffolds. SALSA scaffolds from the stLFR-1 (A) library and the stLFR-4 (B) library were mapped relative to the hg37 reference human genome. 734 million HiC reads obtained from Dixon et al. (29) were also used to generate scaffolds and mapped against hg37 (C). In all cases, only scaffolds covering 5% or more of the chromosome were mapped.

[0015] Fig.11 : LongHap phasing. A full description of the phasing algorithm applied with LongHap can be found in the "Methods and Materials" section.

[0016] Fig.12 : Barcode sequence assembly. Three ligations were necessary to generate approximately 3.6 billion different barcodes. The expected sequences for each step of the barcode assembly are shown as SEQ ID NOs: 1-13.

[0017] Fig.13 : Barcode sequence assembly flow chart.

[0018] Fig.14 : An exemplary flow chart of a barcoding experiment plan.

[0019] Fig.15 : Animation of the hybridization steps.

[0020] Fig.16 : Animation of the ligation and degradation steps. The final steps showing denaturation and C-tailing are optional and are not described here.

[0021] Fig.17 : Production and use of double-chain DNBs.

[0022] Fig.18 : Amplification of long molecules.

[0023] Fig.19 : Random nickase method

[0024] Fig. 20 : The hairpin adapter method.

[0025] Figures 21(A) and 21(B): Figure 21(A) is a schematic diagram of ligation assays on different DNA substrates. The blunt-ended DNA donor is a synthetic partial dsDNA molecule with a dideoxy 3'-end (solid circle) to prevent the adapter from self-ligation. The long arm of the adapter has been 5'-phosphorylated. The DNA acceptor is assembled using 2 or 3 oligonucleotides (black, red and orange lines) to form a nick (no phosphate), a gap (1 or 8nt) or a 36nt 5'-overhang. All strands of the substrate are not phosphorylated and the scaffold strand is protected by 3' dideoxy. Figure 21(B) shows the size migration of the ligation products analyzed using a 6% denaturing polyacrylamide gel. Negative no-ligase controls (lanes 1, 3, 4, 6, 7, 9, 10, 12 and 13) were loaded at 1 volume or half the volume of their corresponding experimental test. If ligation occurs, the substrate size will shift up by 22nt. Red arrows correspond to substrates, blue arrows correspond to adaptor-ligated substrates. M2 = 25 bp DNA ladder from Thermo Fisher (c) Table of expected ligation product sizes, and estimation of ligation efficiency using ImageJ. Ligation efficiency was estimated by dividing the intensity of the ligated product by the total intensity of the ligated and unligated products.

[0026] Figures 22(A) to 22(D): Gel analysis of size migration of ligation products using 6% TBE polyacrylamide gel. Red arrows correspond to substrates, blue arrows correspond to substrates for adaptor ligation: nick (A, left), 5'-overhang (A, right), 1 nt gap (B), 2 nt gap (C) and 3 nt gap (D). M2 = 25 bp DNA ladder from Thermo Fisher. The two adaptor (Adl and Ad2) sequences were compared and the different bases (A or G) at the 5' end of the ligation junction of Ad2 were also examined. **(e) Ligation efficiency table calculated from band intensities using ImageJ.

[0027] Figures 23(A) and 23(B): Figure 23(A) shows a schematic diagram of 3' branch ligation on a DNA / RNA hybrid with a 20 bp complementary region. We tested whether blunt-ended adapters would ligate to the 3'-end of DNA at the 5'-RNA overhang and / or to the 3'-end of RNA at the 5'-DNA overhang. (B) Gel analysis of size migration of ligation products using 6% denaturing polyacrylamide gel. Red arrows correspond to RNA substrates (29nt) and green arrows correspond to DNA substrates (80nt). Blue arrows correspond to adapter-ligated RNA substrates. If ligation occurs, the substrate size will shift up by 20nt. Reactions 1 and 2 were repeated. M2 = 25 bp DNA ladder from Thermo Fisher.

[0028] Figures 24(A) to 24(C): Figure 24(A) is a schematic diagram of transposon insertion followed by 3'-3' branch ligation and PCR amplification using Pr-A (blue arrow) and Pr-B (green arrow). (B) Amplification products after transposon insertion using TnA and / or TnB and / or 3' branch ligation using AdB with primers pr-A, pr-B or both. Products were electrophoresed on a 6% polyacrylamide gel. M1 = ThermoFisher MassRuler low range DNA ladder. (C) Graph of amplification signals using pr-A and pr-B after various transposon insertion and 3' branch ligation conditions.

[0029] Fig.25 : Mid-length mark.

[0030] Figures 26(A) and 26(B): 3' branched ligation of unconventional DNA ends formed at nicks, gaps and overhangs by T4 DNA ligase. (A) Schematic diagram of ligation assays on different DNA acceptor types. The blunt-ended DNA donor is a synthetic partial dsDNA molecule with a dideoxy 3'-end (solid circle) to prevent self-ligation of the DNA donor. The long arm of the donor is 5' phosphorylated. DNA acceptors are assembled using 2 or 3 oligonucleotides to form nicks (no phosphate), gaps (1 or 8 nt) or 36-nt 3'-recessive ends. All strands of the substrate are unphosphorylated and the scaffold strand is protected by 3' dideoxy. (B) Size migration of ligation products of substrates 1, 2, 3 and 4 was analyzed using 6% denaturing polyacrylamide gels. Negative no-ligase controls (lanes 1, 3, 4, 6, 7, 9, 10, 12 and 13) were loaded at 1 volume or half the volume of the corresponding experimental test. If ligation occurs, the substrate size will shift up by 22 nt. The red arrow corresponds to the substrate and the purple arrow corresponds to the donor-ligated substrate. A 25 bp DNA ladder from Thermo Fisher was used. Donor and substrate sequences in Table S1. Table 8 shows the expected sizes of substrates and ligation products in each experimental group and the approximate ligation efficiency. The intensity of each band was estimated using ImageJ and normalized by its expected size. The ligation efficiency was estimated by dividing the normalized intensity of the ligated product by the normalized total intensity of the ligated and unligated products.

[0031] Figures 27(A) to 27(E): Gel analysis of size migration of ligation products using 6% TBE polyacrylamide gel. Red arrows correspond to substrates and purple arrows correspond to substrates for donor ligation: substrate 5 (nicked) (A), substrate 6 (1-nt gap) (B), substrate 7 (2-nt gap) (C), substrate 8 (3-nt gap) (D), and substrate 9 (3'-recessed end) (E). Thermo Fisher's 25 bp DNA ladder was used. Three DNA donors with different bases at the 5' end (T, A or GA) of the ligation junction were examined. Table 9 shows the ligation efficiency calculated based on normalized band intensity using ImageJ.

[0032] Figures 28(A) to 28(D). 3' branched ligation of RNA 3' end in DNA / RNA hybrids. Schematic diagram of 3' branched ligation on DNA / RNA hybrids with 20 bp complementary regions. We tested whether a blunt-ended DNA donor would ligate to the 3' recessed end of DNA and / or the 3' recessed end of RNA. DNA (ON-21) hybridized to the RNA strand (A), while DNA (ON-23) failed to hybridize to the RNA strand (B). Figures 28(C) and (D) show gel analysis of the size migration of the ligation products using a 6% denaturing polyacrylamide gel. The red arrow corresponds to the RNA substrate (29 nt) and the green arrow corresponds to the DNA substrate (80 nt). The purple arrow corresponds to the donor-ligated RNA substrate. If ligation occurs, the substrate size will shift up by 20 nt. (c) Lanes 1 and 2, replicate experiments; lanes 7-10, no ligase control; 10% PEG with T4 DNA ligase was added. (d) Lane 1, no ligase control; Lanes 2, 3 and 8, T4 DNA ligase with 10% PEG; Lanes 4, 5 and 9, T4 RNA ligase 1 with 20% DMSO; Lanes 6, 7 and 10, T4 RNA ligase 2 with 20% DMSO. Thermo Fisher 25 bp DNA ladder was used. It corresponds to Figure 23, but is not completely accurate in many aspects.

[0033] Figures 29(A) to 29(C) show schematics of three transposon tagging methods followed by PCR amplification using Pr-A (blue arrows) and Pr-B (green arrows). Two transposon methods (A); 1 Y transposon tag with 3' gap filling (B); one transposon method with adapter ligation at the 3' gap (C). Figure 29(D) is a graph of amplification signals after purification using pr-A or pr-A and pr-B after various tagging and gap ligation conditions. May correspond to Figure 23, but is not completely accurate in several respects.

[0034] Figures 30(A) to 30(C). Deviation in base distribution of Tn5 gap junction (A), two transposons (B), and conventional TA junction (C). Only the first 20 bases at both ends of the junction are shown; adenine, blue; cytosine, orange; guanine, grey; thymine, yellow; averages and standard deviations of five independent libraries are given. Not present at all.

[0035] Figures 31(A) and 31(B). DNA 3' branch ligation under different addition conditions. (A) Ligation was performed with titrated concentrations of ATP on 5'-overhanging DNA. The experiment was repeated with 0.01 mM ATP (lanes 4 and 5) and 0.005 mM ATP (lanes 6 and 7). Lane 9 is a control without donor. (B) 3' branch ligation of DNA at nick, 1-nt gap, 8-nt gap, 5'-overhang and blunt end with or without SSB and ligase. Red arrows correspond to substrate and purple arrows correspond to donor-ligated substrate. Not present at all.

[0036] Table 1: Phasing and variant calling statistics. Unless otherwise stated, reads were mapped to Hg37 with bait sequences and variants were called using GATK with default settings for all libraries. SNPs from the GIAB high confidence variant calling VCF were used as input for phasing.

[0037] Table 2: Scaffold statistics.

[0038] Table 3: Filtering reduces false positive calls. Final FP calls were calculated by subtracting 1666 from the filtered FPs, except for the STD library, which by definition does not share any of these FPs with the stLFR library, as the stLFR library was made from the GIAB reference.

[0039] Table 4: LongHap SNP and Indel phasing.

[0040] Table 5: Filtering criteria. Various filtering criteria described in the Materials and Methods section were used to remove FPs.

[0041] Table 6: Exemplary sequences. DETAILED DESCRIPTION

[0042] 1.stLFR library process

[0043] 1.1 Introduction

[0044] Here, we describe the implementation of single-tube long fragment read (stLFR) technology (15), an efficient DNA co-barcoding technology that can achieve millions of barcodes in a single tube. See WO 2014 / 145820A2 (2014), incorporated herein by reference for all purposes. This is achieved by using the surface of microbeads to replace compartments (e.g., the wells of a 384-well plate). Each bead carries many copies of a unique barcode sequence, which is transferred to subfragments of each long DNA molecule. These co-barcoded subfragments are then analyzed on a common short-read sequencing instrument (e.g., BGISEQ-500 or equivalent). In this implementation of this method, we used a ligation-based combinatorial barcode generation strategy to create more than 1.8 billion different barcodes through three ligation steps. For a single sample, we use beads with approximately 10-50 million barcodes to capture approximately 10-100 million long DNA molecules in a single tube. Few two beads will share the same barcode because we sampled 10-50 million beads from such a large total barcode library. In addition, with 50 million beads and 10 million long genomic DNA fragments, the vast majority of subfragments of each long DNA fragment were co-barcoded by a unique barcode. This is similar to single-molecule sequencing of long reads and may provide a powerful informatics approach for de novo assembly. Importantly, stLFR is easy to perform and can be achieved with relatively small oligonucleotide inputs to generate barcoded beads. In addition, stLFR uses standard equipment that can be found in almost all molecular biology laboratories and can be analyzed by almost any sequencing strategy. Finally, stLFR replaces standard NGS library preparation methods, requires only 1 ng of DNA, and does not significantly increase the cost of whole genome or whole exome analysis, with a total cost of less than $30 per sample.

[0045] As used herein, "single tube" refers to the analysis of a large number of individual DNA fragments without the need to separate the fragments into separate tubes, containers, aliquots, wells or droplets during the labeling step. Instead, the surface of the beads can replace the compartments.

[0046] The first step of stLFR is to insert hybridization sequences, preferably at regular intervals, along the genomic DNA fragment. Appropriate spacing may vary depending on the application and the desired outcome, but is typically between 100-1500 bp, typically 200-1000 bp. This is accomplished by transposition of the bound DNA sequence. In one embodiment, the transposase is Tn3, Tn5, Tn7, or Mu. Typically, Tn5 transposase is used (see Picelli et al., 2014, incorporated herein by reference for all purposes). The transposed DNA or insert sequence contains a single-stranded region for hybridization ("hybridization sequence") and a double-stranded mosaic sequence (mosaic sequence) that is recognized by the enzyme and capable of undergoing the transposition reaction ( Figure 1A ). This transposition step is accomplished in solution (as opposed to directly attaching the insert sequence to the beads). This allows very efficient incorporation of the hybridizing sequence along the genomic DNA molecule. As previously observed (10), transposases have the property of remaining bound to genomic DNA after the transposition event, effectively leaving the long genomic DNA molecule with the transposon integrated intact.

[0047] After the DNA is treated with, for example, Tn5, it is diluted in hybridization buffer and combined with the cloned barcoded beads. In one approach (Examples below), Clonal barcoded beads. Each bead contains approximately 400,000 capture adapters (also called capture oligonucleotides or capture oligonucleotides), each containing the same barcode sequence. A portion of the capture adapters contain uracil nucleotides to enable destruction of unused adapters in subsequent steps. For example, the capture adapters can be 5-50% uracil, more often 5-50%, and more often 5-20%. The mixture is incubated under optimized temperature and buffer conditions, during which the transposon-inserted DNA is captured onto the beads via the hybridization sequence.

[0048] It has been proposed that genomic DNA in solution forms a ball with two tails sticking out (16). This can enable capture of long DNA fragments toward one end of the molecule, followed by a rolling motion that wraps the genomic DNA molecule around the bead. Each bead has a capture oligonucleotide approximately every 7.8 nm on its surface. This enables very uniform and high rate capture of subfragments. A 100 kb genomic fragment will wrap around a 2.8 um bead approximately 3 times. In our data, the longest fragment size captured was 300 kb, suggesting that larger beads may be required to capture longer DNA molecules.

[0049] In alternative embodiments, parameters such as bead size, oligonucleotide spacing, or the number of different oligonucleotides in each mixture can be varied. For example, the diameter of the beads used can be 1-20 μm, or 2-8 μm, 3-6 μm, or 1-3 μm. For example, the spacing of the barcoded oligonucleotides on the beads can be at least 1 nm, at least 2 nm, at least 3 nm, at least 4 nm, at least 5 nm, at least 6 nm, or at least 7 nm. In some embodiments, the spacing is less than 10 nm (e.g., 5-10 nm), less than 15 nm, less than 20 nm, less than 30 nm, less than 40 nm, or less than 50 nm. In some embodiments, the number of different barcodes used for each mixture can be>1M,>10M,>30M,>100M,>300M, or>1B. As described below, for example, using the methods described herein, a large number of barcodes can be generated for use in the present invention. In some embodiments, the number of different barcodes used for each mixture can be >1M, >10M, >30M, >100M, >300M or >1B, and they are sampled from a pool with at least 10-fold more diversity (e.g., >10M, >0.1B, 0.3B, >0.5B, >1B, >3B, >10B different barcodes from beads).

[0050] Individual barcode sequences are transferred at regular intervals by ligating the 3' end of the capture adaptor to the 5' end of the hybridizing sequence of the transposon insert, mediated by a bridge or splint oligonucleotide (the terms are used interchangeably) having a first region complementary to the capture adaptor and a second region complementary to the hybridizing sequence ( Figure 1A and Fig.15 ). The beads are collected and the DNA / transposase complexes are disrupted, generating subfragments less than 1 kb in size.

[0051] If necessary, sample barcodes can be obtained in this step. Transposons carrying unique barcodes between chimeric sequences and hybridization sequences are used. These can be synthesized in 96, 384 or 1536 plate formats, with each well containing many copies of transposons carrying the same barcode, and each barcode is different between the wells. Using these barcode transposons, different DNA samples can be transposon-inserted in 96, 384 or 1536 plate formats. Samples labeled with sample barcodes can be multiplexed in any manner.

[0052] Due to the large number of beads and the high capture oligonucleotide density per bead, the amount of excess adaptor is four orders of magnitude greater than the amount of product. This large amount of unused adaptor can overwhelm the following steps. To circumvent this, we designed beads with capture oligonucleotides ligated at the 5' end. This enabled the development of an exonuclease strategy that specifically degrades excess unused capture oligonucleotides. See Fig.14 and 16 Uracil DNA glycosylase (UDG) can also be used to degrade excess adaptors.

[0053] In one aspect, the method comprises combining (i) a first fragment of a target nucleic acid and (ii) a population of beads into a single mixture, wherein each bead comprises an oligonucleotide immobilized thereon, the oligonucleotide comprising a tag-containing sequence (or a barcode adapter), wherein each tag-containing sequence comprises a tag sequence, wherein the oligonucleotides immobilized on the same single bead comprise the same tag-containing sequence, and a majority of the beads have different tag sequences. In some embodiments, the DNA fragment is a concatamer of at least 2, at least 10, at least 30, or at least 100 copies of a DNA or cDNA molecule. The length of the nucleic acid monomer can be 0.5 kb to 10 kb, or >1 kb, or >10 kb. In some methods, the sequence of the bases of >50% or >70%>90%, 95%, >99%, 100% of the DNA or cDNA molecules in the mixture is determined.

[0054] 1.1.1 Two transposon methods

[0055] In one approach to obtaining stLFR, two different transposons are used in the initial insertion step, allowing PCR to be performed after exonuclease treatment. However, this approach results in coverage of approximately 50% or less per long DNA molecule because it requires the insertion of two different transposons adjacent to each other to generate a suitable PCR product.

[0056] 1.1.2 Single transposon method using 3' branch ligation

[0057] To achieve the highest coverage per genomic DNA fragment, we used a single transposon in the initial insertion step and added additional adapters by ligation. This noncanonical ligation, called 3' branched ligation, involves covalently attaching the 5' phosphate from the blunt-ended adapter to the recessed 3' hydroxyl group of the genomic DNA ( Figure 1A). Branch connection is described in Example 3 below. See also U.S. Patent Publication No. US2018 / 0044668 and International Application No. WO2016 / 037418, which are incorporated herein by reference for all purposes. See also U.S. Patent Publication No. 2018 / 0044667, which is incorporated by reference for all purposes. Using this method, it is theoretically possible to amplify and sequence all subfragments of the captured genomic molecules.

[0058] Additionally, this ligation step enables sample barcodes to be placed adjacent to genomic sequences for multiplexed sampling. The benefit of using these adapters for sample barcoding is that the barcodes can be placed adjacent to the genomic DNA, allowing the barcode and genomic DNA to be sequenced using the same primers, without the need for additional sequencing primers to read the barcodes. Sample barcoding allows preparations from multiple samples to be pooled prior to sequencing and differentiated by barcode. 3' branched ligation adapters can be synthesized in a 96, 384, or 1536 plate format, with each well containing many copies of the adapter carrying the same barcode, and each barcode being different between wells. After capture on the beads, these adapters can be used for ligation in a 96, 384, or 1536 plate format.

[0059] After this ligation step, PCR is performed and the library is ready to enter any standard next generation sequencing (NGS) workflow. It should be understood that a first primer that hybridizes to a site on the capture oligonucleotide or its complement (see Figure 1A ) and a second primer that hybridizes to a site of the 3' branched ligation adapter or its complement. In the case of the BGISEQ-500, the library is circularized as previously described (17). DNA nanoballs are prepared from the single-stranded circles and loaded onto patterned nanoarrays (17). These nanoarrays are then sequenced on the BGISEQ-500 based on combined probe-anchor synthesis (cPAS) (18-20). After sequencing, the barcode sequences are extracted. Mapping the read data by unique barcodes shows that most reads with the same barcode are clustered in the genomic region that corresponds to the DNA length used during library preparation ( Figure 1B ). This method and the experimental protocol for preparing the beads are described in detail in Examples 1 and 2.

[0060] In some embodiments, >50%, >70%, >80%, >90% or >95% of the barcoded DNA fragments are barcoded with unique barcodes. In some embodiments, >50%, >70%, >80%, >90% of the subfragments in the fragments are linked to barcode oligonucleotides. In some embodiments, an average of >10% or >20%, >40%, >50%, >60% of the subfragments of the long fragments are sequenced.

[0061] 1.2stLFR read coverage and variant calling

[0062] To demonstrate phasing and variant calling of stLFR, we generated four libraries using 1ng (stLFR-1 and stLFR-2) and 10ngs (stLFR-3 and stLFR-4) of DNA from NA12878. The number of beads used varied from 10 million (stLFR-3), 30 million (stLFR-4), and 50 million (stLFR-1 and stLFR-2). Finally, the 3' branch joining method (stLFR-1, stLFR-2, and stLFR-3) and two transposon methods (stLFR-4) were tested. The sequencing depth of stLFR-1 and stLFR-2 reached 336Gb and 660Gb of total base coverage, respectively. We also analyzed the coverage of downsampling. stLFR-3 and stLFR-4 sequenced to a moderate level of 117Gb and 126Gb, respectively. The co-barcoded reads were mapped to build 37 of the human reference genome using BWA-MEM (21). Since stLFR does not require any preamplification step, the distribution of read coverage across the genome was close to a Poisson distribution (Figure 3). Non-replicated coverage was 34–58X, and the number of long DNA molecules per barcode was 1.2–6.8 (Tables 1 and Figure 1C As expected, the stLFR library made from 50 million beads and 1 ng of genomic DNA had the highest single unique barcode co-barcoding rate of over 80% ( Figure 1C These libraries also observed the highest average non-overlapping read coverage of 10.7–12.1% per long DNA molecule and the highest average non-overlapping base coverage of 17.9–18.4% per captured subfragment of each long DNA molecule (Fig. 1d). This coverage is approximately 10-fold higher than that previously demonstrated using 3 ng of DNA and transposons attached to beads (12).

[0063] For each library, variants were called using GATK (22) with default settings. SNP and indel calls were compared to the genome in a bottle (GIAB) (23), and false positive rates (FP) and false negative rates (FN) were determined (Table 1). In addition, we used the same settings in GATK on a standard non-stLFR library made from approximately 1000 times more genomic DNA and also sequenced on the BGISEQ-500 (STD) and on a Chromium library from 10X genome (11). We also compared the accuracy and sensitivity reported by Zhang et al. in a bead haplotype library study (12), which is incorporated herein by reference for all purposes. Compared to the Chromium library, our stLFR method and the method described by Zhang et al. showed lower SNP and Indel FP rates. The FP and FN rates of stLFR were 2 times higher than those of the STD library, and the FN rates were higher or lower than those of the Chromium library depending on the specific stLFR library and filtering criteria. The higher FN rate in the stLFR library compared to the standard library is mainly due to the shorter average insert size ( 300 bp in the standard library). That is, stLFR has much lower FN rates for SNPs and indels than Zhang et al., and much lower FN rates for indels than the Chromium library (Table 1). Overall, our stLFR library performs better than Zhang et al. or the Chromium library for most metrics of variant calling, especially when non-optimized mapping and variant calling methods are used (Table 1, “No filter”).

[0064] One potential problem with using GIAB data to measure FP rates is that we were unable to use the GIAB reference material (NIST RM 8398) due to the small size of the isolated DNA fragments. We therefore used the GM12878 cell line and isolated the DNA using a dialysis-based method that is capable of producing very high molecular weight DNA (see Methods). However, our GM12878 cell line isolates may have many unique somatic mutations compared to the GIAB reference material, resulting in an increased number of FPs in our stLFR libraries. To examine this further, we compared the overlap of single-nucleotide FP variants between the four stLFR libraries and the two non-LFR libraries (Fig. 4a). Overall, 544 FP variants were shared between the six libraries, while 2078 FP variants were unique to the four stLFR libraries. We also compared the stLFR FPs with the Chromium libraries and found that more than half of these shared FPs (1194) were also present in the Chromium libraries (Fig. 4b). Inspection of the read and barcode coverage of these shared variants showed that they were more similar to TP variants (Figs. 5-6). We also examined the distribution of these shared FP variants with 2078 randomly selected variants in the genome (Figure 7a). This analysis showed that 219 variants were found in clusters where two or more of these FPs were within 100bp of each other. However, the distribution of most (90%) variants did not appear to be different from randomly selected variants. In addition, of those FPs shared between stLFR and Chromium libraries, only 41 were found to be clustered (Figure 7a). Ultimately, GIAB called 96 of these variants, but with different zygosity than those called in the stLFR library.

[0065] If we accept the evidence that these shared FP variations are largely real and not present in the GIAB reference material, the FP rate of stLFR may be up to 1859 fewer variations than reported for SNP detection in Table 1. This is still thousands of single nucleotide variations more than the standard BGISEQ-500 library. To further improve the FP rate in the stLFR library, we tested a variety of different filtering strategies to eliminate errors. Ultimately, by applying some filtering criteria based on reference and variant allele ratios and barcode counts (see Examples), we were able to remove 3647-13840 FP variations depending on the library and coverage. Importantly, this was achieved while only increasing the FN rate by 0.10-0.29% in the stLFR library. After this filtering step, we examined the shared FPs between the four stLFR libraries. Filtering removed only 340 shared FP variations, of which 147 were clustered within 100 base pairs of each other and were likely not real (Figure 7b). This further suggests that most of these shared FPs are real variations. Taking these variations into account and the reduction in the number of FP variants after filtering, the FP rate was similar and the FN rate was 2-3 times higher compared to the filtered STD library using SNP calling (Table 3). This increased FN rate was mainly due to the increase in non-unique mapping of partner pairs with short insert lengths in the stLFR library.

[0066] 1.3stLFR phasing performance

[0067] To evaluate the phasing performance of the variants, high-confidence variants from GIAB were phased using the publicly available software package HapCut2 (24). Depending on the library type and the amount of sequence data, more than 99% of all heterozygous SNPs were placed in contigs with N50s ranging from 0.6 to 15.1 Mb (Table 1). The stLFR-1 library with 336 Gb total read coverage (44X unique genome coverage) achieved the highest phasing performance with an N50 of 15.1 Mb. The N50 length appears to be primarily affected by the length and coverage of long genomic fragments. This can be seen from the reduced N50 of stLFR-2, because the DNA used for this sample was slightly older and more fragmented than the material used for stLFR-1 (Table 1, average fragment length 52.5 kb versus 62.2 kb), while the N50 of the 10 ng library (stLFR-3 and 4) was 10 times shorter. Comparison with the GIAB data showed that short and long switch error rates were low and comparable to previous studies (11, 12, 25). stLFR performance was very similar to the Chromium libraries. Since no read data were available for the bead haplotype method by Zhang et al., we could only compare our results with those of a phasing algorithm written and optimized specifically for their data. This showed that the stLFR-1 and stLFR-2 libraries had longer N50s, similar short switch error rates, but higher long switch error rates. stLFR-3 and stLFR-4, which use more DNA, had N50s similar to those of Zhang et al. However, direct comparison is difficult due to differences in DNA input and coverage.

[0068] It should be noted that this phasing result was obtained using a program that was not written for stLFR data. To see if this result could be improved, we developed a phasing program, LongHap, and optimized it specifically for stLFR data. Using GIAB variants, LongHap was able to phase 99% of the SNPs into contigs with an N50 of 18.1Mb (Table 1). Importantly, this increase in contig length was achieved while reducing short and long switching errors (Table 1). LongHap can also phase indels. Applying LongHap to stLFR-1 using GIAB SNPs and indels produced a 23.4Mb N50, but also resulted in an increased switching error rate (Table 4).

[0069] 1.4 Structural variation detection

[0070] Previous studies have shown that long-range information can improve the detection of structural variants (SVs) and have described large deletions (4–155 kb) in NA12878 (11, 12). To demonstrate the ability of stLFR to detect SVs, we examined barcode overlap data from stLFR-1 and stLFR-4 libraries in these regions as previously described (12). In each case, the deletion was observed in the stLFR-1 data, even at lower coverage (Figures 2a and 8). Close inspection of coverage of chromosome 8 The co-barcode sequence reads for the deletion indicated that the deletion was heterozygous and present in a single haplotype (Figure 2b-c). The 10 ng stLFR-4 library also detected most deletions, but the smallest three were difficult to identify due to the lower coverage per fragment of this library (and therefore less barcode overlap).

[0071] To evaluate the performance of stLFR for detecting other types of SVs, we prepared libraries for cell lines from patients with known translocations between chromosomes 5 and 12 (26) and for the cell line GM20759 with a known inversion in chromosome 2 (27). The stLFR libraries were able to identify inversions and translocations in the corresponding cell lines (Fig. 2d). Downsampling the read amount for each library showed that strong signals for translocations were detected even with only 5 Gb of read data ( Finally, examination of the two SVs in the stLFR-1 library did not result in a clear pattern (Fig. 9i-l), indicating that the false positive rate for detecting these types of SVs is low.

[0072] 1.5 Scaffold contigs using stLFR

[0073] stLFR is a powerful method in part because it uses a very large number (e.g., approximately 1.8 billion) of unique barcodes and enables co-barcoding that is unique to each individual long genomic DNA molecule. Such data should be useful for de novo genome assembly and improved scaffolding. To demonstrate how stLFR can be used to improve genome assembly, we used reads from the stLFR-1 and stLFR-4 libraries and SALSA (28), a program designed for chromatin conformation capture (Hi-C) data to support single-molecule, real-time (SMRT) read assembly of NA12878 (29). SALSA was not designed for stLFR data, so it was necessary to modify the stLFR data to resemble a Hi-C structure. This was achieved by selecting read pairs that share the same barcode and are located at the ends of the captured long DNA molecules. These were then marked as read pairs for the SALSA program. Replacing Hi-C data with stLFR data produced excellent scaffolding. Using only 60 million stLFR reads, 1411 contigs were linked to 597 scaffolds with an N50 of 44.7 Mb. These scaffolds cover 2.84 Gb of the genome. These metrics compare very favorably to those generated in the SALSA manuscript using the same contigs and 10x (734 million) Hi-C reads generated from human embryonic stem cells (Table 2). The quality of the stLFR scaffolds was further analyzed by aligning the stLFR scaffolds to build 37 of the human reference genome and comparing them with the dnadiff program (31). In general, the stLFR scaffolds matched the reference genome very well, and the number of breakpoints, translocations, relocations, and inversions was similar to that of scaffolds generated using Hi-C reads (Table 2). The alignment dot plots further demonstrated the high degree of continuity between the stLFR scaffolds and the reference genome (Figure 10).

[0074] 1.6 Discussion

[0075] Here, we describe an efficient whole-genome sequencing library preparation technology, stLFR, that enables subfragments of long genomic DNA molecules to be co-barcoded with individual unique clonal barcodes in a single-tube approach. Using microbeads as miniaturized compartments, a virtually unlimited number of clonal barcodes can be used per sample at negligible cost. Our optimized hybridization-based capture of transposon-inserted DNA on beads, combined with 3′ branch ligation and exonuclease degradation of an extreme excess of capture adapters, successfully barcoded up to ~20% of subfragments in DNA molecules up to 300 kb in length. Importantly, this was achieved without DNA amplification of the original long DNA fragments and without the attendant representation bias. In this way, stLFR addresses the cost and limited co-barcoding capabilities of emulsion-based approaches.

[0076] The quality of variant calling using stLFR is very high and will likely reach the quality of standard WGS methods with further optimization, but with the added benefit that common barcodes enable advanced informatics applications. We demonstrate high-quality, near-complete genome phasing into long contigs with very low error rates, detection of SVs, and scaffolding of contigs enabling de novo assembly applications. All of this is achieved with a library that requires neither specialized equipment nor the added cost of library preparation.

[0077] Due to effective barcoding, we have successfully used as little as 1ng of human DNA (600X genome coverage) to prepare stLFR libraries and obtained high-quality WGS, in which most subfragments were uniquely co-barcoded. Less DNA can be used, but stLFR does not use DNA amplification during the co-barcoding process, so overlapping subfragments are not generated from each individual long DNA molecule. Therefore, as the amount of DNA decreases, the coverage of the entire genome is compromised. In addition, since stLFR currently retains 10-20% of each original long DNA molecule and then performs PCR amplification, a sampling problem is created. This leads to relatively high read duplication rates and increases sequencing costs, but it may be improved. An obvious solution is to remove the PCR step. This will eliminate sampling, but it can also greatly reduce false positive error rates and false negative error rates. In addition, improvements such as optimizing the insertion distance between transposons and increasing the length of sequencing reads to 200 bases at the paired end should be easy to implement and will increase coverage and overall quality. For some applications, such as structural variation detection, it may be necessary to use less DNA and less coverage. As we demonstrate here, only 5 Gb of sequence coverage can faithfully detect inter- and intra-chromosomal translocations, in cases where the duplication rate is negligible. Indeed, stLFR may represent a simple and cost-effective alternative to long mate pair libraries in the clinical setting.

[0078] Additionally, we believe that this type of data could enable complete diploid de novo assembly of a single stLFR library without the need for long physical reads, such as those produced by SMRT or nanopore technologies. An interesting feature of transposon insertion is that it produces an overlap of 9 bases of sequence between adjacent subfragments. Typically, these adjacent subfragments are captured and sequenced, doubling the synthetic length of the read (e.g., for a 200-base read, two adjacent captured subfragments create two 200-base reads, overlapping by 9 bases, or 391 bases). stLFR does not require special equipment like droplet-based microfluidics methods, and the cost per sample is minimal. In this paper, we demonstrate that 50 million beads can be used, but more beads can be used. This will enable many types of cost-effective analyses, in which hundreds of millions of barcodes will be useful. We envision that this inexpensive large-scale barcode can be used for RNA analysis, such as full-length mRNA sequencing of 1000s cells in combination with single-cell technology or deep population sequencing of 16S RNA in microbial samples. Using stLFR, it is also possible to perform phase-angle chromatin mapping or methylation studies using assay for transposase-accessible chromatin (ATAC-seq) (32).

[0079] 1.7 Target nucleic acid

[0080] As used herein, the term "target nucleic acid" (or polynucleotide) or "nucleic acid of interest" refers to any nucleic acid (or polynucleotide) suitable for processing and sequencing by the methods described herein. Nucleic acids can be single-stranded or double-stranded, and can include DNA, RNA or other known nucleic acids. The target nucleic acid can be the nucleic acid of any organism, including but not limited to viruses, bacteria, yeast, plants, fish, reptiles, amphibians, birds and mammals (including but not limited to mice, rats, dogs, cats, goats, sheep, cattle, horses, pigs, rabbits, monkeys and other non-human primates and humans). The target nucleic acid can be obtained from an individual or multiple individuals (i.e., a population). The sample from which the nucleic acid is obtained can contain nucleic acids from cells or even a mixture of organisms, for example: a human saliva sample including human cells and bacterial cells; a mouse xenograft comprising mouse cells and transplanted human tumor cells, etc. The target nucleic acid may not be amplified, or it may be amplified by any suitable nucleic acid amplification method known in the art. The target nucleic acids may be purified according to methods known in the art to remove cellular and subcellular contaminants (lipids, proteins, carbohydrates, nucleic acids other than the nucleic acid to be sequenced, etc.), or they may not be purified, i.e., include at least some cellular and subcellular contaminants, including but not limited to intact cells, which are disrupted to release their nucleic acids for processing and sequencing. The target nucleic acids may be obtained from any suitable sample using methods known in the art. These samples include but are not limited to: tissues, isolated cells or cell cultures, body fluids (including but not limited to blood, urine, serum, lymph, saliva, anal and vaginal secretions, sweat and semen); air samples, agricultural samples, water samples and soil samples, etc. Non-limiting examples of target nucleic acids include "circulating nucleic acids" (CNA), which are nucleic acids circulating in human blood or other body fluids (including but not limited to lymph, cerebrospinal fluid (liquor), ascites, milk, urine, feces and bronchial lavage fluid) and can be distinguished as cell-free (CF) nucleic acids or cell-associated nucleic acids (reviewed in Pinzani et al., Methods 50:302-307, 2010).

[0081] The target nucleic acid can be genomic DNA (e.g., from a single individual), cDNA, and / or can be a composite nucleic acid, including nucleic acids from multiple individuals or genomes. Examples of composite nucleic acids include microbiome, circulating fetal cells in pregnant women's blood (see, e.g., Kavanagh et al., J. Chromatol. B 878: 1905-1911, 2010), tumor cells (CTCs) circulating in cancer patients' blood (see, e.g., Allard et al., Clin Cancer Res. 10: 6897-6904, 2004). Another example is genomic DNA from a single cell or a small amount of cells, such as from a biopsy (e.g., fetal cells biopsied from blastocyst trophectoderm; cancer cells aspirated from a needle of a solid tumor, etc.). Another example is a pathogen in tissue, blood or other body fluids, such as a bacterial cell, a virus or other pathogen. As used herein, the term "composite nucleic acid" refers to a large number of different nucleic acids or polynucleotides. In certain embodiments, the target nucleic acid is genomic DNA; exome DNA (a subset of complete genomic DNA enriched for transcribed sequences, which includes a set of exons in the genome); transcriptome (i.e., a collection of all mRNA transcripts produced in a cell or cell population, or cDNA produced from such mRNA); methylome (i.e., a population of methylation sites and patterns of methylation in a genome); exome (i.e., protein-coding regions of a genome selected by exon capture or enrichment methods; microbiome; a mixture of genomes of different organisms; a mixture of genomes of different cell types of an organism; and other complex nucleic acid mixtures containing a large number of different nucleic acid molecules (examples include, but are not limited to, microbiome, xenografts, solid tumor biopsies containing normal cells and tumor cells, etc.), including subsets of the above-mentioned types of complex nucleic acids. In one embodiment, such a complex nucleic acid has a complete sequence containing at least one gigabase (Gb) (a diploid human genome contains approximately 6Gb of sequence).

[0082] In some cases, the target nucleic acid or the first fragment is a genomic fragment. In some embodiments, the genomic fragment is longer than 10kb, for example, 10-100kb, 10-500kb, 20-300kb or longer than 100kb. The quantity of the DNA (e.g., human genomic DNA) used in a single mixture can be <10ng, <3ng, <1ng, <0.3nm or <0.1ng of DNA. In some cases, the length of the target nucleic acid or the first fragment is 5000 to 100000KB.

[0083] 1.8 Other methods

[0084] Although the working examples described herein use polymerase chain reaction, other nucleic acid amplification methods may be used. It is within the capabilities of those skilled in the art to make modifications suitable for appropriate amplification techniques.

[0085] Fig.17 -B5 shows other methods. Fig.17 The production of double-stranded DNBs is shown, which can be inserted by transposons and captured by stLFR beads. Up to several thousand copies (e.g., 10-10,000 copies, such as 10-1000 copies or 100-1000 copies) can be replicated on the same DNA strand. This can be sequenced by stLFR to provide high coverage of the original molecule. Fig.18 It is demonstrated that when limited amounts of template DNA are available, a limited pre-amplification step can be used prior to stLFR. Fig.19 A method is described in which a low concentration of random nicking enzyme is used, a moderate concentration of Klenow fragment is used, and a high concentration of ligase is used. The concentrations of beads and DNA are suitable for stLFR. When a nick is formed by Klenow and opens to a vacancy, ligation occurs immediately and locks the long fragment to the bead. The nicking is allowed to proceed, and more vacancies are opened to allow more adapters to be attached in the vacancies. Primer extension produces approximately 500 base pair fragments. A second adapter is ligated to the blunt end, and the library is sequenced. Fig. 20 The ligation of hairpin adapters on long DNA is shown, and the use of primers and Ph29 or similar polymerases in the loop to generate concatenated dsDNA before barcoding. In addition to improving the read coverage of each molecule, an interesting result of this method is that the total "length" (number of bases) of each concatemer is similar at the end of the 0.5-3h polymerase reaction, regardless of the starting fragment length. This provides the option of using barcode beads with a binding capacity corresponding to the size of the concatemer, thereby preventing each bead from binding multiple concatemers. This will reduce the number of beads required for each reaction, thereby further reducing costs.

[0086] Figure B5 shows a method for intermediate length labeling. In one method, 96 or more different barcode transposons are connected in groups of 10 or less by a linker portion (e.g., DNA, a long inert molecule such as dextrin or polyethylene glycol (PEG), or a long protein such as keratin or collagen). Hybridization and connection can be used to connect the transposon to the linker DNA. Other methods can be connected to these molecules by chemical bonds or by avidin, and then biotin is connected to the transposon. This achieves two things, it controls the insertion distance between transposons and provides proximity information (10kb or less) of the intermediate read. This is useful for analyzing repetitive sequences (tandem repeats, trinucleotide mapping, etc.). As with other stLFR methods described herein and elsewhere, DNA containing the inserted sequence can be captured on beads. See Joseph C. Mellor et al., "Phased NGS Library Generation via Tethered Synaptic Complexes," seqWell (2017), available on the World Wide Web at (http: / / )seqwell.com / wp-content / uploads / 2017 / 02 / seqWell_LongBow_poster_AGBT2017.pdf (last visited May 16, 2018).

[0087] 1.9 References for Part 1

[0088] 1. K. Zhang et al., Long-range polony haplotyping of individual humanchromosome molecules. Nat Genet 38, 382-387 (2006).

[0089] 2. L.Ma et al., Direct determination of molecular haplotypes bychromosome microdissection. Nat Methods 7, 299-301 (2010).

[0090] 3.JOKitzman et al., Haplotype-resolved genome sequencing ofaGujarati Indian individual.Nat Biotechnol 29,59-63(2011).

[0091] 4.E.K.Suk et al.,A comprehensively molecular haplotype-resolvedgenome of a European individual.Genome Res 21,1672-1685(2011).

[0092] 5.H.C.Fan,J.Wang,A.Potanina,S.R.Quake,Whole-genome molecularhaplotyping of single cells.Nat Biotechnol 29,51-57(2011).

[0093] 6.B.A.Peters et al.,Accurate whole-genome sequencing and haplotypingfrom 10to 20human cells.Nature 487,190-195(2012).

[0094] 7.J.Duitama et al.,Fosmid-based whole genome haplotyping of aHapMaptrio child:evaluation of Single Individual Haplotyping techniques.NucleicAcids Res 40,2041-2053(2012).

[0095] 8.S.Selvaraj,R.D.J,V.Bansal,B.Ren,Whole-genome haplotypereconstruction using proximity-ligation and shotgun sequencing.Nat Biotechnol31,1111-1118(2013).

[0096] 9.V.Kuleshov et al.,Whole-genome haplotyping using long reads andstatistical methods.Nat Biotechnol 32,261-266(2014).

[0097] 10.S.Amini et al.,Haplotype-resolved whole-genome sequencing bycontiguity-preserving transposition and combinatorial indexing.Nat Genet 46,1343-1349(2014).

[0098] 11.G.X.Zheng et al.,Haplotyping germline and cancer genomes withhigh-throughput linked-read sequencing.Nat Biotechnol,(2016).

[0099] 12.F.Zhang et al.,Haplotype phasing of whole human genomes usingbead-based barcode partitioning in a single tube.Nat Biotechnol 35,852-857(2017).

[0100] 13.B.A.Peters,J.Liu,R.Drmanac,Co-barcoded sequence reads fromlong DNAfragments:a cost-effective solution for"perfect genome"sequencing.Frontiersin genetics 5,466(2014).

[0101] 14.R.Drmanac.Nucleic Acid Analysis by Random Mixtures of Non-Overlapping Fragments.WO 2006 / 138284 A2(2006).

[0102] 15.R.Drmanac,Peters,B.A.,Alexeev,A.Multiple tagging of longDNAfragments.WO 2014 / 145820 A2(2014).

[0103] 16.K.Jo,Y.L.Chen,J.J.de Pablo,D.C.Schwartz,Elongation andmigration ofsingle DNA molecules in microchannels using oscillatory shear flows.Lab Chip9,2348-2355(2009).

[0104] 17.R.Drmanac et al.,Human genome sequencing using unchained basereadson self-assembling DNA nanoarrays.Science 327,78-81(2010).

[0105] 18.T.Fehlmann et al.,cPAS-based sequencing on the BGISEQ-500toexplore small non-coding RNAs.Clin Epigenetics 8,123(2016).

[0106] 19.J.Huang et al.,A reference human genome dataset of the BGISEQ-500sequencer.Gigascience 6,1-9(2017).

[0107] 20.S.S.T.Mak et al.,Comparative performance of the BGISEQ-500vsIllumina HiSeq2500 sequencing platforms for palaeogenomicsequencing.Gigascience 6,1-13(2017).

[0108] 21.H.Li,R.Durbin,Fast and accurate short read alignment withBurrows-Wheeler transform.Bioinformatics 25,1754-1760(2009).

[0109] 22.A.McKenna et al.,The Genome Analysis Toolkit:a MapReduceframeworkfor analyzing next-generation DNA sequencing data.Genome Res 20,1297-1303(2010).

[0110] 23.J.M.Zook et al.,Integrating human sequence data sets providesaresource of benchmark SNP and indel genotype calls.Nat Biotechnol 32,246-251(2014).

[0111] 24.P.Edge,V.Bafna,V.Bansal,HapCUT2:robust and accurate haplotypeassembly for diverse sequencing technologies.Genome Res 27,801-812(2017).

[0112] 25.Q.Mao et al.,The whole genome sequences and experimentally phasedhaplotypes of over 100personal genomes.Gigascience 5,1-9(2016).

[0113] 26.Z.Dong et al.,Low-pass whole-genome sequencing in clinicalcytogenetics:a validated approach.Genet Med 18,940-948(2016).

[0114] 27.Z.Dong et al.,Identification of balanced chromosomalrearrangements previously unknown among participants in the 1000GenomesProject:implications for interpretation of structural variation in genomesand the future of clinical cytogenetics.Genet Med,(2017).

[0115] 28.J.Ghurye,M.Pop,S.Koren,D.Bickhart,C.S.Chin,Scaffolding of longread assemblies using long range contact information.BMC Genomics 18,527(2017).

[0116] 29.M.Pendleton et al.,Assembly and diploid architecture of anindividual human genome via single-molecule technologies.Nat Methods 12,780-786(2015).

[0117] 30.J.R.Dixon et al.,Topological domains in mammalian genomesidentified by analysis of chromatin interactions.Nature 485,376-380(2012).

[0118] 31.A.M.Phillippy,M.C.Schatz,M.Pop,Genome assembly forensics:findingthe elusive mis-assembly.Genome biology 9,R55(2008).

[0119] 32.JDBuenrostro,PGGiresi,LCZaba,HYChang,WJGreenleaf,Transposition of native chromatin for fast and sensitive epigenomic profilingof open chromatin,DNA-binding proteins and nucleosome position.Nat Methods10,1213-1218(2013).

[0120] Example

[0121] 2. Example 1: Methods and Materials

[0122] 2.1. Isolation of high molecular weight DNA

[0123] Follow RecoverEase TM A modified version of the DNA Isolation Kit (Agilent Technologies, La Jolla, CA) protocol was used to isolate long genomic DNA from cell lines (1).

[0124] Briefly, approximately one million cells were pelleted and lysed with 500 ul of lysis buffer. After incubation at 4°C for 10 minutes, 20 μL of RNase-IT ribonuclease cocktail in 4 mL of digestion buffer was added directly to the lysed cells and incubated on a heat block at 50°C. After 5 minutes, 4.5 mL of proteinase K solution ( Proteinase K, 0.56% SDS and 0.89X TE), and the mixture was incubated for another 2 hours at 50° C. The genomic DNA was then transferred to dialysis tubing with a molecular weight cutoff of 1000 kD (Spectrum Laboratories, Inc., Rancho Dominguez, CA) and dialyzed against 0.5X TE buffer at room temperature overnight.

[0125] 2.2 Barcoded Bead Construction

[0126] Using three sets of double-stranded barcoded DNA molecules, barcoded beads were constructed by a split-and-pool-based strategy. Fig.12 and 13 . A common adapter sequence containing PCR primer annealing sites was ligated to Dynabeads with a 5' double biotin linker. TMM-280 streptavidin (ThermoFisher, Waltham, MA) magnetic beads. Three sets of 1536 barcode oligonucleotides containing overlapping sequence regions were constructed by Integrated DNA Technologies (IA, Coralville). Ligation was performed in a 384-well plate in a 15 μL reaction containing 50 mM Tris-HCl (pH 7.5), 10 mM MgCl2, 1 mM ATP, 2.5% PEG-8000, 571 units of T4 ligase, 580 pmol barcode oligonucleotides and 65 million M-280 beads. Ligation reactions were incubated on a rotator at room temperature for 1 hour. Between ligations, beads were pooled into a single container by centrifugation, collected to the side of the container using a magnet, and washed once with high salt wash buffer (50 mM Tris-HCl (pH 7.5), 500 mM NaCl, 0.1 mM EDTA, and 0.05% Tween 20) and twice with low salt wash buffer (50 mM Tris-HCl (pH 7.5), 150 mM NaCl, and 0.05% Tween 20). The beads were resuspended in 1X ligation buffer and distributed in a 384-well plate, and the ligation step was repeated.

[0127] Some of the "barcodes" referred to herein are "tripartite barcodes". Tripartate refers to their structure and / or their composition. Fig.12 As shown, a triple barcode can be synthesized by continuously connecting shorter (e.g., 4-20 nucleotide) sequences. In one embodiment, the length of the shorter barcode is 10 bases. As shown, an exemplary structure includes CS1-BC1-CS2-BC2-CS3-BC3-CS4, wherein CS is a constant sequence present on all capture adapters, and the BC sequence is a diversified 10-base barcode as described herein. As shown, a partially double-stranded oligonucleotide having a structure CSa-BC-CSb can be used to anneal to a shorter oligonucleotide (i.e., the complementary sequence of BC (i.e., BC')) to construct a triple barcode.

[0128] On the one hand, the present invention provides a composition comprising beads with capture oligonucleotides comprising attached clonal barcodes, wherein the composition comprises more than 3 billion different barcodes, and wherein the barcode is a triple barcode having the structure 5'-CS1-BC1-CS2-BC2-CS3-BC3-CS4. In some embodiments, CS1 and CS4 are longer than CS2 and CS3. In some embodiments, CS2 and CS3 are 4-20 bases, CS1 and CS4 are 5 or 10 to 40 bases, such as 20-30 bases, and the length of the BC sequence is 4-20 bases (e.g., 10 bases). In some embodiments, CS4 is complementary to a splint oligonucleotide. In some embodiments, the composition comprises a bridge oligonucleotide. In some embodiments, the composition comprises a bridge oligonucleotide, a bead comprising a triple barcode as described above, and a genomic DNA comprising a hybridization sequence having a region complementary to the bridge oligonucleotide.

[0129] 2.3 stLFR using two transposons

[0130] 2 pmol of Tn5-coupled transposon was inserted into 40 ng of genomic DNA in a 60 μL reaction of 10 mM TAPS-NaOH (pH 8.5), 5 mM MgCl2, and 10% DMF at 55°C for 10 min. 1.5 μL of transposon-inserted DNA was transferred to 248.5 μL of hybridization buffer consisting of 50 mM Tris-HCl (pH 7.5), 100 mM MgCl2, and 0.05% 20 composition. 10-50 million barcoded beads are resuspended in the same hybridization buffer. Diluted DNA is added to the barcoded beads, and the mixture is heated to 60°C for 10 minutes with occasional slight mixing. The DNA-bead mixture is transferred to a tube revolver in a laboratory oven and incubated at 45°C for 50 minutes. 500 μL of a connection mixture containing 50mM Tris-HCl (pH 7.8), 10mM DTT, 1mM ATP, 2.5% PEG-8000 and 4000 units of T4 ligase is added directly to the DNA-bead mixture. The ligation reaction is incubated on a rotator at room temperature for 1 hour. 110 μL of 1% SDS is added, and the mixture is incubated at room temperature for 10 minutes to remove the Tn5 enzyme. The beads are collected on the side of the tube via a magnet, washed once with low salt wash buffer, and washed once with NEB2 buffer (New England Biolabs, Ipswich, MA). Excess barcode oligonucleotides were removed using 10 units of UDG (MA, Ipswich, New England Biolabs, MA) and 30 units of APE1 (New England Biolabs, Ipswich, MA) and 40 units of Exonuclease 1 (New England Biolabs, Ipswich, MA) in 100 uL 1X NEB2 buffer. The reaction was incubated at 37°C for 30 minutes. The beads were collected on the side of the tube and washed once with low salt wash buffer and then washed once with 1X PCR buffer (1X PfuCx buffer (Agilent Technologies, La Jolla, CA), 5% DMSO, 1M Betaine, 6mM MgSO4 and 600μM dNTPs). The PCR mixture containing 1X PCR buffer, 400pmol of each primer and 6μL PfuCx enzyme (Agilent Technologies, LaJolla, CA) was heated to 95°C for 3 minutes and then cooled to room temperature. This mixture was used to resuspend the beads, and the combined mixture was incubated at 72°C for 10 minutes, followed by 12 cycles of 10 seconds at 95°C, 30 seconds at 58°C, and then 2 minutes at 72°C.

[0131] 2.4 stLFR with 3' branched ligation adapter

[0132] The method starts with the same hybridization insertion conditions, but only one transposon is used instead of two transposons. As mentioned above, after capture and barcode connection steps, beads are collected on the side of the tube and washed with low salt wash buffer. The adapter digestion mixture of 90 units of exonuclease I (New England Biolabs, Ipswich, MA) and 100 units of exonuclease III (New England Biolabs, Ipswich, MA) in 100 μL 1X TA buffer (Teknova, Hollister, CA) is added to the beads and incubated at 37 ° C for 10 minutes. Terminate the reaction, and remove the Tn5 enzyme by adding 11 μL of 1% SDS. Beads are collected on the side of the tube, and washed once with low salt wash buffer, and then washed once with 1X NEB2 buffer (New England Biolabs, Ipswich, MA). Excess capture oligonucleotide was removed by adding 10 units of UDG (New England Biolabs, Ipswich, MA) and 30 units of APE1 (New England Biolabs, Ipswich, MA) in 100 uL of 1X NEB2 buffer (New England Biolabs, Ipswich, MA) and incubating at 37°C for 30 minutes. The beads were collected on the side of the tube and washed once with high salt wash buffer and once with low salt wash buffer. 300 pmol of the second adapter was ligated to the beads bound to the subfragments with 4000 units of T4 ligase in 100 uL of ligase buffer containing 50 mM Tris-HCl (pH 7.8), 10 mM MgCl2, 0.5 mM DTT, 1 mM ATP, and 10% PEG-8000 at room temperature on a rotator for 2 hours. The beads were collected on the side of the tube and washed once in high salt wash buffer and once in 1X PCR buffer. The PCR mixture and conditions were the same as for the two transposon method described above.

[0133] Exemplary 3' branch ligation adapters include 3' branch ligation adapter-F ( / 5Phos / CTGATGGCGCGAGGGAGGC) and 3' branch ligation adapter-R (TCGCGCCATCA / 3'dd / G) oligonucleotides as shown in Table 6. For example, the adapter F sequence includes a PCR primer annealing sequence. Optionally, a barcode (e.g., a sample barcode) can be included between the 5' phosphate and the sequence shown. In this embodiment, the adapter R sequence is shorter than the primer annealing sequence so that it will melt under conditions of PCR primer annealing.

[0134] 2.5 Sequence Mapping and Variant Calling

[0135] First, the raw read data were demultiplexed by associated barcode sequences using the barcode segmentation tool (available from GitHub https: / / github.com / stLFR / stLFR_read_demux). The barcode-assigned and sheared reads were mapped to the hs37d5 reference genome using BWA-MEM (2). The resulting BAM files were then sorted by chromosome coordinates using SAMtools (3) and duplicates were marked using the picard MarkDuplicate function (http: / / broadinstitute.github.io / picard). Short variant (SNP and indel) calling was performed using HaplotypeCaller in GATK4.0.3.0 (4). The vcf files generated in the above steps were then benchmarked against the Genome in a Bottle (GIAB) high-confidence variant list (ftp: / / ftp-trace.ncbi.nlm.nih.gov / giab / ftp / release / NA12878_HG001 / latest / GRCh37 / ) (5) using the rtgtools vcfeval function (6). After benchmarking, the stLFR library was analyzed using GATK VariantRecalibrator and a Gaussian mixture model was trained using the GIAB truth set. The VCFs were then filtered using GATK ApplyVQSR. In almost all cases, 99.9 tranches were applied to the raw vcfs, with the exception of the 100Gb stLFR-1 library and the STD library, where 100 tranches were applied. We then established and applied further hard filtering criteria based on the GQ scores, the ratio of reference depth to alternative depth, and barcode support listed in Table 5:

[0136] 2.6 Variant phasing using Hapcut2

[0137] SNPs were phased using Hapcut2 (https: / / github.com / vibansal / HapCUT2)(7) using its 10XGenomics data pipeline. The BAM files were first converted to a format with barcode information similar to the 10X Genomics barcode BAM. Specifically, a "BX" field was added to each line to reflect the barcode information of that read. GIAB variants or variants called by GATK for each library were used as input for phasing, and the calculate_haplotype_statistics.py tool of Hapcut2 was used to summarize the phasing results and compare them to the GIAB phased vcf file (5).

[0138] 2.7.LongHap

[0139] The seed extension strategy is used during the phasing process of LongHap. It initially starts with a pair of seeds consisting of the most upstream heterozygous variants in the chromosome. The seed is extended by connecting other downstream candidate variants until no more variants can be added to the extended seed ( Fig.11 ). In this expansion process, candidate variants located at different loci will not be treated equally (i.e., upstream variants have higher priority than downstream variants across the chromosome). Every two heterozygous loci have two possible combinations along two different alleles. Take variants T2 / G2 and G3 / C3 as an example ( Fig.11 ), one combination pattern is T2-G3 and G2-C3, and the other is T2-C3 and G2-G3. The score of each combination is calculated by the number of long DNA fragments spanning the two loci, which is equal to the number of unique barcodes with reads mapped to the two loci. Fig.11 As shown, the final score of the former combination is 3, which is three times that of the latter. Variant T2 / G2 is added to the extended seed and the process is repeated. It is worth noting that if any barcode supports two alleles at a particular locus, it will be ignored when calculating the connection score. This helps to reduce the switching error rate. When conflicts occur when connecting downstream candidate variants, such as Fig.11 As shown in the variant A4 / C4 in , a simple decision will be made by comparing the number of connected loci to allow further expansion of candidate variants. In this case, there are two connected loci in the left scenario and only one in the right scenario. LongHap will choose the left combination mode as the final phasing result.

[0140] 2.8SV detection

[0141] As described above, structural variation was detected by counting shared barcodes between genomic regions. Duplicate reads were first removed. The mapped co-barcode reads were scanned using a sliding window along the genome (the default value was 2kb), and each window recorded how many barcodes were found in this 2kb window, and the Jaccard index was calculated for the shared barcode ratio between window pairs. Structural variation events were identified by the Jaccard index sharing metric between window pairs.

[0142] For each window pair (X, Y) on the genome, the Jaccard index is calculated as follows:

[0143] X=(x1,x2,...x n ); Y = (y1, y2, ...y n )

[0144]

[0145] 2.9 Contig scaffolding using SALSA

[0146] Sequencing reads from the stLFR library were used to support the NA12878 assembly containing 18,903 contigs with an NG50 of 26.83 Mb (9) (contigs downloaded from the NCBI genome website using the scaffolding program SALSA (10)). To mimic the HiC sequence structure suitable for SALSA, stLFR sequence reads were selected from fragments of size >= 5 kb. From each fragment of length >= 5 kb, the "first" and "last" reads were selected to form a read pair. These artificial read pairs were then selected by moving the fragments inward at intervals of 2 kb. These read pairs were then mapped to the NA12878 contigs and scaffolded using SALSA. The resulting scaffolds were then aligned and compared to the hg19 reference genome using nucmer and dnadiff of the MUMmer 4 program (11).

[0147] 2.10 References for Example 1

[0148] 1.I.Agent Technologies,RecoverEase DNA Isolation Kit.Revision C.0,(2015).

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

[0150] 3.H.Li et al.,The Sequence Alignment / Map format andSAMtools.Bioinformatics 25,2078-2079(2009).

[0151] 4.A.McKenna et al.,The Genome Analysis Toolkit:a MapReduce frameworkfor analyzing next-generation DNA sequencing data.Genome Res 20,1297-1303(2010).

[0152] 5.J.M.Zook et al.,Integrating human sequence data sets providesaresource of benchmark SNP and indel genotype calls.Nat Biotechnol 32,246-251(2014).

[0153] 6.J.G.Cleary et al.,Comparing Variant Call Files for PerformanceBenchmarking of Next-Generation Sequencing Variant Calling Pipelines.bioRxiv,(2015).

[0154] 7.P.Edge,V.Bafna,V.Bansal,HapCUT2:robust and accurate haplotypeassembly for diverse sequencing technologies.Genome Res 27,801-812(2017).

[0155] 8.F.Zhang et al.,Haplotype phasing of whole human genomes using bead-based barcode partitioning in a single tube.Nat Biotechnol 35,852-857(2017).

[0156] 9. M. Pendleton et al., Assembly and diploid architecture of anindividual human genome via single-molecule technologies. Nat Methods 12, 780-786 (2015).

[0157] 10.J.Ghurye,M.Pop,S.Koren,D.Bickhart,CSChin,Scaffolding of longread assemblies using long range contact information.BMC Genomics 18,527(2017).

[0158] 11.S.Kurtz et al.,Versatile and open software for comparing largegenomes.Genome biology 5,R12(2004).

[0159] Example 2: Detailed experimental protocol

[0160] 3.1 Materials

[0161] 1Kb Plus DNA Ladder (ThermoFisher, Cat. No. 10787018)

[0162] 100Kd MWCO Biotech CE Dialysis Tubing (Spectrum Labs, Cat. No. 131486)

[0163] 384-well Armadillo PCR plate (ThermoFisher, Cat. No. AB2384)

[0164] AMPure XP beads (Beckman Coulter, Cat. No. A63882)

[0165] APE 1 (10000 units / mL) (New England Biolabs, catalog number M0282L)

[0166] ATP (100 mM) (Teknova, Catalog No. A1210)

[0167] Barcoded bead construction oligonucleotides (IDT) (see Notes)

[0168] Betaine (5M) (Sigma-Aldrich, Catalog No. B0300-5VL)

[0169] BSA (20 mg / mL) (New England Biolabs, Catalog No. B9000S)

[0170] Universal Adaptor Oligonucleotide (IDT)

[0171] DMF (~100%) (Sigma-Aldrich, Catalog No. D4551-250ML)

[0172] DMSO (100%) (Sigma-Aldrich, Catalog No. D9170-5VL)

[0173] dNTPs (25 mM) (ThermoFisher, Cat. No. R1121)

[0174] Dialysis tubing (1000kD MWCO) (Spectrum Laboratories, Inc., Cat. No. 131486)

[0175] DTT (Sigma-Aldrich, Catalog No. 11583786001)

[0176] Dynabeads TM M-280 Streptavidin (ThermoFisher, Cat. No. 60210)

[0177] EDTA (0.5M, pH 8.0) (Sigma-Aldrich, Cat. No. 03690-100ML)

[0178] Exonuclease I (20,000 units / mL) (New England Biolabs, Cat. No. M0293L)

[0179] Exonuclease III (100,000 units / mL) (New England Biolabs, Catalog No. M0206L) Formamide (100%, 250 mL) (Sigma-Aldrich, Catalog No. 47671-250ML-F)

[0180] Glycerol (100%) (Sigma-Aldrich, Product No. G5516-100ML)

[0181] KCl (Sigma-Aldrich, Catalog No. P9333-1KG)

[0182] KH2PO4 (Sigma-Aldrich, Catalog No. 795488-1KG)

[0183] KOH (Sigma-Aldrich, Catalog No. P5958-1KG)

[0184] MgCl2 (1M) (Sigma-Aldrich, Catalog No. 63069-500ML)

[0185] MgSO4 (1M) (Sigma-Aldrich, Product No. M3409-100ML)

[0186] MicroAmp transparent adhesive plate sealing film (ThermoFisher, Cat. No. 4306311)

[0187] NaCl (5M) (ThermoFisher, Cat. No. AM9760G)

[0188] Na2HPO4 (Sigma-Aldrich, Catalog No. S7907-1KG)

[0189] NaOH (10M) (Sigma-Aldrich, Catalog No. 72068-100ML)

[0190] NEB2 buffer (10X) (New England Biolabs, Catalog No. B7002S)

[0191] PEG-8000 (50%) (Rigaku, Cat. No. 1008063)

[0192] Pfu Turbo Cx Hotstart DNA Polymerase (Agilent, Cat. No. 600414)

[0193] Proteinase K, recombinant, PCR grade solution (14-22 mg / mL) (Roche, Catalog No. 03115844001) RiboRuler Low Range RNA Ladder (Thermofisher, Catalog No. SM1831) RNase-IT Ribonuclease Mix (Agilent, Catalog No. 400720)

[0194] SDS (10%) (ThermoFisher, Cat. No. 15553027)

[0195] Sucrose (Sigma-Aldrich, Catalog No. S7903-1KG)

[0196] T4 DNA ligase (2x10 6Unit / mL) (New England Biolabs, Catalog No. M0202M) TA buffer (10X) (Teknova, Catalog No. T0379)

[0197] TAPS-NaOH (1M, pH 8.5) (Boston BioProducts, Cat. No. BB-2375)

[0198] TBE (10X) (ThermoFisher, Cat. No. 15581028)

[0199] TE buffer (10X) (Fisher Scientific, Cat. No. BP24771)

[0200] Tn5 enzyme

[0201] Transposon oligonucleotides (IDT)

[0202] Tris-HCl (1M, pH 7.5) (ThermoFisher, Cat. No. 15567027)

[0203] Tris-HCl (2M, pH 7.8) (Amresco, product number J837-500ML)

[0204] Triton TM X-100 (10%) (Sigma-Aldrich, Catalog No. 93443-100ML)

[0205] 20(10%) (Roche, Cat. No. 11332465001)

[0206] UDG (5,000 units / mL) (New England Biolabs, Cat. No. M0280L)

[0207] 3.2 Equipment

[0208] 2.4L tall polystyrene container (Click Clack, Cat. No. 659030) or equivalent

[0209] DynaMag TM -2 magnet (ThermoFisher, Cat. No. 12321D)

[0210] Easy 50EasySep TM Magnet (Stem Cell Technologies, Cat. No. 18002) or equivalent

[0211] Laboratory oven that can accommodate a tube rotator / rotator

[0212] Magnetic Plate Stirrer

[0213] Medium Magnetic Stirring Bar

[0214] Standard laboratory vortex mixer

[0215] Tetrad PCR Thermal Cycler (Bio-Rad, Cat. No. PTC0240) or equivalent, with a reaction volume of up to 100 μL per well

[0216] Tube Rotator / Spinner (Thermo Fisher, Cat. No. 88881001) or equivalent

[0217] 3.3 Reagent Setup

[0218] Annealing buffer (3X)

[0219] 3 mL of 1M Tris-HCl, pH 7.5

[0220] 6 mL of 5 M NaCl

[0221] 91 mL of sterile dH2O

[0222] Store at room temperature for 1 year.

[0223] 3.4 Buffer D (10X)

[0224] 224 mg KOH

[0225] 50 μL 0.5 M EDTA

[0226] 2.45 mL of sterile dH2O

[0227] Aliquot and store at -20 °C for 1 month.

[0228] 3.5 Coupling buffer (1X)

[0229] 5 mL of 1X TE

[0230] 5 mL of 100% glycerol

[0231] Store at -20°C for 1 year.

[0232] 3.6 Digestion buffer (1X, pH 8.0)

[0233] 1.75 g Na2HPO4

[0234] 0.2 g KCl

[0235] 0.2 g KH2PO4

[0236] 27.4 mL 5M NaCl

[0237] 20 mL of 0.5 M EDTA (pH 8.0)

[0238] 800 mL of sterile dH2O

[0239] Adjust pH to 8.0 with 1 M NaOH.

[0240] Add sterile dH2O to a final volume of 1 L.

[0241] Filter sterilize.

[0242] Store at room temperature for 1 year.

[0243] 3.7 3' branch ligation buffer (3X)

[0244] 6 mL of 50% PEG-8000

[0245] 0.75 mL 2M Tris-HCl (pH 7.8)

[0246] 0.3mL 1M MgCl2

[0247] 0.3mL 0.1M ATP

[0248] 15 μL 1M DTT

[0249] 75 μL 20 mg / mL BSA

[0250] 2.560 mL of sterile dH2O

[0251] Store at -20°C for 1 year.

[0252] 3.8 High Salt Bead Binding Buffer (1X)

[0253] 5 mL of 1M Tris-HCl (pH 7.5)

[0254] 6 mL of 5 M NaCl

[0255] 20 μL 0.5M EDTA

[0256] 88.98 mL of sterile dH2O

[0257] Store at room temperature for 1 year.

[0258] 3.9 High Salt Wash Buffer (1X)

[0259] 5 mL of 1M Tris-HCl, pH 7.5

[0260] 10 mL of 5 M NaCl

[0261] 20 μL 0.5M EDTA

[0262] 0.5mL of 10% 20

[0263] 84.48 mL of sterile dH2O

[0264] Store at room temperature for 1 year.

[0265] 3.10 Hybridization buffer (1X)

[0266] 50 mL of 1M Tris-HCl, pH 7.5

[0267] 100 mL of 1M MgCl2

[0268] 10% of 5mL 20

[0269] 845 mL water

[0270] Store at room temperature for 1 year.

[0271] 3.11 Ligation buffer (10X)

[0272] 25 mL of 50% PEG-8000

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

[0274] 5 mL of 100 mM ATP

[0275] 5 mL of 1M MgCl2

[0276] 2.5 mL of sterile dH2O

[0277] Store at -20°C for 1 year.

[0278] 3.12 Ligation buffer, without MgCl2 (10X)

[0279] 25 mL of 50% PEG-8000

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

[0281] 5 mL of 100 mM ATP

[0282] 5 mL of 1M DTT

[0283] 2.5 mL of sterile dH2O

[0284] Store at -20°C for 1 year.

[0285] 3.13 Low salt wash buffer (1X)

[0286] 5 mL of 1M Tris-HCl, pH 7.5

[0287] 3 mL of 5 M NaCl

[0288] 0.5mL of 10% 20

[0289] 91.5 mL sterile dH2O

[0290] Store at room temperature for 1 year.

[0291] 3.14 Lysis buffer (1X, pH 8.3)

[0292] 0.22 g KCl

[0293] 120g cane sugar

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

[0295] 2 mL of 0.5 M EDTA (pH 8.0)

[0296] 28 mL of 5M NaCl

[0297] 10mL X-100

[0298] 800 mL of sterile dH2O

[0299] Adjust pH to 8.3

[0300] Add sterile dH2O to a final volume of 1 L.

[0301] Filter sterilize.

[0302] Store at 4°C for 1 year.

[0303] 3.15 Transposase Buffer (5X)

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

[0305] 0.25 mL of 1 M MgCl2

[0306] 5 mL of 100% DMF

[0307] 4.25 mL of sterile dH2O

[0308] Store at -20°C for 1 year.

[0309] 3.16PfuCx mixture (2X)

[0310] 2 mL of 10X PfuCx buffer (contains enzyme)

[0311] 0.5 mL of 100% DMSO

[0312] 2 mL of 5 M betaine

[0313] 60 μL 1M MgSO4

[0314] 240 μL of 25 mM dNTPs

[0315] 5.2 mL of sterile dH2O

[0316] 3.16 Barcoded Bead Construction Oligonucleotides

[0317] All barcoded oligonucleotides were synthesized in a 384-well format at a 100 nmol scale by standard desalting methods and delivered by Integrated DNA Technologies (Coralville, IA) at a concentration of 200 μM in 1X TE (pH 8.0). There are a total of 1536 unique barcode oligonucleotides per barcode set, and there are 3 barcode sets. This makes up to 3.6 billion different barcode combinations possible. This may not be necessary for some applications, and fewer barcode combinations can be achieved by ordering fewer oligonucleotide plates. This particular design does require the use of at least one barcode oligonucleotide from each set to produce an appropriate final sequence, however, slight modifications to the 6 base overlap sequences between barcode sets can be made to remove entire barcode sets.

[0318] 3.2 Procedure

[0319] Isolation of high molecular weight DNA from cells

[0320] This method is based on RecoverEase TM DNA Isolation Kit Protocol 26, but using much larger volumes to reduce the viscosity of the resulting solution.

[0321] 1. Precipitate up to 1 × 10 in a 15 or 50 mL conical tube. 7 Disperse nucleated cells (500 x g for 5 minutes). Remove supernatant. Add 500 μL of lysis buffer to the cell pellet and briefly vortex the sample at medium speed for 3-5 seconds, then place the conical tube in the refrigerator for approximately 10 minutes with occasional swirling.

[0322] 2. Prepare proteinase K solution by combining 250 μL of 10% SDS, 250 μL of proteinase K, and 4 mL of 1X TE. Place on a heating block at 50°C and heat briefly.

[0323] 3. Prepare digestion solution by mixing 20 μL of RNase-It Ribonuclease Mix with 4 mL of Digestion Buffer.

[0324] 4. Add approximately 4 mL of the prepared digestion solution to the lysed cells and buffer from step 1 and gently shake the conical tube.

[0325] After 5.5 minutes, place the conical tube in a 50°C heating block and add 4.5 mL of heated proteinase K solution to the free-floating pellet. Gently swirl the conical tube to mix.

[0326] 6. Re-cap the tubes and incubate in a heat block at 50 °C for 2 h, gently swirling the tubes every 30 min.

[0327] 7. Cut about 13 cm of dialysis tubing (its capacity is about 1 mL / cm). Allow to equilibrate in 0.5X TE for 30 minutes. Seal one end with a dialysis clip.

[0328] 8. Pour at least 1 L of 0.5X TE buffer into the dialysis cell.

[0329] 9. Carefully pour the viscous genomic DNA from the conical tube into the open end of the dialysis tubing. Seal the open end of the dialysis tubing with a dialysis clamp. Secure the float to one of the clamps. Place the dialysis tubing with the float into the dialysis tank.

[0330] 10. Dialyze the genomic DNA for 24 to 48 hours at room temperature while gently stirring the buffer with a magnetic stir bar. Change the TE buffer once during the dialysis period to maximize the purity of the recovered DNA.

[0331] 11. After dialysis is complete, remove the dialysis tubing from the TE buffer, remove the float and clamp from the top of the dialysis tubing, and gently pour into a 15mL conical tube. The DNA can be used immediately without shearing.

[0332] 3.3 Barcoded beads

[0333] The barcoded beads were constructed using a split-and-merge strategy of three sets of double-stranded barcoded DNA molecules. The full-length adapter was constructed by continuous ligation ( Fig.12 and Fig.13). The barcode oligonucleotides are provided in 384-well plates (see "Reagent Instructions"). The common adapter oligonucleotides are provided in tubes. Depending on the sequencing technology used, it may be necessary to change the PCR primer sequences within the common adapter oligonucleotides.

[0334] 12. Mix 10 μL of complementary oligonucleotide from each well of the source 384-well plate with 10 μL of 3X Annealing Buffer in a 384-well PCR plate. Mix 30 μL of common adaptor oligonucleotide in one well of an 8-well PCR strip.

[0335] 13. Incubate at 70°C for 3 minutes on a PCR thermal cycler, then slowly increase the temperature to 20°C at a rate of 0.1°C / s. The final concentration of the hybridized barcode oligonucleotide is 66 μM.

[0336] 14. Mix 4.725 mL (157.5 μmol) of hybridized bead linker containing 5' bi-biotin with 3.225 mL of ligation buffer (10X), 460.8 μL (921,600 units) of T4 DNA ligase, and 9.67 mL of dH2O to a total volume of 18.081 mL.

[0337] 15. 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 first barcode plate that was hybridized to each well containing the bead linker ligation mixture. Seal with MicroAmp clear adhesive sealer, vortex, centrifuge, and incubate at room temperature for 1 hour.

[0338] 16. Collect 100 billion (143 mL) of M-280 Streptavidin-coated magnetic beads by transferring 50 mL of beads to an empty 50 mL centrifuge tube. Place the 50 mL tube with beads in the Easy 50 EasySep TM Magnet for 5 minutes to collect the beads to the side of the tube. Carefully remove the supernatant with a pipette. Transfer the second 50 mL of beads to the tube on the magnet. Let stand on the magnet for 5 minutes, then carefully remove the supernatant. Transfer the final 43 mL of beads to the 50 mL tube. Let stand on the magnet for 5 minutes, then carefully remove the supernatant. Wash the beads twice with low salt wash buffer, then thoroughly resuspend in 8 mL of high salt bead binding buffer.

[0339] 17. Dispense 5 μL of beads in high salt bead binding buffer into each well of the plate containing the ligated product. During the dispensing process, occasionally vortex the bead source tube to keep the beads well suspended.

[0340] 18. Seal the plate with MicroAmp clear adhesive sealer, vortex, and incubate on a tube rotator in the "shake" mode at room temperature for 1 hour.

[0341] 19. Centrifuge the plate at 300 x g for 5 seconds to remove beads from the seal, but do not allow pelleting. Remove the seal and add 2.8 μL of 0.1% SDS to each well. Seal the plate again with MicroAmp clear adhesive sealer, vortex briefly and incubate at room temperature for 10 minutes.

[0342] 20. Vortex and centrifuge the plate at 300 x g for 5 seconds to remove the beads from the plate seal. Remove the seal from each plate and invert the plate onto a collection plate. Centrifuge at 500 x g for 2 minutes. Using a 10 ml serological pipette, collect the beads into a new 50 ml tube.

[0343] 21. Collect the beads into Easy 50EasySep TM Place the tube sideways on the magnet for 5 minutes. Discard the supernatant. Wash once with 10 mL of high salt wash buffer, then wash twice with low salt wash buffer. Resuspend the beads in 8 mL of 1X ligation buffer.

[0344] 22. Dispense 5ul of beads into each well of four new 384-well PCR plates. During the dispensing process, vortex the bead source tube occasionally to keep the beads well suspended.

[0345] 23. To ligate the second set of barcodes, bring the total volume of the mixture containing 3.225 mL of ligation buffer (10X), 460.8 μL (921,600 units) of T4 DNA ligase, and 6.33 mL of dH2O to 10.02 mL. Dispense 6.2 μL of the second ligation mixture into each well of four 384-well PCR plates containing beads. Next, add 8.8 μL (580 pmol) from each well of the hybridized second barcode plate to the corresponding well of the 384-well PCR plate containing beads and ligation mixture.

[0346] 24. Repeat steps 18-22.

[0347] 25. To ligate the third set of barcodes, bring the total volume of the ligation mix to 10.02 mL, containing 3.225 mL ligation buffer (10X), 460.8 μL (921600 units) of T4 DNA ligase, and 6.33 mL dH2O. Dispense 6.2 μL of the third ligation mix into each well of four 384-well PCR plates containing beads. Next, add 8.8 μL (580 pmol) from each well of the hybridized third barcode plate to the corresponding well of the 384-well PCR plate containing beads and ligation mix.

[0348] 26. Repeat steps 18-22. The beads can now be stored at 4°C for up to one year. In their current form, the beads are almost completely double stranded and are not yet in the correct form for use in stLFR.

[0349] 27. Count the beads using a hemacytometer and remove 5 million beads for the QC step. Place the tube with beads in the DynaMag TM -2 magnet for 5 minutes. Discard supernatant. Add 5 μL 100% formamide, 4 μL dH2O, and 1 μL 10X loading buffer. Incubate at 95°C for 3 minutes on a PCR thermocycler. Immediately place on ice for 2 minutes. Place the tube with beads on a DynaMag TM -2 magnet for 5 minutes. Collect the supernatant, load onto a 15% TBU gel, and run at 200V for 40 minutes to check the length and quantity of the oligonucleotides. Alternatively, the beads can be checked using a flow cytometer by hybridizing a fluorescently labeled oligonucleotide to the 3' end of the bead adapter sequence. We typically see about 25% of the total streptavidin binding sites have the full-length constructed adapter sequence.

[0350] 3.20stLFR Bead Preparation

[0351] To prepare beads for stLFR, they must first be denatured into single-stranded DNA and then rehybridized with the bridge oligonucleotide.

[0352] 28. Pipette 500 million constructed barcoded beads from step 26 of the previous section into a standard 1.5 mL microcentrifuge tube.

[0353] 29. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant.

[0354] 30. Add 1 mL of the 1X dilution of Buffer D. Vortex briefly and incubate at room temperature for 2 minutes.

[0355] 31. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant.

[0356] 32. Repeat steps 30 and 31 once more.

[0357] 33. Wash once in 1X annealing buffer. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant.

[0358] 34. Mix 36 μL of 100 μM Bridge Oligo, 333.33 μL of Annealing Buffer (3X), and 630.67 μL of dH2O to a final volume of 1 mL. Add the mixture to the beads. Vortex briefly.

[0359] 35. Incubate at 60°C for 5 minutes and at room temperature for 50 minutes.

[0360] 36. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and resuspend in 500 μL low salt wash buffer. The beads are now ready for stLFR and can be stored at 4 °C for 3 months.

[0361] 3.21 Two transposon stLFR experimental protocols

[0362] This protocol utilizes two transposons to generate hybridization sequences and PCR primer sites along the length of a genomic DNA molecule. This is the most simplified and fastest stLFR method, but coverage of each long DNA fragment may be reduced by 50% compared to the 3' branch ligation protocol. For compatibility with sequencing technologies other than BGISEQ-500, it may be necessary to change some of the transposon sequences after the chimeric region. Before ordering these oligonucleotides, check the sequencing primers used. Information on all oligonucleotide sequences is available in the supplementary material.

[0363] 37. Hybridize the captured transposon oligonucleotides by combining 10 μL of Transposon 1T (100 μM), 10 μL of Transposon B (100 μM), 10 μL of Annealing Buffer (3X) into the first well of an 8-well PCR strip tube. Hybridize the non-captured transposon oligonucleotides by combining 10 μL of Transposon 1T (100 μM), 10 μL of Transposon B (100 μM), 10 μL of Annealing Buffer (3X) into the second well of the same PCR strip tube.

[0364] 38. Incubate at 70°C for 3 minutes on a PCR thermocycler, then slowly increase the temperature to 20°C at 0.1°C / s. Combine the two transposons into the third well of the PCR strip tube.

[0365] 39. Couple Tn5 enzyme to the transposon mix by combining 9.6 μL of the mixed transposon with 23.53 μL of Tn5 (13.6 pmol / μL) and 46.87 μL of coupling buffer (1X).

[0366] 40. 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 recommend aliquoting before storage.

[0367] 41. Incorporate the transposon into the long genomic DNA by combining 12 μL of Transposase Buffer (5X), 0.5 μL of the coupled transposon from step 40, and 40 ng of DNA in one well of an 8-well strip tube for a total volume of 60 μL. NOTE: The amount of DNA and the amount of coupled transposon can be adjusted at this step. Because there can be variability between batches, it will be necessary to titrate the amount of Tn5 enzyme used. Again, it is possible to start with less DNA, but for titration purposes it is useful to use 40 ng so that some material can be run on an agarose gel to determine the efficiency of transposon incorporation (see later steps).

[0368] 42. Incubate at 55 °C for 10 minutes.

[0369] 43. Transfer 40 μL of the transposon-incorporated material to one well of a new 8-well tube strip. Add 4 μL of 1% SDS and incubate at room temperature for 10 minutes.

[0370] 44. Load the material from step 43 on a 0.5X TBE 1% agarose gel and run at 150V for 40 minutes. The transposed DNA should be between 200 and 1500bp on the gel. We typically expect to see the brightest portion of the DNA smear around 600bp, this may vary depending on the sequencing technology chosen. We typically load a control that has gone through the same steps but lacks the transposon, Tn5 enzyme, or genomic DNA. If the transposon integration product looks the correct size, proceed to step 45. Otherwise, repeat the above steps but adjust the concentration of the coupled product until the smear is the desired size.

[0371] 45. Dilute 1.5 μL of the remaining product from step 42 with 248.5 μL of 1x hybridization buffer.

[0372] 46. ​​Transfer 50 μL of beads (50 million) from step 36 to a 1.5 mL microcentrifuge tube. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and resuspend in 250 μL of hybridization buffer (1X).

[0373] 47. Heat the DNA and beads separately at 60°C for 30 seconds.

[0374] 48. Add 250 μL of diluted DNA to 250 μL of beads, gently mix by tapping the bottom of the tube with your finger, and continue incubating at 60°C for 10 minutes. Gently mix the tube with your finger every few minutes.

[0375] 49. Place on a tube rotator and incubate in an oven at 45 °C in "shake" mode for 50 minutes.

[0376] 50. Prepare the ligation buffer by combining 100 μL of ligation buffer, 10X MgCl2-free, 2 μL of T4 DNA ligase (2x10 6 Prepare the ligation mix by adding 1 mL of 5% dHO (100 μL / mL) and 398 μL of dH 2 O. Remove the tube from the rotator and add the ligation mix to a total volume of 1 mL.

[0377] 51. Incubate on a tube rotator in "shaking" mode for 1 hour at room temperature.

[0378] 52. Add 110 μL of 1% SDS to the tube and incubate at room temperature for 10 minutes.

[0379] 53. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and wash once with 500 μL low salt wash buffer and once with 500 μL NEB2 buffer (1X).

[0380] 54. Prepare capture oligonucleotide digestion mix by combining 10 μL NEB2 buffer (10X), 2 μL UDG (5,000 U / mL), 3 μL APE1 (10,000 U / mL), 2 μL Exonuclease 1 (20,000 units / mL), and 83 μL dH 2 O. Remove wash buffer and add digestion mix to beads.

[0381] 55. Vortex gently to resuspend the beads and incubate at 37°C for 30 minutes.

[0382] 56. Place on DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and wash once with 500 μL low salt wash buffer and once with 500 μL PfuCx buffer (1X).

[0383] 57. Prepare a PCR premix by adding 150 μL PCR mix (2X), 4 μL PCR primer 1 (100 μM), 4 μL PCR primer 2 (100 μM), 6 μL PfuCx enzyme, and 136 μL dH2O. Preheat the PCR premix at 95°C for 3 minutes. Place in a DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the wash buffer and add the PCR premix to the beads.

[0384] 58. Gently vortex to resuspend the beads and cycle the PCR reaction under the following conditions:

[0385]

[0386] 59. PCR should produce approximately 500 ng of DNA, run 20 ng of the 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.

[0387] 60. Purify the PCR product using 300 μL of Agencourt XP beads according to the manufacturer's protocol. This purified product is now ready for sequencing.

[0388] Single transposon 3' branch ligation stLFR protocol

[0389] This protocol is based on single transposon insertion in a DNA nick and a novel adapter ligation method and can result in higher coverage per fragment, which can be important for certain sequencing strategies such as de novo assembly. This strategy is slightly more expensive due to the addition of additional reagents. It also takes more than 2.5 hours.

[0390] 61. Hybridize the captured transposon oligonucleotides by combining 10 μL Transposon 1T (100 μM), 10 μL Transposon B (100 μM), 10 μL Annealing Buffer (3X) into the first well of an 8-well PCR strip. Hybridize the gap ligation adapters by combining 10 μL Branch T (100 μM), 10 μL Branch B (100 μM), 10 μL Annealing Buffer (3X) into the second well of the same PCR strip.

[0391] 62. Incubate at 70°C for 3 minutes on a PCR thermal cycler, then slowly increase the temperature to 20°C at 0.1°C / s.

[0392] 63. Enzymatically couple the Tn5 to the transposon by combining 9.6 μL of the hybridization captured transposon in step 61 with 23.53 μL of Tn5 (13.6 pmol / μL) and 46.87 μL of coupling buffer (1X).

[0393] 64. Incubate at 30°C for 1 hour. Use immediately or store at -20°C for up to 1 month.

[0394] 65. Execute steps 41-51.

[0395] 66. Place in DynaMag TM -2 magnet for 2 min to collect the beads to the side of the tube. Remove the supernatant and wash once with 500 μL low salt wash buffer.

[0396] 67. Prepare the adapter oligonucleotide digestion mix by combining 10 μL TA buffer (10X), 4.5 μL Exonuclease I (20,000 U / mL), 1 μL Exonuclease III (100,000 U / mL), and 74.5 μL dH 2 O. Remove the wash buffer and add the digestion mix to the beads.

[0397] 68. Vortex gently to resuspend the beads and incubate on a tube rotator at 37 °C for 10 minutes on the "shake" mode.

[0398] 69. Add 11 μL of 1% SDS and incubate at room temperature for 10 minutes.

[0399] 70. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and wash once with 500 μL low salt wash buffer and once with 500 μL NEB2 buffer (1X).

[0400] 71. Prepare capture oligonucleotide digestion mix by combining 10 μL NEB2 buffer (10X), 2 μL UDG (5000 U / mL), 3 μL APE1 (10000 U / mL), and 85 μL dH 2 O. Remove wash buffer and add digestion mix to beads.

[0401] 72. Vortex gently to resuspend the beads and incubate at 37°C for 30 minutes.

[0402] 73. Place in DynaMag TM -2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and wash once with 500 μL high salt wash buffer and once with 500 μL low salt wash buffer (1X).

[0403] 74. Prepare the 3' branch ligation buffer (3X) by combining 33.4 μL of 3' branch ligation buffer (3X), 18 μL of 3' branch ligation adapter (16.7 μM) prepared in step 61, 2 μL of T4 DNA ligase (2x10 6 Prepare the 3' branch ligation mix by adding 46.6 μL of dH 2 O. Remove the wash buffer and add the ligation mix to the beads.

[0404] 75. Vortex gently to resuspend the beads and incubate on a tube rotator at 25 °C for 2 hours on the "shake" mode.

[0405] 76. Place in DynaMag TM-2 magnet for 2 minutes to collect the beads to the side of the tube. Remove the supernatant and wash once with 500 μL high salt wash buffer and once with 500 μL PCR buffer (1X).

[0406] 77. Prepare a 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 dH 2 O. Remove the wash buffer and add the PCR master mix to the beads.

[0407] 78. Gently vortex to resuspend the beads and cycle the PCR reaction under the following conditions:

[0408]

[0409] 79. Perform steps 59-60 above.

[0410] 3.4 Analysis of stLFR data

[0411] 7 The starting point for this process is a FASTQ file. This is the standard format for read data generated by most sequencing technologies. The software we use to deconvolute the barcode information takes a FASTQ file and expects the 42 bases of the barcode and the universal adapter sequence to be appended to the end of the first read. It matches the data of the barcode read to the expected 1536 sequences at each barcode position. The barcoding strategy used by stLFR can error-correct barcodes with single base mismatches. The final output of our software is a FASTQ file with the barcode information appended to the end of the read ID in the format #Barcode1ID_Barcode2ID_Barcode3ID, where BarcodeID is a number from 0-1536. A barcode ID of zero means it does not match any of the expected barcode sequences. We recommend using BWA-mem27 for mapping, GATK28 for variant calling, and HapCUT229 for phasing. We also recommend using a decoy sequence to map to Hg19.

[0412] 3.5 References for Example 2

[0413] 1 Zhang, K. et al. Long-range polony haplotyping of individual humanchromosome molecules. Nat Genet 38, 382-387 (2006).

[0414] 2Ma,L.et al.Direct determination of molecular haplotypes bychromosome microdissection.Nat Methods 7,299-301(2010).

[0415] 3 Kitzman,J.O.et al.Haplotype-resolved genome sequencing of aGujaratiIndian individual.Nat Biotechnol 29,59-63(2011).

[0416] 4 Suk,E.K.et al.A comprehensively molecular haplotype-resolvedgenomeof a European individual.Genome Res 21,1672-1685(2011).

[0417] 5 Fan,H.C.,Wang,J.,Potanina,A.&Quake,S.R.Whole-genomemolecularhaplotyping of single cells.Nat Biotechnol 29,51-57(2011).

[0418] 6 Peters,B.A.et al.Accurate whole-genome sequencing andhaplotypingfrom 10 to 20 human cells.Nature 487,190-195(2012).

[0419] 7 Duitama,J.et al.Fosmid-based whole genome haplotyping of aHapMaptrio child:evaluation of Single Individual Haplotyping techniques.NucleicAcids Res 40,2041-2053(2012).

[0420] 8 Selvaraj,S.,J,R.D.,Bansal,V.&Ren,B.Whole-genomehaplotypereconstruction using proximity-ligation and shotgun sequencing.NatBiotechnol 31,1111-1118(2013).

[0421] 9 Kuleshov,V.et al.Whole-genome haplotyping using long readsandstatistical methods.Nat Biotechnol 32,261-266(2014).

[0422] 10 Amini,S.et al.Haplotype-resolved whole-genome sequencingbycontiguity-preserving transposition and combinatorial indexing.Nat Genet46,1343-1349(2014).

[0423] 11 Zheng,G.X.et al.Haplotyping germline and cancer genomes withhigh-throughput linked-read sequencing.Nat Biotechnol(2016).

[0424] 12 Zhang,F.et al.Haplotype phasing of whole human genomes usingbead-based barcode partitioning in a single tube.Nat Biotechnol 35,852-857(2017).

[0425] 13 Peters,B.A.,Liu,J.&Drmanac,R.Co-barcoded sequence reads fromlongDNA fragments:a cost-effective solution for"perfect genome"sequencing.Frontiers in genetics 5,466(2014).

[0426] 14 Drmanac,R.Nucleic Acid Analysis by Random Mixtures of Non-Overlapping Fragments.WO 2006 / 138284 A2(2006).

[0427] 15 McElwain,M.A.,Zhang,R.Y.,Drmanac,R.&Peters,B.A.LongFragment Read(LFR)Technology:Cost-Effective,High-Quality Genome-WideMolecularHaplotyping.Methods Mol Biol 1551,191-205(2017).

[0428] 16 Schaaf,C.P.et al.Truncating mutations of MAGEL2 cause Prader-Williphenotypes and autism.Nat Genet 45,1405-1408(2013).

[0429] 17 Peters,B.A.et al.Detection and phasing of single base denovomutations in biopsies from human in vitro fertilized embryos by advancedwhole-genome sequencing.Genome Res 25,426-434(2015).

[0430] 18 Ciotlos,S.et al.Whole genome sequence analysis of BT-474usingcomplete Genomics'standard and long fragment readtechnologies.Gigascience 5,8(2016).

[0431] 19 Hellner,K.et al.Premalignant SOX2 overexpression in thefallopiantubes of ovarian cancer patients:Discovery and validationstudies.EBioMedicine10,137-149(2016).

[0432] 20 Mao,Q.et al.The whole genome sequences and experimentallyphasedhaplotypes of over 100 personal genomes.Gigascience 5,1-9(2016).

[0433] 21 Gulbahce,N.et al.Quantitative Whole Genome SequencingofCirculating Tumor Cells Enables Personalized Combination Therapy ofMetastaticCancer.Cancer Res 77,4530-4541(2017).

[0434] 22 Walker,R.F.et al.Clinical and genetic analysis of a raresyndromeassociated with neoteny.Genetics In Medicine(2017).

[0435] 23 Mao,Q.et al.Advanced Whole-Genome Sequencing and Analysis ofFetalGenomes from Amniotic Fluid.Clinical chemistry(2018).

[0436] 24 Drmanac,R.,Peters,B.A.,Alexeev,A.Multiple tagging ofindividuallong DNA fragments.WO 2014 / 145820 A2(2013).

[0437] 25 Picelli,S.et al.Tn5 transposase and tagmentation proceduresformassively scaled sequencing projects.Genome Res 24,2033-2040(2014).

[0438] 26 Agent Technologies,I.RecoverEase DNA Isolation Kit.Revision C.0(2015).

[0439] 27Li,H.&Durbin,R.Fast and accurate short read alignment with Burrows-Wheeler transform.Bioinformatics 25,1754-1760(2009).

[0440] 28McKenna,A.et al.The Genome Analysis Toolkit:a MapReduce frameworkfor analyzing next-generation DNA sequencing data.Genome Res 20,1297-1303(2010).

[0441] 29Edge,P.,Bafna,V.&Bansal,V.HapCUT2:robust and accurate haplotypeassembly for diverse sequencing technologies.Genome Res 27,801-812(2017).

[0442] The BAM file "NA12878_WGS_v2_phased_possorted_bam.bam" from a recent Chromium dataset was downloaded from the 10X Genomics website and processed in the same way as the stLFR library. For filtering the results, we used the VCF file "NA12878_WGS_v2_phased_variants.vcf.gz" from the same Chromium library. This VCF contains data processed through the 10X Genomics optimization pipeline. The fragment sizes of the Chromium library were copied from the 10X Genomics website. 10Genomics calculates fragment sizes using a length-weighted average, which may be larger than the average fragment size. 2 No read data are available yet, as reported by Zhang et al. (12). 3 Data from a standard library processed on BGISEQ-500.

[0443] Table 6 shows exemplary sequences that can be used in the stLFR methods described herein.

[0444] Example 3: 3' branch ligation, a novel method for ligating DNA to the 3'OH terminus of DNA or RNA and its application

[0445] 4.1 Introduction

[0446] This example generally describes 3' branch ligation. In the stLFR embodiments described herein, 3' branch ligation is used to add additional adaptors (3' branch ligation adaptors). See, eg, §1.1.2.

[0447] Ligases join breaks in nucleic acids, which is essential for cell viability and vitality. DNA ligases catalyze the formation of phosphodiester bonds between DNA ends and play a key role in DNA repair, recombination, and replication in vivo. RNA ligases join 5'-phosphoryl (5'PO4) and 3'-hydroxyl (3'OH) RNA ends via phosphodiester bonds and are involved in RNA repair, splicing, and editing. Organisms from all three kingdoms (Bacteria, Archaea, and Eukaryota) can be used as important molecular tools in vitro for applications such as cloning, ligase-based amplification or detection, and synthetic biology.

[0448] One of the most widely used ligases in vitro is the bacteriophage T4 DNA ligase, which is a single 55-kDA polypeptide that requires ATP as an energy source. T4 DNA ligase typically ligates adjacent 5'PO4 and 3'OH termini of double-stranded DNA. In addition to sealing nicks or ligating sticky ends, T4 DNA ligase can also efficiently catalyze blunt-end ligation, which is not seen in all other DNA ligases. Some unusual catalytic properties of this ligase have been previously reported, such as sealing single-stranded gaps in double-stranded DNA, sealing nicks near abasic sites in double-stranded DNA (dsDNA), promoting intramolecular loop formation of partially double-stranded DNA, and ligating DNA strands containing 3' branch extensions. (Nilsson and Magnusson, Nucleic Acids Res 10:1425-1437, 1982; Goffin et al., Nucleic Acids Res 15:8755-8771, 1987; Mendel-Hartvig et al., Nucleic Acids Res. 32:e2, 2004; Western and Rose, Nucleic Acids Res.,19:809-813,1991). Researchers have also observed template-independent ligations mediated by T4 ligases, such as sealing of mispaired nicks in dsDNA (Alexander, 2003, Nucleic Acids Res. 2003 Jun 15; 31(12): 3208-16) and even single-stranded DNA (ssDNA) ligations, although with very low efficiency (H. Kuhn, 2005, FEBS J. 2005 Dec; 272(23): 5991-6000). These results suggest that perfect complementary base pairing at or near the ligation junction is not critically required for some unconventional T4 DNA ligase activities. T4 RNA ligases 1 and 2 are products of T4 bacteriophage genes 63 and 24, respectively. They both require adjacent 5'PO4 and 3'OH ends to successfully hydrolyze ATP to AMP and PPi. The substrates of T4 RNA ligase 1 include single-stranded RNA and DNA, while T4 RNA ligase 2 preferentially seals nicks on dsRNA rather than ligating the ends of ssRNA.

[0449] Here, we demonstrate an unconventional end-joining event mediated by T4 DNA ligase, which we term 3'-branched ligation (3'BL). It can join DNA or DNA / RNA fragments at nicks, single-stranded gaps, or 5'-overhang regions to form branched structures. This report extensively studies various ligation cofactors and activators and optimizes the ligation conditions for this novel ligation. Using our 3'BL protocol, no base pairing is required, and ligation can be completed in more than 90% of cases even in the presence of 1nt gaps. One of its applications is the ligation of adapters to DNA or RNA in NGS library preparation. Several genomic structures that were previously considered unligatable can now become substrates for 3'BL, resulting in high conversion rates of input DNA to adapter-ligated molecules while avoiding chimeras. We demonstrate that 3'BL can be combined with transposon insertion. The directed transposon insertion strategy we propose can theoretically generate 100% sequencing-ready templates. microRNA applications. Our study demonstrates the value of this novel technology for NGS library preparation, as well as the potential to facilitate many other molecular applications, such as radiolabeling the 3' end of RNA.

[0450] 4.2 3' branch ligation, a novel method for joining DNA ends

[0451] Typically, DNA ligation involves the joining of the 5'PO4 and 3'OH DNA ends of sticky or blunt-ended fragments. Sticky-end ligation is generally faster than blunt-end ligation and is less dependent on enzyme concentration. Bacteriophage T4 DNA ligase can catalyze both processes, while bacteriophage T4 DNA ligase uses ATP as a cofactor for energy production and requires Mg. 2+ . T4 DNA ligase has also been reported to ligate specific or degenerate single-stranded oligonucleotides to partially single-stranded substrates via hybridization. Here, we demonstrate unprecedented T4 DNA ligase-mediated ligation that does not require complementary base pairing and can ligate a blunt-ended DNA donor to the 3'OH terminus of a double-stranded DNA acceptor at a nick, a gap, or a 5'-overhang to form a branched structure (Figure 21a). Therefore, we use the term 3'-branched ligation (3'BL) to describe these ligations. The synthetic donor DNA we used contained a blunt-ended duplex end and an ssDNA end. The acceptor substrate contained one of the following structures: a dephosphorylated nick, a 1- or 8-nucleotide (nt) gap, or a 36-nt 5'-overhang. T4 ligase helps to connect the 5'PO4 of the adapter strand to the only ligatable 3'OH of the substrate strand, thereby forming a forked ligation product.

[0452] To optimize ligation efficiency, we extensively tested many factors that affect general ligation efficiency, including adapter::DNA substrate ratio, T4 ligase amount, final ATP concentration, Mg 2+ The concentration, pH, incubation time and presence of different additives such as polyethylene glycol-8000 (PEG-8000) and single-stranded binding protein (SSB) were investigated (Supplementary Figures 1 and 2). We found that the addition of PEG-8000 to a final concentration of 10% significantly increased the ligation efficiency from less than 10% to greater than 90% (Figure 21). The 3' branch ligation also worked over a wide range of ATP concentrations (from 1 uM to 1 mM) and Mg 2+ The ligase concentration (from 3mM to 10mM) of 3'BL is effective. The amount of ligase required for 3'BL is comparable to that of blunt-end ligation. Under our optimized conditions, we used a linker::substrate DNA molar ratio of 10 to 100 and reacted with 1mM ATP, 10mM MgCl2 and 10% PEG-8000 at pH 7.8 at 37°C for one hour. The same linker was ligated to a blunt-ended substrate and no ligase reaction was used as a control. To determine the yield of the ligation product, the reaction was performed on a denaturing polyacrylamide gel (Figure 21b) or a TBE gel (Figure 22a-d). The ratio of product to substrate intensity was used to quantify the ligation efficiency by ImageJ (Figure 21b). The 5'-overhang ligation (lane 11 in Figure 1b) appears to be more than 90% complete, even higher than the blunt-end ligation control (lane 14, 76.9%), indicating that the ligation efficiency of the DNA 5'-overhang is very high. The 1 or 8 nt vacant substrates (lanes 5 and 8) showed good ligation efficiency of about 60%. However, the nick ligation (lane 2) had the lowest efficiency, about 20%. However, if we incubate the nick ligation for 12 hours, the ligation yield can be improved, indicating that the kinetics of the nick ligation reaction are slower (Figure 22).

[0453] We also extended the study to different adapter and substrate sequences (Figure 22). The 5'PO4 end of three different adapters (Ad-T, Ad-A or Ad-GA) contained a T or A or dinucleotide GA before the common CTGCTGA sequence ligation. They were ligated to the 3'OH end of the acceptor template with T at the ligation junction. Overall, higher ligation efficiency (70-90%) than Ad-GA was observed with Ad-T and Ad-A in all cases except nick ligation (Figure 22), indicating some nucleotide preference of T4 DNA ligase at the ligation junction. Despite the adapter and substrate sequences, 5' overhang or 3' branch ligation always showed higher efficiency (60-90%), while the nick ligation efficiency was quite low after 1 hour of incubation. We assume that these differences in ligation efficiency are due to DNA bending, that is, the DNA bends at the beginning of the nick / gap / overhang and exposes the 3'OH group for ligation. Longer ssDNA regions may make the 3' end more accessible to ligation and lead to higher ligation efficiency. We also tested whether similar end-joining events might occur with 5' branched ligations. In contrast, no significant ligation of blunt-ended adapters at gaps or 3' overhangs to 5'PO4 ends was observed, suggesting that the tertiary structure requirements of T4 DNA ligase at the 5' end of the donor might be more stringent than at the 3' end.

[0454] 4.3.3' Branch ligation, a novel method for connecting DNA to RNA

[0455] We further investigated the 3'BL on DNA / RNA hybrids (ON21 / 22), which form 1 DNA and 1 RNA 5'-overhang (Figure 23a). Negative ligation controls included DNA / RNA hybrids, ssDNA or ssRNA oligonucleotides, incubated alone or with adapters (lanes 3, 4 and 5 in Figure 23b). Interestingly, when the DNA / RNA hybrids were incubated with adapters, we saw that the size of the RNA oligonucleotide changed from the original 29nt to 49nt with an efficiency of >90%, indicating that T4 DNA ligase can effectively connect the adapter to RNA. However, the DNA substrate remained unchanged (lanes 1 and 2 in Figure 23b). This indicates that the blunt-ended DNA adapter is connected to the 3' end of the RNA at the 5'-DNA overhang, but the 5'-RNA overhang is connected to the 3' end of the DNA. To confirm the required 5'-overhang structure for 3'BL, we performed the same ligation reaction, replacing the original DNA oligonucleotide (ON21) with another long DNA template (ON23) that is not complementary to ON22 RNA. Not surprisingly, no ligation was observed using ON23 DNA template, indicating that 3'BL can only occur with a 5'-overhang. Our findings suggest that T4 DNA ligase has certain substrate preferences, which may be caused by differences in protein-substrate binding affinity.

[0456] Previous studies have shown that T4 DNA ligase and T4 RNA ligase 2, but not T4 RNA ligase 1, can ligate 5'PO4 DNA ends to juxtaposed 3'OH DNA or RNA ends on RNA / DNA duplexes, but not RNA 3'OH (Bullard 2006, Biochem J 398: 135-144). We performed the same ligation test using T4 RNA ligase 1 and 2 ( Fig.24C ). It appears that R4 RNA ligases 1 and 2 can ligate blunt-ended adaptors to RNA, but the ligation efficiency is very low (<10%).

[0457] 4.4 Construction of directional transposon insertion library

[0458] Since 3' branch ligation has been shown to be useful for efficiently ligating adapters to several genomic structures, we explored its application in NGS workflows. Transposon-based library construction methods are time-saving and consume less input DNA than traditional NGS library preparation. However, using commercial transposon-based library preparation systems, only half of the labeled molecules are flanked by two different adapter sequences, while the labeled DNA is flanked by self-complementary regions that may form stable hairpin structures, which may compromise sequencing quality (Gorbacheva, 2015, Biotechniques Apr; 58(4): 200-202). In addition, PCR-mediated incorporation of adapter sequences is not suitable for whole-genome bisulfite sequencing or PCR-free NGS library construction.

[0459] To overcome these limitations, we developed a new protocol for transposon-based NGS library construction combined with 3'BL. Both Tn5 and MuA transposons work through a "cut and paste" mechanism, in which the transposon adapter sequence end is ligated to the 5' end of the target DNA, creating a 9 bp or 5 bp gap at the 3' end of the genomic DNA, respectively (Figure 24). Then, 3'BL is used to add another adapter sequence to the 3' end of the genomic DNA at the gap to complete the ligation of the directional adapter. We compared the efficiency of the 3'BL method with that of the dual transposon insertion method, which uses two different Tn5-based adapters, TnA and TnB. Human genomic DNA was incubated with either TnA transposome complexes alone or equimolar amounts of TnA and TnB transposome complexes. The TnA transposome fragmentation products were further used in 3'BL with the blunt-end adapter AdB, which shares a common adapter sequence with TnB. PCR amplification using two primers Pr-A and Pr-B designed for TnA and AdB / TnB adapters, respectively, showed similar PCR yields (Figure 24b, lanes 9 and 10, and Figure 24c), indicating that the two methods have the same efficiency. When only one primer specific for TnA or AdB / TnB adapter was used, no significant amplification was observed (Figure 24b and Figure 24c). As expected, due to PCR inhibition, both the 3'-ligation method and the dual transposon insertion method showed significantly higher PCR efficiency (Figure 24b, lanes 3 and 8, and Figure 24c) compared to the transposon insertion reaction using only TnA or TnB transposome complexes alone.

[0460] 4.5 Materials and methods

[0461] 3' branch ligation of double-stranded DNA

[0462] Substrates for 3'BL were prepared by mixing 2 pmol of ON1 or ON9 with 4 pmol of one or two additional oligonucleotides in pH 8 Tris-EDTA (TE) buffer (Life Technologies). Substrates 1 and 5 (nicks): ON1 / 2 / 3 and ON9 / 10 / 11; substrates 2 and 6 (1 bp gap): ON1 / 2 / 4 and ON9 / 10 / 12; substrate 3 (8 bp gap): ON1 / 4 / 5; substrates 4 and 9 (5' overhang): ON1 / 2 and ON9 / 10; substrate 7 (2 bp gap): ON9 / 10 / 13; substrate 8 (3 bp gap): ON9 / 10 / 14; blunt-end control: ON1 / 6 (Figure 1, Supplementary Table 1). The template was ligated to 180 pmol of adaptor (Ad-C: ON7 / 8, Ad-T: ON15 / 16, Ad-A: ON17 / 8 or Ad-GA: ON19 / 20) using 2400 units of T4 ligase (Enzymatics Inc) in 3'BL buffer [0.05 mg / ml BSA (New England Biolabs), 50 mM Tris-Cl pH 7.8 (Amresco), 10 mM MgCl2 (EMD Millipore), 0.5 mM DTT (VWR Scientific), 10% PEG-8000 (Sigma Aldrich) and 1 mM ATP (Sigma Aldrich)]. The ATP concentration was varied from 1 uM to 1 mM, MgCl2 was varied from 1 uM to 1 mM. 2+ Optimization tests were performed for concentration, pH, temperature from 12°C to 42°C, and additives such as PEG-8000 from 2.5% to 10% and SSB from 2.5 to 20 ng / ul. Ligation mixtures were prepared on ice and incubated at 37°C for 1 to 12 hours and then heat-inactivated at 65°C for 15 minutes. Samples were purified using Axygen beads (Corning) and eluted into 40 μL TE buffer. All ligation reactions were performed on 6% TBE or denaturing polyacrylamide gels (Life Technologies) and visualized on an Alpha Imager (Alpha Innotech). The amount of input control loaded was equal to or half the amount of template used for ligation. Ligation efficiency was estimated using ImageJ Software (NIH) by dividing the intensity of the ligated product by the total intensity of the ligated and unligated products.

[0463] 3' branch ligation of DNA / RNA hybrids

[0464] The substrate for 3'BL was composed of 10 pmol ON22 RNA oligonucleotide mixed with 2 pmol ON21 or ON23 DNA oligonucleotide. For T4 DNA ligase-mediated 3'BL, the substrate was incubated with Ad-T (ON15 / 16) in 3'BL buffer as described above and incubated at 37°C for 1 hour. 3'BL was performed using T4 RNA ligase 1 or 2 in their own 1x RNA ligase buffer (NEB) and 20% DMSO. All ligation products were determined on 6% denaturing polyacrylamide gels.

[0465] Construction of directional transposon insertion library

[0466] Transposon oligonucleotides used in this experiment were synthesized by Sangon Biotech. For two transposon experiments using TnA and TnB, TnA, TnB, and MErev oligonucleotides were annealed at a ratio of 1:1:2. For single transposon experiments using tn1, tn1 and MErev were annealed at a ratio of 1:1.

[0467] Transposome assembly was performed by mixing 15pmol of pre-annealed adapter, 7ul Tn5 transposase (Vazemy) and 5.5ul glycerol to obtain 20ul reactants, which were incubated at 30°C for 1 hour. Transposon insertion (Coriell19240) of genomic DNA was performed in 20ul reactants of transposomes comprising 100ng gDNA, TAG buffer (Vazyme) and 2ul assembled. The reaction was incubated at 55°C for 10 minutes, and then 100ul of PB buffer (Qiagen) was added to remove the transposome complex from the tagged DNA, and Agencourt AMPure XP beads (Beckman Coulter) were used for purification. Incubated at 25°C for 1 hour in a reactant containing 100pmol adapter, 600U of T4 DNA ligase (Enzymatics Inc.) and 3'BL buffer, AdB (ONB1, ONB2) was 3' branched and connected to the tag DNA. AMPure XP beads were used for purification of the reaction. PCR amplification of the labeled and gap-ligated DNA was performed in 50ul reactions containing 2ul of labeled or gap-ligated DNA, TAB buffer, 1ul of TruePrep Amplification Enzyme (Vazyme), 200mM dNTPs (Enzymatics Inc.), and 400mM each of primers Pr-A and Pr-B. The labeling reaction was performed as follows: 3 minutes at 72°C; 30 seconds at 98°C; 8 cycles of 10 seconds at 98°C, 30 seconds at 58°C, 2 minutes at 72°C; and 10 minutes extension at 72°C. The gap-ligation reaction was performed using the same procedure, but without the initial extension at 72°C for 3 minutes. The PCR reaction was purified using AMPure XP beads with single-step size selection or by dual fractionation. The purified product was quantified using the Qubit High Sensitivity DNA Kit (Invitrogen).

[0468] Example 4: 3' branch ligation: a novel method to ligate non-complementary DNA to the recessed or internal 3' OH terminus of DNA or RNA

[0469] Nucleic acid ligases are key enzymes that repair breaks in DNA or RNA during synthesis, repair, and recombination. Various molecular tools have been developed using the multiple activities of DNA / RNA ligases. However, other ligase activities remain to be discovered. In this paper, we demonstrated that T4 DNA ligase has an unconventional ability to ligate 5'-phosphorylated blunt-ended double-stranded DNA to DNA breaks at 3' recessed ends, gaps, or nicks to form 3' branched structures. Therefore, this base pairing-independent ligation is called 3' branched ligation (3'BL). In an extensive study of optimal ligation conditions similar to blunt-end ligation, the presence of 10% PEG-8000 in the ligation buffer significantly improved the ligation efficiency. Some nucleotide preferences were observed at the ligation site using different synthetic DNAs, indicating the level of ligation bias for 3'BL. In addition, we found that T4 DNA ligase can efficiently ligate DNA to the 3' end of RNA in DNA / RNA hybrids, while RNA ligase is less efficient in this reaction. These new properties of T4 DNA ligase can be used as a wide range of molecular techniques in many important applications. We performed a proof-of-concept study of a new directional tagging protocol for next-generation sequencing (NGS) library construction that eliminates reverse adapters and allows sample barcodes to be inserted adjacent to genomic DNA. Theoretically, 3'BL after single-transposon tagging can obtain 100% of usable templates, and our empirical data showed that the new approach produced higher yields compared to traditional double-transposon or Y-transposon tagging. We further explored the potential use of 3'BL in the preparation of targeted RNA NGS libraries with reduced structure-based bias and adapter dimer issues.

[0470] 5.1 Introduction

[0471] Ligases repair breaks in nucleic acids, and this activity is essential for cell viability and vitality. DNA ligases catalyze the formation of phosphodiester bonds between DNA ends and play a key role in DNA repair, recombination, and replication in vivo1-3. RNA ligases join 5'-phosphoryl (5'PO4) and 3'-hydroxyl (3'OH) RNA ends via phosphodiester bonds and are involved in RNA repair, splicing, and editing4. Ligases from all three kingdoms (bacteria, archaea, and eukaryotes) can be used as important molecular tools in vitro for applications such as cloning, ligase-based amplification or detection, and synthetic biology5-7.

[0472] One of the most widely used ligases in vitro is bacteriophage T4 DNA ligase, a single 55-kDA polypeptide that requires ATP as an energy source8. T4 DNA ligase typically ligates adjacent 5'PO4 and 3'OH termini of double-stranded DNA. In addition to sealing nicks and ligating sticky ends, T4 DNA ligase can also efficiently catalyze blunt-end ligation, which has not been observed for any other DNA ligase9,10. Some unusual catalytic properties of this ligase have been previously reported, such as sealing single-stranded gaps in double-stranded DNA, sealing nicks near abasic sites in double-stranded DNA (dsDNA), promoting the formation of intramolecular loops in partially double-stranded DNA, and ligating DNA strands containing 3' branch extensions11-13. Researchers have also observed template-independent ligations mediated by T4 ligase, such as sealing of mispaired nicks in dsDNA14 or even single-stranded DNA (ssDNA) ligations, although with very low efficiency15. These results suggest that perfect complementary base pairing at or near the ligation junction is not required for some of the unconventional T4 DNA ligase activities. T4 RNA ligases 1 and 2 are products of T4 bacteriophage genes 63 and 24, respectively. Both require adjacent 5'PO4 and 3'OH termini for successful ligation, with ATP hydrolyzed to AMP and PPi. Substrates for T4 RNA ligase 1 include single-stranded RNA and DNA, whereas T4 RNA ligase 2 preferentially seals nicks on dsRNA rather than ligating the ends of ssRNA16,17.

[0473] Here, we demonstrate an unconventional end-joining event mediated by T4 DNA ligase, which we term 3'-branched ligation (3'BL). This method can join DNA or DNA / RNA fragments at nicks, single-stranded gaps, or 3' recessed ends to form branched structures. This report includes an extensive study of various ligation cofactors and activators, as well as the optimization of ligation conditions for this novel ligation. Using our 3'BL protocol, no base pairing is required, and in most cases, including 1nt gaps, ligation can reach 70-90% completion rates. One application of this method is the ligation of adapters to DNA or RNA during NGS library preparation. Several genomic structures that were previously considered unligatable can now be used as substrates for 3'BL, resulting in high conversion rates of input DNA to adapter-ligated molecules while avoiding chimeras. We demonstrate that 3'BL can be combined with transposon tagmentation to increase library yields. The targeted tagging strategy we propose will theoretically generate 100% sequencing-ready templates. Our studies demonstrate the value of this novel technology for NGS library preparation and the potential to advance many other molecular applications.

[0474] 5.2 Results: 3' branch ligation, a novel method for joining DNA ends

[0475] Typically, DNA ligation involves the joining of the 5'PO4 and 3'OH DNA ends of sticky or blunt-ended fragments. Sticky-end ligation is generally faster and less dependent on enzyme concentration than blunt-end ligation. Both processes can be catalyzed by bacteriophage T4 DNA ligase, which uses ATP as a cofactor for energy production and requires Mg. 2+ 8. T4 DNA ligase has also been reported to ligate specific or degenerate single-stranded oligonucleotides to partially single-stranded substrates by hybridization18,19. Here, we demonstrate an unconventional T4 DNA ligase-mediated ligation that does not require complementary base pairing and can ligate a blunt-ended DNA donor to the 3′OH terminus of a double-stranded DNA acceptor at a 3′ recessed strand, a gap, or a nick (Figure 26a). Therefore, we use the term 3′-branched ligation (3′BL) to describe these ligations. The synthetic donor DNA we used contained a 5′ blunt double-stranded end and a 3′ ssDNA end. The acceptor substrate contained one of the following structures: a dephosphorylated nick, a 1- or 8-nucleotide (nt) gap, or a 3′ 36-nt recessed end (Supplementary Table 1). T4 ligase facilitates the ligation of the 5′PO4 of the donor strand to the uniquely ligatable 3′OH of the acceptor strand, forming a branched ligation product.

[0476] To optimize ligation efficiency, we extensively tested many factors that affect general ligation efficiency, including adapter:DNA substrate ratio, T4 ligase amount, final ATP concentration, Mg 2+ concentration, pH, incubation time and different additives such as polyethylene glycol-8000 (PEG-8000) and single-stranded binding protein (SSB). Adding PEG-8000 to a final concentration of 10% significantly increased the ligation efficiency from less than 10% to greater than 80% (Figures 26 and 27). A wide range of ATP concentrations (from 1 uM to 1 mM) and Mg 2+ The concentration (3mM to 10mM,) is compatible with 3'BL. The amount of ligase required for 3'BL is comparable to that for blunt-end ligation. Under our optimized conditions, we used a 30:100 donor:substrate DNA molar ratio and performed the reaction at 37°C, pH 7.8, with 1mM ATP, 10mM MgCl2, and 10% PEG-8000 for one hour. Ligation of the same donor with a blunt-end substrate and a reaction without ligase were used as positive and negative controls, respectively.

[0477] The ligation donor (Ad-G) is double-stranded at one end (5' phosphorylation and 3' dideoxy protection) and single-stranded at the other end (3' dideoxy protection) (Figure 26). The ligation substrate consists of the same bottom chain (ON1) and different top chains to form the nick, vacancy and protrusion structures. In order to quantify the yield of the ligation product, the reaction products were separated on a 6% denaturing polyacrylamide gel (Figure 26b). The ligation efficiency was calculated as the ratio of product to substrate intensity using ImageJ (Figure 26b). The 3'-recessed ligation (lane 11 in Figure 26b) appears to be about 90% complete, even higher than the flat end ligation control (lane 14, 72.74%), and indicates that the ligation efficiency with the 3'-recessed DNA end is very high. The 1-nt or 8-nt vacancy substrates (lanes 5 and 8) show a good ligation efficiency of about 45%. The nick ligation (lane 2) has the lowest efficiency, about 13%. However, the ligation yield improved when the nick ligation reaction was incubated for a longer time, indicating that the kinetics of the nick ligation reaction were slower.

[0478] We also extended the study to different adapter and substrate sequences (Fig. 27). The 5'PO4 termini of three different adapters (Ad-T, Ad-A or Ad-GA in Supplementary Table 1) contained a single T or A or dinucleotide GA at the ligation junction preceding the consensus CTGCTGA sequence. These 5'PO4 termini were ligated to the 3'OH termini of the acceptor template, respectively, with T added at the ligation junction. Overall, high ligation efficiencies (70-90%) were observed in most cases, except for nick ligation or 3'BL with Ad-GA (Fig. 27), thus indicating that T4 DNA ligase has some nucleotide preference at the ligation junction. Independent of the adapter and substrate sequences, 3' recessed ends or gap ligations consistently showed higher efficiencies (60-90%) in 1 h of incubation, while nick ligation efficiencies were very low. We hypothesize that these differences in ligation efficiencies are due to DNA bending, i.e., DNA bending at the beginning of the nick / gap / overhang and exposing the 3'OH group for ligation. The longer ssDNA region may make the 3' end more accessible for ligation and lead to higher ligation efficiency. We also tested whether similar end-joining events may occur with 5' branched ligations. In contrast to 3'BL, no obvious ligation of blunt-ended adapters to 5'PO4 ends was observed at the vacant or 5' recessed ends. This result suggests that the steric hindrance of T4 DNA ligase at the 5' end of the donor is greater than that at the 3' end.

[0479] 5.3: 3' branch ligation to connect DNA to RNA

[0480] We further studied the 3'BL on DNA / RNA hybrids (ON-21 / ON-23 in Table 3), which form one DNA and one RNA 5'-overhang (Figure 28a). Ligation on DNA / DNA hybrids was used as a positive control, while negative ligation controls included DNA / RNA hybrids, ssDNA or ssRNA oligonucleotides incubated alone or with adapters (lanes 3, 4 and 5 in Figure 28c). Interestingly, when DNA / RNA hybrids were incubated with blunt-ended dsDNA donors, we observed that the size of the RNA oligonucleotide changed from the original 29nt to 49nt after ligation. However, the DNA substrate remained unchanged (lanes 1 and 2 in Figure 28c). This result indicates that the blunt-ended dsDNA donor ligates to the 3'-end of the RNA at the 3'-recessed DNA end, but not to the 3'-end of the DNA at the 3'-recessed RNA end. As a positive control, DNA / DNA hybrids with 3' recessed ends on each side showed a band shift efficiency of nearly 100% toward the larger species on both strands. To confirm that 3' recessed structures are required for 3'BL, we performed the same ligation reaction while replacing the original DNA oligonucleotide (ON-21) with another long DNA template (ON-23) that is not complementary to ON-22 RNA (Figure 28b). Not surprisingly, no ligation was observed using the ON-23 DNA template (lanes 10-13 in Figure 28c). Our findings suggest that T4 DNA ligase can promote 3'BL on DNA / RNA hybrids and that this activity has a certain spatial substrate preference, which may be affected by differences in T4 DNA ligase binding affinity to the substrate.

[0481] Previous studies have reported that T4 DNA ligase and T4 RNA ligase 2, but not T4 RNA ligase 1, can efficiently ligate 5'PO4 DNA ends to juxtaposed 3'OH DNA or RNA ends in order to seal the nick in the DNA / RNA hybrid when the complementary strand is RNA rather than DNA17. Therefore, we performed the same ligation test using T4 RNA ligases 1 and 2 in 20% DMSO (Figure 28d) or 10% PEG. In both tests, T4 RNA ligase 1 and T4 RNA ligase 2 slightly ligated blunt-ended adapters to the 3' end of RNA in the DNA / RNA hybrid. Notably, in the RNA-only control, T4 RNA ligase 2 could ligate blunt-ended dsDNA adapters to ssRNA. In summary, T4 DNA ligase, but not T4 RNA ligase, is able to efficiently ligate blunt-ended dsDNA to the 3' end of RNA via 3'BL.

[0482] 5.4 Construction of directional marker library

[0483] Because 3'BL can be used to efficiently ligate adapters to multiple genomic structures, we explored its use in NGS workflows. Compared with traditional NGS library preparation, transposon-based library construction is fast and consumes less input DNA. However, using commercial transposon-based library preparation systems, only half of the labeled molecules are flanked by two different adapter sequences (Figure 29a), and the labeled DNA is flanked by self-complementary regions that may form stable hairpin structures and impair sequencing quality20. In addition, PCR-mediated incorporation of adapter sequences has not yet been adapted for whole-genome bisulfite sequencing or PCR-free NGS library construction.

[0484] To overcome these limitations, we developed a new protocol for transposon-based NGS libraries by incorporating 3'BL. Both Tn5 and MuA transposons work through a "cut and paste" mechanism in which the transposon adapter sequence end is ligated to the 5' end of the target DNA, creating a 9-bp or 5-bp gap at the 3'-end of the genomic DNA, respectively (Figure 29a). Subsequently, 3'BL can be used to add another adapter sequence to the 3' end of the genomic DNA at the gap to complete the ligation of the directional adapter (Figure 29c). We used the Tn5 transposon in this manuscript to compare the efficiency of the single-tag + 3'BL method (Figure 29c) with the efficiency of the dual-tag method (Figure 29a), which uses two different Tn5-based adapters TnA and TnB, as well as with the efficiency of another directional single-tag strategy using a Y adapter containing two different adapter sequences (Figure 29b). Human genomic DNA was incubated with equimolar amounts of TnA and TnB transposome complexes, or with TnA transposome complexes alone, or with TnY (TnA / B) transposome complexes.

[0485] Only the product of TnA transposome fragmentation is further used as the 3'BL template of the flat-end adapter AdB, which shares a common adapter sequence with TnB. PCR amplification is performed using two primers Pr-A and Pr-B, which are designed to recognize TnA and AdB / TnB adapters, respectively. Quantitative data show that TnA&AdB has the highest efficiency compared to TnA&TnB and TnY (TnA / B) (Figure 29d). When only one primer specific to the TnA adapter was used, no obvious amplification was observed (Figure 29d). As expected, due to PCR inhibition, the TnA-3'BL method, the dual label method, and the TnY method all showed significantly higher PCR efficiency than the label reaction using only TnA or TnB transposome complexes alone (Figure 29d).

[0486] We also sequenced these libraries using BGISEQ-500 and compared the base position deviations between the transposon interference end, the 3'BL end, and the conventional TA connection end (Figure 30). Obviously, the position deviation of the 3'BL end is less than that of the Tn5 end (Figure 30a-b), because the 3'BL end is affected by the transposon interruption and the 3'BL. Because only the first 6nt (positions 1-6) of the 3'BL end showed base deviation, and the deviation was similar but not identical to the hybrid Tn5 end (positions 30-35 after the 9-nt overhang), we believe that the position deviation observed at the 3'BL end is mainly caused by the Tn5 transposon. Therefore, the 3'BL causes the smallest deviation and is similar to the conventional TA connection (Figure 30c).

[0487] 5.5 Discussion

[0488] An important property of T4 DNA ligase is that it can efficiently ligate blunt-ended dsDNA21,22, which has not been observed with other DNA ligases. This ligase has also been reported to mediate some unconventional catalytic events, such as ligating single-stranded gaps or base mismatches in double-stranded DNA11,12, forming stem-loop molecules from partially double-stranded DNA13, or inefficiently ligating ssDNA in a template-independent manner20.

[0489] Here, we demonstrate that T4 DNA ligase catalyzes the ligation of blunt-ended dsDNA to the 3’OH termini of nicked dsDNA and to partially single-stranded double-stranded DNA with nicks or 5’-overhangs. In contrast, ligation to 5’PO4 ends in 5’ recessed ends or gaps was not observed, suggesting that after binding to the 5’PO4 ends of dsDNA adapters, T4 DNA ligase can access the 3’ recessed ends upon DNA bending. Using our 3’BL approach, no base pairing is required, and more than 70% completion rates can be achieved using optimized conditions even with a 1 nt gap. However, different ligation efficiencies were observed for ligating 5’T, A, or GA to 3’T (Figure 2), suggesting some sequence preference at the ligation junction. Despite the well-recognized ligation bias, 23 T4 DNA ligase is commonly used in the adapter addition step during NGS library preparation. The ability of T4 ligase to perform 3’BL has enabled adapter ligation to several genomic structures that were previously thought to be unligatable, thereby improving template utilization. 3’BL can also be combined with transposon tags. The traditional dual-transposon strategy only yields 50% of the labeled molecules suitable for the subsequent amplification step. However, when DNA labeling is performed using one transposon and subsequent 3'BL, an increased yield of molecules with different adapters at each insert end can be obtained (Figure 4). In addition, the labeled 3'BL product can be directly loaded into the Illumina flow cell as a PCR-free WGS library, which is difficult to do using the dual-transposon strategy.

[0490] Other directed transposon protocols have been proposed that use Y transposons composed of two different adapter sequences or replace the unligated strand of a single transposon with a second adapter oligonucleotide followed by gap filling and ligation24. However, these approaches continue to retain the inverted adapter sequence and are unable to insert sample barcodes near genomic DNA as in the tagged-3′BL protocol. According to NGS data, the genomic ends ligated with 3′BL also showed fewer positions with positional base composition bias, and the first 6-nt bias was mild, mainly caused by transposon interruption, indicating that 3′BL has minimal positional bias. Using this new library construction method, Wang et al. successfully achieved highly accurate and complete variant calling in WGS and phased variants nearly perfectly into long contigs with N50 sizes up to 23.4 Mb for long fragment reads (BioRxiv, https: / / doi.org / 10.1101 / 324392).

[0491] In this study, we also investigated 3'BL using templates that form chimeric DNA / RNA duplexes with 5'DNA and 5'RNA overhangs (Figure 3). Unexpectedly, blunt-ended dsDNA was efficiently ligated to the 3' end of RNA rather than DNA, indicating that T4 ligase has a preference for ternary complex formation. If T4 RNA ligase I or II was used to ligate the ends, the ligation efficiency was greatly reduced. Another preliminary but important application of 3'BL with T4 DNA ligase is the enrichment of mRNA or the construction of targeted RNA libraries, especially for miRNAs, which are small regulatory RNAs whose uncontrolled expression leads to many diseases25,26. Therefore, our 3'BL technology can be easily applied to the detection of cancer and Alzheimer's disease using miRNA. Hybridization with DNA probes targeting poly(A) tails or specific miRNA sequences can be used to create DNA-RNA hybrids with DNA 5'-overhangs, which are then ligated to adapter sequences with sample and / or UID barcodes by 3'BL. These universal sequences can then be reverse transcribed to generate cDNA targeting the RNA sequence. The use of 3'BL mediated by T4 DNA ligase may offer some advantages for NGS RNA library construction compared to current miRNA capture technologies. First, hybridization with the DNA strand will prevent the RNA strand from forming secondary structures, thus mitigating bias introduced by other experimental protocols. Second, T4 DNA ligase can achieve efficient adapter addition through 3'BL, thus avoiding intramolecular RNA interactions that RNA ligase can promote. Third, adapter dimers can be effectively eliminated, potentially making unnecessary gel purification unnecessary. This new approach can improve unbiased microRNA expression profiling with a simple and scalable workflow, and therefore, large-scale studies will become more affordable.

[0492] The findings of this study add to the growing understanding of T4 DNA ligase activity. We foresee that 3' branch ligation will become a versatile tool in molecular biology that will drive the development of new DNA engineering methods beyond the described NGS applications.

[0493] 5.6 Materials and methods

[0494] 3' branch ligation of double-stranded DNA

[0495] Substrates for 3'BL were composed of 2 pmol of ON1 or ON9 mixed with 4 pmol of one or two additional oligonucleotides in pH 8 Tris-EDTA (TE) buffer (Life Technologies), as follows: substrates 1 and 5 (nicks), ON-1 / 2 / 3 and ON-9 / 10 / 11; substrates 2 and 6 (1-nt gap), ON1 / 2 / 4 and ON9 / 10 / 12; substrate 3 (8-nt gap), ON1 / 4 / 5; substrates 4 and 9 (5' overhang), ON1 / 2 and ON9 / 10; substrate 7 (2-nt gap), ON9 / 10 / 13; substrate 8 (3-nt gap), ON9 / 10 / 14; blunt-end control, ON1 and ON6 (Figure 26, Supplementary Table 1). The template was ligated to 180 pmol of adaptor (Ad-G: ON7 / 8, Ad-T: ON15 / 16, Ad-A: ON17 / 8, or Ad-GA: ON19 / 20) using 2400 units of T4 ligase (Enzymatics Inc.) in 3'BL buffer [0.05 mg / ml BSA (New England Biolabs), 50 mM Tris-Cl pH 7.8 (Amresco), 10 mM MgCl2 (EMD Millipore), 0.5 mM DTT (VWR Scientific), 10% PEG-8000 (Sigma Aldrich), and 1 mM ATP (Sigma Aldrich)]. The ATP concentration was varied from 1 μM to 1 mM. 2+ Optimization tests were performed by changing the concentration from 3 mM to 10 mM, pH from 3 to 9, temperature from 12°C to 42°C, and adjusting additives (e.g., PEG-8000 from 2.5% to 10%, SSB from 2.5 to 20 ng / μL). The ligation mixture was prepared on ice and incubated at 37°C for 1 to 12 hours, followed by heat inactivation at 65°C for 15 minutes. The samples were purified using Axygen beads (Corning) and eluted into 40 μL TE buffer. All ligation reactions were performed on 6% TBE or denaturing polyacrylamide gels (Life Technologies) and visualized on an Alpha Imager (Alpha Innotech). The amount of input control loaded was equal to or half the amount of template used for ligation. Ligation efficiency was estimated using ImageJ Software (NIH) by dividing the intensity of the ligated product by the total intensity of the ligated and unligated products.

[0496] 3' branch ligation of DNA / RNA hybrids

[0497] The substrate for 3'BL was composed of 10 pmol ON-21 RNA oligonucleotide mixed with 2 pmol ON-21 or ON-23 DNA oligonucleotide. For T4 DNA ligase-mediated 3'BL, the substrate was incubated with Ad-T (ON15 / 16) in 3'BL buffer as described above and incubated at 37°C for 1 hour. 3'BL was performed using T4 RNA ligase 1 or 2 in 1x RNA ligase buffer (NEB) with 20% DMSO or 25% PEG. All ligation products were determined on 6% denaturing polyacrylamide gels.

[0498] Directed marker library construction

[0499] Transposon oligonucleotides used in this experiment were synthesized by Sangon Biotech. For two transposon experiments using TnA / TnB, oligonucleotides of TnA (ON24), TnB (ON25), and MErev (ON26) were annealed in a ratio of 1:1:2. For single transposon experiments using TnA, ON24 and ON26 were annealed in a ratio of 1:1. For Y (TnA & TnB) transposon experiments, ON24 and ON27 were annealed in a ratio of 1:1.

[0500] Transposon assembly was performed by mixing 100 pmol of pre-annealed adapters, 7 μL of Tn5 transposase, and enough glycerol to obtain a total of 20 μL of reactions and incubating them at 30°C for 1 hour. Labeling of genomic DNA (Coriell12878) was performed in a 20 μL reaction containing 100 ng gDNA, TAG buffer (homemade), and 1 μL of assembled transposon. The reaction was incubated at 55°C for 10 minutes; 40 μL of 6M guanidine hydrochloride (Sigma) was then added to remove the transposon complex from the labeled DNA, and the DNA was purified using Agencourt AMPure XP beads (Beckman Coulter). AdB (ON28 and ON29) was ligated to the labeled DNA at 25°C for 1 hour in a reaction containing 100 pmol adapters, 600 U of T4 DNA ligase (Enzymatics Inc.), and 3'BL buffer. The reaction was purified using AMPure XP beads. PCR amplification of labeled and gap-ligated DNA was performed in 50 μL reactions containing 2 μL of labeled or gap-ligated DNA, TAB buffer, 1 μL of TruePrep Amplification Enzyme (Vazyme), 200 mM dNTPs (Enzymatics Inc.), and 400 mM each of primers Pr-A and Pr-B. The labeling reaction was incubated as follows: 72°C for 3 minutes; 98°C for 30 seconds; 8 cycles of: 98°C for 10 seconds, 58°C for 30 seconds, and 72°C for 2 minutes; and 72°C for 10 minutes extension. The gap-ligated reaction was performed using the same procedure, but without the initial extension at 72°C for 3 minutes. PCR reactions using prA (ON30) or prA and prB (ON31) were purified using AMPure XP beads. The purified products were quantified using the Qubit High Sensitivity DNA Kit (Invitrogen).

[0501] 5.6 References for Example 4

[0502] 1.Lehnman,IRDNA ​​ligase:structure,mechanism,and function.Science(80-.).186,790-797(1974).

[0503] 2. Tomkinson, AE & Mackey, ZBStructure and function of mammalian DNAligases. Mutat. Res. Repair 407, 1-9 (1998).

[0504] 3.Timson,D.J.,Singleton,M.R.&Wigley,D.B.DNA ligases in therepair andreplication of DNA.Mutat.Res.Repair 460,301-318(2000).

[0505] 4.Ho,C.K.,Wang,L.K.,Lima,C.D.&Shuman,S.Structure andmechanism of RNAligase.Structure 12,327-339(2004).

[0506] 5.Tomkinson,A.E.,Vijayakumar,S.,Pascal,J.M.&Ellenberger,T.DNAligases:structure,reaction mechanism,and function.Chem.Rev.106,687-699(2006).

[0507] 6.Pascal,J.M.DNA and RNA ligases:structural variations and sharedmechanisms.Curr.Opin.Struct.Biol.18,96-105(2008).

[0508] 7.Shuman,S.DNA ligases:progress and prospects.J.Biol.Chem.284,17365-17369(2009).

[0509] 8.Dickson,K.S.,Burns,C.M.&Richardson,J.P.Determination of thefree-energy change for repair of a DNA phosphodiester bond.J.Biol.Chem.275,15828-15831(2000).

[0510] 9.Cai,L.,Hu,C.,Shen,S.,Wang,W.&Huang,W.Characterizationofbacteriophage T3 DNA ligase.J.Biochem.135,397-403(2004).

[0511] 10. Thermostable DNA Ligase.Available at:http: / / www.epibio.com / enzymes / ligases-kinases-phosphatases / dna-ligases / ampligase-thermostable-dna-ligase?details.

[0512] 11.Nilsson,S.V&Magnusson,G.Sealing of gaps in duplex DNA by T4DNAligase.Nucleic Acids Res.10,1425-1437(1982).

[0513] 12.Goffin,C.,Bailly,V.&Verly,W.G.Nicks 3′or 5′to AP sites ortomispaired bases,and one-nucleotide gaps can be sealed by T4 DNAligase.NucleicAcids Res.15,8755-8771(1987).

[0514] 13.Mendel-Hartvig,M.,Kumar,A.&Landegren,U.Ligase-mediatedconstructionof branched DNA strands:a novel DNA joining activity catalyzed byT4 DNAligase.Nucleic Acids Res.32,e2-e2(2004).

[0515] 14.Alexander,R.C.,Johnson,A.K.,Thorpe,J.A.,Gevedon,T.&Testa,S.M.Canonical nucleosides can be utilized by T4 DNA ligase asuniversaltemplate bases at ligation junctions.Nucleic Acids Res.31,3208-3216(2003).

[0516] 15.Kuhn,H.&Frank-Kamenetskii,M.D.Template-independentligation ofsingle-stranded DNA by T4 DNA ligase.FEBS J.272,5991-6000(2005).

[0517] 16.Ho,C.K.&Shuman,S.Bacteriophage T4 RNA ligase 2(gp24.1)exemplifiesa family of RNA ligases found in all phylogeneticdomains.Proc.Natl.Acad.Sci.99,12709-12714(2002).

[0518] 17.Bullard,D.R.&Bowater,R.P.Direct comparison of nick-joiningactivityof the nucleic acid ligases from bacteriophage T4.Biochem.J.398,135-144(2006).

[0519] 18.Broude,N.E.,Sano,T.,Smith,C.L.&Cantor,C.R.Enhanced DNAsequencingby hybridization.Proc.Natl.Acad.Sci.91,3072-3076(1994).

[0520] 19.Gunderson,K.L.et al.Mutation detection by ligation to complete n-mer DNA arrays.Genome Res.8,1142-1153(1998).

[0521] 20.Gorbacheva,T.,Quispe-Tintaya,W.,Popov,V.N.,Vijg,J.&Maslov,A.Y.Improved transposon-based library preparation for the Ion Torrentplatform.Biotechniques 58,200(2015).

[0522] 21.Sgaramella,V.&Khorana,H.G.CXII.Total synthesis of thestructuralgene for an alanine transfer RNA from yeast.Enzymic joining of thechemicallysynthesized polydeoxynucleotides to form the DNA duplexrepresentingnucleotide sequence 1 to 20.J.Mol.Biol.72,427-444(1972).

[0523] 22.SGARAMELLA,V.&EHRLICH,S.D.Use of the T4 PolynucleotideLigase inThe Joining of Flush-Ended DNA Segments Generated byRestrictionEndonucleases.FEBS J.86,531-537(1978).

[0524] 23.Seguin-Orlando,A.et al.Ligation bias in illumina next-generationDNA libraries:implications for sequencing ancient genomes.PLoS One8,e78575(2013).

[0525] 24.Goryshin,I.,Baas,B.,Vaidyanathan,R.&Maffitt,M.Oligonucleotidereplacement for di-tagged and directional libraries.(2016).

[0526] 25.Bushati,N.&Cohen,S.M.microRNA functions.Annu.Rev.Cell Dev.Biol.23,175-205(2007).

[0527] 26. Mallory, AC & Vaucheret, H. Functions of microRNAs and related smallRNAs in plants. Nat. Genet. 38, S31 (2006).

[0528] While the present invention has been disclosed with reference to certain aspects and embodiments, it will be apparent that other embodiments and variations of the invention may be devised by others skilled in the art without departing from the true spirit and scope of the invention.

[0529] Each publication and patent document cited in this disclosure is incorporated herein by reference for all purposes in the United States of America as if each such publication or document was specifically and individually indicated to be incorporated herein by reference. Citation of publications and patent documents is not intended to indicate that any such document is relevant prior art nor constitutes an admission as to its contents or date.

[0530]

[0531] Table 1B

[0532]

[0533] Table 2. Scaffolding statistics

[0534]

[0535] 1 HiC read pairs from human embryonic stem cells (hESCs) (30) were downloaded and used to scaffold SMRT reads using SALSA (28) and the same process as the stLFR library.

[0536] 2 Results reported by Ghurye et al. using the same HiC read pairs to scaffold SMRT reads using SALSA.

[0537] Table 3

[0538]

[0539] Table 4

[0540]

[0541] Table 5

[0542]

[0543]

[0544]

[0545] Sequence Listing <110> Shenzhen MGI Technology Co., Ltd. Shenzhen BGI Life Sciences Institute <120> Single-tube bead-based DNA co-barcoding for sequencing, haplotype typing, and assembly <130> 092171-1130724-5068-WO <140> PCT / US2019 / 031161 <141> 2019-05-07 <150> 62 / 687,159 <151> 2018-06-19 <150> 62 / 672,501 <151> 2018-05-16 <150> 62 / 668,757 <151> 2018-05-08 <160> 31 <170> PatentIn version 3.5 <210> 1 <211> 28 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 1 aaaaaaaacc cgtagccatg tcgttctg 28 <210> 2 <211> 32 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (7)..(26) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (7)..(26) <223> This region may include 4-20 nucleotides <400> 2 ctctggnnnn nnnnnnnnnn nnnnnntctg cg 32 <210> 3 <211> 26 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 3 ccagagcaga acgacatggc tacggg 26 <210> 4 <211> 60 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (35)..(54) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (35)..(54) <223> This region may include 4-20 nucleotides <400> 4 aaaaaaaacc cgtagccatg tcgttctgct ctggnnnnnn nnnnnnnnnn nnnntctgcg 60 <210> 5 <211> 46 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (1)..(20) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (1)..(20) <223> This region may include 4-20 nucleotides <400> 5 nnnnnnnnnn nnnnnnnnnn ccagagcaga acgacatggc tacggg 46 <210> 6 <211> 60 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (35)..(54) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (35)..(54) <223> This region may include 4-20 nucleotides <400> 6 aaaaaaaacc cgtagccatg tcgttctgct ctggnnnnnn nnnnnnnnnn nnnntctgcg 60 <210> 7 <211> 46 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (1)..(20) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (1)..(20) <223> This region may include 4-20 nucleotides <400> 7 nnnnnnnnnn nnnnnnnnnn ccagagcaga acgacatggc tacggg 46 <210> 8 <211> 80 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (35)..(54) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (35)..(54) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (61)..(80) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (61)..(80) <223> This region may include 4-20 nucleotides <400> 8 aaaaaaaacc cgtagccatg tcgttctgct ctggnnnnnn nnnnnnnnnn nnnntctgcg 60 nnnnnnnnnn nnnnnnnnnn 80 <210> 9 <211> 78 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (7)..(26) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (7)..(26) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (33)..(52) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (33)..(52) <223> This region may include 4-20 nucleotides <400> 9 ggaaggnnnn nnnnnnnnnn nnnnnncgca gannnnnnnn nnnnnnnnnn nnccagagca 60 gaacgacatg gctacggg 78 <210> 10 <211> 121 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (35)..(54) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (35)..(54) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (61)..(80) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (61)..(80) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (87)..(106) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (87)..(106) <223> This region may include 4-20 nucleotides <400> 10 aaaaaaaacc cgtagccatg tcgttctgct ctggnnnnnn nnnnnnnnnn nnnntctgcg 60 nnnnnnnnnn nnnnnnnnnn ccttccnnnn nnnnnnnnnn nnnnnntagc atggactatg 120 g 121 <210> 11 <211> 113 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (16)..(35) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (16)..(35) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (42)..(61) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (42)..(61) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (68)..(87) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (68)..(87) <223> This region may include 4-20 nucleotides <400> 11 ccatagtcca tgctannnnn nnnnnnnnnn nnnnnggaag gnnnnnnnnn nnnnnnnnnn 60 ncgcagannn nnnnnnnnnn nnnnnnncca gagcagaacg acatggctac ggg 113 <210> 12 <211> 121 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (35)..(54) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (35)..(54) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (61)..(80) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (61)..(80) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (87)..(106) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (87)..(106) <223> This region may include 4-20 nucleotides <400> 12 aaaaaaaacc cgtagccatg tcgttctgct ctggnnnnnn nnnnnnnnnn nnnntctgcg 60 nnnnnnnnnn nnnnnnnnnn ccttccnnnn nnnnnnnnnn nnnnnntagc atggactatg 120 g 121 <210> 13 <211> 121 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <221> modified_base <222> (35)..(54) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (35)..(54) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (61)..(80) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (61)..(80) <223> This region may include 4-20 nucleotides <220> <221> modified_base <222> (87)..(106) <223> a, c, t, g, unknown or other <220> <221> misc_feature <222> (87)..(106) <223> This region may include 4-20 nucleotides <400> 13 aaaaaaaacc cgtagccatg tcgttctgct ctggnnnnnn nnnnnnnnnn nnnntctgcg 60 nnnnnnnnnn nnnnnnnnnn ccttccnnnn nnnnnnnnnn nnnnnntagc atggactatg 120 g 121 <210> 14 <211> 27 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 14 aaaaaaaaaa tgtgagccaa ggagttg 27 <210> 15 <211> twenty three <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 15 ccagagcaac tccttggctc aca 23 <210> 16 <211> 43 <212> RNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 16 gcacugacga caugaucacc aaggaucgcc auaguccaug cua 43 <210> 17 <211> 58 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 17 cgatccttgg tgatcatgtc gtcagtgctt gtcttcctaa gatgtgtata agagacag 58 <210> 18 <211> 38 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 18 gcctccctcg cgccatcaga gatgtgtata agagacag 38 <210> 19 <211> 19 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <220> <223> Description of the combined DNA / RNA molecules: Synthetic oligonucleotides <400> 19 ctguctcuta tacacauct 19 <210> 20 <211> 17 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 20 tgtgagccaa ggagttg 17 <210> twenty one <211> 19 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> twenty one gcctccctcg cgccatcag 19 <210> twenty two <211> 38 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: Synthetic primers <400> twenty two gcctccctcg cgccatcaga gatgtgtata agagacag 38 <210> twenty three <211> 34 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: Synthetic primers <400> twenty three cgagaacgtc ttgtgagcca aggagttgct ctgg 34 <210> twenty four <211> 39 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: Synthetic primers <400> twenty four cgtcagtgct tgtcttccta agatgtgtat aagagacag 39 <210> 25 <211> 31 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: Synthetic primers <400> 25 tgatcaccaa ggatcgccat agtccatgct a 31 <210> 26 <211> 40 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: Synthetic primers <400> 26 ctgtctctta tacacatctt aggaagacaa gcactgacga 40 <210> 27 <211> 19 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 27 ctgatggcgc gagggaggc 19 <210> 28 <211> 12 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 28 tcgcgccatc ag 12 <210> 29 <211> 36 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: Synthetic primers <400> 29 caactccttg gctcacacgg agggagcgcg gtagtc 36 <210> 30 <211> 43 <212> RNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 30 gcacugacga caugaucacc aaggaucgcc auaguccaug cua 43 <210> 31 <211> 13 <212> DNA <213> Artificial sequence <220> <223> Description of artificial sequences: synthetic oligonucleotides <400> 31 cccccccccccc 13

Claims

1. A method for preparing a sequencing library for target nucleic acid sequencing without using nanodroplets, the method comprising: (a) transposing an insert sequence into a first fragment of the target nucleic acid in solution by a transposase, wherein the insert sequence comprises a hybridizing sequence, and wherein the transposition generates a nick in the first fragment; (b) combining (i) the first fragment of the target nucleic acid from (a), (ii) the splint oligonucleotide, and (iii) a population of beads into a single mixture or in a single container, wherein each bead comprises a capture oligonucleotide immobilized thereon, the capture oligonucleotide comprising 1) a barcode-containing sequence, wherein the oligonucleotides immobilized on the same single bead comprise the same barcode-containing sequence, while the majority of beads have different barcode-containing sequences, 2) a common sequence that is complementary to at least a portion of a splint oligonucleotide, wherein another portion of the splint oligonucleotide is complementary to at least a portion of the hybridization sequence; 3) first PCR primer annealing site; (c) Ligating the capture oligonucleotide of a single bead to the intervening hybridization sequence of a single first fragment.

2. The method of claim 1, wherein the capture oligonucleotide comprises the first PCR primer annealing site, further comprising connecting a 3' branched linker oligonucleotide to the first fragment at the cut, wherein the connection of the linker oligonucleotide is a 3' branched connection, and wherein the linker oligonucleotide comprises a second PCR primer annealing site.

3. The method of claim 2, wherein the 3' branched ligation adapter oligonucleotide is a blunt-ended adapter, and the 3' branched ligation comprises covalently linking a 5' phosphate from the blunt-ended adapter to a recessed 3' hydroxyl at the nick of the first fragment.

4. The method of claim 2, wherein the first PCR primer annealing site and the second PCR primer annealing site have different sequences.

5. The method of claim 3, wherein the 3' branch ligation adapter oligonucleotide comprises a barcode sequence. The method of claim 1 , wherein in the step of transposing the insert sequence into the first fragment, the transposase remains bound to the first fragment.

7. The method of any one of claims 1-6, comprising removing the transposase from the first fragment, thereby generating subfragments.

8. The method of claim 7, comprising amplifying the subfragments to produce amplicons. 9 . The method of claim 8 , comprising sequencing the amplicons to generate sequence reads, wherein sequence reads having the same barcode sequence are derived from the same first fragment.

10. The method according to any one of claims 1 to 6, wherein the target nucleic acid is genomic DNA.

11. The method of claim 1, wherein two different transposons are inserted into the first fragment.

12. The method of claim 10, wherein the genomic DNA is from a eukaryotic organism.

13. The method of claim 10, wherein the genomic DNA is from a single individual organism.

14. The method according to any one of claims 1 to 6, wherein the target nucleic acid is from a prokaryotic organism.

15. The method according to any one of claims 1 to 6, wherein the target nucleic acid is not amplified prior to step (a).

16. The method of any one of claims 1 to 6, wherein the bead comprises at least 100,000 copies of the capture oligonucleotide.

17. The method of claim 16, wherein more than 10 million beads are used in the reaction.

18. The method of claim 8, wherein prior to said amplifying, (i) at least some of the capture oligonucleotides are enzymatically removed and / or at least some of the splint oligonucleotides are enzymatically removed and / or at least some of the chimeric end (ME) sequences are enzymatically removed.

19. The method according to any one of claims 1 to 6, wherein: At least some of the capture oligonucleotides are enzymatically removed following step (c) using one or more exonucleases.

20. The method of any one of claims 1-6, wherein the capture oligonucleotide and / or splint oligonucleotide comprises uracil.

21. The method of claim 20, wherein the uracil-containing oligonucleotide is degraded by treatment with uracil-DNA glycosylase (UDG).

22. The method of claim 1, wherein the capture oligonucleotide is attached to a bead at the 5' end.

23. The method of any one of claims 1-6, wherein >50% of the first DNA fragments are barcoded with a unique barcode.

24. The method of any one of claims 1-6, wherein >50% of the subfragments in the fragments are linked to barcode oligonucleotides.

25. The method of any one of claims 1-6, wherein the spacing of the barcoded capture oligonucleotides on the beads is <50 nm on average.

26. The method of claim 7, comprising sequencing at least some of the subfragments to generate sequence reads.

27. The method of claim 26, wherein sequence reads having the same barcode sequence are from the same first fragment.

28. The method of claim 27, wherein >10% of the subfragments of the long fragments of the subfragments are sequenced.

29. The method according to claim 9 or 26, comprising the following additional steps: (e) assigning a majority of the sequence reads to the corresponding first fragment; and (f) assembling the sequence reads to generate an assembled sequence of the target.

30. The method of any one of claims 1-6, wherein the majority of the first fragments in step (a) are longer than 20 kb.

31. The method of any one of claims 1-6, wherein the majority of the first fragments in step (a) are longer than 50 kb.

32. The method of any one of claims 1-6, wherein the majority of the first fragments in step (a) are longer than 100 kb.

33. The method of any one of claims 1-6, wherein the majority of the first fragments are 50 kb to 200 kb in length.

34. The method of any one of claims 1-6, wherein a plurality of said first fragments are from a single cell.

35. The method of claim 10, wherein the plurality of genomic DNA fragments are from 1-100 eukaryotic cells.

36. The method of claim 19, wherein the one or more exonucleases comprise exonuclease I or III or both.

37. The method of claim 35, wherein the plurality of genomic DNA fragments are from 3-30 human cells.

38. The method of any one of claims 1-6, wherein the single container or the mixture comprises 5-100 genome equivalents of human DNA.

39. The method of any one of claims 1-6, wherein the single container or the mixture comprises 5-20 genome equivalents of human DNA.

40. The method of any one of claims 1-6, wherein the single container or the mixture comprises 50-1000 genome equivalents of human DNA.

41. The method of claim 1, wherein: The bead population includes at least 100,000 different tag sequences in total.

42. The method according to claim 41, wherein: The bead population includes at least one million different tag sequences in total.

43. The method of claim 7, wherein the subfragments are 200 to 2000 bp in size.

44. The method of claim 8, comprising sequencing at least some of the amplicons to generate sequence reads.

45. The method of any one of claims 1-6, wherein the bead comprises at least 100,000 copies of the capture oligonucleotide.

46. ​​The method of any one of claims 1-6, wherein the plurality of first fragments are from a single cell.

47. The method of any one of claims 1-6, wherein the single container or mixture comprises 5-100 genome equivalents of human DNA.

48. The method of any one of claims 1 to 6, wherein the target nucleic acid is from a heterogeneous population of prokaryotes.

49. The method of claim 5, wherein the barcode sequence is a sample barcode sequence.

Images