Patent Application
20230160004
The Sequence Listing XML associated with this application is provided electronically in XML file format and is hereby incorporated by reference into the specification. The name of the XML file containing the Sequence Listing XML is “NATE-025_C02US_SeqList.xml”. The XML file is 132,673 bytes, created on Jul. 25, 2022, and is being submitted electronically via USPTO Patent Center.
There are currently a variety of methods for nucleic acid sequencing, i.e., the process of determining the precise order of nucleotides within a nucleic acid molecule. Current methods require amplifying a nucleic acid enzymatically, e.g., PCR, and/or by cloning. Further enzymatic polymerizations are required to produce a detectable signal by a light detection means. Such amplification and polymerization steps are costly and/or time-consuming. Thus, there is a need in the art for a method of nucleic acid sequencing that is amplification- and enzyme-free. The present invention addresses these needs.
The present invention provides sequencing probes, methods, kits, and apparatuses that provide enzyme-free, amplification-free, and library-free nucleic acid sequencing that has long-read-lengths and with low error rate. Moreover, the methods, kits, and apparatuses have rapid sample-to-answer capability. These features are particularly useful for sequencing in a clinical setting.
Provided herein are sequencing probes comprising a target binding domain and a barcode domain. The target binding domain and the barcode domain may be operably linked, e.g., covalently linked. A sequencing probe optionally comprises a spacer between the target binding domain and the barcode domain. The spacer can be any polymer with appropriate mechanical properties, for example, a single- or double-stranded DNA spacer (of 1 to 100 nucleotides, e.g., 2 to 50 nucleotides). Non-limiting examples of double-stranded DNA spacers include the sequences covered by SEQ ID NO: 25 to SEQ ID NO: 29.
The target binding domain comprises at least four nucleotides (e.g., 4, 5, 6, 7, 8, 9, 10, 11, 12, or more) and is capable of binding a target nucleic acid (e.g., DNA, RNA, and PNA). The barcode domain comprises a synthetic backbone, the barcode domain having at least a first position which comprises one or more attachment regions. The barcode domain may have one, two, three, four, five, six, seven, eight, nine, ten, eleven, twelve, or more positions; each position having one or more (e.g., one to fifty) attachment regions; each attachment region comprises at least one (i.e., one to fifty, e.g., ten to thirty copies of a nucleic acid sequence(s)) capable of reversibly binding to a complementary nucleic acid molecule (RNA or DNA). Certain positions in a barcode domain may have more attachment regions than other positions; alternately, each position in a barcode domain has the same number of attachment regions. The nucleic acid sequence of a first attachment region determines the position and identity of a first nucleotide in the target nucleic acid that is bound by a first nucleotide of the target binding domain, whereas the nucleic acid sequence of a second attachment region determines the position and identity of a second nucleotide in the target nucleic acid that is bound by a second nucleotide of the target binding domain. Likewise, the nucleic acid sequence of a sixth attachment region determines the position and identity of a sixth nucleotide in the target nucleic acid that is bound by a sixth nucleotide of the target binding domain. In embodiments, the synthetic backbone comprises a polysaccharide, a polynucleotide (e.g., single or double stranded DNA or RNA), a peptide, a peptide nucleic acid, or a polypeptide. The number of nucleotides in a target binding domain equals to or is greater than (e.g., 1, 2, 3, 4, or more) the number of positions in the barcode domain. Each attachment region in a specific position of the barcode domain may include one copy of the same nucleic acid sequence and/or multiple copies of the same nucleic acid sequence. However, an attachment region will include a different nucleic acid sequence than an attachment region in a different position of the barcode domain, even when both attachment regions identify the same type of nucleotide, e.g., adenine, thymine, cytosine, guanine, uracil, and analogs thereof. An attachment region may be linked to a modified monomer, e.g., a modified nucleotide, in the synthetic backbone, thereby creating a branch relative to the backbone. An attachment region may be part of a synthetic backbone’s polynucleotide sequence. One or more attachment regions may be adjacent to at least one flanking single-stranded polynucleotide, that is, an attachment region may be operably linked to a 5′ flanking single-stranded polynucleotide and/or to a 3′ flanking single-stranded polynucleotide. An attachment region with or without one or two flanking single-stranded polynucleotides may be hybridized to a hybridizing nucleic acid molecule lacking a detectable label. A hybridizing nucleic acid molecule lacking a detectable label may be between about 4 and about 20 nucleotides in length, e.g., 12 nucleotides, or longer.
An attachment region may be bound by a complementary nucleic acid comprising a detectable label. Each complementary nucleic acid may comprise a detectable label.
Alternately, an attachment region may be bound by a complementary nucleic acid that is part of a reporter complex (comprising detectable labels). A complementary nucleic acid (either comprising a detectable label or of a reporter complex) may be between about 4 and about 20 nucleotides in length, e.g., about 8, 10, 12, and 14 nucleotides, or more. In a reporter complex, a complementary nucleic acid is linked (directly or indirectly) to a primary nucleic acid molecule. A complementary nucleic acid may be indirectly linked to a primary nucleic acid molecule via a single or double-stranded nucleic acid linker (e.g., a polynucleotide comprising 1 to 100 nucleotides). A primary nucleic acid is hybridized to one or more (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more) secondary nucleic acids. Each secondary nucleic acid is hybridized to one or more (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more) tertiary nucleic acids; the tertiary nucleic acids comprise one or more detectable labels. A or each secondary nucleic acid may comprise a region that does not hybridize to a primary nucleic acid molecule and does not hybridize to a tertiary nucleic acid molecule (an “extra-handle”); this region may be four or more (e.g., about 6 to about 40, e.g., about 8, 10, 12, and 14) nucleotides in length. The region that does not hybridize to a primary nucleic acid molecule and does not hybridize to a tertiary nucleic acid molecule may comprise the nucleotide sequence of the complementary nucleic acid molecule that is linked to the primary nucleic acid molecule. This region may be located near the end of the secondary nucleic acid distal to its end that hybridizes to the primary nucleic acid. By having “extra-handles” comprising the nucleotide sequence of the complementary nucleic acid, the likelihood and speed at which a reporter complex binds to a sequencing probe is greatly increased. In any embodiment or aspect of the present invention, when a reporter complex comprises “extra-handles”, the reporter complex can hybridize to a sequencing probe either via the reporter complex’s complementary nucleic acid or via the “extra-handle.” Thus, for example, the phrase “binding to the first attachment region .. . a first complementary nucleic acid molecule of a first reporter complex” would be understood according to its plain meaning and also understood to mean “binding to the first attachment region ... an ‘extra handle’ of a first reporter complex.”
In embodiments, the terms “barcode domain” and “synthetic backbone” are synonymous.
Provided herein is a method for sequencing a nucleic acid using a sequencing probe of the present invention. The method comprises steps of: (1) hybridizing at least one sequencing probe, of the present invention, to an target nucleic acid that is immobilized (e.g., at one, two, three, four, five, six, seven, eight, nine, ten or more positions) to a substrate; (2) binding to the first attachment region a first complementary nucleic acid molecule (RNA or DNA) which has a detectable label (e.g., a fluorescent label) or a first complementary nucleic acid molecule of a first reporter complex comprising detectable labels (e.g., fluorescent labels); (3) detecting the detectable label(s), and (4) identifying the position and identity of the first nucleotide in the immobilized target nucleic acid. Optionally, the immobilized target nucleic acid is elongated prior to being bound by the probe. The method further comprises steps of: (5) contacting the first attachment region (with or without one or two flanking single-stranded polynucleotides) with a first hybridizing nucleic acid molecule lacking a detectable label, thereby unbinding the first complementary nucleic acid molecule having a detectable label or the first complementary nucleic acid molecule of a first reporter complex comprising detectable labels and binding to, at least, the first attachment region a first hybridizing nucleic acid lacking a detectable label; (6) binding to the second attachment region a second complementary nucleic acid molecule having a detectable label or a complementary nucleic acid molecule of a second reporter complex comprising detectable labels; (7) detecting the detectable label(s); and (8) identifying the position and identity of the second nucleotide in the immobilized target nucleic acid. Steps (5) to (8) are repeated until each nucleotide in the immobilized target nucleic acid and corresponding to the target binding domain has been identified. Steps (5) and (6) may occur concurrently or sequentially. Each (e.g., first, second, third, fourth, fifth, sixth, seventh, eighth, ninth, tenth, or higher) complementary nucleic acid molecule (having a detectable label or part of a reporter complex) has the same nucleic acid sequence as its corresponding (i.e., first, second, third, fourth, fifth, sixth, seventh, eighth, ninth, tenth, or higher) hybridizing nucleic acid molecule lacking a detectable label. The target nucleic acid is immobilized to a substrate by binding a first position and/or second position of the target nucleic acid with a first and/or a second capture probe; each capture probe comprises an affinity tag that selectively binds to a substrate. The first and/or second positions may be at or near a terminus of a target nucleic acid. The substrate can be any solid support known in the art, e.g., a coated slide and microfluidic device (e.g., coated with streptavidin). Other positons which are located distant from a terminus of a target nucleic acid may be selectively bound to the substrate. The nucleic acid may be elongated by applying a force (e.g., gravity, hydrodynamic force, electromagnetic force, flow-stretching, a receding meniscus technique, and combinations thereof) sufficient to extend the target nucleic acid.
Provided herein is a method for sequencing a nucleic acid using one population of probes of the present invention or a plurality of populations of probes of the present invention. The method comprises steps of: (1) hybridizing a first population of sequencing probes (of the present invention) to a target nucleic acid that is immobilized to a substrate (with each sequencing probe in the first population de-hybridizing from the immobilized target nucleic acid under about the same conditions, e.g., level of chaotropic agent, temperature, salt concentration, pH, and hydrodynamic force); (2) binding a plurality of first complementary nucleic acid molecules each having a detectable label or a plurality of first complementary nucleic acid molecules of a plurality of first reporter complexes each complex comprising detectable labels to a first attachment region in each sequencing probe in the first population; (3) detecting the detectable label(s); (4) identifying the position and identity of a plurality of first nucleotides in the immobilized target nucleic acid hybridized by sequencing probes in the first population; (5) contacting each first attachment region of each sequencing probe of the first population with a plurality of first hybridizing nucleic acid molecules lacking a detectable label thereby unbinding the first complementary nucleic acid molecules having a detectable label or of a reporter complex and binding to each first attachment region a first hybridizing nucleic acid molecule lacking a detectable label (6) binding a plurality of second complementary nucleic acid molecules each having a detectable label or a plurality of second complementary nucleic acid molecules of a plurality of second reporter complexes each complex comprising detectable labels to a second attachment region in each sequencing probe in the first population; (7) detecting the detectable label(s); and (8) identifying the position and identity of a plurality of second nucleotides in the immobilized target nucleic acid hybridized by sequencing probes in the first population. In step (9), steps (5) to (8) are repeated until each nucleotide in the immobilized target nucleic acid and corresponding to the target binding domain of each sequencing probe in the first population has been identified. Steps (5) and (6) may occur concurrently or sequentially. Thereby, the linear order of nucleotides is identified for regions of the immobilized target nucleic acid that were hybridized by the target binding domain of sequencing probes in the first population of sequencing probes.
In embodiments, when a plurality of populations (i.e., more than one population) of probes are used, the method further comprises steps of: (10) de-hybridizing each sequencing probe of the first population from the nucleic acid; (11) removing each de-hybridized sequencing probe of the first population; (12) hybridizing at least a second population of sequencing probes of the present invention, where each sequencing probe in the second population de-hybridizes from the immobilized target nucleic acid under about the same conditions and de-hybridizes from the immobilized target nucleic acid under different conditions from the sequencing probes in the first population; (13) binding a plurality of first complementary nucleic acid molecules each having a detectable label or a plurality of first complementary nucleic acid molecules of a plurality of first reporter complexes each complex comprising detectable labels to a first attachment region in each sequencing probe in the second population; (14) detecting the detectable label(s) (15) identifying the position and identity of a plurality of first nucleotides in the immobilized target nucleic acid hybridized by sequencing probes in the second population; (16) contacting each first attachment region of each sequencing probe of the second population with a plurality of first hybridizing nucleic acid molecules lacking a detectable label thereby unbinding the first complementary nucleic acid molecules (having a detectable label or from a reporter complex) and binding to each first attachment region a first hybridizing nucleic acid molecule lacking detectable label; (17) binding a plurality of second complementary nucleic acid molecules each having a detectable label or a plurality of second complementary nucleic acid molecules of a plurality of second reporter complexes each complex comprising detectable labels to a second attachment region in each sequencing probe in the second population; (18) detecting the detectable label(s); (19) identifying the position and identity of a plurality of second nucleotides in the immobilized target nucleic acid hybridized by sequencing probes in the second population; and (20) repeating steps (16) to (19) until the linear order of nucleotides has been identified for regions of the immobilized target nucleic acid that were hybridized by the target binding domain of sequencing probes in the second population of sequencing probes. Steps (16) and (17) may occur concurrently or sequentially.
Each sequencing probe in the second population may de-hybridize from the immobilized target nucleic acid at a different condition (e.g., a higher temperature, higher level of chaotropic agent, higher salt concentration, higher flow rate, and different pH) than the average condition for which the sequencing probes in the first population de-hybridize from the target nucleic acid.
However, when more than two populations of probes are used, then probes in two sequential populations may de-hybridize at different conditions and probes in non-sequential populations may de-hybridize at similar conditions. As an example, probes in a first population and third population may de-hybridize under similar conditions. In embodiments, sequential populations of probes de-hybridized at increasingly more stringent conditions (e.g., higher levels of chaotropic agent, salt concentration, and temperature). For a microfluidic device, using temperature as an example, a first population of probes may remain hybridized at a first temperature but de-hybridize at a second temperature, which is higher than the first. A second population of probes may remain hybridized at the second temperature but de-hybridize at a third temperature, which is higher than the second. In this example, solutions (comprising reagents required by the present method) flowing over a target nucleic acid for initial probe populations are at a lower temperature than solutions flowing over the target nucleic acid for later probe populations.
In some embodiments, after a population of probes has been used, the population of probes is de-hybridized from the target nucleic acid and a new aliquot of the same population of probes is used. For example, after a first population of probes has been hybridized, detected, and de-hybridized, a subsequent aliquot of the first population of probes is hybridized. Alternately, as an example, a first population of probes may be de-hybridized and replaced with a second population of probes; once the second population has been detected and de-hybridized, a subsequent aliquot of the first population of probes is hybridized to the target nucleic acid. Thus, a probe in the subsequent population may hybridize to a region of the target nucleic acid that had been previously sequenced (thereby gaining duplicative and/or confirmatory sequence information) or a probe in the subsequent population may hybridize to a region of the target nucleic acid that had not previously been sequenced (thereby gaining new sequence information). Accordingly, a population of probes may be re-aliquoted when a prior read was unsatisfactory (for any reason) and/or to improve the accuracy of the alignment resulting from the sequencing reads.
The probes hybridizing and de-hybridizing under similar conditions may have similar lengths of their target binding domain, GC content, or frequency of repeated bases and combinations thereof. Relationships between Tm and length of an oligonucleotide are taught, for example, in Sugimoto et al., Biochemistry, 34, 11211-6.
When more than two populations of probes are used, steps, as described for the first and second populations of sequencing probes, are repeated with additional populations of probes (e.g., 10 to 100 to 1000 populations). The number of populations of probes used will depend on a variety of factors, including but not limited to the size of the target nucleic acid, the number of unique probes in each population, the degree of overlap among sequencing probes desired, and the enrichment of probes to regions of interest.
A population of probes may contain extra sequencing probes directed to a specific region of interest in a target nucleic acid, e.g., a region containing a mutation (e.g., a point mutation) or a SNP allele. A population of probes may contain fewer sequencing probes directed to a specific region of less interest in a target nucleic acid.
A population of sequencing probes may be compartmentalized into discrete smaller pools of sequencing probes. The compartmentalization may be based upon predicted melting temperature of the target binding domain in the sequencing probes and/or upon sequence motif of the target binding domain in the sequencing probes. The compartmentalization may be based on empirically-derived rules. The different pools of sequencing probes can be reacted with the target nucleic acid using different reaction conditions, e.g., based on temperature, salt concentration, and/or buffer content. The compartmentalization may be performed to cover target nucleic acid with uniform coverage. The compartmentalization may be performed to cover target nucleic acid with known coverage profile.
The lengths of target binding domains in a population of sequencing probes may be reduced to increase coverage of probes in a specific region of a target nucleic acid. The lengths of target binding domains in a population of sequencing probes may be increased to decrease coverage of probes in a specific region of a target nucleic acid, e.g., to above the resolution limit of the sequencing apparatus.
Alternately or additionally, the concentration of sequencing probes in a population may be increased to increase coverage of probes in a specific region of a target nucleic acid. The concentration of sequencing probes may be reduced to decrease coverage of probes in a specific region of a target nucleic acid, e.g., to above the resolution limit of the sequencing apparatus.
The methods for sequencing a nucleic acid further comprises steps of assembling each identified linear order of nucleotides for each region of the immobilized target nucleic acid, thereby identifying a sequence for the immobilized target nucleic acid. Steps of assembling use a non-transitory computer-readable storage medium with an executable program stored thereon which instructs a microprocessor to arrange each identified linear order of nucleotides, thereby obtaining the sequence of the nucleic acid. Assembling can occur in “real time”, i.e., while data is being collected from sequencing probes rather than after all data has been collected.
The target nucleic acid, i.e., that is sequenced, may be between about 4 and 1,000,000 nucleotides in length. The target may include a whole, intact chromosome or a fragment thereof either of which is greater than 1,000,000 nucleotides in length.
Provided herein are apparatuses for performing a method of the present invention.
Provided herein are kits including sequencing probes of the present invention and for performing methods of the present invention. In embodiments, the kits include a substrate capable of immobilizing a nucleic acid via a capture probe, a plurality of sequencing probes of the present invention, at least one capture probe, at least one complementary nucleic acid molecule having a detectable label, at least one complementary nucleic acid molecule which lacks a detectable label, and instructions for use. In embodiments, the kit comprises about or at least 4096 unique sequencing probes. 4096 is the minimum number of unique probes necessary to include each possible hexameric combination (i.e., for probes each having six attachment regions in the barcode domains). Here, “4096” is achieved since there are four nucleotides options for six positions: 46. For a set of probes having four attachment regions in the barcode domains, only 256 (i.e., 44) unique probes will be needed. For a set of probes having eight nucleotides in their target binding domains, 48 (i.e., 65,536) unique probes will be needed. For a set of probes having ten nucleotides in their target binding domains, 410 (i.e., 1,048,576) unique probes will be needed.
In embodiments, the kit comprises about or at least twenty four distinct complementary nucleic acid molecule having a detectable label and about or at least twenty four distinct hybridizing nucleic acid molecule lacking a detectable label. A complementary nucleic acid may bind to an attachment region having a sequence of one of SEQ ID NO: 1 to 24, as non-limiting examples. Additional exemplary sequences that may be included in a barcode domain are listed in SEQ ID NO: 42 to SEQ ID NO: 81. Indeed, the nucleotide sequence is not limited; preferably it lacks substantial homology (e.g., 50% to 99.9%) with a known nucleotide sequence; this helps avoid undesirable hybridization of a complementary nucleic acid and a target nucleic acid.
Any of the above aspects and embodiments can be combined with any other aspect or embodiment.
Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. In the Specification, the singular forms also include the plural unless the context clearly dictates otherwise; as examples, the terms “a,” “an,” and “the” are understood to be singular or plural and the term “or” is understood to be inclusive. By way of example, “an element” means one or more element. Throughout the specification the word “comprising,” or variations such as “comprises” or “comprising,” will be understood to imply the inclusion of a stated element, integer or step, or group of elements, integers or steps, but not the exclusion of any other element, integer or step, or group of elements, integers or steps. About can be understood as within 10%, 9%, 8%, 7%, 6%, 5%, 4%, 3%, 2%, 1%, 0.5%, 0.1%, 0.05%, or 0.01% of the stated value. Unless otherwise clear from the context, all numerical values provided herein are modified by the term “about.”
Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. The references cited herein are not admitted to be prior art to the claimed invention. In the case of conflict, the present Specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and are not intended to be limiting. Other features and advantages of the invention will be apparent from the following detailed description and claim.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawings will be provided by the Office upon request and payment of the necessary fee.
The above and further features will be more clearly appreciated from the following detailed description when taken in conjunction with the accompanying drawings.
The present invention provides sequencing probes, methods, kits, and apparatuses that provide enzyme-free, amplification-free, and library-free nucleic acid sequencing that has long-read-lengths and with low error rate.
The present invention relates to a sequencing probe comprising a target binding domain and a barcode domain. Non-limiting examples of sequencing probes of the present invention are shown in
Another sequencing probe of the present invention is shown in
The target binding domain has at least four nucleotides, e.g., at least, 4, 5, 6, 7, 8, 9, 10, 11, 12, or more nucleotides. The target binding domain preferable is a polynucleotide. The target binding domain is capable of binding a target nucleic acid.
A probe may include multiple copies of the target binding domain operably linked to a synthetic backbone.
Probes can be designed to control the likelihood of hybridization and/or de-hybridization and the rates at which these occur. Generally, the lower a probe’s Tm, the faster and more likely that the probe will de-hybridize to/from a target nucleic acid. Thus, use of lower Tm probes will decrease the number of probes bound to a target nucleic acid.
The length of a target binding domain, in part, affects the likelihood of a probe hybridizing and remaining hybridized to a target nucleic acid. Generally, the longer (greater number of nucleotides) a target binding domain is, the less likely that a complementary sequence will be present in the target nucleotide. Conversely, the shorter a target binding domain is, the more likely that a complementary sequence will be present in the target nucleotide. For example, there is a 1/256 chance that a four-mer sequence will be located in a target nucleic acid versus a 1/4096 chance that a six-mer sequence will be located in the target nucleic acid. Consequently, a collection of shorter probes will likely bind in more locations for a given stretch of a nucleic acid when compared to a collection of longer probes.
In circumstances, it is preferable to have probes having shorter target binding domains to increase the number of reads in the given stretch of the nucleic acid, thereby enriching coverage of a target nucleic acid or a portion of the target nucleic acid, especially a portion of particular interest, e.g., when detecting a mutation or SNP allele.
However, it may be preferable to have fewer numbers of probes bound to a target nucleic acid since there are occasions when too many probes in a region may cause overlap of their detectable label, thereby preventing resolution of two nearby probes. This is explained as follows. Given that one nucleotide is 0.34 nm in length and given that the lateral (x-y) spatial resolution of a sequencing apparatus is about 200 nm, a sequencing apparatus’s resolution limit is about 588 base pair (i.e., a 1 nucleotide/0.34 nm × 200 nm). That is to say, the sequencing apparatus mentioned above would be unable to resolve signals from two probes hybridized to a target nucleic acid when the two probes are within about 588 base pair of each other. Thus, two probes, depending on the resolution of the sequencing apparatus, will need be spaced approximately 600 bp’s apart before their detectable label can be resolved as distinct “spots”. So, at optimal spacing, there should be a single probe per 600 bp of target nucleic-acid. A variety of software approaches (e.g., utilize fluorescence intensity values and wavelength dependent ratios) can be used to monitor, limit, and potentially deconvolve the number of probes hybridizing inside a resolvable region of a target nucleic acid and to design probe populations accordingly. Moreover, detectable labels (e.g., fluorescent labels) can be selected that provide more discrete signals. Furthermore, methods in the literature (e.g., Small and Parthasarthy: “Superresolution localization methods.” Annu. Rev. Phys Chem., 2014; 65:107-25) describe structured-illumination and a variety of super-resolution approaches which decrease the resolution limit of a sequencing microscope up to 10’s-of-nanometers. Use of higher resolution sequencing apparatuses allow for use of probes with shorter target binding domains.
As mentioned above, designing the Tm of probes can affect the number of probes hybridized to a target nucleic acid. Alternately or additionally, the concentration of sequencing probes in a population may be increased to increase coverage of probes in a specific region of a target nucleic acid. The concentration of sequencing probes may be reduced to decrease coverage of probes in a specific region of a target nucleic acid, e.g., to above the resolution limit of the sequencing apparatus.
The term “target nucleic acid” shall mean a nucleic acid molecule (DNA, RNA, or PNA) whose sequence is to be determined by the probes, methods, and apparatuses of the invention. In general, the terms “target nucleic acid”, “nucleic acid molecule,”, “nucleic acid sequence,” “nucleic acid”, “nucleic acid fragment,” “oligonucleotide” and “polynucleotide” are used interchangeably and are intended to include, but not limited to, a polymeric form of nucleotides that may have various lengths, either deoxyribonucleotides or ribonucleotides, or analogs thereof. Non-limiting examples of nucleic acids include a gene, a gene fragment, an exon, an intron, intergenic DNA (including, without limitation, heterochromatic DNA), messenger RNA (mRNA), transfer RNA, ribosomal RNA, ribozymes, small interfering RNA (siRNA), non-coding RNA (ncRNA), cDNA, recombinant polynucleotides, branched polynucleotides, plasmids, vectors, isolated DNA of a sequence, isolated RNA of a sequence, nucleic acid probes, and primers.
The present methods directly sequence a nucleic acid molecule obtained from a sample, e.g., a sample from an organism, and, preferably, without a conversion (or amplification) step. As an example, for RNA-based sequencing, the present methods do not require conversion of an RNA molecule to a DNA molecule (i.e., via synthesis of cDNA) before a sequence can be obtained. Since no amplification or conversion is required, a nucleic acid sequenced in the present invention will retain any unique base and/or epigenetic marker present in the nucleic acid when the nucleic acid is in the sample or when it was obtained from the sample. Such unique bases and/or epigenetic markers are lost in sequencing methods known in the art.
The target nucleic acid can be obtained from any sample or source of nucleic acid, e.g., any cell, tissue, or organism, in vitro, chemical synthesizer, and so forth. The target nucleic acid can be obtained by any art-recognized method. In embodiments, the nucleic acid is obtained from a blood sample of a clinical subject. The nucleic acid can be extracted, isolated, or purified from the source or samples using methods and kits well known in the art.
A nucleic acid molecule comprising the target nucleic acid may be fragmented by any means known in the art. Preferably, the fragmenting is performed by an enzymatic or a mechanical means. The mechanical means may be sonication or physical shearing. The enzymatic means may be performed by digestion with nucleases (e.g., Deoxyribonuclease I (DNase I)) or one or more restriction endonucleases.
When a nucleic acid molecule comprising the target nucleic acid is an intact chromosome, steps should be taken to avoid fragmenting the chromosome.
The target nucleic acid can include natural or non-natural nucleotides, comprising modified nucleotides, as well-known in the art.
Probes of the present invention may have overall lengths (including target binding domain, barcode domain, and any optional domains) of about 20 nanometers to about 50 nanometers. A probe’s backbone may a polynucleotide molecule comprising about 120 nucleotides.
The barcode domain comprises a synthetic backbone. The synthetic backbone and the target binding domain are operably linked, e.g., are covalently attached or attached via a linker. The synthetic backbone can comprise any material, e.g., polysaccharide, polynucleotide, polymer, plastic, fiber, peptide, peptide nucleic acid, or polypeptide. Preferably, the synthetic backbone is rigid. In embodiments, the backbone comprises “DNA origami” of six DNA double helices (See, e.g., Lin et al, “Submicrometre geometrically encoded fluorescent barcodes self-assembled from DNA.” Nature Chemistry; 2012 Oct; 4(10): 832-9). A barcode can be made of DNA origami tiles (Jungmann et al, “Multiplexed 3D cellular super-resolution imaging with DNA-PAINT and Exchange-PAINT”, Nature Methods, Vol. 11, No. 3, 2014).
The barcode domain comprises a plurality of positions, e.g., one, two, three, four, five, six, seven, eight, nine, ten, or more positions. The number of positions may be less than, equal to, or more than the number of nucleotides in the target binding domain. It is preferable to include additional nucleotides in a target binding domain than number of positions in the backbone domain, e.g., one, two, three, four, five, six, seven, eight, nine, ten, or more nucleotides. The length of the barcode domain is not limited as long as there is sufficient space for at least four positions, as described above.
Each position in the barcode domain corresponds to a nucleotide in the target binding domain and, thus, to a nucleotide in the target nucleic acid. As examples, the first position in the barcode domain corresponds to the first nucleotide in the target binding domain and the sixth position in the barcode domain corresponds to the sixth nucleotide in the target binding domain.
Each position in the barcode domain comprises at least one attachment region, e.g., one to 50, or more, attachment regions. Certain positions in a barcode domain may have more attachment regions than other positions (e.g., a first position may have three attachment regions whereas a second position may have two attachment positions); alternately, each position in a barcode domain has the same number of attachment regions. Each attachment region comprises at least one (i.e., one to fifty, e.g., ten to thirty) copies of a nucleic acid sequence(s) capable of being reversibly bound by a complementary nucleic acid molecule (e.g., DNA or RNA). In examples, the nucleic acid sequence in a first attachment region determines the position and identity of a first nucleotide in the target nucleic acid that is bound by a first nucleotide of the target binding domain. Each attachment region may be linked to a modified monomer (e.g., modified nucleotide) in the synthetic backbone such that the attachment region branches from the synthetic backbone. In embodiments, the attachment regions are integral to a polynucleotide backbone; that is to say, the backbone is a single polynucleotide and the attachment regions are parts of the single polynucleotide’s sequence. In embodiments, the terms “barcode domain” and “synthetic backbone” are synonymous.
The nucleic acid sequence in an attachment region identifies the position and identity of a nucleotide in the target nucleic acid that is bound by a nucleotide in the target binding domain of a sequencing probe. In a probe, each attachment region will have a unique overall sequence. Indeed, each position on a barcode domain can have an attachment region comprising a nucleic acid sequence that encodes one of four nucleotides, i.e., specific to one of adenine, thymine/uracil, cytosine, and guanine. Also, the attachment region of a first position (and encoding cytosine, for example) will include a nucleic acid sequence different from the attachment region of a second position (and encoding cytosine, for example). Thus, to a nucleic acid sequence in an attachment region in a first position that encodes a thymine, there will be no binding of a complementary nucleic acid molecule that identifies an adenine in a target nucleic acid corresponding to the first nucleotide of a target binding domain. Also, to an attachment region in a second position, there will be no binding of a complementary nucleic acid molecule that identifies an adenine in a target nucleic acid corresponding to the first nucleotide of a target binding domain.
Each position on a barcode domain may include one or more (up to fifty, preferably ten to thirty) attachment region; thus, each attachment region may bind one or more (up to fifty, preferably ten to thirty) complementary nucleic acid molecules. As examples, the probe in
Table 1 provides exemplary sequences, for illustration purposes only, for attachments regions for sequencing probes having up to six positions in its barcode domain and detectable labels on complementary nucleic acid that bind thereto.
As seen in Table 1, the nucleic acid sequence of a first attachment region may be one of SEQ ID NO: 1 to SEQ ID NO: 4 and the nucleic acid sequence of a second attachment may be one of SEQ ID NO: 5 to SEQ ID NO: 8. When the first nucleotide in the target nucleic acid is adenine, the nucleic acid sequence of the first attachment region would have the sequence of SEQ ID NO: 1 and when the second nucleotide in the target nucleic acid is adenine, the nucleic acid sequence of the second attachment region would have the sequence of SEQ ID NO: 5.
In embodiments, a complementary nucleic acid molecule may be bound by a detectable label. In alternate embodiments, a complementary nucleic acid is associated with a reporter complex comprising detectable labels.
The nucleotide sequence of a complementary nucleic acid is not limited; preferably it lacks substantial homology (e.g., 50% to 99.9%) with a known nucleotide sequence; this helps avoid undesirable hybridization of a complementary nucleic acid and a target nucleic acid.
An example of the reporter complex useful in the present invention is shown in
In embodiments, a primary nucleic acid molecule may comprise about 90 nucleotides. A secondary nucleic acid molecule may comprise about 87 nucleotides. A tertiary nucleic acid molecule may comprise about 15 nucleotides.
Reporter complexes can be of various designs. For example, a primary nucleic acid molecule can be hybridized to at least one (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more) secondary nucleic acid molecules. Each secondary nucleic acid molecule may be hybridized to at least one (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more) tertiary nucleic acid molecules. Exemplary reporter complexes are shown in
In embodiments, the reporter complexes are “pre-constructed”. That is, each polynucleotide in the complex is hybridized prior to contacting the complex with a probe. An exemplary recipe for pre-constructing five exemplary reporter complexes is shown in
A detectable moiety, label or reporter can be bound to a complementary nucleic acid or to a tertiary nucleic acid molecule in a variety of ways, including the direct or indirect attachment of a detectable moiety such as a fluorescent moiety, colorimetric moiety and the like. One of skill in the art can consult references directed to labeling nucleic acids. Examples of fluorescent moieties include, but are not limited to, yellow fluorescent protein (YFP), green fluorescent protein (GFP), cyan fluorescent protein (CFP), red fluorescent protein (RFP), umbelliferone, fluorescein, fluorescein isothiocyanate, rhodamine, dichlorotriazinylamine fluorescein, cyanines, dansyl chloride, phycocyanin, phycoerythrin and the like. Fluorescent labels and their attachment to nucleotides and/or oligonucleotides are described in many reviews, including Haugland, Handbook of Fluorescent Probes and Research Chemicals, Ninth Edition (Molecular Probes, Inc., Eugene, 2002); Keller and Manak, DNA Probes, 2nd Edition (Stockton Press, New York, 1993); Eckstein, editor, Oligonucleotides and Analogues: A Practical Approach (IRL Press, Oxford, 1991); and Wetmur, Critical Reviews in Biochemistry and Molecular Biology, 26:227-259 (1991). Particular methodologies applicable to the invention are disclosed in the following sample of references: U.S. Pat. Nos. 4,757,141; 5,151,507; and 5,091,519. In one aspect, one or more fluorescent dyes are used as labels for labeled target sequences, e.g., as disclosed by U.S. Pat. Nos. 5,188,934 (4,7-dichlorofluorescein dyes); 5,366,860 (spectrally resolvable rhodamine dyes); 5,847,162 (4,7-dichlororhodamine dyes); 4,318,846 (ether-substituted fluorescein dyes); 5,800,996 (energy transfer dyes); Lee et al. 5,066,580 (xanthine dyes); 5,688,648 (energy transfer dyes); and the like. Labelling can also be carried out with quantum dots, as disclosed in the following patents and patent publications: U.S. Pat. Nos. 6,322,901; 6,576,291; 6,423,551; 6,251,303; 6,319,426; 6,426,513; 6,444,143; 5,990,479; 6,207,392; 2002/0045045; and 2003/0017264. As used herein, the term “fluorescent label” comprises a signaling moiety that conveys information through the fluorescent absorption and/or emission properties of one or more molecules. Such fluorescent properties include fluorescence intensity, fluorescence lifetime, emission spectrum characteristics, energy transfer, and the like.
Commercially available fluorescent nucleotide analogues readily incorporated into nucleotide and/or oligonucleotide sequences include, but are not limited to, Cy3-dCTP, Cy3-dUTP, Cy5-dCTP, Cy5-dUTP (Amersham Biosciences, Piscataway, NJ), fluorescein- 12-dUTP, tetramethylrhodamine-6-dUTP, TEXAS RED™-5-dUTP, CASCADE BLUE™-7-dUTP, BODIPY TMFL-14-dUTP, BODIPY TMR-14-dUTP, BODIPY TMTR-14-dUTP, RHODAMINE GREEN™-5-dUTP, OREGON GREENR™ 488-5-dUTP, TEXAS RED™- 12-dUTP, BODIPY™ 630/650- 14-dUTP, BODIPY™ 650/665- 14-dUTP, ALEXA FLUOR™ 488-5-dUTP, ALEXA FLUOR™ 532-5-dUTP, ALEXA FLUOR™ 568-5-dUTP, ALEXA FLUOR™ 594-5-dUTP, ALEXA FLUOR™ 546- 14-dUTP, fluorescein- 12-UTP, tetramethylrhodamine-6-UTP, TEXAS RED™-5-UTP, mCherry, CASCADE BLUE™-7-UTP, BODIPY™ FL-14-UTP, BODIPY TMR-14-UTP, BODIPY™ TR-14-UTP, RHODAMINE GREEN™-5-UTP, ALEXA FLUOR™ 488-5-UTP, LEXA FLUOR™ 546- 14-UTP (Molecular Probes, Inc. Eugene, OR) and the like. Alternatively, the above fluorophores and those mentioned herein may be added during oligonucleotide synthesis using for example phosphoroamidite or NHS chemistry. Protocols are known in the art for custom synthesis of nucleotides having other fluorophores (See, Henegariu et al. (2000) Nature Biotechnol. 18:345). 2-Aminopurine is a fluorescent base that can be incorporated directly in the oligonucleotide sequence during its synthesis. Nucleic acid could also be stained, a priori, with an intercalating dye such as DAPI, YOYO- 1, ethidium bromide, cyanine dyes (e.g., SYBR Green) and the like.
Other fluorophores available for post-synthetic attachment include, but are not limited to, ALEXA FLUOR™ 350, ALEXA FLUOR™ 405, ALEXA FLUOR™ 430, ALEXA FLUOR™ 532, ALEXA FLUOR™ 546, ALEXA FLUOR™ 568, ALEXA FLUOR™ 594, ALEXA FLUOR™ 647, BODIPY 493/503, BODIPY FL, BODIPY R6G, BODIPY 530/550, BODIPY TMR, BODIPY 558/568, BODIPY 558/568, BODIPY 564/570, BODIPY 576/589, BODIPY 581/591, BODIPY TR, BODIPY 630/650, BODIPY 650/665, Cascade Blue, Cascade Yellow, Dansyl, lissamine rhodamine B, Marina Blue, Oregon Green 488, Oregon Green 514, Pacific Blue, Pacific Orange, rhodamine 6G, rhodamine green, rhodamine red, tetramethyl rhodamine, Texas Red (available from Molecular Probes, Inc., Eugene, OR), Cy2, Cy3, Cy3.5, Cy5, Cy5.5, Cy7 (Amersham Biosciences, Piscataway, NJ) and the like. FRET tandem fluorophores may also be used, including, but not limited to, PerCP-Cy5.5, PE-Cy5, PE-Cy5.5, PE-Cy7, PE-Texas Red, APC-Cy7, PE-Alexa dyes (610, 647, 680), APC-Alexa dyes and the like.
Metallic silver or gold particles may be used to enhance signal from fluorescently labeled nucleotide and/or oligonucleotide sequences (Lakowicz et al. (2003) BioTechniques 34:62).
Other suitable labels for an oligonucleotide sequence may include fluorescein (FAM, FITC), digoxigenin, dinitrophenol (DNP), dansyl, biotin, bromodeoxyuridine (BrdU), hexahistidine (6xHis), phosphor-amino acids (e.g., P-tyr, P-ser, P-thr) and the like. In one embodiment the following hapten/antibody pairs are used for detection, in which each of the antibodies is derivatized with a detectable label: biotin/a-biotin, digoxigenin/a-digoxigenin, dinitrophenol (DNP)/a-DNP, 5-Carboxyfluorescein (FAM)/a-FAM.
Detectable labels described herein are spectrally resolvable. “Spectrally resolvable” in reference to a plurality of fluorescent labels means that the fluorescent emission bands of the labels are sufficiently distinct, i.e., sufficiently non-overlapping, that molecular tags to which the respective labels are attached can be distinguished on the basis of the fluorescent signal generated by the respective labels by standard photodetection systems, e.g., employing a system of band pass filters and photomultiplier tubes, or the like, as exemplified by the systems described in U.S. Pat. Nos. 4,230,558; 4,811,218; or the like, or in Wheeless et al., pgs. 21-76, in Flow Cytometry: Instrumentation and Data Analysis (Academic Press, New York, 1985). In one aspect, spectrally resolvable organic dyes, such as fluorescein, rhodamine, and the like, means that wavelength emission maxima are spaced at least 20 nm apart, and in another aspect, at least 40 nm apart. In another aspect, chelated lanthanide compounds, quantum dots, and the like, spectrally resolvable means that wavelength emission maxima are spaced at least 10 nm apart, and in a further aspect, at least 15 nm apart.
The present invention relates to methods for sequencing a nucleic acid using a sequencing probe of the present invention. Examples of the method are shown in
The method comprises reversibly hybridizing at least one sequencing probe, of the present invention, to a target nucleic acid that is immobilized (e.g., at one, two, three, four, five, six, seven, eight, nine, ten, or more positions) to a substrate.
The substrate can be any solid support known in the art, e.g., a coated slide and a microfluidic device, which is capable of immobilizing a target nucleic acid. In certain embodiments, the substrate is a surface, membrane, bead, porous material, electrode or array. The target nucleic acid can be immobilized onto any substrate apparent to those of skill in the art.
In embodiments, the target nucleic acid is bound by a capture probe which comprises a domain that is complementary to a portion of the target nucleic acid. The portion may be an end of the target nucleic acid or not towards an end.
Exemplary useful substrates include those that comprise a binding moiety selected from the group consisting of ligands, antigens, carbohydrates, nucleic acids, receptors, lectins, and antibodies. The capture probe comprises a binding moiety capable of binding with the binding moiety of the substrate. Exemplary useful substrates comprising reactive moieties include, but are not limited to, surfaces comprising epoxy, aldehyde, gold, hydrazide, sulfhydryl, NHS-ester, amine, thiol, carboxylate, maleimide, hydroxymethyl phosphine, imidoester, isocyanate, hydroxyl, pentafluorophenyl-ester, psoralen, pyridyl disulfide or vinyl sulfone, polyethylene glycol (PEG), hydrogel, or mixtures thereof. Such surfaces can be obtained from commercial sources or prepared according to standard techniques. Exemplary useful substrates comprising reactive moieties include, but are not limited to, OptArray-DNA NHS group (Accler8), Nexterion Slide AL (Schott) and Nexterion Slide E (Schott).
In embodiments, the capture probe’s binding moiety is biotin and the substrate comprises avidin (e.g., streptavidin). Useful substrates comprising avidin are commercially available including TB0200 (Accelr8), SAD6, SAD20, SAD100, SAD500, SAD2000 (Xantec), SuperAvidin (Array-It), streptavidin slide (catalog #MPC 000, Xenopore) and STREPTAVIDINnslide (catalog #439003, Greiner Bio-one).
In embodiments, the capture probe’s binding moiety is avidin (e.g., streptavidin) and the substrate comprises biotin. Useful substrates comprising biotin that are commercially available include, but are not limited to, Optiarray-biotin (Accler8), BD6, BD20, BD100, BD500 and BD2000 (Xantec).
In embodiments, the capture probe’s binding moiety can comprise a reactive moiety that is capable of being bound to the substrate by photoactivation. The substrate could comprise the photoreactive moiety, or the first portion of the nanoreporter could comprise the photoreactive moiety. Some examples of photoreactive moieties include aryl azides, such as N((2-pyridyldithio)ethyl)-4-azidosalicylamide; fluorinated aryl azides, such as 4-azido-2,3,5,6-tetrafluorobenzoic acid; benzophenone-based reagents, such as the succinimidyl ester of 4-benzoylbenzoic acid; and 5-Bromo-deoxyuridine.
In embodiments, the capture probe’s binding moiety can be immobilized to the substrate via other binding pairs apparent to those of skill in the art.
After binding to the substrate, the target nucleic acid may be elongated by applying a force (e.g., gravity, hydrodynamic force, electromagnetic force “electrostretching”, flow-stretching, a receding meniscus technique, and combinations thereof) sufficient to extend the target nucleic acid.
The target nucleic acid may be bound by a second capture probe which comprises a domain that is complementary to a second portion of the target nucleic acid. The portion may be an end of the target nucleic acid or not towards an end. Binding of a second capture probe can occur after or during elongation of the target nucleic acid or to a target nucleic acid that has not been elongated. The second capture probe can have a binding as described above.
A capture probe may comprise or be associated with a detectable label, i.e., a fiducial spot.
The capture probe is capable of isolating a target nucleic acid from a sample. Here, a capture probe is added to a sample comprising the target nucleic acid. The capture probe binds the target nucleic acid via the region of the capture probe that his complementary to a region of the target nucleic acid. When the target nucleic acid contacts a substrate comprising a moiety that binds the capture probe’s binding moiety, the nucleic acid becomes immobilized onto the substrate.
To ensure that a user “captures” as many target nucleic acid molecules as possible from high fragmented samples, it is helpful to include a plurality of capture probes, each complementary to a different region of the target nucleic acid. For example, there may be three pools of capture probes, with a first pool complementary to regions of the target nucleic acid near its 5′ end, a second pool complementary to regions in the middle of the target nucleic acid, and a third pool near its 3′ end. This can be generalized to “n-regions-of-interest” per target nucleic acid. In this example, each individual pool of fragmented target nucleic acid bound to a capture probe comprising or bound to a biotin tag. 1/nth of input sample (where n = the number of distinct regions in target nucleic acid) is isolated for each pool chamber. The capture probe binds the target nucleic acid of interest. Then the target nucleic acid is immobilized, via the capture probe’s biotin, to an avidin molecule adhered to the substrate. Optionally, the target nucleic acid is stretched, e.g., via flow or electrostatic force. All n-pools can be stretched-and-bound simultaneously, or, in order to maximize the number of fully stretched molecules, pool 1 (which captures most 5′ region) can be stretched and bound first; then pool 2, (which captures the middle-of-target region) is then can be stretched and bound; finally, pool 3 is can be stretched and bound.
The number of distinct capture probes required is inversely related to the size of target nucleic acid fragment. In other word, more capture probes will be required for a highly-fragmented target nucleic acid. For sample types with highly fragmented and degraded target nucleic acids (e.g., Formalin-Fixed Paraffin Embedded Tissue) it may be useful to include multiple pools of capture probes. On the other hand, for samples with long target nucleic acid fragments, e.g., in vitro obtained isolated nucleic acids, a single capture probe at a 5′ end may be sufficient.
The region of the target nucleic acid between to two capture probes or after one capture probe and before a terminus of the target nucleic acid is referred herein as a “gap”. The gap is a portion of the target nucleic acid that is available to be bound by a sequencing probe of the present invention. The minimum gap is a target binding domain length (e.g., 4 to 10 nucleotides) and a maximum gap is the majority of a whole chromosome.
An immobilized target nucleic acid is shown in
The number of pools of complementary nucleic acids or reporter complexes is identical to the number of positions in the barcode domain. Thus, for a barcode domain having six positions, six pools will be cycled over the probes.
Alternately, prior to contacting a target nucleic acid with a probe, the probe may be hybridized at its first position to a complementary nucleic acid comprising a detectable label or a reporter complex. Thus, when contacted with its target nucleic acid, the probe is capable of emitting a detectable signal from its first position and it is unnecessary to provide a first pool of complementary nucleic acids or reporter complexes that are directed to the first position on the barcode domain.
In embodiments, the complementary nucleic acids comprising a detectable label or reporter complexes may be removed from the attachment region but not replaced with a hybridizing nucleic acid lacking a detectable label. This can occur, for example, by adding a chaotropic agent, increasing the temperature, changing salt concentration, adjusting pH, and/or applying a hydrodynamic force. In these embodiments fewer reagents (i.e., hybridizing nucleic acids lacking detectable labels) are needed.
In embodiments, the steps shown in
As mentioned above, complementary nucleic acids comprising detectable labels or reporter complexes can be removed from attachment regions but not replaced with hybridizing nucleic acid lacking detectable labels.
If needed, the rate of detectable label exchange can be accelerated by incorporating small single-stranded oligonucleotides that accelerate the rate of exchange of detectable labels (e.g., “Toe-Hold” Probes; see, e.g., Seeling et al., “Catalyzed Relaxation of a Metastable DNA Fuel”; J. Am. Chem. Soc. 2006, 128(37), pp12211-12220).
It is possible to replace the complementary nucleic acids or reporter complexes on a final position on a barcode domain (the sixth position in
Like
Alternately, hybridizing nucleic acid lacking a detectable label and complementary nucleic acids comprising detectable labels/reporter complexes are provided concurrently. This alternate embodiment is shown in
The distribution of probes along a length of target nucleic acid is critical for resolution of detectable signal. As discussed above, the resolution limit for two detectable labels is about 600 nucleotides. Preferably, each sequencing probe in a population of probes will bind no closer than 600 nucleotides from each other. As discussed above, 600 nucleotides is the resolution limit of a typical sequencing apparatus. In this case, a sequencing probe will provide a single read; this is shown in
Randomly, but in part depending on the length of the target binding domain, the Tm of the probes, and concentration of probes applied, it is possible for two distinct sequencing probes in a population to bind within 600 nucleotides of each other. In this case, unordered multiple reads will emit from a single resolution-limited spot; this is shown in
Alternately or additionally, the concentration of sequencing probes in a population may be reduced to decrease coverage of probes in a specific region of a target nucleic acid, e.g., to above the resolution limit of the sequencing apparatus, thereby producing a single read from a resolution-limited spot.
The method further comprising steps of assembling each identified linear order of nucleotides for each region of the immobilized target nucleic acid, thereby identifying a sequence for the immobilized target nucleic acid. The steps of assembling uses a non-transitory computer-readable storage medium with an executable program stored thereon. The program instructs a microprocessor to arrange each identified linear order of nucleotides for each region of the target nucleic acid, thereby obtaining the sequence of the nucleic acid. Assembling can occur in “real time”, i.e., while data is being collected from sequencing probes rather than after all data has been collected.
Any of the above aspects and embodiments can be combined with any other aspect or embodiment as disclosed here in the Summary and/or Detailed Description sections.
In certain exemplary embodiments, the terms “annealing” and “hybridization,” as used herein, are used interchangeably to mean the formation of a stable duplex. In one aspect, stable duplex means that a duplex structure is not destroyed by a stringent wash under conditions such as a temperature of either about 5° C. below or about 5° C. above the Tm of a strand of the duplex and low monovalent salt concentration, e.g., less than 0.2 M, or less than 0.1 M or salt concentrations known to those of skill in the art. The term “perfectly matched,” when used in reference to a duplex means that the polynucleotide and/or oligonucleotide strands making up the duplex form a double stranded structure with one another such that every nucleotide in each strand undergoes Watson-Crick base pairing with a nucleotide in the other strand. The term “duplex” comprises, but is not limited to, the pairing of nucleoside analogs, such as deoxyinosine, nucleosides with 2-aminopurine bases, PNAs, and the like, that may be employed. A “mismatch” in a duplex between two oligonucleotides means that a pair of nucleotides in the duplex fails to undergo Watson-Crick bonding.
As used herein, the term “hybridization conditions,” will typically include salt concentrations of less than about 1 M, more usually less than about 500 mM and even more usually less than about 200 mM. Hybridization temperatures can be as low as 5° C., but are typically greater than 22° C., more typically greater than about 30° C., and often in excess of about 37° C. Hybridizations are usually performed under stringent conditions, e.g., conditions under which a probe will specifically hybridize to its target subsequence. Stringent conditions are sequence-dependent and are different in different circumstances. Longer fragments may require higher hybridization temperatures for specific hybridization. As other factors may affect the stringency of hybridization, including base composition and length of the complementary strands, presence of organic solvents and extent of base mismatching, the combination of parameters is more important than the absolute measure of any one alone.
Generally, stringent conditions are selected to be about 5° C. lower than the Tm for the specific sequence at a defined ionic strength and pH. Exemplary stringent conditions include salt concentration of at least 0.01 M to no more than 1 M Na ion concentration (or other salts) at a pH 7.0 to 8.3 and a temperature of at least 25° C. For example, conditions of 5X SSPE (750 mM NaCl, 50 mM Na phosphate, 5 mM EDTA, pH 7.4) and a temperature of 25-30° C. are suitable for allele-specific probe hybridizations. For stringent conditions, see for example, Sambrook, Fritsche and Maniatis, “Molecular Cloning A Laboratory Manual, 2nd Ed.” Cold Spring Harbor Press (1989) and Anderson Nucleic Acid Hybridization, 1st Ed., BIOS Scientific Publishers Limited (1999). As used herein, the terms “hybridizing specifically to” or “specifically hybridizing to” or similar terms refer to the binding, duplexing, or hybridizing of a molecule substantially to a particular nucleotide sequence or sequences under stringent conditions.
Detectable labels associated with a particular position of a probe can be “readout” (e.g., its fluorescence detected) once or multiple times; a “readout” may be synonymous with the term “basecall”. Multiple reads improve accuracy. A target nucleic acid sequence is “read” when a contiguous stretch of sequence information derived from a single original target molecule is detected; typically, this is generated via multi-pass consensus (as defined below). As used herein, the term “coverage” or “depth of coverage” refers to the number of times a region of target has been sequenced (via discrete reads) and aligned to a reference sequence. Read coverage is the total number of reads that map to a specific reference target sequence; base coverage is the total number of basecalls made at a specific genomic position.
As used in herein, a “hybe and seq cycle” refers to all steps required to detect each attachment region on a particular probe or population of probes. For example, for a probe capable of detecting six positions on a target nucleic acid, one “hybe and seq cycle” will include, at least, hybridizing the probe to the target nucleic acid, hybridizing complementary nucleic acids/reporter complexes to attachment region at each of the six positions on the probe’s barcode domain, and detecting the detectable labels associated with each of the six positions.
The term “k-mer probe” is synonymous with a probe of the present invention.
When two or more sequences from discrete reads are aligned, the overlapping portions can be combined to create a single consensus sequence. In positions where overlapping portions have the same base (a single column of the alignment), those bases become the consensus. Various rules may be used to generate the consensus for positions where there are disagreements among overlapping sequences. A simple majority rule uses the most common base in the column as the consensus. A “multi-pass consensus” is an alignment of all discrete probe readouts from a single target molecule. Depending on the total number of cycles of probe populations/polls applied, each base position within a single target molecules can be queried with different levels of redundancy or overlap; generally, redundancy increases the confidence level of a basecall.
The “Raw Accuracy” is a measure of system’s inherent ability to correctly identify a base. Raw accuracy is dependent on sequencing technology. “Consensus Accuracy” is a measure of system’s ability to correctly identify a base with the use of additional reads and statistical power. “Specificity” refers to the percentage of reads that map to the intended targets out of total reads per run. “Uniformity” refers to the variability in sequence coverage across target regions; high uniformity correlates with low variability. This feature is commonly reported as the fraction of targeted regions covered by ≥20% of the average coverage depth across all targeted regions. Stochastic errors (i.e., intrinsic sequencing chemistry errors) can be readily corrected with ‘multi-pass’ sequencing of same target nucleic acid; given a sufficient number of passes, substantially ‘perfect consensus’ or ‘error-free’ sequencing can be achieved. The methods described herein may be implemented and/or the results recorded using any device capable of implementing the methods and/or recording the results. Examples of devices that may be used include but are not limited to electronic computational devices, including computers of all types. When the methods described herein are implemented and/or recorded in a computer, the computer program that may be used to configure the computer to carry out the steps of the methods may be contained in any computer readable medium capable of containing the computer program. Examples of computer readable medium that may be used include but are not limited to diskettes, CD-ROMs, DVDs, ROM, RAM, non-transitory computer-readable media, and other memory and computer storage devices. The computer program that may be used to configure the computer to carry out the steps of the methods, assemble sequence information, and/or record the results may also be provided over an electronic network, for example, over the internet, an intranet, or other network.
A “Consumable Sequencing Card” (
The sequencing data can be analyzed by any number of standard next-generation-sequencing assemblers (see, e.g., Wajid and Serpedin, “Review of general algorithmic features for genome assemblers for next generation sequencers” Genomics, proteomics & bioinformatics, 10 (2), 58-73, 2012). The sequencing data obtained within a single diffraction limited region of the microscope is “locally-assembled” to generate a consensus sequence from the multiple reads within a diffraction spot. The multiple diffraction spot assembled reads are then mapped together to generate contiguous sequences representing the entire targeted gene set, or a de-novo assembly of entire genome(s).
Additional teaching relevant to the present invention are described in one or more of the following: U.S. 8,148,512, U.S. 7,473,767, U.S. 7,919,237, U.S. 7,941,279, U.S. 8,415,102, U.S. 8,492,094, U.S. 8,519,115, U.S. 2009/0220978, U.S. 2009/0299640, U.S. 2010/0015607, U.S. 2010/0261026, U.S. 2011/0086774, U.S. 2011/0145176, U.S. 2011/0201515, U.S. 2011/0229888, U.S. 2013/0004482, U.S. 2013/0017971, U.S. 2013/0178372, U.S. 2013/0230851, U.S. 2013/0337444, U.S. 2013/0345161, U.S. 2014/0005067, U.S. 2014/0017688, U.S. 2014/0037620, U.S. 2014/0087959, U.S. 2014/0154681, and U.S. 2014/0162251, each of which is incorporated herein by reference in their entireties.
Below is described the timing for steps in the methods of the present invention and as shown in
The present invention requires minimal sample preparation. For example, as shown in
Calculations for an exemplary run are shown in
Binding a population of probes to an immobilized target nucleic acid takes about sixty seconds. This reaction can be accelerated by utilizing multiple copies of the target binding domain on the synthetic backbone. With microfluidic-controlled fluid exchange device, washing away unbound probes takes about a half a second.
Adding a first pool of complementary nucleic acids (comprising a detectable label) and binding them to attachment regions in the first position of the barcode domain takes about fifteen seconds.
Each field of view (FOV) is imaged for four different colors, each color representing a single-base. Fiducial spots placed on a 5′ capture probe or 3′ capture probe (or both) may be helpful for reading only those optical barcodes in-a-line (consistent with the presence of gapped target nucleic acid) between the two locations. Fiducial spots can also be added to each field of view in order to generate equal alignment of images upon successive steps in the sequencing process. All four images can be obtained at a single FOV and then the optical reading device may move to a new FOV, or take all FOV in one color then reimage in a second color. A single FOV can be read in about a half a second. It takes about a half a second to move to a next FOV. Therefore, the time to read “n” FOV’s equals “n” times 1 sec).
The complementary nucleic acids having detectable labels are removed from the first position of the barcode domain by addition of heat or washing with excess of complementary nucleic acids lacking detectable labels. If needed, the rate of detectable label exchange can be accelerated by incorporating small single-stranded oligonucleotides that accelerate the rate of exchange of detectable labels (e.g., “Toe-Hold” Probes; see, e.g., Seeling et al., “Catalyzed Relaxation of a Metastable DNA Fuel”; J. Am. Chem. Soc. 2006, 128(37), pp 12211-12220). A FOV can be reimaged to confirm that all complementary nucleic acids having detectable labels are removed before moving continuing. This takes about fifteen seconds. This step can be repeated until background signal levels are reached.
The above steps are repeated or the remaining positions in the probes’ barcode domain.
The total time to read equals m (bases read) times (15 sec + n FOVs times 1 sec + 15 sec). For example, when the number of positions in the barcode domain is 6 and 20 FOVs, the time to read equals 6 X (30 + 20 + 15) or 390 seconds.
Probes of the first population are de-hybridized. This takes about sixty seconds.
The above steps are repeated for second and subsequent populations of probes. If populations of sequencing probes are organized by melting temperature (Tm), each population of probes will require multiple hybridizations to ensure that each base is covered to required depth (this is driven by error rate). Moreover, by analyzing the hybridization reads during a run, it is possible to recognize each individual gene that is being sequenced well before the entire sequence is actually determined. Hence cycling can be repeated until a particular desired error-frequency (or coverage) is met.
Using the timing described above, together with some gapped-nucleic acid binding density estimates, throughput of a Nanostring (NSTG)-Next Generation Sequencer of the present invention can be estimated.
Net throughput of sequencer is given by: Fractional-Base-Occupancy X <gap-length> X number-of-gaps-per-FOV X number-of-bases-per-optical-barcode / [ 60 sec (hybridizing probes to target nucleic acid) + 0.5 sec (wash) + m: positions in the barcode domain X (15 sec (binding complementary nucleic acids) + nfovsX1 + 15 sec (unbinding complementary nucleic acids)) + 60 sec (de-hybridizing probes to target nucleic acid) ]
Therefore, in an example, a total “cycle” for a single gapped-nucleic acid (adding together from the method shown in
Assuming: 1% occupancy of the gapped-nucleic acid region, 4000 bases per gap, and 5000 gapped nucleic-acid fragments per FOV and an m of 6 and nFOVs of 20 (as described above) yields a net throughput of: 0.01X 4000 X 5000 X 20 = 4,000,000 6-base reads per 510.5 secs = 47,012 \.73 bases/sec.
Therefore, in this example, a net throughput per 24 hours of continuous measurement = 4.062 Gigabases (Gb) per day. Alternate estimates up to 12 Gb per day. See
As shown in
For the claimed invention, an error rate associated with sequencing is related to the free-energy difference between a fully-matched (m+n)-mer and a single-base mismatch (m-1+n)-mer. The sum of m+n is the number of nucleotides in a target binding domain and m represents the number of positions in a barcode domain. An estimate of the selectivity of hybridization can be made using the equation (See, Owczarzy, R. (2005), Biophys. Chem., 117:207-215 and Integrated DNA Technologies website: at the World Wide Web (www) idtdna.com/analyzer/Applications/Instructions/Default.aspx?AnalyzerDefinitions=true#Mismatc hMeltTemp):
where Ka is the association equilibrium constant obtained from predicted thermodynamic parameters,
Theta represents the percent bound of the exact complement and the single base mismatch sequences, which are expected to be annealed to target at the specified hybridization temperature. The T is the hybridization temperature in Kelvins, ΔH° (enthalpy) and ΔS° (entropy) are the melting parameters calculated from the sequence and the published nearest neighbor thermodynamic parameters, R is the ideal gas constant (1.987 cal·K-1mole-1), [strand1/2] is the molar concentration of an oligonucleotide, and the constant of -273.15 converts temperature from Kelvin to degrees of Celsius. The most accurate, nearest-neighbor parameters were obtained from the following publications for DNA/DNA base pairs (See, Allawi,H., SantaLucia, J. Biochemistry, 36, 10581), RNA/DNA base pairs (See, Sugimoto et al., Biochemistry, 34, 11211-6), RNA/RNA base pairs (See, Xia,T. et al., Biochemistry, 37, 14719),
As example of an estimate of the approximate error-rate expected from the NSTG-sequencer follows. For (m + n) equals 8’mer. Consider the following 8-mer barcode and its single-base mismatch.
(region to sequence)
(sequencing optical barcode with perfect match)
(sequencing optical barcode with single-base mismatch (G-T) pairing)
Using the IDT calculator based upon the above equations yields:
At 17.4° C. (the Tm of the perfect match case), (50% / 0.3%) would be the ratio of the correct optical barcode hybridized to that sequence versus the incorrect barcode at the Tm, yielding an estimated error rate for that sequence to be 0.6%.
A very high GC content sequencing calculation yields:
(region to sequence)
(sequencing optical barcode with perfect match)
(sequencing optical barcode with single-base mismatch (G-A) mis-pairing)
At 41.9° C. (the Tm of the perfect match case), (50% / 0.4%) would be the ratio of the correct optical barcode hybridized to that sequence versus the incorrect barcode at the Tm, yielding an estimated error rate for that sequence to be 0.8%.
Examination of a number of 8-mer pairs yields a distribution of error rates, in the range of 0.2% to 1%. While the above calculations will not be identical to the conditions used, these calculations provide an indication that the method of the present invention will have a relatively low intrinsic error rate, when compared to other single-molecule sequencing technologies, such as Pacific Biosciences and Oxford Nanopore Technologies where error rates can be significant (» 10%).
Additional experiments were performed using a target RNA hybridized with barcode and immobilized to the surface of cartridge using normal NanoString gene-expression binding technology (see, e.g., Geiss et al, “Direct multiplexed measurement of gene expression with color-coded probe pairs”; Nature Biotechnology, 26, 317 - 325 (2008)). The ability of a barcode with different target binding domain length and with a perfect match (YGBYGR-2 um optical bar code connected to perfect 10-mer match sequence) to hybridize to RNA-target was measured (
Ability of 10 mer to detect a SNP the real sequence is >15000 counts over background, whilst incorrect sequences are at most > 400 over background. In the presence of correct probe, error rates are expected to be <3% of real sequence. Note that this data is (in essence) a worse-case scenario. Having only a 10-base-pair hybridization sequence attached to a 6.6 Kilobase optical barcode reporter (Gen2 style). No specific condition optimizations were performed. This data, however, does reveal that the NanoString Next-Generation Sequencing approach is capable of resolving single-base pairs of sequence.
The detailed materials and methods utilized in the above study are as follows:
Hybridization Protocol Probe B plus codeset
Target Hybridization protocol (750 mM NaCl)
Materials
This application is a continuation of U.S. Pat. Application No. 16/528,760, filed Aug. 1, 2019, which is a continuation of U.S. Pat. Application No. 14/946,386, filed Nov. 19, 2015, which claims the benefit of U.S. Provisional Application No. 62/082,883, filed Nov. 21, 2014. The contents of each of the aforementioned patent application are incorporated herein by reference in their entireties.
| Number | Date | Country | |
|---|---|---|---|
| 62082883 | Nov 2014 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 16528760 | Aug 2019 | US |
| Child | 17816233 | US | |
| Parent | 14946386 | Nov 2015 | US |
| Child | 16528760 | US |
