Patent Description:
The present filing includes nucleotide/amino acid sequences in form of a sequence listing.

Epigenetics is the study of heritable changes in gene expression that do not involve changes to the underlying primary nucleic acid sequence. Epigenetic change is a regular and natural occurrence but can also be influenced by several factors including age, the environment/lifestyle, and disease state. There is a renewed interest in epigenetics, as epigenic modifications have been associated with a host of disorders including various cancers, mental retardation associated disorders, immune disorders, neuropsychiatric disorders and pediatric disorders.

Epigenetic marks include DNA methylation, histone modifications, and regulatory RNA. Epigenetics is most often studied in multicellular organisms such as humans and plants, where epigenetic marks function in embryogenesis, cellular differentiation, genomic imprinting, and play roles in pathogenesis of diseases such as cancer. In bacteria, histones do not exist; therefore, non-coding RNAs and DNA methylation are the only universal epigenetic marks. The roles of DNA methylation in bacteria range from defense against bacteriophage infection, initiation of DNA replication, DNA repair, and gene regulation. 5N methylcytosine (5mC) is the most common epigenetic mark in higher eukaryotes, whereas 6N methyladenosine (6mA) is the most common epigenetic mark in bacteria. However, both marks are present in eukaryotes and prokaryotes, along with other nucleotide base modifications (See <FIG>). Direct evidence has linked 6mA to several species of eubacteria, mosquitoes, wheat, protists, and indirect evidence links 6mA to several archaea and vertebrates. However, methods to determine the global 6mA patterns using common next generation sequencing (NGS) techniques are not commercially available. Accordingly, the role of 6mA in eukaryotic genomic function may be under appreciated due to a failure to detect 6mA.

In prokaryotes, 6mA not associated with restriction modification systems is produced by the action of two methyl transferases, deoxyadenosine methylase (DAM), which methylates at GATC sites, and DNA methylase N-<NUM>/N-<NUM> domain-containing protein (CcrM), which methylates at the motif GANTC. Studies have shown that 6mA is present in the mitochondrial and chloroplast genomes, and therefore, could play a regulatory role throughout the Tree of Life.

<CIT> describes methods for detecting the presence of methylation at a locus within a population of nucleic acids. <NPL>) describes the characterization of eukaryotic DNA N<NUM>-methyladenine by a highly sensitive restriction enzyme-assisted sequencing. While 5mC can be measured by conventional bisulfite sequencing, global 6mA patterns at the nucleotide level are currently detected only by PacBio Sequencing (<NPL>), requiring a specialized, costly instrument and significant consumables to reach the depth of sequencing necessary for methylated base calls i.e., -$<NUM>,<NUM> for a human genome). Thus there is a need for a new methodology that allows for the rapid and cost effective analysis of the presence of 6mA in a sample of genomic DNA.

In accordance with one embodiment of the present disclosure, applicant provides a novel method, "6mA-Seq", to identify 6mA residues by sequence analysis using an Illumina sequencer or equivalent equipment.

Currently, DNA methylation is one of the most broadly studied and well-characterized epigenetic modifications. Bacteria, like eukaryotes, use methylation to regulate gene expression. However, methylation profiling is not common in bacteria due to a lack of methodology pertinent to study of the dominant epigenetic marker in prokaryotes, 6N methyladenine (<NUM> mA). In accordance with one embodiment of the present disclosure a method is provided for assessing global <NUM> mA profiles in bacterial genomes. The method can be used to monitor changes in methylation patterns of bacterial genomic DNA over time and/or in response to various environmental factors, including for example temperature, availability of nutrients and presence of stimulants or toxins.

