US20150057948A1 - Analysis of measurements of a polymer - Google Patents
Analysis of measurements of a polymer Download PDFInfo
- Publication number
- US20150057948A1 US20150057948A1 US14/379,242 US201314379242A US2015057948A1 US 20150057948 A1 US20150057948 A1 US 20150057948A1 US 201314379242 A US201314379242 A US 201314379242A US 2015057948 A1 US2015057948 A1 US 2015057948A1
- Authority
- US
- United States
- Prior art keywords
- polymer
- feature vector
- measurements
- derived
- nanopore
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6858—Allele-specific amplification
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/48785—Electrical and electronic details of measuring devices for physical analysis of liquid biological material not specific to a particular test method, e.g. user interface or power supply
- G01N33/48792—Data management, e.g. communication with processing unit
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B82—NANOTECHNOLOGY
- B82Y—SPECIFIC USES OR APPLICATIONS OF NANOSTRUCTURES; MEASUREMENT OR ANALYSIS OF NANOSTRUCTURES; MANUFACTURE OR TREATMENT OF NANOSTRUCTURES
- B82Y15/00—Nanotechnology for interacting, sensing or actuating, e.g. quantum dots as markers in protein assays or molecular motors
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/483—Physical analysis of biological material
- G01N33/487—Physical analysis of biological material of liquid biological material
- G01N33/48707—Physical analysis of biological material of liquid biological material by electrical means
- G01N33/48721—Investigating individual macromolecules, e.g. by translocation through nanopores
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/68—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
- G01N33/6803—General methods of protein analysis not limited to specific proteins or families of proteins
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/68—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
- G01N33/6875—Nucleoproteins
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
Definitions
- the present invention relates generally to the field of analysing measurements of a polymer comprising polymer units, for example but without limitation a polynucleotide, made during translocation of the polymer through a nanopore.
- a nanopore measurement is typically made by restricting the flow of material between two pools of solution using a membrane.
- An aperture is provided within that membrane to allow the transfer of material from one pool of solution to another.
- the aperture has at least one dimension on the nanometre scale.
- measurements are made of that material.
- the most commonly used setup relies on the application of an applied potential to drive molecular species through the nanopore.
- An electrode is placed in each solution volume and the solution contains an electrolyte, typically a salt, such as 1 M NaCl.
- the applied potential across the electrodes also drives the electrolyte through the pore and generates a current.
- material passes through the pore it modifies the flow of ions which is directly observed in the current measurement.
- the degree of current block and the duration the material spends in the nanopore are indicative of its identity.
- Branton et al. disclose the use of a polymerase to control the speed of DNA translocation through the nanopore. This elegant solution has been adopted and adapted by many researchers in the field which has led to a number of publications. The basic concept is to provide a ratchet to the motion of the polymer, which could encompass a molecular motor or a molecular brake.
- Phi29 DNA polymerase could function on top of a nanopore (e.g. see Lieberman et al., J Am Chem Soc. 2010 Dec. 22; 132(50):17961-72, 61/402,903).
- the strength of the Phi29 DNAP binding to the template DNA was sufficient to allow multiple enzyme cycles to be performed on top of the nanopore, thus allowing the DNA to be pulled through the nanopore in a ratcheted fashion.
- the paper also revealed that Phi29 DNAP could be used to control DNA motion through the nanopore under conditions where the enzyme motion was inhibited.
- helicase families can be used to provide controlled movement of polynucleotides through a nanopore (e.g. see U.S. 61/549,998 (N115020), U.S. 61/581,332 (N115505), U.S. 61/581,340 (N115506)).
- Helicases have a number of properties that make them suitable for a nanopore system.
- An alternative method of slowing down translocation of a target single stranded DNA is to hybridise additional sections of ssDNA (hyb-DNA) along the length of the target strand.
- the target strand of DNA is rapidly fed through the pore under an applied potential. Once a double stranded section of the strand reaches the constriction of the nanopore, the translocation of the strand is halted, allowing the current to be read with the polymer at a fixed position.
- the hyb-DNA section is un-hybridised by the force of the applied field, and the target DNA strand continues to translocate the nanopore until another hyb-DNA is encountered. In this way, the current signatures for the DNA strand at a number of fixed positions are obtained.
- Derrington et al. propose a method of sequencing a strand of DNA using this approach.
- a measurement system may be created by placing a nanopore in an insulating membrane and measuring voltage-driven ionic transport through the nanopore in the presence of analyte molecules. The controlled movement of polymer through a nanopore results in a number of distinct levels of measurement that are indicative of the polymer sequence.
- the present invention makes use of an appreciation that many applications do not require the exact polymer sequence to be assigned. These include a significant number of diagnostic, clinical, scientific, genetic applications where the desired result can be obtained cheaply, quickly, and to a higher degree of accuracy without resorting to sequence information.
- the present invention involves derivation of a feature vector of time-ordered features representing characteristics of the measurements. Similarity between the derived feature vector and at least one other feature vector is then determined which provides information that is useful in many applications.
- the present invention does not require the assignment of polymer sequence, i.e. there is not necessarily a conversion of the measurement signal into sequence space.
- This provides useful analysis of the polymer in many applications, but reduces the burden on operation of the measurement system, because it is not necessary to resolve every single polymer unit in the sequence.
- This reduction on the constraints of the measurement system also increases the range of measurements systems. This may allow the use of a measurement system that is easier to design or operate, or may allow the use of a measurement system that is specifically adapted to analyse a particular characteristic of the polymer, even without being able to provide complete sequence information.
- An underlying feature of the invention is the conversion of the raw signal, that is the time-ordered series of measurements into a feature vector of time-ordered features.
- the series of measurements are derived as the polymer translocates through the nanopore and so provide information on the overall sequence, even if this is not complete.
- the derivation of the feature vector provides a representation which is also time-ordered but with a reduced data set.
- This feature vector may be thought of as a “signature” of the polymer.
- the feature vector is then compared to at least one other feature vector to determine the similarity.
- the at least one other feature vector may be, for example, a feature vector stored in a memory or another feature vector derived in the same manner. Based on the similarity, characteristics of the polymer may be derived.
- the step of deriving a feature vector may comprise identifying groups of consecutive measurements, and, in respect of each group, deriving values of one or more features that represent characteristics of the measurements of the group.
- the features may comprise: an average of the group of measurements; the period of the group of measurements; a variance of the group of measurements; the distribution of the group of measurements; or any combination thereof.
- the present invention is also applicable to signals with a lesser resolution, such that some k-mers may provide only a single measurement or no measurement at all.
- the derived feature vector may be compared with at least one other feature vector stored in a memory in respect of at least one class.
- the similarity may be determined between between the entirety or part of the derived feature vector and the entirety of the at least one other feature vector stored in the memory, or alternatively between the entirety or part of the derived feature vector and a part of the at least one other feature vector stored in the memory.
- the method may further comprise classifying the polymer from which the derived feature vector is derived as belonging to a said class on the basis of the determined similarity. This provides for identification of the polymer under investigation.
- the at least one other feature vector stored in the memory may be selected depending upon the polymer to be measured, or alternatively a library of plural other feature vectors stored in the memory may be used.
- a combined feature vector may be obtained from two or more feature vectors having overlapping regions wherein the similarity of the derived feature vector is determined between the combined feature vector.
- a non-overlapping region of the combined feature vector may be used to determine similarity between the derived feature vector, for example to identify a particular localised region of the derived feature vector.
- the method may be used to determine similarity between continuous or non-continuous regions of a derived feature vector and one or more feature vectors.
- plural parts of the derived feature vector may be compared all, parts or plural parts of stored feature vectors.
- the derived feature vector may be compared with at least one other feature vector that is a feature vector derived using the same method. This provides for identification of characteristics of plural polymers that are under investigation, relative to each other.
- the method may further comprise identifying clusters of similar feature vectors as a class and classifying the polymers from which the feature vectors are derived as belonging to an identified class.
- the method may further comprise identifying feature vectors that are derived from polymers that are fragments of a common polymer on the basis of similarity in overlapping parts of the feature vectors.
- the method may further comprise counting the numbers of feature vectors belonging to different classes. This provides for analysis of a population of polymers under investigation.
- the method may further comprise identifying localized regions where the derived feature vector is dissimilar to a feature vector in respect of the class in which the polymer is classified as belonging.
- the derived feature vector may be compared to a feature vector stored in a memory and the determination of similarity comprises determining localized regions where the derived feature vector is dissimilar to the at least one other feature vector stored in the memory.
- Such identification of localized regions where the derived feature vector is dissimilar to what is expected provides an analysis technique that is very powerful in many applications where change in relatively small regions of long sequences of polymers is significant.
- One example of such a technique is to identify mutations in a polymer that is a polynucleotide.
- the method may be performed on a series of measurements that has been previously made.
- the method may further comprise: translocating the polymer through a nanopore; and making the continuous series of measurements of the polymer.
- the method of analysing the series of measurements may be used in a method of estimating the presence, absence or amount of a target polymer based on the analysis.
- the polymer may comprise a mixture of two or more polymers and the relative amounts of one or more polymers may be determined.
- the method of estimating the presence, absence or amount of a target polymer may be applied to a polymer analyte in a method comprising: fragmenting the polymer analyte into polymers; and performing the method of estimating on the fragmented polymers.
- the polymer is a polynucleotide, and the polymer units are nucleotides
- the polymer analyte may be fragmented by a restriction enzyme.
- the method of analysing the series of measurements may be applied in a method of determining an alteration in a polymer, comprising: translocating a polymer through a nanopore repeatedly over a period of time; during each translocation, making a continuous series of measurements of the polymer; analysing each series of measurements.
- the step of determining similarity between the derived feature vector and at least one other feature vector may comprise either (a) determining similarity between the derived feature vector derived from each series of measurements and the same at least one other feature vector or (b) determining similarity between all the derived feature vectors derived from the series of measurements.
- the method may be used to determine the presence of a modified base or a point mutation.
- the methods may be used to guide a therapy or diagnosis or to identify an individual.
- the present invention has numerous applications. Some non-limitative examples or applications are as follows.
- This invention can be applied to single molecule label free detection systems for analysis of polymers, for example a nanopore system. It is common for such systems to comprise a recognition element that is influenced by more than one monomer units at a given polymer position. In these systems, extracting the relationship between measurement and polymer sequence may be challenging or resource demanding.
- This invention can be applied to any polymer analysis system where a polymer signature is indicative of a characteristic of that polymer and where the exact polymer sequence does not have to be known to determine said characteristic. Examples include but are not limited to: detection of single nucleotide polymorphisms (SNPs), presence or absence of specific sequences, grouping and counting of polymer sequences, design of labels and biomarkers, and identification of modified or damaged DNA.
- SNPs single nucleotide polymorphisms
- the method may be used for example to determine the presence, absence or amount of a target polymer analyte in a sample.
- the method may be used to measure an amount with respect to a threshold.
- the method may be used to determine the relative amounts of one or more target polymers in a mixture of polymers.
- the method may be used to guide a therapy or diagnosis based upon analysis of a single sample. Alternatively the method may be carried out plural times over a period for example to monitor progression of a disease or improvement of an individual. The method may be used to monitor an efficacy of treatment, for example where used as a theranostic.
- the method may be used in forensic applications for example to detect SNPs in mitochondrial DNA for DNA profiling of individuals, for genetic fingerprinting of individuals, for example by determining the presence of short tandem repeats, variable tandem repeats and the like.
- All the methods may be performed without estimating the sequence of polymer units of the polymer.
- FIG. 1 is a schematic diagram of a measurement system comprising a nanopore
- FIG. 2 is a plot of a signal of an event measured over time by a measurement system
- FIG. 3 is a graph of the frequency distributions of measurements of two different polynucleotides in a measurement system comprising a nanopore;
- FIGS. 4 and 5 are plots of 64 3-mer coefficients and 1024 5-mer coefficients, respectively, against predicted values from a first order linear model applied to sets of experimentally derived current measurements;
- FIG. 6 is a flowchart of a method of analyzing an input signal comprising measurements of a polymer
- FIG. 7 is a flowchart of a state detection step of FIG. 6 ;
- FIGS. 8 and 9 are plots, respectively, of an input signal subject to the state detection step and of the resultant series of measurements
- FIGS. 10 an 11 are flowcharts of examples of the similarity determination step of FIG. 6 ;
- FIG. 12 is plot of feature vectors for three fragments of a sequence identified by their overlap, for Example 2 of the method.
- FIG. 13 is a plot of similarity scores for candidate molecule as compared to all library sequences in Example 2.
- FIG. 14 is a plot of a candidate molecule aligned with best match library molecule in Example 2.
- FIG. 15 is a histogram of classification for 176 candidate molecules in Example 2.
- FIG. 16 is a graph of the feature vector in Example 3 of the method, illustrating the effect of SNPs on molecule 13;
- FIG. 17 is a histogram of classification for 176 candidate molecules with three SNPs in molecule 13 in Example 3;
- FIG. 18 is a graph of the alignment of a measured molecule with the library feature vector in Example 3.
- FIG. 19 is a plot of position-resolved differences between measurements and library feature vectors, illustrating position of SNPs, in Example 3.
- FIG. 20 is a plot of position-resolved differences between measurements and library feature vector without SNPs in Example 3;
- FIG. 21 is a plot of the final alignment of data with consensus landmarks in Example 4 of the method.
- FIG. 22 is a plot of position-resolved differences in candidate molecules 51-60 at approximately position 337 in Example 4.
- FIGS. 23 and 24 are diagrams of trees formed by neighbour joining on alignment similarity scores for a two cluster and a three cluster dataset respectively, in Example 5 of the method;
- FIGS. 25 to 27 are graphs of landmark consensus with final alignment of data for each identified cluster in Example 5;
- FIGS. 28 and 29 are histograms of classifications for the two cluster and three cluster experiment respectively in Example 5;
- FIG. 30 is a diagram of a tree formed by neighbour joining on alignment similarity scores in Example 6 of the method.
- FIG. 31 is a graph of landmark consensus with final alignment of data for each of three fragments in Example 6.
- Polymers that may be applied are as follows.
- the polymer may be a biological polymer.
- the polymer may be natural or synthetic.
- the polymer may be a polynucleotide (or nucleic acid), a polypeptide such as a protein, a polysaccharide, or any other polymer.
- the polymer units may be amino acids that are naturally occurring or synthetic.
- the polymer units may be monosaccharides.
- Polynucleotides that may be applied are as follows.
- a polynucleotide such as a nucleic acid
- the polynucleotide or nucleic acid may comprise any combination of any nucleotides.
- the nucleotides can be naturally occurring or artificial.
- One or more nucleotides in the target polynucleotide can be oxidized or methylated.
- One or more nucleotides in the target polynucleotide may be damaged.
- One or more nucleotides in the target polynucleotide may be modified, for instance with a label or a tag.
- the target polynucleotide may comprise one or more spacers.
- a nucleotide typically contains a nucleobase, a sugar and at least one phosphate group.
- the nucleobase is typically heterocyclic.
- Nucleobases include, but are not limited to, purines and pyrimidines and more specifically adenine, guanine, thymine, uracil and cytosine.
- the sugar is typically a pentose sugar.
- Nucleotide sugars include, but are not limited to, ribose and deoxyribose.
- the nucleotide is typically a ribonucleotide or deoxyribonucleotide.
- the nucleotide typically contains a monophosphate, diphosphate or triphosphate. Phosphates may be attached on the 5′ or 3′ side of a nucleotide.
- Nucleotides include, but are not limited to, adenosine monophosphate (AMP), adenosine diphosphate (ADP), adenosine triphosphate (ATP), guanosine monophosphate (GMP), guanosine diphosphate (GDP), guanosine triphosphate (GTP), thymidine monophosphate (TMP), thymidine diphosphate (TDP), thymidine triphosphate (TTP), uridine monophosphate (UMP), uridine diphosphate (UDP), uridine triphosphate (UTP), cytidine monophosphate (CMP), cytidine diphosphate (CDP), cytidine triphosphate (CTP), 5-methylcytidine monophosphate, 5-methylcytidine diphosphate, 5-methylcytidine triphosphate, 5-hydroxymethylcytidine monophosphate, 5-hydroxymethylcytidine diphosphate, 5-hydroxymethylcytidine triphosphate, cyclic
- the nucleotides are preferably selected from AMP, TMP, GMP, UMP, dAMP, dTMP, dGMP or dCMP.
- the nucleotides may be abasic (i.e. lack a nucleobase).
- the nucleotides may contain additional modifications.
- suitable modified nucleotides include, but are not limited to, 2′ amino pyrimidines (such as 2′-amino cytidine and 2′-amino uridine), 2′-hyrdroxyl purines (such as, 2′-fluoro pyrimidines (such as 2′-fluorocytidine and 2′ fluoro uridine), hydroxyl pyrimidines (such as 5′- ⁇ -P-borano uridine), 2′-O-methyl nucleotides (such as 2′-O-methyl adenosine, 2′-O-methyl guanosine, 2′-O-methyl cytidine and 2′-O-methyl uridine), 4′-thio pyrimidines (such as 4′-thio uridine and 4′-thio cytidine) and nucleotides have modifications of the nucleobase (such as 5-pentynyl-2′-deoxy uridine, 5-(3
- a nucleotide may be abasic (i.e. lack a nucleobase).
- the polynucleotide may be single stranded or double stranded.
- the polynucleotide may comprise one or more double stranded regions and one or more single regions.
- the polynucleotide can be a nucleic acid, such as deoxyribonucleic acid (DNA) or ribonucleic acid (RNA).
- the target polynucleotide can comprise one strand of RNA hybridized to one strand of DNA.
- the polynucleotide may be any synthetic nucleic acid known in the art, such as peptide nucleic acid (PNA), glycerol nucleic acid (GNA), threose nucleic acid (TNA), locked nucleic acid (LNA) or other synthetic polymers with nucleotide side chains.
- PNA peptide nucleic acid
- GNA glycerol nucleic acid
- TAA threose nucleic acid
- LNA locked nucleic acid
- the whole or only part of the target polynucleotide may be characterised using this method.
- the target polynucleotide can be any length.
- the polynucleotide can be at least 10, at least 50, at least 100, at least 150, at least 200, at least 250, at least 300, at least 400 or at least 500 nucleotide pairs in length.
- the polynucleotide can be 1000 or more nucleotide pairs, 5000 or more nucleotide pairs in length or 100000 or more nucleotide pairs in length.
- the target polynucleotide is present in any suitable sample.
- the invention is typically carried out on a sample that is known to contain or suspected to contain the target polynucleotide. Alternatively, the invention may be carried out on a sample to confirm the identity of one or more target polynucleotides whose presence in the sample is known or expected.
- the sample may be a biological sample.
- the invention may be carried out in vitro on a sample obtained from or extracted from any organism or microorganism.
- the organism or microorganism is typically archaean, prokaryotic or eukaryotic and typically belongs to one the five kingdoms: plantae, animalia, fungi, monera and protista.
- the invention may be carried out in vitro on a sample obtained from or extracted from any virus.
- the sample is preferably a fluid sample.
- the sample may be solid or semi-solid in origin which is subsequently treated to provide a fluid sample. Examples of such are faecal, skin, tissue, hair, bone and muscle.
- the sample typically comprises a body fluid of the patient.
- the sample may be chosen for example from urine, blood, plasma, serum, lymph, saliva, interstitial fluid, tears, mucus or amniotic fluid.
- the sample is human in origin, but alternatively it may be from another mammal animal such as from commercially farmed animals such as horses, cattle, sheep or pigs or may alternatively be pets such as cats or dogs.
- a sample of plant origin is typically obtained from a commercial crop, such as a cereal, legume, fruit or vegetable, for example wheat, barley, oats, canola, maize, soya, rice, bananas, apples, tomatoes, potatoes, grapes, tobacco, beans, lentils, sugar cane, cocoa, cotton.
- the sample may be a non-biological sample.
- the non-biological sample is preferably a fluid sample.
- Examples of a non-biological sample include surgical fluids, water such as drinking water, sea water or river water, and industrial samples such as reagents for laboratory tests, samples obtained from the synthesis of a polymer reagent.
- the sample is typically processed prior to being assayed, for example by centrifugation or by passage through a membrane that filters out unwanted molecules or cells, such as red blood cells.
- the sample may be measured immediately upon being taken.
- the sample may also be typically stored prior to assay, preferably below ⁇ 70° C.
- the sample may also be subject to any of the processes, designs, or modifications presented in U.S. 61/490,860.
- Membranes that may be used in a measurement system are as follows.
- the membrane is preferably an amphiphilic layer.
- An amphiphilic layer is a layer formed from amphiphilic molecules, such as phospholipids, which have both hydrophilic and lipophilic properties.
- the amphiphilic layer may be a monolayer or a bilayer.
- the membrane may be a co-block polymer such as disclosed by (Gonzalez-Perez et al., Langmuir, 2009, 25, 10447-10450).
- the membrane may be a lipid bilayer.
- Lipid bilayers are models of cell membranes and serve as excellent platforms for a range of experimental studies.
- lipid bilayers can be used for in vitro investigation of membrane proteins by single-channel recording.
- lipid bilayers can be used as biosensors to detect the presence of a range of substances.
- Suitable amphiphilic layers include, but are not limited to, a planar lipid bilayer, a supported bilayer or a liposome.
- the lipid bilayer is preferably a planar lipid bilayer.
- Suitable lipid bilayers are disclosed in International Application No. PCT/GB08/000,563 (published as WO 2008/102121), International Application No. PCT/GB08/004,127 (published as WO 2009/077734) and International Application No. PCT/GB2006/001057 (published as WO 2006/100484).
- Lipid bilayers are commonly formed by the method of Montal and Mueller (Proc. Natl. Acad. Sci. USA., 1972; 69: 3561-3566), in which a lipid monolayer is carried on aqueous solution/air interface past either side of an aperture which is perpendicular to that interface.
- Montal & Mueller The method of Montal & Mueller is popular because it is a cost-effective and relatively straightforward method of forming good quality lipid bilayers that are suitable for protein pore insertion.
- Other common methods of bilayer formation include tip-dipping, painting bilayers and patch-clamping of liposome bilayers.
- amphiphilic layer is formed as described in International Application No. PCT/GB08/004,127 (published as WO 2009/077734).
- the membrane is a solid state layer.
- a solid-state layer is not of biological origin.
- a solid state layer is not derived from or isolated from a biological environment such as an organism or cell, or a synthetically manufactured version of a biologically available structure.
- Solid state layers can be formed from both organic and inorganic materials including, but not limited to, microelectronic materials, insulating materials such as Si 3 N 4 , Al 2 O 3 , and SiO, organic and inorganic polymers such as polyamide, plastics such as Teflon® or elastomers such as two-component addition-cure silicone rubber, and glasses.
- the solid state layer may be formed from monatomic layers, such as graphene, or layers that are only a few atoms thick.
- Suitable graphene layers are disclosed in International Application No. PCT/US2008/010637 (published as WO 2009/035647).
- the solid state membrane can also support a nanopore derived from biological material, non-limiting examples have been disclosed by Hall et al. (Nat Nanotechnol. 2010 December; 5(12):874-7) and Bell et al. (Nano Lett. 2012 Jan. 11; 12(1):512-7), and International Application No. PCT/US2011/039621 (published as WO/2012/005857).
- the method is typically carried out using (i) an artificial amphiphilic layer comprising a pore, (ii) an isolated, naturally-occurring amphiphilic layer comprising a pore, or (iii) a cell having a pore inserted therein.
- the method is preferably carried out using an artificial amphiphilic layer.
- the bilayer may comprise other transmembrane and/or intramembrane proteins as well as other molecules in addition to the pore. Suitable apparatus and conditions are discussed below.
- the method of the invention is typically carried out in vitro. Nanopores that may be applied are as follows.
- the measurement system comprises a nanopore.
- the measurements are taken during translocation of the polymer through the nanopore.
- the translocation of the polymer through the nanopore generates a characteristic signal in the measured property that may be observed, and may be referred to overall as an “event”.
- the nanopore is a pore, typically having a size broadly speaking of the order of nanometres, that allows the passage of polymers therethrough.
- references to a “pore” mean a nanopore in this sense.
- the nanopore may be a biological pore or a solid state pore.
- a solid state pore is typically an aperture in a solid state layer.
- a solid state pore may be used in combination with additional components which provide an alternative or additional measurement of the polymer such as tunnelling electrodes (Ivanov A P et al., Nano Lett. 2011 Jan. 12; 11(1):279-85), or a field effect transistor (FET) device (International Application WO 2005/124888).
- Solid state pores may be formed by known processes including for example those described in WO 00/79257.
- the nanopore is preferably a transmembrane protein pore.
- a transmembrane protein pore is a polypeptide or a collection of polypeptides that permits hydrated ions to flow from one side of a membrane to the other side of the membrane.
- the transmembrane protein pore is capable of forming a pore that permits hydrated ions driven by an applied potential to flow from one side of the membrane to the other.
- the transmembrane protein pore allows a polymer, such as DNA or RNA, to be moved through the pore.
- the transmembrane protein pore may be a monomer or an oligomer.
- the pore is preferably made up of several repeating subunits, such as 6, 7 or 8 subunits.
- the pore is more preferably a heptameric or octameric pore.
- the transmembrane protein pore typically comprises a barrel or channel through which the ions may flow.
- the subunits of the pore typically surround a central axis and contribute strands to a transmembrane (3-barrel or channel or a transmembrane ⁇ -helix bundle or channel.
- the barrel or channel of the transmembrane protein pore typically comprises amino acids that facilitate interaction with analyte, such as polymers, nucleotides, polynucleotides or nucleic acids. These amino acids are preferably located near a constriction of the barrel or channel.
- the transmembrane protein pore typically comprises one or more positively charged amino acids, such as arginine, lysine or histidine, or aromatic amino acids, such as tyrosine or tryptophan. These amino acids typically facilitate the interaction between the pore and polymers, nucleotides, polynucleotides or nucleic acids.
- Transmembrane protein pores for use in accordance with the invention can be derived from ⁇ -barrel pores or ⁇ -helix bundle pores.
- ⁇ -barrel pores comprise a barrel or channel that is formed from ⁇ -strands.
- Suitable ⁇ -barrel pores include, but are not limited to, ⁇ -toxins, such as ⁇ -hemolysin, anthrax toxin and leukocidins, and outer membrane proteins/porins of bacteria, such as Mycobacterium smegmatis porin (Msp), for example MspA, outer membrane porin F (OmpF), outer membrane porin G (OmpG), outer membrane phospholipase A and Neisseria autotransporter lipoprotein (NalP).
- Msp Mycobacterium smegmatis porin
- OmpF outer membrane porin F
- OmpG outer membrane porin G
- Neisseria autotransporter lipoprotein Neisseria autotransporter lipoprotein
- ⁇ -helix bundle pores comprise a barrel or channel that is formed from ⁇ -helices.
- Suitable ⁇ -helix bundle pores include, but are not limited to, inner membrane proteins and outer membrane proteins, such as WZA and ClyA toxin.
- the transmembrane pore may be derived from Msp or from ⁇ -hemolysin ( ⁇ -HL).
- the transmembrane protein pore is preferably derived from Msp, preferably from MspA. Such a pore will be oligomeric and typically comprises 7, 8, 9 or 10 monomers derived from Msp.
- the pore may be a homo-oligomeric pore derived from Msp comprising identical monomers. Alternatively, the pore may be a hetero-oligomeric pore derived from Msp comprising at least one monomer that differs from the others.
- the pore is derived from MspA or a homolog or paralog thereof.
- a monomer derived from Msp comprises the sequence shown in SEQ ID NO: 2 or a variant thereof.
- SEQ ID NO: 2 is the MS-(B1)8 mutant of the MspA monomer. It includes the following mutations: D90N, D91N, D93N, D118R, D134R and E139K.
- a variant of SEQ ID NO: 2 is a polypeptide that has an amino acid sequence which varies from that of SEQ ID NO: 2 and which retains its ability to form a pore. The ability of a variant to form a pore can be assayed using any method known in the art.
- the variant may be inserted into a lipid bilayer along with other appropriate subunits and its ability to oligomerise to form a pore may be determined.
- Methods are known in the art for inserting subunits into membranes, such as lipid bilayers.
- subunits may be suspended in a purified form in a solution containing a lipid bilayer such that it diffuses to the lipid bilayer and is inserted by binding to the lipid bilayer and assembling into a functional state.
- subunits may be directly inserted into the membrane using the “pick and place” method described in M. A. Holden, H. Bayley. J. Am. Chem. Soc. 2005, 127, 6502-6503 and International Application No. PCT/GB2006/001057 (published as WO 2006/100484).
- a variant will preferably be at least 50% homologous to that sequence based on amino acid identity. More preferably, the variant may be at least 55%, at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 90% and more preferably at least 95%, 97% or 99% homologous based on amino acid identity to the amino acid sequence of SEQ ID NO: 2 over the entire sequence. There may be at least 80%, for example at least 85%, 90% or 95%, amino acid identity over a stretch of 100 or more, for example 125, 150, 175 or 200 or more, contiguous amino acids (“hard homology”).
- Standard methods in the art may be used to determine homology.
- the UWGCG Package provides the BESTFIT program which can be used to calculate homology, for example used on its default settings (Devereux et al (1984) Nucleic Acids Research 12, p38′7-395).
- the PILEUP and BLAST algorithms can be used to calculate homology or line up sequences (such as identifying equivalent residues or corresponding sequences (typically on their default settings)), for example as described in Altschul S. F. (1993) J Mol Evol 36:290-300; Altschul, S F et al (1990) J Mol Biol 215:403-10.
- Software for performing BLAST analyses is publicly available through the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/).
- SEQ ID NO: 2 is the MS-(B1) 8 mutant of the MspA monomer.
- the variant may comprise any of the mutations in the MspB, C or D monomers compared with MspA.
- the mature forms of MspB, C and D are shown in SEQ ID NOs: 15 to 17.
- the variant may comprise the following substitution present in MspB: A138P.
- the variant may comprise one or more of the following substitutions present in MspC: A96G, N 102 E and A138P.
- the variant may comprise one or more of the following mutations present in MspD: Deletion of G1, L2V, E5Q, L8V, D13G, W21A, D22E, K47T, 149H, 168V, D91G, A96Q, N102D, S103T, V104I, S136K and G141A.
- the variant may comprise combinations of one or more of the mutations and substitutions from Msp B, C and D.
- the variant may comprise the mutation L88N.
- the variant of SEQ ID NO: 2 has the mutation L88N in addition to all the mutations of MS-B1 and is called MS-B2.
- the pore used in the invention may be MS-(B2)8 or MS-(B2C)8.
- Amino acid substitutions may be made to the amino acid sequence of SEQ ID NO: 2 in addition to those discussed above, for example up to 1, 2, 3, 4, 5, 10, 20 or 30 substitutions.
- Conservative substitutions replace amino acids with other amino acids of similar chemical structure, similar chemical properties or similar side-chain volume.
- the amino acids introduced may have similar polarity, hydrophilicity, hydrophobicity, basicity, acidity, neutrality or charge to the amino acids they replace.
- the conservative substitution may introduce another amino acid that is aromatic or aliphatic in the place of a pre-existing aromatic or aliphatic amino acid.
- Conservative amino acid changes are well-known in the art and may be selected in accordance with the properties of the 20 main amino acids as defined in Table 2 below. Where amino acids have similar polarity, this can also be determined by reference to the hydropathy scale for amino acid side chains in Table 3.
- One or more amino acid residues of the amino acid sequence of SEQ ID NO: 2 may additionally be deleted from the polypeptides described above. Up to 1, 2, 3, 4, 5, 10, 20 or 30 residues may be deleted, or more.
- Variants may include fragments of SEQ ID NO: 2. Such fragments retain pore forming activity. Fragments may be at least 50, 100, 150 or 200 amino acids in length. Such fragments may be used to produce the pores. A fragment preferably comprises the pore forming domain of SEQ ID NO: 2. Fragments must include one of residues 88, 90, 91, 105, 118 and 134 of SEQ ID NO: 2. Typically, fragments include all of residues 88, 90, 91, 105, 118 and 134 of SEQ ID NO: 2.
- One or more amino acids may be alternatively or additionally added to the polypeptides described above.
- An extension may be provided at the amino terminal or carboxy terminal of the amino acid sequence of SEQ ID NO: 2 or polypeptide variant or fragment thereof.
- the extension may be quite short, for example from 1 to 10 amino acids in length. Alternatively, the extension may be longer, for example up to 50 or 100 amino acids.
- a carrier protein may be fused to an amino acid sequence according to the invention. Other fusion proteins are discussed in more detail below.
- a variant is a polypeptide that has an amino acid sequence which varies from that of SEQ ID NO: 2 and which retains its ability to form a pore.
- a variant typically contains the regions of SEQ ID NO: 2 that are responsible for pore formation. The pore forming ability of Msp, which contains a ⁇ -barrel, is provided by ⁇ -sheets in each subunit.
- a variant of SEQ ID NO: 2 typically comprises the regions in SEQ ID NO: 2 that form ⁇ -sheets.
- One or more modifications can be made to the regions of SEQ ID NO: 2 that form ⁇ -sheets as long as the resulting variant retains its ability to form a pore.
- a variant of SEQ ID NO: 2 preferably includes one or more modifications, such as substitutions, additions or deletions, within its ⁇ -helices and/or loop regions.
- the monomers derived from Msp may be modified to assist their identification or purification, for example by the addition of histidine residues (a hist tag), aspartic acid residues (an asp tag), a streptavidin tag or a flag tag, or by the addition of a signal sequence to promote their secretion from a cell where the polypeptide does not naturally contain such a sequence.
- An alternative to introducing a genetic tag is to chemically react a tag onto a native or engineered position on the pore. An example of this would be to react a gel-shift reagent to a cysteine engineered on the outside of the pore. This has been demonstrated as a method for separating hemolysin hetero-oligomers (Chem Biol. 1997 July; 4(7):497-505).
- the monomer derived from Msp may be labelled with a revealing label.
- the revealing label may be any suitable label which allows the pore to be detected. Suitable labels include, but are not limited to, fluorescent molecules, radioisotopes, e.g. 125I, 35S, enzymes, antibodies, antigens, polynucleotides and ligands such as biotin.
- the monomer derived from Msp may also be produced using D-amino acids.
- the monomer derived from Msp may comprise a mixture of L-amino acids and D-amino acids. This is conventional in the art for producing such proteins or peptides.
- the monomer derived from Msp contains one or more specific modifications to facilitate nucleotide discrimination.
- the monomer derived from Msp may also contain other non-specific modifications as long as they do not interfere with pore formation.
- a number of non-specific side chain modifications are known in the art and may be made to the side chains of the monomer derived from Msp. Such modifications include, for example, reductive alkylation of amino acids by reaction with an aldehyde followed by reduction with NaBH 4 , amidination with methylacetimidate or acylation with acetic anhydride.
- the monomer derived from Msp can be produced using standard methods known in the art.
- the monomer derived from Msp may be made synthetically or by recombinant means.
- the pore may be synthesized by in vitro translation and transcription (IVTT). Suitable methods for producing pores are discussed in International Application Nos. PCT/GB09/001,690 (published as WO 2010/004273), PCT/GB09/001,679 (published as WO 2010/004265) or PCT/GB10/000,133 (published as WO 2010/086603). Methods for inserting pores into membranes are discussed.
- the transmembrane protein pore is also preferably derived from ⁇ -hemolysin ( ⁇ -HL).
- ⁇ -HL ⁇ -hemolysin
- the wild type ⁇ -HL pore is formed of seven identical monomers or subunits (i.e. it is heptameric).
- the sequence of one monomer or subunit of ⁇ -hemolysin-NN is shown in SEQ ID NO: 4.
- the transmembrane protein pore preferably comprises seven monomers each comprising the sequence shown in SEQ ID NO: 4 or a variant thereof.
- Residues 113 and 147 of SEQ ID NO: 4 form part of a constriction of the barrel or channel of ⁇ -HL.
- a pore comprising seven proteins or monomers each comprising the sequence shown in SEQ ID NO: 4 or a variant thereof are preferably used in the method of the invention.
- the seven proteins may be the same (homoheptamer) or different (heteroheptamer).
- a variant of SEQ ID NO: 4 is a protein that has an amino acid sequence which varies from that of SEQ ID NO: 4 and which retains its pore forming ability.
- the ability of a variant to form a pore can be assayed using any method known in the art. For instance, the variant may be inserted into a lipid bilayer along with other appropriate subunits and its ability to oligomerise to form a pore may be determined. Methods are known in the art for inserting subunits into membranes, such as lipid bilayers. Suitable methods are discussed above.
- the variant may include modifications that facilitate covalent attachment to or interaction with the helicase.
- the variant preferably comprises one or more reactive cysteine residues that facilitate attachment to the helicase.
- the variant may include a cysteine at one or more of positions 8, 9, 17, 18, 19, 44, 45, 50, 51, 237, 239 and 287 and/or on the amino or carboxy terminus of SEQ ID NO: 4.
- Preferred variants comprise a substitution of the residue at position 8, 9, 17, 237, 239 and 287 of SEQ ID NO: 4 with cysteine (A8C, T9C, N17C, K237C, S239C or E287C).
- the variant is preferably any one of the variants described in International Application No. PCT/GB09/001,690 (published as WO 2010/004273), PCT/GB09/001,679 (published as WO 2010/004265) or PCT/GB10/000,133 (published as WO 2010/086603).
- the variant may also include modifications that facilitate any interaction with nucleotides.
- the variant may be a naturally occurring variant which is expressed naturally by an organism, for instance by a Staphylococcus bacterium.
- the variant may be expressed in vitro or recombinantly by a bacterium such as Escherichia coli .
- Variants also include non-naturally occurring variants produced by recombinant technology. Over the entire length of the amino acid sequence of SEQ ID NO: 4, a variant will preferably be at least 50% homologous to that sequence based on amino acid identity.
- the variant polypeptide may be at least 55%, at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 90% and more preferably at least 95%, 97% or 99% homologous based on amino acid identity to the amino acid sequence of SEQ ID NO: 4 over the entire sequence.
- Amino acid substitutions may be made to the amino acid sequence of SEQ ID NO: 4 in addition to those discussed above, for example up to 1, 2, 3, 4, 5, 10, 20 or 30 substitutions. Conservative substitutions may be made as discussed above.
- One or more amino acid residues of the amino acid sequence of SEQ ID NO: 4 may additionally be deleted from the polypeptides described above. Up to 1, 2, 3, 4, 5, 10, 20 or 30 residues may be deleted, or more.
- Variants may be fragments of SEQ ID NO: 4. Such fragments retain pore-forming activity. Fragments may be at least 50, 100, 200 or 250 amino acids in length. A fragment preferably comprises the pore-forming domain of SEQ ID NO: 4. Fragments typically include residues 119, 121, 135. 113 and 139 of SEQ ID NO: 4.
- One or more amino acids may be alternatively or additionally added to the polypeptides described above.
- An extension may be provided at the amino terminus or carboxy terminus of the amino acid sequence of SEQ ID NO: 4 or a variant or fragment thereof.
- the extension may be quite short, for example from 1 to 10 amino acids in length. Alternatively, the extension may be longer, for example up to 50 or 100 amino acids.
- a carrier protein may be fused to a pore or variant.
- a variant of SEQ ID NO: 4 is a subunit that has an amino acid sequence which varies from that of SEQ ID NO: 4 and which retains its ability to form a pore.
- a variant typically contains the regions of SEQ ID NO: 4 that are responsible for pore formation.
- the pore forming ability of ⁇ -HL, which contains a ⁇ -barrel, is provided by (3-strands in each subunit.
- a variant of SEQ ID NO: 4 typically comprises the regions in SEQ ID NO: 4 that form ⁇ -strands.
- the amino acids of SEQ ID NO: 4 that form ⁇ -strands are discussed above.
- One or more modifications can be made to the regions of SEQ ID NO: 4 that form ⁇ -strands as long as the resulting variant retains its ability to form a pore. Specific modifications that can be made to the ⁇ -strand regions of SEQ ID NO: 4 are discussed above.
- a variant of SEQ ID NO: 4 preferably includes one or more modifications, such as substitutions, additions or deletions, within its ⁇ -helices and/or loop regions. Amino acids that form ⁇ -helices and loops are discussed above.
- the variant may be modified to assist its identification or purification as discussed above.
- Pores derived from ⁇ -HL can be made as discussed above with reference to pores derived from Msp.
- the transmembrane protein pore is chemically modified.
- the pore can be chemically modified in any way and at any site.
- the transmembrane protein pore is preferably chemically modified by attachment of a molecule to one or more cysteines (cysteine linkage), attachment of a molecule to one or more lysines, attachment of a molecule to one or more non-natural amino acids, enzyme modification of an epitope or modification of a terminus. Suitable methods for carrying out such modifications are well-known in the art.
- the transmembrane protein pore may be chemically modified by the attachment of any molecule. For instance, the pore may be chemically modified by attachment of a dye or a fluorophore.
- any number of the monomers in the pore may be chemically modified.
- One or more, such as 2, 3, 4, 5, 6, 7, 8, 9 or 10, of the monomers is preferably chemically modified as discussed above.
- the molecule (with which the pore is chemically modified) may be attached directly to the pore or attached via a linker as disclosed in International Application Nos. PCT/GB09/001,690 (published as WO 2010/004273), PCT/GB09/001,679 (published as WO 2010/004265) or PCT/GB10/000,133 (published as WO 2010/086603).
- Ratchets that may be used are as follows.
- the translocation of the polymer through the nanopore may be performed in a ratcheted manner. In this case successive k-mers of the polymer are registered with the nanopore. In this manner each measurement is dependent on a particular k-mer. If the registration is held for sufficient time, then a group of plural measurements will be dependent on a particular k-mer. Depending on the nature of the translocation, the period of registration can be unpredictable and may vary in length. Depending on the period of registration, relative to the measurement sampling rate, it might be that there are not plural measurements, or even a signal measurement, that are dependent on every k-mer in the sequence.
- the translocation of the polymer may be controlled by a molecular ratchet that controls the movement of the polymer through the pore.
- the molecular ratchet may be a polymer binding protein.
- the polynucleotide binding protein is preferably a polynucleotide handling enzyme.
- a polynucleotide handling enzyme is a polypeptide that is capable of interacting with and modifying at least one property of a polynucleotide. The enzyme may modify the polynucleotide by cleaving it to form individual nucleotides or shorter chains of nucleotides, such as di- or trinucleotides.
- the enzyme may modify the polynucleotide by orienting it or moving it to a specific position.
- the polynucleotide handling enzyme does not need to display enzymatic activity as long as it is capable of binding the target polynucleotide and controlling its movement through the pore.
- the enzyme may be modified to remove its enzymatic activity or may be used under conditions which prevent it from acting as an enzyme. Such conditions are discussed in more detail below.
- the polynucleotide handling enzyme may be derived from a nucleolytic enzyme.
- the polynucleotide handling enzyme used in the construct of the enzyme is more preferably derived from a member of any of the Enzyme Classification (EC) groups 3.1.11, 3.1.13, 3.1.14, 3.1.15, 3.1.16, 3.1.21, 3.1.22, 3.1.25, 3.1.26, 3.1.27, 3.1.30 and 3.1.31.
- the enzyme may be any of those disclosed in International Application No. PCT/GB10/000,133 (published as WO 2010/086603).
- Preferred enzymes are polymerases, exonucleases, helicases and topoisomerases, such as gyrases.
- Suitable enzymes include, but are not limited to, exonuclease I from E. coli (SEQ ID NO: 8), exonuclease III enzyme from E. coli (SEQ ID NO: 10), RecJ from T. thermophilus (SEQ ID NO: 12) and bacteriophage lambda exonuclease (SEQ ID NO: 14) and variants thereof.
- Three subunits comprising the sequence shown in SEQ ID NO: 14 or a variant thereof interact to form a trimer exonuclease.
- the enzyme is preferably derived from a Phi29 DNA polymerase.
- An enzyme derived from Phi29 polymerase comprises the sequence shown in SEQ ID NO: 6 or a variant thereof.
- a variant of SEQ ID NOs: 6, 8, 10, 12 or 14 is an enzyme that has an amino acid sequence which varies from that of SEQ ID NO: 6, 8, 10, 12 or 14 and which retains polynucleotide binding ability.
- the variant may include modifications that facilitate binding of the polynucleotide and/or facilitate its activity at high salt concentrations and/or room temperature.
- a variant will preferably be at least 50% homologous to that sequence based on amino acid identity. More preferably, the variant polypeptide may be at least 55%, at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 90% and more preferably at least 95%, 97% or 99% homologous based on amino acid identity to the amino acid sequence of SEQ ID NO: 6, 8, 10, 12 or 14 over the entire sequence.
- the variant may differ from the wild-type sequence in any of the ways discussed above with reference to SEQ ID NO: 2.
- the enzyme may be covalently attached to the pore as discussed above.
- the two strategies for single strand DNA sequencing are the translocation of the DNA through the nanopore, both cis to trans and trans to cis, either with or against an applied potential.
- the most advantageous mechanism for strand sequencing is the controlled translocation of single strand DNA through the nanopore under an applied potential.
- Exonucleases that act progressively or processively on double stranded DNA can be used on the cis side of the pore to feed the remaining single strand through under an applied potential or the trans side under a reverse potential.
- a helicase that unwinds the double stranded DNA can also be used in a similar manner.
- sequencing applications that require strand translocation against an applied potential, but the DNA must be first “caught” by the enzyme under a reverse or no potential.
- the single strand DNA exonucleases or single strand DNA dependent polymerases can act as molecular motors to pull the recently translocated single strand back through the pore in a controlled stepwise manner, trans to cis, against the applied potential.
- the single strand DNA dependent polymerases can act as molecular brake slowing down the movement of a polynucleotide through the pore.
- strand sequencing is carried out using a pore derived from Msp and a Phi29 DNA polymerase.
- the method comprises (a) adding the polynucleotide to the solution; (b) allowing the target polynucleotide to interact with a detector in the membrane, which detector comprises a pore derived from Msp and a Phi29 DNA polymerase, such that the polymerase controls the movement of the target polynucleotide through the pore and a proportion of the nucleotides in the target polynucleotide interacts with the pore; and (c) measuring the current passing through the pore during each interaction and thereby determining the sequence of the target polynucleotide, wherein steps (b) and (c) are carried out with a voltage applied across the pore.
- the target polynucleotide When the target polynucleotide is contacted with a Phi29 DNA polymerase and a pore derived from Msp, the target polynucleotide firstly forms a complex with the Phi29 DNA polymerase. When the voltage is applied across the pore, the target polynucleotide/Phi29 DNA polymerase complex forms a complex with the pore and controls the movement of the target polynucleotide through the pore.
- Wild-type Phi29 DNA polymerase has polymerase and exonuclease activity. It may also unzip double stranded polynucleotides under the correct conditions. Hence, the enzyme may work in three modes. This is discussed in more detail below.
- the Phi29 DNA polymerase may comprise the sequence shown in SEQ ID NO: 6 or a variant thereof.
- a variant of SEQ ID NO: 6 is an enzyme that has an amino acid sequence which varies from that of SEQ ID NO: 6 and which retains polynucleotide binding activity.
- the variant must work in at least one of the three modes discussed below. Preferably, the variant works in all three modes.
- the variant may include modifications that facilitate handling of the polynucleotide and/or facilitate its activity at high salt concentrations and/or room temperature.
- a variant will preferably be at least 40% homologous to that sequence based on amino acid identity. More preferably, the variant polypeptide may be at least 50%, at least 55%, at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 90% and more preferably at least 95%, 97% or 99% homologous based on amino acid identity to the amino acid sequence of SEQ ID NO: 6 over the entire sequence. There may be at least 80%, for example at least 85%, 90% or 95%, amino acid identity over a stretch of 200 or more, for example 230, 250, 270 or 280 or more, contiguous amino acids (“hard homology”). Homology is determined as described above. The variant may differ from the wild-type sequence in any of the ways discussed above with reference to SEQ ID NO: 2.
- the salt concentration is typically from 0.15M to 0.6M.
- the salt is preferably KCl.
- the method may be carried out in one of three preferred ways based on the three modes of the Phi29 DNA polymerase. Each way includes a method of proof-reading the sequence.
- the method is preferably carried out using the Phi29 DNA polymerase as a polymerase.
- steps (b) and (c) are carried out in the presence of free nucleotides and an enzyme cofactor such that the polymerase moves the target polynucleotide through the pore against the field resulting from the applied voltage.
- the target polynucleotide moves in the 5′ to 3′ direction.
- the free nucleotides may be one or more of any of the individual nucleotides discussed above.
- the enzyme cofactor is a factor that allows the Phi29 DNA polymerase to function either as a polymerase or an exonuclease.
- the enzyme cofactor is preferably a divalent metal cation.
- the divalent metal cation is preferably Mg 2+ , Mn 2+ , Ca 2+ or Co 2+ .
- the enzyme cofactor is most preferably Mg 2+ .
- the method preferably further comprises (d) removing the free nucleotides such that the polymerase moves the target polynucleotide through the pore with the field resulting from the applied voltage (i.e.
- steps (d) and (e) are also carried out with a voltage applied across the pore.
- the method is preferably carried out using the Phi29 DNA polymerase as an exonuclease.
- steps (b) and (c) are carried out in the absence of free nucleotides and the presence of an enzyme cofactor such that the polymerase moves the target polynucleotide through the pore with the field resulting from the applied voltage.
- the target polynucleotide moves in the 3′ to 5′ direction.
- the method preferably further comprises (d) adding free nucleotides such that the polymerase moves the target polynucleotide through the pore against the field resulting from the applied voltage (i.e.
- steps (d) and (e) are also carried out with a voltage applied across the pore.
- the method is preferably carried out using the Phi29 DNA polymerase in unzipping mode.
- steps (b) and (c) are carried out in the absence of free nucleotides and the absence of an enzyme cofactor such that the polymerase controls the movement of the target polynucleotide through the pore with the field resulting from the applied voltage (as it is unzipped).
- the polymerase acts like a brake preventing the target polynucleotide from moving through the pore too quickly under the influence of the applied voltage.
- the method preferably further comprises (d) lowering the voltage applied across the pore such that the target polynucleotide moves through the pore in the opposite direction to that in steps (b) and (c) (i.e.
- steps (d) and (e) are also carried out with a voltage applied across the pore.
- a helicase is used as a ratchet for the polynucleotide (for example as disclosed in U.S. 61/549,998 (N115020), U.S. 61/581,332 (N115505), U.S. 61/581,340 that are incorporated herein by reference). It has been shown that helicases have a surprisingly high salt tolerance. Helicases can move the target polynucleotide in two directions, namely with or against the field resulting from the applied voltage. Hence, the method may be carried out in one of two preferred modes. Different signals are obtained depending on the direction the target polynucleotide moves through the pore, i.e. in the direction of or against the field.
- Helicases typically move the target polynucleotide through the pore one nucleotide at a time. Helicases can therefore function like a single-base ratchet. This is of course advantageous when sequencing a target polynucleotide because substantially all, if not all, of the nucleotides in the target polynucleotide may be identified using the pore.
- Helicases are capable of controlling the movement of single stranded polynucleotides and double stranded polynucleotides.
- Helicases appear very resistant to the field resulting from applied voltages. Very little movement of the polynucleotide under an “unzipping” condition was observed. This is important because it means that there are no complications from unwanted “backwards” movements when moving polynucleotides against the field resulting from an applied voltage.
- the method comprises: (a) contacting the target polynucleotide with a transmembrane pore and a helicase such that the helicase controls the movement of the target polynucleotide through the pore and nucleotides in the target polynucleotide interact with the pore; and (b) measuring the current passing through the pore during one or more interactions to measure one or more characteristics of the target polynucleotide and thereby characterising the target polynucleotide.
- helicases may work in two modes with respect to the nanopore.
- the two modes are as follows.
- the method is preferably carried out using the helicase such that it moves the target sequence through the pore with the field resulting from the applied voltage.
- the 3′ end of the DNA is first captured in the nanopore, and the enzyme moves the DNA into the nanopore such that the target sequence is passed through the nanopore with the field until it finally translocates through to the trans side of the bilayer.
- the method is preferably carried out such that the enzyme moves the target sequence through the pore against the field resulting from the applied voltage.
- the 5′ end of the DNA is first captured in the nanopore, and the enzyme moves the DNA through the nanopore such that the target sequence is pulled out of the nanopore against the applied field until finally ejected back to the cis side of the bilayer.
- the two modes are as follows.
- the method is preferably carried out using the helicase such that it moves the target sequence through the pore with the field resulting from the applied voltage.
- the 5′ end of the DNA is first captured in the nanopore, and the enzyme moves the DNA into the nanopore such that the target sequence is passed through the nanopore with the field until it finally translocates through to the trans side of the bilayer.
- the method is preferably carried out such that the enzyme moves the target sequence through the pore against the field resulting from the applied voltage.
- the 3′ end of the DNA is first captured in the nanopore, and the enzyme moves the DNA through the nanopore such that the target sequence is pulled out of the nanopore against the applied field until finally ejected back to the cis side of the bilayer.
- Measurement systems that may be used are as follows.
- the methods may be carried out using any apparatus that is suitable for investigating a membrane/pore system in which a pore is inserted into a membrane.
- the method may be carried out using any apparatus that is suitable for transmembrane pore sensing.
- the apparatus comprises a chamber comprising an aqueous solution and a barrier that separates the chamber into two sections.
- the barrier has an aperture in which the membrane containing the pore is formed.
- the methods may be carried out using the apparatus described in International Application No. PCT/GB08/000,562 (WO 2008/102120).
- the methods may involve measuring the current passing through the pore during one or more interactions with the nucleotide(s). Therefore the apparatus may also comprise an electrical circuit capable of applying a potential and measuring an electrical signal across the membrane and pore.
- the methods may be carried out using a patch clamp or a voltage clamp.
- the methods preferably involve the use of a voltage clamp.
- the methods of the invention may involve the measuring of a current passing through the pore during one or more interactions with the nucleotide. Suitable conditions for measuring ionic currents through transmembrane protein pores are known in the art and disclosed in the Example. The method is typically carried out with a voltage applied across the membrane and pore. The voltage used is typically from +2 V to ⁇ 2 V, typically ⁇ 400 mV to +400 mV.
- the voltage used is preferably in a range having a lower limit selected from ⁇ 400 mV, ⁇ 300 mV, ⁇ 200 mV, ⁇ 150 mV, ⁇ 100 mV, ⁇ 50 mV, ⁇ 20 mV and 0 mV and an upper limit independently selected from +10 mV, +20 mV, +50 mV, +100 mV, +150 mV, +200 mV, +300 mV and +400 mV.
- the voltage used is more preferably in the range 100 mV to 240 mV and most preferably in the range of 120 mV to 220 mV. It is possible to increase discrimination between different nucleotides by a pore by using an increased applied potential.
- the methods are typically carried out in the presence of any charge carriers, such as metal salts, for example alkali metal salt, halide salts, for example chloride salts, such as alkali metal chloride salt.
- Charge carriers may include ionic liquids or organic salts, for example tetramethyl ammonium chloride, trimethylphenyl ammonium chloride, phenyltrimethyl ammonium chloride, or 1-ethyl-3-methyl imidazolium chloride.
- the salt is present in the aqueous solution in the chamber. Potassium chloride (KCl), sodium chloride (NaCl) or caesium chloride (CsCl) is typically used. NaCl is preferred.
- the salt concentration may be at saturation.
- the salt concentration may be 3M or lower and is typically from 0.1 to 2.5 M, from 0.3 to 1.9 M, from 0.5 to 1.8 M, from 0.7 to 1.7 M, from 0.9 to 1.6 M or from 1 M to 1.4 M.
- the salt concentration is preferably from 150 mM to 1 M.
- the method is preferably carried out using a salt concentration of at least 0.3 M, such as at least 0.4 M, at least 0.5 M, at least 0.6 M, at least 0.8 M, at least 1.0 M, at least 1.5 M, at least 2.0 M, at least 2.5 M or at least 3.0 M.
- High salt concentrations provide a high signal to noise ratio and allow for currents indicative of the presence of a polymer to be identified against the background of normal current fluctuations.
- the methods are typically carried out in the presence of a buffer.
- the buffer is present in the aqueous solution in the chamber. Any buffer may be used in the method of the invention.
- the buffer is HEPES.
- Another suitable buffer is Tris-HCl buffer.
- the methods are typically carried out at a pH of from 4.0 to 12.0, from 4.5 to 10.0, from 5.0 to 9.0, from 5.5 to 8.8, from 6.0 to 8.7 or from 7.0 to 8.8 or 7.5 to 8.5.
- the pH used is preferably about 7.5.
- the methods may be carried out at from 0° C. to 100° C., from 15° C. to 95° C., from 16° C. to 90° C., from 17° C. to 85° C., from 18° C. to 80° C., 19° C. to 70° C., or from 20° C. to 60° C.
- the methods are typically carried out at room temperature.
- the methods are optionally carried out at a temperature that supports enzyme function, such as about 37° C.
- the method is typically carried out in the presence of free nucleotides or free nucleotide analogues and an enzyme cofactor that facilitate the action of the molecular ratchet or enzyme.
- the free nucleotides may be one or more of any of the individual nucleotides discussed above.
- the free nucleotides include, but are not limited to, adenosine monophosphate (AMP), adenosine diphosphate (ADP), adenosine triphosphate (ATP), guanosine monophosphate (GMP), guanosine diphosphate (GDP), guanosine triphosphate (GTP), thymidine monophosphate (TMP), thymidine diphosphate (TDP), thymidine triphosphate (TTP), uridine monophosphate (UMP), uridine diphosphate (UDP), uridine triphosphate (UTP), cytidine monophosphate (CMP), cytidine diphosphate (CDP), cytidine triphosphate (CTP), cyclic adenosine monophosphate (cAMP), cyclic guanosine monophosphate (cGMP), deoxyadenosine monophosphate (dAMP), deoxyadenosine diphosphate (dADP), deoxyaden
- the free nucleotides are preferably selected from AMP, TMP, GMP, CMP, UMP, dAMP, dTMP, dGMP or dCMP.
- the free nucleotides are preferably adenosine triphosphate (ATP).
- the enzyme cofactor is a factor that allows the enzyme to function.
- the enzyme cofactor is preferably a divalent metal cation.
- the divalent metal cation is preferably Mg 2+ , Mn 2+ , Ca 2+ or Co 2+ .
- the enzyme cofactor is most preferably Mg 2+ .
- the target polymer may be contacted with the molecular ratchet and the pore in any order. In is preferred that, when the target polymer is contacted with the molecular ratchet and the pore, the target polymer firstly forms a complex with the molecular ratchet. When the voltage is applied across the pore, the target polymer/molecular ratchet complex then forms a complex with the pore and controls the movement of the polymer through the pore.
- the nature of the measurements may be as follows.
- a property that depends on the polymer units translocating through the pore may be measured.
- the property may be associated with an interaction between the polymer and the pore. Interaction of the polymer may occur at a constricted region of the pore.
- the measurement system measures the property, producing a measurement that is dependent on the polymer units of the polymer.
- a variety of different types of measurements may be made. This includes without limitation: electrical measurements and optical measurements. Possible electrical measurements include: ion current flow measurements, impedance measurements, tunnelling measurements (Ivanov A P et al., Nano Lett. 2011 Jan. 12; 11(1):279-85), and FET measurements (International Application WO 2005/124888). Optical measurements may be combined with electrical measurements (Soni G V et al., Rev Sci Instrum. 2010 January; 81(1):014301). The measurement may be a transmembrane current measurement such as measurement of ionic current flowing through the pore.
- Electrical measurements may be made using standard single channel recording equipment as describe in Stoddart D et al., Proc Natl Acad Sci, 12; 106(19):7702-7, Lieberman K R et al, J Am Chem Soc. 2010; 132(50):17961-72, and International Application WO-2000/28312.
- electrical measurements may be made using a multi-channel system, for example as described in International Application WO-2009/077734 and International Application WO-2011/067559.
- measurements of more than one property it is possible to use measurements of more than one property.
- one possibility is to use measurements of ion current flow together with measurements of at least one additional property besides ion current flow, for example comprising FET measurements, optical measurements, or both.
- the measurement system may comprise a plurality of pores.
- the apparatus preferably further comprise a plurality of a polymer ratchets.
- the apparatus preferably further comprises instructions for carrying out the method of the invention.
- the apparatus may be any conventional apparatus for polymer analysis, such as an array or a chip. Any of the embodiments discussed above with reference to the methods of the invention are equally applicable to the apparatus of the invention.
- the apparatus is preferably set up to carry out the method of the invention.
- the apparatus may comprises: a sensor device that is capable of supporting the membrane and plurality of pores and being operable to perform polymer characterising using the pores; at least one reservoir for holding material for performing the characterising; a fluidics system configured to controllably supply material from the at least one reservoir to the sensor device; and a plurality of containers for receiving respective samples, the fluidics system being configured to supply the samples selectively from the containers to the sensor device.
- the apparatus may be any of those described in International Application No. No. PCT/GB08/004,127 (published as WO 2009/077734), PCT/GB10/000,789 (published as WO 2010/122293), International Application No. PCT/GB10/002,206 (not yet published) or International Application No. PCT/US99/25679 (published as WO 00/28312), all of which are incorporated herein by reference.
- the apparatus may be a diagnostic device.
- the diagnostic device may be a benchtop or handheld device.
- the device may be operated in conjunction with a cartridge, the cartridge comprising the nanopore assay components and for receiving the fluid sample.
- the cartridge may be housed in the device or otherwise operably connectable with the device.
- the cartridge may be subsequently removed or disconnected from the device in order to clean the cartridge for re-use, or for disposal. Thereafter an unused or cleaned cartridge may be used with the device.
- the cartridge may be an integral part of the device wherein the device is disposable after use.
- the cartridge will typically have a sample application region for receiving a fluid sample.
- the sample application region may be a microfluidic channel or a porous sample pad for example to directly receive a urine sample.
- the size of sample would typically range from 0.25 uL to 10 mL.
- the sample application region may serve to directly receive a sample from a patient, for example a sample of blood obtained with a fingerstick.
- the cartridge may comprise a red blood cell filter for filtering red blood cells.
- the cartridge may comprise dried reagents such as a salt, an anticoagulant, or a buffer.
- the device will typically comprise data input and output ports and a memory for sending or receiving and storing data, such as information in relation to feature vectors, patient ID, and measurement results.
- the device may have wireless connectivity for communicating with a remote server or medical professional.
- the device and cartridge are not restricted to measurement of a particular analyte and may capable of measuring any particular analyte and feature vectors relating to a particular analyte of interest may be uploaded and stored in the memory.
- the measurements are of a property that is associated with an interaction between the polymer and the measurement system.
- measurements that are dependent on larger groups of polymer units i.e. with a “broad” resolution. Such broad resolution may be particularly useful for examining homopolymer regions.
- the measurements are dependent on a k-mer, it is desirable that the measurements are resolvable (i.e. separated) for as many as possible of the possible k-mers. Typically this can be achieved if the measurements produced by different k-mers are well spread over the measurement range and/or have a narrow distribution. This may be achieved to varying extents by different measurement systems. However, it is a particular advantage of the present invention, that it is not essential for the measurements produced by different k-mers to be resolvable.
- FIG. 1 schematically illustrates an example of a measurement system 8 comprising a nanopore that is a biological pore 1 inserted in a biological membrane 2 such as a lipid bilayer.
- a polymer 3 comprising a series of polymer units 4 is translocated through the biological pore 1 as shown by the arrows.
- the polymer 3 may be a polynucleotide in which the polymer units 4 are nucleotides.
- the polymer 3 interacts with an active part 5 of the biological pore 1 causing an electrical property such as the trans-membrane current to vary in dependence on a k-mer inside the biological pore 1 .
- the active part 5 is illustrated as interacting with a k-mer of three polymer units 4 , but this is not limitative. Electrodes 6 arranged on each side of the biological membrane 2 are connected to a measurement circuit 7 that measures the electrical property. Thus the measurements are dependent on the k-mer inside the biological pore 1 .
- a typical type of signal output by a measurement system and which is an input signal to be analysed in accordance with the present invention is a “noisy step wave”, although without limitation to this signal type.
- An example of an input signal having this form is shown in FIG. 2 for the case of an ion current measurement obtained using a measurement system comprising a nanopore.
- This type of input signal comprises an input series of measurements in which successive groups of plural measurements are dependent on the same k-mer.
- the plural measurements in each group are constant, subject to some variance discussed below, and therefore form a “level” in the signal, corresponding to a state of the measurement system.
- the signal moves between a set of levels, which may be a large set. Given the sampling rate of the instrumentation and the noise on the signal, the transitions between levels can be considered instantaneous, thus the signal can be approximated by an idealised step trace.
- the measurements corresponding to each state are constant over the time scale of the event, but for most measurement systems will be subject to variance over a short time scale. Variance can result from measurement noise, for example arising from the electrical circuits and signal processing, notably from the amplifier in the particular case of electrophysiology. Such measurement noise is inevitable due the small magnitude of the properties being measured. Variance can also result from inherent variation or spread in the underlying physical or biological system of the measurement system. Most measurement systems will experience such inherent variation to greater or lesser extents, even in the idealised case that measurement noise is avoided. For any given measurement system, both sources of variation may contribute or one of these noise sources may be dominant.
- each group of measurements may be referred to as a “state”.
- the event consisting of translocation of the polymer through the nanopore may occur in a ratcheted manner.
- the ion current flowing through the nanopore at a given voltage across the nanopore is constant, subject to the variance discussed above.
- each group of measurements is associated with a step of the ratcheted movement.
- Each step corresponds to a state in which the polymer is in a respective position relative to the nanopore.
- the signal may take an arbitrary form.
- the measurements corresponding to k-mers may also be described in terms of a set of emissions and transitions.
- a measurement that is dependent on a particular k-mer may comprise of a series of measurements occurring in a fashion amenable to description by these methods.
- the extent to which a given measurement system provides measurements that are dependent on k-mers and the size of the k-mers may be examined experimentally.
- known polymers may be synthesized and held at predetermined locations relative to the measurement system to investigate from the resultant measurements how the measurements depend on the identity of k-mers that interact with the measurement system.
- One possible approach is to use a set of polymers having identical sequences except for a k-mer at a predetermined position that varies for each polymer of the set.
- the size and identity of the k-mers can be varied to investigate the effect on the measurements.
- FIG. 3 is a frequency distribution of current measurements of two polynucleotides in a measurement system comprising a nanopore.
- every base in the region of the nanopore is a T (labelled polyT)
- N11-TATGAT-N8 11 bases to the left and 8 to the right of a specific fixed 6-mer (having the sequence TATGAT) are allowed to vary.
- This approach can be generalised for any measurement system enabling the location and a minimal k-mer description to be determined.
- Similar methodology may be used to identify location and width of well-approximating k-mers in a general measurement system.
- this is achieved by changing the position of the 6-mer relative to the pore (e.g. by varying the number of Ns before and after) to detect location of the best approximating k-mer and increasing and decreasing the number of fixed bases from 6.
- the value of k can be minimal subject to the spread of values being sufficiently narrow.
- the location of the k-mer can be chosen to minimise peak width.
- FIGS. 4 and 5 are plots of the best model fit against the current measurements. If the data was well described by this model, then the points should closely follow the diagonal line within a typical experimental error (for example 2 pA). This is not the case showing that the data is not well described by this linear model for either set of coefficients.
- the method is illustrated in FIG. 6 and may be computer-implemented in an analysis device 10 illustrated schematically in FIG. 6 .
- the analysis device 10 may be implemented by a computer program executed in a computer apparatus or may be implemented by a dedicated hardware device, or any combination thereof. In either case, the data used by the method is stored in a memory in the analysis device 10 .
- the computer apparatus where used, may be any type of computer system but is typically of conventional construction.
- the computer program may be written in any suitable programming language.
- the computer program may be stored on a computer-readable storage medium (i.e.
- a non-transitory medium which may be of any type, for example: a recording medium which is insertable into a drive of the computing system and which may store information magnetically, optically or opto-magnetically; a fixed recording medium of the computer system such as a hard drive; or a computer memory.
- an input signal 11 that has sufficient time resolution that it comprises a series of measurements (or more generally any number of series, as described further below) of the type described above in which the measurements are time-ordered and comprise successive groups of plural measurements that are dependent on the same k-mer without a priori knowledge of number of measurements in any group.
- FIG. 2 An example of such an input signal 11 is shown in FIG. 2 as previously described.
- a state detection step S 1 the input signal 11 is processed to identify successive groups of measurements.
- the state detection step S 1 may be performed using the method shown in FIG. 7 that looks for short-term increases in the derivative of the input signal 11 as follows.
- step S 1 - 1 the input signal 11 is differentiated to derive its derivative.
- step S 1 - 2 the derivative from step S 1 - 1 is subjected to low-pass filtering to suppress high-frequency noise (which the differentiation tends to amplify).
- step S 1 - 3 the filtered derivative from step S 1 - 2 is thresholded to detect transition points between the groups of measurements, and thereby identify the groups of data.
- the measurements in each identified group are to derive values of one or more features that represent characteristics in respect of each group.
- a single value is derived, for example the mean, but plural values of features that represent the same or different characteristics may be used to increase the information content.
- features include: an average (a mean or a median or other average) of the group of measurements; the period of the group of measurements; a variance of the group of measurements; the distribution of the group of measurements, asymmetry information; the confidence of the measurements; or any combination thereof.
- the values of the features output from step S 2 form a feature vector 12 in which the values are time-ordered in the same order as the groups from which they are derived.
- Step S 2 has the result of providing a representation of the input signal 11 in which the amount of information is reduced, but in which the significant characteristics of the signal are maintained.
- steps S 1 an/or S 2 may alternatively be used in place of steps S 1 an/or S 2 to derive the feature vector 12 of values of one or more features that represent characteristics of the input signal 11 , time-ordered in the same order as the input signal 11 .
- a possible simplification of the state detection step is to use a sliding window analysis whereby one compares the means of two adjacent windows of data.
- a threshold can then be either put directly on the difference in mean, or can be set based on the variance of the data points in the two windows (for example, by calculating Student's t-statistic).
- a particular advantage of these methods is that they can be applied without imposing many assumptions on the data.
- FIG. 8 illustrates an experimentally determined input signal 11 reduced by a moving window t-test.
- FIG. 8 shows the input signal 11 as the light line.
- Levels following state detection are shown overlayed as the dark line.
- FIG. 9 shows the values derived for the entire trace, calculating the level of each state from the mean value between transitions.
- step S 3 the feature vector 12 derived in step S 2 is compared with at least one other feature vector 13 to determine the similarity there between.
- that other feature vector 13 may be one or more feature vectors 14 stored in memory 15 of the analysis device 10 , or alternatively may be one or more feature vectors 12 derived using steps S 1 and S 2 from input signals 11 that are series of measurements of other polymers.
- Step S 3 may be implemented in a variety of manners to derive useful information about the polymer under investigation. Some non-limitative examples of step S 3 are as follows.
- step S 3 the feature vector 12 derived in step S 2 is compared with the other feature vector that is one or more of plural feature vectors 14 stored in a memory 15 of the analysis device 10 in respect of at least one class, as a library.
- step S 3 produces classification data 16 that classifies the polymer from which the derived feature vector 12 is derived as belonging to one of the classes on the basis of the determined similarity.
- similarity may be determined between the entirety or part of the derived feature vector 12 and the entirety of the feature vector 14 stored in the memory 15 , or between the entirety or part of the derived feature vector 12 and a part of the feature vector 14 stored in the memory 15 .
- the method may be repeated on input signals 11 that are series of measurements of other polymers, for example from the same sample. In that case either or both of the following steps S 4 and S 5 may be performed.
- step S 4 the numbers of polymers in each class may be counted. That provides information on the profile of the population of polymers under investigation.
- step S 5 the derived feature vector 12 is compared again with the feature vector feature vector 14 stored in a memory 15 of the class within which the polymer of the derived feature vector 12 is classified as belonging. In this comparison, similarity is again determined, but this time to identify localized regions where the derived feature vector 12 is dissimilar to that feature vector 14 in respect of the class.
- identification of localized regions where the derived feature vector is dissimilar to what is expected provides an analysis technique that is very powerful in many applications where change in relatively small regions of long sequences of polymers is significant.
- One example of such a technique is to identify mutations in a polymer that is a polynucleotide.
- step S 3 the feature vector 13 used for the comparison may be selected from the feature vectors 14 stored in the memory depending upon the polymer to be measured.
- the feature vectors 14 stored in the memory 15 may comprise two or more feature vectors having overlapping regions. In that case, the similarity may be determined in step S 3 with the non-overlapping regions of the feature vectors 14 are used in the determination of similarity with the derived feature vector 12 .
- step S 3 is performed in respect of plural feature vectors 12 derived by performing steps S 1 and S 2 on plural polymers, for example polymers from the same sample or polymers that are fragments of a common polymer.
- step S 3 comprises the following steps.
- step S 3 - 1 the plural derived feature vectors 12 are compared with each other and the similarity therebetween is determined.
- step S 3 - 2 the plural derived feature vectors 12 are clustered on the basis of their similarity. In particular, clusters of similar feature vectors 12 are identified as a class. Step S 3 - 2 produces classification data 16 that classifies the polymers from which each derived feature vector 12 is derived as belonging to one of the classes.
- the classification data 16 may be processed by steps S 4 and/or S 5 as described above.
- step S 3 is performed in respect of plural feature vectors 12 derived by performing steps S 1 and S 2 on plural polymers that are fragments of a common polymer.
- the plural derived feature vectors 12 are compared with each other and the similarity therebetween is determined in overlapping parts of the feature vectors 12 . This allows information on the common polymer to be built up from the input signals of the fragments.
- a fourth example of step S 3 is similar to step S 5 , but involves comparison of the derived feature vector 12 a feature vector 14 stored in a memory 15 . In this comparison, similarity is determined to identify localized regions where the derived feature vector 12 is dissimilar to that feature vector 14 in the memory.
- This fourth example has similar advantages to step S 5 above, but is applicable where the expected type of the polymer is known in advance and so the comparison can be made with a feature vector 14 in respect of that expected type, without needing to classify the derived feature vector 12 first.
- One approach is to modify existing pairwise dynamic programming sequence alignment algorithms e.g. the Needleman-Wunsch algorithm for global alignment or the Smith-Waterman algorithm for local alignment.
- the modifications may include replacing the substitution matrix with a distance measure operating on the feature vector.
- the distance measure may be a measurement of the absolute difference in current between the data points.
- the distance function could also consider multiple measurements at each position e.g. mean and variance of a current measurement.
- Modification may also be made to the to the gap scoring mechanism as are known in the art, for example constant gap penalties, linear gap penalties or affine gap penalties.
- These algorithms output an alignment score that is a function of the two feature vectors, the distance function and the gap penalties.
- the alignment score can be used to determine similarity.
- an alternative approach is to represent the feature vector in terms of shorter sub-vectors, typically comprising consecutive entries in the feature vector. For example, if the feature vector was (1,2,3,4,5) then we could represent it by length 3 sub-vectors to give the new representation ⁇ (1,2,3),(2,3,4),(3,4,5) ⁇ . For our application the sub-vectors are frequently considerably longer (>10) so maintaining much of the time-ordering information.
- the feature sub-vectors are suitably discretized (for example by rounding each number to the nearest 0.1) then exact or partial matches of sub-vectors may be used, and similarity calculated in terms of what proportion of sub-vectors match or partially match. Discretisation also enables integer arithmetic to be used for comparison. Alternatively hash functions may be applied to sub-vectors to give fixed length “fingerprints” (see for instance Karp, R., Rabin, M. (1987) “Efficient randomized pattern matching algorithms”/IBM J. Res. Development 31:249-260.) denoting presence or absence of sub-vectors which can be rapidly compared.
- alignment-based and sub-vector based measures of pairwise similarity treat the pair of feature vectors that are being compared in the same way. The result is that given a pair of feature vectors A and B, the similarity of A to B is equal to the similarity of B to A.
- one of the feature vectors to be compared is a library feature vector
- an alignment can be performed using MEW methods with models constructed in a similar manner to the “forced path” training models described previously (U.S. 61/538,721, GB 1117574.2).
- Algorithms other than Viterbi that are known in the art may also be applied, for example the Forwards-Backwards algorithm.
- this is the total likelihood of the path. The total likelihood is not guaranteed to be equal if we swapped the roles of the two feature vectors, however for classification problems in particular, this is not generally an issue.
- Clustering is performed on input signals 11 from a measured population of polymers, and involves determining the number and/or types of polymer present according to some similarity criteria.
- Hierarchical agglomerative methods are also used for sequence alignment in packages such as CLUSTAL (Higgins, D. G. and Sharp, P. M. (1988). CLUSTAL: a package for performing multiple sequence alignment on a microcomputer. Gene, 73, 237-244.)
- all feature vectors are pairwise aligned with each other such that we have a measure of similarity (or in some cases distance) between each pair of feature vectors.
- These similarity values can be written down as a similarity matrix with the (m,n)th entry containing the similarity of the m′ th to the n′ th feature vector.
- a clustering technique is then used (typically hierarchical agglomerative clustering) based on that similarity matrix.
- Two extremes of agglomerative clustering are single-link (score a pair of clusters during the agglomerative step on the basis of the most similar feature vector pair) and complete-link (score a pair of clusters based on the most dissimilar feature vector pair) clustering.
- the best combination of algorithm to determine similarity and clustering technique is dependent on the nature of the clusters expected for a given application.
- clusters are expected to be made up of feature vectors with overlapping fragments of pairs of feature vectors showing high similarity
- local alignment scores and single-link agglomerative clustering would be one appropriate choice.
- An example of this is shown in Worked Example 2, where sequences 1 and 2 overlap as do sequences 2 and 3. If in our clustering task we wished to identify these as a single cluster amid some other feature vectors, we would be most likely to be successful using a local alignment score to correctly identify the short overlapping regions.
- Single-link clustering would join the sequences into the same cluster because 1 has an overlap with 2 and 2 has an overlap with 3, however complete-link agglomerative clustering would be a poor choice since sequences 1 and 3 have no actual overlap in sequence space and hence are likely to have low similarity in terms of feature vectors.
- clusters are expected to be near-identical across the entire feature vector (for example, where feature vectors have already been identified to begin and end in approximately the same place relative to a known reference, and we are looking to discover classes that vary subtly from that reference) global alignment scores and complete-link agglomerative clustering would be more appropriate.
- the landmark vector can be updated based on many or all possible alignments.
- the landmark vector produced as a result of this process with the feature vectors aligned to it produces a “consensus” of the feature vectors.
- step 1 all pairs of feature vectors may be aligned and the aligned pair with the most states picked, subject to some minimum level of similarity, taking the mean at each position where the states align to generate the initial feature vector.
- Alternatives are possible, for instance just picking the longest feature vector.
- step 2 The pairwise alignment algorithms used in step 2 are described above.
- step 4 p and r can be varied according to the particular situation, mean may be replaced by some other measure of location, and r may be replaced by some other measure of spread.
- This consensus building process provides a multiple alignment algorithm in terms of feature vectors.
- the landmark-aligned states give a fixed length vector representing each feature vector.
- the task for classification is to assign a “query” feature vector to one of m classes for integer m>1.
- Method of solution is dependent on whether the target feature vectors are heterogeneous (mutually dissimilar at a global level) or are homogeneous (all globally similar to each other, with some relatively subtle differences, typically localized, differences), although clearly there are cases that lie between these extremes where mixtures of methods are appropriate.
- the simplest method for class determination is to calculate similarity between the query feature vector and the target feature vectors by one of the methods described above and to assign the query feature vector to the class with the target feature vector of maximum similarity.
- a summary target feature vector may be derived for each class, containing for instance the mean value across target feature vectors in that class, and proceed as before.
- each target feature vector may be treated independently.
- the query feature vector is assigned to the class of the closest target feature vector.
- a re-weighting of statistics to account for the different number of target feature vectors per class is frequently desirable.
- the vector of distances or dissimilarities to target feature vectors can be used as the input to multivariate learning techniques such as multi-class linear discriminant analysis to produce an improved classifier.
- a fixed length vector may be produced from sub-vectors using standard hashing algorithms as described earlier and this used as the input to learning algorithms. More about learning algorithms in the homogeneous case follows.
- the same or similar methods may be applied as in the heterogeneous case, however random variation across the feature vector may well mask the systematic local variations that are of primary interest and provide the key information to correctly discriminate between classes.
- feature vectors may initially be aligned to a common reference feature vector (for instance the landmarks from consensus alignment of the target feature vectors), similarly to the “Consensus Building” case above, and the states aligning to landmarks fed as a fixed length input vector to learning algorithms.
- a common reference feature vector for instance the landmarks from consensus alignment of the target feature vectors
- a decision tree classifier for example, but not limited, to C4.5. Quinlan, J. R. (1993) C4.5: Programs for Machine Learning. Morgan Kaufmann Publishers
- black box methods such as neural networks, random forests and support vector machines may be used to make predictions of class membership, while not necessarily generating interpretable rules.
- Bayesian networks may be implemented, where expert knowledge may also be incorporated.
- sub-vectors as inputs to learning algorithms.
- a fixed length vector may be produced from sub-vectors using standard hashing algorithms as discussed earlier and this used as the input to learning algorithms.
- the sub-vectors themselves may be used directly—for instance with an algorithm searching for sub-vectors that only have near neighbours within-class.
- step S 5 the fourth example of step S 3 .
- an alignment to the target feature vector is performed and then positions that vary between the query feature vector and the target feature vector are identified.
- a reference feature vector is generated (for example the landmarks described in the “Consensus Building” above) from the target feature vectors, and the target feature vectors are aligned to the reference feature vector to gain an idea of location and variability at each position in the reference (for example by calculating the mean and standard deviation of the aligned target feature vectors at that position).
- Localised regions where the query feature vector shows a pattern of values unlikely to be produced in the target class may then be identified, for example by looking at the total likelihood across a number of consecutive reference-aligned states if we assume each distribution to be Gaussian with mean and standard deviation estimated from the target feature vectors.
- the methods may be extended to look at differences between classes of feature vectors as discussed in the homogenous case of classification above.
- classes may be pre-defined, for example they may be DNA samples from patients with and without a particular disease. Alternatively, they may be derived by clustering in the first instance.
- Feature vectors are first “discretised”.
- a transformation is applied to each series of measurments may include any one or combination of the following:
- a first application is in counting molecules against a known library or panel of molecules, which may use a method involving the first example of step S 3 .
- the library comprises the feature vectors 14 stored in the memory 15 .
- Such libraries may be generated for later use, using either supervised or un-supervised learning, based on individual experiments for each molecule or for sets of molecules to learn the feature vectors.
- the fingerprints of these molecules may be known in advance, either from measurements or generated from a model. Given a measurement of a molecule, this can be compared against the known library and the similarity of the molecule to the library members measured. This allows identification of each molecule measured (this identification could be “other”) and quantification of the relative numbers of each type of molecule measured.
- Expression profiles comparing abundance of mRNA transcripts by matching feature vectors. This can be used to measure changes in expression levels. Such gene expression might change during development, disease, treatment for the disease, between one organ and another.
- biomarker miRNAs these are, typically, 20-25-mer RNA oligonucleotides that circulate in blood, and changes in the expression level of groups of these is associated with certain diseases, particularly cancers.
- biomarker miRNAs these are, typically, 20-25-mer RNA oligonucleotides that circulate in blood, and changes in the expression level of groups of these is associated with certain diseases, particularly cancers.
- Foetal copy number variation in circulating blood fragmented foetal DNA circulates in maternal blood. If the foetus has an aneuploidy, e.g. an additional copy of chromosome 21, 18, 11 (the main ones that are not immediately fatal) it would be possible to design capture probes to, for example, exons of the chromosomes of interest, so as to enrich them for pore analysis, and to then compare these to reference feature vectors and to count.
- the main limitation of current methods for this is the inability to distinguish between maternal and foetal chromosomes. There are differences in methylation status between foetal and maternal DNA that are not visible to next gen sequencing that uses PCR, but which would be visible as differences in feature vectors.
- Comparative genomic hybridization changes in copy number of various genomic regions can be altered in tumour cells (and also in foetuses, as described above). For a while, this was identified by comparative genome hybridisation, i.e. where patient/sample is compared to reference by hybridising fragmented genomic DNA to a set of probes on an array. As with foetal testing, feature vector space can be used to profile these copy number changes.
- Viral or bacterial load a measure of the severity of infection. Possibly in conjunction with some form of enrichment, the number of pathogen RNA or DNA copies per ml of blood is measured. It would not have to be done on the whole pathogen genome. Early stage and late stage measurements may be carried out to identify antigenic drift and/or antigenic variation.
- the method may be applications in epidemiology, for example in identification (strain typing) and how a disease is spreading or evolving.
- the method may be used for example to monitor the efficacy of a particular drug therapy or monitor the spread of disease from one area of the body to another or spread of disease between patients.
- Probes There is provided a small population of probes (eg. aptamers to a biomarker panel), some of which attach to a target molecule. Those that didn't bind are separated from those that did, and molecules are counted in the unbound population or in the bound population to quantify the target molecules.
- probes eg. aptamers to a biomarker panel
- the identity of organisms could be determined, e.g. in food or in cultures.
- a second application is in quantification of major population(s) and measurement of “other” present in a sample.
- synthesis of a DNA oligomer can involve Polyacrylamide Gel Electrophoresis (PAGE), High-performance liquid chromatography (HPLC) and mass spectrometry.
- PAGE Polyacrylamide Gel Electrophoresis
- HPLC High-performance liquid chromatography
- mass spectrometry One could measure a sample of the synthesised DNA and determine the signature of the major population present. It is then possible to count the number of molecules in the sample that are different from the major population, possibly indicating errors in synthesis. In particular if the differences occur at particular positions in the feature vectors it may be due to a systematic error that can be rectified by adjusting the synthesis conditions. Any improvements could then be verified by repeating the nanopore measurement.
- a third application is in measurement of modifications/differences at positions and quantification of those modifications/differences within a population of molecules.
- SNP single-nucleotide polymorphisms
- Measurements can be made at known methylation sites.
- the method allows identification of the presence, absence and/or quantification of methylation at those sites.
- the method also allows identification of unknown sites.
- the method allows estimation of “bulk” methylation state of individual molecules, for example whether 100% of the population is 50% modified or 50% of the population is 100% modified, e.g. for use in foetal screening as described above.
- the methylation state of certain genes can be used as a biomarker for cancer
- Another example is identification of splice variants and/or translocation breakpoints. This is similar to the examples described earlier, but one identifies the position where feature vectors stop matching, or where one half of a feature vector maps to one locus and the other half maps elsewhere.
- a fourth application is in identification of the presence and/or absence to a desired confidence of a particular known molecule.
- This method may be used to identify populations of molecules that are related, but not identical to the known molecule to a particular degree of confidence (similar to measurements of homology of DNA or protein sequences), for example in rapidly mutating diseases.
- Another example concerns fusion transcripts, as in splice variants. Detection of specific fusion transcripts is used in cancer diagnosis, e.g. the presence of the Bcl-abl fusion transcript indicates leukaemia.
- Another example concerns diagnosis of NAHR.
- recombination between similar but non-allelic loci results in deletion or duplication of fairly large chunks of genome, with catastrophic consequences for the foetus arising from such a gamete. This will cause a change in copy number of the affected loci (see CGH above), but also results in a fusion of the non-allelic homologs, which would be detectable by looking at PSVs (like SNPs but not the same).
- Another example concerns the case where plural parts of the derived feature vector are compared to plural stored feature vectors.
- the DNA sequence for known protein domains may be used to generate the library feature vectors and the DNA that codes for an unknown protein measured.
- Part of the derived feature vector may be identified with for example a catalytic domain and another part with for example a DNA binding domain.
- the function of the protein may be deduced.
- a fifth application concerns assembly.
- Libraries may be derived experimentally or be generated informatically.
- Examples of the type of library used may include without limitation feature vectors constructed from known DNA sequences using a model, from known protein sequences, from known polymers, feature vectors derived experimentally, feature vectors assembled from overlapping derived feature vectors, feature vectors derived from the consensus of clustered measurements. Libraries may comprise plural related feature vectors, plural unrelated feature vectors, heterogeneous or homogeneous sizes of feature vectors, similar feature vectors with localised differences.
- libraries of feature vectors corresponding to DNA fragments are derived experimentally may use fragments produced systematically, for example by enzymatic fragmentation, or fragments produced randomly, for example by mechanical shearing or through non-selective enzyme action. Randomly fragmented derived feature vectors may preferentially be assembled into larger derived feature vectors for use in a library. Systematically fragmented libraries may preferentially be used as library feature vectors covering similar regions to the fragmentation pattern.
- a library of feature vectors is derived informatically may utilise available databases, for example the NIH Genbank database (Nucleic Acids Research, 2011 January; 39(Database issue):D32-7) which contains publically available DNA sequences.
- a model may be used derived from a training process such as that used previously (U.S. 61/538,721, GB 1117574.2, N114722).
- the library may be reduced to those sequences that are of interest for the particular application, for example the library may be reduced to coding regions of the human genome.
- Example 1 concerns data acquisition in a typical nanopore experiment with the following experimental conditions:
- Buffered solution 1 M NaCl, 100 mM Hepes pH 8.0, 1 mM ATP, 1 mM MgCl 2 , 1 mM DTT, 10 mM Potassium Ferrocyanide (II), 10 mM Potassium Ferricyanide (III), Pt electrodes.
- Enzyme Helicase 100 nM
- Helicase ATPase activity was initiated as required by the addition of divalent metal (1 mM MgCl 2 ) and NTP (1 mM ATP) to the cis compartment. Experiments were carried out at a constant potential of +120 mV.
- the analyte DNA samples used in this study are shown as ANA ID NO. 1-19.
- Example 2 concerns identification and quantification of particular DNA molecules from a panel of DNA molecules. This example describes the process of identification of DNA molecules in a solution from a pre-determined library of feature vectors.
- the library was constructed by taking 18 approximately 400 mer sequences (ANA ID NO 1 to 18), each overlapping the previous sequence by approximately 100 bases from a 5 kilobase genome (PhiX174). For example, ANA ID NO 2 will share 100 bases with ANA ID NO 1 and 100 bases with ANA ID NO 3). These sequences contain a sequence at the beginning and a sequence at the end, common to all strands and not part of the larger genome. The overlapping sequences allow a demonstration of identification of different molecules, even in the presence of large similar regions.
- the library feature vectors are constructed for the mean current by considering a model of the current associated with each 5 mer position (1024 values). The determination of this type of model has been disclosed previously (for example in U.S. 61/538,721, GB 1117574.2, N114722).
- FIG. 12 Feature Vectors for sequences 1, 2 and 3 are shown in FIG. 12 , which illustrates overlapping sections. Common ends (as described above) of each sequence have been removed for this illustration.
- Candidate molecule feature vectors were obtained as follows. Candidate molecules were acquired using the experimental methods describe above and in Example 1. Candidates are reduced to feature vectors consisting of the mean current between identified transitions as described previously.
- An example candidate belonging to one of the sequences (ANA ID NO 1 to 18) was considered. This molecule was compared against the library (ANA ID NO 1 to 18) using an alignment algorithm as described above. The output score from the alignment is used as a measurement of similarity to each of the library members.
- Example 3 concerns measurement of single-nucleotide polymorphisms (SNPs) in a DNA fragment.
- FIG. 19 shows the difference in current between Viterbi aligned library and candidate feature vectors.
- the three SNPs are visible, in the case of 335 and 357 at several positions as several of the measured features are affected by each single change (i.e. a single change to sequences affects several adjacent kmers).
- Example 4 concerns identification of a major population and measurement of a sub-population that is similar but different.
- This example is worked through with simulated data.
- a set of 60 feature vectors (of mean current) is simulated of ANA ID NO 13.
- Ten of the simulations also contain a SNP.
- Gaussian noise with standard deviation of 1 pA is added to each value and 5% of values within each vector are deleted at random. Apart from simulating the data, no more knowledge of the sequence is used.
- FIG. 21 shows the final output of this process with all the data aligned to the consensus. We clearly see the region where the SNP is contained at approx. position 337.
- Example 3 Performing the same analysis as for Example 3 one can see the SNP usually identified in molecules 51-60, as shown in FIG. 22 .
- Example 5 concerns identification of a number of populations, generation of a library and a relative count.
- FIGS. 28 and 29 show the counts against the three clusters for the two cluster and three cluster experiment. We see we have correctly quantified the mixtures in each experiment.
- Example 6 concerns assembly of a larger library feature vector from smaller feature vectors.
- This example uses simulated data from the overlapping sequences S 1 -S 18 as described above. However, to illustrate the assembly process we remove the sequences at the start and end, common to all sequences (as described in Example 2) such that the sequences overlap without any mis-matched regions (as was shown in FIG. 12 ). Since the sequences are guaranteed to be overlapping a relatively simple method can be used. Were this not the case one could use more complex assembly algorithms adapted from those known in the art as described above.
- a tree by neighbour joining on pairwise alignment scores was constructed, similar to Example 5. However since relatively large non-similar regions were expected, there was used a scoring function that does not penalise gaps at the beginning or end of the alignment as strongly as those within the alignment.
- the tree is shown in FIG. 30 . Here it can be seen that all the sequences have similar relation to two other sequences, representing the ⁇ 100 base overlap each sequence shares with the sequences either side.
- consensus landmarks for the aligned sequences are constructed with the output landmarks from a pair of sequences acting as the feature vector where that pair is joined to another sequence.
- the output of the process is a fully assembled feature vector.
- the original data was aligned to the assembled features for illustration. Alignments for three fragments are shown in FIG. 31 , wherein the overlaps can be clearly seen.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Organic Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biophysics (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Biotechnology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Hematology (AREA)
- Urology & Nephrology (AREA)
- Nanotechnology (AREA)
- General Physics & Mathematics (AREA)
- Microbiology (AREA)
- Pathology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- General Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Bioinformatics & Computational Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Biology (AREA)
- Cell Biology (AREA)
- Optics & Photonics (AREA)
- Databases & Information Systems (AREA)
- Human Computer Interaction (AREA)
- Electrochemistry (AREA)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/379,242 US20150057948A1 (en) | 2012-02-16 | 2013-02-18 | Analysis of measurements of a polymer |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261599573P | 2012-02-16 | 2012-02-16 | |
PCT/GB2013/050381 WO2013121224A1 (en) | 2012-02-16 | 2013-02-18 | Analysis of measurements of a polymer |
US14/379,242 US20150057948A1 (en) | 2012-02-16 | 2013-02-18 | Analysis of measurements of a polymer |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/GB2013/050381 A-371-Of-International WO2013121224A1 (en) | 2012-02-16 | 2013-02-18 | Analysis of measurements of a polymer |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/245,306 Continuation US11959906B2 (en) | 2012-02-16 | 2019-01-11 | Analysis of measurements of a polymer |
Publications (1)
Publication Number | Publication Date |
---|---|
US20150057948A1 true US20150057948A1 (en) | 2015-02-26 |
Family
ID=47749886
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/379,242 Abandoned US20150057948A1 (en) | 2012-02-16 | 2013-02-18 | Analysis of measurements of a polymer |
US16/245,306 Active 2034-07-14 US11959906B2 (en) | 2012-02-16 | 2019-01-11 | Analysis of measurements of a polymer |
US18/591,644 Pending US20240264143A1 (en) | 2012-02-16 | 2024-02-29 | Analysis of measurements of a polymer |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/245,306 Active 2034-07-14 US11959906B2 (en) | 2012-02-16 | 2019-01-11 | Analysis of measurements of a polymer |
US18/591,644 Pending US20240264143A1 (en) | 2012-02-16 | 2024-02-29 | Analysis of measurements of a polymer |
Country Status (10)
Country | Link |
---|---|
US (3) | US20150057948A1 (de) |
EP (2) | EP3736339B1 (de) |
JP (1) | JP6226888B2 (de) |
KR (1) | KR102106499B1 (de) |
CN (1) | CN104321441B (de) |
AU (1) | AU2013220179B2 (de) |
BR (1) | BR112014020211A2 (de) |
CA (1) | CA2866587C (de) |
IN (1) | IN2014DN06795A (de) |
WO (1) | WO2013121224A1 (de) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018034745A1 (en) * | 2016-08-18 | 2018-02-22 | The Regents Of The University Of California | Nanopore sequencing base calling |
US9984147B2 (en) | 2008-08-08 | 2018-05-29 | The Research Foundation For The State University Of New York | System and method for probabilistic relational clustering |
US10131943B2 (en) | 2012-12-19 | 2018-11-20 | Oxford Nanopore Technologies Ltd. | Analysis of a polynucleotide via a nanopore system |
WO2019204720A1 (en) * | 2018-04-20 | 2019-10-24 | The Regents Of The University Of California | Nucleic acid sequencing methods and computer-readable media for practicing same |
US20200176082A1 (en) * | 2018-11-28 | 2020-06-04 | Oxford Nanopore Technologies Ltd. | Analysis of nanopore signal using a machine-learning technique |
US10689697B2 (en) | 2014-10-16 | 2020-06-23 | Oxford Nanopore Technologies Ltd. | Analysis of a polymer |
US10783268B2 (en) | 2015-11-10 | 2020-09-22 | Hewlett Packard Enterprise Development Lp | Data allocation based on secure information retrieval |
US10823721B2 (en) | 2016-07-05 | 2020-11-03 | Quantapore, Inc. | Optically based nanopore sequencing |
WO2020236995A1 (en) * | 2019-05-22 | 2020-11-26 | Icahn School Of Medicine At Mount Sinai | Method for de novo detection, identification and fine mapping of multiple forms of nucleic acid modifications |
CN112703256A (zh) * | 2018-09-04 | 2021-04-23 | 牛津纳米孔科技公司 | 用于测定聚合物序列的方法 |
US11080301B2 (en) | 2016-09-28 | 2021-08-03 | Hewlett Packard Enterprise Development Lp | Storage allocation based on secure data comparisons via multiple intermediaries |
US11514289B1 (en) * | 2016-03-09 | 2022-11-29 | Freenome Holdings, Inc. | Generating machine learning models using genetic data |
US11921103B2 (en) | 2011-09-23 | 2024-03-05 | Oxford Nanopore Technologies Plc | Method of operating a measurement system to analyze a polymer |
US11959906B2 (en) | 2012-02-16 | 2024-04-16 | Oxford Nanopore Technologies Plc | Analysis of measurements of a polymer |
Families Citing this family (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0724736D0 (en) | 2007-12-19 | 2008-01-30 | Oxford Nanolabs Ltd | Formation of layers of amphiphilic molecules |
GB201202519D0 (en) | 2012-02-13 | 2012-03-28 | Oxford Nanopore Tech Ltd | Apparatus for supporting an array of layers of amphiphilic molecules and method of forming an array of layers of amphiphilic molecules |
GB201313121D0 (en) | 2013-07-23 | 2013-09-04 | Oxford Nanopore Tech Ltd | Array of volumes of polar medium |
GB201418512D0 (en) | 2014-10-17 | 2014-12-03 | Oxford Nanopore Tech Ltd | Electrical device with detachable components |
US10302972B2 (en) | 2015-01-23 | 2019-05-28 | Pacific Biosciences Of California, Inc. | Waveguide transmission |
JP6800862B2 (ja) | 2015-02-05 | 2020-12-16 | プレジデント・アンド・フェロウズ・オブ・ハーバード・カレッジ | 流体通路を含むナノポアセンサ |
GB201508669D0 (en) | 2015-05-20 | 2015-07-01 | Oxford Nanopore Tech Ltd | Methods and apparatus for forming apertures in a solid state membrane using dielectric breakdown |
EP3332032B1 (de) | 2015-08-06 | 2023-12-27 | Pacific Biosciences Of California, Inc. | Integrierte schaltungen zur selektiven adressierung spärlich angeordneter, elektronischer messvorrichtungen |
EP4019543A1 (de) | 2016-03-02 | 2022-06-29 | Oxford Nanopore Technologies plc | Mutante pore |
GB201611770D0 (en) | 2016-07-06 | 2016-08-17 | Oxford Nanopore Tech | Microfluidic device |
MX2019000037A (es) * | 2016-07-06 | 2019-07-10 | Guardant Health Inc | Metodos para perfilado de fragmentoma de acidos nucleicos libres de celula. |
GB201619930D0 (en) | 2016-11-24 | 2017-01-11 | Oxford Nanopore Tech | Apparatus and methods for controlling insertion of a membrane channel into a membrane |
US20200251184A1 (en) * | 2016-12-16 | 2020-08-06 | Osaka University | Classification analysis method, classification analysis device, and storage medium for classification analysis |
GB201707122D0 (en) | 2017-05-04 | 2017-06-21 | Oxford Nanopore Tech Ltd | Pore |
GB201707138D0 (en) * | 2017-05-04 | 2017-06-21 | Oxford Nanopore Tech Ltd | Machine learning analysis of nanopore measurements |
KR102620326B1 (ko) * | 2017-09-15 | 2024-01-02 | 일루미나, 인코포레이티드 | 서열-검출 시스템 |
GB2569630B (en) | 2017-12-21 | 2022-10-12 | Sharp Life Science Eu Ltd | Droplet Interfaces in Electro-wetting Devices |
SG11201903333SA (en) | 2017-12-29 | 2019-08-27 | Clear Labs Inc | Automated priming and library loading services |
EP3745857A1 (de) | 2018-02-02 | 2020-12-09 | Bayer Aktiengesellschaft | Bekämpfung resistenter schadorganismen |
US20210205814A1 (en) | 2018-05-24 | 2021-07-08 | Oxford Nanopore Technologies Limited | Droplet interfaces in electro-wetting devices |
EP3887046A1 (de) | 2018-11-28 | 2021-10-06 | Oxford Nanopore Technologies Limited | Messsystem und betriebsverfahren |
CN113574381A (zh) | 2019-03-12 | 2021-10-29 | 牛津纳米孔科技公司 | 纳米孔感测装置、组件和操作方法 |
CN110628597B (zh) * | 2019-10-19 | 2023-04-25 | 上海新微技术研发中心有限公司 | 拉曼光谱法生物分子测序方法 |
US11651281B2 (en) * | 2020-05-18 | 2023-05-16 | International Business Machines Corporation | Feature catalog enhancement through automated feature correlation |
CN113030548A (zh) * | 2021-03-09 | 2021-06-25 | 南京邮电大学 | 基于局域阈值的纳米孔检测电流信号特征信息提取方法及装置 |
GB202103605D0 (en) * | 2021-03-16 | 2021-04-28 | Oxford Nanopore Tech Ltd | Alignment of target and reference sequences of polymer units |
US20240240248A1 (en) | 2021-05-19 | 2024-07-18 | Oxford Nanopore Technologies Plc | Methods for complement strand sequencing |
WO2023029044A1 (zh) * | 2021-09-06 | 2023-03-09 | 百图生科(北京)智能技术有限公司 | 用于单细胞测序的方法、装置、设备、介质和程序产品 |
CN113837036B (zh) * | 2021-09-09 | 2024-08-02 | 成都齐碳科技有限公司 | 生物聚合物的表征方法、装置、设备及计算机存储介质 |
CN114908358B8 (zh) * | 2022-04-14 | 2024-06-21 | 孔确(成都)科技有限公司 | 一种两亲性分子层的制备方法及装置 |
GB202215442D0 (en) | 2022-10-19 | 2022-11-30 | Oxford Nanopore Tech Plc | Analysis of a polymer |
WO2024099985A1 (en) | 2022-11-10 | 2024-05-16 | Bayer Aktiengesellschaft | Targeted crop protection product application based on genetic profiles |
GB202216905D0 (en) | 2022-11-11 | 2022-12-28 | Oxford Nanopore Tech Plc | Novel pore monomers and pores |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6128587A (en) * | 1997-01-14 | 2000-10-03 | The Regents Of The University Of California | Method and apparatus using Bayesian subfamily identification for sequence analysis |
Family Cites Families (77)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5795782A (en) * | 1995-03-17 | 1998-08-18 | President & Fellows Of Harvard College | Characterization of individual polymer molecules based on monomer-interface interactions |
JPH11178575A (ja) * | 1997-12-22 | 1999-07-06 | Hitachi Ltd | Dna塩基配列解析装置、方法及び記録媒体 |
US6267872B1 (en) | 1998-11-06 | 2001-07-31 | The Regents Of The University Of California | Miniature support for thin films containing single channels or nanopores and methods for using same |
WO2000039333A1 (en) | 1998-12-23 | 2000-07-06 | Jones Elizabeth Louise | Sequencing method using magnifying tags |
JP2003502166A (ja) | 1999-06-22 | 2003-01-21 | プレジデント・アンド・フェローズ・オブ・ハーバード・カレッジ | 固体状態の次元的特徴の制御 |
AU2002239284A1 (en) | 2000-11-27 | 2002-06-03 | The Regents Of The University Of California | Methods and devices for characterizing duplex nucleic acid molecules |
US20020197618A1 (en) | 2001-01-20 | 2002-12-26 | Sampson Jeffrey R. | Synthesis and amplification of unstructured nucleic acids for rapid sequencing |
JP2002325581A (ja) * | 2001-04-27 | 2002-11-12 | Adgene Co Ltd | 核酸溶解曲線及び核酸解離曲線を用いた未知あるいは既知核酸変異検出法及び表示法 |
JP2003256435A (ja) | 2002-03-06 | 2003-09-12 | Fujitsu Ltd | 配列データ統合処理方法、配列データ統合処理装置及び配列データ統合処理プログラム |
US7220405B2 (en) | 2003-09-08 | 2007-05-22 | E. I. Du Pont De Nemours And Company | Peptide-based conditioners and colorants for hair, skin, and nails |
JP2005176730A (ja) | 2003-12-19 | 2005-07-07 | Hitachi Ltd | cDNA配列をゲノム配列にマッピングする方法 |
US20050136408A1 (en) | 2003-12-19 | 2005-06-23 | May Tom-Moy | Methods and systems for characterizing a polymer |
US7279337B2 (en) | 2004-03-10 | 2007-10-09 | Agilent Technologies, Inc. | Method and apparatus for sequencing polymers through tunneling conductance variation detection |
WO2006028508A2 (en) * | 2004-03-23 | 2006-03-16 | President And Fellows Of Harvard College | Methods and apparatus for characterizing polynucleotides |
KR100614827B1 (ko) | 2004-05-06 | 2006-08-25 | 재단법인서울대학교산학협력재단 | 양방향 은닉 마코프 모델을 이용한 완숙한마이크로알엔에이 위치예측방법 및 이를 구현하기 위한컴퓨터 프로그램을 기록한 저장매체 |
WO2005124888A1 (en) | 2004-06-08 | 2005-12-29 | President And Fellows Of Harvard College | Suspended carbon nanotube field effect transistor |
US20060019259A1 (en) | 2004-07-22 | 2006-01-26 | Joyce Timothy H | Characterization of biopolymers by resonance tunneling and fluorescence quenching |
US20130071837A1 (en) | 2004-10-06 | 2013-03-21 | Stephen N. Winters-Hilt | Method and System for Characterizing or Identifying Molecules and Molecular Mixtures |
US20060086626A1 (en) | 2004-10-22 | 2006-04-27 | Joyce Timothy H | Nanostructure resonant tunneling with a gate voltage source |
GB0505971D0 (en) | 2005-03-23 | 2005-04-27 | Isis Innovation | Delivery of molecules to a lipid bilayer |
WO2007065025A2 (en) | 2005-11-29 | 2007-06-07 | Wisconsin Alumni Research Foundation | Method of dna analysis using micro/nanochannel |
WO2007117832A2 (en) | 2006-03-12 | 2007-10-18 | Applera Corporation | Methods of detecting target nucleic acids |
US8116988B2 (en) | 2006-05-19 | 2012-02-14 | The University Of Chicago | Method for indexing nucleic acid sequences for computer based searching |
US7731826B2 (en) | 2006-08-17 | 2010-06-08 | Electronic Bio Sciences, Llc | Controlled translocation of a polymer in an electrolytic sensing system |
US8003319B2 (en) | 2007-02-02 | 2011-08-23 | International Business Machines Corporation | Systems and methods for controlling position of charged polymer inside nanopore |
NZ579083A (en) | 2007-02-20 | 2012-07-27 | Oxford Nanopore Tech Ltd | Lipid bilayer sensor system |
JP5646987B2 (ja) | 2007-04-04 | 2014-12-24 | ザ リージェンツ オブ ザ ユニバーシティ オブ カリフォルニア | ナノポアを使用するための組成物、デバイス、システム、及び方法 |
WO2009035647A1 (en) | 2007-09-12 | 2009-03-19 | President And Fellows Of Harvard College | High-resolution molecular graphene sensor comprising an aperture in the graphene layer |
GB2453377A (en) * | 2007-10-05 | 2009-04-08 | Isis Innovation | Transmembrane protein pores and molecular adapters therefore. |
GB0724736D0 (en) | 2007-12-19 | 2008-01-30 | Oxford Nanolabs Ltd | Formation of layers of amphiphilic molecules |
US8628940B2 (en) | 2008-09-24 | 2014-01-14 | Pacific Biosciences Of California, Inc. | Intermittent detection during analytical reactions |
CA2730068A1 (en) * | 2008-07-07 | 2010-01-14 | Oxford Nanopore Technologies Limited | Base-detecting pore |
WO2010004265A1 (en) | 2008-07-07 | 2010-01-14 | Oxford Nanopore Technologies Limited | Enzyme-pore constructs |
ES2781858T3 (es) | 2008-09-22 | 2020-09-08 | Univ Washington | Nanoporos MSP y procedimientos relacionados |
WO2010053820A1 (en) | 2008-10-29 | 2010-05-14 | Trustees Of Boston University | Sequence preserved dna conversion |
US8452546B1 (en) | 2008-11-07 | 2013-05-28 | Electronic Biosciences, Inc. | Method for deducing a polymer sequence from a nominal base-by-base measurement |
GB0820927D0 (en) * | 2008-11-14 | 2008-12-24 | Isis Innovation | Method |
EP2391655B1 (de) | 2009-01-30 | 2017-10-11 | Oxford Nanopore Technologies Limited | Hybridisierungslinker |
JP5843614B2 (ja) | 2009-01-30 | 2016-01-13 | オックスフォード ナノポア テクノロジーズ リミテッド | 膜貫通配列決定における核酸構築物のためのアダプター |
GB0905140D0 (en) | 2009-03-25 | 2009-05-06 | Isis Innovation | Method |
WO2010117470A2 (en) | 2009-04-10 | 2010-10-14 | Pacific Biosciences Of California, Inc. | Nanopore sequencing devices and methods |
KR20110138286A (ko) | 2009-04-20 | 2011-12-26 | 옥스포드 나노포어 테크놀로지즈 리미티드 | 지질 이분자층 센서 어레이 |
AU2010307229B2 (en) | 2009-09-18 | 2016-02-25 | President And Fellows Of Harvard College | Bare single-layer graphene membrane having a nanopore enabling high-sensitivity molecular detection and analysis |
CN102741430B (zh) | 2009-12-01 | 2016-07-13 | 牛津楠路珀尔科技有限公司 | 生化分析仪器、用于进行生化分析的第一模块以及相关方法 |
US8748091B2 (en) | 2009-12-18 | 2014-06-10 | The Board Of Trustees Of The University Of Illinois | Characterizing stretched polynucleotides in a synthetic nanopassage |
US8324914B2 (en) | 2010-02-08 | 2012-12-04 | Genia Technologies, Inc. | Systems and methods for characterizing a molecule |
EP4268944A3 (de) | 2010-02-23 | 2024-03-20 | University of Washington | Analytsequenzierung mit nanoporen |
CN103154729B (zh) | 2010-06-08 | 2015-01-07 | 哈佛大学校长及研究员协会 | 具有由石墨烯支持的人工脂质膜的纳米孔装置 |
US20140051068A1 (en) | 2010-09-07 | 2014-02-20 | The Regents Of The University Of California | Control of dna movement in a nanopore at one nucleotide precision by a processive enzyme |
WO2012109483A2 (en) | 2011-02-09 | 2012-08-16 | Life Technologies Corporation | A method of analysis of genetic markers |
US9751915B2 (en) | 2011-02-11 | 2017-09-05 | Oxford Nanopore Technologies Ltd. | Mutant pores |
WO2012135658A2 (en) | 2011-03-30 | 2012-10-04 | Noblegen Biosciences, Inc. | Sequence preserved dna conversion for optical nanopore sequencing |
WO2012138357A1 (en) | 2011-04-04 | 2012-10-11 | President And Fellows Of Harvard College | Nanopore sensing by local electrical potential measurement |
EP3633370B1 (de) | 2011-05-27 | 2024-05-01 | Oxford Nanopore Technologies plc | Verfahren und vorrichtung zur bestimmung der anwesenheit, abwesenheit oder der eigenschaften eines analyten |
US9145623B2 (en) | 2011-07-20 | 2015-09-29 | Thermo Fisher Scientific Oy | Transposon nucleic acids comprising a calibration sequence for DNA sequencing |
EP2737084B1 (de) | 2011-07-25 | 2017-10-18 | Oxford Nanopore Technologies Limited | Haarnadelschleifenverfahren für doppelsträngige polynukleotidsequenzierung mit transmembranporen |
CA2849624C (en) | 2011-09-23 | 2021-05-25 | Oxford Nanopore Technologies Limited | Analysis of a polymer comprising polymer units |
EP2987870B1 (de) | 2011-10-21 | 2020-02-19 | Oxford Nanopore Technologies Limited | Verfahren zur kennzeichnung eines zielpolynukleotids mit einer transmembranpore und eines molekularen motors |
BR112014016112B8 (pt) | 2011-12-29 | 2022-12-27 | Oxford Nanopore Tech Ltd | Método e aparelho para caracterizar um polinucleotídeo alvo, kit, e, uso de uma helicase |
WO2013098561A1 (en) | 2011-12-29 | 2013-07-04 | Oxford Nanopore Technologies Limited | Method for characterising a polynucelotide by using a xpd helicase |
CN107828877A (zh) | 2012-01-20 | 2018-03-23 | 吉尼亚科技公司 | 基于纳米孔的分子检测与测序 |
JP6312607B2 (ja) | 2012-02-16 | 2018-04-18 | ザ リージェンツ オブ ザ ユニバーシティ オブ カリフォルニア | 酵素仲介タンパク質トランスロケーションのためのナノポアセンサー |
WO2013121224A1 (en) | 2012-02-16 | 2013-08-22 | Oxford Nanopore Technologies Limited | Analysis of measurements of a polymer |
AU2013246702B2 (en) | 2012-04-10 | 2017-06-15 | Oxford Nanopore Technologies Limited | Mutant lysenin pores |
US10612083B2 (en) | 2012-04-19 | 2020-04-07 | University Of Washington Through Its Center For Commercialization | Methods and compositions for generating reference maps for nanopore-based polymer analysis |
WO2013185137A1 (en) | 2012-06-08 | 2013-12-12 | Pacific Biosciences Of California, Inc. | Modified base detection with nanopore sequencing |
WO2014013259A1 (en) | 2012-07-19 | 2014-01-23 | Oxford Nanopore Technologies Limited | Ssb method |
BR112015001054A8 (pt) | 2012-07-19 | 2022-08-02 | Oxford Nanopore Tech Ltd | Helicase monomérica que é membro da superfamília 1 ou superfamília 2, helicase hel308 isolada, método para caracterizar um polinucleotídeo alvo, e, sensor para a caracterização de um polinucleotídeo alvo |
AU2013291765C1 (en) | 2012-07-19 | 2019-08-08 | Oxford Nanopore Technologies Limited | Enzyme construct |
GB201313121D0 (en) | 2013-07-23 | 2013-09-04 | Oxford Nanopore Tech Ltd | Array of volumes of polar medium |
CN104936682B (zh) | 2012-10-26 | 2017-12-15 | 牛津纳米孔技术公司 | 微滴界面 |
GB201222928D0 (en) | 2012-12-19 | 2013-01-30 | Oxford Nanopore Tech Ltd | Analysis of a polynucleotide |
CA3157586A1 (en) | 2013-11-26 | 2015-06-04 | Illumina, Inc. | Compositions and methods for polynucleotide sequencing |
EP3120277A1 (de) | 2014-03-21 | 2017-01-25 | Oxford Nanopore Technologies Limited | Analyse eines polymers aus mehrdimensionalen messungen |
GB201408652D0 (en) | 2014-05-15 | 2014-07-02 | Oxford Nanopore Tech Ltd | Model adjustment during analysis of a polymer from nanopore measurements |
US9925679B2 (en) | 2014-05-19 | 2018-03-27 | I+D+M Creative, Llc | Devices and methods for assisting with slicing items |
JP6709213B2 (ja) | 2014-10-16 | 2020-06-10 | オックスフォード ナノポール テクノロジーズ リミテッド | ポリマーの解析 |
-
2013
- 2013-02-18 WO PCT/GB2013/050381 patent/WO2013121224A1/en active Application Filing
- 2013-02-18 EP EP20168745.6A patent/EP3736339B1/de active Active
- 2013-02-18 BR BR112014020211A patent/BR112014020211A2/pt not_active Application Discontinuation
- 2013-02-18 AU AU2013220179A patent/AU2013220179B2/en active Active
- 2013-02-18 KR KR1020147025898A patent/KR102106499B1/ko active IP Right Grant
- 2013-02-18 CA CA2866587A patent/CA2866587C/en active Active
- 2013-02-18 EP EP13706058.8A patent/EP2814980B8/de active Active
- 2013-02-18 US US14/379,242 patent/US20150057948A1/en not_active Abandoned
- 2013-02-18 IN IN6795DEN2014 patent/IN2014DN06795A/en unknown
- 2013-02-18 JP JP2014557122A patent/JP6226888B2/ja active Active
- 2013-02-18 CN CN201380020403.6A patent/CN104321441B/zh active Active
-
2019
- 2019-01-11 US US16/245,306 patent/US11959906B2/en active Active
-
2024
- 2024-02-29 US US18/591,644 patent/US20240264143A1/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6128587A (en) * | 1997-01-14 | 2000-10-03 | The Regents Of The University Of California | Method and apparatus using Bayesian subfamily identification for sequence analysis |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9984147B2 (en) | 2008-08-08 | 2018-05-29 | The Research Foundation For The State University Of New York | System and method for probabilistic relational clustering |
US11921103B2 (en) | 2011-09-23 | 2024-03-05 | Oxford Nanopore Technologies Plc | Method of operating a measurement system to analyze a polymer |
US11959906B2 (en) | 2012-02-16 | 2024-04-16 | Oxford Nanopore Technologies Plc | Analysis of measurements of a polymer |
US10131943B2 (en) | 2012-12-19 | 2018-11-20 | Oxford Nanopore Technologies Ltd. | Analysis of a polynucleotide via a nanopore system |
US11085077B2 (en) | 2012-12-19 | 2021-08-10 | Oxford Nanopore Technologies Ltd. | Analysis of a polynucleotide via a nanopore system |
US11401549B2 (en) | 2014-10-16 | 2022-08-02 | Oxford Nanopore Technologies Plc | Analysis of a polymer |
US10689697B2 (en) | 2014-10-16 | 2020-06-23 | Oxford Nanopore Technologies Ltd. | Analysis of a polymer |
US10783268B2 (en) | 2015-11-10 | 2020-09-22 | Hewlett Packard Enterprise Development Lp | Data allocation based on secure information retrieval |
US11514289B1 (en) * | 2016-03-09 | 2022-11-29 | Freenome Holdings, Inc. | Generating machine learning models using genetic data |
US10823721B2 (en) | 2016-07-05 | 2020-11-03 | Quantapore, Inc. | Optically based nanopore sequencing |
WO2018034745A1 (en) * | 2016-08-18 | 2018-02-22 | The Regents Of The University Of California | Nanopore sequencing base calling |
US11080301B2 (en) | 2016-09-28 | 2021-08-03 | Hewlett Packard Enterprise Development Lp | Storage allocation based on secure data comparisons via multiple intermediaries |
WO2019204720A1 (en) * | 2018-04-20 | 2019-10-24 | The Regents Of The University Of California | Nucleic acid sequencing methods and computer-readable media for practicing same |
CN112703256A (zh) * | 2018-09-04 | 2021-04-23 | 牛津纳米孔科技公司 | 用于测定聚合物序列的方法 |
US20200176082A1 (en) * | 2018-11-28 | 2020-06-04 | Oxford Nanopore Technologies Ltd. | Analysis of nanopore signal using a machine-learning technique |
WO2020236995A1 (en) * | 2019-05-22 | 2020-11-26 | Icahn School Of Medicine At Mount Sinai | Method for de novo detection, identification and fine mapping of multiple forms of nucleic acid modifications |
Also Published As
Publication number | Publication date |
---|---|
EP2814980B1 (de) | 2020-04-15 |
US20240264143A1 (en) | 2024-08-08 |
EP3736339A1 (de) | 2020-11-11 |
CN104321441A (zh) | 2015-01-28 |
JP2015509710A (ja) | 2015-04-02 |
CN104321441B (zh) | 2016-10-19 |
US11959906B2 (en) | 2024-04-16 |
KR102106499B1 (ko) | 2020-05-04 |
AU2013220179B2 (en) | 2018-11-08 |
EP2814980B8 (de) | 2020-06-10 |
IN2014DN06795A (de) | 2015-05-22 |
AU2013220179A1 (en) | 2014-09-04 |
BR112014020211A2 (pt) | 2017-07-04 |
WO2013121224A1 (en) | 2013-08-22 |
KR20140125874A (ko) | 2014-10-29 |
EP2814980A1 (de) | 2014-12-24 |
CA2866587C (en) | 2022-07-26 |
CA2866587A1 (en) | 2013-08-22 |
EP3736339B1 (de) | 2022-07-27 |
US20190154655A1 (en) | 2019-05-23 |
JP6226888B2 (ja) | 2017-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20240264143A1 (en) | Analysis of measurements of a polymer | |
US11085077B2 (en) | Analysis of a polynucleotide via a nanopore system | |
US11401549B2 (en) | Analysis of a polymer | |
Deamer et al. | Three decades of nanopore sequencing | |
JP5328362B2 (ja) | 細孔を使用する方法 | |
US20230039783A1 (en) | Protein and Peptide Fingerprinting and Sequencing by Nanopore Translocation of Peptide-Oligonucleotide Complexes | |
Vilgis et al. | Sequencing in precision medicine | |
CN118120017A (zh) | 纳米孔测量信号分析 | |
Pasquarelli et al. | Genome Sequencing | |
WO2020049293A1 (en) | Method for determining a polymer sequence |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: OXFORD NANOPORE TECHNOLOGIES LIMITED, UNITED KINGD Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:REID, STUART WILLIAM;CLARKE, JAMES ANTHONY;WHITE, JAMES;AND OTHERS;SIGNING DATES FROM 20140828 TO 20140901;REEL/FRAME:033745/0737 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |