CYTOSINE NUCLEOTIDES

[0001] CROSS-REFERENCE TO RELATED APPLICATION

[0002] This application claims the benefit of U.S. Provisional Application Serial No. 63/412,228, filed on September 30, 2022, which is incorporated by reference herein in its entirety.

[0003] SEQUENCE LISTING

[0004] This application contains a Sequence Listing electronically submitted via EFS-Web to the United States Patent and Trademark Office as an XML file entitled "0531_002463W001_SL.xml" having a size of 470 kilobytes and created on September 25, 2023. The information contained in the Sequence Listing is incorporated by reference herein.

[0005] FIELD

[0006] Embodiments of the present disclosure relate to preparing nucleic acids for sequencing, or other applications. In particular, embodiments of the proteins, methods, compositions, and kits provided herein relate to mapping of methylation status by sequencing libraries and other methods.

[0007] BACKGROUND

[0008] Modified DNA cytosines, including 5-methylcytosine (5mC) and 5 -hydroxymethyl cytosine (5hmC), are a well-studied epigenetic modification that play fundamental roles in human development and disease. Its genome-wide distribution differs between tissue types, and between healthy and diseased states. In recent years, 5mC has also gained prominence as a tool for clinical diagnostics: its distribution in cell-free DNA (cfDNA) - obtained from a liquid biopsy - can be used for the tissue-specific prediction of early-stage cancer or monitoring of cancer recurrence or remission after treatment. As a result, there has been an intense focus on developing methods for mapping modified DNA cytosines at single base resolution, with minimal loss of sample DNA quantity, quality, and complexity. Current methods for mapping modified DNA cytosines, however, exhibit limitations including (i) degradation of sample DNA due to prolonged chemical treatment at non-neutral pH and high temperatures, (ii) loss of sample DNA complexity due to conversion of unmethylated DNA bases to uracil, resulting in low complexity genome mapping, (iii) multi-step conversion, requiring both enzymes and chemical treatment, and (iv) for antibody-based 5mC detection, resolution of detection is limited to ~15Obp, precluding the identification of its exact location in the genome.

[0009] 5-hydroxymethylcytosine (5hmC) is an oxidized derivative of the widely studied epigenetic modification 5-methylcytosine (5mC). Increasing evidence supports the biological importance of 5hmC in diverse developmental processes in mammals, such as neurogenesis. As such, there is widespread interest in determining the localization of 5hmC in DNA from healthy and diseased patients. The majority of methods for mapping 5- hydroxymethylcytosine (5hmC) require bisulfite treatment, which results in significant DNA loss and damage. Recent methods for mapping 5hmC have been developed, such as oxBS-seq, TAB-seq, and ACE-seq, but some include bisulfite treatment, and all involve multiple steps using different enzymes.

[0010] SUMMARY OF THE APPLICATION

[0011] The present disclosure provides proteins, methods, compositions, and kits for determining the methylation status of DNA and RNA. The methods and compositions include a CpG binding protein. In one or more embodiments, the composition and methods include a cytidine deaminase. In one or more embodiments, the composition and methods include a ten-eleven translocase (TET) protein.

[0012] In one aspect, the present disclosure presents a one-step enzymatic method using a CpG binding protein and cytidine deaminases. In one embodiment, the cytidine deaminase is a wild-type cytidine deaminase. Wild-type cytidine deaminases act on unmodified cytosines and convert them to U and also act on modified cytosines and convert them to T. As a result, wild-type cytidine deaminases generally cannot be used to distinguish unmodified cytosines from modified cytosines. [0013] In one embodiment, the cytidine deaminase is an altered cytidine deaminase that selectively acts on certain modified cytosines of target nucleic acids and converts them to thymidine (T). The altered cytidine deaminases described herein can include a residual activity of also converting unmethylated cytosines to uracil (U), and the inventors determined that inclusion of CpG binding protein will reduce this residual activity. The methods described herein permit the use of wild-type cytidine deaminases in the mapping of modified cytosines and improves the ability of using altered cytidine deaminases in the mapping of modified cytosines.

[0014] The present disclosure also provides proteins, methods, compositions, and kits for mapping 5-hydroxymethylcytosine (5hmC) nucleotides present in DNA and RNA. The compositions and methods include a CpG binding protein and a cytidine deaminase. Current methods for mapping methylation status of 5hmC nucleotides include a step of modifying or blocking 5hmC nucleotides. For instance, the ACE-seq method (Schutsky et al., Nature biotechnology, 10.1038/nbt.4204. 8 Oct. 2018, doi : 10.1038/nbt.4204) blocks 5hmC by conversion to 5ghmC using the enzyme p-glucosyltransferase (PGT). Unlike current methods for mapping methylation status of 5hmC nucleotides, the methods presented herein do not require modifying or blocking 5hmC nucleotides. Instead, the present disclosure presents one-step, fully enzymatic methods using altered cytidine deaminases that selectively act on certain modified cytosines of target nucleic acids and converts them to uracil (U) or thymidine (T) but do not act on 5hmC, 5-formylcytosine (5fC), or 5-carboxycytosine (5-caC). The altered cytidine deaminases described herein circumvent the limitations of currently available methods for mapping 5hmC nucleotides because (i) harsh chemical treatments that cause significant loss of DNA and RNA are not required, (ii) the conversion is a single-step enzymatic reaction, and (iii) the process results in detection of 5hmC at single base resolution.

[0015] In some embodiments, the compositions and methods includes altered cytidine deaminases. In one embodiment, an altered cytidine deaminase includes amino acid substitution mutations at positions functionally equivalent to (Tyr/Phe)130 and Tyrl32 in a wild-type APOBEC3A protein. In another embodiment, an altered cytidine deaminase includes an amino acid substitution mutation at a position functionally equivalent to (Tyr/Phe)130 in a wild-type AP0BEC3A protein, where the substitution mutation is (Tyr/Phe)130Trp. The (Tyr/Phe)130 of an altered cytidine deaminase can be Tyrl30, and the wild-type AP0BEC3A protein is SEQ ID NO:3. The present disclosure also includes a polynucleotide encoding an altered cytidine deaminase.

[0016] The present disclosure also provides compositions that include an altered cytidine deaminase and a CpG binding protein. In one embodiment, a composition can further include at least one of (i) a sample including DNA including at least one modified cytosine, where the modified cytosine is 5-methyl cytosine (5mC), 5 -hydroxymethyl cytosine (5hmC), 5-formyl cytosine (5fC), 5-carboxy cytosine (5caC), or a combination thereof; or (ii) a buffer having a pH that is lower than 7; or (iii) combinations thereof.

[0017] Also provided by the present disclosure are methods. In one embodiment, a method includes providing a sample of DNA suspected of including single- stranded DNA including at least one 5-methyl cytosine (5mC), at least one 5 -hydroxymethyl cytosine (5hmC), at least one 5-formyl cytosine (5fC), at least one 5-carboxy cytosine (5CaC), or a combination thereof; contacting the single-stranded DNA with a CpG binding protein and a cytidine deaminase (e.g., a wild type or altered cytidine deaminase) under conditions suitable for (i) conversion of 5-methylcytosine (5mC) to thymidine (T) by deamination at a greater rate than conversion of cytosine (C) to uracil (U) by deamination, to result in converted single-stranded DNA, or (ii) conversion of C to U by deamination and 5mC to T by deamination at a greater rate than conversion of 5-hydroxymethyl cytosine (5hmC) to 5- hydroxymethyl uracil (5hmU) by deamination; and processing the converted singlestranded DNA to produce a sequencing library.

[0018] In another embodiment, a method includes providing a sample of DNA suspected of including double-stranded DNA including at least one 5-methyl cytosine (5mC), at least one 5-hydroxymethyl cytosine (5hmC), at least one 5-formyl cytosine (5fC), at least one 5- carboxy cytosine (5caC), or a combination thereof; processing the double-stranded DNA to produce a sequencing library; denaturing the sequencing library to result in a singlestranded DNA; contacting the single-stranded DNA with a CpG binding protein and a cytidine deaminase (e g., a wild type or altered cytidine deaminase) i) conversion of 5- methyl cytosine (5mC) to thymidine (T) by deamination at a greater rate than conversion of cytosine (C) to uracil (U) by deamination, or (ii) conversion of C to U by deamination and 5mC to T by deamination at a greater rate than conversion of 5 -hydroxymethyl cytosine (5hmC) to 5 -hydroxymethyl uracil (5hmU) by deamination, to result in converted singlestranded DNA; and converting the converted single-stranded DNA to a converted doublestranded DNA sequencing library.

[0019] In one embodiment, a method can include detecting the location of a modified cytosine in a target nucleic acid. The method can include (a) contacting target nucleic acids suspected of including at least one modified cytosine with a CpG binding protein and an altered cytidine deaminase of claim 1 or 2 to produce converted nucleic acids including at least one converted cytosine; and (b) detecting the at least one converted cytosine in the converted nucleic acids of (a). The detecting can include sequencing the converted nucleic acids or hybridizing nucleic acid probes to the converted nucleic acids.

[0020] In embodiments where the detecting includes sequencing the converted nucleic acids, the method can further include comparing the sequence of the converted nucleic acids with an untreated reference sequence to determine which cytosines in the target nucleic acids are modified.

[0021] In embodiments where the detecting includes hybridizing the converted nucleic acids to nucleic acid probes, then the nucleic acid probes can be present on an analyte array, and the method can further include sequencing the hybridized converted nucleic acids. In another embodiment the method can further include amplifying the converted nucleic acid, where the nucleic acid probes include two primers for amplification of a predetermined sequence, where the primers anneal to regions of converted nucleic acids including at least one converted cytosine with a greater affinity than to the regions of converted nucleic acids where at least one cytosine is not a converted cytosine, and where the presence of an amplified product is indicative of a modified cytosine in the target nucleic acid. In another embodiment the method can further include cleaving a single stranded DNA (ssDNA) reporter substrate by a CRISPR-based system, where the ssDNA reporter substrate includes a fluorophore and a quencher, where the presence of fluorescence is indicative of a modified cytosine in the target nucleic acid. In another embodiment the converted nucleic acids can be present in a fixed cell, where the nucleic acid probes include a fluorescent labeled probe, and where the nucleic acid probes anneal to a predetermined sequence of converted nucleic acids including at least one converted cytosine with a greater affinity than to the regions of converted nucleic acids where at least one cytosine is not a converted cytosine, where the presence of cell-associated fluorescence is indicative of a modified cytosine in the target nucleic acid.

[0022] In one embodiment of detecting the location of a modified cytosine in a target nucleic acid, the target nucleic acids can be obtained from a subject, and the detecting can include obtaining a pattern of cytosine modification in the converted nucleic acids. In some embodiments, the method can further include comparing the pattern of cytosine modification in the converted nucleic acids with the pattern of cytosine modification in a reference nucleic acid. For instance, the subject can be one that has or is at risk of having a disease or condition, and the reference nucleic acid can be from a normal subject. In one embodiment, a pattern of cytosine modification is linked in-cis to a coding region that is correlated with a disease or condition. For instance, the pattern of cytosine modification is linked in-cis to a coding region, where the coding region in the reference nucleic acid is transcriptionally active or transcriptionally inactive. The comparing can further include determining if the pattern of cytosine modification of the converted nucleic acid indicates the coding region is transcriptionally active or transcriptionally inactive in the subject. The transcription of the coding region can be correlated with a disease or condition.

[0023] For any method disclosed herein that includes discrete steps, the steps may be conducted in any feasible order. And, as appropriate, any combination of two or more steps may be conducted simultaneously.

[0024] The above summary of the present disclosure is not intended to describe each disclosed embodiment or every implementation of the present disclosure. The description that follows more particularly exemplifies illustrative embodiments. In several places throughout the application, guidance is provided through lists of examples, which examples can be used in various combinations. In each instance, the recited list serves only as a representative group and should not be interpreted as an exclusive list.

[0025] In another aspect, the present disclosure presents a method that includes oxidation of modified cytosines with, for instance, a ten-eleven translocase (TET) protein. The inventors have found that the presence of unmethylated cytosines inhibit the oxidation, and the inclusion of CpG binding protein reduces the inhibition.

[0026] Terms used herein will be understood to take on their ordinary meaning in the relevant art unless specified otherwise. Several terms used herein, and their meanings are set forth below.

[0027] Terms used herein will be understood to take on their ordinary meaning in the relevant art unless specified otherwise. Several terms used herein and their meanings are set forth below.

[0028] As used herein, the terms "organism," "subject," are used interchangeably and refer to microbes (e.g., prokaryotic or eukaryotic) animals and plants. An example of an animal is a mammal, such as a human.

[0029] As used herein, the term "target nucleic acid," is intended as a semantic identifier for the nucleic acid in the context of a method or composition or kit set forth herein and does not necessarily limit the structure or function of the nucleic acid beyond what is otherwise indicated. Reference to a nucleic acid such as a target nucleic acid includes both singlestranded and double-stranded nucleic acids, and both DNA and RNA, unless indicated otherwise. The term library refers to the collection of target nucleic acids containing known common sequences, such as a universal sequence or adapter, at their 3' and 5' ends.

[0030] As used herein, the term "adapter" and its derivatives, e.g., universal adapter, refers generally to any linear oligonucleotide which can be attached to a target nucleic acid. An adapter can be single- stranded or double-stranded DNA, or can include both doublestranded and single-stranded regions. An adapter can include a universal sequence that is substantially identical, or substantially complementary, to at least a portion of a primer, for example a universal primer; an index (also referred to herein as a barcode or tag) to assist with downstream error correction, identification, or sequencing; and/or a unique molecular identifier. In some embodiments, the adapter is substantially non-complementary to the 3' end or the 5' end of any target sequence present in the sample. In some embodiments, suitable adapter lengths are in the range of about 6-100 nucleotides, about 12-60 nucleotides, or about 15-50 nucleotides in length. For instance, the terms "adaptor" and "adapter" are used interchangeably.

[0031] As used herein, the term "universal," when used to describe a nucleotide sequence, refers to a region of sequence that is common to two or more nucleic acid molecules where the molecules also have regions of sequence that differ from each other. A universal sequence that is present in different members of a collection of nucleic acids can be used as, for instance, a "landing pad" in a subsequent step to anneal a nucleotide sequence that can be used as a primer for addition of another nucleotide sequence, such as an index, to a target nucleic acid. A universal sequence that is present in different members of a collection of nucleic acids can allow capture of multiple different nucleic acids using a population of universal capture nucleic acids, e.g., capture oligonucleotides that are complementary to a portion of the universal sequence, e.g., a universal capture sequence. Non-limiting examples of universal capture sequences include sequences that are identical to or complementary to P5 and P7 primers. Similarly, a universal sequence present in different members of a collection of molecules can allow the replication (e.g., sequencing) or amplification of multiple different nucleic acids using a population of universal primers that are complementary to a portion of the universal sequence, e.g., a universal anchor sequence. In one embodiment universal anchor sequences are used as a site to which a universal primer (e.g., a sequencing primer for read 1 or read 2) anneals for sequencing. A capture oligonucleotide or a universal primer therefore includes a sequence that can hybridize specifically to a universal sequence.

[0032] The terms "P5" and "P7" may be used when referring to a universal capture sequence or a capture oligonucleotide. The terms "P5 ¹ " (P5 prime) and "P7 ¹ " (P7 prime) refer to the complement of P5 and P7, respectively. It will be understood that any suitable universal capture sequence or a capture oligonucleotide can be used in the methods presented herein, and that the use of P5 and P7 are exemplary embodiments only. Uses of capture oligonucleotides such as P5 and P7 or their complements on flow cells are known in the art, as exemplified by the disclosures of WO 2007/010251, WO 2006/064199, WO 2005/065814, WO 2015/106941, WO 1998/044151, and WO 2000/018957 , which are incorporated by reference as to P5 and P7 and their uses. For example, any suitable forward amplification primer, whether immobilized or in solution, can be useful in the methods presented herein for hybridization to a complementary sequence and amplification of a sequence. Similarly, any suitable reverse amplification primer, whether immobilized or in solution, can be useful in the methods presented herein for hybridization to a complementary sequence and amplification of a sequence. One of skill in the art will understand how to design and use primer sequences that are suitable for capture and/or amplification of nucleic acids as presented herein.

[0033] As used herein, the term "primer" and its derivatives refer generally to any nucleic acid that can hybridize to a target sequence of interest. Typically, the primer functions as a substrate onto which nucleotides can be polymerized by a polymerase or to which a polynucleotide can be ligated; in some embodiments, however, the primer can become incorporated into the synthesized nucleic acid strand and provide a site to which another primer can hybridize to prime synthesis of a new strand that is complementary to the synthesized nucleic acid molecule. In some embodiments, the primer can be used for hybridization to a predetermined sequence, for instance a predetermined sequence that includes one or more nucleotides that identify the location of a modified cytosine. In one embodiment, a “primer” includes a sequence present in a guide RNA used with a CRISPR-based system to hybridize to a predetermined sequence. The primer can include any combination of nucleotides or analogs thereof. In some embodiments, the primer is a single-stranded oligonucleotide or polynucleotide.

[0034] The terms "polynucleotide," "oligonucleotide" and “nucleic acid” are used interchangeably herein to refer to a polymeric form of nucleotides of any length, and may include ribonucleotides, deoxyribonucleotides, analogs thereof, or mixtures thereof. The terms should be understood to include, as equivalents, analogs of either DNA, RNA, cDNA, or antibody-oligo conjugates made from nucleotide analogs and to be applicable to single stranded (such as sense or antisense) and double stranded polynucleotides. The term as used herein also encompasses cDNA, that is complementary or copy DNA produced from a RNA template, for example by the action of reverse transcriptase.

[0035] As used herein, the term "protein" refers broadly to a polymer of two or more amino acids joined together by peptide bonds. The term "protein" also includes molecules which contain more than one polypeptide joined by disulfide bonds, ionic bonds, or hydrophobic interactions, or complexes of polypeptides that are joined together, covalently or noncovalently, as multimers (e.g., dimers, tetramers). Thus, the terms peptide, oligopeptide, and polypeptide are all included within the definition of protein and these terms are used interchangeably. It should be understood that these terms do not connote a specific length of a polymer of amino acids, nor are they intended to imply or distinguish whether the protein is produced using recombinant techniques, chemical or enzymatic synthesis, or is naturally occurring.

[0036] A protein or polynucleotide can be isolated. An "isolated" protein or polynucleotide is one that has been removed from a cell. For instance, an isolated protein is a polypeptide that has been removed from the cytoplasm or from the membrane of a cell, and many of the proteins, nucleic acids, and other cellular material of its natural environment are no longer present. A protein or polynucleotide produced by chemical or enzymatic synthesis is understood to be isolated.

[0037] As used herein, an "index" (also referred to as an "index region," "index adaptor," "tag," or a "barcode") refers to a unique nucleic acid tag that can be used to identify a sample or source of the nucleic acid material, or a compartment in which a target nucleic acid was present. The index can be present in solution or on a solid-support, or attached to or associated with a solid-support and released in solution or compartment. When nucleic acid samples are derived from multiple sources, the nucleic acids in each nucleic acid sample can be tagged with different nucleic acid tags such that the source of the sample can be identified. Any suitable index or set of indexes can be used, as known in the art and as exemplified by the disclosures of U.S. Pat. No. 8,053,192, PCT Publication No. WO 05/068656, and U.S. Pat. Publication No. 2013/0274117. In some embodiments, an index can include a six-base Index 1 (i7) sequence, an eight-base Index 1 (i7) sequence, an eight- base Index 2 (i5e) sequence, a ten-base Index 1 (i7) sequence, or a ten-base Index 2 (i5) sequence from Illumina, Inc. (San Diego, CA).

[0038] As used herein, the term "amplicon," when used in reference to a nucleic acid, means the product of copying the nucleic acid, wherein the product has a nucleotide sequence that is the same as or complementary to at least a portion of the nucleotide sequence of the nucleic acid. An amplicon can be produced by any of a variety of amplification methods that use the nucleic acid, or an amplicon thereof, as a template including, for example, polymerase extension, polymerase chain reaction (PCR), rolling circle amplification (RCA), ligation extension, or ligation chain reaction. An amplicon can be a nucleic acid molecule having a single copy of a particular nucleotide sequence (e.g. a PCR product) or multiple copies of the nucleotide sequence (e g. a concatameric product of RCA). A first amplicon of a target nucleic acid is typically a complementary copy. Subsequent amplicons are copies that are created, after generation of the first amplicon, from the target nucleic acid or from the first amplicon. A subsequent amplicon can have a sequence that is substantially complementary to the target nucleic acid or substantially identical to the target nucleic acid.

[0039] As used herein, "amplify", "amplifying" or "amplification reaction" and their derivatives, refer generally to any action or process whereby at least a portion of a nucleic acid molecule is replicated or copied into at least one additional nucleic acid molecule. The additional nucleic acid molecule optionally includes sequence that is substantially identical or substantially complementary to at least some portion of the template nucleic acid molecule. The template nucleic acid molecule can be single-stranded or double-stranded and the additional nucleic acid molecule can independently be single-stranded or doublestranded. Amplification is typically the exponential replication of a nucleic acid molecule. In some embodiments, such amplification can be performed using isothermal conditions; in other embodiments, such amplification can include thermocycling. In some embodiments, the amplification is a multiplex amplification that includes the simultaneous amplification of a plurality of target sequences in a single amplification reaction. In some embodiments, "amplification" includes amplification of at least some portion of DNA and RNA based nucleic acids alone, or in combination. The amplification reaction can include any of the amplification processes known to one of ordinary skill in the art. In some embodiments, the amplification reaction includes polymerase chain reaction (PCR).

[0040] As used herein, the term "polymerase chain reaction" ("PCR") refers to the method of Mullis U.S. Pat. Nos. 4,683,195 and 4,683,202, which describe a method for increasing the concentration of a segment of a polynucleotide of interest in a mixture of genomic DNA without cloning or purification. This process for amplifying the polynucleotide of interest consists of introducing a large excess of two oligonucleotide primers to the DNA mixture containing the desired polynucleotide of interest, followed by a series of thermal cycling in the presence of a DNA polymerase. The two primers are complementary to their respective strands of the double stranded polynucleotide of interest. The mixture is denatured at a higher temperature first and the primers are then annealed to complementary sequences within the polynucleotide of interest molecule. Following annealing, the primers are extended with a polymerase to form a new pair of complementary strands. The steps of denaturation, primer annealing and polymerase extension can be repeated many times (referred to as thermocycling) to obtain a high concentration of an amplified segment of the desired polynucleotide of interest. The length of the amplified segment of the desired polynucleotide of interest (amplicon) is determined by the relative positions of the primers with respect to each other, and therefore, this length is a controllable parameter. By virtue of repeating the process, the method is referred to as PCR. Because the desired amplified segments of the polynucleotide of interest become the predominant nucleic acid sequences (in terms of concentration) in the mixture, they are said to be "PCR amplified". In a modification to the method discussed above, the target nucleic acid molecules can be PCR amplified using a plurality of different primer pairs, in some cases, one or more primer pairs per target nucleic acid molecule of interest, thereby forming a multiplex PCR reaction.

[0041] As used herein, "amplification conditions" and its derivatives, generally refers to conditions suitable for amplifying one or more nucleic acid sequences. In some embodiments, the amplification conditions can include isothermal conditions or alternatively can include thermocycling conditions, or a combination of isothermal and thermocycling conditions. In some embodiments, the conditions suitable for amplifying one or more nucleic acid sequences include polymerase chain reaction (PCR) conditions. Typically, the amplification conditions refer to a reaction mixture that is sufficient to amplify nucleic acids such as one or more target sequences flanked by a universal sequence, or target specific primers, or to amplify an amplified target sequence flanked by one or more adapters. Generally, the amplification conditions include a catalyst for amplification or for nucleic acid synthesis, for example a polymerase; a primer that possesses some degree of complementarity to the nucleic acid to be amplified; and nucleotides, such as deoxyribonucleotide triphosphates (dNTPs) to promote extension of the primer once hybridized to the nucleic acid. The amplification conditions can require hybridization or annealing of a primer to a nucleic acid, extension of the primer and a denaturing step in which the extended primer is separated from the nucleic acid sequence undergoing amplification. Typically, but not necessarily, amplification conditions can include thermocycling; in some embodiments, amplification conditions include a plurality of cycles where the steps of annealing, extending and separating are repeated. Typically, the amplification conditions include cations such as Mg ²⁺ or Mn ²⁺ and can also include various modifiers of ionic strength.

[0042] As defined herein "multiplex amplification" refers to selective and non-random amplification of two or more target sequences within a sample using at least one targetspecific primer. In some embodiments, multiplex amplification is performed such that some or all of the target sequences are amplified within a single reaction vessel. The "plexy" or "plex" of a given multiplex amplification refers generally to the number of different target-specific sequences that are amplified during that single multiplex amplification. In some embodiments, the plexy can be about 12-plex, 24-plex, 48-plex, 96- plex, 192-plex, 384-plex, 768-plex, 1536-plex, 3072-plex, 6144-plex or higher. It is also possible to detect the amplified target sequences by several different methodologies (e.g., gel electrophoresis followed by densitometry, quantitation with a bioanalyzer or quantitative PCR, hybridization with a labeled probe; incorporation of biotinylated primers followed by avidin-enzyme conjugate detection; incorporation of ³² P -labeled deoxynucleotide triphosphates into the amplified target sequence). [0043] As used herein, the term "amplification site" refers to a site in or on an array where one or more amplicons can be generated. An amplification site can be further configured to contain, hold or attach at least one amplicon that is generated at the site.

[0044] As used herein, the term "array,” “analyte array,” and “microarray” are used interchangeably and refer to a population of sites that can be differentiated from each other according to relative location. Different molecules that are at different sites of an array can be differentiated from each other according to the locations of the sites in the array. An individual site of an array can include one or more molecules of a particular type. For example, a site can include a single target nucleic acid molecule having a particular sequence or a site can include several nucleic acid molecules having the same sequence (and/or complementary sequence, thereof). The sites of an array can be different features located on the same substrate. Exemplary features include without limitation, droplets, wells in a substrate, beads (or other particles) in or on a substrate, projections from a substrate, ridges on a substrate or channels in a substrate. The sites of an array can be separate substrates each bearing a different molecule. Different molecules attached to separate substrates can be identified according to the locations of the substrates on a surface to which the substrates are associated or according to the locations of the substrates in a liquid or gel. Exemplary arrays in which separate substrates are located on a surface include, without limitation, those having beads in wells.

[0045] As used herein, the term "compartment" is intended to mean an area or volume that separates or isolates something from other things. Exemplary compartments include, but are not limited to, vials, tubes, wells, droplets, boluses, beads, vessels, surface features, flow cell, or areas or volumes separated by physical forces such as fluid flow, magnetism, electrical current or the like. In one embodiment, a compartment is a well of a multi-well plate, such as a 96- or 384-well plate. As used herein, a droplet may include a hydrogel bead, which is a bead for encapsulating one or more nuclei or cell, and includes a hydrogel composition. In some embodiments, the droplet is a homogeneous droplet of hydrogel material or is a hollow droplet having a polymer hydrogel shell. Whether homogenous or hollow, a droplet may be capable of encapsulating one or more nuclei or cells. In some embodiments, the droplet is a surfactant stabilized droplet. In some embodiments, a single cell or Nuclei is present per compartment. In some embodiments, two or more cells or Nuclei are present per compartment. In some embodiments, each compartment contains a compartment-specific index. In some embodiments, the index is in solution or attached or associated with a solid-phase in each compartment.

[0046] The term "flow cell" as used herein refers to a chamber comprising a solid surface across which one or more fluid reagents can be flowed. Examples of flow cells and related fluidic systems and detection platforms that can be readily used in the methods of the present disclosure are described, for example, in Bentley et al., Nature 456:53-59 (2008), WO 04/018497; US 7,057,026; WO 91/06678; WO 07/123744; US 7,329,492; US 7,211,414; US 7,315,019; US 7,405,281; and US 2008/0108082.

[0047] As used herein, the term "clonal population" refers to a population of nucleic acids that is homogeneous with respect to a particular nucleotide sequence. The homogenous sequence is typically at least 10 nucleotides long, but can be even longer including for example, at least 50, 100, 250, 500 or 1000 nucleotides long. A clonal population can be derived from a single target nucleic acid or template nucleic acid. Typically, all of the nucleic acids in a clonal population will have the same nucleotide sequence. It will be understood that a small number of mutations (e.g., due to amplification artifacts) can occur in a clonal population without departing from clonality.

[0048] As used herein, a "pattern of cytosine modification," also referred to as a "methylation profile," refers to the pattern with which both methylation and unmethylation of cysteines is distributed in the genome of a cell or an organism. A “pattern” is inclusive of both modified cytosines and non-modified cytosines. The pattern can be defined in several distribution dimensions: by organ, by tissue, by status of disease or pathological condition (e g., cancer, neurophysiological), by genome segment (e.g., chromosome or genetic coordinates on a chromosome), by gene, by CpG island, a group of cytosines, or by the site of a modified cytosine. A pattern of cytosine modification can have a known correlation with a disease or pathological condition, or correlation of a pattern of cytosine modification with a disease or pathological condition can be identified using methods described herein. A pattern of cytosine modification can be present at a specific locus (e.g., location) in a genome, and that specific location can be a single modified cytosine or a set of modified cytosines, e.g., a CpG island. A pattern of cytosine modification can be identified by using a predetermined sequence, e.g., a method of using an cytidine deaminase can be designed and practiced with the intent of determining a pattern of cytosine modification, for instance, the methylation status of one of more specific cytosines, the methylation status of one or more specific cytosines present at a specific location of a genome, or the combination thereof.

[0049] As used herein, the term "each," when used in reference to a collection of items, is intended to identify an individual item in the collection but does not necessarily refer to every item in the collection unless the context clearly dictates otherwise.

[0050] As used in this specification and the appended claims, the term "or" is generally employed in its sense including "and/or" unless the content clearly dictates otherwise. The term "and/or" means one or all of the listed elements or a combination of any two or more of the listed elements. The use of "and/or" in some instances does not imply that the use of "or" in other instances may not mean "and/or."

[0051] Unless otherwise specified, "a," "an," "the," and "at least one" are used interchangeably and mean one or more than one.

[0052] As used in this specification and the appended claims, the term "or" is generally employed in its sense including "and/or" unless the content clearly dictates otherwise. The term "and/or" means one or all of the listed elements or a combination of any two or more of the listed elements. The use of "and/or" in some instances does not imply that the use of "or" in other instances may not mean "and/or."

[0053] The words "preferred" and "preferably" refer to embodiments of the disclosure that may afford certain benefits, under certain circumstances. However, other embodiments may also be preferred, under the same or other circumstances. Furthermore, the recitation of one or more preferred embodiments does not imply that other embodiments are not useful, and is not intended to exclude other embodiments from the scope of the disclosure. [0054] As used herein, "have," "has," "having," "include," "includes," "including," "comprise," "comprises," "comprising" or the like are used in their open ended inclusive sense, and generally mean "include, but not limited to," "includes, but not limited to," or "including, but not limited to."

[0055] It is understood that wherever embodiments are described herein with the language "have," "has," "having," "include," "includes," "including," "comprise," "comprises," "comprising" and the like, otherwise analogous embodiments described in terms of "consisting of and/or "consisting essentially of are also provided. The term "consisting of means including, and limited to, whatever follows the phrase "consisting of." That is, "consisting of indicates that the listed elements are required or mandatory, and that no other elements may be present. The term "consisting essentially of indicates that any elements listed after the phrase are included, and that other elements than those listed may be included provided that those elements do not interfere with or contribute to the activity or action specified in the disclosure for the listed elements.

[0056] Conditions that are "suitable" for an event to occur, such as converting 5 methyl cytosine to thymidine by deamination, or "suitable" conditions are conditions that do not prevent such events from occurring. Thus, these conditions permit, enhance, facilitate, and/or are conducive to the event.

[0057] As used herein, "providing" in the context of a protein, sample of DNA or RNA, or composition means making the protein, sample of DNA or RNA, or composition, purchasing the protein, sample of DNA or RNA, or composition, or otherwise obtaining the protein, sample of DNA or RNA, or composition.

[0058] Reference throughout this specification to "one embodiment," "an embodiment," "certain embodiments," or "some embodiments," etc., means that a particular feature, configuration, composition, or characteristic described in connection with the embodiment is included in at least one embodiment of the disclosure. Thus, the appearances of such phrases in various places throughout this specification are not necessarily referring to the same embodiment of the disclosure. Furthermore, the particular features, configurations, compositions, or characteristics may be combined in any suitable manner in one or more embodiments.

[0059] While polynucleotide sequences encoding an cytidine deaminase are described herein as DNA sequences, it is understood that the complements, reverse sequences, and reverse complements of the DNA sequences can be easily determined by the skilled person. It is also understood that the sequences described herein as DNA sequences can be converted from a DNA sequence to an RNA sequence by replacing each thymidine nucleotide with a uracil nucleotide.