In accordance with one embodiment a method is provided for monitoring global methylation patterns in genomic DNAs recovered from organisms or cell populations. More particularly, the method is directed to assessing the methylated state of genomic sequences associated with one or more target restriction enzyme recognition sites as further specified in the appended claims. In one embodiment the method of detecting methylated nucleic acid residues comprises the steps of first obtaining a library of genomic DNA and then subjecting the library of genomic DNA to restriction enzymatic digestion with either <NUM>) enzymes that cut only at sites when the DNA is methylated (leaving a library enriched for unmethylated regions) or enzymes that cut only unmethylated regions (leaving a library enriched for methylated sequences). The DNA sequences of the restriction enzyme cleaved library are then analyzed using Next Generation Sequencing (NGS) to determine the sequences remaining in the restriction enzyme cleaved library, wherein comparison of the sequences identified by the Next Generation Sequencing step to a reference library of all available target restriction enzyme recognition sites present in relevant genome reveals those sites that were methylated in the analyzed sample.

In one embodiment the method of detecting methylated nucleic acid residues in a sample comprising genomic DNA comprises the steps of first obtaining a library of genomic DNAs and then contacting the library with a methylation specific restriction enzyme that cleaves said target nucleic acid recognition site only when said target nucleic acid recognition site is unmethylated, to produce a set of digested genomic nucleic acid molecules enriched in methylated sequences. In one embodiment the methylation sensitive restriction enzyme is selected from the restriction enzymes listed in <FIG>. In one embodiment the restriction enzyme used to digest the genomic DNA is selected from the group consisting of DpnI, DpnII and MboI.

The DNA sequences of the restriction enzyme cleaved library is then analyzed using Next Generation Sequencing to determine the sequences remaining in the restriction enzyme cleaved library. Comparison of the sequences identified by the Next Generation Sequencing step with a reference genomic sequences (said reference genomic sequence comprising all of the respective genomic sequence comprising the the restriction enzyme recognition site) reveals the unmethylated sequence in the analyzed genome as sequences missing relative in the restriction enzyme cleaved library relative to the reference genomic sequences. The sequences detected in the restriction enzyme cleaved library that comprise a target restriction enzyme recognition site are identified as methylated sequences.

In accordance with one embodiment the genomic DNA to be analyzed for the presence of methylated nucleotides (e.g., 6mA) is processed prior to enzymatic directions. More particularly, a PCR-Free library of genomic DNA is prepared suitable for Next Generation Sequencing. In one embodiment the preparation of the PCR-free library comprises the steps of isolating genomic DNA from a cell/ organism without an amplification step, fragmenting the isolated genomic DNA into DNA sequences less than <NUM> Kb in length, and typically having an average size of about <NUM> bp to about <NUM> bp, and the ligating the fragments of the genomic DNA to adapters that include all necessary components for sequencing primer annealing and attachment to a flow-cell surface for conducting Next-Generation Sequencing.

In one embodiment the step of restriction enzyme digestion of the genomic DNA is conducted using two different restriction enzymes that cleave the same recognition sequence but where one enzyme is sensitive to methylation and the other is not. The initial library of genomic DNA is divided into a first and second pool of genomic DNA wherein the first pool of genomic DNA is digested with a restriction enzyme that cleaves its target nucleic acid recognition site only when said target nucleic acid recognition site is unmethylated to produce a first set of digested genomic nucleic acid molecules (enriched in methylated DNAs), and second pool of genomic DNA is digested with a restriction enzyme that cleaves its target nucleic acid recognition site regardless of the methylated state of the target recognition site to produce a second set of digested genomic nucleic acid molecules. The sequence of the digested nucleic acids of the first and second digested genomic nucleic acid molecules is then determined, typically by using Next Generation Sequencing techniques. The nucleic acid sequences of the first and second digested genomic nucleic acid molecules are compared, and sequences present in the first set of digested genomic nucleic acid molecules relative to the second set of digested genomic nucleic acid molecules are identified as methylated sequences. In one embodiment the restriction enzyme used to digest the first pool of genomic DNA is selected from the group consisting of DpnI, DpnII and MboI and the restriction enzyme used to digest the second pool of genomic DNA is Sau3AI.

