Patent Application
20230357838
A Sequence Listing is provided herewith as a Sequence Listing XML, “NEB-451-CIP.xml” created on May 24, 2023, and having a size of 350 KB. The contents of the Sequence Listing XML are incorporated by reference herein in their entirety.
In many organisms, cytosine in the genome can be covalently modified to, for example, 5-methylcytosine (5mC) or 5-hydroxymethylcytosine (5hmC). These epigenetic changes are believed to play a role in a wide variety of phenomena, including gene expression. Global or regional changes of DNA methylation are among the earliest events known to occur in cancer. The identification of methylation profiles in humans is a key step in studying disease processes and is increasingly used for diagnostic purposes.
Current methods for identifying modified cytosine include a deamination step in which cytosines are converted to uracils, leaving the modified cytosines undeaminated. Uracils in these deaminated DNA molecules are copied into thymines during amplification and, after sequencing the amplification products, each of the modified cytosines in the starting sequences can be readily identified as a “C” in the sequenced amplification product, whereas each of the cytosines appear as a “T” in the sequenced amplification product.
DNA may be deaminated chemically (using, e.g., bisulfite; see Frommer et al PNAS 1992 89: 1827-1831) or enzymatically using a DNA deaminase (e.g., APOBEC3A, see, e.g., Sun et al, Genome Res. 2021 31: 291-300 and Vaisvila et al Genome Res. 2021 31: 1280-1289). However, both of these approaches require a single-stranded substrate. As such, current workflows for analyzing modified cytosines typically involve a denaturation step. It would be desirable to eliminate the denaturation step from current workflows. Use of DNA deaminases having particular specificities, such as a bias for deaminating cytosines in a particular sequence context (e.g., the “CpG” context, the most common context for mammalian cytosine methylation) and/or selectivity for deaminating or not deaminating particular modifications, may further simplify such workflows as well as enable other genome analysis and engineering tools.
The present disclosure relates, in some embodiments, to deaminases having one or more desirable properties including, for example, cytosine deaminases that are active on DNA substrates.
These enzymes may deaminate cytosines in a double-stranded DNA substrate (e.g., without denaturing the DNA). Double-stranded DNA deaminases may deaminate cytosines in single-stranded DNA, in addition to deaminating cytosines in double-stranded DNA. Cytosines adjacent to guanines (“CG”) may be deaminated by disclosed deaminases as well as, not as well as, or better than cytosines in other sequence contexts (“CH”, H=A, C, T). Double-stranded DNA deaminase compositions may comprise a deaminase and, optionally, a buffer, one or more enzymes that alter the deamination susceptibility of one or more modified cytosines (e.g., a TET methylcytosine dioxygenase and/or a DNA beta-glucosyltransferase).
The present disclosure relates, in some embodiments, to methods for deaminating double-stranded DNA substrates. For example, deaminating a double-stranded DNA may comprise contacting the double-stranded DNA substrate and a double-stranded DNA deaminase to deaminate cytosines in the double-stranded substrate, for example, without denaturing the substrate or otherwise using any agents that unwind or otherwise separate the strands of the substrate (e.g., a gyrase or a helicase), to produce deamination products. In some embodiments, a double-stranded DNA deaminase may be used to deaminate cytosines in a single-stranded substrate, which may be preceded by separating the strands of the substrate. In some embodiments, methods may include sequencing at least one strand of the product of a deamination reaction (which is a deaminated double-stranded DNA molecule referred to herein as a “deamination product”) to produce sequence reads. A method may include amplifying a deamination product to produce an amplification product and then sequencing the amplification product to produce sequence reads. Disclosed cytosine deaminases may deaminate cytosines without deaminating modified cytosines (e.g., 5mC, 5hmC, 5fC, 5caC, 5ghmC, N4mC) also present in a DNA substrate or may both deaminate cytosines and deaminate one or more modified cytosines in a substrate. Accordingly, the positions of modified cytosines (e.g., 5mC or 5hmC) in a double-stranded DNA substrate can be identified by analysis of sequence reads. Some of the double-stranded DNA deaminases do not deaminate N4mC, but can deaminate other modified cytosines, others do not deaminate 5mC, and 5hmC, others do not deaminate 5hmC but can deaminate 5mC, others do not deaminate 5ghmC but can deaminate 5mC and/or 5hmC, and others that do not deaminate 5fC and 5caC but can deaminate 5mC and 5hmC (see, for example, Table 3). As such, the positions of one or more modified cytosines may be determined in a double-stranded substrate by contacting the substrate with a deaminase having a selected specificity and, optionally, pre-treating the substrate with one or more enzymes that alter the deamination susceptibility of one or more modified cytosines. For example, a method may include pre-treating the double-stranded DNA substrate with: (a) a TET methylcytosine dioxygenase and DNA beta-glucosyltransferase or (b) a TET methylcytosine dioxygenase but not DNA beta-glucosyltransferase. These enzymes modify 5mC and/or 5hmC in double-stranded nucleic acids to make those residues resistant to certain double-stranded DNA deaminases. In some embodiments, a method may include contacting a double-stranded DNA deaminase with a double-stranded nucleic acid not contacted (previously or concurrently) with a TET methylcytosine dioxygenase or a DNA beta-glucosyltransferase, for example, where the double-stranded DNA deaminase does not deaminate 5mC and/or 5hmC. In some embodiments, methods may include base editing and other genome engineering approaches.ln some embodiments, the double-stranded DNA substrate may comprise at least one N4mC or pyrrolo-dC. N4mC is found in prokaryotes and archaea. As such, in some embodiments, a double-stranded DNA substrate may be prokaryotic or archaeal. In some embodiments, a double-stranded DNA substrate may be made by ligating a hairpin adapter to a double-stranded fragment of DNA to produce a ligation product, enzymatically generating a free 3′ end in a double-stranded region of the hairpin adapter in the ligation product, and extending the free 3′ end in a dCTP-free reaction mix that comprises a strand-displacing or nick-translating polymerase, dGTP, dATP, dTTP and modified dCTP. In this method, the modified dCTP is incorporated into the new strand, to produce a double-stranded nucleic acid that has modified Cs. Enzymes and kits for performing the methods are also provided including, for example, a double-stranded DNA deaminase and a reaction buffer.
The file of this patent contains at least one drawing executed in color. Copies of this patent with color drawing(s) will be provided by the Patent and Trademark Office upon request and payment of the necessary fee.
Illumina short-read sequencing-based detection method (Example 3). Total area of the circles corresponds to total activity and the relative sizes of colored sectors show relative activity on the indicated substrates. The inner-most ring shows relative deamination activity on unmodified cytosines in double-stranded DNA (blue sectors) compared to single-stranded DNA (red sectors). The middle ring shows activity on 5-methylated cytosine in double-stranded DNA. The outermost ring shows activity on 5-hydroxymethylated cytosine in double-stranded DNA. Enzyme names are colored according to their phylogenetic family.
The present disclosure provides double-stranded DNA deaminases, variants, ancestors, fusions, compositions, systems, apparatus, methods, and workflows for deaminating double-stranded DNA (in duplex form, without denaturation). Applications of these deaminases include, for example, EM-seq, methyl-SNP-seq, and N4mC detection, among others.
Aspects of the present disclosure can be understood in light of the provided descriptions, figures, sequences, embodiments, section headings, and examples, none of which should be construed as limiting the entire scope of the present disclosure in any way. Accordingly, the innovations set forth herein should be construed in view of the full breadth and spirit of the disclosure.
Each of the individual embodiments described and illustrated herein has discrete components and features which can be readily separated from or combined with the components and/or features of any of the other several embodiments without departing from the scope or spirit of the present teachings. Any recited method can be carried out in the order of events recited or in any other order which is logically possible. Unless otherwise expressly stated to be required herein, each component, feature, and method step disclosed herein is optional and the disclosure contemplates embodiments in which each optional element may be expressly excluded.
Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs. Still, certain terms are defined herein with respect to embodiments of the disclosure and for the sake of clarity and ease of reference.
Sources of commonly understood terms and symbols may include: standard treatises and texts such as Kornberg and Baker, DNA Replication, Second Edition (W.H. Freeman, New York, 1992); Lehninger, Biochemistry, Second Edition (Worth Publishers, New York, 1975); Strachan and Read, Human Molecular Genetics, Second Edition (Wiley-Liss, New York, 1999); Eckstein, editor,
Oligonucleotides and Analogs: A Practical Approach (Oxford University Press, New York, 1991); Gait, editor, Oligonucleotide Synthesis: A Practical Approach (IRL Press, Oxford, 1984); Singleton, et al., Dictionary of Microbiology and Molecular biology, 2d ed., John Wiley and Sons, New York (1994), and Hale & Markham, the Harper Collins Dictionary of Biology, Harper Perennial, N.Y. (1991) and the like. As used herein and in the appended claims, the singular forms “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise. For example, the term “a protein” refers to one or more proteins, i.e., a single protein and multiple proteins. Optional elements may be expressly excluded where exclusive terminology is used, such as “solely,” “only”, in connection with the recitation of the optional elements or when a negative limitation is specified.
Numeric ranges are inclusive of the numbers defining the range. All numbers should be understood to encompass the midpoint of the integer above and below the integer i.e., the number 2 encompasses 1.5-2.5. The number 2.5 encompasses 2.45-2.55 etc. When sample numerical values are provided, each alone may represent an intermediate value in a range of values and together may represent the extremes of a range unless specified.
In the context of the present disclosure, “buffer” and “buffering agent” refer to a chemical entity or composition that itself resists and, when present in a solution, allows such solution to resist changes in pH when such solution is contacted with a chemical entity or composition having a higher or lower pH (e.g., an acid or alkali). Examples of suitable non-naturally occurring buffering agents that may be used in disclosed compositions, kits, and methods include HEPES, MES, MOPS, TAPS, tricine, and Tris. Additional examples of suitable buffering agents that may be used in disclosed compositions, kits, and methods include ACES, ADA, BES, Bicine, CAPS, carbonic acid/bicarbonic acid, CHES, citric acid, DIPSO, EPPS, histidine, MOPSO, phosphoric acid, PIPES, POPSO, TAPS, TAPSO, and triethanolamine.
In the context of the present disclosure, “deaminase substrate” refers to a polynucleotide (e.g., a DNA) molecule that optionally may be exclusively double-stranded, partially double-stranded and partially single-stranded, or exclusively single-stranded. A deaminase substrate may comprise one or more cytosines, one or more modified cytosines, one or more adenines, one or more modified adenines, or combinations thereof. A DNA substrate may comprise one or more adapters. As described in Example 10, such adapters may contain modified nucleotides that are not deaminated during the deamination step. Adapters that do not contain modified nucleotides may be used, so long as base pairing is sufficient to allow the adapters to attach to cognate binding partners as required for a particular method. Additionally, adapters containing modified nucleotides are not required, for example, when the adapters are attached after the deamination step.
In the context of the present disclosure, “double-stranded DNA deaminase” refers to a hydrolyase that deaminates cytosines in double-stranded DNA to uracils and/or deaminates adenines in double-stranded DNA to hypoxanthines. A double-stranded DNA deaminase may deaminate cytosines and/or adenines in double-stranded DNA as well as or better than it deaminates cytosines and/or adenines, respectively, in single-stranded DNA. For example, a double-stranded DNA deaminase may deaminate cytosines double-stranded DNA, but not deaminate cytosines in single-stranded DNA. A double-stranded DNA deaminase may be modification sensitive. For example, a double-stranded DNA deaminase may deaminate an unmodified cytosine or adenine in double-stranded DNA, but not deaminate one or more corresponding modified cytosines or adenines.
In the context of the present disclosure, “duplex” and “double stranded” refer to any conformation of a polynucleotide in which two polynucleotide strands (e.g., separate molecules or spatially separated portions of a single molecule) are arranged anti parallel to one another in a helix with complementary bases of each strand paired with one another (e.g., in Watson-Crick base pairs). Paired bases may be stacked relative to one another to permit pi electrons of the bases to be shared.
Duplex stability, in part, may be related to the ratio of complementary bases to mismatches (if any) in the two strands, ratio of pairs with three hydrogen bonds (e.g., G:C) to pairs with two hydrogen bonds (e.g., A:T, A:U) in the duplex, and the length of the strands with higher ratios and longer strands generally associated with higher stability. Duplex stability, in part, may be related to ambient conditions including, for example, temperature, pH, salinity, and/or the presence, concentration and identity of any buffer(s), denaturant(s) (e.g., formamide), crowding agent(s) (e.g., PEG), detergent(s) (e.g., SDS), surfactant(s), polysaccharide(s) (e.g., dextran sulfate), chelator(s) (e.g., EDTA), and nucleic acid(s) (e.g., salmon sperm DNA). A duplex polynucleotide may comprise one or more unpaired bases including, for example, a mismatched base, a hairpin loop, a single-stranded (5′ and/or 3′) end.
Duplex polynucleotides (e.g., double-stranded DNA deaminase substrates) may have any desired length. For example, a duplex polynucleotide may have a length of 50 nucleotides, 10-200 nucleotides, 80-400 nucleotides, 50-500 nucleotides, 500 nucleotides, 1 kb, 2 kb, 5 kb or 10 kb.
Duplex polynucleotides may have any desired number of mismatched or unpaired nucleotides, for example, 1 per 100 nucleotides, 2 per 100 nucleotides, 3 per 100 nucleotides, 5 per 100 nucleotides, or 10 per 100 nucleotides.
In the context of the present disclosure, “fusion protein” refers to a protein composed of two or more polypeptide components that are un-joined in their native state. Fusion proteins may be a combination of two, three or four or more different proteins. For example, a fusion protein may comprise two naturally occurring polypeptides that are not joined in their respective native states. A fusion protein may comprise two polypeptides, one of which is naturally occurring and the other of which is non-naturally occurring. The term polypeptide is not intended to be limited to a fusion of two heterologous amino acid sequences. A fusion protein may have one or more heterologous domains added to the N-terminus, C-terminus, and or the middle portion of the protein. If two parts of a fusion protein are “heterologous”, they are not part of the same protein in its natural state. Examples of fusion proteins include proteins comprising a double-stranded DNA deaminase fused to a protein such as albumin, another enzyme (e.g., an endonuclease), an antibody, a binding domain suitable for immobilization such as maltose binding domain (MBP), a histidine tag (“His-tag”), a chitin binding domain, an alpha mating factor or a SNAP-Tag® (New England Biolabs, Ipswich, MA (see for example U.S. Pat. Nos. 7,939,284 and 7,888,090)), a DNA-binding domain (e.g., the DNA binding domain of a transcription factor, a non-specific DNA-binding domain (e.g., Sso7d), or a specific DNA binding domain (e.g., BD09; see, for example, U.S. Pat. No. 9,963,687), or a methyl binding domain (MBD), with the deaminase optionally positioned closer to the N-terminus or closer to the C-terminus than the other component(s). A binding peptide may be used to improve solubility or yield of the deaminase during the production of the protein reagent. Other examples of fusion proteins include fusions of a deaminase and a heterologous targeting sequence, a linker, an epitope tag, a detectable fusion partner, such as a fluorescent protein, (3-galactosidase, luciferase and/or functionally similar peptides. Components of a fusion protein may be joined by one or more peptide bonds, disulfide linkages, and/or other covalent bonds.
In the context of the present disclosure, “modified cytosine” refers to any covalent modification of cytosine including naturally occurring and non-naturally occurring modifications. Modified cytosines include, for example, 1-methylcytosine (1mC), 2-O-methylcytosine (m2C), 3-ethylcytosine (e3C), 3,N4 -ethylenocytosine (SC), 3-methylcytosine (3mC), 4-methylcytosine (4mC), 5-carboxylcytosine (5CaC), 5-formylcytosine (5fC), 5-hydroxymethylcytosine (5hmC), 5-methylcytosine (5mC), N4 -methylcytosine (N4mC), and pyrrolo-cytosine (pyrrolo-C). 5-carboxylcytosine (5caC) is the final oxidized derivative of 5-methylcytosine (5mC). 5mC is oxidized to 5-hydroxymethylcytosine (5hmC) which is then oxidized to 5-formylcytosine (5fC) then 5caC. Additional examples of modified nucleotides may be found at https://dnamod.hoffmanlab.org.
In the context of the present disclosure, “non-naturally occurring” refers to a polynucleotide, polypeptide, carbohydrate, lipid, or composition that does not exist in nature. Such a polynucleotide, polypeptide, carbohydrate, lipid, or composition may differ from naturally occurring polynucleotides polypeptides, carbohydrates, lipids, or compositions in one or more respects. For example, a polymer (e.g., a polynucleotide, polypeptide, or carbohydrate) may differ in the kind and arrangement of the component building blocks (e.g., nucleotide sequence, amino acid sequence, or sugar molecules). A polymer may differ from a naturally occurring polymer with respect to the molecule(s) to which it is linked. For example, a “non-naturally occurring” protein may differ from naturally occurring proteins in its secondary, tertiary, or quaternary structure, by having a chemical bond (e.g., a covalent bond including a peptide bond, a phosphate bond, a disulfide bond, an ester bond, and ether bond, and others) to a polypeptide (e.g., a fusion protein), a lipid, a carbohydrate, or any other molecule. Similarly, a “non-naturally occurring” polynucleotide or nucleic acid may contain one or more other modifications (e.g., an added label or other moiety) to the 5′- end, the 3′ end, and/or between the 5′- and 3′-ends (e.g., methylation) of the nucleic acid. A “non-naturally occurring” composition may differ from naturally occurring compositions in one or more of the following respects: (a) having components that are not combined in nature; (b) having components in concentrations not found in nature; (c) omitting one or components otherwise found in naturally occurring compositions; (d) having a form not found in nature, e.g., dried, freeze dried, crystalline, aqueous; and (e) having one or more additional components beyond those found in nature (e.g., buffering agents, a detergent, a dye, a solvent or a preservative).
With reference to an amino acid, “position” refers to the place such amino acid occupies in the primary sequence of a peptide or polypeptide numbered from its amino terminus to its carboxy terminus. A position in one primary sequence may correspond to a position in a second primary sequence, for example, where the two positions are opposite one another when the two primary sequences are aligned using an alignment algorithm (e.g., BLAST (Journal of Molecular Biology. 215 (3): 403-410) using default parameters (e.g., expect threshold 0.05, word size 3, max matches in a query range 0, matrix BLOSUM62, Gap existence 11 extension 1, and conditional compositional score matrix adjustment) or custom parameters). An amino acid position in one sequence may correspond to a position within a functionally equivalent motif or structural motif that can be identified within one or more other sequence(s) in a database by alignment of the motifs. Analogously, with reference to a nucleotide, “position” refers to the place such nucleotide occupies in the nucleotide sequence of an oligonucleotide or polynucleotide numbered from its 5′ end to its 3′ end.
All publications, patents, and patent applications mentioned in this specification are herein incorporated by reference to the same extent as if each individual publication, patent, or patent application was specifically and individually indicated to be incorporated by reference. Reagents referenced in this disclosure may be made using available materials and techniques, obtained from the indicated source, and/or obtained from New England Biolabs, Inc. (Ipswich, MA).
Double-stranded DNA Deaminases The present disclosure relates to naturally occurring and non-naturally occurring double-stranded DNA deaminases. A non-naturally occurring double-stranded DNA deaminase may relate to, but differ from, a naturally occurring protein. Naturally-occurring proteins often include a deaminase as a single domain of a larger, multi-domain structure with the deaminase domain positioned at the most C-terminal end. Non-naturally occurring double-stranded DNA deaminases may constitute truncated versions of a naturally-occurring protein, in which cases, the non-naturally occurring double-stranded
DNA deaminases may have a high degree of identity to a portion of a naturally-occurring sequence, but lack, for example, structural and/or functional domains or sub-units of the corresponding naturally-occurring proteins. A non-naturally occurring double-stranded DNA deaminase may have any number of insertions, deletions, or substitutions relative to a naturally occurring enzyme. For example, a non-naturally occurring double-stranded DNA deaminase may have less than 100% identity, less than 99% identity, less than 98% identity, less than 90% identity, less than 85% identity, less than 80% identity, less than 70% identity, less than 60% identity, less than 50% identity, less than 40% identity, less than 30% identity, or less than 20% identity to a naturally occurring enzyme. Non-naturally occurring double-stranded DNA deaminases may include expression and/or purification tags. Non-naturally occurring double-stranded DNA deaminase disclosed herein may have an amino acid sequence that is at least 80% identical (e.g., at least 90% identical, at least 95% identical or at least 98% identical or at least 99% identical to) the C-terminal deaminase domain of a naturally-occurring protein, wherein the double-stranded DNA deaminase possesses a double-stranded DNA deaminase activity and does not comprise the N-terminus of the corresponding naturally-occurring protein (if any). In some embodiments, a non-naturally occurring double-stranded DNA deaminase lacks at least 10, at least 20, at least 50 or at least 100 of the N-terminal amino acids of the corresponding naturally-occurring protein. In some embodiments, a double-stranded DNA deaminase is no more than 300 amino acids in length, e.g., no more than 200 amino acids in length or no more than 150 amino acids in length.
According to some embodiments, a double-stranded DNA deaminase may comprise an amino acid sequence having at least 80%, at least 85%, at least 88% identical, at least 90%, at least 92%, at least 93%, at least 95%, at least 96%, at least 97%, at least 98% or at least 99% identity to any of SEQ ID NOS: 1-152. In some embodiments, a double-stranded DNA deaminase may be encoded by a nucleic acid sequence that, when transcribed, translated, and/or processed, results in an amino acid sequence having at least 80%, at least 85%, at least 90%, at least 93%, at least 96%, at least 97%, at least 98% or at least 99% identity to any of SEQ ID NOS: 1-152. A double-stranded DNA deaminase may have an amino acid sequence at least 90% (e.g., at least 95%, at least 98%, at least 99%) identical to any of SEQ ID NOS:
In some embodiments, a double-stranded DNA deaminase may be a fusion protein. For example, a double-stranded DNA deaminase may have a purification tag (e.g., a His tag or the like) at either end. In some embodiments, a double-stranded DNA deaminase may be fused to a DNA binding protein (e.g., the DNA binding domain of a transcription factor) or the protein component of a nucleic acid-guided endonuclease (e.g., a catalytically dead Cas9 (dCas9) or a Cas9 nickase (nCas9) or TALEN (transcription activator-like effector nucleases)) so that the fusion protein can affect site-specific C to T substitutions in a genome. Example methods of “base editing” are described in, for example, Komor et al (Nature 533: 420-424), among other publications.
A double-stranded DNA deaminase optionally may deaminate cytosine, but not adenine (a “dsDNA cytosine deaminase”), deaminate adenine, but not cytosine (a “dsDNA adenine deaminase”), or deaminase both adenine and cytosine (appreciating that one may be a better substrate than the other under otherwise equivalent conditions). A double-stranded DNA deaminase may be modification sensitive. For example, a double-stranded DNA deaminase may deaminate cytosine, but not deaminate one or more modified cytosines in double stranded DNA. For example, a double-stranded DNA deaminase may deaminate cytosine, but not deaminate 5mC or N4mC or it may deaminate C and 5mC, but not 5hmC, 5ghmC or N4mC.
The present disclosure provides double-stranded DNA deaminase compositions including, for example, reaction mixtures. According to some embodiments, deaminase compositions may comprise (a) a double-stranded DNA deaminase and (b) a double-stranded DNA. A deaminase composition may comprise, for example, a deaminase variant (e.g., having an amino acid sequence at least 80% identical to one or more of SEQ ID NOS:1-152). A double-stranded DNA deaminase composition may be free of one or more other catalytic activities. For example, a double-stranded DNA deaminase composition may be free of nucleases that cleave dsDNA, free of nucleases that cleave ssDNA, free of polymerase activity, free of DNA modification activity, and/or free of protease activity, in each case, under desired test conditions (e.g., conditions of time, temperature, pH, salinity, model substrate and/or others), for example, conditions intended to replicate conditions of a specific use of the double-stranded DNA deaminase composition or intended to represent conditions for a range of uses.
In some embodiments, double-stranded DNA deaminases and compositions comprising one or more double-stranded DNA deaminase may have any desirable form including, for example, a liquid, a gel, a film, a powder, a cake, and/or any dried or lyophilized form. A double-stranded DNA deaminase composition may comprise a double-stranded DNA deaminase and a support or matrix, for example, a film, gel, fabric, or bead comprising, for example, a magnetic material, agarose, polystyrene, polyacrylamide, and/or chitin.
In some embodiments, a reaction mix may comprise: a double-stranded DNA substrate that comprises cytosines and a double-stranded DNA deaminase. A double-stranded DNA substrate may comprise cytosines and at least one modified cytosine, e.g., a 5fC, 5CaC, 5mC, 5hmC, N4mC or pyrrolo-C. A double-stranded DNA substrate may be eukaryotic DNA (e.g., plant or animal) or bacterial. In some embodiments, the double-stranded DNA substrate may be mammalian, e.g., from a human. In some embodiments, the double-stranded DNA substrate may be human cfDNA. The reaction mix may additionally comprise one or more of a TET methylcytosine dioxygenase (e.g., TET2) and a DNA beta-glucosyltransferase, as described herein and/or a ligase, a polymerase, a proteinase K, and/or a thermolabile proteinase K. A reaction mix may be free of unwinding agents (e.g., gyrases, topoisomerases, single-stranded DNA binding proteins, or helicases) and/or free of denaturants.
The present disclosure provides methods for identifying the type and/or position of modified nucleotides in, for example, DNA using a deaminase. In some embodiments, a method may comprise providing a double-stranded DNA substrate of any desired length. For example, a double-stranded DNA substrate may have a length of 50 nucleotides, 10-200 nucleotides, 80-400 nucleotides, 50-500 nucleotides, 500 nucleotides, 1 kb, 2 kb, 5 kb or 10 kb. A double-stranded DNA substrate, in some embodiments, may be a fragment of genomic DNA, organelle DNA, cDNA, or other DNAs of interest and can be or arise from any desired source (e.g., human, non-human mammal, plants, insects, microbial, viral, or synthetic DNA). A DNA substrate may be prepared, in some embodiments by extracting (e.g., genomic DNA) from a biological sample and, optionally, fragmenting it. In some embodiments, fragmenting DNA may comprise mechanically fragmenting the DNA (e.g., by sonication, nebulization, or shearing) or enzymatically fragmenting the DNA (e.g., using a double stranded DNA “dsDNA” fragmentation mix). Examples of enzymes for fragmentation include NEBNext® Fragmentase®, Ultrashear, and FS systems (New England Biolabs, Ipswich MA)), among others. In some embodiments, DNA for deamination may already be fragmented (e.g., as is the case for FFPE samples and circulating cell-free DNA (cfDNA)).
According to some embodiment, a method may include polishing DNA ends (e.g., the ends of fragmented DNA). For example, DNA ends may be contacted with (a) a proofreading polymerase to excise 3′ overhanging nucleotides, if any, (b) a proofreading and/or non-proofreading polymerase to fill in 5′ overhangs, if any, and/or (c) a polynucleotide kinase (PNK) to phosphorylate unphosphorylated 5′ ends, if any. In some embodiments, a method may comprise contacting DNA ends (e.g., blunt ends) with a non-proofreading polymerase to add an untemplated A-tail (e.g., a single base overhang comprising adenine) to the 3′ end. Methods may include, according to some embodiments, ligating one or more adapters to DNA ends. Adapters may comprise one or more sample tags, unique molecular identifiers (UMIs), modified nucleotides, primer sequences (e.g., for sequencing). In some embodiments, adapters may comprise cytosines (or adenines) that are not substrates for the deaminase to be used. If desired, polishing products and/or ligation products may be cleaned up, for example, to separate polishing products or ligation products, as applicable, from enzymes, unreacted nucleotides and/or adapters.
In some embodiments, a method may comprise contacting (a) a deaminase substrate and (b) a glucosyltransferase (e.g., T4-BGT) and/or Ten-eleven translocation (TET) dioxygenase to produce a modified deaminase substrate. BGT may glucosylate 5hmC to form 5ghmC. TET may oxidize 5mC to 5caC. If subsequently treated with sodium bisulfite or Apolipoprotein B mRNA editing enzyme subunit 3A (APOBEC3A), all Cs except 5ghmC in the modified deaminase substrate would be deaminated. Deaminases disclosed herein may obviate the need to denature the DNA prior to deamination (e.g., as with APOBEC3A) and may provide methylation sensitivities.
A method may comprise contacting a double-stranded DNA substrate that comprises cytosines and a double-stranded DNA deaminase to produce a deamination product that comprises deaminated cytosines. A double-stranded DNA substrate may further comprise one or more modified cytosines, e.g., one or more modified cytosines selected from 5fC, 5CaC, 5mC, 5hmC, N4mC and pyrrolo-C, 4mC, EC, 3mC, e3C, m2C, and 1mC. A double-stranded DNA deaminase substrate does not need to be denatured before or during deamination. As such, methods can be practiced in the absence of a denaturation step. In some embodiments, deamination methods may comprise contacting a double-stranded DNA substrate comprising cytosines and a double-stranded DNA deaminase to produce a reaction mix to produce a deamination product comprising deaminated cytosines.
Deamination methods may further comprise amplifying the deamination product to produce an amplification product, thereby copying any deaminated Cs in the original strand to Ts in the amplification product. Deamination methods may further comprise ligating an asymmetric (or “Y”) adapter, e.g., an Illumina P5/P7 adapter, onto the deamination product and amplifying the deaminated product using primers complementary to sequences in the adapter. In some embodiments, a method may comprise sequencing a deamination product, or amplifying a deamination product to produce amplification products and sequencing the amplification products, in each case, to produce sequence reads. Deamination products and/or amplification products may be sequenced using any suitable system including Illumina's reversible terminator method (see, e.g., Shendure et al, Science 2005 309: 1728). In some embodiments, a deaminated product may be sequenced directly, without amplification, for example, by nanopore or PacBio sequencing. A sequencing step may result in at least 10,000, at least 100,000, at least 500,000, at least 1M, at least 10M, at least 100M, at least 1B or at least 106 sequence reads per reaction. In some cases, the reads may be paired-end reads. A method may comprise analyzing sequence reads to identify a modified cytosine in the double-stranded DNA substrate, where a modified cytosine can be identified as a “C” because it is deaminase-resistant.
Double-stranded DNA deaminases that are “blocked” by or do not deaminate modified cytosines (e.g., 5mC, 5hmC, 5ghmC, N4mC) may be used in a variety of “EM-seq”-like workflows for the analysis of modified cytosines (e.g., see
Workflows for example deamination methods are shown in
As illustrated, the double-stranded DNA deaminase can be added to the reaction without any clean-up, denaturation or addition of unwinding agents.
As illustrated in
As illustrated in
In some embodiments, a double-stranded DNA substrate may comprise at least one N4mC (N4-methyl-cytosine) which is a cytosine modification that is resistant to some double-stranded DNA deaminases. Double-stranded DNA deaminases useful for detecting N4mC may have an amino acid sequence that is at least 90% identical to the amino acid sequence of any of SEQ ID NOS:1-28. For example, double-stranded DNA deaminases useful for detecting N4mC may have an amino acid sequence that is at least 90% identical to the amino acid sequence of any of CseDa01 (SEQ ID NO:3) and LbDa01 (SEQ ID NO:19) double-stranded DNA deaminases. In these embodiments, the double-stranded DNA substrate may be or comprise prokaryotic or archaeal DNA.
In some embodiments, the double-stranded DNA deaminase may be used in a “methyl-SNP-seq” workflow (see, e.g., Yan et al, Genome Res. 2022; gr.277080.122). For example, a method may comprise; (a) ligating a hairpin adapter to a double-stranded fragment of DNA to produce a ligation product, (b) enzymatically generating a free 3′ end in a double-stranded region of the hairpin adapter in the ligation product; and (c) extending the free 3′ end in a dCTP-free reaction mix that comprises a strand-displacing or nick-translating polymerase, dGTP, dATP, dTTP and modified dCTP to produce the double-stranded DNA substrate, as described in U.S. Provisional Application Ser. No. 63/399,970, filed on Aug. 22, 2022, which application is incorporated by reference herein. Examples of modified dCTPs include 5mdCTP, pyrrolo-dCTP, and N4mdCTP among other modified dCTPs that can be incorporated by a polymerase. Deaminases may have an amino acid sequence that is at least 90% identical to the amino acid sequence of any of MGYPDa20 (SEQ ID NO: 11), NsDa01 (SEQ ID NO: 27), AshDa01 (SEQ ID NO: 40).
Other current workflows for sequencing genomic DNA employ multi-copy sequencing strategies to simultaneously obtain both genetic bases (e.g., G, A T, C) and epigenetic bases (e.g., methylated bases), for example to better distinguish between genetic C-to-T mutations and epigenetically modified cytosine. These workflows (e.g., Yan et al, ibid, and Fullgrabe, et al., Nat Biotechnol (2023). https://doi.org/10.1038/s41587-022-01652-0) involve copying strands of genomic DNA fragments to generate two copies of the strand on a single DNA strand. The multi-copy strands are then deaminated, followed by sequencing. This results in obtaining two sequence reads of the same strand of the genomic DNA fragment (and if desired four reads because both strands of the genomic DNA fragment may be subjected to this process); bioinformatics tools are used to discern whether Cs arose from modified C or from a mutation in the genomic DNA, as well as to identify errors arising from sequencing or amplification steps. Such workflows involve linking together the two strands of the genomic DNA, e.g., using a hairpin; breaking that linkage to synthesize the copy, thereby creating the multi-copy strand. Then, in either order, deaminating the multi-copy strand, and adding sequencing primers to the multi-copy strands to obtain reads of the original and copied sequences. The sequences are determined using rules based on the selected deamination process. A double-stranded DNA deaminase described herein may be used to reduce the complexity of such workflows (see, e.g, Example 10). A double-stranded DNA deaminase described herein may also be used to sequence genetic and epigenetic bases using a standard sequencing workflow by adding a deamination step, without the need for making a multi-copy strand (e.g., see
According to some embodiments, a double-stranded DNA deaminase composition may comprise a double-stranded DNA deaminase and, optionally, any of (including one or more of) a buffering agent (e.g., a storage buffer, a reaction buffer), an excipient, a salt (e.g., NaCl, MgCll2, CaCl2), a protein (e.g., albumin, an enzyme), a stabilizer, a detergent (for example, ionic, non-ionic, and/or zwitterionic detergents (e.g., octoxinol, polysorbate 20)), a polynucleotide, a cell (e.g., intact, digested, or any cell-free extract), a biological fluid or secretion (e.g., mucus, pus), an aptamer, a crowding agent, a sugar (e.g., a mono, di, tri, tetra, or higher saccharide), a starch, cellulose, a glass-forming agent (e.g., for lyophilization), a lipid, an oil, aqueous media, a support (e.g., a bead) and/or (non-naturally occurring) combinations thereof. Combinations may include for example, two or more of the listed components (e.g., a salt and a buffer) or a plurality of a single listed component (e.g., two different salts or two different sugars). Examples of proteins that may be included in a double-stranded DNA deaminase composition include one or more enzymes that alter the deamination susceptibility of one or more modified cytosines (e.g., a TET methylcytosine dioxygenase and/or a DNA beta-glucosyltransferase).
The present disclosure relates, in some embodiments, to a deaminase kit comprising a double-stranded DNA deaminase. A kit may comprise any of the components described herein. A double-stranded DNA deaminase composition or kit may include, for example, double-stranded DNA deaminase and, optionally, a storage buffer (e.g., comprising a buffering agent and comprising or lacking glycerol), and/or a reaction buffer. A reaction buffer for a deaminase composition or a deaminase kit may be in concentrated form, and the buffer may include one or more additives (e.g., glycerol), one or more salts (e.g. KCl), one or more reducing agents, EDTA, one or more detergents, one or more non-ionic surfactants, one or more ionic (e.g. anionic or zwitterionic) surfactants, and/or crowding agents. A kit comprising dNTPs may include one, two, three of all four of dATP, dTTP, dGTP and dCTP. A kit may further comprise one or more modified nucleotides.
One or more components of a kit may be included in one container for a single step reaction, or one or more components may be contained in one container, but separated from other components for sequential use or parallel use. For example, a kit may comprise two components in a single tube (e.g., a deaminase and a storage buffer) and all other components in separate, individual tubes, in each case, with the contents provided in any desired form (e.g., liquid, dried, lyophilized). One tube in a kit may contain a mastermix, for example, for receiving and amplifying a DNA (e.g., a deaminated DNA). For example, a double-stranded DNA deaminase may be deposited in the cap of a tube while components for transcribing a template nucleic acid are deposited in the body of the tube. As desired, for example, upon completion of the deamination reaction, the tube may be tapped, shaken, turned, spun, or otherwise moved to contact the deposited double-stranded DNA deaminase with the deamination reaction mixture. A kit may include a double-stranded DNA deaminase and the reaction buffer in a single tube or in different tubes and, if included in a single tube, the double-stranded DNA deaminase and the buffer may be present in the same or separate locations in the tube. For example, a kit may comprise a double-stranded DNA deaminase, as described above, and a reaction buffer (e.g., a 5× or 10× buffer). The contents of a kit may be formulated for use in a desired method or process. In some embodiments, the kit may further comprise (a) a TET methylcytosine dioxygenase (e.g., TET2) and a DNA beta-glucosyltransferase or (b) a TET methylcytosine dioxygenase and no DNA beta-glucosyltransferase. In some embodiments, a kit does not contain either a TET methylcytosine dioxygenase or DNA beta-glucosyltransferase. In some embodiments, a kit further comprises a modified dCTP selected from 5hmdCTP, 5fdCTP, 5cadCTP, 5mdCTP, pyrrolo-dCTP and N4mdCTP and/or a strand-displacing or nick translating polymerase. In some embodiments, a kit may additionally comprise a ligase, a polymerase, a proteinase K, and/or a thermolabile proteinase K. A double-stranded DNA deaminase may be lyophilized or in a buffered storage solution that contains glycerol.
As would be apparent to those having the benefit of the present disclosure, a double-stranded DNA deaminase may be used in a variety of genome analysis methods, particularly methods whose goal is to identify the position and/or identity of one or more modified cytosines and/or determine the methylation status of a cytosine. In other embodiments, a double-stranded DNA deaminase can be a component of a fusion protein for based editing, i.e., generating site-specific C to T substitutions in a genome.
The present disclosure further relates to embodiments disclosed in U.S. Provisional Application No. 63/264,513 including all of the following:
Embodiment 1. A polypeptide comprising at least 90% sequence identity with any of SEQ ID NOs: 1-8, not including 100% identity to SEQ ID NO: 3.
Embodiment 2. The polypeptide according to embodiment 1, comprising at least 90% sequence identity with any of SEQ ID NOs: 1-3 not including 100% identity to SEQ ID NO: 3.
Embodiment 3. The polypeptide according to embodiment 1, comprising at least 90% sequence identity with any of SEQ ID NOs: 1 or 2.
Embodiment 4. The polypeptide according to any of embodiments 1-3, capable of deaminating cytosine in double stranded DNA (dsDNA) with no sequence bias.
Embodiment 5. The polypeptide according to any of embodiments 1-3, capable of deaminating cytosine in single stranded DNA (ssDNA) with no sequence bias.
Embodiment 6. The polypeptide of any of embodiments 1-5, comprising a fusion protein.
Embodiment 7. The polypeptide of any of embodiments 1-6, wherein the polypeptide is lyophilized.
Embodiment 8. The polypeptide of any of embodiments 1-7, wherein the polypeptide is immobilized on a substrate.
Embodiment 9. The polypeptide of any of embodiments 1-8, wherein the polypeptide is combined with one or more reagents in a mixture wherein one or more reagents in the mixture comprises a second polypeptide.
Embodiment 10. The polypeptide of embodiment 9, wherein the second polypeptide is selected from the group consisting of a ligase, a polymerase, a methylcytosine (mC) dioxygenase, DNA glucosyltransferase, a Proteinase K, and a Thermolabile Proteinase K.
Embodiment 11. The polypeptide of any of embodiments 9-10, wherein the one or more reagents in the mixture further comprises a reversible inhibitor of the deaminase.
Embodiment 12. The polypeptide of any of embodiments 1-11, wherein the mixture further comprises DNA.
Embodiment 13. A method for methylome analysis comprising
Embodiment 14. The method according to embodiment 13, wherein prior to (a) adding to the reaction mixture, a methylcytosine (mC) dioxygenase to the genomic DNA for converting mC to hydroxymethylcytosine (hmC).
Embodiment 15. The method according to any of embodiments 13-14, wherein prior to (a) adding a hydroxymethylcytosine (hmC) modifying reagent to the reaction mixture.
Embodiment 16. The method according to any of embodiments 13-15, wherein (b) further comprises inactivating the DNA deaminase with a Proteinase K or Thermolabile Proteinase K.
Embodiment 17. The method according to any of embodiments 13-16, wherein (b) further comprises amplifying the DNA containing the converted cytosines.
Embodiment 18. The method according to any of embodiments 13-17, further comprising sequencing the amplified DNA.
Embodiment 19. The method according to any of embodiments 13-18, further comprising determining the location of methylcytosine (mC) in genomic DNA.
Embodiment 20. A kit comprising a deaminase capable of deaminating cytosine in double stranded DNA (dsDNA) and optionally single stranded DNA (ssDNA) with no sequence bias.
Embodiment 21. The kit according to embodiment 20, further comprising a methyl dioxygenase in a separate container from the dixoygenase.
Embodiment 22. The kit according to embodiment 20 or 21, further comprising a hydroxymethylcytosine (hmC) modifying enzyme in the same container with the dioxygenase or in a different container.
Embodiment 23. A method for deaminating a double-stranded nucleic acid, the method comprising:
contacting:
Embodiment 24. The method according to Embodiment 23, wherein the double-stranded DNA substrate further comprises a modified cytosine.
Embodiment 25. The method according to Embodiment 24, wherein the modified cytosine is a 5fC, 5CaC, 5mC, 5hmC, N4mC, 5ghmC, or pyrrolo-C.
Embodiment 26. The method according to Embodiment 23, wherein the method further comprises:
Embodiment 27. The method according to Embodiment 26, wherein the method further comprises:
Embodiment 28. The method according to Embodiment 23, wherein the double-stranded DNA substrate is eukaryotic or bacterial DNA.
Embodiment 29. The method according to Embodiment 23, wherein the double-stranded DNA substrate is human cfDNA.
Embodiment 30. The method according to Embodiment 23, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 90% identical to any of SEQ ID NOS: 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 15, 16, 19, 24, 26, 27, 28, 33, 40, 49, 50, 63, 95, 96, 97, and 99.
Embodiment 31. The method according to Embodiment 23, wherein the double-stranded DNA substrate is pre-treated with a TET methylcytosine dioxygenase and DNA beta-glucosyltransferase.
Embodiment 32. The method according to Embodiment 31, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 90% identical to any of the SEQ ID NOS for MGYPDa829 (SEQ ID NO: 96), MGYPDa06 (SEQ ID NO: 4), CrDa01 (SEQ ID NO: 12), AvDa02 (SEQ ID NO: 2), CsDa01 (SEQ ID NO: 9), LbsDa01 (SEQ ID NO: 10), FIDa01 (SEQ ID NO: 8), MGYPDa26 (SEQ ID NO: 7), MGYPDa23 (SEQ ID NO: 6), chimera_10 (SEQ ID NO: 97) and AncDa04 (SEQ ID NO: 95).
Embodiment 33. The method according to Embodiment 23, wherein the double-stranded DNA substrate is pre-treated with a TET methylcytosine dioxygenase but not DNA beta-glucosyltransferase.
Embodiment 34. The method according to Embodiment 33, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 90% identical to any of the SEQ ID NOS for CseDa01 (SEQ ID NO: 3) and LbDa02 (SEQ ID NO: 1).
Embodiment 35. The method according to Embodiment 23, wherein the double-stranded DNA substrate is not pre-treated with either a TET methylcytosine dioxygenase or DNA beta-glucosyltransferase.
Embodiment 36. The method according to Embodiment 23, wherein the double-stranded DNA substrate comprises at least one N4mC.
Embodiment 37. The method according to Embodiment 36, wherein the double-stranded DNA substrate is bacterial DNA.
Embodiment 38. The method according to Embodiment 36, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 90% identical to any of the SEQ ID NOS for MGYPDa20 (SEQ ID NO: 11), NsDa01 (SEQ ID NO: 27), and AshDa01 (SEQ ID NO: 40).
Embodiment 39. The method according to Embodiment 23, further comprising:
Embodiment 40. The method according to Embodiment 39, wherein the modified dCTP is 5mdCTP, pyrrolo-dCTP, 5hmdCTP or N4-mdCTP.
Embodiment 41. The method according to Embodiment 39, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 90% identical to any of the SEQ ID NOS for MGYPDa20 (SEQ ID NO: 11), NsDa01 (SEQ ID NO: 27), AshDa01 (SEQ ID NO:40).
Embodiment 42. An enzyme comprising an amino acid sequence that is at least 80% identical to the C-terminal deaminase domain of a naturally-occurring protein, wherein the enzyme:
Embodiment 43. The enzyme according to Embodiment 42, wherein the enzyme is no more than 300 amino acids in length.
Embodiment 44. The enzyme according to Embodiment 42, wherein the enzyme is at least 80% identical to any of SEQ ID NOS: 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 15, 16, 19, 24, 26, 27, 28, 33, 40, 49, 50, 63, 95, 96, 97, and 99.
Embodiment 45. The enzyme according to Embodiment 42, wherein the enzyme is fused with a catalytically dead Cas9 (dCas9) or a nicking Cas9 (nCas9) or Transcription activator-like effector nucleases (TALEN).
Embodiment 46. A kit comprising:
Embodiment 47. The kit according to Embodiment 46, wherein the kit further comprises:
Embodiment 48. The kit according to Embodiment 46, wherein the kit is free of TET methylcytosine dioxygenase and DNA beta-glucosyltransferase.
Embodiment 49. The kit according to Embodiment 46, wherein the kit further comprises a modified dCTP selected from 5mdCTP, pyrrolo-dCTP, 5hmdCTP and N4-mdCTP.
Embodiment 50. A reaction mix comprising:
Embodiment 51. The reaction mix according to Embodiment 50, wherein the double-stranded DNA substrate comprises cytosines and at least one modified cytosine.
Embodiment 52. The reaction mix according to Embodiment 50, wherein the modified cytosine is a 5fC, 5caC, 5mC, 5hmC, N4mC or pyrrolo-C.
Embodiment 53. The reaction mix according to Embodiment 50, wherein the double-stranded DNA substrate comprises eukaryotic or bacterial DNA.
Embodiment 54. The reaction mix according to Embodiment 50, wherein the double-stranded DNA substrate is human cfDNA.
Embodiment 55. The reaction mix according to Embodiment 50, wherein the deaminase has an amino acid sequence that is at least 90% identical to any of SEQ ID NOS: 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 15, 16, 19, 24, 26, 27, 28, 33, 40, 49, 50, 63, 95, 96, 97, and 99.
Embodiment 56. A method for sequencing, comprising: contacting a single-stranded DNA substrate comprising a genomic DNA fragment with a double-stranded DNA deaminase to produce a deamination product; sequencing the deamination product, or amplifying the deamination product to produce amplification products and sequencing the amplification products, in each case, to produce sequence reads, wherein the double-stranded DNA deaminase is an enzyme of Embodiment 42.
Embodiment 57. A method for sequencing, comprising:
Embodiment 58. The method of Embodiment 57, wherein the double-stranded DNA deaminase has sequence bias for cytosine in a CpG context.
Embodiment 59. The method of Embodiment 58, wherein the double-stranded DNA deaminase is modification sensitive.
Embodiment 60. The method of Embodiment 59, wherein the double-stranded DNA deaminase does not deaminate one or more of 5fC, 5CaC, 5mC, 5hmC, N4mC, or 5ghmC.
Embodiment 61. The method of Embodiment 58, wherein the double-stranded DNA deaminase is not modification sensitive.
Embodiment 62. The method of Embodiment 58, wherein the double-stranded DNA substrate or the genomic fragment is not pre-treated with either a TET methylcytosine dioxygenase or DNA beta-glucosyltransferase.
Embodiment 63. The method of Embodiment 58, wherein the double-stranded DNA substrate or the genomic DNA fragment is pre-treated with a TET methylcytosine dioxygenase, and optionally is pre-treated with a DNA beta-glucosyltransferase.
Embodiment 64. The method of Embodiment 59, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 80% identical to any of SEQ ID NOS: 40, 62, 63, 65, 67, 71, 110, 112, 114, and 117.
Embodiment 65. The method of Embodiment 61, wherein the double-stranded DNA deaminase has an amino acid sequence that is at least 80% identical to any of SEQ ID NOS: 47, 49, 50, 55, 58, 59, 70, 76, 106, 107, 163 and 164.
Embodiment 66. The method of Embodiment 57, wherein the double-stranded DNA substrate further comprises a genomic fragment linked to an adapter.
Embodiment 67. The method of Embodiment 66, wherein the adapter comprises a primer.
Embodiment 68. The method of Embodiment 57, wherein the strands of the double-stranded DNA substrate are not linked together by an adapter.
Embodiment 69. The method of Embodiment 57, wherein the deamination product is double-stranded.
Embodiment 70. The method of Embodiment 57, wherein the double-stranded DNA substrate is not a multi-copy strand.
Embodiment 71. The method of Embodiment 57, further comprising analyzing the sequence reads to identify a modified cytosine in the double-stranded DNA substrate.
Embodiment 72. The method of Embodiment 71, wherein a reference sequence is not used for the analyzing.
Embodiment 73. The method of Embodiment 71, wherein the modified cytosine is one or more of 5fC, 5CaC, 5mC, 5hmC, N4mC, or 5ghmC.
Embodiment 74. The method of Embodiment 71, wherein the modified cytosine is 5hmC. and the double-stranded DNA deaminase has an amino acid sequence that is at least 80% identical to any of SEQ ID NOS: 4, 5, 10, 13, 16, 96, 99 and 106.
Embodiment 75. A method for deaminating a nucleic acid, the method comprising: contacting:
Embodiment 76. The method of Embodiment 75, wherein the DNA substrate further comprises a modified cytosine.
Embodiment 77. The method of Embodiment 76, wherein the modified cytosine is a 5fC, 5CaC, 5mC, 5hmC, N4mC, 5ghmC, or pyrrolo-C.
Embodiment 78. An enzyme comprising an amino acid sequence that is at least 80% identical to any of SEQ ID NOS: 21, 40, 47, 49, 50, 55, 58, 59, 62, 63, 65, 67, 70, 71, 76, 106, 107, 110, 112, 114, 117, 163, and 164.
Embodiment 79. The enzyme of Embodiment 78, wherein the enzyme is fused with a DNA binding domain.
Embodiment 80. The enzyme of Embodiment 78, wherein the DNA binding domain is selected from a Cas9 domain, a Cas12 domain, a transcription activator-like effector nuclease (TALEN domain), a zinc finger (ZF) domain, a transcription activator-like effector (TALE) domain, an Sso7d domain, and a methyl binding domain (MBD) domain.
Embodiment 81. The enzyme of Embodiment 78, wherein the enzyme is no more than 300 amino acids in length.
Embodiment 82. A method for sequencing, comprising:
Embodiment 83. A kit comprising:
Embodiment 84. The kit of Embodiment 83, wherein the kit further comprises:
Embodiment 85. The kit of Embodiment 83, wherein the kit is free of TET methylcytosine dioxygenase and DNA beta-glucosyltransferase.
Embodiment 86. A reaction mix comprising:
Embodiment 87. The reaction mix of Embodiment 86, wherein the DNA substrate comprises cytosines and at least one modified cytosine.
Embodiment 88. The reaction mix of Embodiment 87, wherein the modified cytosine is a 5fC, 5caC, 5mC, 5hmC, N4mC or pyrrolo-C.
Embodiment 89. A method for base editing comprising:
Embodiment 90. The method of Embodiment 89, wherein the DNA binding domain is selected from a Cas9 domain, a Cas12 domain, a transcription activator-like effector nuclease (TALEN domain), a zinc finger (ZF) domain, a transcription activator-like effector (TALE) domain, and a methyl binding domain (MBD) domain.
Embodiment 91. The method of Embodiment 90, wherein the fusion protein further comprises a guide RNA complementary to at least a portion of the targeted sequence.
Embodiment 92. The method of Embodiment 89 wherein the fusion protein comprises an enzyme at is at least 80% identical to any of SEQ ID NOS:1-152.
Candidate DNA deaminase genes first were codon-optimized and then flanking sequences were added to each end, specifically, sequences containing T7 promoter at 5′ end and T7 terminator at 3′ end. These sequences were ordered as liner gBlocks from Integrated DNA Technologies (Coralville, IA, USA). Template DNA for in vitro protein synthesis was generated with Phusion® Hot Start Flex DNA Polymerase using gBlocks as template and flanking primers. The PCR products were purified using
Monarch PCR and DNA Cleanup kit (New England Biolabs, Inc., Ipswich, MA, USA). DNA concentration was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA). 100-400 ng PCR fragments were used as template DNA to synthesize analytic amounts of DNA deaminases using PURExpress In Vitro Protein Synthesis kit (New England Biolabs, Inc., Ipswich, MA, USA) following manufacturer's recommendations.
To test the activity of in vitro expressed DNA deaminases, a 2 ul aliquot of PURExpress sample was mixed with 300 ng of ΦX174 Virion DNA (ssDNA substrate) or ΦX174 RF I DNA (dsDNA substrate) in buffer containing 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100 and incubated for 1 h at 37° C. The deaminated ΦX174 DNA was purified using Monarch PCR and DNA Cleanup kit (New England Biolabs, Inc., Ipswich, MA, USA). DNA concentration was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA). 150 ng of deaminated DNAs were digested to nucleosides with the Nucleoside Digestion Mix (New England Biolabs, Inc., Ipswich, MA, USA) following manufacturer's recommendations. LC-MS/MS analysis was performed by injecting digested DNAs on an Agilent 1290 Infinity II UHPLC equipped with a G7117A diode array detector and a 6495C triple quadrupole mass detector operating in the positive electrospray ionization mode (+ESI). UHPLC was carried out on a Waters XSelect HSS T3 XP column (2.1×100 mm, 2.5 um) with a gradient mobile phase consisting of methanol and 10 mM aqueous ammonium acetate (pH 4.5). MS data acquisition was performed in the dynamic multiple reaction monitoring (DMRM) mode. Each nucleoside was identified in the extracted chromatogram associated with its specific MS/MS transition: dC [M+H]+at m/z 228.1112.1; dU [M+H]+ at m/z 229.1113.1; dmC [m+Fi]+ at m/z 242.1126.1; and dT [m+Fi]+ at m/z 243.1127.1. External calibration curves with known amounts of the nucleosides were used to calculate their ratios within the samples analyzed.
50 ng of E. coli C2566 genomic DNA was combined with control modified DNA's:
E. coli C2566
Then the DNA was transferred to a Covaris microTUBE (Covaris, Woburn, MA, USA) and sheared to 300 bp using the Covaris S2 instrument. The 50 μL of sheared material was transferred to a PCR strip tube to begin library construction. NEBNext DNA Ultra II Reagents (New England Biolabs, Ipswich, MA,
USA) were used according to the manufacturer's instructions for end repair, A-tailing, and adaptor ligation using an Illumina-compatible adapter. The ligated samples were mixed with 110 μL of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library was eluted in 17 μL of water.
The DNA was then deaminated in 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100, using 1μI of dsDNA deaminase synthesized as described above with an incubation time of 1 hour at 37° C. After deamination reaction, 1μI of Thermolabile Proteinase K (New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. 5p.M of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (New England Biolabs, Ipswich, MA, USA) were added to the DNA and PCR amplified. The PCR reaction samples were mixed with 50 ul of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library was eluted in 15 ul of water. The libraries were analyzed and quantified by High sensitivity DNA analysis using a chip inserted into an Agilent Bioanalyzer 2100. The whole-genome libraries were sequenced using the Illumina NextSeq platform. Pair-end sequencing of 150 cycles (2×75 bp) was performed for all the sequencing runs. Base calling and demultiplexing were carried out with the standard Illumina pipeline. Results of CseDa01 are shown in
50 ng of NA12878 genomic DNA was combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA and made up to 50 ul with 5 mM Tris pH=8.0. DNA was prepared according to Example 3 and the library was eluted in 29 ul of water. DNA was oxidized in a 50 ul reaction volume containing 50 mM Tris HCl pH 8.0, 1 mM DTT, 5 mM Sodium-L-Ascorbate, 20 mM a-KG, 2 mM ATP, 50 mM Ammonium Iron (II) sulfate hexahydrate, 0.04 mM UDG-glucose (NEB, Ipswich, MA), 16 μg mTET2, 10 U T4-BGT (NEB, Ipswich, MA). The reaction was initiated by adding Fe (II) solution to a final reaction concentration of 40 μM and then incubated for 1h at 37° C. The DNA was then deaminated, using 1 ul of MGYPDa829 dsDNA deaminase with an incubation time of 3 hour at 37° C. After deamination reaction, 1 ul of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. and 15 min at 60° C. At the end of the incubation, DNA was purified using 70 ul of resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. The sample was eluted in 16 ul water and 15 ul was transferred to a new tube. 1 μM of NEBNext Unique Dual Index Primers and 25 ul NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) were added to the DNA and PCR amplified. The libraries were analyzed and quantified with an Agilent Bioanalyzer 2100 DNA analyzer. The whole-genome libraries were sequenced, and analyzed as described below.
Raw reads were first trimmed by the Trim Galore software to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads due to adapter/quality trimming were also removed during this process. The trimmed read sequences were C to T converted and were then mapped to a composite reference sequence including the human genome (GRCh38) and the complete sequences of lambda and pUC19 controls using the Bismark program with default Bowtie2 setting (Langmead and Salzberg 2012). The aligned reads were then subjected to two post-processing QC steps: 1, alignment pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded; 2, reads that aligned to the human genome and contained excessive cytosines in non-CpG context (e.g., more than 3 in 75 bp) were removed because they are likely resulted from conversion errors. The numbers of T's (converted not methylated) and C's (unconverted modified) of each covered cytosine position were then calculated from the remaining good quality alignments using Bismark methylation extractor, and the methylation level was calculated as # of C/(# of C+# of T).
1500 ng of oligonucleotides
(ACACCCATCACATTTACAC(5caC)GGGAAAGAGTTGAATGTAGAGTTGG; SEQ ID NO: 157) or ACACCCATCACATTTACAC(5fC)GGGAAAGAGTTGAATGTAGAGTTGG; SEQ ID NO:158 with one modified cytosine (5caC or 5fC) were treated with CseDa01 DNA deaminase for 4 h in buffer containing 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100 and incubated for 1 h at 37° C. The deaminated oligonucleotides were purified using Monarch PCR and DNA Cleanup kit (New England Biolabs, Inc., Ipswich, MA, USA). DNA concentration was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA). 1500 ng of deaminated DNAs were digested to nucleosides with the Nucleoside
Digestion Mix (New England Biolabs, Inc., Ipswich, MA, USA) following manufacturer's recommendations. UHPLC-MS analysis was performed using an Agilent 1290 Infinity II UHPLC equipped with G7117A Diode Array Detector and 6135 XT MS Detector, on a Waters XSelect HSS T3 XP column (2.1×100 mm, 2.5 μm) with the gradient mobile phase consisting of methanol and 10 mM ammonium acetate buffer (pH 4.5). The identity of each peak was confirmed by MS. The relative abundance of each nucleoside was determined by the integration of each peak at 260 nm or their respective UV absorption maxima. Results are shown in
50 ng of NA12878 genomic DNA was combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA and made up to 50 μL with 5 mM Tris pH=8.0. DNA was prepared according to Example 3 and the library was eluted in 29 μL of water. DNA was oxidized in a 50 μL reaction volume containing 50 mM Tris HCl pH 8.0, 1 mM DTT, 5 mM Sodium-L-Ascorbate, 20 mM a-KG, 2 mM ATP, 50 mM Ammonium Iron (II) sulfate hexahydrate, and 16 μg mTET2. The reaction was initiated by adding Fe (II) solution to a final reaction concentration of 40 μM and then incubated for 1 h at 37° C. The DNA was then deaminated, using 1 μl of CseDa01 dsDNA deaminase with an incubation time of 3 hour at 37° C. After deamination reaction, 1 μl of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. and 15 min at 60° C. At the end of the incubation, DNA was purified using 70 μL of resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. The sample was eluted in 16 μL water and 15 μL was transferred to a new tube. 1 μM of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) were added to the DNA and PCR amplified. The libraries were analyzed and quantified with an Agilent Bioanalyzer 2100 DNA analyzer. The whole-genome libraries were sequenced, and analyzed as described below. Raw reads were first trimmed by the Trim Galore software to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads due to adapter/quality trimming were also removed during this process. The trimmed read sequences were C to T converted and were then mapped to a composite reference sequence including the human genome
(GRCh38) and the complete sequences of lambda and pUC19 controls using the Bismark program with default Bowtie2 setting (Langmead and Salzberg 2012). The aligned reads were then subjected to two post-processing QC steps: 1, alignment pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded; 2, reads that aligned to the human genome and contained excessive cytosines in non-CpG context (e.g., more than 3 in 75 bp) were removed because they are likely resulted from conversion errors. The numbers of T's (converted not methylated) and C's (unconverted modified) of each covered cytosine position were then calculated from the remaining good quality alignments using Bismark methylation extractor, and the methylation level was calculated as # of C/(# of C+# of T).
To test the activity of CseDa01 DNA deaminase in TET2 buffer a 2μI of PURExpress sample was mixed with 300 ng of ΦX174 Virion DNA (ssDNA substrate) or ΦX174 RF I DNA (dsDNA substrate) in buffer containing 50 mM Tris HCl pH 8.0, 1 mM DTT, 5 mM Sodium-L-Ascorbate, 20 mM a-KG, 2 mM ATP, 50 mM Ammonium Iron (II) sulfate hexahydrate, 0.04 mM, and incubated for 1 h at 37° C. The deaminated ΦX174 DNA was purified using Monarch PCR and DNA Cleanup kit (New England Biolabs, Inc., Ipswich, MA, USA). DNA concentration was quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Inc., Waltham, MA, USA). 150 ng of deaminated DNAs were digested to nucleosides with the Nucleoside Digestion Mix (New England Biolabs, Inc., Ipswich, MA, USA) following manufacturer's recommendations. LC-MS/MS analysis was performed by injecting digested DNAs on an Agilent 1290 Infinity II UHPLC equipped with a G7117A diode array detector and a 6495C triple quadrupole mass detector operating in the positive electrospray ionization mode (+ESI). UHPLC was carried out on a Waters XSelect HSS T3 XP column (2.1×100 mm, 2.5 μm) with a gradient mobile phase consisting of methanol and 10 mM aqueous ammonium acetate (pH 4.5). MS data acquisition was performed in the dynamic multiple reaction monitoring (DMRM) mode. Each nucleoside was identified in the extracted chromatogram associated with its specific MS/MS transition: dC [m+H]+ at m/z 228.1→112.1; dU [M+H]+ at m/z 229.1→113.1; dmC [m+H]+ at m/z 242.1→126.1; and dT [m+H]+ at m/z 243.1→127.1. External calibration curves with known amounts of the nucleosides were used to calculate their ratios within the samples analyzed. Results are shown in
50 ng of E. coli C2566 genomic DNA was combined with 2 ng unmethylated lambda, phage XP12 (all cytosines are 5-methylcytosines) and T4 phage DNA (all cytosines are 5-hydroxymethyl cytosines) control DNAs and made up to 50 μL with 10 mM Tris, pH 8.0. Then the DNA was prepared according to Example 3 with a sheared size of 240-290 bp and a library elution volume of 15 μL of water. The DNA was then deaminated in 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100, using 1μI of a modification-sensitive dsDNA deaminase (e.g., MGYPDa20 or NsDa01) synthesized as described above with an incubation time of 1 hour at 37° C. After deamination reaction, 1μI of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. 1μM of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) were added to the DNA and PCR amplified. The PCR reaction samples were mixed with 50 μL of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library was eluted in 15 μL of water. The libraries were analyzed and quantified by High sensitivity DNA analysis using a chip inserted into an Agilent Bioanalyzer 2100. The whole-genome libraries were sequenced using the Illumina NextSeq platform. Pair-end sequencing of 150 cycles (2×75 bp) was performed for all the sequencing runs. Base calling and demultiplexing were carried out with the standard Illumina pipeline. Raw reads were first trimmed by the Trim Galore to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads owing to adapter/quality trimming were also removed during this process. The trimmed read sequences were C-to-T converted and were then mapped to a composite reference sequence including the E. coli C2566 genome and the complete sequences of lambda, phage XP12, and T4 controls using the Bismark program with the default Bowtie 2 setting.
The first 5 bp at the 5′ end of R2 reads were removed to reduce end-repair errors and aligned read pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded. Next deamination events (C->T) were called by comparing the remaining good alignment sequences to the reference sequences using Bismark methylation extractor program. The 20 bp flanking sequences (10 bp upstream and 10 bp downstream) of all the covered cytosines from the individual genomes were then extracted and the cytosines sites were divided into different groups based on their deamination rates (>=90%, >=50%, >=25% or <=10%). Flanking sequences of each cytosine group were used to make sequence logo using WebLogo 3 to infer deamination sequence preference. Results are shown in
50 ng of NA12878 genomic DNA was combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA and made up to 50 μL with 5 mM Tris pH=8.0. DNA was prepared according to Example 3 and the library was eluted in 17 μL of molecular grade water. The DNA was then deaminated in 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100, using 1μI of MGYPDa20 dsDNA deaminase with an incubation time of 3 hours at 37° C. Other modification sensitive deaminases may be substituted (e.g., see Table 3). After deamination reaction, 1μI of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. 5 μM of NEBNext Unique Dual Index Primers, 20 μM deaminated DNA and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) were combined and PCR amplified. The PCR reaction samples were mixed with 50 μl of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library was eluted in 15 μL of water. The libraries were analyzed and quantified by High sensitivity DNA analysis using a chip inserted into an Agilent Bioanalyzer 2100. The whole-genome libraries were sequenced using the Illumina NextSeq platform and analyzed as described below. Raw reads were first trimmed by the Trim Galore software to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads due to adapter/quality trimming were also removed during this process. The trimmed read sequences were C to T converted and were then mapped to a composite reference sequence including the human genome (GRCh38) and the complete sequences of lambda and pUC19 controls using the Bismark program with default Bowtie2 setting (Langmead and Salzberg 2012). The aligned reads were then subjected to two post-processing QC steps: 1, alignment pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded; 2, reads that aligned to the human genome and contained excessive cytosines in non-CpG context (e.g., more than 3 in 75 bp) were removed because they are likely resulted from conversion errors. The numbers of T's (converted not methylated) and C's (unconverted modified) of each covered cytosine position were then calculated from the remaining good quality alignments using Bismark methylation extractor, and the methylation level was calculated as # of C/(# of C+# of T).
For whole human genome methyl-SNP-seq sequencing 4μg of NA12878 gDNA and 40 ng of unmethylated lambda DNA as spiked in to monitor the deamination efficiency were used. The genomic DNA was fragmented using 250 bp sonication protocol using a Covaris S2 sonicator. Two technical replicates were set up. The fragmented gDNA was end repaired and dA-tailed (NEB Ultra II E7546 module), then ligated to the custom hairpin adapter using NEB ligase master mix (NEB, M0367). The incomplete ligation product (fragment having only one or no adaptor ligated) was removed using two exonucleases (NEB exolll and NEB exoVIl). Two nick sites were created at the uracil positions in the hairpin adapters at both ends after being treated with UDG and EndoVIll. The nick sites were translated towards 3′ terminus by DNA polymerase I in the presence of dATP, dGTP, dTGP and 5-methyl-dCTP. The nick translation causes double stranded DNA break when DNA polymerase I encounters the other nick on the opposite strand. The resulting fragments have one end ligated to a hairpin adapter and blunt end on the other side. The blunt end was dA-tailed and ligated with methylated Illumina adapter. The ligated product was deaminated at 37° C. for 3 h with double stranded DNA deaminase MGYPDa20. The deaminated DNA product was amplified using NEBNext Q5U Master Mix (NEB, M0597). The resulting indexed library was used for Illumina sequencing. The human Methyl-SNP-seq libraries were sequenced using an Illumina Novaseq 6000 sequencer for 100 bp paired end reads.
50 ng of Paenibacillus species JDR-2 (CCGG target sequence) and Salmonella enterica FDAARGOS_312 (CACCGT target sequence) DNAs were combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA and made up to 50 μL with 5 mM Tris pH=8.0. DNA was prepared according to Example 3 with a sheared size of240-290 bp and an elution volume of 15 μL of water. The DNA was then deaminated in 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100, using 1μI of CseDa01 dsDNA deaminase synthesized as described above with an incubation time of 1 hour at 37° C. After deamination reaction, 1 μl of Thermolabile Proteinase K (P8111S, New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. 1 μM of NEBNext Unique Dual Index Primers and 25 μl NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) were added to the DNA and PCR amplified. The PCR reaction samples were mixed with 50 ul of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library was eluted in 15 μl of water. The libraries were analyzed and quantified by High sensitivity DNA analysis using a chip inserted into an Agilent Bioanalyzer 2100. The whole-genome libraries were sequenced using the Illumina NextSeq platform. Pair-end sequencing of 150 cycles (2×75 bp) was performed for all the sequencing runs. Raw reads were first trimmed by the Trim Galore to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads owing to adapter/quality trimming were also removed during this process. The trimmed read sequences were C-to-T converted and were then mapped to the reference sequence and the complete sequences of lambda and pUC19 controls using the Bismark program with the default Bowtie 2 setting. The first 5 bp at the 5′ end of R2 reads were removed to reduce end-repair errors and aligned read pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded. Next deamination events (C->T) were called by comparing the remaining good alignment sequences to the reference sequences using Bismark methylation extractor program. An N4mC modified site is called when it is largely un-deaminated (C->T conversion rate <=20%). The flanking 20 bp sequences of all the called N4mC sites were extracted and a sequence logo using WebLogo 3 was generated. Results are shown in
50 ng of NEB1569 Thermus species M and NEB 394 Acinetobacter species H genomic DNAs was combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA and made up to 50 μl with 5 mM Tris pH=8.0. Then the DNA was prepared according to Example 3 with a sheared size of 240-290 bp and a library elution volume of 15 μl of water. The DNA was then deaminated in 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100, using 1 μl of dsDNA deaminase synthesized as described above with an incubation time of 1 hour at 37° C. After deamination reaction, 1 μl of Thermolabile Proteinase K (P8111S, New England Biolabs, Ipswich, MA) was added and incubated additional 30 min at 37° C. 1 μM of NEBNext Unique Dual Index Primers and 25 μl NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) were added to the DNA and PCR amplified. The PCR reaction samples were mixed with 50 μl of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library was eluted in 15 μl of water. The libraries were analyzed and quantified by High sensitivity DNA analysis using a chip inserted into an Agilent Bioanalyzer 2100. The whole-genome libraries were sequenced using the Illumina NextSeq platform. Pair-end sequencing of 150 cycles (2×75 bp) was performed for all the sequencing runs. Base calling and demultiplexing were carried out with the standard Illumina pipeline. Raw reads were first trimmed by the Trim Galore to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads owing to adapter/quality trimming were also removed during this process. The trimmed read sequences were C-to-T converted and were then mapped to a composite reference sequence including the NEB1569 Thermus species M and NEB 394 Acinetobacter species H and the complete sequences of lambda and pUC19 controls using the Bismark program with the default Bowtie 2 setting. The first 5 bp at the 5′ end of R2 reads were removed to reduce end-repair errors and aligned read pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded. Next deamination events (C->T) were called by comparing the remaining good alignment sequences to the reference sequences using Bismark methylation extractor program. The N4mC modification is called from the CseDa01 deaminase-treated library. An N4mC modified site is called when it is largely un-deaminated (C->T conversion rate <=20%). For 5mC modification detection, a differential methylation analysis was conducted between the MGYPDa20 deaminase-treated library (detect both N4mC and 5mC) and the CseDa01 deaminase-treated library (detect only N4mC) of the same sample to identify modified sites (i.e., 5mC) that are only detected in the MGYPDa20 library. The differentially methylated sites were called by a logistic regression method with SLIM corrected Q value <=0.01, and methylation difference >=80% using the Methylkit program. To identify methyltransferase recognition sequences, the 9 bp flanking sequences were extracted, including 4 bp upstream and 4 bp downstream of all the modified sites, and the unique 9 bp sequences were clustered using a hierarchical linkage method based on the difference between each pair of sequences. A sequence logo was generated using WebLogo 3 for each cluster representing a distinct methyltransferase recognition motif.
A list of HMMER3 (Eddy, S. R. Accelerated Profile HMM Searches. PLOS Comput. Biol. 7, e1002195 (2011)) cytosine deaminase sequence profiles was curated. 29 profiles came from the CDA clan (CL0109) from the Pfam (Mistry, J. et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 49, D412—D419 (2021)) database (excluding the TM1506, Lpxl_C, FdhD-NarQ, and AICARFT_IMPCHas, which do not encode deaminases), 17 profiles were built from multiple sequence alignments (MSAs) of deaminase families defined by lyer et al. (Nucleic Acids Res. 39, 9473-9497, 2011), and one profile was built from a multiple sequence alignment found in Zhang et al. (Biol. Direct 7, 18, 2012).
Some candidate sequences were selected directly from the MSAs listed in lyer et al. (2011), and Zhang et al. (2012). Others were selected from hmmsearch hits of the profiles described above against six different databases: UniProt, Mgnify, IMG/VR, IMG/M, wastewater treatment plant metagenomes, and GenBank (respectively, The UniProt Consortium. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 49, D480—D489 (2021); Mitchell, A. L. et al. MGnify: the microbiome analysis resource in 2020. Nucleic Acids Res. 48, D570—D578 (2020); Paez-Espino, D. et al. IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses. Nucleic Acids Res. 45, gkw1030 (2017); Chen, I.-M. A. et al. The IMG/M data management and analysis system v.6.0: new tools and advanced capabilities. Nucleic Acids Res. 49, D751—D763 (2021); Singleton, C. M. et al. Connecting structure to function with the recovery of over 1000 high-quality metagenome-assembled genomes from activated sludge using long-read sequencing. Nat. Commun. 12, 2009 (2021); and Da, B. et al. GenBank. Nucleic Acids Res. 41, (2013)).
Most of the deaminases tested were found as fusions to larger proteins, for example as parts of polymorphic toxin systems. To determine the boundaries of the deaminase domain, AlphaFold2 (Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 1-11 (2021) doi:10.1038/s41586-021-03819-2) structural predictions were generated and visualized. N-terminal truncation sites were generally selected at several amino acids before helix 1 of the deaminase domain.
For convenience, each screened sequence was given a short name. The names are arbitrary, but relate somehow to the database or species of origin for the sequence. Da=deaminase, MGYP=Mgnify protein, Hm=hot metagenome, VR=IMG/VR, WWTP=waste water treatment plant, chimera=chimeric sequence, Anc=ancestral sequence reconstruction. Other prefixes are mostly two or three letters drawn from the name of the source organism or the source environment of the metagenome data. Some sequences also have prefixes or suffixes of the form extN #, extC #, d #, Cd #, which indicate, respectively, N-terminal extensions, C-terminal extensions, N-terminal deletions, and C-terminal deletions of the indicated number of residues, compared to the candidate with the un-affixed name.
Amino acid sequence alignments were all calculated using MAFFT (v7.490) (Katoh, K. & Standley, D. M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 30, 772-780 (2013)) using globalpair mode. Trees were generated using raxml-ng (v. 1.1)(Kozlov, A. M., Darriba, D., Flouri, T., Morel, B. & Stamatakis, A. RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 35, 4453-4455 (2019)). Ancestral sequence reconstructions were built from phylogenetic trees using raxml-ng (v. 1.1)
Assay results for 29 deaminases are shown in Table 1 below, in which APOBEC3A (a single-stranded DNA deaminase) served as a negative control. The other 28 deaminases (double-stranded DNA deaminases) in the table all have significant activity on a double-stranded DNA substrate.
Double-stranded DNA deaminases disclosed herein may be used in many methods, processes, and workflows including, for example, the applications shown in Table 2 below. Deamination products may contain one or more modified cytosines, for example, where the substrate dsDNA included such modified cytosines and the operative deaminase does not or only poorly deaminases such modified cytosines. Each of the listed methods/applications may further comprise (a)(i) sequencing the deamination products and/or (ii) amplifying (e.g., by PCR) the deamination products to produce amplification products and sequencing the amplification products, in each of (a)(i) and (a)(ii), to produce sequence reads, and (b) optionally determining the kind and/or position of modified cytosines in the dsDNA substrate from the sequence reads.
Screening results for over 100 deaminases are shown in Table 3 below, in which APOBEC3A (a single-stranded DNA deaminase) served as a negative control. Many were observed to have double-stranded DNA deaminase activity under the conditions tested. Relatedness of the enzymes tested is illustrated in
The names and SEQ ID NOS of certain double-stranded DNA deaminases disclosed herein are shown in Table 4 along with the corresponding names included in U.S. Provisional Application No. 63/264,513 filed Nov. 24, 2021.
Combine 50 ng of NA12878 genomic DNA with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNAs and make up to 50 μl with 5 mM Tris pH=8.0. Shear the DNA to about 300 bp using any method. For example, transfer DNA to a Covaris microtube (Covaris, Woburn, MA) and shear according to the manufacturer's protocol.
Add 50 ng of sheared DNA to a PCR strip tube to begin library construction. Use NEBNext DNA Ultra II Reagents (NEB, Ipswich, MA) according to the manufacturer's instructions for end repair, A-tailing, and adaptor ligation of the custom made Pyrollo-dC adaptor, where all dC's are replaced with Pyrollo-dC: ACACTCTTTCCCTACACGACGCTCTTCCGATCT (SEQ ID NO:165) and [Phos]GATCGGAAGAGCACACGTCTGAACTCCAGTCA (SEQ ID NO:166).. EM-seq adaptor (E7120S/L, NEB, Ipswich, MA) or any other desired adapter with or without replacement of dCs with modified dCs may also be used. Purify the adapter ligated DNA using 1× NEBNext Sample Purification Beads according to the manufacturer's instructions.
Combine the adapter ligated DNA with 5X TET2 reaction buffer and the mTET2 and T4-BGT enzymes in a 50 μl reaction volume. Initiate 5mC oxidation reaction by adding Fe (II) solution (Ammonium Iron (II) sulfate hexahydrate) to a final reaction concentration of 50 μM and then incubate for 1 h at 37° C. Add 1 μl of Thermolabile Proteinase K (P8111S, New England Biolabs, Ipswich, MA) and incubate additional 30 min at 37° C. and 10 min at 60° C.
Purify DNA using 1X NEBNext Sample Purification Beads according to the manufacturer's protocol and elute in 17 μl water. Add 2 μl of 10X deaminase buffer, 1 μl of CbDa01 and incubate for 3 h at 37° C. Incubation time may be shortened or extended depending on factors such as temperature, enzyme concentration, etc. The deamination reaction may be stopped in a variety of ways (e.g., enzymatically, separation step, etc.). For example, add 1 μl of Thermolabile Proteinase K (P8111S, New England Biolabs, Ipswich, MA) to each tube and incubate for an additional 30 min at 37° C. and 10 min at 60° C. Mix with 5 μl of H20, 15 μl of deamination reaction, 5 μl of NEBNext Unique Dual Index Primers and 25 μl NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) and amplify using EM-seq protocol (6 PCR cycles). After PCR reaction, purify amplified DNA using 1X resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. Use a Bioanalyzer or TapeStation to determine the size distribution and concentration of the libraries. Data analysis may be conducted as described in Example 22.
Combine 50 ng of NA12878 genomic DNA with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNAs and make up to 50 μl with 5 mM Tris pH=8.0. Shear the DNA to about 300 bp using any method. For example, transfer DNA to a Covaris microtube (Covaris, Woburn, MA) and shear according to the manufacturer's protocol. Add 50 ng of sheared DNA to a PCR strip tube to begin library construction. Use NEBNext DNA Ultra II Reagents (NEB, Ipswich, MA) according to the manufacturer's instructions for end repair, A-tailing, and adaptor ligation of the custom made Pyrollo-dC adaptor, where all dC's are replaced with Pyrollo-dC: ACACTCTTTCCCTACACGACGCTCTTCCGATCT (SEQ ID NO:165) and [Phos]GATCGGAAGAGCACACGTCTGAACTCCAGTCA (SEQ ID NO:166). EM-seq adaptor (E7120S/L, NEB, Ipswich, MA) or any other desired adapter with or without replacement of dCs with modified dCs may also be used. Purify the adapter ligated DNA using 1× NEBNext Sample Purification Beads according to the manufacturer's instructions and elute in 17 μl of water.
Add 2 μl of 10X deaminase buffer and 1 μl of the RhDa01 deaminase and incubate for 3 h at 37° C. Incubation time may be shortened or extended depending on factors such as temperature, enzyme concentration, etc. After deamination reaction, add 1μI of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) to each tube and incubated additional 30 min at 37° C. and 10 min at 60° C. Note that the deamination reaction may be stopped in a variety of alternative ways (e.g., enzymatically, separation step, etc.). Mix with 5μI of H2O, 15 μL of deamination reaction, 5μI of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) and amplify using EM-seq protocol (6 PCR cycles). After PCR reaction, purify amplified DNA using 1X resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. Use a Bioanalyzer or TapeStation to determine the size distribution and concentration of the libraries. Data analysis may be conducted as described in Example 22.
Combine 5˜10 pg of NA12878 genomic DNA or DNA from a single cell with 0.2% of CpG methylated pUC19 and 2% of unmethylated lambda control DNAs and make up to 50 μL with 5 mM Tris pH=8.0. Shear the DNA to about 300 bp using any method. For example, transfer DNA to a Covaris microtube (Covaris, Woburn, MA) and shear according to the manufacturer's protocol. Transfer the sheared DNA to a PCR strip tube to begin library construction. Use NEBNext DNA Ultra II Reagents (NEB, Ipswich, MA) according to the manufacturer's instructions but reduce the reaction volumes to half for end repair, A-tailing, and adaptor ligation of the custom made Pyrollo-dC: ACACTCTTTCCCTACACGACGCTCTTCCGATCT (SEQ ID NO:165) and [Phos]GATCGGAAGAGCACACGTCTGAACTCCAGTCA (SEQ ID NO:166). EM-seq adaptor (E7120S/L, NEB, Ipswich, MA) or any other desired adapter with or without replacement of dCs with modified dCs may also be used. Purify the adapter ligated DNA using 1× NEBNext Sample Purification Beads according to the manufacturer's instructions and elute in 16 μl of water.
Add 2 μL of 10X deaminase buffer, 1 μL of 10 ng carrier DNA (any not adapter ligated DNA could be used as carrier DNA), and 1 μL the RhDa01 enzyme. Mix and incubate for 3 h at 37° C. Incubation time could be shortened or extended depending on enzyme concentration. For deaminases acting on ssDNA, the DNA substrate could be denatured using heat or any chemical denaturing agent. The deamination reaction may be stopped in a variety of ways (e.g., enzymatically, separation step, etc.). For example, add 1 μL of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) to each tube and incubated additional 30 min at 37° C. and 10 min at 60° C.). Mix with 5 μL of H20, 15 pi of deamination reaction, 5 μL of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) and amplify using EM-seq protocol (1215 PCR cycles). After PCR reaction, purify amplified DNA twice using 1X resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. Use a Bioanalyzer or TapeStation to determine the size distribution and concentration of the libraries.
The library was sequenced on an Illumina NovaSeq machine and resulted in 154 million 2x100 bp paired-end reads. Data analysis of 5mC and genetic base detection was conducted using the method described in the Example 22. About 95% of the reads were mapped to the human reference genome using a standard DNA sequencing aligner Bowtie2. The final average sequencing coverage of the human genome after removing PCR duplicates and muti-mapped reads is about 6×. The base-resolution methylation results of this method agree well with of the published EM-seq method, with a high Pearson correlation value of 0.87 for the CpG sites that have a minimum 5× coverage. The CpG island methylation levels are also highly correlated between the two methods (Pearson correlation=0.95). This method also produced SNP calling results even at a low sequencing depth of 6× average coverage. Using chromosome 21 as an example and benchmarking against the variants obtained using the NA12878 whole-genome sequencing data set (WGS, performed by the JIMB NIST project), more than 90% of the SNPs detected by this method are identified in the reference dataset.
Combine 50 ng of NA12878 genomic DNA with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNAs and make up to 50 μL with 5 mM Tris pH=8.0. Shear the DNA to about 300 bp using any method. For example, transfer DNA to a Covaris microtube (Covaris, Woburn, MA) and shear according to the manufacturer's protocol. Add 50 ng of sheared DNA to a PCR strip tube to begin library construction. Use NEBNext DNA Ultra II Reagents (NEB, Ipswich, MA) according to the manufacturer's instructions for end repair, A-tailing, and adaptor ligation of the custom made Pyrollo-dC adaptor, where all dC's are replaced with Pyrollo-dC: ACACTCTTTCCCTACACGACGCTCTTCCGATCT (SEQ ID NO:165) and [Phos]GATCGGAAGAGCACACGTCTGAACTCCAGTCA (SEQ ID NO:166). EM-seq adaptor (E7120S/L, NEB, Ipswich, MA) or any other desired adapter with or without replacement of dCs with modified dCs may also be used. Purify the adapter ligated DNA using 1× NEBNext Sample Purification Beads according to the manufacturer's instructions.
Combine the adapter ligated DNA with T4-BGT enzyme in T4-BGT buffer (NEB, Ipswich, MA) in a 50 μL reaction volume, and incubate for 1 h at 37° C. Reaction time may be adjusted according to substrate quantity. After deamination, stop the reaction (e.g., add 1 μL of Thermolabile Proteinase K (P8111S, New England Biolabs, Ipswich, MA) and incubate for an additional 30 min at 37° C. and 10 min at 60° C.).
Purify DNA using 1× NEBNext Sample Purification Beads according to the manufacturer's protocol and elute in 17 μL water. Add 2 μL of 10X deaminase buffer, 1 μL of CbDa01 (SEQ ID NO: 50), AcDa01 (SEQ ID NO: 49), Sso7d_GGGVTS_AcDa01 (SEQ ID NO: 163) or Sso7d_LSGLSDDKLKEI_AcDa01 (SEQ ID NO: 164), and incubate for a sufficient time for deamination, e.g., 3 h at 37° C. Incubation time may be shortened or extended depending on factors such as temperature, enzyme concentration, etc. For deaminases acting on ssDNA, the DNA substrate may be denatured using heat, enzymatic methods, or chemical methods. The deamination reaction may be stopped in a variety of ways (e.g., enzymatically, separation step, etc.). For example, add 1 μL of Thermolabile Proteinase K (P8111S, New England Biolabs, Ipswich, MA) to each tube and incubate for an additional 30 min at 37° C. and 10 min at 60° C. Mix with 5 μL ofH20, 15 μL of deamination reaction, 5 μL of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) and amplify using EM-seq protocol (6 PCR cycles). After PCR reaction, purify amplified DNA using 1X resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. Use a Bioanalyzer or TapeStation to determine the size distribution and concentration of the libraries. Data analysis may be conducted as described in Example 22.
50 ng of NA12878 genomic DNA is combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA. Deamination is performed in a suitable buffer (e.g., 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100), using 1 μL of dsDNA deaminase having bias for CpG (see, e.g., Table 3) with an incubation time, for example, of 3 hours at 37° C. Enzyme amount, temperature, and incubation time could be adjusted depending on deaminase activity. As is described herein, dsDNA deaminases are active on both ssDNA and dsDNA. For deaminases acting on ssDNA, the DNA substrate may be denatured using heat, enzymatic methods, or chemical methods. The deamination reaction may be stopped in a variety of ways (e.g., enzymatically, separation step, etc.). For example, add 1 μL of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) and incubate for an additional 30 min at 37° C. The DNA samples are mixed with purification beads and cleaned up according to the manufacturer's instructions. DNA is fragmented to about 300 bp. The 50 μL of sheared material is transferred to a PCR strip tube to begin library construction. NEBNext DNA Ultra II Reagents (NEB, Ipswich, MA) are used according to the manufacturer's instructions for end repair, A-tailing, and adaptor ligation. Any adapter sequence could be used. The ligated samples are mixed with 110 pi of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. 5 μL of NEBNext Unique Dual Index Primers, 20 μL of deaminated DNA and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) are combined and PCR amplified. The PCR reaction samples are mixed with 50 μL of resuspended NEBNext Sample Purification Beads and cleaned up according to the manufacturer's instructions. The library is eluted in 15 μL of water. The libraries are analyzed and quantified by High sensitivity DNA analysis using a chip inserted into an Agilent Bioanalyzer 2100. The whole-genome libraries are sequenced, e.g., using the Illumina NextSeq platform. Data analysis of 5mC and genetic base detection may be conducted as described in Example 22.
For 5mC detection, 200 ng of human genomic DNA is oxidized by incubating with 16 μg of TET2 for 30 min at 37° C. followed 30-min incubation with BGT in the same buffer at 37° C. For 5hmC detection, genomic DNA is glucosylated by incubating with BGT enzyme for 2 h at 37° C. Modification protected genomic DNA is incubated for an additional 30 min with Proteinase K at 37° C. and subsequently purified using a genomic DNA purification kit. Purified DNA is deaminated with 2 μL of CbDa01 CpG dsDNA deaminase in 100 μL reaction volume for 3 hours. Incubation time may be shortened or extended depending on factors such as temperature, enzyme concentration, etc. For deaminases acting on ssDNA, the DNA substrate may be denatured using heat, enzymatic methods, or chemical methods. The deamination reaction may be stopped in a variety of ways (e.g., enzymatically, separation step, etc.). For example, add 1 μL of Thermolabile Proteinase K and incubate for an additional 30 min at 37° C. and 10 min at 60° C. Targeted genomic regions are amplified from the purified deaminated DNA using custom designed primers. After purification of PCR product, the long amplicons are used to prepare a PacBio SMRT sequencing (Pacific Biosciences) library following the “amplicon template preparation and sequencing” protocol, and the library is sequenced on a PacBio machine following manufacturer's instruction. Circular Consensus Sequences (CCS) are extracted from the raw data and converted into FASTQ file using SMRT Link program. 5mC and sequence analysis may be conducted as described in Example 22.
For 5mC detection, 200 ng of human genomic DNA is oxidized by incubating with 16 μg of TET2 for 30 min at 37° C. followed 30-min incubation with BGT in the same buffer at 37° C. For 5hmC detection, genomic DNA is glucosylated by incubating with BGT enzyme for 2 h at 37° C. Modification protected genomic DNA was incubated for an additional 30 min with Proteinase K at 37° C. and subsequently purified using a genomic DNA purification kit. Purified DNA is deaminated with 2 μL of CbDa01 CpG deaminase in 100 μL reaction volume for 3 hours. Incubation time may be shortened or extended depending on factors such as temperature, enzyme concentration, etc. For deaminases acting on ssDNA, the DNA substrate could be denatured using heat or any chemical denaturing agent. The deamination reaction may be stopped in a variety of ways (e.g., enzymatically, separation step, etc.). For example, add 1 μL of Thermolabile Proteinase K and incubate additional 30 min at 37° C. and 10 min at 60° C. Targeted genomic regions are amplified from the purified deaminated DNA using custom designed primers. After purification of PCR product, the long amplicons are used to prepare a Nanopore sequencing library using the 1D Native barcoding genomic DNA kit, and the library is sequenced on a MinION flow cell. Raw reads are base called using the Guppy base caller. Analysis of 5mC and sequence may be conducted as described in Example 22.
Raw sequencing reads are quality trimmed to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads due to adapter/quality trimming are also removed during this process. The trimmed read sequences are then mapped to a composite reference sequence including the human genome (GRCh38) and the complete sequences of lambda and pUC19 controls using a standard sequence alignment tool e.g., Bowtie2. Alignment pairs that shared the same alignment start positions (5′ ends) are regarded as PCR duplicates and are discarded. The alignments are used for SNP and other genetic variants detection (except for variant detection of Cs in CpG context) using a standard variant calling analysis pipeline, such as GATK.
For variant and methylation analysis of Cs in CpG context, C to T conversion events in CpG context are detected and summarized in a strand-specific manner using the samtools by the mpileup function. The G to A conversion events in the G positions that pair with the Cs of CpGs on the opposite strand are also detected and summarized using the samtools mpileup program. For each cytosine position in a CpG context, the C->T conversion rate and paired opposite strand G->A conversion rate is compared. If the two conversion rates are not statistically different, then all the C->T conversions are considered as a results of genetic variant, and no deamination on this cytosine position. If the C->T conversion rate is significantly larger than the G->A conversion rate, then the difference is considered as a result of deamination of unmodified cytosines. The methylation level of a cytosine in CpG context is calculated as #of methylated unconverted Cs/(#of methylated unconverted Cs+#of unmodified and converted Cs), whereas #of unmodified and converted Cs=number of C->T conversion events—number of G->A conversion events on the opposite strand.
This example describes detecting 5hmC using a dsDNA deaminase having activity on unmodified and methylated dsDNA and is blocked by 5ghmC (e.g., MGYPDa829, LbsDa01, MGYPDa01, PeDa01, MGYPDa06, or any of these deaminases fused to Sso7d, or a combination). 50 ng of NA12878 genomic DNA is combined with 0.1 ng of CpG methylated pUC19 and 1 ng of unmethylated lambda control DNA and made up to 50 μL with 5 mM Tris pH=8.0. DNA according to Example 3 and the library is eluted in 29 μL of water. The adapter ligated DNA is combined with T4-BGT enzyme in T4-BGT buffer (NEB, Ipswich, MA) in a 50 μL reaction volume, and incubated for 1 h at 37° C. Reaction time may be adjusted, for example, according to substrate quantity. 1 μL of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) is added and incubated additional 30 min at 37° C. and 10 min at 60° C. The DNA is purified using 1× NEBNext Sample Purification Beads according to the manufacturer's protocol and eluted in 17 μL water. The DNA is then deaminated, using, e.g., 1 μL of MGYPDa829 dsDNA deaminase and 2 μL 5× deamination buffer, incubation time of 3 hour at 37° C. After deamination reaction, the reaction is stopped, e.g., using 1 μL of Thermolabile Proteinase K (P81115, New England Biolabs, Ipswich, MA) and incubation for an additional 30 min at 37° C. and 15 min at 60° C. At the end of the incubation, DNA is purified using 70 μL of resuspended NEBNext Sample Purification Beads according to the manufacturer's protocol. The sample is eluted in 16 μL water and 15 μL is transferred to a new tube. 1 μM of NEBNext Unique Dual Index Primers and 25 μL NEBNext Q5U Master Mix (M0597, New England Biolabs, Ipswich, MA) are added to the DNA and PCR amplification is performed. The libraries are analyzed and quantified with an Agilent Bioanalyzer 2100 DNA analyzer.
Raw reads are first trimmed by the Trim Galore software to remove adapter sequences and low-quality bases from the 3′ end. Unpaired reads due to adapter/quality trimming are also removed during this process. The trimmed read sequences were C to T converted and were then mapped to a composite reference sequence including the human genome (GRCh38) and the complete sequences of lambda and pUC19 controls using the Bismark program with default Bowtie2 setting (Langmead and Salzberg 2012). The aligned reads were then subjected to two post-processing QC steps: 1, alignment pairs that shared the same alignment start positions (5′ ends) were regarded as PCR duplicates and were discarded; 2, reads that aligned to the human genome and contained excessive cytosines in non-CpG context (e.g., more than 3 in 75 bp) were removed because they are likely resulted from conversion errors. The numbers of T's (converted not methylated) and C's (unconverted modified) of each covered cytosine position were then calculated from the remaining good quality alignments using Bismark methylation extractor, and the 5hmC level was calculated as #of C/(#of C+#of T).
Multiple deaminases can be combined in the same mixture to achieve sequence specificities not accessible from a single deaminase. For example, a C proceeded by C or T can be selectively deaminated by a mixture of MGYPDa917 (SEQ ID NO: 48) and NoDa01 (SEQ ID NO: 39). C followed by a G or C can be selectively deaminated by a mixture of XcDa01 (SEQ ID NO: 68), MGYPDa21 (SEQ ID NO: 64), and AcDa01 (SEQ ID NO: 49). A C followed by a T or G can be selectively deaminated by a mixture of PdDa01 (SEQ ID NO: 60) and CbDa01 (SEQ ID NO: 50). Enzymes in all three of the described mixtures are blocked by the 5ghmC modification, and in combination with TET2 and BGT would be suitable for selectively mapping C modifications in their target contexts. Libraries can be constructed, and data analyzed as in Example 4. Three base context preferences for individual deaminases is described in
For base editing, a dsDNA deaminase (e.g., CseDa01 or other deaminase from Table 3) may be split into two (or more) inactive subdomains. The breakpoint of the deaminase domain is selected such that when brought together, the deaminase domain is competent for cytosine deamination of the cytidine (the “DD” refers to this deaminase domain below). Each subdomain of the DD is genetically fused to a TALE (Transcription activator-like effector) protein with the N- to C-terminal arrangement as follows: bpNLS-TALELeft-DDN-TERM-UGI; and bpNLS-TALERight-DDC-TERM-UGI, where DD is the deaminase domain, UGI is the uracil glycosidase inhibitor from Bacillus subtilis bacteriophage, bpNLS is a bi-partite nuclear localization signal, and DDN-TERM and DDC-TERM denote the N-terminal and C-terminal subdomains of the DD. The TALE protein pair are designed to target the WTAP (Wilms tumor 1 associated protein) gene locus. Codon sequences are optimized for mammalian expression. DNA constructs encoding the TALE-base editors are placed into mammalian expression plasmids where transcription is directed by the CMV immediate-early promoter enhancer. mRNA cleavage and polyadenylation are directed by the bovine growth hormone polyadenylation signal. The plasmids are co-electroporated into HEK293 cells using a Lonza nucleofector 4D. 48 hours post electroporation, genomic DNA is extracted from the cells and the WTAP locus is PCR amplified using primers spanning the site targeted by the TALE base editor pair. Amplified products are deep sequenced, and reads are analyzed using CRISPResso2. C to T mutations are measured in a quantification window measured centered between the two TALE binding sites.
A deaminase domain (DD) of HcDa01 or other deaminase from Table 3 is genetically fused to a S. pyogenes Cas9 nickase variant D10A to encode a polypeptide with the N- to C-terminal arrangement as follows: bpNLS-DD-Cas9(D10A)-UGI-UGI-bpNLS, where DD is the deaminase domain, UGI is the uracil glycosidase inhibitor from Bacillus subtilis bacteriophage, and bpNLS is a bi-partite nuclear localization signal. Codon sequences are optimized for mammalian expression. The DNA construct encoding the genetic fusion is placed into a mammalian expression plasmid where transcription was directed by the cytomegalovirus (CMV) enhancer and the chicken beta-actin promoter. mRNA processing is directed by the chimeric intron (chicken/rabbit beta-globin) at the 5′-end of the transcript and cleavage and poly(A) tailing is directed from the rabbit beta-globin polyadenylation signal. Single guide RNAs targeting the WTAP locus are expressed from a separate DNA plasmid in which transcription of the sgRNA is directed from the U6 promoter. The sgRNA targets the sequence 5′-GGATTTAAGTGTAAATGTAC-3′ (SEQ ID NO:168). Plasmids are co-transfected into HEK293 cells. After 48 hours, genomic DNA is extracted and the WTAP locus is PCR amplified using primers spanning the site targeted by the transfected sgRNA and Cas9 base editor fusion. Amplified products are deep sequenced, and reads are analyzed using CRISPResso2. C to T mutations are measured in a quantification window measured −10 bp relative to the 3′-end of the hybridizing region of the sgRNA.
MapR may be performed, as described (Yan and Sarma, 2020 (DOI: https://doi.org/10.1002/cpmb.113, PM ID: 31943854); Yan et al., 2019 (DOI: https://doi.org/10.1016/j.celrep.2019.09.052), with the exception that RNase A may be omitted from the stop buffer. Following DNA extraction, the DNA sample is enzymatically deaminated with the ssDNA specific DNA deaminase activity of a deaminase set forth herein, e.g., HcDa01, followed by separation or removal from the reaction and/or inactivation by any means (e.g., heat, chemical, or specific or non-specific enzymatic degradation such as proteinase K digestion at 60° C. for 10 minutes). The DNA sample is purified by column purification and the eluted product used as a template for second-strand synthesis, e.g., using reagents from the NEBNext Ultra II Directional RNA Library Prep Kit (NEB E7760) following manufacturer's instructions. Libraries are PCR amplified with dual-index barcode primers for Illumina sequencing (NEB E7600) using Q5 DNA polymerase (NEB M0491) and purified. Uracil DNA glycosylase (NEB) are added to the PCR amplification mix to degrade dUTP-containing molecules and remove adaptor hairpins. Sequencing is performed on a NextSeq 500 instrument (Illumina) with 38×2 paired-end cycles.
A dsDNA deaminase with low sequence context preference, such as CseDa01, LbDa02, BaDa01, MGYPDa01, MGYPDa20, MGYPDa06, CrDa01, AvDa01, or AvDa02, is added to a dsDNA substrate, such as a plasmid, genome, or amplicon, containing the mutagenesis target, to cause base mutations resulting from deamination of one or more bases in the dsDNA substrate. Alternatively, a dsDNA deaminase with stronger sequence preference may be used to bias the mutagenesis towards or away from specific parts of the target sequence (see, for example, context preferences set forth in
To test activity of dsDNA deaminases on 5fC and 5caC modified DNA, a set of three oligonucleotide substrates (40 bp, listed below) containing four modified cytosines were prepared. The set of oligonucleotides included preferable deamination sites for eleven DNA deaminases described herein from five representative clades. In each reaction, the modified oligonucleotide (dcaC or dfC) was mixed with the control oligonucleotide (C only) in a ratio of 1:1 (800 ng+800 ng) to monitor deamination of cytosine to uracil. After incubation for 5 h at 37° C. in—reaction buffer containing 50 mM Bis-Tris pH 6.0, 0.1% Triton X-100 with different DNA deaminases, the oligonucleotide substrates were purified using Monarch PCR and DNA Cleanup kit, digested to nucleosides with the Nucleoside Digestion Mix (NEB, Ipswich MA) and the reaction products were quantified with LC-MS/MS. The results are shown in Table 5.
Oligonucleotides used to test deaminase activity on 5fC and 5caC were as follows:
“5caC” in Table 5—each of the four Cs are substituted with 5caC
“5fC” in Table 5 —each of the four Cs are substituted with 5fC
The NGS assay described in Example 3 was used to survey sequence preference among selected dsDNA deaminases. Deaminase recognition sequences sites can extend beyond the nCn context, with preferences for sequences of various lengths and compositions (
This application is a continuation-in-part of U.S. application Ser. No. 18/058,115, filed on Nov. 22, 2022, which claims priority to provisional application Ser. No. 63/264,513, filed on Nov. 24, 2021, which applications are incorporated by reference herein in their entirety.
| Number | Date | Country | |
|---|---|---|---|
| 63264513 | Nov 2021 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 18058115 | Nov 2022 | US |
| Child | 18323143 | US |