[0060] Throughout this disclosure, various aspects of the disclosure can be presented in a range format. It should be understood that the description in range format is merely for convenience and brevity and should not be construed as an inflexible limitation on the scope of the disclosure. Accordingly, the description of a range should be considered to have specifically disclosed all the possible subranges as well as individual numerical values within that range. For example, description of a range such as from 1 to 6 should be considered to have specifically disclosed subranges such as from 1 to 3, from 1 to 4, from 1 to 5, from 2 to 4, from 2 to 6, from 3 to 6 etc., as well as individual numbers within that range, for example, 1, 2, 2.7, 3, 4, 4.5, 5, 5.3, and 6. This applies regardless of the breadth of the range.

[0061] In the description herein particular embodiments may be described in isolation for clarity. Unless otherwise expressly specified that the features of a particular embodiment are incompatible with the features of another embodiment, certain embodiments can include a combination of compatible features described herein in connection with one or more embodiments.

[0062] For any method disclosed herein that includes discrete steps, the steps may be conducted in any feasible order. And, as appropriate, any combination of two or more steps may be conducted simultaneously.

[0063] The above summary of the present disclosure is not intended to describe each disclosed embodiment or every implementation of the present disclosure. The description that follows more particularly exemplifies illustrative embodiments. In several places throughout the application, guidance is provided through lists of examples, which examples can be used in various combinations. In each instance, the recited list serves only as a representative group and should not be interpreted as an exclusive list.

[0064] BRIEF DESCRIPTION OF THE FIGURES

[0065] The following detailed description of illustrative embodiments of the present disclosure may be best understood when read in conjunction with the following drawings.

[0066] FIG. 1A-F shows the deamination scheme of APOBEC3A. Cytosine (C), 5-methylcytosine (5mC), and 5-hydroxymethylcytosine (5hmC) nucleobases in single-stranded DNA are well characterized substrates of APOBEC3A. FIG. 1A shows conversion of C to uracil (U) by APOBEC3A. FIG. IB shows conversion of 5mC to thymidine (T) by APOBEC3A. FIG. 1C shows conversion of 5hmC to 5 -hydroxymethyl uracil (5hmU) by APOBEC3A. " denotes p _ie connection of the nucleobases to a DNA molecule. FIG. 1D-F shows the result of treating a DNA sample with a wild-type APOBEC3A enzyme (FIG. ID), an example of one-step detection of 5mC using an altered cytidine deaminase described herein (FIG. IE), and an example of one-step detection of 5hmC using a different altered cytidine deaminase described herein (FIG. IF). The top strand of FIG. 1D-F shows C, 5mC, and/or 5hmC bases, and the bottom strand of FIG. ID, IE, and IF underlines the changed bases. In FIG. ID, 5mC nucleobases are marked with CH3, 5hmC nucleobases are marked with CH2-OH, 5-hydroxymethyl uracil nucleobases are designated with small case "u" and uracil nucleobases are designated with capital "U."

[0067] FIG. 2 shows cytosine and examples of modified cytosine nucleobases in DNA. The denotes the connection of the nucleobases to a DNA molecule.

[0068] FIG. 3 is a schematic showing alignment of cytidine deaminase amino acid sequences using the Clustal O algorithm. An "*" (asterisk) indicates positions which have a single, fully conserved residue between all cytidine deaminases. A (colon) indicates conservation between groups of strongly similar properties as below - roughly equivalent to scoring > 0.5 in the Gonnet PAM 250 matrix. A (period) indicates conservation between groups of weakly similar properties as below - roughly equivalent to scoring =< 0.5 and > 0 in the Gonnet PAM 250 matrix. The amino acids marked with " ^A " show the ZDD motif SEQ ID NO: 12 (e.g., above amino acids 70 to 106 of sp|P3194111-199). The amino acids marked with " ^A " and "#" show the ZDD motif SEQ ID NO: 13 (e.g., above amino acids 70 to 153 of sp|P31941|l-199). sp|P3194111-199 is a human APOBEC3A, SEQ ID NO:3; XP 045219544.1 is an APOBEC3A from Macaca fascicularis, SEQ ID NO: 19; AER45717.1 is an APOBEC3A from Pongo pygmaeus, SEQ ID NO:20; XP_003264816.1 is an APOBEC3A from Nomascus leucogenys, SEQ ID NO:21;

PNI48846.1 is an APOBEC3A from Pan troglodytes, SEQ ID NO:22; and ADO85886.1 is an APOBEC3A from Gorilla gorilla, SEQ ID NO:23.

[0069] FIG. 4 shows a schematic of restriction enzyme (SirM) -based assay for deamination by cytidine deaminase. X is H in the context of C, and X is methyl in the context of 5mC, and X is hydroxymethyl in the context of 5hmC. “Matched oligos” refers to complete complementarity between the two oligonucleotides; “mismatched oligos” refers to incomplete complementarity between the two oligonucleotides. A mismatch results in cleavage of the double stranded oligonucleotides by b'wal.

[0070] FIG. 5A-B shows a positive control experiment of N l-based assay using synthesized oligonucleotides. FIG. 5A shows sequences of synthesized oligonucleotides. oLB1609, SEQ ID NO:24; 0LBI6IO, SEQ ID NO:25; 0LBI611, SEQ ID NO:26; oLB1612, SEQ ID NO:27; OLB1679, SEQ ID NO:28; oJT1910, SEQ ID NO:90; and oJT1911, SEQ ID NO:91. FIG. 5B shows visualization of results of 5 I digestion. “Matched” refers to complete complementarity between the two oligonucleotides; “mismatched” refers to incomplete complementarity between the two oligonucleotides. A mismatch results in cleavage of the double stranded oligonucleotides by biral.

[0071] FIG. 6 shows a positive control deamination experiment using a commercially available APOBEC3A enzyme.

[0072] FIG. 7A-B shows a SDS-PAGE panel of APOBEC3A(Y130X) proteins. FIG. 7A shows SDS-PAGE analysis of purified APOBEC3A(Y130X) mutant proteins. FIG. 7B shows SDS-PAGE analysis of purified APOBEC3A(Y130A_Y132H) mutant proteins. [0073] FTG. 8 shows results of APOBEC3A(Y130X) deamination end-point assay panel using Svral assay readout.

[0074] FIG. 9 shows a bar graph representation of APOBEC3A(Y130X) deaminase activity.

[0075] FIG. 10 shows deaminase activity of wild type (NEB APOBEC) and mutant APOBEC variants on C, 5mC, and 5hmC substrates. Percent deamination values were determined from a Swal restriction enzyme assay and quantified as in FIG. 4. C deamination activity was measured in two independent experiments corresponding to the left and right panels.

[0076] FIG. 11A-F shows APOBEC3A(Y130A) time course C and 5mC deamination reaction.

FIG. 11A-C shows the reaction at 37°C, and FIG. 11D-F shows the reaction at 22°C. FIG. HA and FIG. 11D) Y130A deamination Swal-based assay at 37°C and 22°C, respectively, visualized by 15% Urea-PAGE with FAM filter. FIG. 11B and FIG. HE) Graphical representation of Y130A time course C and 5mC deamination reaction at 37°C and 22°C, respectively. FIG. 11C and FIG. 11F) Table depicting percentage deamination at 37°C and 22°C, respectively, at various time points.

[0077] FIG. 12 shows preliminary Michaelis-Menten kinetics of Y130A on C and 5mC oligonucleotide substrates.

[0078] FIG. 13 shows DNA oligonucleotide substrates for evaluating deaminase activity of a double mutant cytidine deaminase. Set(A), SEQ ID NO:29; Set(B), SEQ ID NO:30, 31, 32, and 36, respectively.

[0079] FIG. 14 shows percent deamination at each NCN motif in the DNA oligo substrate (A) after incubation (37°C, 6 hour reaction) with APOBEC3A mutants. This metric is calculated as the percentage of OT (cytosine to thymidine) mutations at each position as determined by DNA sequencing.

[0080] FIG. 15 shows percent deamination at each NCN motif in the DNA oligo substrate (A) after incubation (37°C, 1 hour reaction) with APOBEC3A mutants. This metric is calculated as the percentage of C>T mutations at each position. [0081] FTG. 16 shows percent deamination at each NCpGN motif in DNA oligo substrate (B) after incubation with APOBEC3 A mutants. Four different DNA oligos were mixed together as a substrate for APOBEC3A deamination. Percent deamination was calculated as the percentage of OT mutations in each NCpGN motif. Methylated and unmethylated forms of each NCpGN motif were assayed (32 sites in total).

[0082] FIG. 17A-B shows the time course of APOBEC3A(Y130W) deamination against C, 5mC, and 5hmC-containing substrates. The Swal restriction enzyme assay was performed to measure deaminase activity of APOBEC3A(Y130W). Percent deamination is calculated as the ratio of the cut band intensity to the (uncut + cut) band intensity. FIG.17A, gel images and FIG. 17B, quantification of band intensities.

[0083] FIG. 18A-B shows comparison of APOBEC3A deaminase activity against C, 5mC, and oxidized derivatives. The NM ZI restriction enzyme assay was performed with a 90 min reaction time to measure deaminase activities of APOBEC3A wild type and Y130W mutant enzymes. A no-protein control was included to account for potential degradation of oligonucleotide substrates and account for non-specific activity of Swal during the course of the assay. Percent deamination was calculated as the ratio of the cut band intensity to the (uncut + cut) band intensity. FIG. 18A, gel images and FIG. 18B, quantification of band intensities, subtracted from the corresponding no-protein control lanes.

[0084] FIG. 19 shows a method for 5hmC detection using deaminase-based sequencing. Treatment of a polynucleotide (SEQ ID NO: 110) with pGT results in a polynucleotide (SEQ ID NO: 113). Treatment of the polynucleotide SEQ ID NO: 110 with a selective deaminase results in the polynucleotide of SEQ ID NO: 11, which is read as SEQ ID NO: 112 during sequencing.

[0085] FIG. 20A-B shows comparison of different reaction conditions and the resulting methylation reporting on pUC19 (CG methylated) and Lambda (fully unmethylated) with the altered cytidine deaminase having Y130A and Y132H. FIG. 20A is the methylation level of methylated pUC19 and unmethylated lambda DNA using the altered cytidine deaminase under different concentrations. FIG. 20B is the methylation level of methylated pUC19 and unmethylated lambda DNA using the altered cytidine deaminase under different buffers.

[0086] FIG. 21 shows impact of RNAse A on deamination of methylated pUC19 and unmethylated lambda as determined by sequencing.

[0087] FIG. 22A-B shows activity of APOBEC Y130A-Y132H on 5hmC. (FIG. 22A) Strategy for construction of control oligo; (FIG. 22B) Observed methylation level on control oligo substrate. mC, 5mC; hmC, 5hmC.

[0088] FIG. 23 shows analysis of regional methylation in CpG islands indicate that deaminase- based assay produces the expected methylation profile. mC-deaminase-seq, APOBEC3A Y130A-Y132H.

[0089] FIG. 24 shows performance of SNV and indel calling with and without methylation conversion with APOBEC3A Y130A-Y132H. mC-deaminase-seq, APOBEC3A Y130A- Y132H.

[0090] FIG. 25A-F shows visualization of DMRs identified by EM-Seq™ and the altered cytidine deaminase assay in ZNF154 gene. Representation of methylation level across this region is shown in (FIG. 25A) HCC2218-Normal, EM-Seq™ conversion, (FIG. 25B) HCC2218- Tumor, EM-Seq™ conversion, (FIG. 25C) HCC2218-Normal, altered cytidine deaminase assay, (FIG. 25D) HCC2218-Tumor, altered cytidine deaminase assay, (FIG. 25E) differentially methylated regions (DMRs) called between Tumor/Normal samples using EM-Seq™ data, and (FIG. 25F) differentially methylated regions (DMRs) called between Tumor/Normal sample s using altered cytidine deaminase data.

[0091] FIG. 26 shows recall and precision of DMRs in HCC1187 Tumor/Normal and HCC2218 Tumor/Normal paired genomes. DMRs from each workflow were compared to DMRs called by EM-Seq™, which was used as the truth set. BiSulfite identifies a majority of the DMRs that EM-Seq™ identifies. Additionally, the mC-selective deamination protocols described here in Method A, Method B, and Method C are able to identify most of the DMRs identified by EM-Seq™. mC-deaminase-seq, APOBEC3A Y130A-Y132H. [0092] FTG. 27A-B shows methylation levels detected in promoter regions using Method A and EM-Seq™. (FIG. 27A) Methylation level at H3K36me3 regions, which are expected to be hypermethylated, (FIG. 28B) Methylation levels at H3K27ac regions, which are expected to be hypomethylated. Dotted traces, methylation levels detected in promoter regions using Method A; solid traces, methylation levels detected in promoter regions using EM-Seq™.

[0093] FIG. 28 shows tumor signal for 0% spike-in of tumor DNA into normal DNA vs a 10% spike-in of tumor DNA into normal DNA. Methylation level of HCC2218 normal DNA at individual CpG sites within the PanSeer cancer panel were assessed to create a baseline. Methylation level at individual CpG sites were then assessed in separate replicates of HCC2218 normal DNA and a 10% spike-in of HCC2218 tumor DNA into HCC2218 normal DNA. Tumor signal indicates the fraction of CpGs that have a significantly different methylation level compared to the background.

[0094] FIG. 29 shows a diagram depicting the benefits of 5mC>T conversion for enrichment.

[0095] FIG. 30A-C shows different workflows for enrichment of methyl -converted libraries.

(FIG. 30A) Hybridization after conversion and amplification, in which specialized probe designs shown in FIG. 29 are typically employed. (FIG. 30B) Illumina Methyl Capture EPIC Workflow where standard probe designs are employed prior to bisulfite conversion, requiring higher DNA inputs. (FIG. 30C) Workflow for data presented for use with an altered cytidine deaminase (mC-deaminase) with enrichment.

[0096] FIG. 31A-B shows enrichment performance of altered cytidine deaminase (mC-deaminase- seq) libraries. (FIG. 31A) Read enrichment performance of altered cytidine deaminase libraries as compared to libraries with no methylation conversion. (FIG. 31B) Correlation of regional methylation levels in CpG islands measured in the unenriched (WGS) sample as compared to the enriched sample.

[0097] FIG. 32 shows a schematic of qPCR-based detection of 5mC in a genomic locus of interest. Selective 5mC deamination by APOBEC3A(Y130A/Y132H) results in a DNA template that is fully complementary to qPCR primers, and can therefore be amplified. No deamination by APOBEC3A(Y130A/Y132H) will be observed in unmethylated substrates, resulting in qPCR primers being mismatched to the target site, therefore compromising amplification. “Selective 5mC -> T deamination” refers to the result of incubating the ssDNA substrate with the altered cytidine deaminase AP0BEC3A Y130A/Y132H, and the underlined thymidine nucleotides are the result of selective deamination of 5mC. “Non- selective 5mC -> T and C-> U deamination” refers to the result of incubating the ssDNA substrate with a wild type AP0BEC3A, and the uracil and underlined thymidine nucleotides are the result of non-selective deamination by wild type AP0BEC3A.

[0098] FIG. 33 demonstrates qPCR assay using purified APOBEC proteins show decrease in Cq value after treatment of methylated ssDNA substrate with Y130A Y132H. NEB APOBEC, APOBEC protein from New England Biolabs; Y130A, Y130A/Y132H, Y130A/E72A, and Y130A/Y132/ E72A are the substitution mutations present in four APOBEC3A proteins. The E72A mutation abrogates APOBEC activity and is used as a negative control to show observed Cq value differences are due to the mutant APOBEC enzymes.

[0099] FIG. 34 shows CRISPR-Casl2 mediated detection of 5mC. Altered cytidine deaminase - mediated conversion of 5mC to T restores full complementarity of the Casl2 guide RNA for its DNA target, allowing the Casl2-guide RNA protein complex to engage the converted substrate. This activates the collateral cleavage activity of Casl2, resulting in cleavage of the reporter ssDNA containing a fluorophore and quencher. Liberation of the fluorophore leads to increased fluorescence, which is measured in a standard fluorimeter. F, fluorophore; q, quencher.

[00100] FIG. 35A-C shows incubation of substrate DNA with APOBEC3A(Y130A/Y132H) results in high 5mC deamination and low but detectable C deamination. The ssDNA oligo substrate detailed in FIG. 32 was treated with the enzyme and subsequently analyzed by Illumina sequencing. The % methylation at each C or 5mC site in the oligo is calculated as the % conversion to T. In this experiment, varying concentrations of APOBEC3A(Y130A/Y132H) were tested at different reaction temperatures and times. An elevated level of C deamination was observed at 25C, with increased enzyme concentration and reaction time. Ih, 3h, and 6h refer to the hours of incubation time; 0.75uM, 1 ,5uM, and 4uM refer to the micromolar amounts of the enzyme used in each reaction; and 25C, 30C, and 37C refer to temperature of the reaction. In each histogram the 17 bars on the left are C deamination, and the 16 bars on the right are 5mC deamination. The nucleotide triplets on the X-axis of the histograms for each micromolar amount of enzyme are the following: ACT, ACT, CCT, TCT, ACA, GCA, CCA, GCT, GCC, ACG, ACG, ACT, TCG, CCT, TCG, GCA, ACG, TCA, TCC, ACG, ACA, ACG, GCC, ACG, TCG, CCA, GCA, GCG, TCG, GCC, ACA, GCG, TCA.

[00101] FIG. 36 shows detection of 5mC in fixed cell or tissue preparations using FISH. Fixed biological samples are permeabilized and denatured to render DNA accessible to an altered cytidine deaminase. Enzymatic deamination selectively converts 5mC to T. Methylation events are detected using fluorescent probes that are specific to the converted DNA sequence.

[00102] FIG. 37A to 37N show amino acid sequences of SEQ ID NOs:3, 16, 17, 37-66, and 90-93.

[00103] FIG. 38 shows a schematic illustration of generation of false positives by wild-type APOBEC in the absence of CpG binding protein and the reduction of false positives in the presence of CpG binding protein.

[00104] FIG. 39 shows a schematic illustration of the inhibition of TET protein activity by unmethylated CpG dinucleotides and the reduction of inhibition by CpG binding protein.

[00105] FIG. 40 shows cytosine and examples of modified cytosine nucleobases in DNA, and the iterative conversion of 5mC to 5caC catalyzed by TET protein.

[00106] The schematic drawings are not necessarily to scale. Like numbers used in the figures refer to like components, steps and the like. However, it will be understood that the use of a number to refer to a component in a given figure is not intended to limit the component in another figure labeled with the same number. In addition, the use of different numbers to refer to components is not intended to indicate that the different numbered components cannot be the same or similar to other numbered components.

[00107] DETAILED DESCRIPTION [00108] Described herein are methods for mapping modified cytosine nucleotides, such as 5- methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), at single base resolution. The methods include the use of proteins, referred to herein as "CpG binding proteins," that selectively bind to unmethylated CpG dinucleotides.

[00109] In one aspect, a CpG binding protein is used with methods that include a cytidine deaminase, such as an APOBEC protein. Wild-type cytidine deaminases cannot be used for distinguishing unmethylated cytosines from methylated cytosines because they convert 5mC to thymidine (T) and cytosine (C) to uracil (U). Altered cytidine deaminases described herein preferentially convert 5mC to thymidine (T) compared to cytosine (C) to uracil (U), but the conversion of C to U can lead to false positives (FIG. 38, left panel). The use of a CpB binding protein is expected to reduce the conversion of cytosine to uracil but have no detectable effect on the conversion of 5mC to thymidine (FIG. 38, right panel). Thus, a wild-type cytidine deaminase or an altered cytidine deaminase can be used in one aspect of the present disclosure.

[00110] In another aspect, a CpG binding protein is used with methods that include a ten-eleven translocase (TET) protein. TET proteins catalyze the iterative conversion of 5mC to 5hmC, 5hmC to 5fC, and 5fC to 5caC (FIG. 39, left panel), however in some embodiments not all 5fC is converted to 5caC resulting in a mixture of 5fC and 5caC. The inventor found that unmethylated sequences, such as unmethylated cytosine, inhibits the conversion of methylated cytosine nucleotides, which leads to decreased beta values after sequencing (FIG. 39, left panel). The use of a CpG binding protein is expected to reduce the inhibition.

[00111] A. CpG binding proteins

[00112] Proteins that selectively bind to unmethylated CpG dinucleotides are referred to herein as "CpG binding proteins." In one embodiment, a CpG binding protein described herein is a naturally occurring protein, often referred to as a "CXXC domain." A CpG binding protein typically includes eight conserved cysteine residues which define two Zn ²⁺ -chelating units, and the residues TGHQ (SEQ ID NO:59) which are predicted to play a role in CpG binding (Ko et al., Nature. 2013 May 2; 497(7447): 122-126. doi: 10.1038/naturel2052). An example of a CpG binding protein is the IDAX CXXC domain (ICD) protein. ICD proteins useful in the methods disclosed herein include, but are not limited to, an ICD protein disclosed at Genbank accession number NP 079488.2 (SEQ ID NO:57, FIG. 37G). Other examples of CpG binding proteins include, but are not limited to, the amino acid sequence disclosed at Genbank accession number NP_057547 (SEQ ID NO:58, FIG.

37G), and the amino acid sequence disclosed at UNIPROT accession number K7EKZ6 (SEQ ID NO:63, FIG. 37L). CpG binding proteins useful in the methods disclosed herein are highly conserved. Other examples of CpG binding proteins include, but are not limited to, the CXXC domains of class I and class II (Xu et al., 2018, Structure 26, 85-95). An example of class I is the CXXC domain of protein CFP1, and examples of class II include the CXXC domain of proteins MLL1, MLL2, FBXL19, KDM2A, and KDM2B, and the 3rd CXXC domain of MBD1. Unless the context indicates otherwise, a "CpG binding protein" does not catalyze the oxidation of 5mC to 5hmC.

[00113] Whether a protein has the activity of preferably binding to unmethylated CpG dinucleotides can be determined by measuring the binding of a candidate CpG binding protein to a polynucleotide containing a single unmethylated CpG dinucleotide and comparing it to the binding of the candidate to polynucleotides containing methylated CpG or no CpG (TpG) (see, for instance, Ko et al., Nature. 2013 May 2; 497(7447): 122-126. doi: 10.1038/nature!2052). The activity of preferably binding to unmethylated CpG dinucleotides compared to methylated CpG dinucleotides is referred to herein as "CpG binding protein activity." A protein that preferably binds unmethylated CpG dinucleotides compared to methylated CpG dinucleotides has CpG binding protein activity.

[00114] The present disclosure also includes CpG binding proteins that have an amino acid sequence that is structurally similar to a reference CpG binding protein disclosed herein. In one embodiment, a reference CpG binding protein is one that includes the amino acid sequence of SEQ ID NO:57, or SEQ ID NO:58.

[00115] As used herein, a CpG binding protein may be "structurally similar" to a reference CpG binding protein if the amino acid sequence of the CpG binding protein possesses a specified amount of sequence similarity and/or sequence identity compared to the reference CpG binding protein.

[00116] Structural similarity of two amino acid sequences can be determined by aligning the residues of the two sequences (for example, a candidate CpG binding protein and a reference CpG binding protein described herein) to optimize the number of identical amino acids along the lengths of their sequences; gaps in either or both sequences are permitted in making the alignment in order to optimize the number of identical amino acids, although the amino acids in each sequence must nonetheless remain in their proper order. A candidate CpG binding protein is the CpG binding protein being compared to the reference CpG binding protein. A candidate CpG binding protein that has structural similarity with a reference CpG binding protein and CpG binding protein activity is included in the scope of the present disclosure.

[00117] Unless modified as otherwise described herein, a pair-wise comparison analysis of amino acid sequences can be conducted, for instance, by the local homology algorithm of Smith & Waterman, Adv. Appl. Math. 2:482 (1981), by the homology alignment algorithm of Needleman & Wunsch, J. Mol. Biol. 48:443 (1970), by the search for similarity method of Pearson & Lipman, Proc. Nat'l. Acad. Sci. USA 85:2444 (1988), by computerized implementations of these algorithms (GAP, BESTFIT, FASTA, and TFASTA in the Wisconsin Genetics Software Package, Genetics Computer Group, 575 Science Dr., Madison, Wis.), or by visual inspection (see generally Current Protocols in Molecular Biology, Ausubel et al., eds., Current Protocols, a joint venture between Greene Publishing Associates, Inc. and John Wiley & Sons, Inc., supplemented through 2004). One example of an algorithm that is suitable for determining structural similarity is the BLAST® algorithm, which is described in Altschul et al., J. Mol. Biol. 215:403-410 (1990). The BLAST® algorithm can be used to calculate percent sequence identity and percent sequence similarity between two sequences. Software for performing BLAST® analyses is publicly available through the National Center for Biotechnology Information.

[00118] In the comparison of two amino acid sequences, structural similarity may be referred to by percent "identity" or may be referred to by percent "similarity." "Identity" refers to the presence of identical amino acids. "Similarity" refers to the presence of not only identical amino acids but also the presence of conservative substitutions. Thus, in one embodiment the amino acid sequence of a cytidine deaminase protein having sequence similarity to a reference sequence may include conservative substitutions of amino acids present in that reference sequence.

[00119] A conservative substitution for an amino acid in a protein may be selected from other members of the class to which the amino acid belongs. For example, it is well-known in the art of protein biochemistry that an amino acid belonging to a grouping of amino acids having a particular size or characteristic (such as charge, hydrophobicity, or hydrophilicity) can be substituted for another amino acid without altering the activity of a protein, particularly in regions of the protein that are not directly associated with biological activity. For example, amino acids having a non-polar side chain include alanine, glycine, isoleucine, leucine, methionine, phenylalanine, proline, tryptophan, and valine; amino acids having a hydrophobic side chain include glycine, alanine, valine, leucine, isoleucine, proline, phenylalanine, methionine, and tryptophan; amino acids having a polar side chain include arginine, asparagine, aspartic acid, glutamine, glutamic acid, histidine, lysine, serine, cysteine, tyrosine, and threonine; and amino acids having an uncharged side chain include glycine, serine, cysteine, asparagine, glutamine, tyrosine, and threonine.

[00120] Thus, as used herein, reference to a CpG binding protein as described herein, such as reference to the amino acid sequence of one or more SEQ ID NOs described herein can include a protein with at least 80%, at least 85%, at least 86%, at least 87%, at least 88%, at least 89%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99% amino acid sequence similarity to the reference CpG binding protein.

[00121] Alternatively, as used herein, reference to a CpG binding protein as described herein, such as reference to the amino acid sequence of one or more SEQ ID NOs described herein can include a protein with at least 80%, at least 85%, at least 86%, at least 87%, at least 88%, at least 89%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99% amino acid sequence identity to the reference CpG binding protein.

[00122] A CpG binding protein can include additional amino acids as a fusion. In one embodiment, additional amino acids can have an activity that includes binding to single-stranded nucleic acids, such as single-stranded DNA. An example of a suitable amino acid sequence includes a single stranded binding protein. Examples of single stranded binding proteins include, but are not limited to E. colt SSB, T4 gp32, T7 gene 2.5 SSB, phage phi 29 SSB, RB69 bacteriophage gp32 protein. In another embodiment, additional amino acids can have an activity that includes increasing the solubility of a CpG binding protein. Examples of additional amino acids can increase the solubility of a CpG binding protein include, but are not limited to, a SUMO tag, a maltose binding protein (MBP) tag, or glutathione S- transferase (GST) tag. The additional amino acids fused to a CpG binding protein can present at the N-terminus or the C-terminus of a CpG binding protein, or at both the N- terminus and the C-terminus of a CpG binding protein. Reference to a CpG binding protein as described herein can include a CpG binding protein fusion.

[00123] B. Improvement of cytidine deaminase-based methods with CpG binding protein

[00124] One aspect of the present disclosure is methods, compositions, and kits related to the use of a CpG binding protein and a cytidine deaminase for detecting methylated nucleotides. The methods of using a CpG binding protein and a cytidine deaminase provided by the present disclosure can be easily integrated into essentially any application that includes sequencing library preparation, such as whole genome, accessible (e g., ATAC), and conformational state (e.g., HiC). It can be particularly useful in essentially any application using low input DNA or RNA such as, but not limited to, single cell combinatorial indexing (sci) methods like sci-WGS-seq, sci-MET-seq, and sci-ATAC-seq, sci-RNA-seq, and cell free DNA- based methods. Instead of using destructive methods for identifying methylated cytosines, integrating a cytidine deaminase into sequencing library production results in the more efficient enzyme-catalyzed conversion of modified cytosines during generation of target nucleic acids, thereby permitting more efficient detection of modified cytosines. A multitude of sequencing library methods are known to a skilled person that can be used in the construction of whole-genome or targeted libraries.

[00125] In general, methods for using a CpG binding protein and cytidine deaminase, such as an altered cytidine deaminase, include contacting target nucleic acids, e.g., DNA or RNA, with the CpG binding protein and the cytidine deaminase, under conditions suitable for conversion of modified cytidines, such as 5mC, to thymidine. The target DNA can be contacted with a CpG binding protein before adding a cytidine deaminase, or both CpG binding protein and cytidine deaminase can be added to the target DNA at the same time. Because amplification of DNA does not preserve the modification status of cytidine (e.g., the methylation status of 5mC is not retained), use of a CpG binding protein and a cytidine deaminase typically occurs before amplification of target DNA. Target nucleic acids can be contacted with a CpG binding protein and a cytidine deaminase at essentially any time in a method before an amplification, and the target nucleic acids can be single-stranded or double-stranded. In one embodiment, the target nucleic acids are single- stranded. For instance, target nucleic acids can be contacted with a CpG binding protein and a cytidine deaminase while the nucleic acids are inside a cell or a nucleus, after isolation of genomic or cell free DNA or mRNA, before or after fragmentation, or before or after tagmentation. The skilled person will recognize that target nucleic acids can be contacted with a CpG binding protein and a cytidine deaminase after addition of a universal sequence and/or an adapter, provided the universal sequence and/or an adapter is not added by amplification.

[00126] In some embodiments, a cytidine deaminase, such as an altered cytidine deaminase, as presented herein can be used to differentiate between 5-methyl cytosine (5mC) and 5- hydroxymethyl cytosine (5hmC). In such an embodiment, a sample of DNA suspected of including single-stranded DNA comprising at least one 5-methyl cytosine (5mC) or 5- hydroxymethyl cytosine (5hmC) is modified to prevent a cytidine deaminase from converting 5hmC to thymidine. Methods for blocking deaminase activity are known in the art, and any one of a number of methods can be used to protect 5hmC from deaminase activity. As one example, target DNA can be treated to modify 5hmC but not 5mC such that 5hmC is an unsuitable substrate for cytidine deaminase activity. In a specific example, a glucosyltransferase enzyme can be used to glucosylate 5hmC but not 5mC. Glucosyltransferase enzymes are known to those of skill in the art, and include, for example, p-glucosyltransferase (PGT or BGT). By way of example, the enzyme T4 p-glucosyltransferase is commercially available (PGT, NEB) and can be used for modification of 5hmC. Methods for using a PGT to glucosylate 5hmC are known in the art and can be used in conjunction with the use of cytidine deaminase enzymes as presented here. For example, a sample of DNA can be treated with a PGT to glucosylate 5hmC in the sample DNA prior to treating the DNA with the cytidine deaminase enzyme. By treating the sample DNA with a PGT, 5hmC is protected from the deaminase activity of the cytidine deaminase enzyme. Thus, 5mC will be detected in downstream readout, such as sequencing, PCR, array, and the like, as a thymidine. In contrast, any protected 5hmC sites will be detected as cytosine in the same readout. Enzymes, buffers, and conditions for performing glucosylation of 5hmC are known in the art, as exemplified by the methods disclosed in Schutsky et al., Nature biotechnology, 10.1038/nbt.4204. 8 Oct. 2018, doi: 10.1038/nbt.4204. A specific example of modification of 5hmC to protect it from cytidine deaminase activity is provided below in Example 9.