In describing and claiming the invention, the following terminology will be used in accordance with the definitions set forth below.

The term "about" as used herein means greater or lesser than the value or range of values stated by <NUM> percent, but is not intended to designate any value or range of values to only this broader definition. Each value or range of values preceded by the term "about" is also intended to encompass the embodiment of the stated absolute value or range of values.

As used herein, the term "purified" and like terms relate to the isolation of a molecule or compound in a form that is substantially free of contaminants normally associated with the molecule or compound in a native or natural environment. As used herein, the term "purified" does not require absolute purity; rather, it is intended as a relative definition. The term "purified nucleic acid" is used herein to describe a nucleic acid which has been separated from other compounds including, but not limited to polypeptides, lipids and carbohydrates.

The term "isolated" requires that the referenced material be removed from its original environment (e.g., the natural environment if it is naturally occurring). For example, a naturally-occurring nucleic acid present in a living animal is not isolated, but the same nucleic acid, separated from some or all of the coexisting materials in the natural system, is isolated.

As used herein the terms "restriction endonuclease" or "restriction enzyme" are used interchangeably and encompass proteins that are able to cleave a double stranded DNA sequence at or near a specific sequence of nucleotides.

As used herein the term "restriction enzyme recognition site" defines locations on a DNA molecule containing a specific sequence of nucleotides that are recognized by the individual restriction enzyme and result in the cleavage of the sequence between two nucleotides within its recognition site, or somewhere nearby.

As used herein the term "methylation sensitive restriction enzyme" encompasses restriction enzymes whose cleavage is blocked or inhibited when the restriction enzyme recognition site is methylated by the cognate methylase.

In accordance with one embodiment, the present disclosure is directed to a method for analyzing global methylation patterns in genomes, particularly in bacterial genomes or the genomes of eukaryotic organelles as further specified in the appended claims. Restriction enzymes that select for digestion of DNA at sites with differential methylation enable one to remove either unmethylated or methylated DNA and then compare the library that remains to a full genome to identify regions of methylation on a genome wide basis. Such methods can be used to establish correlations between certain methylation patterns and disease states or conditions, and/or determine the impact of various environmental factors on the methylated state of the genome and their corresponding impact on gene expression.

Restriction endonucleases are known that will only cleave their target recognition site when the DNA is in a unmethylated state. Advantageously, as disclosed herein such restriction enzymes can be used to detect the presence of methylated nucleic acids in genomic sequences, and more significantly analyze the methylated state of DNA sequences at genomic level.

<FIG> provides a list of restriction enzymes that are methylation sensitive. Restriction enzyme cleavage is blocked or substantially inhibited when the recognition sequence is methylated by the cognate methylase. More particularly, methylation of nucleic acid bases at or near the restriction enzyme recognition site can block cleavage, leave cleavage unaffected, or slow the rate or extent of cleavage. A restriction enzyme database, REBASE, is known to those skilled in the art for providing more detailed information regarding methylation sensitive restriction enzymes.

In accordance with the present disclosure, methylation sensitive restriction enzymes can be used to determine and monitor methylation patterns of genomic DNAs on a global level. In particular, one can track how methylation of genomic DNA sequences is altered by exposure to various environmental factors or genetic background. The methods disclosed herein can be used to analyze genomic DNA isolated from any cell, including organelle genomic DNA from chloroplasts and/or mitochondria of eukaryotic cells. In one embodiment the methylation state of genomic DNA isolated from bacterial cells is analyzed using the methods disclosed herein.

