Contiguity Preserving Transposition
Embodiments provided herein relate to methods and compositions for preparing an immobilized library of barcoded DNA fragments of a target nucleic acid, identifying genomic variants, determining the contiguity information, phasing information, and methylation status of the target nucleic acid.
This application is a continuation of Ser. No. 15/519,482, filed on Apr. 14, 2017, which is a national stage entry of International Patent Application No. PCT/US2015/056040, filed Oct. 16, 2015, which claims priority to U.S. Provisional Patent Application No. 62/065,544 filed on Oct. 17, 2014, U.S. Provisional Patent Application No. 62/157,396 filed on May 5, 2015, and U.S. Provisional Patent Application No. 62/242,880 filed on Oct. 16, 2015, which are hereby incorporated by reference in their entirety.
SEQUENCE LISTINGThe instant application contains a Sequence Listing which has been submitted electronically in ASCII format and is hereby incorporated by reference in its entirety. Said ASCII copy, created on Oct. 26, 2018, is named SequenceListing.txt and is 2,492 bytes in size.
FIELD OF THE INVENTIONEmbodiments of the present invention relate to sequencing nucleic acids. In particular, embodiments of the methods and compositions provided herein relate to preparing nucleic acid templates and obtaining sequence data therefrom.
BACKGROUND OF THE INVENTIONThe detection of specific nucleic acid sequences present in a biological sample has been used, for example, as a method for identifying and classifying microorganisms, diagnosing infectious diseases, detecting and characterizing genetic abnormalities, identifying genetic changes associated with cancer, studying genetic susceptibility to disease, and measuring response to various types of treatment. A common technique for detecting specific nucleic acid sequences in a biological sample is nucleic acid sequencing.
Nucleic acid sequencing methodology has evolved significantly from the chemical degradation methods used by Maxam and Gilbert and the strand elongation methods used by Sanger. Today several sequencing methodologies are in use which allow for the parallel processing of nucleic acids all in a single sequencing run. As such, the information generated from a single sequencing run can be enormous.
SUMMARY OF THE INVENTIONIn one aspect, described herein are methods of preparing a library of barcoded DNA fragments of a target nucleic acid. The methods include contacting a target nucleic acid with a plurality of transposome complexes, each transposome complex includes: transposons and transposases, in which the transposons comprise transferred strands and non-transferred strands. At least one of the transposons of the transposome complex comprises an adaptor sequence capable of hybridizing to a complementary capture sequence. The target nucleic acid is fragmented into a plurality of fragments and inserting plurality of transferred strands to the 5′ end of at least one strand of the fragments while maintaining the contiguity of the target nucleic acid. The plurality of fragments of the target nucleic acid are contacted with a plurality of solid supports, each of the solid supports in the plurality comprising a plurality of immobilized oligonucleotides, each of the oligonucleotides comprising a complementary capture sequence and a first barcode sequence, and wherein the first barcode sequence from each solid support in the plurality of the solid supports differs from the first barcode sequence from other solid supports in the plurality of solid supports. The barcode sequence information is transferred to the target nucleic acid fragments, thereby producing an immobilized library of double-stranded fragments wherein at least one strand is 5′-tagged with the first barcode such that at least two fragments of the same target nucleic acid receives identical barcode information.
In one aspect, described herein are methods for determining contiguity information of a target nucleic acid sequence The methods include contacting the target nucleic acid with a plurality of transposome complexes, each transposome complex comprising: transposons and transposases, in which the transposons comprise transferred strands and non-transferred strands, in which at least one of the transposons of the transposome complex comprise an adaptor sequence capable of hybridizing to a complementary capture sequence. The target nucleic acid is fragmented into a plurality of fragments and plurality of transferred strands is inserted into the plurality of fragments while maintaining the contiguity of the target nucleic acid. The plurality of fragments of the target nucleic acid is contacted with a plurality of solid supports. Each of the solid supports in the plurality comprising a plurality of immobilized oligonucleotides, each of the oligonucleotides comprising a complementary capture sequence and a first barcode sequence, and wherein the first barcode sequence from each solid support in the plurality of the solid supports differs from the first barcode sequence from other solid supports in the plurality of solid supports. The barcode sequence information is transferred to the target nucleic acid fragments such that at least two fragments of the same target nucleic acid receive identical barcode information. The sequence of the target nucleic acid fragments and the barcode sequences are determined. The contiguity information of the target nucleic acid are determined by identifying the barcode sequences. In some embodiments, the transposases of transposome complexes are removed after transposition and subsequent hybridization of the adaptor sequences of the transposon to the complimentary capture sequence. In some embodiments, the transposases are removed by SDS treatment. In some embodiments, the transposases are removed by proteinase treatment.
In one aspect, described herein are methods for simultaneously determining phasing information and methylation status of a target nucleic acid sequence. The methods include contacting the target nucleic acid with a plurality of transposome complexes, each transposome complex includes transposons and transposases, in which the transposons comprise transferred strands and non-transferred strands, wherein at least one of the transposons of the transposome complex comprise an adaptor sequence capable of hybridizing to a complementary capture sequence. The target nucleic acid is fragmented into a plurality of fragments and plurality of transferred strands is inserted into the target nucleic acid fragments while maintaining the contiguity of the target nucleic acid. The plurality of fragments of the target nucleic acid are contacted with a plurality of solid supports, each of the solid supports in the plurality comprising a plurality of immobilized oligonucleotides, each of the oligonucleotides comprising a complementary capture sequence and a first barcode sequence, and wherein the first barcode sequence from each solid support in the plurality of the solid supports differs from the first barcode sequence from other solid supports in the plurality of solid supports. The barcode sequence information is transferred to the target nucleic acid fragments such that at least two fragments of the same target nucleic acid receive identical barcode information. The target nucleic acid fragments comprising barcodes are subjected to bisulfite treatment, thereby generating bisulfite treated target nucleic acid fragments comprising barcodes. The sequence of the bisulfite treated target nucleic acid fragments and the barcode sequences are determined. The contiguity information of the target nucleic acid is determined by identifying the barcode sequences.
In one aspect, described herein are methods of preparing an immobilized library of tagged DNA fragments. The methods include providing a plurality of solid supports having transposome complexes immobilized thereon, in which the transposome complexes are multimeric and the transposome monomeric units of the same transposome complex are linked to each other, and wherein said transposome monomeric units comprise a transposase bound to a first polynucleotide, said first polynucleotide comprising (i) a 3′ portion comprising a transposon end sequence, and (ii) a first adaptor comprising a first barcode. A target DNA is applied to the plurality of solid supports under conditions whereby the target DNA is fragmented by the transposome complexes, and the 3′ transposon end sequence of the first polynucleotide is transferred to a 5′ end of at least one strand of the fragments; thereby producing an immobilized library of double-stranded fragments wherein at least one strand is 5′-tagged with the first barcode.
In one aspect, described herein are methods of preparing a sequencing library for determining the methylation status of a target nucleic acid. The methods include fragmenting the target nucleic acid into two or more fragments. A first common adaptor sequence is incorporated into the 5′-end of the fragments of the target nucleic acid, wherein the adaptor sequence comprises a first primer binding sequence and an affinity moiety, wherein the affinity moiety in one member of the binding pair. The target nucleic acid fragments are denatured. The target nucleic acid fragments are immobilized on a solid support, in which the solid support comprises other member of the binding pair and the immobilization of the target nucleic acid is by binding of the binding pair. The immobilized target nucleic acid fragments are subjected to bisulfite treatment. A second common adaptor sequence is incorporated to the bisulfite treated immobilized target nucleic acid fragments, wherein the second common adaptor comprises a second primer binding site. The bisulfite treated target nucleic acid fragments immobilized on solid support is amplified thereby producing a sequencing library for determining the methylation status of a target nucleic acid.
In one aspect, described herein are methods of preparing a sequencing library for determining the methylation status of a target nucleic acid. The methods include providing a plurality of solid support comprising immobilized transposome complexes immobilized thereon. The transposome complexes comprise transposons and transposases, in which the transposons comprise transferred strands and non-transferred strands. The transferred strand comprises (i) a first portion at the 3′-end comprising the transposase recognition sequence, and (ii) a second portion located 5′ to the first portion comprising a first adaptor sequence and first member of a binding pair. The first member of the binding pair binds to a second member of the binding pair on the solid support, thereby immobilizes the transposon to the solid support. The first adaptor also comprises a first primer binding sequence. The non-transferred strand comprises (i) a first portion at the 5′-end comprising the transposase recognition sequence and (ii) a second portion located 3′ to the first portion comprising a second adaptor sequence, in which the terminal nucleotide at the 3′-end is blocked. The second adaptor also comprises a second primer binding sequence The target nucleic acid is contacted with the plurality of solid support comprising immobilized transposome complexes. The target nucleic acid is fragmented into a plurality of fragments and plurality of transferred strands are inserted to the 5′ end of at least one strand of the fragments, thereby immobilizing the target nucleic acid fragments to the solid support. The 3′-end of the fragmented target nucleic acid is extended with a DNA polymerase. The non-transferred strand is ligated to the 3′-end of the fragmented target nucleic acid. The immobilized target nucleic acid fragments are subjected to bisulfite treatment. The 3′-end of the immobilized target nucleic acid fragments damaged during the bisulfite treatment is extended by using a DNA polymerase such that the 3′-end of the immobilized target nucleic acid fragments comprise a homopolymeric tail. A second adaptor sequence is introduced to the 3′-end of the immobilized target nucleic acid fragments damaged during the bisulfite treatment. The bisulfite treated target nucleic acid fragments immobilized on solid support are amplified using a first and a second primer, thereby producing a sequencing library for determining the methylation status of a target nucleic acid.
In one aspect, disclosed herein are methods of preparing a sequencing library for determining the methylation status of a target nucleic acid. The methods include a. contacting the target nucleic acid with transposome complexes, in which the transposome complexes comprise transposons and transposases. The transposons comprise transferred strands and non-transferred strands. The transferred strand includes (i) a first portion at the 3′-end comprising the transposase recognition sequence, and (ii) a second portion located 5′ to the first portion comprising a first adaptor sequence and first member of a binding pair, wherein the first member of the binding pair binds to a second member of the binding pair. The non-transferred strand includes (i) a first portion at the 5′-end comprising the transposase recognition sequence and (ii) a second portion located 3′ to the first portion comprising a second adaptor sequence, in which the terminal nucleotide at the 3′-end is blocked, and wherein the second adaptor comprises a second primer binding sequence. The target nucleic acid is fragmented into a plurality of fragments and inserting plurality of transferred strands to the 5′ end of at least one strand of the fragments, thereby immobilizing the target nucleic acid fragments to the solid support. The target nucleic acid fragments comprising the transposon end are contacted with the plurality of solid support comprising second member of the binding pair, wherein binding of the first member of the binding pair to the second member of the binding pair immobilizes the target nucleic acid to the solid support. The 3′-end of the fragmented target nucleic acid is extended with a DNA polymerase. The non-transferred strand is ligated to the 3′-end of the fragmented target nucleic acid. The immobilized target nucleic acid fragments are subjected to bisulfite treatment. The 3′-end of the immobilized target nucleic acid fragments damaged during the bisulfite treatment is extended by using a DNA polymerase such that the 3′-end of the immobilized target nucleic acid fragments comprise a homopolymeric tail. A second adaptor sequence is introduced to the 3′-end of the immobilized target nucleic acid fragments damaged during the bisulfite treatment. The bisulfite treated target nucleic acid fragments immobilized on solid support are amplified using a first and a second primer, thereby producing a sequencing library for determining the methylation status of a target nucleic acid.
In some embodiments, the terminal nucleotide at the 3′-end of the second adaptor is blocked by a member selected from the group consisting of a dideoxy nucleotide, a phosphate group, thiophosphate group, and an azido group.
In some embodiments, affinity moieties can be members of a binding pair. In some cases, the modified nucleic acids may comprise a first member of a binding pair and the capture probe may comprise a second member of the binding pair. In some cases, capture probes may be immobilized to a solid surface and the modified nucleic acid may comprise a first member of a binding pair and the capture probe may comprise a second member of the binding pair. In such cases, binding the first and second members of the binding pair immobilizes the modified nucleic acid to the solid surface. Examples of binding pair include, but are not limited to biotin-avidin, biotin-streptavidin, biotin-neutravidin, ligand-receptor, hormone-receptor, lectin-glycoprotein, oligonucleotide-complementary oligonucleotide, and antigen-antibody.
In some embodiments, the first common adaptor sequence is incorporated to the 5′-end fragments of the target nucleic acid by one-sided transposition. In some embodiments, the first common adaptor sequence is incorporated to the 5′-end fragments of the target nucleic acid by ligation. In some embodiments, incorporating the second common adaptor sequence into the bisulfite treated immobilized target nucleic acid fragments includes (i) extending the 3′-end of the immobilized target nucleic acid fragments using terminal transferase to comprise a homopolymeric tail; (ii) hybridizing an oligonucleotide comprising a single stranded homopolymeric portion and a double stranded portion comprising the second common adaptor sequence, wherein the ingle stranded homopolymeric portion is complementary to the homopolymeric tail; and (iii) ligating the second common adaptor sequence to the immobilized target nucleic acid fragments, thereby incorporating the second common adaptor sequence into the bisulfite treated immobilized target nucleic acid fragments.
In some embodiments, the target nucleic acid is from a single cell. In some embodiments, the target nucleic acid is from a single organelle. In some embodiments, the target nucleic acid is genomic DNA. In some embodiments, the target nucleic acid is cross-linked to other nucleic acids. In some embodiments, target nucleic acid is from formalin fixed paraffin embedded (FFPE) sample. In some embodiments, the target nucleic acid is cross-linked with proteins. In some embodiments, the target nucleic acid is cross-linked with DNA. In some embodiments, the target nucleic acid is histone protected DNA. In some embodiments, histones are removed from the target nucleic acid. In some embodiments, the target nucleic acid is cell free tumor DNA. In some embodiments, the cell free tumor DNA is obtained from placental fluid. In some embodiments, the cell free tumor DNA is obtained from plasma. In some embodiments, the plasma is collected from whole blood using a membrane separator comprising a collection zone for the plasma. In some embodiments, the collection zone for the plasma comprises transposome complexes immobilized on solid support. In some embodiments, the target nucleic acid is cDNA. In some embodiments, the solid support is a bead. In some embodiments, the plurality of solid supports are plurality of beads and wherein the plurality of beads are of different sizes.
In some embodiments, a single barcode sequence is present in the plurality of immobilized oligonucleotides on each individual solid support. In some embodiments, different barcode sequences are present in the plurality of immobilized oligonucleotides on each individual solid support. In some embodiments, the transferring of the barcode sequence information to the target nucleic acid fragments is by ligation. In some embodiments, transferring of the barcode sequence information to the target nucleic acid fragments is by polymerase extension. In some embodiments, the transferring of the barcode sequence information to the target nucleic acid fragments is by both ligation and polymerase extension. In some embodiments, the polymerase extension is by extending the 3′-end of the non-ligated transposon strand with a DNA polymerase using the ligated immobilized oligonucleotide as a template. In some embodiments, at least a portion of the adaptor sequences further comprise a second barcode sequence.
In some embodiments, the transposome complexes are multimeric, and wherein the adaptor sequences of the transposons of each monomeric unit are different from the other monomeric unit in the same transposome complex. In some embodiments, the adaptor sequence further comprises a first primer binding sequence. In some embodiments, the first primer binding site has no sequence homology to the capture sequence or to the complement of the capture sequence. In some embodiments, the immobilized oligonucleotides on the solid support further comprise a second primer binding sequence.
In some embodiments, the transposome complexes are multimeric, and the transposome monomeric units are linked to each other in the same transposome complex. In some embodiments, the transposase of a transposome monomeric unit is linked to the transposase of another transposome monomeric unit of the same transposome complex. In some embodiments, the transposons of a transposome monomeric unit are linked to transposons of another transposome monomeric unit of the same transposome complex. In some embodiments, the transposase of a transposome monomeric unit is linked to the transposase of another transposome monomeric unit of the same transposome complex by covalent bond. In some embodiments, the transposases of one monomeric unit is linked to the transposase of another transposome monomeric unit of the same transposome complex by di-sulfide bond. In some embodiments, the transposons of a transposome monomeric unit are linked to transposons of another transposome monomeric unit of the same transposome complex by covalent bond.
In some embodiments, the contiguity information of a target nucleic acid sequence is indicative of haplotype information. In some embodiments, the contiguity information of a target nucleic acid sequence is indicative of genomic variants. In some embodiments, the genomic variants are selected from the group consisting of deletions, translocations, interchromosomal gene fusions, duplications, and paralogs. In some embodiments, the oligonucleotides immobilized on the solid support comprise a partially double stranded region and a partially single stranded region. In some embodiments, the partially single stranded region of the oligonucleotide comprises the second barcode sequence and the second primer binding sequence. In some embodiments, the target nucleic acid fragments comprising the barcodes are amplified prior to determining the sequence of the target nucleic acid fragments. In some embodiments, subsequent amplification are carried out in a single reaction compartment prior to determining the sequence of the target nucleic acid fragments. In some embodiments, a third barcode sequence is introduced to the target nucleic acid fragments during the amplification.
In some embodiments, the methods may further include combining the target nucleic acid fragments comprising the barcodes from plurality of first set of reaction compartments into a pool of target nucleic acid fragments comprising the barcodes; redistributing the pool of target nucleic acid fragments comprising the barcodes to a plurality of second set of reaction compartments; and introducing a third barcode in to the target nucleic acid fragments by amplifying the target nucleic acid fragments in the second set of reaction compartments prior to sequencing.
In some embodiments, the methods may further include pre-fragmenting the target nucleic acid prior to contacting the target nucleic acid with transposome complexes. In some embodiments, the pre-fragmenting the target nucleic acid is by a method selected from the group consisting of sonication and restriction digestion.
Embodiments of the present invention relate to sequencing nucleic acids. In particular, embodiments of the methods and compositions provided herein relate to preparing nucleic acid templates and obtaining sequence data therefrom.
In one aspect, the present invention relate to methods of tagmenting (fragmenting and tagging) target nucleic acid on a solid support for the construction of a tagmented target nucleic acid library. In one embodiment, the solid support is a bead. In one embodiment, the target nucleic acid is DNA.
In one aspect, the present invention relate to methods and compositions of solid-support, transposase-based methods that can derive contiguity information of a target nucleic acid. In some embodiments, the compositions and the methods can derive assembly/phasing information.
In one aspect, the present invention relate to methods and compositions to derive contiguity information by means of capturing contiguously-linked, transposed, target nucleic acid onto a solid support.
In one aspect the compositions and methods disclosed herein relate to analysis of genomic variants. Exemplary genomic variants include but are not limited to deletions, inter chromosomal translocations, duplications, paralogs, interchromosomal gene fusions. In some embodiments, the compositions and methods disclosed herein relate to determining phasing information of the genomic variants.
In one aspect, the compositions and methods disclosed herein relate to phasing specific regions of the target nucleic acid. In one embodiment, the target nucleic acid is DNA. In one embodiment, the target nucleic acid is genomic DNA. In some embodiments, the target nucleic acid is RNA. In some embodiments, the RNA is mRNA. In some embodiments, the target nucleic acid is complimentary DNA (cDNA). In some embodiments, target nucleic acid is from a single cell. In some embodiments, target nucleic acid is from circulating tumor cells. In some embodiments, target nucleic acid is cell free DNA. In some embodiments, target nucleic acid is cell free tumor DNA. In some embodiments, target nucleic acid is from formalin fixed paraffin embedded tissue samples. In some embodiments, target nucleic acid is cross-linked target nucleic acid. In some embodiments, target nucleic acid is cross-linked to proteins. In some embodiments, target nucleic acid is cross-linked to nucleic acid. In some embodiments, target nucleic acid is histone-protected DNA. In some embodiments, histone-protected DNA is precipitated from a cell lysate using antibodies to histones and the histones are removed.
In some aspects, indexed libraries are created from the target nucleic acid using the clonally indexed beads. In some embodiments, the tagmented target nucleic acid, while the transposase is still bound to the target DNA can be captured using the clonally indexed beads. In some embodiments, specific capture probes are used to capture the specific region of interest in the target nucleic acid. The captured regions of the target nucleic acid can be washed at various stringencies and optionally amplified, followed by sequencing. In some embodiments, the capture probe may be biotinylated. The complex of the biotinylated capture probes hybridized to the specific regions of the indexed target nucleic acids can be separated by using streptavidin beads. Exemplary scheme of targeted phasing is shown in
In some aspects, the compositions and methods disclosed herein can be used phasing exomes. In some embodiments, exons, promoters can be enriched. Markers, for example, heterozygous SNPs between exonic regions, can aid in phasing the exons, especially when the distance between exons is large. Exemplary exome phasing is shown in
In one aspect, the compositions and methods disclosed herein can be used for phasing and simultaneous methylation detection. Methylation detection through bisulfite conversion (BSC) is challenging as the BSC reaction is harsh on DNA, fragmenting the DNA and therefore removing contiguity/phasing information. Also, methods disclosed in the present application has an additional advantage because no additional purification steps are required in contrast to those required in traditional BSC approaches, thereby improving the yield.
In one aspect, the compositions and methods disclosed herein can be used to prepare different size libraries in single assay. In some embodiment, different sizes of clonally indexed beads can be used to prepare different size libraries.
In one embodiment, transposons may comprise sequencing primer binding sites. Exemplary sequences of sequence binding sites include, but are not limited to AATGATACGGCGACCACCGAGATCTACAC (P5 sequence) (SEQ ID NO: 1) and CAAGCAGAAGACGGCATACGAGAT (P7 sequence) (SEQ ID NO: 2). In some embodiments, the transposons may be biotinylated.
At a step 110 of
At a step 115 of
At a step 120 of
At step 110 of method 100, a plurality of biotinylated P5 transposons 210a and a plurality of P7 transposons 210b are generated. P5 transposons 210a and P7 transposons 210b are biotinylated.
At step 115 of method 100, P5 transposons 210a and P7 transposons 210b are mixed with transposase Tn5 215 to form a plurality of assembled transposomes 220.
At step 120 of method 100, transposomes 220 are bound to a bead 225. Bead 225 is a streptavidin coated bead. Transposomes 220 are bound to bead 225 via a biotin-streptavidin binding complex.
In one embodiment, a mixture of transposomes may be formed on a solid support such as bead surface as shown in
In another embodiment, the transposomes may be bound to any solid surface, such as the walls of a microfuge tube.
In another embodiment of forming a mixture of transposome complexes on a bead surface, oligonucleotides are first bound to a bead surface prior to transposome assembly. FIG. 10 illustrates a flow diagram of an example of a method 1000 of forming transposome complexes on a bead surface. Method 1000 includes, but is not limited to, the following steps.
At a step 1010, P5 and P7 oligonucleotides are bound to a bead surface. In one example, the P5 and P7 oligonucleotides are biotinylated and the bead is a streptavidin coated bead. This step is also shown pictorially in schematic diagram 1100 of
At a step 1015, complementary mosaic end (ME′) oligonucleotides are hybridized to the bead-bound P5 and P7 oligonucleotides. This step is also shown pictorially in schematic diagram 1200 of
At a step 1020, transposase enzyme is added to the bead-bound oligonucleotides to form a mixture of bead-bound transposome complexes. This step is also shown pictorially in schematic diagram 1300 of
The length of bridged molecules 1415 is independent of the quantity of beads 1120 with transposome complexes 1310 bound thereon used in a tagmentation reaction. Similarly, adding more or less DNA 1410 in a tagmentation reaction does not alter the size of the final tagmented product, but may affect the yield of the reaction.
In one example, bead 1120 is a paramagnetic bead. In this example, purification of the tagmentation reaction is readily achieved by immobilizing beads 1120 with a magnet and washing. Therefore, tagmentation and subsequent PCR amplification may be performed in a single reaction compartment (“one-pot”) reaction.
In one aspect, the present invention relate to methods and compositions of transposase-based methods that can derive contiguity information of a target nucleic acid on a solid support. In some embodiments, the compositions and the methods can derive assembly/phasing information. In one embodiment, the solid support is a bead. In one embodiment, the target nucleic acid is DNA. In one embodiment, the target nucleic acid is genomic DNA. In some embodiment, the target nucleic acid is RNA. In some embodiments, the RNA is mRNA. In some embodiments, the target nucleic acid is complimentary DNA (cDNA).
In some embodiments, transposons may be immobilized as dimers to solid-support such as beads, followed by the binding of transposase to the transposons to form transposomes.
In some embodiments, particularly related to formation of transposomes on solid-phases by solid-phase immobilized transposons and addition of transposase, two transposons may be immobilized in close proximity (preferably fixed distance) to one another in a solid support. There are several advantages to this approach. First, the two transposons will always be immobilized simultaneously, with preferably an optimum linker length and orientation of the two transposons to form transposomes efficiently. Second, transposome formation efficiency will not be a function of transposon density. Two transposons will always be available with the right orientation and distance between them to form transposomes. Third, with random immobilized transposons on surfaces, various distances are created between transposons, therefore only a fraction has the optimum orientation and distance to form transposomes efficiently. As a consequence, not all transposons are converted into transposomes and solid-phase immobilized non-complexed transposons will be present. These transposons are susceptible as a target to transposition as the ME-part is double-stranded DNA. This could result in a reduction of transposition efficiency or creates undesired side products. Thus, transposomes may be prepared on solid support, which can subsequently be used to derive contiguity information through tagmentation and sequencing. An exemplary scheme is illustrated in
In some embodiments, transposomes can be pre-assembled and then immobilized on a solid-support. In some embodiments, the transposons comprise unique indexes, barcodes, and amplification primer binding sites. Transposase can be added in solution comprising transposons to form transposome dimers, which can be immobilized on a solid support. In one embodiment, multiple bead sets can be generated in which each set has the same index derived from the immobilized transposons thus generating indexed beads. Target nucleic acid can be added to each set of indexed beads as shown in
In some embodiments, target nucleic acid can be added to each set of indexed beads, tagmented and subsequent PCR amplification may be performed separately.
In some embodiments, target nucleic acid, indexed beads, and transposomes can be combined in droplets such that a number of droplets contain a single bead with one or more DNA molecules and adequate transposomes.
In some embodiments, the indexed beads can be pooled, target nucleic acid can be added to the pool, tagmented and subsequent PCR amplification may be performed in a single reaction compartment (“one-pot”).
In one aspect, the present invention relate to methods and compositions to derive contiguity information by means of capturing contiguously-linked, transposed, target nucleic acid onto a solid support. In some embodiments, contiguity preserving transposition (CPT) is carried out on the DNA, but the DNA is kept intact (CPT-DNA), thus making contiguously-linked libraries. Contiguity information can be preserved by the use of transposase to maintain the association of template nucleic acid fragments adjacent in the target nucleic acid. The CPT-DNA can be captured by hybridization of complimentary oligonucleotides having unique indexes or barcodes and immobilized on solid support, e.g., beads (
Advantageously, such use of transposomes to maintain physical proximity of fragmented nucleic acids increases the likelihood that fragmented nucleic acids from the same original molecule, e.g., chromosome, will receive the same unique barcode and index information from the oligonucleotides immobilized on a solid support. This will result in a contiguously-linked sequencing library with unique barcodes. The contiguously-linked sequencing library can be sequenced to derive contiguous sequence information.
In some embodiments, the oligonucleotide immobilized on a solid support can comprise a partially double stranded structure such that one strand is immobilized to the solid support and the other strand is partially complementary to the immobilized strand resulting in a Y-adaptor.
In some embodiments, the Y-adaptor immobilized on the solid surface is linked to the contiguously linked tagmented DNA by ligation and gap filling and shown in
In some embodiments, Y-adaptor is formed through hybridization capture of CPT-DNA with the probe/index on the solid support such as beads.
In some embodiments, free transposomes may be separated from CPT-DNA. In some embodiments, the separation of the free transposomes is by size exclusion chromatography. In one embodiment, the separation may be achieved by MicroSpin S-400 HR Columns (GE Healthcare Life Sciences, Pittsburgh, Pa.).
Capturing contiguously-linked, transposed, target nucleic acid onto a solid support through hybridization has several unique advantages. First, the method is based on hybridization and not transposition. Intramolecular hybridization rate>>intermolecular hybridization rate. Thus, chances of contiguously-transposed libraries on a single target DNA molecule to wrap around a uniquely indexed bead is much higher than having two or more different single target DNA molecule to wrap around a uniquely indexed bead. Second, DNA transposition and barcoding of the transposed DNA occur in two separate steps. Third, the challenges associated with active transposome assembly on beads and surface density optimization of transposons on solid-surfaces can be avoided. Fourth, self-transposition products can be removed by column purification. Fifth, as contiguously linked, transposed, DNA contains gaps, the DNA is more flexible and therefore puts less of a burden on transposition density (insert size) compared to immobilizing transposome on bead methods. Sixth, the method can be used with combinatorial barcoding schemes. Seventh, it is easy to covalently-link indexed oligos to the beads. Thus, there is less chance for index exchange. Eight, the tagmentation and subsequent PCR amplification may be multiplexed and can be performed in a single reaction compartment (“one-pot”) reaction eliminating the need to carryout individual reactions for each index sequences.
In some embodiments, a plurality of unique barcodes throughout the target nucleic acid may be inserted during transposition. In some embodiments, each barcode includes a first barcode sequence and a second barcode sequence, having a fragmentation site disposed therebetween. The first barcode sequence and second barcode sequence can be identified or designated to be paired with one another. The pairing can be informative so that a first barcode is associated with a second barcode. Advantageously, the paired barcode sequences can be used to assemble sequencing data from the library of template nucleic acids. For example, identifying a first template nucleic acid comprising a first barcode sequence and a second template nucleic acid comprising a second barcode sequence that is paired with the first indicates that the first and second template nucleic acids represent sequences adjacent to one another in a sequence representation of the target nucleic acid. Such methods can be used to assemble a sequence representation of a target nucleic acid de novo, without the requirement of a reference genome.
In one aspect, the present invention relate to methods and compositions to generate shotgun sequence library of a specific DNA fragment.
In one embodiment, clonal indexed beads are generated with immobilized oligonucleotide sequences: random or specific primer and unique indexes. Target nucleic acid is added to the clonal indexed beads. In some embodiments, the target nucleic acid is DNA. In one embodiment, the target DNA is denatured. The target DNA hybridizes with primers comprising unique indexes immobilized on the solid surface (e.g., bead) and subsequently with other primers with the same index. The primers on the bead amplify the DNA. One or more further rounds of amplification may be carried out. In one embodiment, the amplification may be carried out by whole genome amplification using bead immobilized primers with a 3′ random n-mer sequence. In a preferred embodiment, the random n-mer contains pseudocomplementary bases (2-thiothymine, 2-amino dA, N4-ethyl cytosine, etc.) to prevent primer-primer interaction during amplification (Hoshika, S; Chen, F; Leal, N A; Benner, S A, Angew. Chem. Int. Ed. 49(32) 5554-5557 (2010).
The methods of the above embodiments have several advantages. Intra-molecular amplification on a bead is much faster than inter-bead amplification. Thus, the products on a bead will have the same index. A shotgun library of a specific DNA fragment can be created. Random primers amplify the template at random locations and therefore a shotgun library with the same index can be generated from a specific molecule and the sequence information can be assembled using the indexed sequence. A significant advantage of the methods of the above embodiments is that the reactions can be multiplexed in a single reaction (one pot reaction) and will not require using many individual wells. Many index clonal beads can be prepared so many different fragments can be uniquely labeled, and discrimination can be made to the parental alleles for same genomic regions. With a high number of indexes, the chance that the DNA copy of the father and copy of the mother will receive the same index for the same genomic region is low. The method takes advantage of the fact that intra reactions are much faster than inter, the beads basically generate a virtual partition in a larger physical compartment.
In some embodiments of all of the above aspect of the inventions, the method may be used for cell free DNA (cfDNA) in cfDNA assays. In some embodiments, the cfDNA is obtained from plasma, placental fluids.
In one embodiment, the plasma can be obtained from undiluted whole blood using membrane based, sedimentation assisted plasma separator (Liu et al. Anal Chem. 2013 Nov. 5; 85(21):10463-70). In one embodiment, the collection zone of the plasma of the plasma separator may comprise solid support comprising transposomes. The solid support comprising transposomes may capture the cfDNA from the isolated plasma as it is separated from the whole blood and can concentrate the cfDNA and/or tagment the DNA. In some embodiments, the tagmentation will further introduce unique barcodes to allow subsequent demultiplexing after sequencing of the pool of libraries.
In some embodiments, the collection zone of the separator may comprise PCR master mix (primers, nucleotides, buffers, metals) and polymerase. In one embodiment, the master mix can be in dry form such that it will be reconstituted as the plasma comes out of the separator. In some embodiments the primers are random primers. In some embodiments, the primers can be specific primers for a particular gene. PCR amplification of the cfDNA will result in the generation of library directly from the separated plasma.
In some embodiments, the collection zone of the separator may comprise RT-PCR master mix (primers, nucleotides, buffers, metals), reverse transcriptase and polymerase. In some embodiments the primers are random primers or oligo dT primers. In some embodiments, the primers can be specific primers for a particular gene. The resulting cDNA can be used for sequencing. Alternatively, the cDNA can be treated with transposomes immobilized on a solid support for sequence library preparation.
In some embodiments, the plasma separator may comprise barcodes (1D or 2D barcodes). In some embodiments, the separation device may comprise blood collection device. This would result in direct delivery of the blood to the plasma separator and library prep device. In some embodiments, the device may comprise a downstream sequence analyzer. In some embodiments, sequence analyzer is a single use sequencer. In some embodiments, the sequencer is capable of queuing samples before sequencing in a batch. Alternatively, the sequencer may have random access capability, where samples are delivered to their sequencing area.
In some embodiments, the collection zone for plasma may comprise silica substrates, such that the cell free DNA is concentrated.
Simultaneous Phasing and Methylation Detection
The 5-methyl Cytosine (5-Me-C) and 5-hydroxymethyl Cytosine (5-hydroxy-C), also known as epi modifications play an important role in cellular metabolism, differentiation and cancer development. Inventors of the present application has surprisingly and unexpectedly found that phasing and simultaneous methylation detection is possible using the methods and compositions of the present application. The present methods will allow to combine CPT-seq on beads (indexed contiguity linked libraries) with DNA methylation detection. For example, individual libraries generated on beads can be treated with bisulfite, converting non-methylated Cs, but not methylated Cs to Us, allowing the detection of 5-Me-C. Through additional phasing analysis using heterozygous SNPs, epi-medication-phasing blocks can be established multi megabase range.
In some embodiments, the size of the DNA analyzed can be about hundred bases to about multi mega bases. In some embodiments, the size of the DNA analyzed can be about 100, 200, 300, 400, 500, 600, 700, 800, 900, 1000, 1200, 1300, 1500, 2000, 3000, 3500, 4000, 4500, 5000, 5500, 6000, 6500, 7000, 7,500, 8000, 8500, 9000, 9500, 10,000, 10,500, 11,000, 11,500, 12,000, 12500, 13000, 14000, 14500, 15000, 15500, 16000, 16500, 17000, 17,500, 18,000, 18,500, 19,000, 19,500, 20,000, 20,500, 21,000, 21,500, 22,000, 22,500, 23,000, 23,500, 24,000, 24,500, 25,000, 25,500, 26,000, 26,500, 27,000, 27,500, 28,000, 28,500, 29,500, 30,000, 30,500, 31,000, 31,500, 32,000, 33,000, 34,000, 35,000, 36,000, 37,000, 38,000, 39,000, 40,000, 42,000, 45,000, 50,000, 55,000, 60,000, 65,000, 70,000, 75,000, 80,000, 85,000, 90,000, 95,000, 100,000, 110,000, 120,000, 130,000, 140,000, 150,000, 160,000, 170,000, 180,000, 200,000, 225,000, 250,000, 300,000, 350,000, 400,000, 450,000, 500,000, 550,000, 600,000, 650,000, 700,000, 750,000, 800,000, 850,000, 900,000, 1,000,000, 1,250,000, 1,500,000, 2,000,000, 2,500,000, 3,000,000, 4,000,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 15,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 75,000,000, 100,000,000 or more bases.
Other epi-modifications like 5-hydroxy-C, DNA oxidation products, DNA alkylation products, histone-foot printing etc. can also be analyzed in the context of phasing using the disclosed methods and compositions of the present application.
In some embodiments, DNA is first transformed into indexed-linked libraries on a solid-support. Individual indexed libraries, much smaller than the original DNA, are less prone to fragmentation since the individual libraries are smaller. Even if a small fraction of indexed libraries are lost, phasing information is still maintained across the long span of the indexed DNA molecule. For example, if a 100 kb molecule in traditional bisulfite conversion (BSC) is fragmented in half the contiguity is now restricted to 50 kb. In the methods disclosed herein, a 100 kb library is first indexed and even if a fraction of individual libraries are lost, contiguity is still at ˜100 kb (except in the unlikely event when all libraries lost are from one end of the DNA molecule. Also, methods disclosed in the present application has an additional advantage because no additional purification steps are required in contrast to those required in traditional bisulfite conversion approaches, thereby improving the yield. In the methods of the present application, the beads are simply washed after bisulfite conversion. Additionally, while DNA is bound to a solid phase, buffer exchanges can be readily performed with minimal loss of DNA (indexed libraries) and reduced hands on time.
Exemplary scheme of simultaneous phasing and methylation detection is shown in
One challenge to performing bisulfite conversion of DNA bound to a solid phase, such as streptavidin magnetic beads is that extended treatment of bead bound DNA with sodium bisulfite at high temperatures damages both the DNA and the beads. To help ameliorate DNA damage, carrier DNA (i.e. Lambda DNA) is added to the reaction mixture prior to bisulfite treatment. Even in presence of carrier DNA, it has been estimated that approximately 80% of starting DNA is lost. As a result, CPTSeq contiguity blocks have fewer members than those in the traditional CPTSeq protocol.
Therefore, several strategies are proposed herein to improve DNA yield of the Epi-CPTSeq protocol. The first strategy relies on decreasing library insert size by more densely populating transposome complexes to the streptavidin beads. By decreasing library size, a smaller proportion of library elements are degraded by bisulfite treatment.
The second strategy to improve DNA yield of the Epi-CPTSeq protocol is enzymatic recovery of broken library elements. The purpose of the recovery strategy is to add the 3′ common sequence necessary for library amplification back to the bead bound library elements that became digested and lost their 3′ portion during bisulfite treatment. After the addition of the 3′ common sequence these elements can now be PCR amplified and sequenced.
In an alternative workflow (
Both of the described workflows rely on a controllable TdT tailing reaction recently developed and described in US Patent Application Publication 20150087027. A common sequencing adapter can also be added to the 3′ end of broken library elements by a recently introduced ssDNA template switching activity of MMLV RT. In short, MMLV RT and a template switch oligo (TS_oligo) are added to damaged DNA (
As a part of the third strategy, an Epicentre's EpiGenome kit “post-bisulfite conversion” library construction method can be used to rescue library elements which lost their common sequences at the 3′ end during bisulfite conversion. As shown in
Next, DNA is denatured (e.g. incubation at high heat) and bound to a solid support. If biotin is used as a capture tag on CS1, for example, DNA can be bound using streptavidin magnetic beads (as pictured). Once bound to the solid support buffer exchanges can be readily made.
In the next step, bisulfite conversion of ssDNA is performed. In the single stranded form, DNA should be readily accessible for bisulfite conversion; up to 95% conversion efficiencies have been observed using a modified version of Promega's Methyl Edge BSC kit (
After bisulfite conversion, a second common sequence is covalently attached to the 3′ end of ssDNA attached to solid support. Several methods have been described above to covalently attach oligos to ssDNA. Using the TdT attenuator/adapter ligation method, ligation efficiencies of >95% have been achieved. As a result, final library yields using the proposed MethylSeq workflow should be greater than existing methods.
In the final step, PCR is performed to amplify the library and remove it from the solid support. PCR primers can be designed to add additional commons sequences, such as sequencing adapters, to the ends of the MethylSeq library.
Preparation of Different Size Libraries in a Single Assay
The accuracy of the assembly of genomes is contingent on the use of different length scale technologies. For example, shotgun (100's of bp)-matepair (˜3 Kb) to -Hi-C (Mb-scale) are all methods that sequentially improve assemblies and contig lengths. The challenge is that multiple assays are required to accomplish this, making the multi-layered approach cumbersome and costly. The compositions and methods disclosed herein can address multiple length scales in a single assay.
In some embodiments, library preparation can be achieved in a single assay using differentially sized solid support, for example, beads. Each bead size will generate a specific library size or range of sizes, with the physical size of the bead determining the library size. The various sized beads all have unique clonal indices that are transferred to the library. As such, different sizes libraries are generated with each different library scale-length uniquely indexed. The various length-scale libraries are prepared simultaneously in the same physical compartment, reducing cost and improving overall work flow. In some embodiments, each specific solid support size, for example, bead size receives a unique index. In some other embodiments, multiple different indexes of the same solid support size, for example, bead size are also prepared so multiple DNA molecules can be index partitioned for that size range.
In some embodiments, the size of the libraries generated are about 50, 75, 100, 150, 200, 250, 300, 350, 400, 500, 600, 700, 800, 900, 1000, 1200, 1300, 1500, 2000, 3000, 3500, 4000, 4500, 5000, 5500, 6000, 6500, 7000, 7,500, 8000, 8500, 9000, 9500, 10,000, 10,500, 11,000, 11,500, 12,000, 12500, 13000, 14000, 14500, 15000, 15500, 16000, 16500, 17000, 17,500, 18,000, 18,500, 19,000, 19,500, 20,000, 20,500, 21,000, 21,500, 22,000, 22,500, 23,000, 23,500, 24,000, 24,500, 25,000, 25,500, 26,000, 26,500, 27,000, 27,500, 28,000, 28,500, 29,500, 30,000, 30,500, 31,000, 31,500, 32,000, 33,000, 34,000, 35,000, 36,000, 37,000, 38,000, 39,000, 40,000, 42,000, 45,000, 50,000, 55,000, 60,000, 65,000, 70,000, 75,000, 80,000, 85,000, 90,000, 95,000, 100,000, 110,000, 120,000, 130,000, 140,000, 150,000, 160,000, 170,000, 180,000, 200,000, 225,000, 250,000, 300,000, 350,000, 400,000, 450,000, 500,000, 550,000, 600,000, 650,000, 700,000, 750,000, 800,000, 850,000, 900,000, 1,000,000, 1,250,000, 1,500,000, 2,000,000, 2,500,000, 3,000,000, 4,000,000, 5,000,000, 6,000,000, 7,000,000, 8,000,000, 9,000,000, 10,000,000, 15,000,000, 20,000,000, 30,000,000, 40,000,000, 50,000,000, 75,000,000, 100,000,000 or more bases.
In some embodiments, multiple length scale libraries discussed above can be used in the assembly of pseudogenes, paralogs etc. instead of having one large length scale. In some embodiments, multiple length scale libraries are prepared simultaneously in a single assay. The advantage is that at least one length-scale will link a unique region with only the pseudo-gene and or gene, but not both. As such, variants detected with this length-scale can uniquely assign the variant to either the gene or the pseudo-gene. The same holds true for copy number variants, paralogs etc. The strength of assembly is the use of different length scales. Using the methods disclosed herein different length scale indexed linked libraries can be generated in a single assay instead of individual, different library preparations for different length scales.
Analysis of Genomic Variants
The compositions and methods disclosed herein relate to analysis of genomic variants. Exemplary genomic variants include but are not limited to deletions, inter chromosomal translocations, duplications, paralogs, interchromosomal gene fusions. In some embodiments, the compositions and methods disclosed herein relate to determining phasing information of the genomic variants. The table below shows exemplary interchromosomal gene fusions.
Table 2 shows exemplary deletions in chromosome 1:
In some embodiments, target nucleic acid can be fragmented prior to exposing it to transposomes. Exemplary fragmentation methods include, but are not limited to sonication, mechanical shearing, and restriction digestion. Fragmentation of target nucleic acid prior to tagmentation (fragmentation and tagging) is advantageous for assembly/phasing of pseudogenes (e.g., CYP2D6). Long islands (>30 kb) of indexed linked reads will span the pseudogenes A and A′ as shown in
Linked Transposomes
In some embodiments, transposases are multimeric in a transposome complex, e.g., they form dimers, tetramers etc. in a transposome complex. Inventors of the present application have surprisingly and unexpectedly found that linking the monomer transposases in multimeric transposome complex or linking the transposon ends of a transposome monomer in multimeric transposome complex has several advantages. First, the linking of the transposases or the transposons leads to the complexes that are more stable and a large fraction is in an active state. Second, lower concentrations of transposomes can potentially be used in the fragmentation by transposition reaction. Third, the linking leads to lower exchange of the mosaic ends (ME) of transposome complexes, thus less mixing of barcodes or adaptor molecules. Such swapping of ME ends are possible if the complexes fall apart and reform, or in case where transposomes are immobilized on solid support by streptavidin/biotin, the streptavidin/biotin interaction can break and reform, or when there is a possible contamination. Inventors of the present application noted that there is a significant swap or exchange of ME ends under various reaction conditions. In some embodiments, the exchange can be as high as 15%. The exchange is pronounced in high salt buffer and the exchange is reduced in glutamate buffer.
In some embodiments, the transposase subunits in the transposome complex can be linked to each other by covalent and non-covalent means. In some embodiments, transposase monomers can be linked before making the transposome complex (before addition of the transposons). In some embodiments, transposase monomers can be linked after transposome formation.
In some embodiments, native amino acid residues may be substituted with Cysteine (Cys) amino acids at the multimeric interface to promote disulfide bond formation. For example, in Tn5 transposase, Asp468, Tyr407, Asp461, Lys459, Ser458, Gly462, Ala466, Met470 may be substituted with Cys to promote disulfide bond between the monomer subunits and shown in
In some embodiments, transposome multimer complexes can be covalently linked to solid support. Exemplary solid supports include but are not limited to nanoparticles, beads, flow cell surfaces, column matrices. In some embodiments, solid surfaces may be coated with amine groups. Modified transposase with amino acid residues substituted with cysteine can be chemically cross-linked to such amine groups using an amine-to-sulfhydryl crosslinker (i.e., succinimidyl-4-(N-maleimidomethyl)cyclohexane-1-carboxylate (SMCC)). Exemplary scheme is shown in
In some embodiments, transposase gene can be modified to express multimeric protein in a single polypeptide. For example, Tn5 or Mos-1 genes can be modified to express two Tn5 or Mos-1 proteins in a single polypeptide. Similarly Mu transposase gene can be modified to encode four mu transposase units in a single polypeptide.
In some embodiments, the transposon ends of a transposome monomer unit can be linked to form a linked transposome multimer complex. Linking the transposon ends allow insertion of primer sites, sequencing primers, amplification primers or any role DNA can play into gDNA without fragmenting the target DNA. Insertion of such functionality are advantages in haplotyping assays or junction tagging assays in which information needs to be extracted from intact molecules or in which sub-sampling are important. In some embodiments, transposon ends of Mu transposomes can be linked to a “looped” Mu transposase/transposon configuration. Since Mu is a tetramer, various configurations are possible but not limited by linking either R2UJ and/or R1UJ with R2J and/or R1J. In these configurations R2UJ and R1UJ can/are not connected with R2J and R1J, respectively.
As used herein the term “transposon” means a double-stranded DNA that exhibits only the nucleotide sequences (the “transposon end sequences”) that are necessary to form the complex with the transposase or integrase enzyme that is functional in an in vitro transposition reaction. A transposon forms a “complex” or a “synaptic complex” or a “transposome complex” or a “transposome composition with a transposase or integrase that recognizes and binds to the transposon, and which complex is capable of inserting or transposing the transposon into target DNA with which it is incubated in an in vitro transposition reaction. A transposon exhibits two complementary sequences consisting of a “transferred transposon sequence” or “transferred strand” and a “non-transferred transposon sequence,” or “non transferred strand”. For example, one transposon that forms a complex with a hyperactive Tn5 transposase (e.g., EZ-Tn5™ Transposase, EPICENTRE Biotechnologies, Madison, Wis., USA) that is active in an in vitro transposition reaction comprises a transferred strand that exhibits a “transferred transposon sequence” as follows:
and a non-transferred strand that exhibits a “non-transferred transposon sequence” as follows:
The 3′-end of a transferred strand is joined or transferred to target DNA in an in vitro transposition reaction. The non-transferred strand, which exhibits a transposon sequence that is complementary to the transferred transposon end sequence, is not joined or transferred to the target DNA in an in vitro transposition reaction. In some embodiments, the transposon sequences may comprise one or more of the following: a barcode, an adaptor sequence, a tag sequence, a primer binding sequence, a capture sequence, unique molecular identifier (UMI) sequence.
As used herein the term “adaptor” means a nucleic acid sequence that can comprise a barcode, a primer binding sequence, a capture sequence, a sequence complementary to a capture sequence, unique molecular identifier (UMI) sequence, an affinity moiety, restriction site.
As used herein the term “contiguity information” refers to a spatial relationship between two or more DNA fragments based on shared information. The shared aspect of the information can be with respect to adjacent, compartmental and distance spatial relationships. Information regarding these relationships in turn facilitates hierarchical assembly or mapping of sequence reads derived from the DNA fragments. This contiguity information improves the efficiency and accuracy of such assembly or mapping because traditional assembly or mapping methods used in association with conventional shotgun sequencing do not take into account the relative genomic origins or coordinates of the individual sequence reads as they relate to the spatial relationship between the two or more DNA fragments from which the individual sequence reads were derived. Therefore, according to the embodiments described herein, methods of capturing contiguity information may be accomplished by short range contiguity methods to determine adjacent spatial relationships, mid-range contiguity methods to determine compartmental spatial relationships, or long range contiguity methods to determine distance spatial relationships. These methods facilitate the accuracy and quality of DNA sequence assembly or mapping, and may be used with any sequencing method, such as those described above.
Contiguity information includes the relative genomic origins or coordinates of the individual sequence reads as they relate to the spatial relationship between the two or more DNA fragments from which the individual sequence reads were derived. In some embodiments, contiguity information includes sequence information from non-overlapping sequence reads.
In some embodiments, the contiguity information of a target nucleic acid sequence is indicative of haplotype information. In some embodiments, the contiguity information of a target nucleic acid sequence is indicative of genomic variants.
As used herein the term “maintaining the contiguity of the target nucleic acid” in the context of fragmenting a nucleic acid means maintaining the order of the nucleic acid sequence of the fragments from the same target nucleic acid.
As used herein the term “at least a portion” and/or grammatical equivalents thereof can refer to any fraction of a whole amount. For example, “at least a portion” can refer to at least about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 99%, 99.9% or 100% of a whole amount.
As used herein the term “about” means +/−10%.
As used herein, the term “sequencing read” and/or grammatical equivalents thereof can refer to a repetitive process of physical or chemical steps that is carried out to obtain signals indicative of the order of monomers in a polymer. The signals can be indicative of an order of monomers at single monomer resolution or lower resolution. In particular embodiments, the steps can be initiated on a nucleic acid target and carried out to obtain signals indicative of the order of bases in the nucleic acid target. The process can be carried out to its typical completion, which is usually defined by the point at which signals from the process can no longer distinguish bases of the target with a reasonable level of certainty. If desired, completion can occur earlier, for example, once a desired amount of sequence information has been obtained. A sequencing read can be carried out on a single target nucleic acid molecule or simultaneously on a population of target nucleic acid molecules having the same sequence, or simultaneously on a population of target nucleic acids having different sequences. In some embodiments, a sequencing read is terminated when signals are no longer obtained from one or more target nucleic acid molecules from which signal acquisition was initiated. For example, a sequencing read can be initiated for one or more target nucleic acid molecules that are present on a solid phase substrate and terminated upon removal of the one or more target nucleic acid molecules from the substrate. Sequencing can be terminated by otherwise ceasing detection of the target nucleic acids that were present on the substrate when the sequencing run was initiated. Exemplary methods of sequencing are described in U.S. Pat. No. 9,029,103, which is incorporated herein by reference in its entirety.
As used herein, the term “sequencing representation” and/or grammatical equivalents thereof can refer to information that signifies the order and type of monomeric units in the polymer. For example, the information can indicate the order and type of nucleotides in a nucleic acid. The information can be in any of a variety of formats including, for example, a depiction, image, electronic medium, series of symbols, series of numbers, series of letters, series of colors, etc. The information can be at single monomer resolution or at lower resolution. An exemplary polymer is a nucleic acid, such as DNA or RNA, having nucleotide units. A series of “A,” “T,” “G,” and “C” letters is a well-known sequence representation for DNA that can be correlated, at single nucleotide resolution, with the actual sequence of a DNA molecule. Other exemplary polymers are proteins having amino acid units and polysaccharides having saccharide units.
Solid Support
Throughout this application, solid support and solid surface are used interchangeably. In some embodiments, the solid support or its surface is non-planar, such as the inner or outer surface of a tube or vessel. In some embodiments, the solid support comprises microspheres or beads. By “microspheres” or “beads” or “particles” or grammatical equivalents herein is meant small discrete particles. Suitable bead compositions include, but are not limited to, plastics, ceramics, glass, polystyrene, methylstyrene, acrylic polymers, paramagnetic materials, thoria sol, carbon graphite, titanium dioxide, latex or cross-linked dextrans such as Sepharose, cellulose, nylon, cross-linked micelles and Teflon, as well as any other materials outlined herein for solid supports may all be used. “Microsphere Detection Guide” from Bangs Laboratories, Fishers Ind. is a helpful guide. In certain embodiments, the microspheres are magnetic microspheres or beads. In some embodiments, the beads can be color coded. For example, MicroPlex® Microspheres from Luminex, Austin, Tex. may be used.
The beads need not be spherical; irregular particles may be used. Alternatively or additionally, the beads may be porous. The bead sizes range from nanometers, i.e. about 10 nm, to millimeters in diameter, i.e. 1 mm, with beads from about 0.2 micron to about 200 microns being preferred, and from about 0.5 to about 5 micron being particularly preferred, although in some embodiments smaller or larger beads may be used. In some embodiments, beads can be about 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.5, 2, 2.5, 2.8, 3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10, 10.5, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 150, or 200 μm in diameter.
Transposomes
A “transposome” comprises an integration enzyme such as an integrase or transposase, and a nucleic acid comprising an integration recognition site, such as a transposase recognition site. In embodiments provided herein, the transposase can form a functional complex with a transposase recognition site that is capable of catalyzing a transposition reaction. The transposase may bind to the transposase recognition site and insert the transposase recognition site into a target nucleic acid in a process sometimes termed “tagmentation”. In some such insertion events, one strand of the transposase recognition site may be transferred into the target nucleic acid. In one example, a transposome comprises a dimeric transposase comprising two subunits, and two non-contiguous transposon sequences. In another example, a transposome comprises a transposase comprises a dimeric transposase comprising two subunits, and a contiguous transposon sequence.
Some embodiments can include the use of a hyperactive Tn5 transposase and a Tn5-type transposase recognition site (Goryshin and Reznikoff, J. Biol. Chem., 273:7367 (1998)), or MuA transposase and a Mu transposase recognition site comprising R1 and R2 end sequences (Mizuuchi, K., Cell, 35: 785, 1983; Savilahti, H, et al., EMBO J., 14: 4893, 1995). An exemplary transposase recognition site that forms a complex with a hyperactive Tn5 transposase (e.g., EZ-Tn5™ Transposase, Epicentre Biotechnologies, Madison, Wis.) comprises the following 19b transferred strand (sometimes “M” or “ME”) and non-transferred strands: 5′ AGATGTGTATAAGAGACAG 3′ (SEQ ID NO: 3), 5′ CTGTCT CTTATACACATCT 3′ (SEQ ID NO: 4), respectively. ME sequences can also be used as optimized by a skilled artisan.
More examples of transposition systems that can be used with certain embodiments of the compositions and methods provided herein include Staphylococcus aureus Tn552 (Colegio et al., J. Bacteriol., 183: 2384-8, 2001; Kirby C et al., Mol. Microbiol., 43: 173-86, 2002), Ty1 (Devine & Boeke, Nucleic Acids Res., 22: 3765-72, 1994 and International Publication WO 95/23875), Transposon Tn7 (Craig, N L, Science. 271: 1512, 1996; Craig, N L, Review in: Curr Top Microbiol Immunol., 204:27-48, 1996), Tn/O and IS10 (Kleckner N, et al., Curr Top Microbiol Immunol., 204:49-82, 1996), Mariner transposase (Lampe D J, et al., EMBO J., 15: 5470-9, 1996), Tc1 (Plasterk R H, Curr. Topics Microbiol. Immunol., 204: 125-43, 1996), P Element (Gloor, G B, Methods Mol. Biol., 260: 97-114, 2004), Tn3 (Ichikawa & Ohtsubo, J Biol. Chem. 265:18829-32, 1990), bacterial insertion sequences (Ohtsubo & Sekine, Curr. Top. Microbiol. Immunol. 204: 1-26, 1996), retroviruses (Brown, et al., Proc Natl Acad Sci USA, 86:2525-9, 1989), and retrotransposon of yeast (Boeke & Corces, Annu Rev Microbiol. 43:403-34, 1989). More examples include IS5, Tn10, Tn903, IS911, Sleeping Beauty, SPIN, hAT, PiggyBac, Hermes, TcBuster, AeBuster1, Tol2, and engineered versions of transposase family enzymes (Zhang et al., (2009) PLoS Genet. 5:e1000689. Epub 2009 Oct. 16; Wilson C. et al (2007) J. Microbiol. Methods 71:332-5).
More examples of integrases that may be used with the methods and compositions provided herein include retroviral integrases and integrase recognition sequences for such retroviral integrases, such as integrases from HIV-1, HIV-2, SIV, PFV-1, RSV.
Barcodes
Generally, a barcode can include one or more nucleotide sequences that can be used to identify one or more particular nucleic acids. The barcode can be an artificial sequence, or can be a naturally occurring sequence generated during transposition, such as identical flanking genomic DNA sequences (g-codes) at the end of formerly juxtaposed DNA fragments. In some embodiments, the barcodes are artificial sequences that are absent in the target nucleic acid sequence and can be used to identify one or more target nucleic acid sequences.
A barcode can comprise at least about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 or more consecutive nucleotides. In some embodiments, a barcode comprises at least about 10, 20, 30, 40, 50, 60, 70 80, 90, 100 or more consecutive nucleotides. In some embodiments, at least a portion of the barcodes in a population of nucleic acids comprising barcodes is different. In some embodiments, at least about 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95%, 99% of the barcodes are different. In more such embodiments, all of the barcodes are different. The diversity of different barcodes in a population of nucleic acids comprising barcodes can be randomly generated or non-randomly generated.
In some embodiments, a transposon sequence comprises at least one barcode. In some embodiments, such as transposomes comprising two non-contiguous transposon sequences, the first transposon sequence comprises a first barcode, and the second transposon sequence comprises a second barcode. In some embodiments, a transposon sequence comprises a barcode comprising a first barcode sequence and a second barcode sequence. In some of the foregoing embodiments, the first barcode sequence can be identified or designated to be paired with the second barcode sequence. For example, a known first barcode sequence can be known to be paired with a known second barcode sequence using a reference table comprising a plurality of first and second bar code sequences known to be paired to one another.
In another example, the first barcode sequence can comprise the same sequence as the second barcode sequence. In another example, the first barcode sequence can comprise the reverse complement of the second barcode sequence. In some embodiments, the first barcode sequence and the second barcode sequence are different. The first and second barcode sequences may comprise a bi-code.
In some embodiments of compositions and methods described herein, barcodes are used in the preparation of template nucleic acids. As will be understood, the vast number of available barcodes permits each template nucleic acid molecule to comprise a unique identification. Unique identification of each molecule in a mixture of template nucleic acids can be used in several applications. For example, uniquely identified molecules can be applied to identify individual nucleic acid molecules, in samples having multiple chromosomes, in genomes, in cells, in cell types, in cell disease states, and in species, for example, in haplotype sequencing, in parental allele discrimination, in metagenomic sequencing, and in sample sequencing of a genome.
Exemplary barcode sequences include, but are not limited to TATAGCCT, ATAGAGGC, CCTATCCT, GGCTCTGA, AGGCGAAG, TAATCTTA, CAGGACGT, and GTACTGAC.
Primer Sites
In some embodiments, a transposon sequence can include a “sequencing adaptor” or “sequencing adaptor site”, that is to say a region that comprises one or more sites that can hybridize to a primer. In some embodiments, a transposon sequence can include at least a first primer site useful for amplification, sequencing, and the like. Exemplary sequences of sequence binding sites include, but are not limited to AATGATACGGCGACCACCGAGATCTACAC (SEQ ID NO: 1) (P5 sequence) and CAAGCAGAAGACGGCATACGAGAT (SEQ ID NO: 2) (P7 sequence).
Target Nucleic Acids
A target nucleic acid can include any nucleic acid of interest. Target nucleic acids can include DNA, RNA, peptide nucleic acid, morpholino nucleic acid, locked nucleic acid, glycol nucleic acid, threose nucleic acid, mixed samples of nucleic acids, polyploidy DNA (i.e., plant DNA), mixtures thereof, and hybrids thereof. In a preferred embodiment, genomic DNA or amplified copies thereof are used as the target nucleic acid. In another preferred embodiment, cDNA, mitochondrial DNA or chloroplast DNA is used. In some embodiments, the target nucleic acid is mRNA.
In some embodiments, target nucleic acid is from a single cell or from fractions of a single cell. In some embodiments, the target nucleic acid is from a single organelle. Exemplary single organelle includes but is not limited to single nuclei, single mitochondria, and a single ribosome. In some embodiments, target nucleic acid is from formalin fixed paraffin embedded (FFPE) sample. In some embodiments, target nucleic acid is cross-linked nucleic acid. In some embodiments, the target nucleic acid is cross-linked with protein. In some embodiments, the target nucleic acid is cross-linked DNA. In some embodiments, the target nucleic acid is histone protected DNA. In some embodiments, histones are removed from the target nucleic acid. In some embodiments, target nucleic acid is from nucleosomes. In some embodiments, target nucleic acid is from nucleosomes from which nuclear proteins are removed.
A target nucleic acid can comprise any nucleotide sequence. In some embodiments, the target nucleic acid comprises homopolymer sequences. A target nucleic acid can also include repeat sequences. Repeat sequences can be any of a variety of lengths including, for example, 2, 5, 10, 20, 30, 40, 50, 100, 250, 500 or 1000 nucleotides or more. Repeat sequences can be repeated, either contiguously or non-contiguously, any of a variety of times including, for example, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15 or 20 times or more.
Some embodiments described herein can utilize a single target nucleic acid. Other embodiments can utilize a plurality of target nucleic acids. In such embodiments, a plurality of target nucleic acids can include a plurality of the same target nucleic acids, a plurality of different target nucleic acids where some target nucleic acids are the same, or a plurality of target nucleic acids where all target nucleic acids are different. Embodiments that utilize a plurality of target nucleic acids can be carried out in multiplex formats so that reagents are delivered simultaneously to the target nucleic acids, for example, in one or more chambers or on an array surface. In some embodiments, the plurality of target nucleic acids can include substantially all of a particular organism's genome. The plurality of target nucleic acids can include at least a portion of a particular organism's genome including, for example, at least about 1%, 5%, 10%, 25%, 50%, 75%, 80%, 85%, 90%, 95%, or 99% of the genome. In particular embodiments the portion can have an upper limit that is at most about 1%, 5%, 10%, 25%, 50%, 75%, 80%, 85%, 90%, 95%, or 99% of the genome.
Target nucleic acids can be obtained from any source. For example, target nucleic acids may be prepared from nucleic acid molecules obtained from a single organism or from populations of nucleic acid molecules obtained from natural sources that include one or more organisms. Sources of nucleic acid molecules include, but are not limited to, organelles, cells, tissues, organs, or organisms. Cells that may be used as sources of target nucleic acid molecules may be prokaryotic (bacterial cells, for example, Escherichia, Bacillus, Serratia, Salmonella, Staphylococcus, Streptococcus, Clostridium, Chlamydia, Neisseria, Treponema, Mycoplasma, Borrelia, Legionella, Pseudomonas, Mycobacterium, Helicobacter, Erwinia, Agrobacterium, Rhizobium, and Streptomyces genera); archeaon, such as crenarchaeota, nanoarchaeota or euryarchaeotia; or eukaryotic such as fungi, (for example, yeasts), plants, protozoans and other parasites, and animals (including insects (for example, Drosophila spp.), nematodes (e.g., Caenorhabditis elegans), and mammals (for example, rat, mouse, monkey, non-human primate and human). Target nucleic acids and template nucleic acids can be enriched for certain sequences of interest using various methods well known in the art. Examples of such methods are provided in Int. Pub. No. WO2012/108864, which is incorporated herein by reference in its entirety. In some embodiments, nucleic acids may be further enriched during methods of preparing template libraries. For example, nucleic acids may be enriched for certain sequences, before insertion of transposomes after insertion of transposomes and/or after amplification of nucleic acids.
In addition, in some embodiments, target nucleic acids and/or template nucleic acids can be highly purified, for example, nucleic acids can be at least about 70%, 80%, 90%, 95%, 96%, 97%, 98%, 99%, or 100% free from contaminants before use with the methods provided herein. In some embodiments, it is beneficial to use methods known in the art that maintain the quality and size of the target nucleic acid, for example isolation and/or direct transposition of target DNA may be performed using agarose plugs. Transposition can also be performed directly in cells, with population of cells, lysates, and non-purified DNA.
In some embodiments, target nucleic acid may be obtained from a biological sample or a patient sample. The term “biological sample” or “patient sample” as used herein includes samples such as tissues and bodily fluids. “Bodily fluids” may include, but are not limited to, blood, serum, plasma, saliva, cerebral spinal fluid, pleural fluid, tears, lactal duct fluid, lymph, sputum, urine, amniotic fluid, and semen. A sample may include a bodily fluid that is “acellular.” An “acellular bodily fluid” includes less than about 1% (w/w) whole cellular material. Plasma or serum are examples of acellular bodily fluids. A sample may include a specimen of natural or synthetic origin (i.e., a cellular sample made to be acellular).
In some embodiments of the above disclosed methods, target nucleic acid can be fragmented (e.g., by sonication, by restriction digestion, other mechanical means) prior to exposing the target nucleic acid to the transposomes.
The term “Plasma” as used herein refers to acellular fluid found in blood. “Plasma” may be obtained from blood by removing whole cellular material from blood by methods known in the art (e.g., centrifugation, filtration, and the like).
Unless otherwise specified, the terms “a” or “an” mean “one or more” throughout this application.
When the terms “for example”, “e.g.”, “such as”, “include”, “including” or variations thereof are used herein, these terms will not be deemed to be terms of limitation, and will be interpreted to mean “but not limited to” or “without limitation.”
The following Examples provide illustrative embodiments and do not in any way limit the inventions provided herein.
EXAMPLES Example 1—DNA Cluster Yield from the Bead-Based Tagmentation ProcessDNA cluster yield from the bead-based tagmentation process of
The reproducibility of the bead-based tagmentation process of
The insert size of pool 1 and the insert size of pool 2 are shown in
The reproducibility of total number of reads and percent reads aligned for the experiment described in
In one application, the bead-based tagmentation process may be used in an exome enrichment assay which includes a tagmentation step, e.g., Illumina's Nextera® Rapid Capture Enrichment protocol. In the current exome enrichment assay (i.e., Illumina's Nextera® Rapid Capture Enrichment protocol), solution-based tagmentation (Nextera) is used to fragment the genomic DNA. Gene specific primers are then used to pull down specific gene fragments of interest. Two enrichment cycles are performed and fragments pulled down are then enriched by PCR and sequenced.
To evaluate the use of the bead-based tagmentation process in the exome enrichment assay, human NA12878 DNA was tagmented using 25, 50, 100, 150, 200, and 500 ng of input DNA. A control library (NA00536) was prepared from 50 ng input DNA according to the standard protocol. Each DNA input had a different index (unique identifier). Ten cycles of PCR using enhanced polymerase mastermix (EPM) were used to match standard methods and to ensure a sufficient amount of fragments were present for pulldown. The amplification protocol was 3 minutes at 72° C., 30 seconds at 98° C., followed by 10 cycles of 10 seconds at 98° C., 30 seconds at 65° C., and 1 minute at 72° C. The samples were then held at 10° C. The samples were then processed through the exome enrichment pulldown process and sequenced.
Example 5—Insert Size in a Control and Bead-Based Tagmented Library in the Exome Enrichment AssayTwo techniques may be used to optimize the insert size distribution. In one example, an SPRI clean-up may be used to remove fragments that are too small or too large. SPRI clean-up is a process of removing fragments that are larger or smaller than the desired size, by selective DNA precipitation based on size and either retention of the precipitated or non-precipitated DNA as desired (i.e., a first step is to precipitate only DNA that is larger than the desired size and retain the soluble smaller fragments). The smaller fragments are then further precipitated and this time the very small fragments that are not wanted (still in solution) are removed and the precipitated DNA is retained, washed and then resolubilized to give a desired size range of DNA. In another example, the spacing of active transposomes on the bead surface may be used to control the insert size distribution. For example, gaps on the bead surface may be filled with inactive transposomes (e.g., transposomes with inactive transposons).
Contiguity of the bead-based tagmentation process was assessed. Table 3 shows the number of times 0, 1, 2, or 3 reads occur within a 1000 bp windows sharing an index. Beads were generated with 9 different indexed transposomes and used to tagment a small amount of human DNA. Reads were generated, aligned, and analyzed for the number of reads within a 1000 bp or 10 Kb window that shared the same index. Some reads within a small window sharing an index may be generated by chance and a prediction of how many times this is likely to occur is given in the “Random” row of Table 3 and Table 4. The numbers in the “Bead” row show the actual number of 1000 bp (Table 3) or 10 Kb (Table 4) windows that share an index. As shown in Table 3 and Table 4, the actual number of times the same index was found within 1000 bp or 10 Kb window is significantly greater than expected in the random case. “0” windows show all the times a particular 1000 bp window had no indexed reads mapping to it. The number is largest here because only a very small amount of the human genome was sequence and most windows have no reads aligning to them. “1” is the number of times just one read maps to a 1000 bp (or 10 Kb) window; “2” the number of times 2 reads share an index within a 1000 bp (or 10 KB) window, etc. This data suggests that in over 1400 cases the same piece of DNA (over 10 Kb) is being tagmented by the same bead at least twice and up to 5 times, out of about 15000 tagmentation events. Since the fragments share an index, they are unlikely to be there by chance, but are coming from the same bead.
Table 4 shows the number of reads (up to 5) within a 10 kb windows sharing an index.
Following transposition, the reaction mixture comprising CPT-DNA and free transposomes were subjected to column chromatography using Sephacryl S-400 and Sephacryl S-200 size exclusion chromatography and shown in
Densities of capture probes A7 and B7 were optimized on 1 μm beads and the results were shown in
Transposomes were prepared by mixing transposons having A7′ and B7′ capture sequences, that are complementary to A7 and B7 capture sequences on beads, with hyperactive Tn5 transposase. High molecular weight genomic DNA is mixed with the transposomes to generate CPT-DNA. Separately, beads are prepared with immobilized oligonucleotides: P5-A7, P7-B7, or P5-A7+P7-B7, where P5 and P7 are primer binding sequences and A7 and B7 are capture sequences complementary to A7′ and B7′ sequences respectively. Beads comprising P5-A7 alone, P7-B7 alone, P5-A7+P7-B7, or a mixture of P5-A7 and P7-B7 beads are treated with CPT-DNA and ligase was added to the reaction mixture to determine the efficiency of the hybridization of the immobilized oligos to the transposed DNA. The results are shown in
Several sets of transposomes were prepared. In one set, hyperactive Tn5 transposase is mixed with transposon sequences Tnp1 with 5′ biotin to prepare transposome 1. In another set, Tnp2 having unique index2 with 5′ biotin to prepare a transposome 2. In another set, hyperactive Tn5 transposase is mixed with transposon sequences Tnp3 with 5′ biotin to for transposome 3. In another Tnp4 having unique index 4 and 5′-biotin to prepare a transposome 4. Each of transposome 1&2 and transposome 3&4 are mixed separately with streptavidin beads to generate bead set 1 and bead set 2. The two set of beads are then mixed together and incubated with genomic DNA and tagmentation buffer to promote tagmentation of the genomic DNA. This is then followed by PCR amplification of the tagmented sequences. The amplified DNA is sequenced to analyze the insertion of the index sequences. If tagmentation is confined to the beads, majority of fragments will be coded with Tnp1/Tnp2 and Tnp3/Tnp4 indexes. If there is intra-molecular hybridization, the fragments may be coded with Tnp1/Tnp4, Tnp2/Tnp3, Tnp1/Tnp3, and Tnp2/Tnp4 indexes. Sequencing results after 5 and 10 cycles of PCR were shown in
Ninety six indexed transposome bead sets are prepared. Individual indexed transposomes were prepared by mixing transposon comprising an oligonucleotide comprising a Tn5 mosaic end sequence (ME) at the 5′-end and index sequence. Individually indexed transposomes were immobilized on beads through streptavidin-biotin interaction. Transposomes on beads were washed and all 96 individually indexed transposomes on beads were pooled. Oligonucleotides complimentary to the ME sequence and comprising an index sequence is annealed to the immobilized oligonucleotide creating transposons with unique indexes. The ninety six clonal indexed transposome bead sets are combined and incubated with high molecular weight (HMW) genomic DNA in presence of Nextera tagmentation buffer in a single tube.
The beads are washed and the transposase are removed by treating the reaction mixture with 0.1% SDS. The tagmented DNA is amplified with indexed primers and sequenced with PE HiSeq flow cell v2 using TruSeq v3 cluster kit and sequencing data are analyzed.
Clusters or islands of reads are observed. A plot of the nearest neighbor distances between the reads for each sequence shows essentially to major peaks, one from within the cluster (proximal) and another from between clusters (distal). A schematic of the method and the results are shown in
Transposomes are first assembled in solution by mixing a first oligonucleotide having ME′ sequence, a second oligonucleotide having ME-barcode-P5/P7 sequence, and Tn5 transposase. In first set, the first oligonucleotide having ME′ sequence is biotinylated at the 3′-end. In second case the oligonucleotide having ME-barcode-P5/P7 sequence is biotinylated at the 5′-end. To various concentrations (10 nM, 50 nM, and 200 NM) of each of the resulting transposome sets streptavidin beads are added such that the transposomes are immobilized on the streptavidin beads. The beads are washed and HMW genomic DNA is added and tagmentation is carried out. In some cases, the tagmented DNA is treated with 0.1% SDS and in other cases the tagmented DNA are untreated. The tagmented DNA is PCR amplified for 5-8 cycles and sequenced. The schematic is shown in
As shown in
Various amounts of target HMW DNA was added to clonally indexed beads with 50 mM Tn5: Transposon density and incubated for 15 or 60 min at 37 degree C. or for 60 min at room temperature. The transposomes comprised oligonucleotides with 3′-biotin. The tagmentation was carried out, the reaction mixture was treated with 0.1% SDS, and PCR amplified. The amplified DNA was sequenced.
Island size and distribution using solution based and bead based methods are compared. In a solution based approach, 96 transposomes each with unique index in the transposons are assembled in a 96 well plate. HMW genomic DNA is added, and the tagmentation reaction is carried out. The reaction product is treated with 0.1% SDS and PCR amplified. The amplified products were sequenced.
In a bead based approach, 96 transposomes each with unique index in the transposons are assembled in a 96 well plate. The oligonucleotides comprised 3′-end biotin. Streptavidin beads are added to each of the 96 well plate and incubated such that the transposomes are immobilized on the streptavidin beads. The beads are individually washed and pooled, HMW genomic DNA is added, and the tagmentation reaction is carried out in a single reaction vessel (one pot). The reaction product is treated with 0.1% SDS and PCR amplified. The amplified products were sequenced.
In the negative control, all 96 transposon sequences, each with unique index, are mixed together first. The oligonucleotides comprised 3′-end biotin. Transposomes are prepared from the individually mixed indexed transposons. Streptavidin beads are added to the mixture. HMW genomic DNA is added, and the tagmentation reaction is carried out. The reaction product is treated with 0.1% SDS and PCR amplified. The amplified products were sequenced.
The number of intra island reads is plotted versus the island size. The results as shown in
Detection of 60 kb Heterozygous Deletion
The sequencing data are extracted as fastq files and go through the demultiplexing process to generate individual fastq file for each barcode. The fastq files from the CPT sequencing are demultiplexed according to their indexes and aligned to the reference genome with the duplicates removed. The chromosomes are scanned by 5 kb/1 kb window, in which the number of the indexes showing any reads within the scanning window is recorded. Statistically for heterozygous deletion region only half amount of DNA is available for the library generation compared to its neighboring regions, therefore the number of indexes should be roughly half as its neighbors' as well. The NA12878 chr1 60 kb heterozygous deletion are shown in
Detection of Gene Fusion
The fastq files from the CPT sequencing are demultiplexed according to their index and aligned to the reference genome with the duplicates removed. Chromosomes are scanned in 2 kb window. Each 2 kb window is a 36864 vector in which each element records how many reads from a unique index have been found in this 2 kb window. For every 2 kb window pair (X,Y) across the genome, the weighted-Jaccard index is calculated. This index indicates the de facto distance between (X,Y) in the sample. Those indexes are displayed as the heatmap shown in
Detection of Deletions
The fastq files from the CPT sequencing are demultiplexed according to their index and aligned to the reference genome with the duplicates removed. Chromosomes are scanned in 1 kb window.
Bisulfite Conversion Efficiency Optimization
Conversion was assessed at the ME (mosaic element region) and gDNA region for index linked CPT-Seq libraries on beads. Promega's MethylEdge Bisulfite Conversion system was optimized to improve efficiency.
ME sequences were analyzed to determine efficiency of bisulfite conversion treatments and shown in
Expected sequencing read structure after sequencing BSC converted CPT-seq on beads libraries observed. Percent base metrics displayed with the IVC plot in
Whole-genome indexed linked CPT-seq libraries were enriched.
Enrichment statistics for HLA region is shown below:
To evaluate the exchange of the mosaic ends (ME) of transposome complexes, bead with different indices were prepared. After mixing, index exchange was determined by sequencing the libraries and reporting the indices for each library. % “swapped” was calculated as (D4+D5+E3+E5+f4)/(sum of all 96) and shown in the
Streptavidin magnetic beads were loaded with 1×, 6×, and 12× concentrations of TsTn5 transposome complex. The Epi-CPT seq protocol was performed for each bead type. The final PCR product was loaded on the Agilent BioAnalyzer for analysis and shown in
After bisulfite conversion, DNA becomes damaged, resulting in loss of the common sequences (CS2) needed for PCR amplification. DNA fragments CPTSeq and Epi-CPTSeq (Me-CPTSeq) libraries were analyzed by BioAnalyzer. Due to DNA damage during bisulfite conversion, the Epi-CPTSeq library has 5-fold lower yield and a smaller library size distribution compared to the CPTSeq library as shown in
Feasibility of the DNA end-recovery by Terminal transferase (TdT) mediated ligation was tested. Briefly, 5 pmoles of ssDNA template was incubated with TdT (10/50 U), Attenuator/adapter duplex (0/15/25 pmoles) and DNA Ligase (0/10 U) were incubated for 15 μm at 37 C. DNA products of extension/ligation were analyzed on a TBE-Urea gel and the results were shown in
Feasibility of the DNA end-recovery by Terminal transferase (TdT) mediated ligation was tested for sodium bisulfate converted bead bound library and shown in
Results of Methyl-CPTSeq assay are presented in
Claims
1.-58. (canceled)
59. A composition comprising a plurality of solid supports,
- wherein each of the solid supports comprises a plurality of immobilized oligonucleotides; and
- wherein each of the oligonucleotides on a given solid support comprises (a) a complementary capture sequence that is capable of hybridizing to at least a portion of an adaptor sequence on a transposon; (b) a first barcode sequence, wherein the first barcode sequence from each solid support in the plurality of solid supports differs from the first barcode sequence from other solid supports in the plurality of solid supports; and (c) a primer binding site,
- wherein (a), (b) and (c) are the same for each of the oligonucleotides on the given solid support.
60. The composition of claim 59, wherein the immobilized oligonucleotides further comprise
- (a) a second barcode sequence, or
- (b) a second barcode sequence and a third barcode sequence.
61. The composition of claim 59, wherein the solid support is a microsphere, bead, or particle.
62. The composition of claim 59, wherein the primer binding site of each immobilized oligonucleotide is a P5 or P7 primer binding site or a complement thereof.
63. The composition of claim 59, wherein an immobilized oligonucleotide on a given solid support comprises a partially double-stranded structure and a partially single-stranded structure.
64. The composition of claim 63, wherein one strand of the double-stranded structure of the immobilized oligonucleotide is immobilized on the solid support and a second strand of the double-stranded structure is partially complementary to the strand immobilized on the solid support and comprises the complementary capture sequence.
65. The composition of claim 64, wherein the complementary capture sequence is hybridized to the adaptor sequence of the transposon.
66. The composition of claim 65, wherein the complementary capture sequence is ligated to the adaptor sequence of the transposon.
67. The composition of claim 65, wherein the transposon comprises a transferred strand and a non-transferred strand and the adaptor sequence is on the transferred strand of the transposon.
68. The composition of claim 67, further comprising a fragment of a target nucleic acid ligated to the 3′ end of the transferred strand of the transposon.
69. The composition of claim 67, further comprising a transposase bound to the transposon that forms a transposome complex.
70. The composition of claim 69, wherein the transposase is a Tn5 transposase.
71. The composition of claim 69, wherein the transposome complex is bound to a target nucleic acid.
72. The composition of claim 71, wherein the target nucleic acid comprises a plurality of fragments and the transposome complex is bound to at least one fragment of the target nucleic acid.
73. The composition of claim 72, wherein the transferred strand is inserted into the 5′ end of at least one strand of the fragments while maintaining the contiguity of the target nucleic acid.
74. The composition of claim 63, wherein an immobilized oligonucleotide on a given solid support within the plurality of solid supports comprises a single-stranded structure.
75. The composition of claim 74, wherein the immobilized oligonucleotide is immobilized to the solid support via a linker at the 5′ end of the oligonucleotide and the complementary capture sequence is at the 3′ end of the immobilized oligonucleotide.
76. The composition of claim 74, wherein the complementary capture sequence of the immobilized oligonucleotide on a given solid support is hybridized to the adaptor sequence of a transposon.
77. The composition of claim 76, wherein the complementary capture sequence is ligated to the adaptor sequence of the transposon.
78. The composition of claim 76, wherein the transposon comprises a transferred strand and a non-transferred strand and the adaptor sequence is on the non-transferred strand of the transposon.
79. The composition of claim 76, wherein the transposon comprises a transferred strand and a non-transferred strand and the adaptor sequence is on the transferred strand of the transposon.
80. The composition of any one of claim 74, further comprising a fragment of a target nucleic acid ligated to the 3′ end of the transferred strand of the transposon.
81. The composition of claim 76, further comprising a transposase bound to the transposon to form a transposome complex.
82. The composition of claim 81, wherein the transposase is Tn5 transposase.
83. The composition of claim 81, wherein the transposome complex is bound to a target nucleic acid.
84. The composition of claim 83, wherein the target nucleic acid comprises a plurality of fragments and the transposome complex is bound to at least one fragment of the target nucleic acid.
85. The composition of claim 84, wherein the transposon comprises a transferred strand and a non-transferred strand, and wherein the transferred strand is inserted into the 5′ end of at least one strand of the at least one fragment while maintaining the contiguity of the target nucleic acid.
86. The composition of claim 59, wherein the adaptor sequence of the transposon comprises a barcode sequence.
87. The composition of claim 59, wherein the adaptor sequence of the transposon comprises a primer binding site.
88. The composition claim 87, wherein the primer binding site of the transposon is a P5 or P7 primer binding site or a complement thereof.
Type: Application
Filed: Oct 29, 2018
Publication Date: Feb 14, 2019
Inventors: Frank J. Steemers (San Diego, CA), Kevin L. Gunderson (San Diego, CA), Fan Zhang (San Diego, CA), Jason Richard Betley (Essex), Niall Anthony Gormley (Essex), Wouter Meuleman (San Diego, CA), Jacqueline Weir (Essex), Avgousta Ioannou (Essex), Gareth Jenkins (Essex), Rosamond Jackson (Essex), Natalie Morrell (Essex), Dmitry K. Pokholok (San Diego, CA), Steven J. Norberg (San Diego, CA), Molly He (San Diego, CA), Amirali Kia (San Diego, CA), Igor Goryshin (Madison, WI), Rigo Pantoja (San Diego, CA)
Application Number: 16/173,202