[00127] In one embodiment, the cytidine deaminase is a naturally occurring wild-type cytidine deaminase, such as a wild-type APOBEC3A. Wild-type cytidinedeaminases, such as APOBEC3A, deaminates cytosine (C), 5 methyl cytosine (5mC), and 5-hydroxymethyl cytosine (5hmC) efficiently in single-stranded DNA (FTG. 1A-1C). Treatment of DNA, such as genomic DNA, with wild type APOBEC3A results in the conversion of both C and 5mC to uracil and thymidine (FIG. 1A and FIG. IB), respectively, and reduces the complexity of the DNA to three bases for sequencing (FIG. ID. The use of a CpG binding protein with a wild-type cytidine deaminase binds to unmethylated cytosine and reduces the conversion of unmethylated cytosine to uracil, but has no detectable effect on the conversion of 5mC to thymidine. This allows one to distinguish between methylated and non-methylated cytosines, and has the advantage of maintaining the complexity of the DNA at four bases for sequencing.

[00128] In another embodiment, the cytidine deaminase is a non-naturally occurring altered cytidine deaminase. Altered cytidine deaminases described herein include mutated APOBEC proteins. For example, point mutations in human APOBEC3A proteins were produced in previous analyses and the ability of the mutant AP0BEC3A proteins to convert cytosine to uracil was determined. Modification of the tyrosine residue at position 130 to alanine (Y130A) consistently resulted in an APOBEC protein with no activity (see FIG. 6c of Bulliard et al., 2011, J Virol., 85(4): 1765-1776, and FIG. 5a of Shi et al., 2017, Nat Struct Mol Biol., 24(2): 131-139). Proceeding contrary to Bulliard and Shi, the inventors made the surprising and unexpected discovery that certain mutations at position 130 of APOBEC3A alter the enzyme's rate of deamination on 5mC compared to C substrates.

[00129] As described herein, the cognate tyrosine (Y) at position 130 was individually mutated to all possible canonical amino acid substitutions, including A, C, D, E, F, G, H, I, K, L, M, N, P, Q, R, S, T, V, and W, and evaluated for activity on C, 5mC, and 5hmC substrates. For example, an APOBEC3A mutant containing a tyrosine to alanine point mutation in position 130 (Y130A) was found to preferentially deaminate 5mC instead of C (5mC was converted to T at a greater rate than C was converted to U), and an APOBEC3A mutant containing a tyrosine to leucine point mutation in position 130 (Y130L) was found to preferentially deaminate C instead of 5mC (C was converted to U at a greater rate than 5mC was converted to T). The deamination of 5mC to T leads to C to T mutations which can be identified by standard sequencing methods. As a result, in one embodiment the treatment of DNA with an altered cytidine deaminase of the present disclosure preferentially converts 5mC to thymidine (FIG. IE). Analysis of the sample DNA after treatment with the modified cytidine deaminase described herein, for example, by sequencing of the sample DNA, and optional comparison to a reference (e.g., reference sequence) permits easy identification of C to T point mutations, and these point mutations are inferred as 5mC positions.

[00130] In another example, an APOBEC3A mutant containing a tyrosine to tryptophan point mutation in position 130 (Y130W) maintained the ability to deaminate C and 5mC to U and T, respectively, but lost the ability to deaminate 5hmC, 5fC, and 5caC. When the sequence of the nucleic acid exposed to this APOBEC3A mutant is determined, C and 5mC are deaminated to U and T, respectively, and are read out as T by a sequencer. On the other hand, 5hmC is not deaminated by the APOBEC3A mutant and is read out as a C. (FIG. IF). 5fC and 5caC are also not deaminated by this AP0BEC3A mutant, and therefore cannot be distinguished from 5hmC. However, the abundance of 5fC and 5caC is orders of magnitude lower than 5hmC in human genomic DNA, approaching the detection limit of mass spectrometers used for such measurements (Ito et al., 2011, Science 333, 1300-1303 (2011); Wagner et al., 2015, Angew. Chem. Int. Edn Engl. 54, 12511-12514 (2015); Bachman et al., 2015, Nat. Chem. Biol. 11, 555-557 (2015)). Therefore, signals from 5fC and 5caC should be insignificant compared to 5hmC.

[00131] Sequencing of the sample DNA after treatment with a CpG binding protein and altered cytidine deaminase and optional comparison to a reference sequence permits easy identification of C to T point mutations, and these point mutations are inferred as 5mC positions.

[00132] B.l. Cytidine deaminases

[00133] The compositions, kits, and methods may include an altered cytidine deaminase. The present disclosure provides three types of altered cytidine deaminases. One type of altered cytidine deaminase preferentially deaminates 5mC instead of C (i.e., converts 5mC to T at a greater rate than converting C to U) and is referred to herein as having “cytosinedefective deaminase activity.” A second type of altered cytidine deaminase preferentially deaminates C instead of 5mC (i.e., converts C to U at a greater rate than converting 5mC to T) and is referred to herein as having “5mC-defective deaminase activity.” A third type of altered cytidine deaminase preferentially deaminates C and 5mC to U and T, respectively, and has significantly reduced deamination of 5hmC, 5fC, and 5caC. The third type is referred to herein as having “5hmC-defective deaminase activity.” Unless the context indicates otherwise, reference to an altered cytidine deaminase includes altered cytidine deaminases having cytosine-defective deaminase activity, altered cytidine deaminases having 5mC-defective deaminase activity, and altered cytidine deaminases having 5mC- defective deaminase activity.

[00134] Cytidine deaminases described herein include apolipoprotein B mRNA editing enzymes, catalytic polypeptide-like (APOBEC) and activation induced cytidine deaminase (AID). Wild-type APOBEC and AID cytidine deaminases have the activity of deaminating cytidine (C) of DNA and/or RNA to form uridine (U). An altered cytidine deaminase of the present disclosure has an altered rate of deamination of 5mC, 5mC, and/or 5hmC verses C substrates when compared to the wild-type enzyme. Unless the context indicates otherwise, reference to a "cytidine deaminase" includes both wild-type cytidine deaminases, altered cytidine deaminases, and wild-type cytidine deaminases and altered cytidine deaminases structurally similar to a reference cytidine deaminase. Structural similarity is described herein. An "altered cytidine deaminase" and "recombinant cytidine deaminase" refers to any of the altered cytitidine deaminases described herein that comprise one or more changes from the reference (i.e., wildtype) amino acid sequence that provide the unexpected property of an altered deamination profde, e.g., alters its ability to preferentially deaminate one form of cytosine over another.

[00135] Whether a protein has cytidine deaminase activity may be determined by in vitro assays. One example of an in vitro assay is based on digestion with the restriction enzyme Ai l (see Example 1). A protein that can deaminate 5mC to thymidine has cytidine deaminase activity.

[00136] An altered cytidine deaminase that preferentially deaminates 5mC instead of C (i.e., has cytosine-defective deaminase activity) can have a catalytic efficiency that is at least 10- fold, at least 50-fold, or at least 100-fold higher on 5mC than C substrates. In one embodiment, an altered cytidine deaminase that preferentially deaminates 5mC instead of C can have a catalytic efficiency that is no greater than 1500-fold higher on 5mC than C substrates.

[00137] An altered cytidine deaminase that preferentially deaminates C instead of 5mC (i.e., has 5mC -defective deaminase activity) can have a catalytic efficiency that is at least 10-fold, at least 50-fold, or at least 100-fold higher on C than 5mC substrates. In one embodiment, an altered cytidine deaminase that preferentially deaminates C instead of 5mC can have a catalytic efficiency that is no greater than 1500-fold higher on C than 5mC substrates.

[00138] When compared to a wild type cytidine deaminase, an altered cytidine deaminase that deaminates C and 5mC to U and T, respectively, and has significantly reduced deamination of 5hmC (i.e., has 5hmC-defective deaminase activity), the deamination of 5hmC by an altered cytidine deaminase disclosed herein is reduced by at least 80%, at least 90%, or at least 99% compared to the wild type cytidine deaminase. In one embodiment, the deamination of 5hmC by an altered cytidine deaminase disclosed herein is undetectable using an assay such as the .S'liz/I-based assay described herein.

[00139] In certain embodiments, an altered cytidine deaminase of the present disclosure is based on a wild-type cytidine deaminase, e.g., a member of the APOBEC protein family. An altered cytidine deaminase of the present disclosure that is "based on" a member of the APOBEC protein family means the altered cytidine deaminase is an APOBEC protein that includes one or more of the substitution mutations described herein as compared to a reference APOBEC sequence. An altered cytidine deaminase of the present disclosure that is "based on" a member of the APOBEC protein family can also include conservative and/or nonconservative mutations as described herein.

[00140] The APOBEC protein family includes subfamilies AID, APOBEC 1, APOBEC2, APOBEC3 (including 3A, 3B, 3C, 3D, 3F, 3G, 3H), and APOBEC4. A wild-type cytidine deaminase can be a member of the APOBEC family of proteins or have structural similarity to a member of the APOBEC family of proteins. An altered cytidine deaminase of the present disclosure can be based on a member of the AID subfamily, the APOBEC 1 subfamily, the APOBEC2 subfamily, the APOBEC3 subfamily (e.g., the 3A subfamily, the 3B subfamily, the 3C subfamily, the 3D subfamily, the 3F subfamily, the 3G subfamily, or the 3H subfamily), or the APOBEC4 subfamily. A member of the APOBEC protein family can be from a vertebrate, such as a mammal. Examples of mammals include, but are not limited to, rodents, primates, rabbit, bovine (e.g., cow), porcine (e.g., pig), and equine (e.g., horse). An example of a primate is a human and a chimpanzee.

[00141] The APOBEC protein family is a member of the large cytidine deaminase superfamily that contains a canonical zinc-dependent deaminase (ZDD) signature motif embedded within a core cytidine deaminase fold. This fold includes a five-stranded mixed beta (b)-sheet surrounded by six alpha (a)-helices with the order al-bl-b2-a2-b3-a3-b4-a4-b5-a5-a6 (Salter et al., Trends Biochem Sci. 201641(7):578-594. doi: 10.1016/j tibs.2016.05.001 ; Salter et al., Trends Biochem. Sci. 2018, 43(8):606-622

31 doi.org/10.1016/j tibs.2018.04.013). Each cytidine deaminase domain core structure of APOB EC proteins contains a highly conserved spatial arrangement of the catalytic center residues of a zinc-binding motif H- [P/A/V]-E-Xp3-28|-P-C-X[2-4]-C (SEQ ID NO: 12) (referred to herein as the ZDD motif, where X is any amino acid, and the subscript range of numbers after X refers to the number of amino acids) (Salter et al., Trends Biochem Sci. 2016 41 (7): 578— 594. doi : 10.1016/j .tibs.2016.05.001). Without intending to be limited by theory, the H and two C residues coordinate a Zn atom, and the E residue polarizes a water molecule near the Zn-atom for catalysis (Chen et al., 2021, Viruses, 13:497, doi.org/ 10.3390/v 13030497).

[00142] Some members of the APOB EC protein family, e.g., the AID subfamily, the APOB EC 1 subfamily, the APOBEC2 subfamily, the APOBEC3A subfamily, the APOBEC3C subfamily, the APOBEC3H subfamily, and the APOBEC4 subfamily, include one copy of the ZDD motif. Other members of the APOBEC protein family, e.g., the APOBEC3B subfamily, the APOBEC3D subfamily, the APOBEC3F subfamily, and the APOBEC3G subfamily, include two copies of the ZDD motif, but often only the C-terminal copy is active (Salter et al., Trends Biochem Sci. 2016 41 (7):578— 594. doi: 10.1016/j .tibs.2016.05.001). Thus, a cytidine deaminase disclosed herein includes one or two ZDD motifs. In one embodiment, a wild-type cytidine deaminase that is a member of the APOBEC3 A subfamily or an altered cytidine deaminase based on a member of the APOBEC3A subfamily includes the following ZDD motif: HXEX24SW(S/T)PCX[2- 4]CX6FX8LX5R(L/I)YX[8-ii]LX2LX[io]M (SEQ ID NO: 13) (where X is any amino acid, and the subscript number or range of numbers after X refers to the number of amino acids) (Salter et al., Trends Biochem Sci. 201641(7)578-594. doi: 10.1016/j .tibs.2016.05.001).

[00143] In one embodiment, a wild-type cytidine deaminase that is a member of the following subfamilies or an altered cytidine deaminase based on a member of the following subfamilies APOBEC3A, APOBEC3B, APOBEC3C, APOBEC3D, APOBEC3F, and APOBEC 3 G, and can include one or more highly conserved sites that are part of the active site and within the ZDD motif SEQ ID NO: 12. The sites include tryptophan at position 98 and serine or threonine at position 99 (Kouno et al., 2017, Nat. Comm, 8:15024, DOI: 10.1038/ncommsl5024). [00144] In addition to the ZDD motif, a member of the APOBEC protein family also includes other highly conserved residues that are part of the active site but not present as part of the ZDD motif SEQ ID NO: 12. A member ofthe APOBEC3A subfamily, APOBEC3B subfamily, APOBEC3C subfamily, APOBEC3D subfamily, APOBEC3F subfamily, and APOBEC3G subfamily typically includes one or more of the following highly conserved sites that are part of the active site: arginine at position 28; histidine, asparagine, or arginine at position 29; serine or threonine, preferably threonine, at position 31; asparagine or aspartic acid at position 57; tyrosine or phenylalanine at position 130; asparagine or tyrosine at position 131; asparagine, tyrosine, or phenylalanine, preferably tyrosine, at position 132; and arginine or lysine at position 189 (Kouno et al., 2017, Nat. Comm, 8: 15024, DOI: 10.1038/ncommsl5024).

[00145] An altered cytidine deaminase of the present disclosure includes a substitution mutation at one or more residues when compared to a reference cytidine deaminase. A substitution mutation can be at the same position or a functionally equivalent position compared to the reference cytidine deaminase. Reference cytidine deaminases and functionally equivalent positions are described in detail herein. The skilled person will readily appreciate that an altered cytidine deaminase described herein is not naturally occurring.

[00146] A reference cytidine deaminase can be a member of the APOBEC protein family. Essentially any known member of the APOBEC protein family can be a reference cytidine deaminase. The skilled person can easily identify members of each of the subfamilies by using a publicly available database such as the Protein database available at the National Center for Biotechnology Information (ncbi.nlm.nih.gov/protein) and searching for APOBEC 1, APOBEC2, APOBEC3A, APOBEC3B, APOBEC3C, APOBEC3D, APOBEC3F, APOBEC3G, APOBEC3H, APOBEC4, or, when identifying members of the AID family, Activation-induced cytidine deaminase. A wild-type reference cytidine deaminase has the activity of binding single-stranded DNA (ssDNA) and deaminating a cytosine present on the ssDNA to convert it to uracil. In one embodiment, a wild-type reference cytidine deaminase has the activity of binding single-stranded RNA (ssRNA) and deaminating a cytosine present on the ssRNA to convert it to uracil. Methods for determining whether a protein binds ssDNA or ssRNA and deaminates a cytosine present are known to the skilled person.

[00147] In one embodiment, a deaminase has an amino acid sequence that is based on a reference sequence that is a member of the APOBEC protein family, includes a ZDD motif H- [P/A/V]-E-X[23-28]-P-C-X[2-4]-C (SEQ ID NO: 12), and an altered cytidine deaminase further includes at least one substitution mutation disclosed herein. Optionally, a cytidine deaminase includes other active site residues disclosed herein. Non-limiting examples of reference cytidine deaminase proteins, which can be used as reference cytidine deaminase proteins, are shown in the following table.

[00148] Table 1 . Examples of members of the APOBEC protein subfamilies. UniProt, database of protein sequence and functional information, available at uniprot.org; GenBank, collection of nucleotide sequences and their protein translations, available at ncbi.nlm.nih.gov/protein/.

[00149] In one embodiment, a wild type or altered cytidine deaminase has an amino acid sequence that is based on a reference sequence that is a member of the AP0BEC3A subfamily, and includes a ZDD motif HXEX24SW(S/T)PCX[2-4]CX6FX8LX5R(L/I)YX[ ₈ -ii]LX2LX[io]M (SEQ ID NO: 13) (where X is any amino acid, and the subscript number or range of numbers after X refers to the number of amino acids), and an altered cytidine deaminase further includes at least one substitution mutation disclosed herein. In one embodiment, the substitution mutation is a substitution mutation at the underlined tyrosine, such as a substitution mutation to alanine (A). Optionally, the altered cytidine deaminase includes other active site residues disclosed herein. In one embodiment, the amino acid sequence of a wild type or altered cytidine deaminase includes the amino acids of a member of the AP0BEC3A subfamily: X[i6-26]-GRXXTXLCYXV-Xi ₅ -GXXXN-Xi2-HAEXXF-Xi4- YXXTWXXSWSPC- X[2-4]-CA-X5-FL-X7-LXIXXXR(L/I)Y-X8-GLXXLXXXG-X5-M- X4-FXXCWXXFV-X6-FXPW-X13-LXXI- X _L2 -6] (SEQ ID NO: 14) (where X is any amino acid, and the subscript number or range of numbers after X refers to the number of amino acids), or a subset thereof, and an altered cytidine deaminase further includes at least one substitution mutation disclosed herein. In one embodiment, the substitution mutation is a substitution mutation at the underlined tyrosine, such as a substitution mutation to alanine (A) or to tryptophan (W).

[00150] In one embodiment, the amino acid sequence of a wild type or altered cytidine deaminase includes a subset of the amino acids of a member of the APOBEC3A subfamily: X26- GRXXTXLCYXV-X15-G-X16-HAEXXF-X14-YXXTWXXSWSPC-X4-CA-X5-FL-X7 - LXIFXXR(L/I)Y-X8-GLXXLXXXG-X ₅ -M-X4-FXXCWXXFV-X6-FXPW-X13-LXXI-X6 (SEQ ID NO: 15) (where X is any amino acid, and the subscript number after X refers to the number of amino acids present), or a subset thereof, and an altered cytidine deaminase further includes at least one substitution mutation disclosed herein. In one embodiment, the substitution mutation is a substitution mutation at the underlined tyrosine (Y), such as a substitution mutation to alanine (A) or to tryptophan (W). [00151] A substitution mutation can be at the same position or a functionally equivalent position compared to a reference cytidine deaminase. By "functionally equivalent" it is meant that the altered cytidine deaminase has the amino acid substitution at the amino acid position in a reference cytidine deaminase that has the same functional role in both the reference cytidine deaminase and the altered cytidine deaminase.

[00152] In general, functionally equivalent substitution mutations in two or more different cytidine deaminases occur at homologous amino acid positions in the amino acid sequences of the cytidine deaminases. Hence, use herein of the term "functionally equivalent" also encompasses mutations that are "positionally equivalent" or "homologous" to a given mutation, regardless of whether or not the particular function of the mutated amino acid is known. It is possible to identify the locations of functionally equivalent and positionally equivalent amino acid residues in the amino acid sequences of two or more different cytidine deaminases on the basis of sequence alignment and/or molecular modelling. An example of a sequence alignment to identify positionally equivalent and/or functionally equivalent residues is set forth in FIG. 3. For example, the residues in the members of the APOBEC3A subfamily in FIG. 3 that are vertically aligned are considered positionally equivalent as well as functionally equivalent to the corresponding residue in the human APOBEC3A amino acid sequence. Thus, for example, as shown in FIG. 3, the tyrosine at residue 130 of the APOBEC3 A proteins of Homo sapiens, Pongo pygmaeus, Nomascus leucogenys, Pan troglodytes, and Gorilla gorilla and the tyrosine at residue 133 of the APOBEC3A protein from Macaca fascicularis are functionally equivalent and positionally equivalent. The skilled person can easily identify functionally equivalent residues in cytidine deaminases.

[00153] In one embodiment, a cytidine deaminase, such as an altered cytidine deaminase has an amino acid sequence that is structurally similar to a reference cytidine deaminase disclosed herein. In one embodiment, a reference cytidine deaminase is one that includes the amino acid sequence of a sequence listed in Table 1, SEQ ID NO: 14, or SEQ ID NO:15.

[00154] As used herein, a cytidine deaminase, such as an altered cytidine deaminase may be "structurally similar" to a reference cytidine deaminase if the amino acid sequence of the altered cytidine deaminase possesses a specified amount of sequence similarity and/or sequence identity compared to the reference cytidine deaminase.

[00155] Structural similarity of two amino acid sequences can be determined by aligning the residues of the two sequences (for example, a candidate altered cytidine deaminase and a reference cytidine deaminase described herein) to optimize the number of identical amino acids along the lengths of their sequences; gaps in either or both sequences are permitted in making the alignment in order to optimize the number of identical amino acids, although the amino acids in each sequence must nonetheless remain in their proper order. A candidate altered cytidine deaminase is the cytidine deaminase being compared to the reference cytidine deaminase. A candidate altered cytidine deaminase that has structural similarity with a reference cytidine deaminase and cytidine deaminase activity is an altered cytidine deaminase.

[00156] Structural similarity of two amino acid sequences can be determined as described herein. A candidate cytidine deaminase (for example, an altered cytidine deaminase) that has structural similarity with a reference cytidine deaminase and cytidine deaminase activity is a cytidine deaminase. In one embodiment the amino acid sequence of a cytidine deaminase protein having sequence similarity to a reference sequence may include conservative substitutions of amino acids present in that reference sequence. Conservative substitutions are described herein.

[00157] Thus, as used herein, reference to a cytidine deaminase as described herein, such as reference to the amino acid sequence of one or more SEQ ID NOs of a cytidine deaminase described herein, such as SEQ ID NOs: 1-15 and 19-23, can include a protein with at least at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 86%, at least 87%, at least 88%, at least 89%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99% amino acid sequence similarity to the reference cytidine deaminase. Examples of altered cytidine deaminases having similarity with a reference amino acid sequence includes those having, for instance, at least 80%, at least 85%, at least 90%, or at least 95% similarity with SEQ ID NO: 16 (FIG. 37A) and having an alanine at amino acid 130. Other examples of altered cytidine deaminases having similarity with a reference amino acid sequence includes those having, for instance, at least 80%, at least 85%, at least 90%, or at least 95% similarity or identity with SEQ ID NO: 17 (FIG. 37A) and having an alanine at amino acid 130 and a histidine at amino acid 132.

[00158] Alternatively, as used herein, reference to a cytidine deaminase as described herein, such as reference to the amino acid sequence of one or more SEQ ID NOs described herein, such as SEQ ID NOs: 1-15 and 19-23, can include a protein with at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 86%, at least 87%, at least 88%, at least 89%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99% amino acid sequence identity to the reference cytidine deaminase. -Examples of altered cytidine deaminases having identity with a reference amino acid sequence includes those having, for instance, at least 80%, at least 85%, at least 90%, or at least 95% similarity with SEQ ID NO: 16 (FIG. 37A) and having an alanine (A) at amino acid 130. Other examples of altered cytidine deaminases having identity with a reference amino acid sequence includes those having, for instance, at least 80%, at least 85%, at least 90%, or at least 95% similarity or identity with SEQ ID NO: 17 (FIG. 37A) and having an alanine (A)at amino acid 130 and a histidine (H) at amino acid 132.

[00159] B.2. Substitution mutations in a cytidine deaminase

[00160] An altered cytidine deaminase of the present disclosure includes a substitution mutation at a position functionally equivalent to tyrosine at position 130 (Y130) in a member of the APOBEC3A subfamily (for instance, SEQ ID NO:3). Accordingly, an alignment can be produced using a member of the APOBEC3A subfamily (for instance, SEQ ID NO:3) and another candidate cytidine deaminase from the APOBEC3A subfamily or a different APOBEC subfamily. In one embodiment, the candidate is selected from APOPEC subfamilies APOBEC 1 or AID. An example of an algorithm that can be used to produce an alignment is Clustal O. In some APOBEC family proteins, the wild-type residue at a position functionally equivalent to Y130 is phenylalanine (F). [00161] In another embodiment, an altered cytidine deaminase of the present disclosure includes a substitution mutation at a position functionally equivalent to the tyrosine (Y) of ZDD motif HXEX24SW(S/T)PCX[2.4]CX6FX8LX5R(L/I)YX[8-ii]LX2LX[io]M (SEQ ID NO: 13) in a member of the APOBEC family, such as a member of the APOBEC3A subfamily. The underlined tyrosine (Y) of SEQ ID NO: 13 is the position functionally equivalent to the tyrosine amino acid 130 of the APOBEC3A protein SEQ ID NO:3.

[00162] In one embodiment, the substitution mutation at a position functionally equivalent to Y130 increases cytidine deaminase activity and preferentially acts on 5mC compared to cytosine (i.e., has cytosine-defective deaminase activity). The substitution mutation can be a mutation to alanine (A), glycine (G), phenylalanine (P), histidine (H), glutamine (Q), methionine (M), asparagine (N), lysine (K), valine (V), aspartic acid (D), glutamic acid (E), serine (S), cysteine (C), proline (P), or threonine (T). For example, the altered cytidine deaminase can comprise SEQ ID NO:64 (FIG. 37L), wherein X is selected from A, G, F, H, Q, M, N, K, V, D, E, S, C, P or T (and is not Y), or can comprise SEQ ID NO:65 (FIG. 37L), wherein Z is selected from A, G, F, H, Q, M, N, K, V, D, E, S, C, P or T (and is not Y), preferably, in one embodiment, X or Z is A or L. In an exemplary aspect of this embodiment, the substitution mutation at a position functionally equivalent to Y130 is a mutation to alanine (A) (e.g., SEQ ID NO: 16). Specific examples of altered cytidine deaminases having increased activity and preferentially acting on 5mC compared to cytosine include SEQ ID NO: 16 or a sequence having at least 90%, at least 95%, at least 98%, at least 99% sequence identity to SEQ ID NO: 16 (FIG. 37A) and comprising Y130A.

[00163] An altered cytidine deaminase of the present disclosure having cytosine-defective deaminase activity (i.e., converts 5mC to T at a greater rate than converting C to U) optionally includes a second substitution mutation at a position two, three, four, or five amino acids on the C-terminal side of the Y130 position, or functionally equivalent to the Y130 position. In one embodiment, the second mutation is a tyrosine (Y), tryptophan (W), cysteine (C), histidine (H), or phenylalanine (F) at a position two, three, four, or five amino acids on the C-terminal side of the Y130 position, or functionally equivalent to the Y130 position. In one embodiment, the second mutation is at a position functionally equivalent to tyrosine at position 132 (Y132) in a member of the APOBEC3A subfamily (for instance, SEQ ID N0:3). An APOBEC protein, such as an APOBEC3A protein, containing substitution mutations at both the first site, at a position functionally equivalent to Y130, and the second site, at a position two, three, four, or five amino acids on the C-terminal side of the Y130 position, increases the preferential activity to act on 5mC compared to the same APOBEC protein, such as an AP0BEC3A protein, containing one substitution mutation at Y130. In one embodiment, the substitution mutation at the second position is an amino acid having a positively charged side chain and selected from arginine (R), histidine (H), lysine (K), or a polar side chain selected from glutamine (Q). In one embodiment, the substitution mutation at the second position is histidine (H), such as Y132 to histidine. The double mutant containing both first and second mutations can be any substitution mutation at a position functionally equivalent to Y130 described herein and any second substitution mutation at a position two, three, four, or five amino acids on the C-terminal side of the Y130 position described herein, in any combination. . For example, the altered cytidine deaminase can be SEQ ID NO: 3, 15, or 16 and have a substitution at Y130 and Y132, or the position functionally equivalent to Y130 and Y132 as described herein. One example of an altered cytidine deaminase is SEQ ID NO:66 (FIG. 37M) comprising Y130X and Y132Z, where X is selected from (A), (L), or (W) (preferably (A)), and Z is selected from (R), (H), (L), or (Q), preferably (H). This encompasses examples including, but not limited to, for example Y130A and Y132R, Y130A and Y132H, Y130A and Y132L, Y130A and Y132Q, Y130L and Y132R, Y130L and Y132H, Y130L and Y132L, Y130L and Y132Q, Y130W and Y132R, Y130W and Y132H, Y130W and Y132L, Y130W and Y130Q, or any suitable combinations therein. In one embodiment, the double mutant includes substitution mutations Y130A and Y132H. Specific examples of altered cytidine deaminases having both substitution mutations and preferentially acting on 5mC compared to the APOBEC protein having just the single substitution mutation at cytosine include SEQ ID NO: 17 or a sequence having at least 90%, at least 95%, at least 98%, at least 99% sequence identity to SEQ ID NO: 17 and comprising Y130A and Y132H.

[00164] The person of ordinary skill in the art can confirm the preferential 5mC deaminase activity of the arginine, glutamine, histidine and lysine substitution mutations at the second position in the double mutants described above. For example, double mutants can be constructed to create an altered cytidine deaminase having a first substitution mutation at a position functionally equivalent to Y130 and a second arginine, glutamine, histidine or lysine substitution mutation at the tyrosine position two amino acids on the C-terminal side of the Y130 position, and then evaluated for deamination of C residues in one assay and deamination of 5mC residues in a second assay. Using an assay such as the .S'liv/I-based assay described herein, the ratio of 5mC deamination and C deamination can be compared to identify those double mutants that preferentially deaminate 5mC compared to C. One of ordinary skill in the art could similarly test double mutants having a tyrosine at a position three, four or five positions C-terminal to the position functionally equivalent to Y130 and confirm that a substitution mutation at that position to arginine, glutamine, histidine or lysine, in combination with a mutation at the position functionally equivalent to Y130 (such as Y130A), as double mutants that preferentially deaminate 5mC compared to C.

[00165] Some embodiments presented herein relate to substitution mutations that result in 5mC- defective deaminase activity (i.e., converts C to U at a greater rate than converting 5mC to T). In one embodiment, the substitution mutation at a position functionally equivalent to Y130 increases cytidine deaminase activity and preferentially acts on cytosine compared to 5mC and is a mutation to an amino acid having a non-polar side chain or a hydrophobic side chain, such as leucine (L) or tryptophan (W). In an exemplary aspect of this embodiment, the substitution mutation at a position functionally equivalent to Y130 is a mutation to leucine. Other examples of mutations that result in increased preferential deamination activity on cytosine compared to 5mC include a single mutant with Y132P, and double mutants with a substitution mutation at Y130V and Y132H, or Y130W and Y132H. Specific examples of altered cytidine deaminases having increased cytidine deaminase activity and preferentially acts on cytosine compared to 5mC include SEQ ID NO: 18 or a sequence having at least 90%, at least 95%, at least 98%, at least 99% sequence identity to SEQ ID NO: 18 and comprising Y130L.

[00166] In one embodiment, the substitution mutation is at a position functionally equivalent to Y130 that results in 5hmC-defective deaminase activity (i.e., preferentially deaminates C and 5mC to U and T, respectively, and has significantly reduced deamination of 5hmC). In an exemplary aspect of this embodiment, the substitution mutation at a position functionally equivalent to Y130 is a mutation to an amino acid having a non-polar side chain or a hydrophobic side chain, such as tryptophan (W). Specific examples of altered cytidine deaminases having the ability to deaminate C and 5mC to U and T, respectively, but reduced ability to deaminate 5hmC, preferably no detectable ability to deaminate 5hmC include SEQ ID NO:92 where Y130 is W (FIG. 37N) or a sequence having at least 90%, at least 95%, at least 98%, at least 99% sequence identity to SEQ ID NO:92 and comprising Y130W.

[00167] A cytidine deaminase, such as an altered cytidine deaminase described herein can include additional mutations. Typically, additional mutations do not unduly alter the activity of the altered cytidine deaminase. One or more additional mutations can be a conservative mutation, a non-conservative mutation, or a combination of both conservative mutations and non-conservative mutations. A cytidine deaminase, such as an altered cytidine deaminase described herein can be a truncated protein. A truncated protein is a fragment of a cytidine deaminase of the present disclosure that retains the ability to deaminate 5mC to thymidine. A truncated cytidine deaminase can include a deletion of 1 to 13 amino acids on the N-terminal end of the protein, a deletion of 1 to 3 amino acids on the C-terminal end of the protein, or a combination thereof

[00168] An altered cytidine deaminase of the present application may be an altered cytidine deaminases described in International Patent Application No. PCT/US2023/017846, entitled “ALTERED CYTIDINE DEAMINASES AND METHODS OF USE,” filed on April 7, 2023. Other altered cytosine deaminases can be used in the practice of this invention, including, but not limited to, for example, the hyperactive altered forms of APOBEC, for instance, those described in U.S. Patent No. 10,961,525 entitled “Hyperactive AID/ APOBEC and HMC Dominant TET enzymes”, the content of the altered APOBECs incorporated herein by reference in their entirety (for example, but not limited to chimeric deaminases comprising a first domain of APOBEC 4B (A3B) and the catalytic domain of APOBECC3A (A3 A), A3B protein having amino acid mutations at D196H, T197I, delta(206-210), Ins(206)GIG, R212H, Q213K, W228S, I230K, M235R, C239H, E241K, E342K, Y343H, Y350D, R351H, and E363D). Further, other altered APOBEC proteins are described in U.S. Patent No. 9,896,726 entitled “Methods and compositions for discrimination between cytosine and modifications thereof, and for methylome analysis”, the altered cytosine deaminases incorporated by reference herein. The listing of wild-type and altered cytosine deaminases that can be used in the present disclosure herein is not exhaustive and it is contemplated others may be used and contemplated with this disclosure.

[00169] Reaction conditions suitable for conversion of modified cytidines, such as conversion of 5mC to thymidine, by a wild type or an altered cytidine deaminase described herein include, but are not limited to, a substrate of target nucleic acid that is single-stranded (ss) DNA or RNA suspected of including at least one modified cytidine, buffer, pH, temperature of the reaction, time of the reaction, and concentration of the cytidine deaminase and/or ss DNA or RNA substrate. In one embodiment, double-stranded (ds) DNA can be denatured and exposed to a wild type or altered cytidine deaminase. Methods for denaturing dsDNA are known and routine, and include heat treatment or chemical treatment, such as NaOH, formamide, dimethyl sulfoxide (DMSO), N,N- dimethylformamide (DMF), or a combination thereof.