There are three common DNA methylases that are present in bacterial cells such as E. coli: the Dam methylase, Dcm methylase and EcoKI methylase. The methylase encoded by the dam gene (Dam methylase) transfers a methyl group from S-adenosylmethionine (SAM) to the N6 position of the adenine residues in the sequence GATC. The Dcm methylase (encoded by the dcm gene; referred to as the Mec methylase in earlier references) methylates the internal (second) cytosine residues in the sequences CCAGG and CCTGG at the C5 position. The EcoKI methylase, M. EcoKI, modifies adenine residues in the sequences AAC(N6)GTGC and GCAC(N6)GTT. EcoKI sites (~<NUM> site per <NUM> kb) are much less common than Dam sites (~<NUM> site per <NUM> bp) or Dcm sites (~<NUM> site per <NUM> bp) in DNA of random sequence (GC=AT).

In one embodiment genomic DNA is isolated from a cell and subjected to enzymatic digestion using a methylation sensitive endonuclease. In one embodiment the isolated genomic DNA is first processed and a library of the genomic DNA is prepared from the processed genomic DNA. In accordance with one embodiment, the genomic DNA is first fractionated to reduce the average size of the genomic DNA prior to preparation of the library. The genomic DNA can be fractionated using any standard technique known to those skilled in the art to reduce the size of the average genomic DNA to less than <NUM> Kb in length. In one embodiment the DNA is fractionated by ultasonification and/or nebulization, including for example the use of the Covaris Adaptive Focused Acoustics (AFA) technology (Covaris, Inc. , <NUM> Gill Street, Unit H, Woburn, Massachusetts, <NUM>-<NUM> USA) to generate fragments having an average size of about <NUM> bp to about <NUM> bp, or about <NUM> bp to about <NUM> bp. Ultasonification shearing can be used to generate double-stranded DNA (dsDNA) fragments with <NUM>' or <NUM>' overhangs (see <CIT>).

The fractionated DNA is then linked to the appropriate adapters to create a library of fractionated genomic DNA sequences. Preferably the creation of the genomic libraries is conducted in the absence of a PCR amplification step (i.e., a PCR-free library). In accordance with one embodiment the fragmented genomic DNA is subjected to a repair step wherein the overhangs resulting from fragmentation are converted into blunt ends. In one embodiment a <NUM>' to <NUM>' exonuclease activity removes the <NUM>' overhangs, and a <NUM>' to <NUM>' polymerase activity completes the <NUM>' overhangs. In a final step the fragments of genomic DNA are ligated to adapters that include all necessary components for sequencing primer annealing and attachment to a flow-cell surface for conducting Next-Generation Sequencing. See <NPL>; and Illumina TruSeq DNA PCR-Free Illumina, Inc. <NUM> Illumina Way, San Diego, CA <NUM> USA).

In accordance with one embodiment PCR-free genomic libraries are prepared from the cells whose genomic DNA will be assessed for methylation patterns. PCR amplification is commonly used in generating libraries for Next-Generation Sequencing (NGS) to efficiently enrich and amplify sequenceable DNA fragments. However, it introduces bias in the representation of the original complex template DNA. Accordingly, in one embodiment of the present invention libraries of genomic DNA will be prepared in the absence of a PCR amplification step and the digestion with methylation sensitive restriction enzymes will be conducted on the PCR-free library components.

In one embodiment a method is provided for analyzing 6mA methylation patterns in genomic DNA, including genomic DNA isolated from prokaryotics. The method comprises obtaining a PCR-Free library of genomic sequences and subjecting the library to at least one restriction enzyme digestion (cleaving DNA but only when the recognition site is unmethylated (e.g., lacking a 6mA residue). Optionally a second enzymatic digestion can be conducted wherein the second enzyme cleaves the same sites as the first but irregardless of whether the site is methylated or unmethylated. The digested library can optionally be amplified by PCR to enrich for uncut sequences. The enzymatic digested library sequence can be sequenced (typically by NGS) and compared to the corresponding reference sequences (comprising all known target recognition sites for the relevant genome) to determine where methylation was or was not, depending on the enzyme used. In one embodiment this analysis is conducted on bacterial DNA or other DNA to assess N6-methyladenosine patterns over time and/or and in relation to exposure to different environmental factors.

Bacteria primarily use 6mA for epigenetic regulation of gene expression. In accordance with one embodiment methods are provided for detecting 6mA in genomic DNA using Illumina platforms. In one embodiment, whole genome sequencing libraries will be prepared, polymerase chain reaction-free, from DNA from each organism or cell line. These libraries are treated with a mix of enzymes that remove library molecules that contain methylated adenosines (such as 6mA), leaving behind only molecules that do not have a methylated adenosine. The libraries will then be amplified to enrich the 6mA libraries for complete molecules and sequenced on the appropriate Illumina sequencing platform. Adapters and low-quality reads will be trimmed with Trimmomatic software prior to analysis of the data. Data analysis includes alignment to the reference sequence by BWA-MEM software and visualization through IGV genome viewer. A bed file of each genome will be produced with the coordinates of possible adenosine methylation sites. Coverage at these possible methylation sites will be determined and compared to the average coverage across the genome. Regions where reads are absent or more than two standard deviations below the mean in coverage represent regions where adenosine bases were not methylated. Analysis of methylated versus unmethylated sites in this manner will provide a global map of 6mA across the genome.

A list of methylation sensitive and insensitive restriction enzymes is provided in <FIG>. In accordance with one embodiment the methylation sensitive restriction enzyme used in the present invention is selected from the group consisting of DpnI, DpnII and MboI. MboI, DpnI and DpnII each recognize and cleave at the recognition sequence GATC and both enzymes are blocked by dam methylation, but not blocked by dcm methylation. HinfI recognizes and cleaves at the recognition sequence GANTC, and cleavage is not blocked by either dam methylation, or dcm methylation. HpaII and MspI recognize and cleave at the recognition sequence CCGG recognition site and can be used for C5-methylcytosine detection.

The global analysis of 6mA in species beyond prokaryotes will revolutionize the field of epigenetics, demonstrating the impact of 6mA in response to environmental changes. State-of-the-art analysis methods overlook 6mA, although it is known to regulate processes in both prokaryotes and eukaryotes. It is now apparent that 6mA is an important regulator in all organisms but requires a lower limit of detection than possible back in the <NUM>. The 6mA modification can regulate gene expression in specific organelles, including mitochondria which is targeted by environmental contaminants such as PFOA and causes a large impact on the lives of multicellular eukaryotes.

In accordance with the present disclosure a kit is provided for analyzing DNA for the presence of 6mA. The kit comprises one or more methylation sensitive restriction enzymes, optionally selected from the group consisting of DpnI, DpnII and MboI. The kit may contain additional reagents for preparing PCR-free libraries.

In accordance with a first embodiment, a method of detecting methylated nucleic acid residues, optionally a N6-methyladenosine (6mA) modification, in one or more target nucleic acid recognition sites present in a sample of genomic DNA is provided, wherein said method comprising the steps of.

In accordance with another embodiment the library of genomic DNA of the first embodiment is prepared by isolating genomic DNA from a cell without a PCR amplification step;.

In accordance with another embodiment, the method of the above embodiments is provided wherein the library of genomic DNA is contacted with a methylation specific restriction enzyme selected from the group consisting of DpnI, DpnII and MboI, optionally the library of genomic DNA is contacted with two or more methylation sensitive restriction enzymes.

In accordance with another embodiment, the method of any one of the above embodiments is provided wherein the genomic DNA is isolated from a prokaryote.

In accordance with another embodiment, the method of any one of the above embodiments is provided wherein the analyzing step comprises determining the nucleic acid sequence of said digested genomic nucleic acid molecules and comparing those nucleic acid sequences to a reference set of nucleic acids that represent all available target nucleic acid recognition sites present in said sample of genomic DNA, wherein sequences missing from the digested genomic nucleic acid molecules relative to said reference set represents a unmethylated target sequence in said sample of genomic DNA.

In accordance with a second embodiment, a method of detecting methylated nucleic acid residues, optionally a N6-methyladenosine (6mA) modification, in one or more target nucleic acid recognition sites present in genomic DNA is provided wherein, said method comprising the steps of.

In accordance with another embodiment, the method of the second embodiment is provided wherein the methylation specific restriction enzyme is selected from the group consisting of DpnI, DpnII and MboI.

In accordance with another embodiment, the method of the second embodiments above is provided wherein the genomic DNA fragments further comprise sequences complementary to sequences linked to a solid support.

In accordance with another embodiment, the method of any one of the second embodiments above is provided wherein the library of genomic DNA is contacted with two or more methylation sensitive restriction enzymes.

In accordance with another embodiment, the method of any one of the second embodiments above is provided wherein the genomic DNA is isolated from a prokaryote.

In accordance with another embodiment, the method of any one of the second embodiments above is provided wherein the analyzing step comprises determining the nucleic acid sequence of said digested genomic nucleic acid molecules and comparing those nucleic acid sequences to a reference set of nucleic acids that represent all available target nucleic acid recognition sites present in said sample of genomic DNA, wherein sequences missing from the digested genomic nucleic acid molecules relative to said reference set represents a unmethylated target sequence in said sample of genomic DNA.

In accordance with another embodiment, the method of any one of the second embodiments above is provided wherein the library of genomic DNA of step a) is divided into a first and second pool of genomic DNA, and said method comprises the steps of.