[00170J Wild-type cytidine deaminases typically function at near-neutral pH, e.g., pH 7. Altered cytidine deaminases described herein can have increased activity below neutral pH. In some embodiments, the pH of a reaction that includes an altered cytidine deaminase described herein can be no greater than pH 6.9, no greater than pH 6.7, no greater than pH 6.5, no greater than pH 6.3, no greater than pH 6.1, no greater than pH 6.0, no greater than pH 5.9, no greater than pH 5.7, no greater than pH 5.5, no greater than pH 5.3, no greater than pH 5.2, or no greater than pH 5.1. In some embodiments, the pH of a reaction that includes an altered cytidine deaminase described herein can be at least pH 5.1, at least pH 5.3, at least pH 5.5, at least pH 5.7, at least pH 5.9, at least pH 6.1, at least pH 6.3, at least pH 6.5, at least pH 6.7, or at least pH 6.9. In some embodiments, the pH of a reaction that includes an altered cytidine deaminase described herein can be no greater than pH 7.5, no greater than pH 7.3, or no greater than pH 7.1. Examples of ranges of pH in a reaction include at least 5.1 to no greater than 6.9, at least 5.1 to no greater than 6.5, at least 5.1 to no greater than 6.3, or at least 5.1 to no greater than 6.1. The activity of an altered cytidine deaminase to deaminate a 5mC oligonucleotide substrate can be an increased catalytic activity that is at least 10-fold greater, at least 50-fold greater, or at least 100-fold greater when comparing activity (see Example 3).

[00171] It is expected that a cytidine deaminase, such as an altered cytidine deaminase can function in essentially any buffer. Examples of useful buffers include, but are not limited to: a citrate buffer, such as the citrate buffer available from Thermo Fisher Scientific (Cat. No. #005000); sodium acetate buffer, Bis Tris-Propane HC1; and Tris-HCl Tris. Examples of other buffers include, but are not limited to, Bicine, DIPSO, glycylglycine, HEPES, imidazole, malonate, MES, MOPS, PB, phosphate, PIPES, SPG, succinate, TAPS, TAPSO, trincine. In some embodiments a reducing agent such as dithiothreitol (DTT) can be present. In some embodiments a divalent cation is not included

[00172] A deamination reaction can occur at a temperature of 25°C to 37°C, such as 37°C. Some altered cytidine deaminases described herein preferentially deaminate a modified cytosine to thymidine at a faster rate than deamination of cytosine to uracil. Thus, in some embodiments the time of reaction can be used to maximize the difference of deamination of modified cytosine versus deamination of cytosine. In one embodiment, the reaction can proceed for at least 15 minutes, at least 30 minutes, at least 45 minutes, at least 60 minutes, at least 90 minutes, at least 120 minutes, or at least 150 minutes, and for no greater than 15 minutes, no greater than 30 minutes, no greater than 45 minutes, no greater than 60 minutes, no greater than 90 minutes, no greater than 120 minutes, no greater than 150 minutes, or no greater than 180 minutes. In one embodiment, a deamination reaction can include an altered cytidine deaminase at a concentration from at least 0.5 micromolar (pM) to no greater than 5 pM. For instance, the concentration of the enzyme can be at least 0.5, at least 1 pM, at least 2 pM, at least 3 pM, at least 4 pM, or 5 pM, and/or no greater than 5 pM, no greater than 4 pM, no greater than 3 pM, no greater than 2 pM, no greater than 1 pM, or 0.5 pM. In one embodiment, a deamination reaction can include nucleic acids at a concentration of at least 1 picomolar (pM) to no greater than 2 pM. For instance, the concentration of nucleic acids can be at least 1 pM, at least 3 pM, at least 6 pM, at least 10 pM, at least 100 pM, at least 1 nanomolar (nM), at least 40 nM, at least 400 nM, at least 500 nM, at least, 600 nM, at least 700 nM, at least 800 nM, at least 900 nM, or 1 pM, and/or no greater than 1 pM, no greater than 900 nM, no greater than 800 nM, no greater than 700 nM, no greater than 600 nM, no greater than 500 nM, no greater than 400 nM, no greater than 40 nM, no greater than 1 nM, no greater than 100 nM, no greater than 6 nM, or no greater than 3 nM.

[00173]

[00174] In one embodiment, a deamination reaction can include an RNAse. RNase A has been implicated in increasing activity of cytidine deaminases (Bransteitter et al., Proceedings of the National Academy of Sciences of the United States of America 100, no. 7 (2003): 4102-7. doi.org/10.1073/pnas.0730835100). When activity of an altered cytidine deaminase of the present disclosure was determined in the presence of RNAse A the opposite was observed. When RNAse A was included in the reaction, an altered cytidine deaminase having cytosine-defective deaminase activity (i.e., converts 5mC to T at a greater rate than converting C to U) had reduced activity, and the reduced activity was more pronounced for off-target cytidine deamination. Thus, RNAse A resulted in greater selectivity for deamination of 5mC compared to C. An RNAse A can be included in a deamination reaction at a concentration from at least 1 microgram/milliliter (ug/ml) to no greater than 20 pM. For instance, the concentration of RNAse A can be at least 1 ug/ml , at least 2 ug/ml, at least 3 ug/ml, at least 4 ug/ml, 5 ug/ml, 6 ug/ml, 7 ug/ml, 8 ug/ml, or 9 ug/ml, and/or no greater than 50 ug/ml, no greater than 40 ug/ml, no greater than 30 ug/ml, no greater than 20 ug/ml, no greater than 19 ug/ml, no greater than 18 ug/ml, no greater than 17 ug/ml, no greater than 16 ug/ml, no greater than 15 ug/ml, no greater than 14 ug/ml, no greater than 13 ug/ml, no greater than 12 ug/ml, or no greater than 11 ug/ml. In one embodiment, the concentration of RNAse A is from 2 ug/ml to 13 ug/ml, or from 5 ug/ml to 10 ug/ml. [00175] C. Improvement of TET-based methods with CpG binding protein

[00176] One aspect of the present disclosure is methods, compositions, and kits related to use of a CpG binding protein and a ten-eleven translocase (TET) protein for detecting methylated nucleotides. The methods of using a CpG binding protein and a TET protein can be easily integrated into essentially any application that includes sequencing library preparation, such as whole genome, accessible (e.g., ATAC), and conformational state (e.g., HiC). It can be particularly useful in essentially any application using low input DNA or RNA such as, but not limited to, single cell combinatorial indexing (sci) methods like sci-WGS-seq, sci-MET-seq, and sci-ATAC-seq, sci-RNA-seq, and cell free DNA-based methods. Instead of using destructive methods for identifying methylated cytosines, integrating a TET protein into sequencing library production results in the more efficient enzyme-catalyzed conversion of modified cytosines during generation of target nucleic acids, thereby permitting more efficient detection of modified cytosines. A multitude of sequencing library methods are known to a skilled person that can be used in the construction of whole-genome or targeted libraries (see, for instance, Sequencing Methods Review, available on the world wide web at illumina.com/content/dam/illumina- marketing/documents/products/research_reviews/sequencing-met hods-review.pdf).

[00177] In one embodiment, the TET protein is a naturally occurring wild-type TET protein including, but not limited to, those described by Tahiliani (Science. 2009 May 15; 324(5929): 930-935. doi: 10.1126/science.H70116) and U.S. Published Application 2011/0236894. Wild-type TET protein has the activity of catalyzing the transfer of an oxygen molecule to the N5 methyl group on 5mC resulting in the formation of 5hmC. Wild-type TET further catalyzes the oxidation of 5hmC to 5-formylcytosine (5fC), and the oxidation of 5fC to form 5-carboxylcytosine (5caC) to result in 5caC or a mixture of 5fC and 5caC (FIG. 40). The inventor has found that unmethylated sequences, such as unmethylated cytosine, inhibit the iterative conversion of 5mC, 5hnC, and 5fC to 5caC (FIG. 2), which leads to decreased beta values after sequencing. The use of a CpG binding protein with a wild-type TET protein is expected to reduce the inhibitory effect of unmethylated nucleotides and lead to improved beta values, which will result in improved quality metrics in sequencing reactions. [00178] In general, methods for using a CpG binding protein and TET protein include contacting target nucleic acids, e.g., DNA or RNA, with the CpG binding protein and the enzyme, under conditions suitable for conversion of modified cytidines, such as 5mC, to 5fC and 5caC. The target nucleic acids can be contacted with a CpG binding protein before adding a TET protein, or both CpG binding protein and TET protein can be added to the target nucleic acids at the same time. Because amplification of DNA does not preserve the modification status of cytidine (e.g., the methylation status of 5mC is not retained), use of a CpG binding protein and a TET protein typically occurs before amplification of target nucleic acid. Target nucleic acids can be contacted with a CpG binding protein and a TET protein at essentially any time in a method before an amplification. For instance, target nucleic acids can be contacted with a CpG binding protein and a TET protein while the nucleic acids are inside a cell or a nucleus, after isolation of genomic or cell free DNA or mRNA, before or after fragmentation, or before or after tagmentation. The skilled person will recognize that target nucleic acids can be contacted with a CpG binding protein and a TET protein after addition of a universal sequence and/or an adapter, provided the universal sequence and/or an adapter is not added by amplification.

[00179] In one embodiment, the method further includes reducing the CpG binding protein/TET- treated DNA or RNA with a borane derivative. Examples of suitable borane derivatives include, but are not limited to, pyridine borane or 2-pi coline borane (pic- BEE), borane, sodium borohydride, sodium cyanoborohydride, and sodium triacetoxyborohydride. Borane reduction results in converting 5fC and 5caC to dihydrouracil (DHU), a uracil derivative that is recognized by both DNA and RNA polymerases as thymine (T). Thus, when CpG binding protein/TET and borane treated nucleic acids are sequenced, methylated cytosine nucleotides are identified as thymine and unmethylated nucleotides are identified as cytosine. This type of method of using CpG binding protein/TET followed by borane reduction is based on a method referred to as TET Assisted Pyridine borane Sequencing (TAPS, U.S. Published Patent Application No. 2022/0213543). Sequencing of the sample nucleic acid after treatment with a CpG binding protein/TET and borane reduction and optional comparison to a reference sequence permits easy identification of C to T point mutations, and these point mutations are inferred as sites having a methylated cytosine 5mC, 5hmC, 5fC, or 5caC. [00180] Optionally, the method can further include treatment of target nucleic acids with glucosyltransferase prior to the CpG binding protein/TET treatment or CpG binding protein/cytidine deaminase treatment. A glucosyltransferase enzyme can be used to glucosylate 5hmC but not 5mC. Glucosyltransferase enzymes are known to those of skill in the art, and include, for example, p-glucosyltransferase (PGT or BGT). A non-limiting example is T4 p-glucosyltransferase, which is commercially available (NEB). Enzymes, buffers, and conditions for performing glucosylation of 5hmC are known in the art, as exemplified by the methods disclosed in Schutsky et al., Nature biotechnology, 10.1038/nbt.4204. 8 Oct. 2018, doi: 10.1038/nbt.4204. A specific example of modification of 5hmC to protect it from cytidine deaminase activity is provided below in Example 9.

[00181] In one embodiment, the method includes treating the target nucleic acid with a glucosyltransferase prior to CpG binding protein/TET treatment. Glucosyltransferase treatment (e.g., using BGT) adds a glucose to 5hmC, converting it to 5gmC and thereby protecting 5hmC from TET-mediated oxidation and borane reduction. As a result, only 5mC is oxidized by TET and reduced by borane. This type of method of using BGT treatment, followed by CpG binding protein/TET treatment and borane reduction is based on a method referred to as TAPS-beta or TAPSB (U.S. Published Patent Application No. 2022/0213543). Sequencing of the sample nucleic acid after treatment with a BGT followed by CpG binding protein/TET and borane reduction and optional comparison to a reference sequence permits easy identification of C to T point mutations, and these point mutations are inferred as sites having a 5mC. 5hmC sites can be identified by comparison with the sequence data from the method using CpG binding protein/TET followed by borane reduction.

[00182] In one embodiment, the method includes treating the target nucleic acid with a glucosyltransferase prior to CpG binding protein/cytidine deaminase treatment. By treating the target nucleic acid with a glucosyltransferase (e.g., 0GT), 5hmC is protected from the deaminase activity of the cytidine deaminase enzyme. Thus, 5mC will be detected in downstream readout, such as sequencing, PCR, array, and the like, as a thymidine. In contrast, any protected 5hmC sites will be detected as cytosine in the same readout.TET enzymes useful in the methods described herein include human TET1, TET2, and TET3; murine Tetl, Tet2, and Tet3; Naegleria TET (NgTET); Coprinopsis cinerea TET (CcTET) and derivatives or analogues thereof. An example of a Tet protein is Tetl . A specific example of a Tetl protein is depicted at SEQ ID NO:60 (FIG. 37H, Genbank accession number NP 085128. Another example of a Tet protein is Tet2. A specific example of a Tet2 protein is depicted at SEQ ID NO:61 (FIG. 371, Genbank accession number NP 001120680). Another example of a Tet protein is Tet3. A specific example of a Tet3 protein is depicted at SEQ ID NO:62 (FIG. 37J, Genbank accession number NP_001274420).

[00183] In one embodiment, a TET protein has an amino acid sequence that is structurally similar to a reference TET protein disclosed herein. In one embodiment, a reference TET protein is one that includes the amino acid sequence of SEQ ID NO:60-62.

[00184] As used herein, a TET protein may be "structurally similar" to a reference TET protein if the amino acid sequence of the TET protein possesses a specified amount of sequence similarity and/or sequence identity compared to the reference TET protein. Structural similarity of two amino acid sequences can be determined as described herein. A candidate TET protein that has structural similarity with a reference TET protein and TET oxidation activity is included in the scope of the present disclosure. In one embodiment the amino acid sequence of a TET protein having sequence similarity to a reference sequence may include conservative substitutions of amino acids present in that reference sequence. Conservative substitutions are described herein.

[00185] Thus, as used herein, reference to a TET protein as described herein, such as reference to the amino acid sequence of one or more SEQ ID NOs of a TET protein described herein, can include a protein with at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 86%, at least 87%, at least 88%, at least 89%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99% amino acid sequence similarity to the reference TET protein.

[00186] Alternatively, as used herein, reference to a TET protein as described herein, such as reference to the amino acid sequence of one or more SEQ ID NOs described herein can include a protein with at least 60%, at least 65%, at least 70%, at least 75%, at least 80%, at least 85%, at least 86%, at least 87%, at least 88%, at least 89%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99% amino acid sequence identity to the reference TET protein.

[00187] D. Polynucleotides encoding CpG binding proteins, cytidine deaminases, and TET proteins

[00188] CpG binding proteins, cytidine deaminases, and TET proteins described herein also may be identified in terms of the polynucleotide that encodes the protein. Thus, this disclosure provides polynucleotides that encode a CpG binding protein, cytidine deaminase, or TET protein described herein or hybridize, under standard hybridization conditions, to a polynucleotide that encodes a CpG binding protein, cytidine deaminase, or TET protein described herein, and the complements of such polynucleotide sequences.

[00189] A polynucleotide as described herein can include any polynucleotide that encodes a CpG binding protein, cytidine deaminase, or TET protein described herein. Thus, the nucleotide sequence of the polynucleotide may be deduced from the amino acid sequence that is to be encoded by the polynucleotide. A CpG binding protein, cytidine deaminase, or TET protein can be encoded by multiple codons, and certain translation systems (e.g., prokaryotic or eukaryotic cells) often exhibit codon bias, e.g., different organisms often prefer one of the several synonymous codons that encode the same amino acid. As such, polynucleotides presented herein are optionally "codon optimized," meaning that the polynucleotides are synthesized to include codons that are preferred by the particular translation system being employed to express the protein. For example, when it is desirable to express the protein in a bacterial cell (or even a particular strain of bacteria), the polynucleotide can be synthesized to include codons most frequently found in the genome of that bacterial cell, for efficient expression of CpG binding protein, cytidine deaminase, or TET protein. A similar strategy can be employed when it is desirable to express the CpG binding protein, cytidine deaminase, or TET protein in a eukaryotic cell, e.g., the nucleic acid can include codons preferred by that eukaryotic cell.

[00190] A polynucleotide described herein may also, advantageously, be included in a suitable expression vector to express the CpG binding protein, cytidine deaminase, or TET protein encoded therefrom in a suitable host. Incorporation of cloned DNA into a suitable expression vector for subsequent transformation of a host cell and subsequent selection of the transformed cells is well known to those skilled in the art as provided in Sambrook et al. (1989), Molecular cloning: A Laboratory Manual, Cold Spring Harbor Laboratory. Suitable host cells include, but are not limited to, E. coli and S. cerevisiae.

[00191] Such an expression vector includes a vector having a polynucleotide described herein operably linked to heterologous regulatory sequences, such as promoter regions, that are capable of effecting expression of said DNA fragments. The term "operably linked" refers to a juxtaposition wherein the components described are in a relationship permitting them to function in their intended manner. Such vectors may be transformed into a suitable host cell to provide for the expression of CpG binding protein, cytidine deaminase, or TET protein.

[00192] The nucleic acid molecule may encode a mature protein or a protein having a pro-sequence, including that encoding a leader sequence on the preprotein which is then cleaved by the host cell to form a mature protein. The vectors may be, for example, plasmid, virus or phage vectors provided with an origin of replication, and optionally a promoter for the expression of said nucleotide and optionally a regulator of the promoter. The vectors may contain one or more selectable markers, such as, for example, an antibiotic resistance gene.

[00193] Regulatory elements required for expression include promoter sequences to bind RNA polymerase and to direct an appropriate level of transcription initiation and also translation initiation sequences for ribosome binding. For example, a bacterial expression vector may include a promoter such as the lac promoter and for translation initiation the Shine- Dalgarno sequence and the start codon AUG. Similarly, a eukaryotic expression vector may include a heterologous or homologous promoter for RNA polymerase II, a downstream polyadenylation signal, the start codon AUG, and a termination codon for detachment of the ribosome. Such vectors may be obtained commercially or be assembled from the sequences described by methods well known in the art.

[00194] Transcription of DNA encoding an altered cytidine deaminase may be optimized by including an enhancer sequence in the vector. Enhancers are cis-acting elements of DNA that act on a promoter to increase the level of transcription. Vectors will also generally include origins of replication in addition to the selectable markersE, Making and isolating CpG binding proteins, cytidine deaminases, and TET proteins

[00195] Generally, polynucleotides encoding a CpG binding protein, cytidine deaminase, or a TET protein described herein can be made by cloning, recombination, in vitro synthesis, in vitro amplification and/or other available methods. A variety of recombinant methods can be used for expressing an expression vector that encodes a CpG binding protein, cytidine deaminase, or a TET protein described herein. Methods for making recombinant polynucleotides, expression, and isolation of expressed products are well known and described in the art.

[00196] Polynucleotides encoding wild type cytidine deaminases can be obtained from a source and subjected to mutagenesis to introduce one or more substitution mutations described herein. In general, any available mutagenesis procedure can be used for making a cytidine deaminase, CpG binding protein, or TET protein described herein. Procedures that can be used include, but are not limited to site-directed point mutagenesis, in vitro or in vivo homologous recombination, oligonucleotide-directed mutagenesis, mutagenesis by total gene synthesis, and many others known to a person skilled in the art.

[00197] Additional useful references for mutation, recombinant, and in vitro nucleic acid manipulation methods (including cloning, expression, PCR, and the like) include Berger and Kimmel, Guide to Molecular Cloning Techniques, Methods in Enzymology volume 152 Academic Press, Inc., San Diego, Calif. (Berger); Kaufman et al. (2003) Handbook of Molecular and Cellular Methods in Biology and Medicine Second Edition Ceske (ed) CRC Press (Kaufman); The Nucleic Acid Protocols Handbook Ralph Rapley (ed) (2000) Cold Spring Harbor, Humana Press Inc (Rapley); Chen et al. (ed) PCR Cloning Protocols, Second Edition (Methods in Molecular Biology, volume 192) Humana Press; and in Viljoen et al. (2005) Molecular Diagnostic PCR Handbook Springer, ISBN 1402034032.

[00198] In addition, many kits are commercially available for the purification of plasmids or other relevant nucleic acids from cells. An isolated polynucleotide can be further manipulated to produce other polynucleotides, used to transfect or transform cells, incorporated into related vectors and introduced into cells for expression, and/or the like. Typical cloning vectors contain transcription and translation terminators, transcription and translation initiation sequences, and promoters useful for regulation of the expression of the particular target nucleic acid. The vectors optionally comprise generic expression cassettes containing at least one independent terminator sequence, sequences permitting replication of the cassette in eukaryotes, or prokaryotes, or both, (e.g., shuttle vectors) and selection markers for both prokaryotic and eukaryotic systems. Vectors are suitable for replication and integration in prokaryotes, eukaryotes, or both.

[00199] Other useful references, e.g. for cell isolation and culture (e.g., for subsequent nucleic acid isolation) include Freshney (1994) Culture of Animal Cells, a Manual of Basic Technique, third edition, Wiley-Liss, New York and the references cited therein; Payne et al. (1992) Plant Cell and Tissue Culture in Liquid Systems John Wiley & Sons, Inc. New York, N.Y.; Gamborg and Phillips (eds) (1995) Plant Cell, Tissue and Organ Culture; Fundamental Methods Springer Lab Manual, Springer-Verlag (Berlin Heidelberg New York); and Atlas and Parks (eds) The Handbook of Microbiological Media (1993) CRC Press, Boca Raton, Fla. Construction of vectors containing a nucleic acid encoding a CpG binding protein, cytidine deaminase, or a TET protein described herein employs standard ligation techniques known in the art. See, e.g., Sambrook et al, Molecular Cloning: A Laboratory Manual., Cold Spring Harbor Laboratory Press (1989) or Ausubel, R.M., ed. Current Protocols in Molecular Biology (1994).

[00200] A variety of protein isolation and detection methods are known and can be used to isolate a CpG binding protein, cytidine deaminase, or a TET protein, e.g., from recombinant cultures of cells expressing a CpG binding protein, cytidine deaminase, or a TET protein described herein. A variety of protein isolation and detection methods are well known in the art, including, e.g., those set forth in R. Scopes, Protein Purification, Springer-Verlag, N.Y. (1982); Deutscher, Methods in Enzymology Vol. 182: Guide to Protein Purification, Academic Press, Inc. N.Y. (1990); Sandana (1997) Bioseparation of Proteins, Academic Press, Inc.; Bollag et al. (1996) Protein Methods, 2nd Edition Wiley-Liss, NY; Walker (1996) The Protein Protocols Handbook Humana Press, NJ, Harris and Angal (1990) Protein Purification Applications: A Practical Approach IRL Press at Oxford, Oxford, England; Harris and Angal Protein Purification Methods: A Practical Approach IRL Press at Oxford, Oxford, England; Scopes (1993) Protein Purification: Principles and Practice 3rd Edition Springer Verlag, NY; Janson and Ry den (1998) Protein Purification: Principles, High Resolution Methods and Applications, Second Edition Wiley-VCH, NY; and Walker (1998) Protein Protocols on CD-ROM Humana Press, NJ; and the references cited therein. Additional details regarding protein purification and detection methods can be found in Satinder Ahuja ed., Handbook of Bioseparations, Academic Press (2000).

[00201] F. Target nucleic acids

[00202] The target nucleic acids contacted with a CpG binding protein, a cytidine deaminase, a TET protein, or any combination thereof and used in the methods, compositions, and kits provided herein may be essentially any nucleic acid of known or unknown sequence. Sequencing may result in determination of the sequence of the whole or a part of the target molecule. In one embodiment, target nucleic acids can be processed into templates suitable for amplification by the placement of universal amplification sequences, e.g., sequences present in a universal adaptor, at the ends of each target fragment.

[00203] Target nucleic acids useful in the methods of the present disclosure are described herein. The target nucleic acid may have or be suspected of having one or more modified cytidines. A modified cytidine present on a substrate single-stranded (ss) DNA or RNA includes, but is not limited to, 5-methyl cytosine (5mC), 5 -hydroxymethyl cytosine (5hmC), 5-formyl cytosine (5fC), and 5-carboxy cytosine (5CaC) (FIG. 2). In one embodiment, the modified cytidine is 5-methyl cytosine. In one embodiment, the modified cytidine is 5 -hydroxymethyl cytosine. Methods that use double stranded target DNA for generating a sequencing library can be modified to include denaturation to convert the double stranded target DNA to ssDNA. In some embodiments, dsDNA that is used in a tagmentation reaction or for adapter attachment can be denatured and then treated with an cytidine deaminase. Conditions for denaturation are known and routine. In those embodiments where ssDNA is contacted with an cytidine deaminase and subsequently used in a process that requires dsDNA, e.g., addition of universal adapters by tagmentation or ligation, the ssDNA can be converted to dsDNA using routine methods. [00204] In some embodiments, an altered cytidine deaminase as presented herein can be used to differentiate between 5-methyl cytosine (5mC) and 5 -hydroxymethyl cytosine (5hmC). In such an embodiment, a sample ofDNA suspected of including single-stranded DNA comprising at least one 5-methyl cytosine (5mC) or 5 -hydroxymethyl cytosine (5hmC) is modified to prevent an altered cytidine deaminase from converting 5hmC to thymidine. Methods for blocking deaminase activity are known in the art, and any one of a number of methods can be used to protect 5hmC from deaminase activity. As one example, target DNA can be treated to modify 5hmC but not 5mC such that 5hmC is an unsuitable substrate for cytidine deaminase activity. In a specific example, a glucosyltransferase enzyme can be used to glucosylate 5hmC but not 5mC.

[00205] In some embodiments, a method that includes using a CpG binding protein and a cytidine deaminase can be used in in vitro diagnostic (IVD) approaches for profiling methylation in a locus-specific manner. Current methods for methylation biomarker detection typically include digestion of genomic DNA with methylation-sensitive enzymes and then quantitative PCR (qPCR) at a locus of interest to quantify the extent of restriction enzyme digestion, and therefore the percent methylation at that site. This is followed by mismatchsensitive qPCR of bisulfite-treated DNA, where 5mC is read out as a lack of 5mC>T conversion. These methods, however, have drawbacks. The recognition site of the methylsensitive restriction enzyme must be present in the methylated region of the target locus. Bisulfite treatment requires a large quantity of starting DNA and results in conversion to a low complexity genome (unmethylated cytosines - which represent the majority of cytosines in the genome - are converted to U and read as T). This reduced complexity of the genomic template constrains the design of qPCR primers that hybridize specifically to the locus of interest. During bisulfite conversion, DNA is intrinsically damaged or lost, which can hinder downstream analysis. DNA damage decreases coverage uniformity of the genome, which can lead to bias coverage. Furthermore, incomplete bisulfite conversion has the potential to adversely affect results, since it can exaggerate DNA methylation levels (Sam et al., PLoS One. 2018; 13(6); Ehrich et al., Nucleic Acids Res. Oxford University Press; 2007;35: e29). [00206] 5mC to T conversion by a cytidine deaminase described herein obviates the need for restriction enzymes or bisulfite treatment, and preserves DNA complexity. The resulting modifications of one or more cytosines can be detected using established in vitro diagnostic (IVD) approaches for profiling methylation in a locus-specific manner. Examples of approaches include detection of 5mC loci via amplification, e.g., quantitative PCT (qPCR) (Example 21), detection of 5mC loci using a CRISPR-based system, e.g., CRISPR-Casl2 (Example 22), spatial detection of 5mC using molecular cytogenic methods, e.g., fluorescence in situ hybridization (FISH) (Example 23), and array-based detection of 5mC (Example 24). In one embodiment, in vitro diagnostic (IVD) approaches for profiling methylation in a locus-specific manner use one or more primers to anneal to a predetermined sequence that may include one or more modified cytosines. After treatment of target nucleic acids with an cytidine deaminase, the modified cytosines present in the target nucleic acids are converted as described herein (e.g., 5mC is converted to T), and primers can be easily designed to anneal with higher affinity to a predetermined sequence when it includes nucleotides resulting from the deaminase treatment (e.g., a T nucleotide where a 5mC was present prior to treatment). For example, as described in Example 21, primers used for an amplification bind with greater affinity to a nucleic acid that includes T nucleotides where 5mC nucleotides were present prior to treatment. The annealing of a primer to a predetermined sequence that includes the expected 5mC to T conversion(s) allows one to infer the location of a modified cytosine in the untreated target nucleic acid. A primer that binds with greater affinity to a nucleic acid that includes T nucleotides where 5mC nucleotides were present prior to treatment can include at least 1, at least 2, at least 3, at least 4 or at least 5 nucleotides that will base-pair with a nucleotide that results from conversion of 5mC to T, i.e., an adenine (A), and when amplification is used, then a second primer for the reverse strand that has a T instead of guanine (G).

[00207] In some embodiments, target nucleic acids obtained from a subject can be treated with a CpG binding protein and a cytidine deaminase and/or a TEET to result in converted nucleic acids, and a pattern of cytosine modification can be identified in the converted nucleic acids. The pattern of cytosine modification can optionally be compared with the pattern of cytosine modification in a reference nucleic acid. In embodiments where a pattern of cytosine modification correlates with a disease or condition, the method can be used in diagnostic or prognostic applications. For instance, the subject can have or be at risk of having a disease or condition, and the reference nucleic acid can be from a normal subject, e.g., a subject that does not have and is not at risk for the disease or condition. The pattern of cytosine modification can be associated with a disease or condition (e.g., the target nucleic acid can be a predetermined sequence), and identification in the subject of a pattern of cytosine modification associated with a disease or condition can indicate the subject has or is at risk of having the disease or condition. For instance, a pattern of cytosine modification can be linked in-cis to a coding region that is correlated with a disease or condition and identification of that pattern, or absence of that pattern, in the subject can be used for diagnosis or prognosis. In one embodiment, the coding region can be one that is transcriptionally active or transcriptionally inactive in a reference nucleic acid. The comparison of the converted nucleic acid to the reference nucleic acid can include determining if the pattern of cytosine modification of the converted nucleic acid indicates the coding region is transcriptionally active or transcriptionally inactive in the subject. When that coding region is associated with a disease or condition, the status of transcriptional activity can be used for diagnosis or prognosis.

[00208] Comparison of a pattern of cytosine modification in a subject can also be used in identifying changes in a pattern of cytosine modification in a subject over time. For instance, a subject can have a disease or conditions and is undergoing treatment, or a subject had a disease or condition and is cured (e.g., the subject was treated and no signs of the disease or condition are present) or in remission (e.g., the subject was treated and signs of the disease or condition are reduced). Target nucleic acids from the subject at different times, e.g., before treatment started, during treatment, after treatment is stopped, can be compared and a pattern of cytosine modification of a sequence, e.g., a predetermined sequence compared and used to determine the progress of a treatment or the status of the disease or condition in the subject.

[00209] In some embodiments where detection of 5mC nucleotides uses amplification, the use of a polymerase that disfavors uracil can aid in reducing the amplification of treated target nucleic acids that include spurious C to U conversion that may result from use of an cytidine deaminase. B-family polymerases are known to exhibit “uracil read-ahead” function which causes stalling of the polymerase at uracil residues (Greagg et al., 1999, PNAS USA; 96(16):9045-50). Examples of B-family polymerases that disfavor uracil include archaeal B-family polymerases from Pyrococcus furiosus (Pfu), Thermococcus kodakarensis (KOD), Thermococcus litoralis (Tli/Vent), Pyrococcus woesei (Pwo), and Thermococcus fumicolans (Tfu). Other examples of uracil-disfavoring polymerases include PhusionTM, Q5®, and Kapa HiFiTM. In other embodiments where amplification of nucleic acids containing uracil nucleotides is desired, the use of a uracil tolerant polymerase can be used. Examples of uracil-tolerant polymerases include PhusionUTM, Q5U®, KapaUTM, Taq, and Dpo4.

[00210] Target nucleic acids are typically derived from primary nucleic acids present in a sample, such as a biological sample. The primary nucleic acids may originate as DNA or RNA. DNA primary nucleic acids may originate in double-stranded DNA (dsDNA) form (e.g., genomic DNA, genomic DNA fragments, cell-free DNA, and the like) from a sample or may originate in single-stranded form from a sample. RNA primary nucleic acids may be mRNA or non-coding RNA, e.g., microRNA or small interfering RNA. The precise sequence of the polynucleotide molecules from a primary nucleic acid sample is generally not material to the disclosure and may be known or unknown.

[00211] The primary nucleic acid molecules may represent the entire genetic complement of an organism, e.g., genomic DNA molecules which include both intron and exon sequences, as well as non-coding regulatory sequences such as promoter and enhancer sequences. The primary nucleic acid molecules may represent the entire genetic complement of specific cells of an organism, e.g., from tumor cells, where the genomic DNA molecules which include both intron and exon sequences, as well as non-coding regulatory sequences such as promoter and enhancer sequences. In one embodiment, particular subsets of genomic DNA can be used, such as, for example, particular chromosomes, DNA associated with open chromatin, DNA associated with closed chromatin, or one or more specific sequences such as a region of a specific gene (e.g., targeted sequencing). In one embodiment, the primary nucleic acid molecules may represent a particular subset of DNA, e.g., DNA having a specific sequence that anneals with a primer such as one used for targeted sequencing or target enrichment. In one embodiment, a particular subset of DNA can be used, such as cell-free DNA, which can include DNA of the subject including DNA from normal cells, DNA from diseased cells such as tumor cells, and/or DNA from fetal cells.

[00212] The primary nucleic acid molecules may represent the entire transcriptome of cells of an organism, e.g., mRNA molecules. The primary nucleic acid molecules may represent the entire transcriptome of specific cells of an organism, e.g., from tumor cells or for instance the cells of a tissue. In one embodiment, the primary nucleic acid molecules may represent a particular subset of mRNA, e g., mRNA having a specific sequence that anneals with a primer such as one used for targeted sequencing or target enrichment.

[00213] A sample, such as a biological sample, can include nucleic acid molecules obtained from biopsies, tumors, scrapings, swabs, blood, mucus, urine, plasma, semen, hair, laser capture micro-dissections, surgical resections, and other clinical or laboratory obtained samples. In some embodiments, the sample can be an epidemiological, agricultural, forensic or pathogenic sample. In some embodiments, the sample can include cultured cells. In some embodiments, the sample can include nucleic acid molecules obtained from an animal such as a human or mammalian source. In another embodiment, the sample can include nucleic acid molecules obtained from a non-mammalian source such as a plant, bacteria, virus, or fungus. In some embodiments, the source of the nucleic acid molecules may be an archived or extinct sample or species.

[00214] Additional non-limiting examples of sources of biological samples can include whole organisms as well as a sample obtained from a patient. The biological sample can be obtained from any biological fluid or tissue and can be in a variety of forms, including liquid fluid and tissue, solid tissue, and preserved forms such as dried, frozen, and fixed forms. The sample may be of any biological tissue, cells or fluid. Such samples include, but are not limited to, sputum, blood, serum, plasma, fecal, mucosal, intestinal fluid, skin (e.g. punch biopsy), blood cells (e.g., white cells), ascitic fluid, urine, saliva, tears, sputum, vaginal fluid (discharge), washings obtained during a medical procedure (e.g., pelvic or other washings obtained during biopsy, endoscopy or surgery), tissue, nipple aspirate, core or fine needle biopsy samples, cell-containing body fluids, peritoneal fluid, and pleural fluid, or cells therefrom, and free floating nucleic acids such as cell-free circulating DNA. Biological samples may also include sections of tissues such as frozen or fixed sections taken for histological purposes or micro-dissected cells or extracellular parts thereof. In some embodiments, the sample can be a blood sample, such as, for example, a whole blood sample. In another example, the sample is an unprocessed dried blood spot (DBS) sample. In yet another example, the sample is a formalin-fixed paraffin-embedded (FFPE) sample. In yet another example, the sample is a saliva sample. In yet another example, the sample is a dried saliva spot (DSS) sample.

[00215] Exemplary biological samples from which target nucleic acids can be derived include, for example, those from a eukaryote, for instance a mammal, such as a rodent, mouse, rat, rabbit, guinea pig, ungulate, horse, sheep, pig, goat, cow, cat, dog, primate, human or nonhuman primate; a plant, such as Arabidopsis thaliana, corn, sorghum, oat, wheat, rice, canola, or soybean; an algae, such as Chlamydomonas reinhardtii, a nematode such as Caenorhabditis elegans', an insect, such as Drosophda melanogaster , mosquito, fruit fly, honey bee or spider; a fish, such as zebrafish; a reptile; an amphibian, such as a frog or Xenopus laevis,' a Dictyostelium discoideum, a fungi, such as Pneumocystis carinii, Takifiigu rubripes, yeast, Saccharamoyces cerevisiae, or Schizosaccharomyces pombc, or a protozoan such as Plasmodium falciparum. Target nucleic acids can also be derived from a prokaryote such as a bacterium, Escherichia coli, Staphylococcus ox Mycoplasma pneumoniae,' an archaeon; a virus such as Hepatitis C virus or human immunodeficiency virus; or a viroid. Target nucleic acids can be derived from a homogeneous culture or population of organisms described herein or alternatively from a collection of several different organisms, for example, in a community or ecosystem.

[00216] In some embodiments, a biological sample includes tissue that is processed to obtain the desired primary nucleic acids. In some embodiments, cells are used to obtain the desired primary nucleic acids. In some embodiments, nuclei are used to obtain the desired primary nucleic acids. The method can further include dissociating cells, and/or isolating nuclei from cells. Methods for isolating cells and nuclei from tissue are available (WO 2019/236599). [00217] In some embodiments, nucleic acids present in tissue, in cells, or in isolated nuclei can be processed depending on the desired read-out. For instance, nucleic acids can be fixed during processing, and useful fixation methods are available (WO 2019/236599). Fixation can be useful to preserve a sample or maintain contiguity of analytes from a sample, a cell, or a nucleus. Fixation methods preserve and stabilize tissue, cell, and nucleus morphology and architecture, inactivate proteolytic enzymes, strengthen samples, cells, and nuclei so they can withstand further processing and staining, and protect against contamination. Examples of methods where fixation can be useful include, but are not limited to, whole genome sequencing of isolated nuclei and chromosome conformation capture methods such as Hi-C. Common methods of fixation include perfusion, immersion, freezing, and drying (Srinivasan et al., Am J Pathol. 2002 Dec; 161(6): 1961— 1971. doi: 10.1016/S0002-9440(10)64472-0). In some embodiments such as whole genome sequencing, isolated nuclei can be processed to dissociate nucleosomes from DNA while leaving the nuclei intact, and methods for generating nucleosome-free nuclei are available (WO 2018/018008).

[00218] In some embodiments, primary nucleic acids in bulk, e.g., from a plurality of cells, can be used to produce a sequencing library as described herein. In other embodiments, individual cells or nuclei can be used as sources of primary nucleic acids to obtain sequence information from single cells and nuclei. Many different single cell library preparation methods are known in the art, including, but not limited to, Drop-seq, Seq-well, and single cell combinatorial indexing ("sci-") methods. . Companies providing single cell products and related technologies include, but are not limited to, Illumina, 10X Genomics, Takara Biosciences, BD Biosciences, Bio-Rad Laboratories, Icellbio, Isoplexis, CellSee, NanoSelect, and Dolomite Bio. Sci-seq is a methodological framework that employs splitpool barcoding to uniquely label the nucleic acid contents of large numbers of single cells or nuclei. Typically, the number of nuclei or cells can be at least two. The upper limit is dependent on the practical limitations of equipment (e.g., multi-well plates, number of indexes) used in other steps of the methods as described herein. The number of nuclei or cells that can be used is not intended to be limiting and can number in the billions. [00219] The target nucleic acids used in the methods, compositions, and kits of the present disclosure can be derived by fragmentation. Random fragmentation refers to the fragmentation of a polynucleotide molecule from a primary nucleic acid sample in a nonordered fashion by enzymatic, chemical, or mechanical methods. Such fragmentation methods are known in the art and use standard methods (Sambrook and Russell, Molecular Cloning, A Laboratory Manual, third edition). Moreover, random fragmentation is designed to produce fragments irrespective of the sequence identity or position of nucleotides comprising and/or surrounding the break. In one embodiment, the random fragmentation is by mechanical means such as nebulization or sonication to produce fragments of about 50 base pairs in length to about 1500 base pairs in length, still more particularly 50-700 base pairs in length, yet more particularly 50-400 base pairs in length. Most particularly, the method is used to generate smaller fragments of from 50-150 base pairs in length

[00220] Fragmentation of polynucleotide molecules by mechanical means (nebulization, sonication, and Hydroshear, for example) results in fragments with a heterogeneous mix of blunt and 3'- and 5'-overhanging ends. It is therefore desirable to repair the fragment ends using methods or kits (such as the Lucigen DNA terminator End Repair Kit) known in the art to generate ends that are optimal for insertion, for example, into blunt sites of cloning vectors. In a particular embodiment, the fragment ends of the population of nucleic acids are blunt ended. More particularly, the fragment ends are blunt ended and phosphorylated. The phosphate moiety can be introduced via enzymatic treatment, for example, using polynucleotide kinase.

[00221] In a particular embodiment, the target fragment sequences are prepared with single overhanging nucleotides by, for example, activity of certain types of DNA polymerase such as Taq polymerase or Klenow exo-minus polymerase which has a non-template-dependent terminal transferase activity that adds a single deoxynucleotide, for example, deoxyadenosine (A) to the 3' ends of a DNA molecule, for example, a PCR product. Such enzymes can be used to add a single nucleotide 'A' to the blunt ended 3' terminus of each strand of the double-stranded target fragments. Thus, an 'A' could be added to the 3' terminus of each end repaired strand of the double-stranded target fragments by reaction with Taq or Klenow exo-minus polymerase, while the universal adapter polynucleotide construct could be a T-construct with a compatible 'T' overhang present on the 3' terminus of each region of double-stranded nucleic acid of the universal adapter. This end modification also prevents self-ligation of both vector and target such that there is a bias towards formation of target nucleic acids having a universal adapter at each end.

[00222] In one embodiment, fragmentation can be accomplished using a process often referred to as tagmentation. Tagmentation uses a transposome complex and combines into a single step fragmentation and ligation to add universal adapters (WO 2016/130704). A transposome complex is a transposase bound to a transposase recognition site and can insert the transposase recognition site into a target nucleic acid in a process sometimes termed "tagmentation." In some such insertion events, one strand of the transposase recognition site may be transferred into the target nucleic acid. Such a strand is referred to as a "transferred strand." In one embodiment, a transposome complex includes a dimeric transposase having two subunits, and two non-contiguous transposon sequences. In another embodiment, a transposase includes a dimeric transposase having two subunits, and a contiguous transposon sequence.

[00223] Some embodiments can include the use of a hyperactive Tn5 transposase and a Tn5-type transposase recognition site (Goryshin and Reznikoff, J. Biol. Chem., 273:7367 (1998)), or MuA transposase and a Mu transposase recognition site comprising R1 and R2 end sequences (Mizuuchi, K., Cell, 35: 785, 1983; Savilahti, H, et al., EMBO J., 14: 4893, 1995). Tn5 Mosaic End (ME) sequences can also be used by a skilled artisan.

[00224] Examples of transposon sequences useful with the methods and compositions described herein are provided in U.S. Patent Application Pub. No. 2012/0208705, U.S. Patent Application Pub. No. 2012/0208724 and Int. Patent Application Pub. No. WO 2012/061832. In some embodiments, a transposon sequence includes a first transposase recognition site and a second transposase recognition site.

[00225] Some transposome complexes useful herein include a transposase having two transposon sequences. In some such embodiments, the two transposon sequences are not linked to one another, in other words, the transposon sequences are non-contiguous with one another. Examples of such transposomes are known in the art (see, for instance, U.S. Patent Application Pub. No. 2010/0120098).

[00226] In one embodiment, tagmentation is used to produce target nucleic acids that include different universal sequences at each end. This can be accomplished by using two types of transposome complexes, where each transposome complex includes a different nucleotide sequence that is part of the transferred strand.

[00227] A population of target nucleic acids can have an average strand length that is desired or appropriate for a particular application of the methods, compositions, or kits set forth herein. For example, the average strand length can be less than about 100,000 nucleotides, 50,000 nucleotides, 10,000 nucleotides, 5,000 nucleotides, 1,000 nucleotides, 500 nucleotides, 100 nucleotides, or 50 nucleotides. Alternatively or additionally, the average strand length can be greater than about 10 nucleotides, 50 nucleotides, 100 nucleotides, 500 nucleotides, 1,000 nucleotides, 5,000 nucleotides, 10,000 nucleotides, 50,000 nucleotides, or 100,000 nucleotides. The average strand length for a population of target nucleic acids can be in a range between a maximum and minimum value set forth herein. It will be understood that amplicons generated at an amplification site (or otherwise made or used herein) can have an average strand length that is in a range between an upper and lower limit selected from those exemplified above.

[00228] In some cases, a population of target nucleic acids can be produced under conditions or otherwise configured to have a maximum length for its members. For example, the maximum length for the members that are used in one or more steps of a method set forth herein or that are present in a particular composition can be less than 100,000 nucleotides, less than 50,000 nucleotides, less than 10,000 nucleotides, less than 5,000 nucleotides, less than 1,000 nucleotides, less than 500 nucleotides, less than 100 nucleotides, or less than 50 nucleotides. Alternatively or additionally, a population of target nucleic acids can be produced under conditions or otherwise configured to have a minimum length for its members. For example, the minimum length for the members that are used in one or more steps of a method set forth herein or that are present in a particular composition can be more than 10 nucleotides, more than 50 nucleotides, more than 100 nucleotides, more than 500 nucleotides, more than 1,000 nucleotides, more than 5,000 nucleotides, more than 10,000 nucleotides, more than 50,000 nucleotides, or more than 100,000 nucleotides. The maximum and minimum strand length for target nucleic acids in a population can be in a range between a maximum and minimum value set forth above. It will be understood that amplicons generated at an amplification site (or otherwise made or used herein) can have maximum and/or minimum strand lengths in a range between the upper and lower limits exemplified above.

[00229] In some embodiments, a sample can be enriched for sequences of interest, such as, for example, a predetermined sequence. For example, a subset of genes or regions of the genome are isolated and sequenced , or a subset of genes or regions of the genome are interrogated by other methods, such as a locus-specific in vitro diagnostic method. A predetermined sequence can be, for instance, one that can have a pattern of cytosine modification.

[00230] In some embodiments, target enrichment works by capturing genomic regions of interest by hybridization to target-specific probes that can be used to physically separate target DNA that has hybridized to bait probes from all other DNA in solution, which are then washed away. For example, some methods of enrichment use biotinylated probes, which are then isolated by magnetic pulldown with streptavidin-coated magnetic particles. In another example, some methods of enrichment use analyte arrays, also known as microarrays, that allow for the hybridization of predetermined sequences.

[00231] Enrichment can occur, for example, prior to treatment to detect a modified cytosine. In such embodiments, enriching a nucleic acid of interest, or a fragment thereof, such as enriching nucleic acids in a sample, may include any suitable enrichment techniques. In some embodiments, enrichment of nucleic acids may include enrichment through molecular inversion probes, in solution capture, pulldown probes, bait sets, hybrid capture, endonuclease digestion, DNase I hypersensitivity, and selective circularization. Preferably, enrichment prior to treatment to detect a modified cytosine does not include amplification. Enrichment can be achieved through negative selection of nucleic acids by eliminating undesired material. This sort of enrichment includes 'footprinting' techniques or 'subtractive' hybrid capture. During the former, the target sample is safe from nuclease activity through the protection of protein or by single and double stranded arrangements. During the latter, nucleic acids that bind 'bait' probes are eliminated.

[00232] In some embodiments, enriching after treatment to detect a modified cytosine can comprise amplification using target-specific primers in, for instance, standard PCR ir multiplex PCR. In some embodiments, amplification is performed subsequent to another form of enrichment. Typically, however, in embodiments where amplification is used for enrichment, the amplification step occurs after treatment with a CpG binding proteins and either a cytidine deaminase or a TET protein, to preserve methylation status of the target DNA. In some such embodiments, amplification can include PCR amplification or genome-wide amplification.

[00233] In some embodiments, enrichment can occur after treatment with an altered cytidine deaminase. Typically, methods used to identity methylated cytosines result in the loss of DNA complexity due to conversion of unmethylated DNA bases to uracil, resulting in lower complexity genome and limits the use of sequences that specifically hybridize to a predetermined sequence. Accordingly, typical methods for identifying methylated cytosines are more difficult to use in methods that include enrichment, such as hybridenrichment sequencing and amplicon-based targeted sequencing, after conversion of methylated cytosines. In contrast, because of (i) the 5mC to T conversion by altered cytidine deaminases and (ii) only a small percentage of cytosines are methylated and expected to be converted by an altered cytidine deaminase. Examples of enrichment-based methods that can be used after treatment of a target nucleic acid with an altered cytidine deaminase include but are not limited to analyte arrays, use of primers for selective amplification, CRISPR-Cas systems, and molecular cytogenic techniques such as FISH. Examples of arrays include, for instance, methylation arrays for interrogation of selected methylation sites across a genome (e.g., the Infinium Methyl ationEPIC BeadChip, Illumina).

[00234] The primary nucleic acid molecules include at least one modified cytosine. A modified cytidine present on a substrate single-stranded (ss) DNA or RNA includes, but is not limited to, 5-methyl cytosine (5mC), 5-hydroxymethyl cytosine (5hmC), 5-formyl cytosine (5fC), and 5-carboxy cytosine (5CaC) (FIG. 2). In one embodiment, the modified cytidine is 5-methyl cytosine. Methods that use double-stranded target DNA for generating a sequencing library can be modified to include denaturation to convert the ds target DNA to ss DNA. In some embodiments, dsDNA that is used in a tagmentation reaction or for adapter attachment can be denatured and then treated with a CpG binding protein and a cytidine deaminase or a TET protein. Conditions for denaturation are routine and include thermal or chemical denaturation. In those embodiments where ssDNA is contacted with a CpG binding protein and a cytidine deaminase or a TET protein and subsequently used in a process that requires dsDNA, e.g., addition of universal adapters by tagmentation or ligation, the ss DNA can be converted to ds DNA using routine methods.

[00235] G. Attachment of Universal Adapters

[00236] In some embodiments, a target nucleic acid used in a method, composition, or kit described herein can include a universal adapter attached to each end. A target nucleic acid having a universal adapter at each end can be referred to as a "modified target nucleic acid." Methods for attaching a universal adapter to each end of a target nucleic acid used in a method described herein are known to the person skilled in the art. The attachment can be through tagmentation using transposase complexes (WO 2016/130704), or through standard library preparation techniques using ligation (U.S. Pat. Pub. No. 2018/0305753). Attachment of a universal adapter to the ends of a target nucleic acid can occur before or after treatment of the target nucleic acid with a cytidine deaminase.

[00237] In some embodiments, adapters compatible with the methods and compositions described herein may include additional regions not generally included in sequencing adapters. For some helicases of interest, certain motifs may preferably improve recruitment of the helicase to a target nucleic acid. Examples of motifs that may be included in an adapter sequence include Y-shaped adapter structures, partial duplexes, ssDNA bubbles, forks, blunt duplexes, three-strand junctions, and others. In some embodiments, an adapter may include more than one motif to recruit a helicase. Either the 3' end or the 5' end may include a motif. In some embodiments, both the 3' and the 5' end include a motif. [00238] In one embodiment, double-stranded target nucleic acids from a sample, e.g., a fragmented sample that has been contacted with a cytidine deaminase and converted from singlestranded to double-stranded nucleic acids, are treated by first ligating identical universal adaptor molecules to the 5' and 3' ends of the double-stranded target nucleic acids. In one embodiment, the universal adapters are "matched" adapters of Y-adapters because the two strands of the adaptors are formed by annealing complementary polynucleotide strands. In one embodiment, the universal adapters used in the method of the disclosure are referred to as "mismatched" adaptors because the adaptors include a region of sequence mismatch, i.e., they are not formed by annealing fully complementary polynucleotide strands. The general features of mismatched adaptors are further described in Gormley et al., U.S. Pat. No. 7,741,463, and Bignell et al., U.S. Pat. No. 8,053,192,). The universal adaptor typically includes universal capture binding sequences that aid in immobilizing the target nucleic acids on an array for subsequent sequencing, and universal primer binding sites useful for the sequencing. In another embodiment, a sample is contacted with a cytidine deaminase and helicase or fusion protein as described herein. Double-stranded target nucleic acids from the sample are then subjected to tagmentation with a transposome complex that inserts a universal adapter, or sequences that can be used to add a universal adapter, into the target double-stranded nucleic acids.

[00239] A universal adapter can optionally include at least one index. An index can be used as a marker characteristic of the source of particular target nucleic acids on an array (U.S. Pat. No. 8,053,192). Generally, the index is a synthetic sequence of nucleotides that is part of the universal adapter which is added to the target nucleic acids as part of the library preparation step. Accordingly, an index is a nucleic acid sequence which is attached to each of the target molecules of a particular sample, the presence of which is indicative of, or is used to identify, the sample or source from which the target molecules were isolated.

[00240] Preferably an index may be up to 20 nucleotides in length, more preferably 1-10 nucleotides, and most preferably 4-6 nucleotides in length. A four nucleotide index gives a possibility of multiplexing 256 samples on the same array, a six base index enables 4096 samples to be processed on the same array. [00241] The precise nucleotide sequence of the universal adapters is generally not material to the disclosure and may be selected by the user such that the desired sequence elements are ultimately included in the common sequences of the plurality of different modified target nucleic acids, for example, to provide for the universal capture binding sequences for immobilizing the target nucleic acids on an array for subsequent sequencing, and binding sites for particular sets of universal amplification primers and/or sequencing primers. Additional sequence elements may be included, for example, to provide binding sites for sequencing primers which will ultimately be used in sequencing of target nucleic acids in the library, sequencing of an index, or products derived from amplification of the target nucleic acids in the library, for example on a solid support.

[00242] In order to prepare a library of deaminase-treated DNA for analysis using a sequencing platform, it may be useful to make additional modifications to the target DNA, either prior to or after treatment with cytidine deaminase. In some embodiments, single-stranded deaminase-treated DNA is prepared for sequencing using a single-stranded library preparation method, as is known in the art. Such methods include, but are not limited to, template switching based second strand synthesis, adapters containing a single-stranded splint overhang, and the like. Reagents for performing single-stranded library preparation methods are commercially available. Examples include xGen ssDNA & Low-Input DNA Library Prep Kit (Integrated DNA Technologies catalog number 10009859), previously sold as Accel-NGS (Swift Biosciences), NGS Single Stranded DNA Library Prep Kit (BioDynami catalog number 30082). Another example includes single-reaction singlestranded library (SRSLY) as set forth in Troll et al., BMC Genomics 20, 1023 (2019).

[00243] In some embodiments, library preparation modifications are made to double-stranded target DNA prior to treatment with cytidine deaminase. Methods for library preparation of double-stranded DNA template are known in the art, and include Y-adaptor ligation, transposome-based tagmentation, and the like. It will be appreciated by those of skill in the art that methods of double-strand library preparation often include one or more amplification steps using for example, PCR. In such methods, the amplification step may be deferred until after cytidine deaminase treatment, to preserve the methylation status of the template strand. For example, in Y-adapter ligation methods, the Y-adapters can be ligated to the double-stranded template, after which the adapter-ligated template DNA is denatured and treated with cytidine deaminase as described elsewhere herein. Following treatment with cytidine deaminase, the resulting treated single-strand DNA molecules can be amplified using PCR, bridge amplification, and other methods as are commonly known in the art.

[00244] H. Preparation of Immobilized Samples for Sequencing

[00245] The library of modified target nucleic acids, e.g., target nucleic acids having universal adapters at each end, can be prepared for sequencing. Methods for attaching modified target nucleic acids to a substrate are known in the art. In one embodiment, modified fragments are enriched using a plurality of capture oligonucleotides having specificity for the modified fragments, and the capture oligonucleotides can be immobilized on a surface of a solid substrate such as a flow cell or a bead. For instance, capture oligonucleotides can include a first member of a universal binding pair, and where a second member of the binding pair is immobilized on a surface of a solid substrate. Likewise, methods for amplifying immobilized target nucleic acids include, but are not limited to, bridge amplification and exclusion amplification (also referred to as kinetic exclusion amplification (KEA). Methods for immobilizing and amplifying prior to sequencing are described in, for instance, Bignell et al. (US 8,053,192), Gunderson et al. (W02016/130704), Shen et al. (US 8,895,249), and Pipenburg et al. (US 9,309,502).

[00246] A pooled sample can be immobilized in preparation for sequencing. Sequencing can be performed as an array of single molecules or can be amplified prior to sequencing. The amplification can be carried out using one or more immobilized primers. The immobilized primer(s) can be, for instance, a lawn on a planar surface, or on a pool of beads. The pool of beads can be isolated into an emulsion with a single bead in each "compartment" of the emulsion. At a concentration of only one template per "compartment," only a single template is amplified on each bead.

[00247] The term "solid-phase amplification" as used herein refers to any nucleic acid amplification reaction carried out on or in association with a solid support such that all or a portion of the amplified products are immobilized on the solid support as they are formed. In particular, the term encompasses solid-phase polymerase chain reaction (solid-phase PCR) and solid phase isothermal amplification which are reactions analogous to standard solution phase amplification, except that one or both of the forward and reverse amplification primers is/are immobilized on the solid support. Solid phase PCR covers systems such as emulsions, where one primer is anchored to a bead and the other is in free solution, and colony formation in solid phase gel matrices wherein one primer is anchored to the surface, and one is in free solution.

[00248] In some embodiments, the solid support comprises a patterned surface. A "patterned surface" refers to an arrangement of different regions in or on an exposed layer of a solid support. For example, one or more of the regions can be features where one or more amplification primers are present. The features can be separated by interstitial regions where amplification primers are not present. In some embodiments, the pattern can be an x- y format of features that are in rows and columns. In some embodiments, the pattern can be a repeating arrangement of features and/or interstitial regions. In some embodiments, the pattern can be a random arrangement of features and/or interstitial regions. Exemplary patterned surfaces that can be used in the methods and compositions set forth herein are described in U.S. Pat. Nos. 8,778,848, 8,778,849 and 9,079,148, and U.S. Pat. Appl. Pub. No. 2014/0243224.

[00249] In some embodiments, the solid support includes an array of wells or depressions in a surface. This may be fabricated as is generally known in the art using a variety of techniques, including, but not limited to, photolithography, stamping techniques, molding techniques and micro-etching techniques. As will be appreciated by those of skill in the art, the technique used will depend on the composition and shape of the array substrate.

[00250] The features in a patterned surface can be wells in an array of wells (e.g. microwells or nanowells) on glass, silicon, plastic or other suitable solid supports with patterned, covalently-linked gel such as poly(N-(5-azidoacetamidylpentyl)acrylamide-co-acrylamide) (PAZAM, see, for example, US Pub. No. 2013/184796, WO 2016/066586, and WO 2015/002813). The process creates gel pads used for sequencing that can be stable over sequencing runs with a large number of cycles. The covalent linking of the polymer to the wells is helpful for maintaining the gel in the structured features throughout the lifetime of the structured substrate during a variety of uses. However, in many embodiments the gel need not be covalently linked to the wells. For example, in some conditions silane free acrylamide (SFA, see, for example, US Pat. No. 8,563,477) which is not covalently attached to any part of the structured substrate, can be used as the gel material.

[00251] In particular embodiments, a structured substrate can be made by patterning a solid support material with wells (e.g. microwells or nanowells), coating the patterned support with a gel material (e.g. PAZAM, SFA, or chemically modified variants thereof, such as the azidolyzed version of SFA (azido-SFA)) and polishing the gel coated support, for example via chemical or mechanical polishing, thereby retaining gel in the wells but removing or inactivating substantially all of the gel from the interstitial regions on the surface of the structured substrate between the wells. Primer nucleic acids can be attached to gel material. A solution of modified target nucleic acids can then be contacted with the polished substrate such that individual modified target nucleic acids will seed individual wells via interactions with primers attached to the gel material; however, the target nucleic acids will not occupy the interstitial regions due to absence or inactivity of the gel material. Amplification of the modified target nucleic acids will be confined to the wells since absence or inactivity of gel in the interstitial regions prevents outward migration of the growing nucleic acid colony. The process can be conveniently manufactured, being scalable and utilizing conventional micro- or nanofabrication methods.

[00252] Although the disclosure encompasses "solid-phase" amplification methods in which only one amplification primer is immobilized (the other primer usually being present in free solution), in one embodiment the solid support is provided with both the forward and the reverse primers immobilized. In practice, there will be a plurality of identical forward primers and/or a plurality of identical reverse primers immobilized on the solid support, since the amplification process requires an excess of primers to sustain amplification. References herein to forward and reverse primers are to be interpreted accordingly as encompassing a plurality of such primers unless the context indicates otherwise. [00253] As will be appreciated by the skilled reader, any given amplification reaction requires at least one type of forward primer and at least one type of reverse primer specific for the template to be amplified. However, in certain embodiments the forward and reverse primers may include template-specific portions of identical sequence, and may have entirely identical nucleotide sequence and structure (including any non-nucleotide modifications). In other words, it is possible to carry out solid-phase amplification using only one type of primer, and such single-primer methods are encompassed within the scope of the disclosure. Other embodiments may use forward and reverse primers which contain identical template-specific sequences but which differ in some other structural features. For example, one type of primer may contain a non-nucleotide modification which is not present in the other.

[00254] Primers for solid-phase amplification are preferably immobilized by single point covalent attachment to the solid support at or near the 5' end of the primer, leaving the templatespecific portion of the primer free to anneal to its cognate template and the 3' hydroxyl group free for primer extension. Any suitable covalent attachment means known in the art may be used for this purpose. The chosen attachment chemistry will depend on the nature of the solid support, and any derivatization or functionalization applied to it. The primer itself may include a moiety, which may be a non-nucleotide chemical modification, to facilitate attachment. In a particular embodiment, the primer may include a sulphur- containing nucleophile, such as phosphorothioate or thiophosphate, at the 5' end. In the case of solid-supported polyacrylamide hydrogels, this nucleophile will bind to a bromoacetamide group present in the hydrogel. A more particular means of attaching primers and templates to a solid support is via 5' phosphorothioate attachment to a hydrogel comprised of polymerized acrylamide and N-(5-bromoacetamidylpentyl) acrylamide (BRAPA), as described in Int. Pub. No. WO 05/065814.

[00255] Certain embodiments of the disclosure may make use of solid supports that include an inert substrate or matrix (e.g. glass slides, polymer beads, etc.) which has been "functionalized," for example by application of a layer or coating of an intermediate material including reactive groups which permit covalent attachment to biomolecules, such as polynucleotides. Examples of such supports include, but are not limited to, polyacrylamide hydrogels supported on an inert substrate such as glass. In such embodiments, the biomolecules (e.g. polynucleotides) may be directly covalently attached to the intermediate material (e.g. the hydrogel), but the intermediate material may itself be non-covalently attached to the substrate or matrix (e.g. the glass substrate). The term "covalent attachment to a solid support" is to be interpreted accordingly as encompassing this type of arrangement.

[00256] The pooled samples may be amplified on beads wherein each bead contains a forward and reverse amplification primer. In one embodiment, a library of modified target nucleic acids is used to prepare clustered arrays of nucleic acid colonies, analogous to those described in U.S. Pub. No. 2005/0100900, U.S. Pat. No. 7,115,400, WO 00/18957 and WO 98/44151 by solid-phase amplification and more particularly solid phase isothermal amplification. The terms "cluster" and "colony" are used interchangeably herein to refer to a discrete site on a solid support including a plurality of identical immobilized nucleic acid strands and a plurality of identical immobilized complementary nucleic acid strands. The term "clustered array" refers to an array formed from such clusters or colonies. In this context, the term "array" is not to be understood as requiring an ordered arrangement of clusters.

[00257] The term "solid phase" or "surface" is used to mean either a planar array wherein primers are attached to a flat surface, for example, glass, silica or plastic microscope slides or similar flow cell devices; beads, wherein either one or two primers are attached to the beads and the beads are amplified; or an array of beads on a surface after the beads have been amplified.

[00258] Clustered arrays can be prepared using either a process of thermocycling, as described in WO 98/44151, or a process whereby the temperature is maintained as a constant, and the cycles of extension and denaturing are performed using changes of reagents. Such isothermal amplification methods are described in patent application numbers WO 02/46456 and U.S. Pub. No. 2008/0009420.

[00259] It will be appreciated that any of the amplification methodologies described herein or generally known in the art may be used with universal or target-specific primers to amplify immobilized DNA fragments. Suitable methods for amplification include, but are not limited to, the polymerase chain reaction (PCR), strand displacement amplification (SDA), transcription mediated amplification (TMA) and nucleic acid sequence-based amplification (NASBA), as described in U.S. Pat. No. 8,003,354. The above amplification methods may be employed to amplify one or more nucleic acids of interest. For example, PCR, including multiplex PCR, SDA, TMA, NASBA and the like may be utilized to amplify immobilized DNA fragments. In some embodiments, primers directed specifically to the polynucleotide of interest are included in the amplification reaction.

[00260] Other suitable methods for amplification of polynucleotides may include oligonucleotide extension and ligation, rolling circle amplification (RCA) (Lizardi et al., Nat. Genet. 19:225-232 (1998)) and oligonucleotide ligation assay (OLA) (See generally U.S. Pat. Nos. 7,582,420, 5,185,243, 5,679,524 and 5,573,907; EP 0 320 308 Bl; EP 0 336 731 Bl; EP 0 439 182 Bl; WO 90/01069; WO 89/12696; and WO 89/09835) technologies. It will be appreciated that these amplification methodologies may be designed to amplify immobilized DNA fragments. For example, in some embodiments, the amplification method may include ligation probe amplification or oligonucleotide ligation assay (OLA) reactions that contain primers directed specifically to the nucleic acid of interest. In some embodiments, the amplification method may include a primer extension-ligation reaction that contains primers directed specifically to the nucleic acid of interest. As a non-limiting example of primer extension and ligation primers that may be specifically designed to amplify a nucleic acid of interest, the amplification may include primers used for the Golden Gate assay (Illumina, Inc., San Diego, CA) as exemplified by U.S. Pat. No. 7,582,420 and 7,611,869.

[00261] DNA nanoballs can also be used in combination with methods, systems, compositions and kits as described herein. Methods for creating and using DNA nanoballs for genomic sequencing can be found at, for example, US patents and publications U.S. Pat. No. 7,910,354, 2009/0264299, 2009/0011943, 2009/0005252, 2009/0155781, 2009/0118488 and as described in, for example, Drmanac et al. (2010, Science 327(5961): 78-81). Briefly, following production of modified target nucleic acids, the modified target nucleic acids are circularized and amplified by rolling circle amplification (Lizardi et al., 1998. Nat. Genet. 19:225-232; US 2007/0099208 Al). The extended concatemeric structure of the amplicons promotes coiling creates compact DNA nanoballs. The DNA nanoballs can be captured on substrates, preferably to create an ordered or patterned array such that distance between each nanoball is maintained thereby allowing sequencing of the separate DNA nanoballs. In some embodiments such as those used by Complete Genomics (Mountain View, Calif.), consecutive rounds of adapter addition, amplification, and digestion are carried out prior to circularization to produce head to tail constructs having several target nucleic acids separated by adapter sequences.

[00262] Exemplary isothermal amplification methods that may be used in a method of the present disclosure include, but are not limited to, Multiple Displacement Amplification (MDA) as exemplified by, for example Dean et al., Proc. Natl. Acad. Sci. USA 99:5261-66 (2002) or isothermal strand displacement nucleic acid amplification exemplified by, for example U.S. Pat. No. 6,214,587. Other non-PCR-based methods that may be used in the present disclosure include, for example, strand displacement amplification (SDA) which is described in, for example Walker et al., Molecular Methods for Virus Detection, Academic Press, Inc., 1995; U.S. Pat. Nos. 5,455,166, and 5,130,238, and Walker et al., Nucl. Acids Res. 20:1691-96 (1992) or hyper-branched strand displacement amplification which is described in, for example Lage et al., Genome Res. 13:294-307 (2003). Isothermal amplification methods may be used with, for instance, the strand-displacing Phi 29 polymerase or Bst DNA polymerase large fragment, 5 '->3 ' exo- for random primer amplification of genomic DNA. The use of these polymerases takes advantage of their high processivity and strand displacing activity. High processivity allows the polymerases to produce fragments that are 10-20 kb in length. As set forth above, smaller fragments may be produced under isothermal conditions using polymerases having low processivity and strand-displacing activity such as Klenow polymerase. Additional description of amplification reactions, conditions and components are set forth in detail in the disclosure of U.S. Patent No. 7,670,810.

[00263] In some embodiments, isothermal amplification can be performed using kinetic exclusion amplification (KEA), also referred to as exclusion amplification (ExAmp). A nucleic acid library of the present disclosure can be made using a method that includes a step of reacting an amplification reagent to produce a plurality of amplification sites that each includes a substantially clonal population of amplicons from an individual target nucleic acid that has seeded the site. In some embodiments, the amplification reaction proceeds until a sufficient number of amplicons are generated to fill the capacity of the respective amplification site. Filling an already seeded site to capacity in this way inhibits target nucleic acids from landing and amplifying at the site thereby producing a clonal population of amplicons at the site. In some embodiments, apparent clonality can be achieved even if an amplification site is not filled to capacity prior to a second target nucleic acid arriving at the site. Under some conditions, amplification of a first target nucleic acid can proceed to a point that a sufficient number of copies are made to effectively outcompete or overwhelm production of copies from a second target nucleic acid that is transported to the site. For example, in an embodiment that uses a bridge amplification process on a circular feature that is smaller than 500 nm in diameter, it has been determined that after 14 cycles of exponential amplification for a first target nucleic acid, contamination from a second target nucleic acid at the same site will produce an insufficient number of contaminating amplicons to adversely impact sequencing-by-synthesis analysis on an Illumina sequencing platform.

[00264] In some embodiments, amplification sites in an array can be, but need not be, entirely clonal. Rather, for some applications, an individual amplification site can be predominantly populated with amplicons from a first modified target nucleic acid and can also have a low level of contaminating amplicons from a second modified target nucleic acid. An array can have one or more amplification sites that have a low level of contaminating amplicons so long as the level of contamination does not have an unacceptable impact on a subsequent use of the array. For example, when the array is to be used in a detection application, an acceptable level of contamination would be a level that does not impact signal to noise or resolution of the detection technique in an unacceptable way. Accordingly, apparent clonality will generally be relevant to a particular use or application of an array made by the methods set forth herein. Exemplary levels of contamination that can be acceptable at an individual amplification site for particular applications include, but are not limited to, at most 0.1%, 0.5%, 1%, 5%, 10% or 25% contaminating amplicons. An array can include one or more amplification sites having these exemplary levels of contaminating amplicons. For example, up to 5%, 10%, 25%, 50%, 75%, or even 100% of the amplification sites in an array can have some contaminating amplicons. It will be understood that in an array or other collection of sites, at least 50%, 75%, 80%, 85%, 90%, 95% or 99% or more of the sites can be clonal or apparently clonal.

[00265] In some embodiments, kinetic exclusion can occur when a process occurs at a sufficiently rapid rate to effectively exclude another event or process from occurring. Take, for example, the making of a nucleic acid array where sites of the array are randomly seeded with modified target nucleic acids from a solution and copies of the modified target nucleic acids are generated in an amplification process to fill each of the seeded sites to capacity. In accordance with the kinetic exclusion methods of the present disclosure, the seeding and amplification processes can proceed simultaneously under conditions where the amplification rate exceeds the seeding rate. As such, the relatively rapid rate at which copies are made at a site that has been seeded by a first target nucleic acid will effectively exclude a second nucleic acid from seeding the site for amplification. Kinetic exclusion amplification methods can be performed as described in detail in the disclosure of U.S. Pat. Appl. Pub. No. 2013/0338042.

[00266J Kinetic exclusion can exploit a relatively slow rate for initiating amplification (e.g. a slow rate of making a first copy of a modified target nucleic acids) vs. a relatively rapid rate for making subsequent copies of the modified target nucleic acids (or of the first copy of the modified target nucleic acids). In the example of the previous paragraph, kinetic exclusion occurs due to the relatively slow rate of modified target nucleic acids seeding (e.g. relatively slow diffusion or transport) vs. the relatively rapid rate at which amplification occurs to fill the site with copies of the modified target nucleic acid seed. In another exemplary embodiment, kinetic exclusion can occur due to a delay in the formation of a first copy of a modified target nucleic acid that has seeded a site (e.g. delayed or slow activation) vs. the relatively rapid rate at which subsequent copies are made to fill the site. In this example, an individual site may have been seeded with several different modified target nucleic acids (e.g. several modified target nucleic acids can be present at each site prior to amplification). However, first copy formation for any given modified target nucleic acid can be activated randomly such that the average rate of first copy formation is relatively slow compared to the rate at which subsequent copies are generated. In this case, although an individual site may have been seeded with several different modified target nucleic acids, kinetic exclusion will allow only one of those to be amplified. More specifically, once a first modified target nucleic acid has been activated for amplification, the site will rapidly fill to capacity with its copies, thereby preventing copies of a second modified target nucleic acid from being made at the site.

[00267] In one embodiment, the method is carried out to simultaneously (i) transport modified target nucleic acids to amplification sites at an average transport rate, and (ii) amplify the modified target nucleic acids that are at the amplification sites at an average amplification rate, wherein the average amplification rate exceeds the average transport rate (U.S. Pat. No. 9,169,513). Accordingly, kinetic exclusion can be achieved in such embodiments by using a relatively slow rate of transport. For example, a sufficiently low concentration of modified target nucleic acids can be selected to achieve a desired average transport rate, with lower concentrations resulting in slower average rates of transport. Alternatively or additionally, a high viscosity solution and/or presence of molecular crowding reagents in the solution can be used to reduce transport rates. Examples of useful molecular crowding reagents include, but are not limited to, polyethylene glycol (PEG), ficoll, dextran, or polyvinyl alcohol. Exemplary molecular crowding reagents and formulations are set forth in U.S. Pat. No. 7,399,590, which is incorporated herein by reference. Another factor that can be adjusted to achieve a desired transport rate is the average size of the target nucleic acids.

[00268] An amplification reagent can include further components that facilitate amplicon formation, and in some cases increase the rate of amplicon formation. An example is a recombinase. A recombinase can facilitate amplicon formation by allowing repeated invasion/extension. More specifically, a recombinase can facilitate invasion of a modified target nucleic acid by the polymerase and extension of a primer by the polymerase using the modified target nucleic acid as a template for amplicon formation. This process can be repeated as a chain reaction where amplicons produced from each round of invasion/extension serve as templates in a subsequent round. The process can occur more rapidly than standard PCR since a denaturation cycle (e.g. via heating or chemical denaturation) is not required. As such, recombinase-facilitated amplification can be carried out isothermally. It is generally desirable to include ATP, or other nucleotides (or in some cases non-hydrolyzable analogs thereof) in a recombinase-facilitated amplification reagent to facilitate amplification. A mixture of recombinase and single-stranded binding (SSB) protein is particularly useful as SSB can further facilitate amplification. Exemplary formulations for recombinase-facilitated amplification include those sold commercially as TwistAmp kits by TwistDx (Cambridge, UK). Useful components of recombinase- facilitated amplification reagent and reaction conditions are set forth in US 5,223,414 and US 7,399,590.

[00269] Another example of a component that can be included in an amplification reagent to facilitate amplicon formation and in some cases to increase the rate of amplicon formation is a helicase. A helicase can facilitate amplicon formation by allowing a chain reaction of amplicon formation. The process can occur more rapidly than standard PCR since a denaturation cycle (e.g. via heating or chemical denaturation) is not required. As such, helicase-facilitated amplification can be carried out isothermally. A mixture of helicase and single-stranded binding (SSB) protein is particularly useful as SSB can further facilitate amplification. Exemplary formulations for helicase-facilitated amplification include those sold commercially as IsoAmp kits from Biohelix (Beverly, MA). Further, examples of useful formulations that include a helicase protein are described in US 7,399,590 and US 7,829,284.

[00270] Yet another example of a component that can be included in an amplification reagent to facilitate amplicon formation and in some cases increase the rate of amplicon formation is an origin binding protein.

[00271] I. Methods of Sequencing

[00272] Following attachment of modified target nucleic acids to a surface, the sequence of the immobilized and amplified modified target nucleic acids is determined. Sequencing can be carried out using any suitable sequencing technique, and methods for determining the sequence of immobilized and amplified modified target nucleic acids, including strand resynthesis, are known in the art and are described in, for instance, Bignell et al. (US 8,053,192), Gunderson et al. (W02016/130704), Shen et al. (US 8,895,249), and Pipenburg et al. (US 9,309,502). [00273] The methods described herein can be used in conjunction with a variety of nucleic acid sequencing techniques. Particularly applicable techniques are those wherein nucleic acids are attached at fixed locations in an array such that their relative positions do not change and wherein the array is repeatedly imaged. Embodiments in which images are obtained in different color channels, for example, coinciding with different labels used to distinguish one nucleotide base type from another are particularly applicable. In some embodiments, the process to determine the nucleotide sequence of a modified target nucleic acid can be an automated process. Preferred embodiments include sequencing-by-synthesis ("SBS") techniques.

[00274] SBS techniques generally involve the enzymatic extension of a nascent nucleic acid strand through the iterative addition of nucleotides against a template strand. In traditional methods of SBS, a single nucleotide monomer may be provided to a target nucleotide in the presence of a polymerase in each delivery. However, in the methods described herein, more than one type of nucleotide monomer can be provided to a target nucleic acid in the presence of a polymerase in a delivery.

[00275] In one embodiment, a nucleotide monomer includes locked nucleic acids (LNAs) or bridged nucleic acids (BNAs). The use of LNAs or BNAs in a nucleotide monomer increases hybridization strength between a nucleotide monomer and a sequencing primer sequence present on an immobilized modified target nucleic acid.

[00276] SBS techniques can use nucleotide monomers that have a terminator moiety or those that lack any terminator moieties. Methods using nucleotide monomers lacking terminators include, for example, pyrosequencing and sequencing using y-phosphate-labeled nucleotides, as set forth in further detail herein. In methods using nucleotide monomers lacking terminators, the number of nucleotides added in each cycle is generally variable and dependent upon the template sequence and the mode of nucleotide delivery. For SBS techniques that use nucleotide monomers having a terminator moiety, the terminator can be effectively irreversible under the sequencing conditions used as is the case for traditional Sanger sequencing which utilizes dideoxynucleotides, or the terminator can be reversible as is the case for sequencing methods developed by Solexa (now Illumina, Inc.). [00277] SBS techniques can use nucleotide monomers that have a labeled moiety or those that lack a labeled moiety. Accordingly, incorporation events can be detected based on a characteristic of the label, such as fluorescence of the label; a characteristic of the nucleotide monomer such as molecular weight or charge; a byproduct of incorporation of the nucleotide, such as release of pyrophosphate; or the like. In embodiments where two or more different nucleotides are present in a sequencing reagent, the different nucleotides can be distinguishable from each other, or alternatively the two or more different labels can be indistinguishable under the detection techniques being used. For example, the different nucleotides present in a sequencing reagent can have different labels and they can be distinguished using appropriate optics as exemplified by the sequencing methods developed by Solexa (now Illumina, Inc.).

[00278] Preferred embodiments include pyrosequencing techniques. Pyrosequencing detects the release of inorganic pyrophosphate (PPi) as particular nucleotides are incorporated into the nascent strand (Ronaghi, M., Karamohamed, S., Pettersson, B., Uhlen, M. and Nyren, P. (1996) "Real-time DNA sequencing using detection of pyrophosphate release." Analytical Biochemistry 242(1), 84-9; Ronaghi, M. (2001) "Pyrosequencing sheds light on DNA sequencing." Genome Res. 11(1), 3-11; Ronaghi, M., Uhlen, M. and Nyren, P. (1998) "A sequencing method based on real-time pyrophosphate." Science 281(5375), 363; U.S. Pat. Nos. 6,210,891 ; 6,258,568 and 6,274,320). In pyrosequencing, released PPi can be detected by being immediately converted to adenosine triphosphate (ATP) by ATP sulfurase, and the level of ATP generated can be detected via luciferase-produced photons. The nucleic acids to be sequenced can be attached to features in an array and the array can be imaged to capture the chemiluminescent signals that are produced due to incorporation of a nucleotides at the features of the array. An image can be obtained after the array is treated with a particular nucleotide type (e.g. A, T, C or G). Images obtained after addition of each nucleotide type will differ with regard to which features in the array are detected. These differences in the image reflect the different sequence content of the features on the array. However, the relative locations of each feature will remain unchanged in the images. The images can be stored, processed and analyzed using the methods set forth herein. For example, images obtained after treatment of the array with each different nucleotide type can be handled in the same way as exemplified herein for images obtained from different detection channels for reversible terminator-based sequencing methods.

[00279] In another exemplary type of SBS, cycle sequencing is accomplished by stepwise addition of reversible terminator nucleotides containing, for example, a cleavable or photobleachable dye label as described, for example, in WO 04/018497 and U.S. Pat. No. 7,057,026. This approach is being commercialized by Solexa (now Illumina Inc.), and is also described in WO 91/06678 and WO 07/123,744. The availability of fluorescently- labeled terminators in which both the termination can be reversed and the fluorescent label cleaved facilitates efficient cyclic reversible termination (CRT) sequencing. Polymerases can also be co-engineered to efficiently incorporate and extend from these modified nucleotides.

[00280] In some reversible terminator-based sequencing embodiments, the labels do not substantially inhibit extension under SBS reaction conditions. However, the detection labels can be removable, for example, by cleavage or degradation. Images can be captured following incorporation of labels into arrayed nucleic acid features. In particular embodiments, each cycle involves simultaneous delivery of four different nucleotide types to the array and each nucleotide type has a spectrally distinct label. Four images can then be obtained, each using a detection channel that is selective for one of the four different labels. Alternatively, different nucleotide types can be added sequentially and an image of the array can be obtained between each addition step. In such embodiments, each image will show nucleic acid features that have incorporated nucleotides of a particular type. Different features will be present or absent in the different images due the different sequence content of each feature. However, the relative position of the features will remain unchanged in the images. Images obtained from such reversible terminator- SBS methods can be stored, processed and analyzed as set forth herein. Following the image capture step, labels can be removed and reversible terminator moieties can be removed for subsequent cycles of nucleotide addition and detection. Removal of the labels after they have been detected in a particular cycle and prior to a subsequent cycle can provide the advantage of reducing background signal and crosstalk between cycles. Examples of useful labels and removal methods are set forth herein. [00281] In particular embodiments some or all of the nucleotide monomers can include reversible terminators. In such embodiments, reversible terminators/cleavable fluorophores can include fluorophores linked to the ribose moiety via a 3' ester linkage (Metzker, Genome Res. 15: 1767-1776 (2005)). Other approaches have separated the terminator chemistry from the cleavage of the fluorescent label (Ruparel et al., Proc Natl Acad Sci USA 102: 5932-7 (2005)). Ruparel et al. described the development of reversible terminators that used a small 3' allyl group to block extension, but could easily be deblocked by a short treatment with a palladium catalyst. The fluorophore was attached to the base via a photocl eavable linker that could easily be cleaved by a 30 second exposure to long wavelength UV light. Thus, either disulfide reduction or photocleavage can be used to cleave a linker. Another approach to reversible termination is the use of natural termination that ensues after placement of a bulky dye on a dNTP. The presence of a charged bulky dye on the dNTP can act as an effective terminator through steric and/or electrostatic hindrance. The presence of one incorporation event prevents further incorporations unless the dye is removed. Cleavage of the dye removes the fluorophore and effectively reverses the termination. Examples of modified nucleotides are also described in U.S. Pat. Nos. 7,427,673, and 7,057,026.

[00282] Additional exemplary SBS systems and methods which can be used with the methods and systems described herein are described in U.S. Patent Application Publication Nos. 2007/0166705, 2006/0188901, 2006/0240439, 2006/0281109, 2012/0270305, and 2013/0260372, U.S. Patent No. 7,057,026, PCT Publication No. WO 05/065814, U.S. Patent Application Publication No. 2005/0100900, and PCT Publication Nos. WO 06/064199 and WO 07/010,251.

[00283] Some embodiments can use detection of four different nucleotides using fewer than four different labels. For example, SBS can be performed using methods and systems described in the incorporated materials of U.S. Pub. No. 2013/0079232. As a first example, a pair of nucleotide types can be detected at the same wavelength, but distinguished based on a difference in intensity for one member of the pair compared to the other, or based on a change to one member of the pair (e.g., via chemical modification, photochemical modification or physical modification) that causes apparent signal to appear or disappear compared to the signal detected for the other member of the pair. As a second example, three of four different nucleotide types can be detected under particular conditions while a fourth nucleotide type lacks a label that is detectable under those conditions, or is minimally detected under those conditions (e.g., minimal detection due to background fluorescence, etc.). Incorporation of the first three nucleotide types into a nucleic acid can be determined based on presence of their respective signals and incorporation of the fourth nucleotide type into the nucleic acid can be determined based on absence or minimal detection of any signal. As a third example, one nucleotide type can include label(s) that are detected in two different channels, whereas other nucleotide types are detected In no more than one of the channels. The aforementioned three exemplary configurations are not considered mutually exclusive and can be used in various combinations. An exemplary embodiment that combines all three examples, is a fluorescent-based SBS method that uses a first nucleotide type that is detected in a first channel (e.g. dATP having a label that is detected in the first channel when excited by a first excitation wavelength), a second nucleotide type that is detected in a second channel (e.g. dCTP having a label that is detected in the second channel when excited by a second excitation wavelength), a third nucleotide type that is detected in both the first and the second channel (e.g. dTTP having at least one label that is detected in both channels when excited by the first and/or second excitation wavelength) and a fourth nucleotide type that lacks a label that is not, or minimally, detected in either channel (e.g. dGTP having no label).

[00284] Further, as described in U.S. Pub. No. 2013/0079232, sequencing data can be obtained using a single channel. In such so-called one-dye sequencing approaches, the first nucleotide type is labeled but the label is removed after the first image is generated, and the second nucleotide type is labeled only after a first image is generated. The third nucleotide type retains its label in both the first and second images, and the fourth nucleotide type remains unlabeled in both images.

[00285] Some embodiments can use sequencing by ligation techniques. Such techniques use DNA ligase to incorporate oligonucleotides and identify the incorporation of such oligonucleotides. The oligonucleotides typically have different labels that are correlated with the identity of a particular nucleotide in a sequence to which the oligonucleotides hybridize. As with other SBS methods, images can be obtained following treatment of an array of nucleic acid features with the labeled sequencing reagents. Each image will show nucleic acid features that have incorporated labels of a particular type. Different features will be present or absent in the different images due the different sequence content of each feature, but the relative position of the features will remain unchanged in the images. Images obtained from ligation-based sequencing methods can be stored, processed and analyzed as set forth herein. Exemplary SBS systems and methods which can be utilized with the methods and systems described herein are described in U.S. Pat. Nos. 6,969,488, 6,172,218, and 6,306,597.

[00286] Some embodiments can use nanopore sequencing (Deamer, D. W. & Akeson, M. "Nanopores and nucleic acids: prospects for ultrarapid sequencing." Trends Biotechnol. 18, 147-151 (2000); Deamer, D. and D. Branton, "Characterization of nucleic acids by nanopore analysis", Acc. Chem. Res. 35:817-825 (2002); Li, J., M. Gershow, D. Stein, E. Brandin, and J. A. Golovchenko, "DNA molecules and configurations in a solid-state nanopore microscope" Nat. Mater. 2:611-615 (2003)). In such embodiments, the modified target nucleic acid passes through a nanopore. The nanopore can be a synthetic pore or biological membrane protein, such as a-hemolysin. As the modified target nucleic acid passes through the nanopore, each base-pair can be identified by measuring fluctuations in the electrical conductance of the pore. (U.S. Pat. No. 7,001,792; Soni, G. V. & Meller, "A. Progress toward ultrafast DNA sequencing using solid-state nanopores." Clin. Chem. 53, 1996-2001 (2007); Healy, K. "Nanopore-based single-molecule DNA analysis." Nanomed. 2, 459-481 (2007); Cockroft, S. L., Chu, J., Amorin, M. & Ghadiri, M. R. "A singlemolecule nanopore device detects DNA polymerase activity with single-nucleotide resolution." J. Am. Chem. Soc. 130, 818-820 (2008)). Data obtained from nanopore sequencing can be stored, processed and analyzed as set forth herein. In particular, the data can be treated as an image in accordance with the exemplary treatment of optical images and other images that is set forth herein.

[00287] Some embodiments can use methods involving the real-time monitoring of DNA polymerase activity. Nucleotide incorporations can be detected through fluorescence resonance energy transfer (FRET) interactions between a fluorophore-bearing polymerase and y-phosphate-labeled nucleotides as described, for example, in U.S. Pat. Nos. 7,329,492 and 7,211,414, or nucleotide incorporations can be detected with zero-mode waveguides as described, for example, in U.S. Pat. No. 7,315,019, and using fluorescent nucleotide analogs and engineered polymerases as described, for example, in U.S. Pat. No. 7,405,281 and U.S. Pub. No. 2008/0108082. The illumination can be restricted to a zeptoliter-scale volume around a surface-tethered polymerase such that incorporation of fluorescently labeled nucleotides can be observed with low background (Levene, M. J. et al. "Zero-mode waveguides for single-molecule analysis at high concentrations." Science 299, 682-686 (2003); Lundquist, P. M. et al. "Parallel confocal detection of single molecules in real time." Opt. Lett. 33, 1026-1028 (2008); Korlach, J. et al. "Selective aluminum passivation for targeted immobilization of single DNA polymerase molecules in zero-mode waveguide nano structures." Proc. Natl. Acad. Sci. USA 105, 1176-1181 (2008)). Images obtained from such methods can be stored, processed and analyzed as set forth herein.

[00288] Some SBS embodiments include detection of a proton released upon incorporation of a nucleotide into an extension product. For example, sequencing based on detection of released protons can use an electrical detector and associated techniques that are commercially available from Ion Torrent (Guilford, CT, a Life Technologies subsidiary) or sequencing methods and systems described in U.S. Pub. Nos. 2009/0026082;

2009/0127589; 2010/0137143; and 2010/0282617. Methods set forth herein for amplifying target nucleic acids using kinetic exclusion can be readily applied to substrates used for detecting protons. More specifically, methods set forth herein can be used to produce clonal populations of amplicons that are used to detect protons.

[00289] The above SBS methods can be advantageously carried out in multiplex formats such that multiple different modified target nucleic acids are manipulated simultaneously. In particular embodiments, different modified target nucleic acids can be treated in a common reaction vessel or on a surface of a particular substrate. This allows convenient delivery of sequencing reagents, removal of unreacted reagents and detection of incorporation events in a multiplex manner. In embodiments using surface-bound target nucleic acids, the modified target nucleic acids can be in an array format. In an array format, the modified target nucleic acids can be typically bound to a surface in a spatially distinguishable manner. The modified target nucleic acids can be bound by direct covalent attachment, attachment to a bead or other particle or binding to a polymerase or other molecule that is attached to the surface. The array can include a single copy of a modified target nucleic acid at each site (also referred to as a feature) or multiple copies having the same sequence can be present at each site or feature. Multiple copies can be produced by amplification methods such as, bridge amplification or emulsion PCR as described in further detail herein.

[00290] The methods set forth herein can use arrays having features at any of a variety of densities including, for example, at least about 10 features/cm ² , 100 features/ cm ² , 500 features/ cm ² , 1,000 features/ cm ² , 5,000 features/ cm ² , 10,000 features/ cm ² , 50,000 features/ cm ² , 100,000 features/ cm ² , 1,000,000 features/ cm ² , 5,000,000 features/ cm ² , or higher.

[00291] An advantage of the methods set forth herein is that they provide for rapid and efficient detection of a plurality of cm ² , in parallel. Accordingly, the present disclosure provides integrated systems capable of preparing and detecting nucleic acids using techniques known in the art such as those exemplified herein. Thus, an integrated system of the present disclosure can include fluidic components capable of delivering amplification reagents and/or sequencing reagents to one or more immobilized modified target nucleic acids, the system including components such as pumps, valves, reservoirs, fluidic lines and the like. A flow cell can be configured and/or used in an integrated system for detection of target nucleic acids. Exemplary flow cells are described, for example, in US Pat. No. 8,241,573 and US Pat. No. 8,951,781. As exemplified for flow cells, one or more of the fluidic components of an integrated system can be used for an amplification method and for a detection method. Taking a nucleic acid sequencing embodiment as an example, one or more of the fluidic components of an integrated system can be used for an amplification method set forth herein and for the delivery of sequencing reagents in a sequencing method such as those exemplified above. Alternatively, an integrated system can include separate fluidic systems to carry out amplification methods and to carry out detection methods. Examples of integrated sequencing systems that are capable of creating amplified nucleic acids and also determining the sequence of the nucleic acids include, without limitation, the MiSeqTM platform (Illumina, Inc., San Diego, CA) and devices described in US Pat. No. 8,951,781.

[00292] While the embodiments presented herein are generally described using a sequencing platform (such as a sequencing by synthesis platform) as a readout, one of ordinary skill in the art will recognize that nucleic acids modified by the cytidine deaminases, TET proteins, glucosyltransferase enzymes, and/or borane derivatives described herein can also be detected using any other suitable readout methodology. For example, the location and identity of modified cytosines can be assessed using a microarray. Any of a variety of analyte arrays (also referred to as “microarrays”) known in the art can be used in a method or system set forth herein. A typical array contains analytes, each having an individual probe or a population of probes. In the latter case, the population of probes at each analyte is typically homogenous having a single species of probe. For example, in the case of a nucleic acid array, each analyte can have multiple nucleic acid molecules each having a common sequence. However, in some implementations the populations at each analyte of an array can be heterogeneous. Similarly, protein arrays can have analytes with a single protein or a population of proteins typically, but not always, having the same amino acid sequence. The probes can be attached to the surface of an array for example, via covalent linkage of the probes to the surface or via non-covalent interaction(s) of the probes with the surface. In some implementations, probes, such as nucleic acid molecules, can be attached to a surface via a gel layer as described, for example, in U.S. patent application Ser. No. 13/784,368 and US Pat. App. Pub. No. 2011/0059865 Al.

[00293] Example arrays include, without limitation, a BeadChip Array available from Illumina, Inc. (San Diego, Calif.) or others such as those where probes are attached to beads that are present on a surface (e.g. beads in wells on a surface) such as those described in U.S. Pat. Nos. 6,266,459; 6,355,431; 6,770,441; 6,859,570; or 7,622,294; or PCT Publication No. WO 00/63437. Further examples of commercially available microarrays that can be used include, for example, an Affymetrix® GeneChip® microarray or other microarray synthesized in accordance with techniques sometimes referred to as VLSIPS™ (Very Large Scale Immobilized Polymer Synthesis) technologies. A spotted microarray can also be used in a method or system according to some implementations of the present disclosure. An example spotted microarray is a CodeLink™ Array available from Amersham Biosciences. Another microarray that is useful is one that is manufactured using inkjet printing methods such as SurePrint™ Technology available from Agilent Technologies.

[00294] In a specific embodiment, nucleic acids in a sample can be treated with a CpG binding protein, and a cytidine deaminase as presented herein can be used to convert 5-methyl cytosine (5mC) to thymidine (T) by deamination as described herein, such as by providing a sample of DNA suspected of including single-stranded DNA including at least one 5- methyl cytosine (5mC), at least one 5 -hydroxymethyl cytosine (5hmC), at least one 5- formyl cytosine (5fC), at least one 5-carboxy cytosine (5CaC), or a combination thereof; contacting the DNA with CpG binding protein and the cytidine deaminase under conditions suitable for conversion of 5 methylcytosine (5mC) to thymidine (T) by deamination at a greater rate than conversion of cytosine (C) to uracil (U) by deamination, to result in converted single-stranded DNA, wherein 5mC, 5hmC, 5fC, and/or 5CaC are converted to T.

[00295] In a specific embodiment, an altered cytidine deaminase of the present disclosure can be used to detect 5hmC as described herein, such as by providing a sample of DNA suspected of including single-stranded DNA that has at least one 5 -hydroxymethyl cytosine (5hmC); contacting the DNA with the altered cytidine deaminase under conditions suitable for conversion of unmodified cytosine to uracil and 5mC to thymidine and no detectable conversion of 5hmC to 5hmU.

[00296] In another specific embodiment, nucleic acids in a sample can be treated with a CpG binding protein, a TET protein, and a borane derivative, and optionally a BGT protein can be used to convert 5-methyl cytosine (5mC) to 5-carboxylcytosine (5caC) as described herein, such as by providing a sample of DNA suspected of including DNA including at least one 5-methyl cytosine (5mC), at least one 5 -hydroxymethyl cytosine (5hmC), at least one 5-formyl cytosine (5fC), at least one 5-carboxy cytosine (5CaC), or a combination thereof; contacting the DNA with the CpG binding protein, TET protein, and a borane derivative, and optionally BGT protein under suitable conditions to result in converted DNA.

[00297] The converted DNA can then be processed as needed to facilitate hybridization to a microarray. For example, the converted DNA can be amplified. Any one of a number of amplification methods as are known in the art can be performed. For example, wholegenome amplification or amplification using universal primers that hybridize to a common region in the converted DNA, such as an adaptor sequence, can be used. Additionally or alternatively, the converted DNA can be fragmented. Fragmentation can be performed prior to or following amplification, or in the absence of amplification. Any one of a number of fragmentation methods as are known in the art can be performed. As one example, fragmentation can be performed using an enzymatic process, such as a restriction endonuclease or other enzyme capable of cleaving the converted DNA. As another example, fragmentation can be performed using mechanical means, such as shearing using, for example a sonication device such as those supplied by Covaris. The fragmented converted DNA can then be precipitated and/or resuspended in a buffer suitable for hybridization to a microarray. Following hybridization, the methylation state of regions of interest, such as a specific CpG locus or loci, can be interrogated at specific locations on the microarray. Methods of preparing converted DNA for microarray analysis are known in the art. One example of such methods is described in the Methylation Protocol Guide for the Infinium HD Assay from Illumina (San Diego, CA). Whereas such a protocol guide may describe use of a microarray designed for interrogation of bisulfite-converted DNA, it will be understood that array features, specifically probe sequences, can be specifically designed for DNA that has not been bisulfite converted. As an example, a commercially available microarray such as the Infinium Methyl ationEPIC BeadChip (Illumina) is specifically designed to hybridize with DNA fragments with reduced complexity, as found in bisulfite converted DNA, where most if not all cytosines are converted to thymidine. Thus, for example, the same CpG sites can be interrogated in non-bisulfite-converted DNA by using a microarray including probes designed to hybridize to the same regions of native, non-bisulfite-converted DNA. One of skill in the art could readily obtain such a microarray. In one embodiment, a custom array could be designed using the manifest for an array such as the Infinium Methyl ationEPIC BeadChip, by using the “Forward Sequence” to identify a probe sequence including native DNA sequence that covers a similar or identical sequence region for the allele-specific probe sequences, which are designed to hybridize to DNA sequences where most or all cytosines have been converted to thymidine. Using such an array designed to hybridize to native (non-bisulfite-converted) DNA sequences, the methodologies and analysis methods described in the Methylation Protocol Guide for the Infinium HD Assay from Illumina (San Diego, CA) could be followed to identify methylated CpG sites in the sample DNA.

[00298] J. Compositions

[00299] The present disclosure also provides compositions that include a CpG binding protein as described herein. In one embodiment a composition includes a CpG binding protein and a cytidine deaminase, such as an altered cytidine deaminase as described herein. In another embodiment a composition includes a CpG binding protein and a TET protein. In another embodiment a composition includes a CpG binding protein, a TET protein, and a cytidine deaminase. The composition can include one or more additional other components. For example, the other component can include a single-stranded or double-stranded DNA or RNA substrate that includes, or is suspected of including, at least one modified cytosine, such as a 5-methyl cytosine, a 5 -hydroxymethyl cytosine, a 5-formyl cytosine (5fC), a 5- carboxy cytosine (5CaC), or a combination thereof. In another example, a single-stranded or double-stranded DNA or RNA substrate can be one including one or more known modified cytosine, e.g., a single-stranded DNA or RNA substrate that can be used as a control to measure conversion efficiency. In another example, the other component can include a buffer having a pH that is described herein. In another example, the other components can include a buffer described herein, such as a citrate buffer, a sodium acetate buffer, or a Bis-Tris buffer. In another example, the other component can include a reductant, including but not limited to, DTT and/or TCEP, as well as Zn. In one embodiment, a composition includes a CpG binding protein that is a fusion protein as described herein.

[00300] A composition can also include a polynucleotide encoding an CpG binding protein, a TET protein, a cytidine deaminase, or any combination thereof. The polynucleotide can be present in a vector, such as a plasmid or virus vector. A vector that includes the polynucleotide can be present in a host cell, such as E. coli.

[00301] K Kits

[00302J The present disclosure also provides kits for determining the methylation status of DNA or RNA. In one embodiment, a kit includes at least CpG binding protein and at least one cytidine deaminase described herein and one or more other components in a suitable packaging material in an amount sufficient for at least one reaction. In another embodiment, a kit includes at least CpG binding protein and at least one TET protein as described herein, optionally a borane derivative, and one or more other components in a suitable packaging material in an amount sufficient for at least one reaction. Examples of other components include a positive control polynucleotide, such as a single-stranded DNA or a double-stranded DNA including one or more known modified cytosines for use in measuring conversion efficiency, or a negative control polynucleotide, such as a singlestranded DNA a double-stranded DNA including unmodified cytosines. Another component can be a glucosyltransferase, such as T4-beta glucosyltransferase. Optionally, other reagents such as buffers and solutions needed to use the altered cytidine deaminase and nucleotide solution are also included. Instructions for use of the packaged components are also typically included.

[00303] As used herein, the phrase "packaging material" refers to one or more physical structures used to house the contents of the kit. The packaging material is constructed by known methods, preferably to provide a sterile, contaminant-free environment. The packaging material has a label which indicates that the components can be used for determining the methylation status of DNA or RNA. In addition, the packaging material contains instructions indicating how the materials within the kit are employed to practice a reaction with a CpG binding protein and a cytidine deaminase or a TET protein. As used herein, the term "package" refers to a solid matrix or material such as glass, plastic, paper, foil, and the like, capable of holding within fixed limits the polypeptides. "Instructions for use" typically include a tangible expression describing the reagent concentration or at least one assay method parameter, such as the relative amounts of reagent and sample to be admixed, maintenance time periods for reagent/sample admixtures, temperature, buffer conditions, and the like.

[00304] The invention is defined in the claims. However, below there is provided a non-exhaustive listing of non-limiting exemplary aspects. Any one or more of the features of these aspects may be combined with any one or more features of another example, embodiment, or aspect described herein.

[00305] Exemplary Aspects

[00306] Aspect 1 is a composition comprising a cytidine deaminase and a CpC binding protein.

[00307] Aspect 2 is the composition of aspect 1, wherein the composition further comprises singlestranded DNA, double-stranded DNA, or both single-stranded DNA and double-stranded DNA.

[00308] Aspect 3 is the composition of aspect 1, wherein the sample comprises genomic DNA or cell free DNA.

[00309] Aspect 4 is the composition of aspect 3, wherein the genomic DNA is from a single cell or is a mixture from a plurality of cells.

[00310] Aspect 5 is the composition of any preceding aspect, wherein the composition further comprises cells or isolated nuclei.

[00311] Aspect 6 is the composition of any preceding aspect, wherein the cytidine deaminase is a member of the AID subfamily, the APOB EC 1 subfamily, the APOBEC2 subfamily, the APOBEC3A subfamily, the APOBEC3B subfamily, the APOBEC3C subfamily, the APOBEC3D subfamily, the APOBEC3F subfamily, the APOBEC3G subfamily, the APOBEC3G subfamily, the APOBEC3H subfamily, or the APOBEC4 subfamily.

[00312] Aspect 7 is the composition of any preceding aspect, wherein the cytidine deaminase is an altered cytidine deaminase comprising an amino acid substitution mutation at a position functionally equivalent to (Tyr/Phe)130 in a wild-type APOBEC3A protein [00313] Aspect 8 is the composition of any preceding aspect, wherein the wherein the (Tyr/Phe)130 is Tyrl30, and the wild-type APOBC3A protein is SEQ ID NO:3.

[00314] Aspect 9 is the composition of any preceding aspect, wherein the substitution mutation at the position functionally equivalent to Tyrl30 comprises a mutation to Ala, Vai, or Trp.

[00315] Aspect 10 is the composition of any preceding aspect, wherein the altered cytidine deaminase comprises amino acid substitution mutations at positions functionally equivalent to (Tyr/Phe)130 and Tyrl32 in a wild-type APOBEC3A protein.

[00316] Aspect 11 is the composition of any preceding aspect, wherein the substitution mutation at the position functionally equivalent to (Tyr/Phe)130 comprises a mutation to Ala, Vai, or Trp and the substitution mutation at the position functionally equivalent to Tyrl32 comprises a mutation to Arg, His, Leu, or Gin, or a combination thereof..

[00317] Aspect 12 is the composition of any preceding aspect, wherein the altered cytidine deaminase is derived from a member of the AID subfamily, the APOBEC1 subfamily, the APOBEC2 subfamily, the APOBEC3A subfamily, the APOBEC3B subfamily, the APOBEC3C subfamily, the APOBEC3D subfamily, the APOBEC3F subfamily, the APOBEC3G subfamily, the APOBEC3G subfamily, the APOBEC3H subfamily, or the APOBEC4 subfamily.

[00318] Aspect 13 is the composition of any preceding aspect, wherein the cytidine deaminase or the altered cytidine deaminase comprises a ZDD motif H- [P/A/V]-E-X[23-28]-P-C-X[2- 4]-C (SEQ ID NO: 12).

[00319] Aspect 14 is the composition of any preceding aspect, wherein the cytidine deaminase or the altered cytidine deaminase is a member of the APOBEC3A subfamily and comprises a ZDD motif HXEX24SW(S/T)PCX[2-4]CX6FX8LX5R(L/I)YX[8-1 l]LX2LX[10]M (SEQ ID N0: 13).

[00320] Aspect 15 is the composition of any preceding aspect, wherein the CpC binding protein is a ID AX CXXC domain protein. [00321] Aspect 16 is the composition of any preceding aspect, wherein the CpC binding protein comprises at least 80% identity to SEQ ID NO:57 or SEQ ID NO:58.

[00322] Aspect 17 is the composition of any preceding aspect, wherein the CpG binding protein is a fusion protein.

[00323] Aspect 18 is the composition of any preceding aspect, wherein the fusion protein comprises a domain comprising the activity of binding to a single-stranded nucleic acid.

[00324] Aspect 19 is the composition of any preceding aspect, wherein the fusion protein comprises a domain comprising the activity of increasing the solubility of the CpG binding protein.

[00325] Aspect 20 is a composition comprising a means for converting a 5 methylcytosine to a thymine and a means for binding to unmethylated CpG dinucleotides.

[00326] Aspect 21 is a method for preparing a sequencing library suitable for identifying methylated nucleotides comprising providing a sample of DNA suspected of comprising single-stranded or double-stranded DNA comprising at least one 5-methyl cytosine (5mC); contacting the single-stranded or double-stranded DNA with a CpG binding protein to result in a mixture, wherein the CpG binding protein binds unmethylated CpG dinucleotides present in the single-stranded or double-stranded DNA; adding to the mixture a cytidine deaminase or an altered cytidine deaminase, processing the converted singlestranded or double-stranded DNA to produce a sequencing library.

[00327] Aspect 22 is the method of any preceding aspect, further comprising, before the contacting, converting double-stranded DNA to single-stranded DNA.

[00328] Aspect 23 is the method of any preceding aspect, wherein the contacting and the adding occur at simultaneously.

[00329] Aspect 24 is a method for preparing a sequencing library suitable for identifying methylated nucleotides comprising providing a sample of DNA suspected of comprising double-stranded DNA comprising at least one 5-methyl cytosine (5mC); processing the double-stranded DNA to produce a sequencing library; denaturing the sequencing library to result in single-stranded DNA; contacting the single- stranded DNA with a CpG binding protein to result in a mixture, wherein the CpG binding protein binds unmethylated CpG dinucleotides present in the single-stranded DNA; contacting the mixture with a cytidine deaminase or an altered cytidine deaminase, converting the converted single-stranded DNA to a converted double-stranded DNA sequencing library.

[00330] Aspect 25 is the method of any preceding aspect, wherein the CpG binding protein is an ID AX CXXC domain.

[00331] Aspect 26 is the method of any preceding aspect, wherein the CpG binding protein is a fusion protein.

[00332] Aspect 27 is the method of any preceding aspect, wherein the CpC binding protein comprises at least 80% identity to SEQ ID NO:57 or SEQ ID NO:58.

[00333] Aspect 28 is the method of any preceding aspect, wherein the sample is a biological sample.

[00334] Aspect 29 is the method of any preceding aspect, wherein the biological sample comprises cell-free DNA.

[00335] Aspect 30 is the method of any preceding aspect, wherein the biological sample comprises a fluid selected from blood or serum.

[00336] Aspect 31 is the method of any preceding aspect, wherein the sample comprises single cells or isolated nuclei.

[00337] Aspect 32 is the method of any preceding aspect, wherein the biological sample comprises a tissue.

[00338] Aspect 33 is the method of any preceding aspect, wherein the tissue comprises tumor tissue.

[00339] Aspect 34 is the method of any preceding aspect, wherein the sample comprises genomic DNA. [00340] Aspect 35 is the method of any preceding aspect, wherein the genomic DNA is from a single cell or is a mixture from a plurality of cells.

[00341] Aspect 36 is the method of any preceding aspect, wherein the cytidine deaminase is a member of the AID subfamily, the APOB EC 1 subfamily, the APOBEC2 subfamily, the APOBEC3A subfamily, the APOBEC3B subfamily, the APOBEC3C subfamily, the APOBEC3D subfamily, the APOBEC3F subfamily, the APOBEC3G subfamily, the APOBEC3G subfamily, the APOBEC3H subfamily, or the APOBEC4 subfamily.

[00342] Aspect 37 is the method of any preceding aspect, wherein the cytidine deaminase is an altered cytidine deaminase comprising an amino acid substitution mutation at a position functionally equivalent to (Tyr/Phe)130 in a wild-type APOBEC3A protein

[00343] Aspect 38 is the method of any preceding aspect, wherein the substitution mutation at the position functionally equivalent to Tyrl30 comprises a mutation to Ala, Vai, or Trp.

[00344] Aspect 39 is the method of any preceding aspect, wherein the altered cytidine deaminase comprises amino acid substitution mutations at positions functionally equivalent to (Tyr/Phe)130 and Tyrl32 in a wild-type APOBEC3A protein.

[00345] Aspect 40 is the method of any preceding aspect, wherein the substitution mutation at the position functionally equivalent to Tyrl30 comprises a mutation to Ala, Vai, or Trp and the substitution mutation at the position functionally equivalent to Tyrl32 comprises a mutation to Arg, His, Leu, or Gin, or a combination thereof..

[00346] Aspect 41 is the method of any preceding aspect, wherein the cytidine deaminase or the altered cytidine deaminase comprises a ZDD motif H- [P/A/V]-E-X[23-28]-P-C-X[2-4]-C (SEQ ID NO: 12).

[00347] Aspect 42 is the method of any preceding aspect, wherein the cytidine deaminase or the altered cytidine deaminase is a member of the APOBEC3A subfamily and comprises a ZDD motif HXEX24SW(S/T)PCX[2-4]CX6FX8LX5R(L/I)YX[8-1 l]LX2LX[10]M (SEQ ID N0: 13). [00348] Aspect 43 is a composition comprising a ten-eleven translocase (TET); and a CpC binding protein.

[00349] Aspect 44 is the composition of any preceding aspect, wherein the composition further comprises double-stranded DNA.

[00350] Aspect 45 is the composition of any preceding aspect, wherein the CpG binding protein is an ID AX CXXC domain protein.

[00351] Aspect 46 is the composition of any preceding aspect, 45 wherein the CpG binding protein is a fusion protein.

[00352] Aspect 47 is the composition of any preceding aspect, wherein the fusion protein comprises a domain comprising the activity of binding to a single-stranded nucleic acid.

[00353] Aspect 48 is the composition of any preceding aspect, wherein the fusion protein comprises a domain comprising the activity of increasing the solubility of the CpG binding protein.

[00354] Aspect 49 is a composition comprising a means for converting 5mC to 5caC; and a means for binding to unmethylated CpG dinucleotides.

[00355] Aspect 50 is a method for identifying methylated nucleotides comprising: providing a sample of DNA suspected of comprising double-stranded DNA comprising at least one 5- methyl cytosine (5mC); contacting the double-stranded DNA with a CpG binding protein, wherein the CpG binding protein binds unmethylated CpG dinucleotides present in the double-stranded DNA; and exposing the double-stranded DNA to a ten-eleven translocase (TET), wherein the TET oxidizes 5mC to 5-carboxy cytosine (5CaC).

[00356] Aspect 51 is a method for reducing the inhibition of oxidation of 5mC by a ten-eleven translocase (TET), the method comprising: providing a sample of DNA suspected of comprising double-stranded DNA comprising at least one 5-methyl cytosine (5mC); contacting the double-stranded DNA with a CpG binding protein, wherein the CpG binding protein binds unmethylated CpG dinucleotides present in the double-stranded DNA; and exposing the DNA to a ten-eleven translocase (TET), wherein the TET oxidizes 5mC to 5- carboxy cytosine (5CaC). [00357] Aspect 52 is the method of any preceding aspect, wherein the CpG binding protein is an ID AX CXXC domain.

[00358] Aspect 53 is the composition of any preceding aspect, wherein the CpG binding protein is a fusion protein.

[00359] Aspect 54 is the composition of any preceding aspect, wherein the fusion protein comprises a domain comprising the activity of binding to a single-stranded nucleic acid.

[00360] Aspect 55 is the composition of any preceding aspect, wherein the fusion protein comprises a domain comprising the activity of increasing the solubility of the CpG binding protein.

[00361] Aspect 56 is the method of any preceding aspect, wherein the contacting and the exposing occur at the same time or the exposing occurs before the contacting.

[00362] Aspect 57 is the method of any preceding aspect, further comprising contacting the doublestranded DNA with a borane reducing agent after the exposing to result in a reduced double-stranded DNA, wherein the borane reducing agent converts 5CaC to dihydrouracil (DHU).

[00363] Aspect 58 is the method of any preceding aspect, wherein the borane reducing agent is pyridine borane or 2-picoline borane (pic- BH3).

[00364] Aspect 59 is the method of any preceding aspect, further comprising contacting the doublestranded DNA with a glucosyltransferase that glucosylates 5hmC and does not glucosylates 5mC.

[00365] Aspect 60 is the method of any preceding aspect, wherein the glucosyltransferase comprises p-glucosyltransferase.

[00366] EXAMPLES

[00367] The present disclosure is illustrated by the following examples. It is to be understood that the particular examples, materials, amounts, and procedures are to be interpreted broadly in accordance with the scope and spirit of the disclosure as set forth herein. [00368] Example 1

[00369] Experimental assay for cytidine deaminase activity

[00370] Deamination by cytidine deaminases was monitored using a gel-based assay with the restriction enzyme rral (FIG. 4).

[00371] Deamination of C to U or 5mC to T results in a mismatched DNA substrate, which can be cleaved by Swal. These products are visualized as two distinct species on a denaturing polyacrylamide gel.

[00372] Upon incubation with APOBEC3 A wildtype or its mutants with 5' fluorescein amidite (FAM) labeled oligonucleotide substrate containing C or 5mC, deamination and no deamination occurs depending on the enzyme substrate preference. Introduction of a complementary oligonucleotide thereafter may result in complete matched or mismatched base pairing. If there is a mismatch, Sirc/I will cleave the double- stranded oligonucleotide. The original or cleaved 5'-FAM labeled oligonucleotide can be visualized by the extent of migration using 15% Urea-PAGE with FAM filter. Although FAM is used here, essentially any label can be used, and either the 5' or the 3' end can be labeled.

[00373] This assay was adapted and modified from Schutsky et al., Nucleic Acid Research, 45, 7655-7665, 2017. doi: 10.1093/nar/gkx345. Modifications to Schutsky et al. included the following. Instead of performing DNA precipitation and redissolving the DNA substrate into Swal compatible buffer, 1 uL of the altered cytidine deaminase APOBEC3A(Y130A) deamination reaction mixture was aliquoted into 9 uL Swal compatible buffer for restriction enzyme digestion for our gel assay. Appropriate controls were performed to determine the Swal restriction enzyme digestion efficiency was not compromised by the APOBEC reaction buffer. Instead of introducing 1.5-fold excess complementary strand prior to overnight Nvi zzI restriction enzyme digestion, 3 -fold excess complementary strand was introduced. Instead of running the pre-heated 20% acrylamide/Tris-Borate-EDTA (TBE)/urea gel reported by Schutsky et al., the gel run was performed at room temperature with a 15% acrylamide/Tris-Borate-EDTA (TBE)/urea gel and observed good resolution between cut (deaminated) and uncut (unreacted) oligo substrates. [00374] The Swal assay was first validated using FAM-labeled DNA oligonucleotides that contain a C, 5mC, U, or T residue (FIG. 5A). Oligonucleotides purchased from Integrated DNA Technologies (IDT) and visualized by the 15% Urea-PAGE with FAM filter. Synthesized oligonucleotides. oLB 1609 contains substrate C, and 0LBI6IO contains its corresponding deaminated product U. 0LBI611 contains substrate 5mC, and oLB1612 contains its corresponding deaminated product T. oLB 1679 is the complementary oligonucleotide for OLB1609, 0LBI6IO, 0LBI611, and OLB1612 (FIG. 5A). As shown in FIG. 5B, annealing of oLB1609/oLB1679 and oLB 1611/oLB 1679 resulted in complete base pairing of the oligonucleotide. Addition of Anv did not result in cleavage of this substrate. However, a single base mismatched substrate formed by annealing of oLB1610/oLB1679 (U/G mismatched) or oLB1612/oLB1679 (T/G mismatched) resulted in a cleavage product upon addition of Swal. As a specific cleavage product was observed in the presence of U or T, but not C or 5mC, the Swal assay could serve as a readout of C to U and 5mC to T deamination by APOBEC3A.

[00375] As an additional control, C and 5mC-containing DNA oligonucleotides oLB1609 (C oligo) and oLB1612 (5mC oligo) were incubated with APOBEC3A enzyme purchased from New England Biolabs (NEBNext® Enzymatic Methyl-seq Kit (Catalog # E7120)), which was reported to deaminate both C and 5mC efficiently (FIG. 6). The deamination reaction mixture (5 pL, 2 pL and 1 pL) were then directly added to Ntv/l assay buffer containing SM’M, resulting in a total volume of 10 pL. 1 pL of deamination reaction was sufficient to observe the cut band, allowing quantification of the extent of deamination. Treatment with NEB APOBEC3A and subsequent digestion by Swal resulted in formation of a cleavage product with similar mobility to both U and T-containing substrates in FIG. 5A. These data further support the notion that Swal cleavage can be used to monitor the deaminase activity of APOBEC3A and other cytidine deaminases.

[00376] Example 2

[00377] Purification of APOBEC3A(Y130X) mutant proteins

[00378] The impact of all possible amino acid substitutions at position 130 of APOBEC3A on the deaminase activity of this enzyme was systematically assessed. To this end, 19 different His-tagged APOBEC3 A constructs were cloned, each encoding a different amino acid at position 130 relative to the wild type tyrosine. The corresponding proteins were expressed in BL21(DE3) cells, purified using Ni-NTA agarose beads, and desal ted/concentrated using spin columns to storage buffer (50mM Tris pH 7.5, 200mM NaCl, 5%(v/v) glycerol, 0.01% (v/v) Tween-20, 0.5mM DTT). This yielded APOBEC3A(Y130X) mutant protein preparations with 80-85% purity, as judged by SDS-PAGE analysis (FIG. 7A).

[00379] Example 3

[00380] DNA deaminase activity of APOBEC3A(Y130X) mutant proteins

[00381] The deaminase activity of all purified APOBEC3A(Y130X) proteins was then analyzed using the AMY/I assay, with a 37°C/2 hour reaction time and NEB APOBEC3A as positive control (FIG. 8). 10-20 pM final concentration of Y130X recombinant enzymes were incubated with oLB1609 (C oligo, top panel) and oLB1612 (5mC oligo, bottom panel) at 37°C for 2 hours. NEB APOBEC3A enzyme was purchased from NEBNext® Enzymatic Methyl-seq Kit (Catalog # E7120). Wild type APOBEC3A deaminated 5mC and C substrates to completion, consistent with previous literature. Different mutants exhibited a wide range of reactivities towards 5mC and C substrates, with some showing preference towards either substrate. Remarkably, APOBEC3A(Y130A) (first box) deaminated 5mC substrates almost completely (94.2%), while it deaminated the corresponding C substrate to a minor extent (29.4%). Other mutants, such as APOBEC3A(Y130P) and APOBEC3A(Y130T), also exhibited more complete deamination of the 5mC than C substrate, albeit to a lesser extent than APOBEC3A(Y130A). In contrast, APOBEC3A(Y130L) (second box) deaminated approximately half of the C substrate (56%), but almost none of the 5mC substrate (6.8%). The deaminase activity of all APOBEC3A(Y130X) mutants is quantified and summarized in FIG. 8, FIG. 9, and FIG.

10

[00382] Because these AIIY/I assays were performed as a single endpoint measurement (2 hour), it could be possible that the respective deamination reactions had already saturated. A time course analysis of APOBEC3A(Y130A) deaminase activity was therefore performed. The extent of C and 5mC deamination was monitored at 0, 5, 10, 30, 60 and 120 minutes by incubation of -10-20 pM of APOBEC(Y130A) with 500nM C and 5mC oligonucleotide substrate (FIG. 11 A, 11B, 11C, 11D, HE, HF, ). A greater difference in the extent of 5mC versus C deamination was observed at t < 30 min.

[00383] The kinetics of deamination by wild type APOBEC3A and mutant APOBEC3A(Y130A) were quantitively compared. The initial deamination reaction velocity was measured at a range of DNA substrate concentrations and used to construct Michaelis-Menten curves for 5mC and C substrates, respectively. The resulting Km and Kcat values were then derived from these data. The catalytic efficiency of APOBEC3A(Y130A) was -100-fold higher on 5mC than C substrates (FIG. 12), corroborating the endpoint YI /I assays shown in FIG. 8, FIG. 9, and FIG. 10

[00384] Example 4

[00385] Purification of APOBEC3A(Y130A-Y132H) double mutant protein

[00386] APOBEC3A(Y130A-Y132H) protein was expressed in BL21(DE3) cells, purified using Ni-NTA agarose beads, and desal ted/ concentrated using spin columns to storage buffer (50mM Tris pH 7.5, 200mM NaCl, 5%(v/v) glycerol, 0.01% (v/v) Tween-20, 0.5mM DTT). This yielded APOBEC3A(Y130A-Y132H) mutant protein preparations with 90- 95% purity, as judged by SDS-PAGE analysis (FIG. 7B).

[00387] Example 5

[00388] DNA deaminase activity of APOBEC3A(Y130A-Y132H) double mutant protein

The deaminase activity of purified APOBEC3A(Y130A-Y132H) double mutant protein was then analyzed using the ,S'n aI assay, with a 37°CZ 2 hour reaction time and NEB APOBEC3A as positive control. The conditions used were the same as described in Example 3 with the exception that the S I assay used reaction conditions of 40 mM sodium acetate pH 5.2, 37°C for 1 hour to 16 hours. The DNA substrates are shown in FIG. 13. After the deaminase reaction the deaminated oligo substrates were PCR- amplified, sequenced, and the number of C and 5mC deamination events per read were counted The DNA oligonucleotide substrates used for experiments with

I l l APOBEC3A(Y130A-Y132H) are shown in FIG. 13. APOBEC3A(Y130A-Y132H) exhibited higher levels of deamination at all methylated sites compared to unmethylated sites. This was consistent across both CpG and non-CpG contexts, and was robust to variation in reaction time (FIG. 14, 15). The difference in deamination level between methylated and unmethylated sites was markedly higher for APOBEC3A(Y130A-Y132H) than APOBEC3A(Y130A), indicating that APOBEC3A(Y130A-Y132H) achieves better discrimination of methylated sites than APOBEC3A(Y130A). In addition, APOBEC3A(Y130A-Y132H) deaminated methylated sites more efficiently than unmethylated sites across all xCpGx motifs (FIG. 16).

[00389] Example 6

[00390] Purification of APOBEC3A(Y130W) mutant protein

[00391] Recombinant human His-tagged APOBEC3A(Y130W) protein was expressed in Escherichia coli BL21(DE3) cells, purified using Ni-NTA affinity chromatography, and desal ted/ concentrated using spin columns to storage buffer (50mM Tris pH 7.5, 200mM NaCl, 5%(v/v) glycerol, 0.01% (v/v) Tween-20, 0.5mM DTT).

[00392] Example 7

[00393] DNA deaminase activity of APOBEC3A(Y130W) mutant protein

[00394] Using the Yit/I assay, we measured the activity of APOBEC3A(Y130W) on C, 5mC, and 5hmC substrates across a 90 min time interval. While the majority of C and 5mC was deaminated, no detectable activity was observed for 5hmC. (FIG. 17A and 17B).

[00395] To better understand the substrate preferences of APOBEC3A(Y130W), we assayed deaminase activity on C, 5mC, and all oxidized derivatives of 5mC (FIG. 18A and 18B). The Swal restriction enzyme assay was performed with a 90 min reaction time to measure deaminase activities of APOBEC3A wild type and Y130W mutant enzymes. A no-protein control was included to account for potential degradation of oligonucleotide substrates and account for non-specific activity of Swal during the course of the assay. APOBECY130A(Y130W) showed no detectable deamination of 5hmC, 5fC, and 5caC, while wild type APOBEC retained significant deaminase activity towards 5hmC and showed residual activity on 5fC and 5caC. Both enzymes deaminated C and 5mC efficiently.

[00396] Example 8

[00397] Detection of 5mC

[00398] The following example generally follows the methods described in Schutsky et al. Nature biotechnology, 10.1038/nbt.4204. 8 Oct. 2018, doi: 10.1038/nbt.4204. Genomic DNA (gDNA) from an organism suspected of having 5hmC is provided. 20 ng of the gDNA mixture is then treated as follows.

[00399] In order to provide single-stranded DNA and facilitate deamination activity, 1 pL of DMSO is added and the sample is denatured at 95 °C for 5 min and snap cooled by transfer to a PCR tube rack pre-incubated at -80 °C. Before thawing, reaction buffer is overlaid to a final concentration of 20 mM MES pH 6.0 + 0.1 % Tween, and an altered cytidine deaminase described herein is added to the sample to a final concentration of 5 pM in a total volume of 10 pL. The deamination reaction is incubated under linear ramping temperature conditions from 4-50 °C over 2 hrs.

[00400] After deamination, the sample is prepared for Illumina sequencing using the Accel Methyl- NGS kit (Swift Biosciences). Specifically, the reactions are purified using Zymo Oligo Clean and Concentrator Kit and eluted in 15 pL elution buffer (10 mM Tris pH 8.0). The Accel-NGS Methyl-Seq kit (Swift Biosciences) is then used for library preparation of single-stranded DNA according to manufacturer’s instructions. After purification of the library, 1-5 ng amplified DNA is run on an Agilent Bioanalyzer High Sensitivity DNA Chip to confirm proper library fragment sizes. The resulting ACE-Seq library is sequenced at 1.9 pM with single-end mode on a NextSeq 500 sequencer (Illumina) using the NextSeq 500/500 High Output kit v2 (150 cycles). [00401] Sequencing data from the sample is analyzed to detect 5mC at 500 CpG alleles suspected of 5mC modification. Those CpG sites that sequence as thymidine are characterized as likely containing 5mC modifications in the original sample.

[00402] Example 9

[00403] Detection of 5hmC

[00404] The following example generally follows the methods described in Schutsky et al., Nature biotechnology, 10.1038/nbt.4204. 8 Oct. 2018, doi:10.1038/nbt.4204. Genomic DNA (gDNA) from an organism suspected of having 5hmC is provided. Two aliquots of 20 ng each are provided. One of the aliquiots of 20 ng of the gDNA mixture is glucosylated using UDP-glucose and T4 P-glucosyltransferase (PGT, NEB). Specifically, 0.5 pL of 10X Cutsmart Buffer (NEB), 0.1 pL of 50X UDP-Glucose, 0.5 pL of T4 pGT (NEB) are combined with 20ng of gDNA and water is added for a total volume of 5pL. A control reaction is assembled using the other 20ng aliquot of gDNA, and the reaction components are mixed as described above, but with 0.5 pL water added in place of T4 PGT. The reactions are incubated at 37 °C for 1 hr.

[00405] In order to provide single-stranded DNA and facilitate deamination activity, 1 pL of DMSO is added and the samples are denatured at 95 °C for 5 min and snap cooled by transfer to a PCR tube rack pre-incubated at -80 °C. Before thawing, reaction buffer is overlaid to a final concentration of 20 mM MES pH 6.0 + 0.1 % Tween, and an altered cytidine deaminase is added to each sample to a final concentration of 5 pM in a total volume of 10 pL. The deamination reactions are incubated under linear ramping temperature conditions from 4-50 °C over 2 hrs.

[00406] After deamination, the samples are prepared for Illumina sequencing using the Accel Methyl-NGS kit (Swift Biosciences). Specifically, the reactions are purified using Zymo Oligo Clean and Concentrator Kit and eluted in 15 pL elution buffer (10 mM Tris pH 8.0). The Accel-NGS Methyl-Seq kit (Swift Biosciences) is used for library preparation according to manufacturer’s instructions. After purification of the libraries, 1-5 ng amplified DNA is run on an Agilent Bioanalyzer High Sensitivity DNA Chip to confirm proper library fragment sizes. The resulting ACE-Seq libraries are sequenced at 1 .9 pM with single-end mode on a NextSeq 500 sequencer (Illumina) using the NextSeq 500/500 High Output kit v2 (150 cycles).

[00407] Sequencing data from the sample and control are analyzed to detect 5hmC and 5mC at 500 CpG alleles suspected of 5hmC modification. Specifically, any CpG sites that are sequenced as cytosine in the control sample are compared to the same CpG sites in the sample treated with [3GT. Those sites that sequence as thymidine are characterized as likely containing 5hmC modifications in the original sample.

[00408] Example 10

[00409] Detection of 5mC and 5hmC

[00410] The following example generally follows the methods described in Schutsky et al. Nature biotechnology, 10.1038/nbt.4204. 8 Oct. 2018, doi: 10.1038/nbt.4204. Genomic DNA (gDNA) from an organism suspected of having 5mC and 5hmC is provided. Two aliquots of 20 ng each are provided. Each 20 ng aliquot of the gDNA mixture is then treated as follows.

[00411] In order to provide single-stranded DNA and facilitate deamination activity, 1 pL of DMSO is added and the sample is denatured at 95 °C for 5 min and snap cooled by transfer to a PCR tube rack pre-incubated at -80 °C. Before thawing, reaction buffer is overlaid to a final concentration of 20 mM MES pH 6.0 + 0.1 % Tween. To one of the aliquots, and the altered cytidine deaminase (Y130W) described herein is added to the sample to a final concentration of 5 pM in a total volume of 10 pL. To the other aliquot, wild type cytidine deaminase is added to the sample to a final concentration of 5 pM in a total volume of 10 pL. The deamination reactions are incubated under linear ramping temperature conditions from 4-50 °C over 2 hrs.

[00412] After deamination, the samples are prepared for Illumina sequencing using the Accel Methyl-NGS kit (Swift Biosciences). Specifically, the reactions are purified using Zymo Oligo Clean and Concentrator Kit and eluted in 15 pL elution buffer (10 mM Tris pH 8.0). The Accel-NGS Methyl-Seq kit (Swift Biosciences) is then used for library preparation of single-stranded DNA according to manufacturer’s instructions. After purification of the library, 1-5 ng amplified DNA is run on an Agilent Bioanalyzer High Sensitivity DNA Chip to confirm proper library fragment sizes. The resulting APOBEC-coupled epigenetic sequencing (ACE-seq) library is sequenced at 1.9 pM with single-end mode on a NextSeq 500 sequencer (Illumina) using the NextSeq 500/500 High Output kit v2 (150 cycles).

[00413] Sequencing data from the sample is analyzed to detect 5mC at 500 CpG alleles suspected of 5mC and/or 5hmC modification. Those CpG sites that sequence predominantly as thymidine in the wild type sample but that sequence predominantly as cytosine in the Y130W sample are characterized as likely containing 5hmC modifications in the original sample.

[00414] Example 11

[00415] Detection of 5mC and 5hmC

[00416] Human genomic DNA is combined with fully unmethylated lambda control DNA (New England Biolabs) and enzymatically CpG methylated pUC19 control DNA and mechanically sheared to give fragments of approximately ~300bp. This sheared DNA (10- lOOng) is then subjected to end-repair, A-tailing and adapter ligation according to standard Illumina library preparation procedures. The sample is then split into 2 aliquots to enable differential treatment.

[00417] One of the aliquots of the adapter-ligated DNA is glucosylated using UDP-glucose and T4 P-glucosyltransferase (PGT, NEB). Specifically, the sample is treated with 0.5 U/uL of T4 PGT (NEB) in IX CutSmart Buffer (NEB) with 40 pM UDP-Glucose at 37 °C for 3 hr. A control reaction is assembled using the other aliquot of the DNA sample, omitting the T4 PGT enzyme. Subsequently, the samples are optionally SPRI purified and then denatured via incubation in 0.02 N sodium hydroxide at 50°C for 10 minutes. Subsequently, ssDNA samples are enzymatically deaminated in 50 mM Bis-Tris (pH 6.5), 10 pg/mL RNAse A with the cytidine deaminase (200nM) for 25 minutes at 37C. The libraries are then PCR amplified using unique-dual indexing primers and Q5U (New England Biolabs) using 9 cycles of PCR. Samples are sequenced on a NovaSeq6000 and analysis is performed with DRAGEN Methylation Pipeline. The methylation calls between the two treatment conditions (+/- T4 PGT) are compared in order to determine which sites were hydroxymethylated. Bases that are detected as methylated in both sample treatments are assigned mC, whereas bases that are detected as methylated in the -PGT condition and unmethylated in the + PGT condition are assigned hmC (FIG. 19).

[00418] Example 12

[00419] Generalized Method for Creation of Methylation Sequencing libraries

[00420] Sheared DNA or cfDNA (10-200ng) was first subjected to end-repair, A-tailing and adapter ligation according to standard library preparation procedures. Suitable adapters for this method can have unmodified C’s or pyrrolo-C modifications to disfavor deamination of the adapter sequence. The adapter ligated DNA was then denatured to ssDNA. Subsequently, this ssDNA sample was enzymatically deaminated in a buffered solution with an engineered deaminase (APOBEC3A-Y130A-Y132H, 50 nM-1000nM) for 5 minutes to 3 hours at incubation temperatures ranging between 20°C to 55°C. Deaminated libraries were optionally SPRI purified before PCR amplification using unique-dual indexing primers, using either a uracil tolerant polymerase (for instance, Q5U®, KapaU™) or a uracil intolerant polymerase (for instance, Q5 HiFi®, KAPA HiFi™) using 9 to 12 cycles of PCR. Libraries were then sequenced on a NextSeq550 or NovaSeq 6000 and analysis was performed with DRAGEN Methylation Pipeline.

[00421] Methods for denaturation

[00422] A variety of methods for denaturation known to those skilled in the art. These include, but are not limited to: (i) heating in the presence of NaOH or high pH buffer at a moderate temperature (e.g., 0.02 N sodium hydroxide at 50°C for 10 minutes); (ii) heating to high temperatures (e.g., 95°C for 10 minutes); (iii) heating in the presence of DMF (e.g., 50% DMF at 95°C for 10 minutes); (iv) heating in the presence of formamide (e.g., 50% formamide at 95°C for 10 minutes); and (v) heating in the presence of DMSO (e.g., 50% DMSO at 95°C for 10 minutes). [00423] Suitable buffers for deamination

[00424] Typically, the deamination buffer is added to the denatured DNA sample, and thus any additives to promote denaturation will be present at a lower (diluted) concentration in the deamination reaction.

[00425] Several buffer systems have been identified in which deamination can occur. The buffer 50 mM Bis Tris, pH 6.5 can be used, and other pH levels between 5-7.5 may be feasible. Additionally, other buffer strengths (concentrations) may be feasible. The skilled person will recognize that if NaOH is used for denaturation, adequate buffer strength/pH can be used to ensure that the final pH is within the ideal range for the deaminase. Alternative buffer systems (e g., tcichemicals.com/US/en/c/10367), including MES, may also provide good results.

[00426] In some embodiments, other buffer systems have been shown to have lower performance when used for the deamination of a genomic DNA (gDNA) sample. Examples include citrate buffer and sodium acetate buffer (FIG. 20). These buffer systems may have lower performance due to sodium content, which is believed to be detrimental to activity of the altered deaminase. APOBEC3A is a zinc-dependent enzyme (Marx et al., Scientific Reports 5, no. December (2015): 1-9. https://doi.org/10.1038/srepl 8191), and thus the presence of metal cations (e.g., sodium, magnesium) may generally be detrimental. Another explanation is that metal cations can impact the formation of secondary structure in nucleic acids (Einert et al., Biophysical Journal 100, no. 11 (2011): 2745-53, doi.org/10.1016/j.bpj.2011.04.038), which may also impact activity of the deaminase on the substrate, since APOBEC3A cannot deaminate Cs in double stranded DNA such as hairpins. Therefore, it would be expected that any buffer system can be used that would maintain pH in the desired range, and in some embodiments buffers with minimal sodium can be used.

[00427] Selectivity of altered cytidine deaminases on genomic DNA

[00428] The altered cytidine deaminase maintains some activity for unmethylated cytidine deamination, thus it can be useful to tune the activity to maintain the desired selectivity for methylated cytosines. For example, high concentrations of the enzyme, high incubation times, or buffers that promote high activity of the enzyme may lead to undesired levels of unmethylated cytidine deamination. Typically, enzyme concentration may be from 50 nM to lOOOnM), the reaction can occur from 5 minutes to 3 hours at incubation temperatures ranging from 20°C to 55°C. In particular, it may be helpful to tune enzyme concentration based on the purification method of the deaminase and how much activity a given enzyme preparation contains.

[00429] As an example, libraries were prepared, deaminated according to the general methods described above, and sequenced. Variation of the concentration of APOBEC -Y130A- Y132H led to observable differences in the level of activity and level of off -target cytidine deamination (Fig. 20A). Furthermore, although one may expect a commercial buffer for APOBEC activity to work well for this application, testing of the APOBEC buffer included in the EM-seq™ kit (New England Biolabs) led to an undesirable high level of activity on the cytidine substrate, resulting in elevated conversion of unmethylated C nucleobases in lambda DNA (FIG. 20B).

[00430] Methods for determining suitable conditions for deamination

[00431] To determine whether conditions were suitable for selective deamination on genomic DNA, a DNA mixture containing NA12878 (human) DNA, fully CpG methylated pUC19, and fully unmethylated lambda was prepared. After adapter ligation, the samples were subjected to various deamination conditions and prepared as a sequencing library, using the general methods described above. The methylation level of pUC19 and lambda was used to select conditions with desired 5mC activity and selectivity.

[00432] Example 13

[00433] Specific Examples of Conditions used for Preparation and Deamination of Libraries by Altered Cytidine Deaminase

[00434] In all methods described below in Examples 13-20, the altered cytidine deaminase refers specifically to APOBEC3A-Y130A-Y132H. A variety of human genomic DNA samples were tested for different application assessments. [00435] Method A: Human genomic DNA was combined with fully unmethylated lambda control DNA and enzymatically CpG methylated pUC19 control DNA and mechanically sheared to give fragments of approximately ~300bp. This sheared DNA (10-100ng) was then subjected to end-repair, A-tailing, and adapter ligation according to standard Illumina library preparation procedures. The adapter ligated DNA was denatured via incubation in 0.02 N sodium hydroxide at 50°C for 10 minutes. Subsequently, ssDNA samples were enzymatically deaminated in 50 mM Bis-Tris (pH 6.5), 10 pg/mL RNAse A with the cytidine deaminase (200nM) for 25 minutes at 37°C. The libraries were then PCR amplified using unique-dual indexing primers and Q5U (New England Biolabs) using 9 cycles of PCR. Samples were sequenced on aNovaSeq6000 and analysis was performed with DRAGEN Methylation Pipeline.

[00436] Method B: Human genomic DNA was combined with fully unmethylated lambda control DNA and enzymatically CpG methylated pUC19 control DNA and mechanically sheared to give fragments of approximately ~300bp. This sheared DNA (10-100ng) was then subjected to end-repair, A-tailing, and adapter ligation according to standard Illumina library preparation procedures. The adapter ligated DNA was denatured via incubation in 0.02 N sodium hydroxide at 50°C for 10 minutes. Subsequently, ssDNA samples were enzymatically deaminated in 50 mM Bis-Tris (pH 6.5), 10 pg/mL RNAse A with the altered cytidine deaminase (200nM) for 15 minutes at 37°C. The libraries were then PCR amplified using unique-dual indexing primers and Q5U (New England Biolabs) using 9 cycles of PCR. Samples were sequenced on a NovaSeq6000 and analysis was performed with DRAGEN Methylation Pipeline.

[00437] Method C: Human genomic DNA was combined with fully unmethylated lambda control DNA and enzymatically CpG methylated pUC19 control DNA and mechanically sheared to give fragments of approximately ~300bp. This sheared DNA (10-100ng) was then subjected to end-repair, A-tailing, and adapter ligation according to standard Illumina library preparation procedures. The adapter ligated DNA was denatured via incubation in 0.02 N sodium hydroxide at 50°C for 10 minutes. Subsequently, ssDNA samples were enzymatically deaminated in 50 mM Bis-Tris (pH 6.5), 10 pg/mL RNAse A with the altered cytidine deaminase (200nM) for 15 minutes at 37°C. The libraries were then PCR amplified using unique-dual indexing primers and Q5 HiFi (New England Biolabs) using 9 cycles of PCR. Samples were sequenced on aNovaSeq6000 and analysis was performed with DRAGEN Methylation Pipeline.

[00438] Method D: Human genomic DNA was combined with fully unmethylated lambda control DNA and enzymatically CpG methylated pUC19 control DNA and mechanically sheared to give fragments of approximately ~300bp. This sheared DNA (10-100ng) was then subjected to end-repair, A-tailing and adapter ligation according to standard Illumina library preparation procedures. The adapter ligated DNA was denatured via incubation in 0.02 N sodium hydroxide at 50°C for 10 minutes. Subsequently, ssDNA samples were enzymatically deaminated in 50 mM Bis-Tris (pH 6.5) with the altered cytidine deaminase (200nM) for 25 minutes at 37°C. The libraries were then PCR amplified using unique-dual indexing primers and Q5U (New England Biolabs) using 9 cycles of PCR. Samples were sequenced on a NovaSeq6000 and analysis was performed with DRAGEN Methylation Pipeline.

[00439J Use ofRNAse A

[00440] Literature suggests that RNAses may increase the activity of cytidine deaminases by removing contaminating RNA (Bransteitter et al ., Proceedings of the National Academy of Sciences of the United States of America 100, no. 7 (2003): 4102-7. https://doi.org/10.1073/pnas.0730835100). However, testing ofRNAse A showed, contrary to this hypothesis, that RNAse A reduced activity of the altered cytidine deaminase, with a more pronounced impact to the reduction of off-target cytidine deamination, thus leading to the result of more 5mC selectivity (FIG. 21).

[00441] Example 14

[00442] Deamination of 5hmC

[00443] To assess 5hmC deamination, oligos were designed with C, 5mC, or 5hmC modifications in defined contexts. These oligos were assembled together using ligation, and oligos were built such that the resulting ligated fragment would contain handles needed for subsequent amplification and sequencing (FIG. 22A). The assembled control oligo was spiked into an adapter-ligated DNA library, and treated according to Method A (Example 12). Analysis of the reported methylation from this control oligo showed that 5mC was the most preferred substrate for APOBEC Y130A-Y132H, however there was still significant activity on 5hmC (FIG. 22B).

[00444] Example 15

[00445] Determination of Methylation on CpG islands

[00446] Using NA12878 gDNA, libraries were prepared and deaminated according to Method A (Example 13). A comparative dataset was also generated using EM-Seq™ conversion (New England Biolabs). To assess the methylation performance on the human genome, regional methylation values for CpG islands across the human genome were calculated using methylpy. The per-region methylation values were plotted against the per-region methylation values derived from the EM-seq™ dataset (FIG. 23). The regional analysis demonstrates high correlation between the use of the Y130A-Y132H deaminase and EM- Seq™, a commercial methylation detection method. This indicates that the 5mC-selective deamination assay described herein can effectively detect and report methylation in CpG islands.

[00447] Example 16

[00448] Variant calling on methylated libraries

[00449] Libraries from human genome samples NA12878, NA24385, and NA24631 were prepared and deaminated according to Method A (Example 13). A comparative dataset was also generated without conversion, by proceeding directly to PCR after adapter ligation. After running variant calling analysis, comparison to truth sets for each genome using hap.py showed that the libraries subjected to the methylation conversion (the altered cytidine deaminase) provided good SNV/indel calling performance, approaching the performance of the no conversion controls (FIG. 24). While SNV/indel calling is discussed here, other types of variant calling, including copy number variation (CNV) and short tandem repeats (STR), and structural variants (SV) are also feasible. [00450] Example 17

[00451] Differentially methylated region (DMR) Calling

[00452] Differential methylation analysis is commonly employed for comparison of methylomes across different diseases, tissues, and cell types (Chen et al., Briefings in Functional Genomics 15, no. 6 (2016): 485-90. https://doi.org/10.1093/bfgp/elw018). Libraries were prepared and deaminated according to Method A (Example 13) using genomic DNA isolated from HCC2218 tumor and normal cell types (CRL-2363D, CRL-2343D, ATCC) and HCC1187 tumor and normal cell types (CRL-2323D, CRL-2322D, ATCC). Comparative dataset was also generated using EM-Seq™ conversion (New England Biolabs) and Bisulfite conversion (EZ DNA Methylation-Gold Kit - Zymo Research). For differential methylation analysis, HOME, a program for identifying DMRs, was used (Srivastava et al., BMC Bioinformatics 20, no. 1 (2019): 1-15, doi.org/10.1186/sl2859- 019-2845-y). Comparisons of the methylation between tumor and normal samples were performed for each methylation assay. As shown in FIG. 25A, 25B, 25C, 25D, 25E, and 25F, the assay was able to detect a relevant differentially methylated region in the ZNF-154 gene (Almeida et al., BMC Cancer 19, no. 1 (2019): 1-12. https://doi.org/10.1186/sl2885- 019-5403-0). Furthermore, quantitative assessment of DMR calling performance using Methods A, B, and C (Example 13) showed that the altered cytidine deaminase-based methods detected expected DMRs with high precision and recall (FIG. 26).

[00453] Example 18

[00454] Methylation at promoters

[00455] The methylation status of promoter regions can be correlated to both histone status as well as gene expression activity. In order to assess whether the altered cytidine deaminase assay could be used to probe the methylation status of promoters, libraries from human genome samples NA12878 were prepared and deaminated according to Method A (Example 13). Comparative libraries were also generated with EM-seq™ conversion. The methylation level at known histone marker sites, including H3K36me3 and H3K27ac were quantified from the resulting data. H3K36me3 sites are expected to be inactive promoters with high histone and DNA methylation, while H3K27ac sites are expected to be active promoters with high histone acetylation and low DNA methylation. The altered cytidine deaminase assay was able to report these methylation trends, as expected (FIG. 27A and 27B).

[00456] Example 19

[00457] Detecting tumor sample in a background of normal sample

[00458] Libraries from human genome samples with either HCC2218 normal DNA or a 10% spikein of HCC2218 tumor DNA into a background of HCC2218 normal DNA were made with appropriate adapters. Assessment of methylation level for the 10% spike-in samples were made in comparison to established methylation methods on CpGs within regions also targeted by the PanSeer cancer panel (Chen, et al. Nature Communications 11, no. 1 (2020): 1-10, https://doi.org/10.1038/s41467-020-17316-z). Libraries subjected to the altered cytidine deaminase reflected methylation changes between the HCC2218 normal DNA and a 10% spike-in of HCC2218 tumor DNA into HCC2218 normal DNA. Methods A, B, and C (Example 13) produced higher tumor signal in the 10% tumor spike-in sample compared to a sample with no tumor DNA spiked in (FIG. 28). This reflects the ability of the altered cytidine deaminase to detect low levels of tumor DNA, which could be applied to cell free DNA (cfDNA) samples for early cancer detection (Jamshidi et al., Cancer Cell 40, no. 12 (2022): 1537-1549.el2. doi.org/10.1016/j.ccell.2022.10.022.) or minimal residual disease (MRD) testing (Jin et al., Proceedings of the National Academy of Sciences of the United States of America 118, no. 5 (2021): 1-8, doi.org/10.1073/pnas.2017421118.).

[00459] Example 20

[00460] Enrichment

[00461] Most on-market methylation assays convert C>T, leading to extensive C>T conversion and low complexity sequences. This makes enrichment very challenging, because probes must be designed for all possible methylation states/strands (www.twistbioscience.com/sites/default/files/resources/2020- 02/AppNote_Methylation- singles.pdf). This makes hybrid-enrichment sequencing of methyl -converted libraries expensive and more difficult to optimize. An alternative is amplicon-based targeted sequencing; however, advantages of hybridization-based targeted sequencing include the ability to easily sequence both the un-enriched (whole genome sequence (WGS)) and enriched sample and target sequencing on multiple contiguous regions of interest (Singh et al., Diagnostics (Basel, Switzerland) 12, no. 7 (June 24, 2022): 1539. doi.org/10.3390/diagnosticsl2071539.). Because of the 5mOT conversion enabled by the altered cytidine deaminase, enrichment of methyl-converted libraries can be accomplished using less complex panels while still maintaining high quality methylation information (FIG. 29). Standard probe designs prepared for non-converted libraries may be used, as only a small percentage of cytosines are methylated and are expected to be converted. Furthermore, hybridization probes tolerate mismatches (Paskey et al., BMC Genomics, 20(1). doi.org/10.1186/S12864-019-5543-2), and thus some methylation conversion can be tolerated by standard panels. In the case of highly CpG dense, methylated regions of the genome, it may be beneficial to design extra probes that account for mismatches to optimize performance.

[00462] In order to assess enrichment performance, the Illumina Methyl Capture EPIC panel, a probe panel typically employed prior to bisulfite-based conversion, was utilized. This panel is designed to target unconverted DNA due to the challenges with enrichment of converted DNA, as discussed above. However, the Illumina Methyl Capture EPIC kit requires high inputs (500ng) due to enrichment upstream of the conversion step. In contrast, in the example below, samples were converted, amplified, and then enriched, enabling lower inputs (FIG. 30A-C).

[00463] Libraries containing NA 12878 genomic DNA were prepared according to Method A (Example 13). As a control, libraries were also prepared that were not deaminated (no conversion control). Following conversion and amplification, libraries were sequenced to generate whole genome methylation control data. Samples were also enriched according to the protocols set forth in the RNA Prep with Enrichment Reference Guide (support, illumina. com/ content/ dam/illumina- support/documents/documentation/chemistry_documentation/illu mina_prep/RNA/illumina- rna-prep-reference-guide-1000000124435-03.pdf) using the Illumina Methyl Capture EPIC panel. Read enrichment was assessed using Picard CollectHsMetrics, showing that the deaminated libraries achieved similar read enrichment compared to non-converted libraries (FIG. 31A). CpG island methylation was assessed using the methods described in Example 15. Comparison of sequencing data from unenriched libraries to the enriched libraries showed that the methylation levels reported with and without enrichment showed good agreement, suggesting enrichment can be used to target and report methylation data (FIG. 31B)

[00464] Example 21

[00465] Altered cytidine deaminase-mediated detection of 5mC loci via qPCR

[00466] 5mC to T conversion can be read out by PCR using mismatch-sensitive primers. Upon 5mC to T conversion, these primers can anneal without mismatches to the substrate DNA. An unmethylated DNA substrate is not converted by an altered cytidine deaminase; accordingly, these substrates exhibit mismatches within primer annealing sequences, resulting in poor PCR amplification. Measuring amplification with respect to standard curves using qPCR would yield a quantitative readout of percent methylation at a target locus.

[00467] A 78nt ssDNA oligo containing 5mC and C was incubated with APOBEC3A(Y130A/Y132H) at 37°C for 15min, allowing 5mC deamination to occur (FIG. 32). Following this, the deamination reaction was heat-inactivated at 95°C for 5 min, and a small aliquot was directly added to a qPCR master mix containing PCR buffer, dNTPs, primers, and Q5 hot start DNA polymerase. The resulting mixture was cycled in a standard qPCR instrument. (Note: that Q5 hot start DNA polymerase was used for the qPCR assay as it is selective against uracil-containing templates, thus providing an additional layer of discrimination against spurious C to U conversion by APOBEC3A(Y130A/Y132H).) The Cq value obtained with APOBEC3A(Y130A/Y132H) and a slightly-less selective variant (Y130A) was significantly lower when compared to wild type APOBEC3A [which is non-selective for 5mC], Importantly, when a catalytically inactive variant, APOBEC3A(Y130A/Y132H_E72A), was used for deamination, Cq values match non-enzyme controls (BSA) confirming that the observed qPCR detection of 5mC is dependent on the evolved APOBEC variant.

[00468] Because the qPCR primers are selective for 5mC-converted DNA, significant qPCR amplification can be detected even though only a small fraction of substrate is deaminated. A 15 min reaction time was sufficient to observe a significant decrease in Cq value for APOBEC3A(Y130A/Y132H) (FIG. 33). Based on this initial success, we expect that alternative iterations of this detection methodology are compatible with any number of DNA amplification-based diagnostic approaches including: (1) LAMP, (2) Recombinase- polymerase amplifications, (3) other isothermal amplification methodologies.

[00469] Example 22

[00470] Detection of 5mC loci using an altered cytidine deaminase and CRISPR-Casl2

[00471] CRISPR-Cas systems have emerged as promising tools for the rapid and specific quantification of nucleic acid sequences. Techniques such as DETECTR (Chen et al., Science. 2018; 360: 436-439) and CDetection (Teng et al., Genome Biol. 2019; 20: 132) use Casl2-family enzymes to detect single-nucleotide polymorphisms in analyte DNA with high sensitivity. Casl2-family enzymes bind and cleave DNA substrates, dictated by complementarity to a guide RNA that is complexed with Casl2. Mismatches between the guide RNA and target DNA inhibit this process. Importantly, Cast 2 orthologs exhibit ‘collateral cleavage’ activity: target DNA binding results in the activation of a highly processive, non-specific nuclease activity, leading to the cleavage of ssDNA in trans. This collateral cleavage activity can be visualized using a separate ssDNA reporter substrate containing a fluorophore and quencher. RNA-guided engagement of target DNA by Cast 2 results in trans-cleavage of the ssDNA reporter. Subsequent liberation of the fluorophore from quencher can be visualized with a fluorimeter. Because the trans-cleavage activity of Casl2 is highly catalytic, sub-attomolar concentrations of target DNA can be detected (Teng et al., Genome Biol. 2019; 20: 132).

[00472] The collateral cleavage activity of Cast 2 can be harnessed for the sensitive detection of 5mC (FIG. 34). The workflow involves: 1) analyte DNA denaturation, 2) 5mC to T conversion using an altered cytidine deaminase, and 3) application of Casl2-guide RNA complex to the sample. The guide RNA is fully complementary to 5mC to T converted DNA, therefore only allowing Casl2 to engage analyte DNA if it contained 5mC. The ensuing activation of Casl2-mediated trans-cleavage results in cleavage of reporter ssDNA and increased fluorescence. This effect is measured using standard plate reader instruments.

[00473] Casl2-family enzymes exhibit different requirements for target DNA engagement, being specific for either dsDNA or ssDNA. The use of ssDNA-specific Casl2b or Casl2f is advantageous for the above workflow as it would bypass the requirement for target DNA rehybridization prior to Casl2-gRNA application. Casl2b/f-gRNA and APOBEC could then be combined in a one-pot reaction, as Casl2b/f will engage the target ssDNA after APOBEC-mediated 5mC to T conversion has occurred. This will accelerate the development of this workflow as a point-of-care (POC) diagnostic for DNA methylation.

[00474] Example 23

[00475J Spatial detection of 5mC using Fluorescence In Situ Hybridization (FISH)

[00476] FISH is routinely employed to detect the abundance and spatial localization of DNA sequences of interest in fixed cell and tissue sections using fluorescently labeled probes. It is increasingly used in IVD applications to diagnose cytogenic abnormalities such as chromosome translocations, deletions, and copy number variation. For example, deletion of chromosome 5q is associated with poor prognosis in acute myeloid leukemia, and has thus been assayed as a biomarker using FISH probes (molecular, abbott/us/ en/products/ oncology/vy si s-egr 1 -fi sh-probe-kit).

[00477] While the detection of DNA sequences by FISH is well-established, extending this technique to detect DNA modifications such as 5mC remains challenging due to technical difficulties with in situ bisulfite treatment of cells or tissue sections. Altered cytidine deaminase-mediated 5mC to T conversion enables such a workflow due to the mild conditions required for APOBEC enzymatic activity: for instance, 37°C at pH 5.2 - 6.5. Furthermore, we have discovered that the 5mC-selective APOBEC3A(Y130A/Y132H) enzyme is active in citrate buffer pH 6.0 (FIG. 35A, 35B, and 35C) buffer is commonly used in FISH protocols to increase overall probe signal intensity (Yu et al., Exp Ther Med. 2021; 22:1480). The envisaged APOBEC -FISH protocol involves: 1) permeabilization and denaturation of fixed cells or tissue sections, 2) 5mC to T conversion using APOBEC, and 3) hybridization of probes that are specific to the deaminated DNA sequence (FIG. 36).

[00478] Example 24

[00479] Array-based detection of 5mC

[00480] ILMN methylation array products typically involve the selective hybridization of bisulfite- converted DNA to bead-based probes. Using an altered cytidine deaminase for direct 5mC to T conversion lowers DNA input requirements and simplifies probe design as overall 4- base genome complexity is maintained.

[00481] In the case of 4 base genomes, array probes are more specific, so quantification of methylated samples is expected to be accurate than quantification of unmethylated samples. Hybridization efficiency for 4-base genome is better than degenerated lower complexity genome due to the preservation of complexity between the DNA and the probe, and can produce better resolutions.

[00482] Example 25

[00483] Expression of a CpG binding protein in soluble form

[00484] In this example, the coding region for a CXXC4 protein is cloned into a vector for recombinant protein expression. Suitable CXCC4 proteins include CXXC4 from H. sapiens (UniProt:K7EKZ6, SEQ ID NO:63, FIG. 37K). Multiple versions of this vector are created in which CXXC4 is fused to a soluble partner, including a SUMO tag, maltose binding protein (MBP) tag, or glutathione S-transferase (GST) tag. Vectors are created with each tag attached to the C-terminus, the N-terminus, or both termini of the CXXC4. Each protein fusion is expressed and purified. Yield, stability, and solubility of each protein fusion is measured to determine which soluble partners and fusion positions are compatible with soluble expression of CXXC4.

[00485] Example 26

[00486] Determination of the dissociation constant for CpG regions

[00487] In this example, the affinity for CpG regions of each fusion protein expressed and purified in Example 8 is measured using fluorescence polarization. Each fusion protein is provided with double-stranded DNA (dsDNA) or single-stranded DNA (ssDNA) including CpG regions with either 5mC or unmethylated C (a total of four possible substrates). The ssDNA and dsDNA substrates are fluorescently labeled.

[00488] The affinity of each fusion protein for each substrate is measured using fluorescence polarization. The fluorescent label on each DNA substrate exhibits an increase in polarization upon binding to a fusion protein. Higher affinity binding to a fusion protein results in less free rotation of the DNA substrate, which is measured as an increase in fluorescence polarization. The degree of fluorescence polarization is used to calculate a dissociation constant (Kd) for each substrate-fusion protein combination.

[00489] From this experiment, it is learned whether certain fusion proteins display preference for ssDNA or dsDNA as a substrate. It is also learned whether fusion proteins display preference for CpG regions including 5mC or unmethylated C. Fusion proteins that display high affinity and a preference for CpG regions including unmethylated C in ssDNA are determined to be good candidates for use with methods that include a cytidine deaminase or a TET protein.

[00490] Example 27

[00491] Test of CpG binding proteins in methods that include a cytidine deaminase or a TET protein

[00492] In this example, each fusion protein identified as a good candidate for in methods that include a cytidine deaminase or a TET protein in Example 9 are tested in either a cytidine deaminase or a TET protein.

[00493] Samples including target nucleic acids with CpG regions with unmethylated C, 5mC, or both are prepared. Each fusion protein is added to the samples before treatment with cytidine deaminase or TET and incubated to allow binding of the fusion protein to CpG regions. Samples with increasing concentrations of each fusion protein are tested. Sample preparation proceeds and samples are sequenced using standard protocols for each method.

[00494] The level of 5mC conversion to 5fC and 5caC and the level of unmethylated C to U conversion in samples with each fusion protein at each concentration are determined.

[00495] For each fusion protein, the concentration of protein required to maximally decrease unmethylated C conversion to U is determined. As conversion of 5mC to 5fC and 5caC is desired in some embodiments, fusion protein concentrations with decreased levels of 5mC conversion are determined to be undesirable. The highest measured level of 5mC conversion with minimal unmethylated C conversion of each fusion protein is compared. The fusion protein that decreased unmethylated C conversion the most with the lowest impact on 5mC conversion is determined to be the best candidate.

[00496] The complete disclosure of all patents, patent applications, and publications, and electronically available material (including, for instance, nucleotide sequence submissions in, e.g., GenBank and RefSeq, and amino acid sequence submissions in, e.g., SwissProt, PIR, PRF, PDB, and translations from annotated coding regions in GenBank and RefSeq) cited herein are incorporated by reference in their entirety. Supplementary materials referenced in publications (such as supplementary tables, supplementary figures, supplementary materials and methods, and/or supplementary experimental data) are likewise incorporated by reference in their entirety. In the event that any inconsistency exists between the disclosure of the present application and the disclosure(s) of any document incorporated herein by reference, the disclosure of the present application shall govern. The foregoing detailed description and examples have been given for clarity of understanding only. No unnecessary limitations are to be understood therefrom. The disclosure is not limited to the exact details shown and described, for variations obvious to one skilled in the art will be included within the disclosure defined by the claims.

[00497] Unless otherwise indicated, all numbers expressing quantities of components, molecular weights, and so forth used in the specification and claims are to be understood as being modified in all instances by the term "about." Accordingly, unless otherwise indicated to the contrary, the numerical parameters set forth in the specification and claims are approximations that may vary depending upon the desired properties sought to be obtained by the present disclosure. At the very least, and not as an attempt to limit the doctrine of equivalents to the scope of the claims, each numerical parameter should at least be construed in light of the number of reported significant digits and by applying ordinary rounding techniques.

[00498] Notwithstanding that the numerical ranges and parameters setting forth the broad scope of the disclosure are approximations, the numerical values set forth in the specific examples are reported as precisely as possible. All numerical values, however, inherently contain a range necessarily resulting from the standard deviation found in their respective testing measurements.

[00499] All headings are for the convenience of the reader and should not be used to limit the meaning of the text that follows the heading, unless so specified.