In accordance with another embodiment, the method of any one of the second embodiments above is provided wherein the target nucleic acid recognition site comprises a nucleic acid sequence of GATC or GANTC, optionally wherein the target nucleic acid recognition site comprises a nucleic acid sequence of GATC.

In accordance with a third embodiment, a method of detecting genomic sequences comprising N6-methyladenosine (6mA) is provided wherein said method comprises the steps of.

In accordance with another embodiment, the method of the third embodiment is provided wherein the genomic DNA of said library is contacted with two or more methylation sensitive restriction enzymes.

In accordance with another embodiment, the method of any one of the third embodiments above is provided wherein the methylation sensitive restriction enzyme is selected from the group consisting of MboI, DpnI and DpnII.

In accordance with another embodiment, the method of any one of the third embodiments above is provided wherein the genomic DNA is prokaryotic DNA.

In accordance with another embodiment, the method of any one of the third embodiments above is provided wherein the library of genomic DNA fragments is prepared by.

In accordance with another embodiment, the method of any one of the third embodiments above is provided wherein the library of genomic DNA of step a) is divided into a first and second pool of genomic DNA, and said method comprises the steps of.

In accordance with another embodiment, the method of any one of the third embodiments above is provided wherein the methylation sensitive restriction enzyme is selected from the group consisting of DpnI, DpnII and MboI.

In accordance with another embodiment, the method of any one of the third embodiments above is provided wherein the first restriction enzyme is selected from the group consisting of DpnI, DpnII and MboI and the second restriction enzyme is Sau3AI.

A proof of concept study was conducted using two strains of Escherichia coli. One strain was derived from E. coli K-<NUM> with the genotype dam-/dcm-, causing it to be unable to methylate adenosine or cytosine. The other strain was E. coli ATCC <NUM>, which is wild-type dam+/dcm+, enabling it to methylate its genome at both cytosine and adenosine. PCR-free Illumina sequencing libraries were prepared from DNA from both strains. These libraries were then treated with enzymes that remove library sequences that contained unmethylated adenosines leaving behind only sequences that contained potential 6mA sites. The dam-/dcm- strain had <NUM> possible 6mA sites of which none were methylated according to the 6mA sequencing analysis, whereas the dam+/dcm+ strain had <NUM>% of its <NUM> possible 6mA sites methylated (Table <NUM>). This data demonstrates the utility of the 6mA-Seq method for detection of differential 6mA.

This method does not require any specialized equipment, utilizes a protocol that is very similar to whole genome sequencing on the Illumina platform, and costs the same as whole genome sequencing. Therefore, the method is easy to use by any lab with access to an Illumina sequencer. The current method should be expandable to detect differential 6mA methylomes across the Tree of Life as the enzymes used to target 6mA will work on any 6mA methylated DNA with the same efficiency.

Two Escherichia coli strains, E coli ATCC <NUM> with both dam and dcm activity and E coli K <NUM> with genotype dam-/dcm- were cultured in different media. DNA from each strain was extracted for further analysis using methylation specific restriction enzymes and next generation sequencing. Data analysis consisted in evaluating methylation patterns to find unmethylated locations as absence of sequencing reads (See Table <NUM>).

<NUM> GATC sites demonstrated hyper or hypomethylation in E. coli ATCC25922. <NUM> GATC loci had consistent hypomethylation in all four growth conditions (position start: <NUM>, <NUM>, and <NUM>), whereas <NUM> were consistently hypermethylated (position start: <NUM> and <NUM>). <NUM> GATC sites (position start: <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, and <NUM>) were consistently hypermethylated in M9 minimal media regardless of carbon source and were not in LB. <NUM> GATC sites were hypomethylated only in LB (position start: <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, and <NUM>). LB shared <NUM> hypomethylated GATC sites with M9-glycerol and M9-glucose (position start: <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, and <NUM>).

The data demonstrates the utility and high sensitivity of <NUM> mA Seq for detection of differentially methylated loci in bacteria. While this method is limited by the identification of methylation at the restriction enzyme motifs, the analysis can be expanded by the use of additional methylation sensitive enzymes to digest the library of genomic sequences, thus analyzing methylation at motifs beyond GATC, which will increase the diversity of epigenetic signatures that can be identified by this technique.

As demonstrated <NUM> mA Seq analytical technique readily identifies sites with differential methylation patterns. E coli which cannot methylate adenine has no coverage at GATC sites, whereas E coli with DAM methylase has <NUM>% of sites methylated
<NUM> mA hyper and hypomethylation status changes in different growth conditions at loci across the genome.

LB grown E coli cultures are very different in methylation pattern from those in M9 minimal media, with only <NUM> GATC sites consistently hyper or hypomethylated between all conditions tested. The <NUM> hypomethylated sites correspond to a biosynthetic threonine ammonia lyase, a tRNA Ala and a location approximately <NUM> bp from this tRNA Ala. The two hypermethylated sites correspond to IS <NUM> family transposase and
ClbS/DfsB family four helix bundle protein.

In each culture condition there were unique E coli hyper or hypomethylated sites compared to other conditions tested.

Claim 1:
A method of detecting methylated nucleic acid residues in one or more target nucleic acid recognition sites present in genomic DNA, said method comprising the steps of
a) obtaining a library of genomic DNA fragments, wherein said genomic DNA fragments comprise isolated genomic DNAs that have not been subjected to PCR, have been fragmented to have an average size of about <NUM> bp to about <NUM> bp and have been further modified by covalent linkage of an adapter sequences to each DNA fragment wherein the adapter sequence comprises a DNA sequencing primer;
b) contacting said library with a methylation specific restriction enzyme that cleaves said target nucleic acid recognition site only when said target nucleic acid recognition site is unmethylated to produce a restriction enzyme cleaved library comprising a set of digested genomic nucleic acid molecules and uncut methylated nucleic acid molecules; and
c) analyzing the restriction enzyme cleaved library using next generation sequencing to determine the nucleic acid sequence of said uncut methylated genomic nucleic acid molecules and comparing those nucleic acid sequences to a reference set of nucleic acids that represent all available target nucleic acid recognition sites present in said library of genomic DNA fragments wherein sequences missing from the restriction enzyme cleaved library relative to said reference set of nucleic acids represent unmethylated target sequence in said sample of genomic DNA.