Patent Description:
There is growing appreciation for the profound ways in which the human microbiome can impact our health, but the comprehensive characterization of these microbial populations remains difficult. Amplicon sequencing of the <NUM> rRNA gene provides a culture-free means of identifying many of the taxa present in a metagenomic sample, but the phylogenetic resolution of this technique is limited and the microbial genomic architecture outside of this single gene is left unexamined or only inferred indirectly. Whole metagenome shotgun sequencing provides access to all the genomic features of the constituent organisms, including bacterial and archaeal chromosomes, plasmids, transposons, and even bacteriophages with a phylogenetic resolution extending up to the strain level. However, multiple technical challenges hinder the interpretation of metagenomic sequencing data collected by short read next-generation sequencing (NGS) methods.

NGS data typically consists of millions of reads that are <<NUM> bp in length, providing considerable depth of sequencing but limited ability to resolve both complex repeats and similar sequences that exist in multiple genomes. This presents significant challenges for de novo metagenomic assembly and interpretation of the resulting thousands of small assembled sequences (called contigs), which relies heavily on either reference-based annotation methods or segregation into putative taxa through a process known as metagenomic binning. Unsupervised (reference-free) methods have the potential to identify novel species, unlike supervised binning methods that require existing references to train classification algorithms. Several reference-free methods attempt to bin metagenomic reads prior to de novo assembly by using k-mer frequency metrics to assess sequence composition profiles or by tracking k-mer covariance across multiple samples. These methods do not depend on the results of a de novo assembly, but the binning resolution is limited by the information content found in short reads from standard NGS technologies.

Owing to the limited information content in short reads, most reference-free binning methods instead utilize the longer sequences of assembled contigs. Composition-based contig binning approaches not only rely on a successful de novo assembly, but also often fail to segregate sequences when the sample contains multiple high-similarity bacterial genomes. Differential coverage (or coverage covariance) methods, which partition sequences based on their similar abundance profiles across multiple samples, provide a powerful means of binning sequences in projects studying a large number of complex samples. However, they sometimes fail to untangle genomes of organisms that share similar abundances across samples and cannot effectively bin independently replicating mobile genetic elements (MGE), such as plasmids, transposons, bacteriophages, and Group I and II introns, which can have dramatically different abundance levels from their host chromosome(s). An alternative approach involves using Hi-C chromosomal interaction maps to link assembled contigs, including MGEs, but these methods are also limited by difficulties in distinguishing between closely related organisms due to high sequence similarity and uneven Hi-C link densities.

The information content of DNA is not limited to the primary nucleotide sequence (A, C, G and T), but is also conveyed by chemical modifications of individual nucleotides, including DNA methylation. In the bacterial (and archaeal) kingdom, DNA methylation is catalyzed by DNA methyltransferases (MTases) that apply methyl groups to DNA bases in a highly sequence-specific manner, causing certain sequence motifs to be nearly <NUM>% methylated while the other motifs remain non-methylated Single-molecule, real-time (SMRT) sequencing of native (amplification-free) DNA makes it possible to detect methylated bases and motifs in prokaryotic genomes. A recent survey of <NUM> diverse bacterial and archaeal genomes found DNA methylation in <NUM>% of genomes across a wide diversity of methylated motifs (<NUM> distinct motifs; averaging three motifs per organism). Importantly, the genetic contents of a cell (chromosomes and extrachromosomal DNA elements) all share the same set of methylation motifs. yet these motifs often differ dramatically across species and strains. The primary reason for such widespread diversity of methylated motifs is horizontal gene transfer (HGT) by mobile genetic elements. Since MTases are often shuttled by HGT, the process plays a crucial role in reconfiguring the bacterial methylomes. Additionally, mutation events can occur in the target recognition domain of MTase genes and thereby modify the sequence motif targeted for methylation, providing a route to further diversification of bacterial methylomes.

This raises the possibility of using SMRT sequencing to access DNA methylation in these communities, which essentially provides an orthogonal data dimension (endogenous epigenetic barcode) that can be leveraged for genome segregation in support of complementary features like coverage and sequence composition.

Whole metagenome shotgun sequencing is a comprehensive approach for characterizing complex microbial communities. However, significant challenges arise in the analysis of metagenomic sequences, often stemming from the presence of highly similar bacterial strains with varying relative abundances. Although a number of metagenomic binning methods have been developed that use features capturing sequence composition, organism abundance, and chromosome organization, many applications still suffer from insufficient discriminative power to distinguish among closely related species and strains with high sequence similarity Single-molecule long-read sequencing technologies enabled the comprehensive detection of DNA methylation events in bacteria, a rich dimension of discriminative features beyond DNA sequences that have not yet been exploited in metagenomic analyses.

The foregoing discussion is presented solely to provide a better understanding of nature of the problems confronting the art and should not be construed in any way as an admission as to prior art nor should the citation of any reference herein be construed as an admission that such reference constitutes "prior art" to the instant application. <CIT> concerns nucleic acid characteristics as guides for sequence assembly. <CIT> concerns methods and compositions for sequencing modified nucleic acids. <NPL>, page e0144612, concerns DNA Methylation Assessed by SMRT Sequencing Is Linked to Mutations in Neisseria meningitidis Isolates. <NPL>, concerns single molecule-level detection and long read-based phasing of epigenetic variations in bacterial methylomes.

A novel approach is presented for binning metagenomic sequences that leverages long reads from a single-molecule long-read sequencing technology and, for the first time, utilizes the DNA methylation signatures inferred from these reads to resolve individual reads and assembled contigs into species- and even strain-level clusters. This novel methylation-based binning approach also enables the mapping of mobile genetic elements (e.g., plasmids, transposons, including retrotransposons, DNA transposons, and insertion sequences, bacteriophages, group I introns, and group II introns) to their host species directly in a microbiome sample.

A novel approach is described to identify the DNA methylation patterns present in metagenomic data using read-level polymerase kinetics of SMRT reads and demonstrate how to exploit this data to derive a sequence-independent, endogenous epigenetic barcode that improves the resolution of metagenomic binning. Because the methylated motifs often differ even between closely-related species and strains, the methylation patterns (sets of motifs) present in SMRT reads and their assembled contigs offer a means for better differentiating sequences from taxonomical groups with high sequence similarity.

In one embodiment, an approach for organizing assembled contigs into taxon-specific clusters using DNA methylation profiles is described, and its complementarity with existing binning approaches that rely on sequence composition and coverage-covariance features is demonstrated.

To complement contig-level DNA methylation-based binning, an approach has been developed and applied to leverage the long read lengths of SMRT sequencing to directly bin individual single-molecule reads using sequence composition and DNA methylation profiles, facilitating the detection of low-abundance organisms and resolving multi-strain de novo assemblies into isolated single-strain assemblies.

In one aspect of the description, a method of deconvoluting genomes of prokaryotic organisms in a microbiome sample is provided, said method comprising the steps of:.

The method of the invention is defined in appended claim <NUM>.

In some embodiments, two or more of the prokaryotic organisms in the microbiome sample have high sequence similarity. In some embodiments, two or more of the prokaryotic organisms in the microbiome sample have an average nucleotide identity of greater than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, or than about <NUM>%.

Microbiome samples of any size or complexity are within the scope to be analyzed by the methods of the present invention. In one embodiment, the microbiome sample analyzed by the methods of the present invention comprises greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>, or greater than <NUM>,<NUM> prokaryotic host organisms.

In one embodiment the methylated nucleotides are selected from N<NUM>-methyladenine, N<NUM>-methylcytosine, and <NUM>-methylcytosine and combinations thereof.

Any prokaryotic organisms known to those skilled in the art are within the scope of the present invention. In one non-limiting embodiment, the prokaryotic organisms are bacterial organisms, archaeal organisms, and combinations thereof. In some non-limiting embodiments, the prokaryotic organisms are bacterial organisms, bacterial species, or strains of bacterial species. In other non-limiting embodiments, the prokaryotic organisms are archaeal organisms, archaeal species, or strains of archaeal species.

In some non-limiting embodiments, the bacterial organisms comprise organisms of bacterial orders Bacteroidales, Bacillales, Bifidobacteriales, Burkholderiales, Clostridiales, Cytophagales, Eggerthallales, Enterobacterales, Erysipelotrichales, Flavobacteriales, Lactobacillales, Rhizobiales, or Verrucomicrobiales, and combinations thereof.

In some non-limiting embodiments, the bacterial organisms are strains of Bacteroides dorei, Bacteroides fragilis, Bacteroides thetaiotaomicron, Bifidobacterium breve, Bifidobacterium longum, Alistipes finegoldii, or Alistipes shahii.

Microbiome samples analyzed by the methods of the invention can be obtained from any source known to those skilled in the art. In one non-limiting embodiment, the microbiome sample is obtained from soil, air, water (including, without limitation, marine water, fresh water, and rain water), sediment, oil, and combinations thereof. In another non-limiting embodiment, the microbiome sample is obtained from a subject selected from a protozoa, an animal (e.g., a mammal, e.g., human), or a plant. The subject (e.g., a mammal, e.g., a human) can be of any age (e.g., infant, child, adolescent, adult, or elderly.

In some embodiments, the subject is at a genetic risk for development a disease, e.g. diabetes mellitus, e.g., type I diabetes mellitus. In other embodiments, the subject may be at a risk of having, or have a bacterial infection, e.g., pneumonia infection.

Any single-molecule sequencing technology can be used in the methods of the present invention. In some embodiments, sequencing nucleic acids of the prokaryotic organisms is accomplished using a single-molecule real time (SMRT) technology or nanopore (e.g., Oxford Nanopore) sequencing technology.

In some embodiments of the above method, the nucleic acid methylation profile is a DNA methylation profile.

In some embodiments, the method described above comprises further steps. In one embodiment, the method described above further comprises the step of combining the methylation profiles of step (d) with other sequence features of the nucleic acids of the prokaryotic organisms in the microbiome sample prior to separating the assembled contigs and/or the single molecule reads into bins.

In one embodiment, the method described above comprises other sequence features, such as k-mer frequency profiles and coverage profiles across multiple samples.

In another embodiment, the method described above further comprises the step of combining contig binning assignments from other tools, such as cross-coverage and composition-based binning tools, with methylation scores in each bin, resulting in detection of methylated motifs in each bin and assignment of bin-level methylation scores in the microbiome sample.

In another embodiment, the method described above further comprises the step of aligning the single molecule reads to the contigs assembled from single molecule reads of the nucleic acids of step (a) prior to the step of assigning a methylation score.

Detailed embodiments of the present invention are disclosed herein; however, it is to be understood that the disclosed embodiments are merely illustrative of the invention that may be embodied in various forms, as defined in the appended claims. In addition, each of the examples given in connection with the various embodiments of the invention is intended to be illustrative, and not restrictive. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a representative basis for teaching one skilled in the art to variously employ the present disclosure.

As used in this specification and the appended claims, the singular forms "a", "an", and "the" include plural references unless the context clearly dictates otherwise. Thus, for example, a reference to "a method" includes one or more methods, and/or steps of the type described herein and/or which will become apparent to those persons skilled in the art upon reading this disclosure.

The terms "treat" or "treatment" of a state, disorder or condition include: (<NUM>) preventing, delaying, or reducing the incidence and/or likelihood of the appearance of at least one clinical or sub-clinical symptom of the state, disorder or condition developing in a subject that may be afflicted with or predisposed to the state, disorder or condition but does not yet experience or display clinical or subclinical symptoms of the state, disorder or condition; or (<NUM>) inhibiting the state, disorder or condition, i.e., arresting, reducing or delaying the development of the disease or a relapse thereof or at least one clinical or sub-clinical symptom thereof; or (<NUM>) relieving the disease, i.e., causing regression of the state, disorder or condition or at least one of its clinical or sub-clinical symptoms. The benefit to a subject to be treated is either statistically significant or at least perceptible to the patient or to the physician.

In one aspect of the description a methodology is provided that enables DNA methylation signatures in unamplified prokaryotic genomes to be profiled by SMRT sequencing and serve as endogenous epigenetic barcodes that present a rich, yet unexplored, dimension of discriminative features capable of providing high resolution metagenomic analyses.

In another aspect of the description, methylation profiles are exploited as a general discriminative feature to segregate assembled contigs, and this methodology is superior to existing methods based on sequence composition profiles and coverage covariance.

Furthermore, in yet another embodiment, it is disclosed how the proposed single molecule read-level binning of long SMRT reads can be used to address multiple challenges in metagenomic de novo assembly, such as assisting in the identification of low-abundance organisms and simplifying de novo metagenome assembly of multiple co-existing strains with high sequence similarity.

Sequence binning by DNA methylation profiles enables multiple other applications. First, methylation profiling can be a tool to track the transmission of plasmids and bacteriophages across geographical locations, time points or conditions, such as antibiotic treatment. Because the methylation signature of a plasmid or phage reflects the most recent bacterial host in which it replicated, transmission events can be detected by comparing the methylation profile of a specific plasmid or phage (and the bacterial community) between two conditions. Second, aside from serving as endogenous epigenetic barcodes for metagenomic binning, bacterial DNA methylation events also plays an important role in the regulation of gene expression and pathogenicity. While existing methods require a clonal sample for methylation analysis, the proposed approach opens up the study of DNA methylations dynamics and epigenetic regulation to the vast research space of uncultured bacteria. Finally, de novo detection of methylation motifs in a metagenomic community also holds promise for the discovery of novel MTases and restriction enzymes, expanding the repertoire of enzymes available for use in biomedical research.

This study focuses on one of the three forms of DNA methylations 6mA (N<NUM>-methyladenine) because it is the most abundant DNA methylation in prokaryotes and it has strong signal-to-noise ratio in SMRT polymerase kinetics. Other less prevalent types of DNA methylation in bacteria, such as N<NUM>-methylcytosine (4mC, medium-to-high signal) and <NUM>-methylcytosine (5mC, low-to-medium signal) are also within the scope of the present invention. As single-molecule long-read sequencing technologies continue to mature, generating larger yields and longer reads, the longer read lengths will provide more robust composition and methylation signatures that can be leveraged to more effectively segregate metagenomic reads, while also leading to even longer contigs with higher quality.

Though the present embodiments focus on SMRT sequencing, the binning framework of the invention applies generally to other third-generation technology, for example Oxford Nanopore. By integrating the features of second- and third-generation sequencing with complementary approaches, like Hi-C intrachromosomal maps, contig coverage covariance or single cell techniques, practitioners in the microbiome and metagenomics arts will gain a much more complete understanding of both the genomic and epigenomic landscape of complex microbial communities.

In some embodiments, the prokaryotic organisms in the microbiome sample do not have high sequence similiarity. In some embodiments, two or more of the prokaryotic organisms in the microbiome sample have high sequence similarity. In some embodiments, two or more of the prokaryotic organisms in the microbiome sample have an average nucleotide identity of greater than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, than about <NUM>%, or than about <NUM>%.

Any methylated nucleotides are within the scope of the methods of the present invention. In one embodiment the methylated nucleotides are selected from, without limitation, N<NUM>-methyladenine, N<NUM>-methylcytosine, and <NUM>-methylcytosine and combinations thereof.

Microbiome samples for use with the methods provided herein can be of any type that includes a microbial community comprising prokaryotic organisms. Prokaryotic organisms include, without limitation, bacterial organisms and archaeal organisms. The sample can include microorganisms from one or more domains. For example, in one embodiment, the sample comprises a heterogeneous population of bacteria and/or archaea.

Any prokaryotic organisms known to those skilled in the art are within the scope of the present disclosure. In one non-limiting embodiment, the prokaryotic organisms are bacterial organisms, archaeal organisms, and combinations thereof. In some non-limiting embodiments, the prokaryotic organisms are bacterial organisms, bacterial species, or strains of bacterial species. In other non-limiting embodiments, the prokaryotic organisms are archaeal organisms, archaeal species, or strains of archaeal species.

In one implementation, microbiome samples for use with the methods provided herein encompass, without limitation, samples obtained from the environment, including soil (e.g., rhizosphere), air, water (e.g., marine water, fresh water, rain water, wastewater sludge), sediment, oil, an extreme environmental sample (e.g., acid mine drainage, hydrothermal systems) and combinations thereof. In the case of marine or freshwater samples, the sample can be from the surface of the body of water, or any depth of the body of water, e.g., a deep sea sample. In one embodiment, the water sample is an ocean, a sea, a river, or a lake sample.

In one embodiment, the sample is a soil sample (e.g., bulk soil or rhizosphere sample). It has been estimated that <NUM> gram of soil contains tens of thousands of bacterial taxa, and up to <NUM> billion bacteria cells as well as about <NUM> million fungal hyphae (<NPL>). Bacteria, archaea, actinomycetes, fungi, algae, protozoa and viruses are all found in soil. Soil microorganism community diversity has been implicated in the structure and fertility of the soil microenvironment, nutrient acquisition by plants, plant diversity and growth, as well as the cycling of resources between above- and below-ground communities. Accordingly, assessing the microbial contents of a soil sample over time provides insight into microorganisms associated with an environmental metadata parameter such as nutrient acquisition and/or plant diversity.

The soil sample in one embodiment is a rhizosphere sample, i.e., the narrow region of soil that is directly influenced by root secretions and associated soil microorganisms. As plants secrete many compounds into the rhizosphere, analysis of the organism types in the rhizosphere may be useful in determining features of the plants which grow therein.

In another embodiment, the sample is a marine or fresh water sample. Ocean water contains up to one million microorganisms per milliliter and several thousand microbial types. These numbers may be an order of magnitude higher in coastal waters with their higher productivity and higher load of organic matter and nutrients. Marine microorganisms are crucial for the functioning of marine ecosystems; maintaining the balance between produced and fixed carbon dioxide; production of more than <NUM>% of the oxygen on Earth through marine phototrophic microorganisms such as Cyanobacteria, diatoms and pico- and nanophytoplankton; providing novel bioactive compounds and metabolic pathways; ensuring a sustainable supply of seafood products by occupying the critical bottom trophic level in marine foodwebs. Organisms found in the marine environment include viruses, bacteria, archaea and some eukarya. Marine bacteria are important as a food source for other small microorganisms as well as being producers of organic matter. Archaea found throughout the water column in the ocean are pelagic Archaea and their abundance rivals that of marine bacteria.

In another embodiment, the sample comprises a sample from an extreme environment, i.e., an environment that harbors conditions that are detrimental to most life on Earth. Organisms that thrive in extreme environments are called extremophiles. Though the domain Archaea contains well-known examples of extremophiles, the domain bacteria can also have representatives of these microorganisms. Extremophiles include: acidophiles which grow at pH levels of <NUM> or below; alkaliphiles which grow at pH levels of <NUM> or above; anaerobes such as Spinoloricus Cinzia which does not require oxygen for growth; cryptoendoliths which live in microscopic spaces within rocks, fissures, aquifers and faults filled with groundwater in the deep subsurface; halophiles which grow in about at least <NUM> concentration of salt; hyperthermophiles which thrive at high temperatures (about <NUM>-<NUM>. ) such as found in hydrothermal systems; hypoliths which live underneath rocks in cold deserts; lithoautotrophs such as Nitrosomonas europaea which derive energy from reduced mineral compounds like pyrites and are active in geochemical cycling; metallotolerant organisms which tolerate high levels of dissolved heavy metals such as copper, cadmium, arsenic and zinc; oligotrophs which grow in nutritionally limited environments; osmophiles which grow in environments with a high sugar concentration; piezophiles (or barophiles) which thrive at high pressures such as found deep in the ocean or underground; psychrophiles/cryophiles which survive, grow and/or reproduce at temperatures of about -<NUM>. or lower; radioresistant organisms which are resistant to high levels of ionizing radiation; thermophiles which thrive at temperatures between <NUM>-<NUM>. ; xerophiles which can grow in extremely dry conditions. Polyextremophiles are organisms that qualify as extremophiles under more than one category and include thermoacidophiles (prefer temperatures of <NUM>-<NUM>. and pH between <NUM> and <NUM>). The Crenarchaeota group of Archaea includes the thermoacidophiles.

In another implementation, microbiome samples for use with the methods provided herein encompass, without limitation, samples obtained from a subject, e.g., an animal subject, a protozoa subject, or a plant subject. The subject can be, for example, a human, mammal, primate, bovine, porcine, canine, feline, rodent (e.g., mouse or rat), or bird. In one embodiment, the animal subject is a mammal, e.g. a human. In one embodiment, the human subject is an adult, a child, an adolescent, an adult, or an elderly person.

In one embodiment the sample obtained from an animal subject is a body fluid. In another embodiment, the sample obtained from an animal subject is a tissue sample. Non-limiting samples obtained from an animal subject include tooth, perspiration, fingernail, skin, hair, feces, urine, semen, mucus, saliva, and gastrointestinal tract samples. The human microbiome comprises the collection of microorganisms found on the surface and deep layers of skin, in mammary glands, saliva, oral mucosa, conjunctiva and gastrointestinal tract. The microorganisms found in the microbiome include bacteria, fungi, protozoa, viruses and archaea. Different parts of the body exhibit varying diversity of microorganisms. The quantity and type of microorganisms may signal a healthy or diseased state for an individual. The number of bacteria taxa are in the thousands, and viruses may be as abundant. The bacterial composition for a given site on a body varies from person to person, not only in type, but also in abundance or quantity.

In the methods provided herein the one or more prokaryotic organisms can be of any type. For example, the one or more prokaryotic organisms can be from the domain Bacteria, Archaea, a combination thereof. Bacteria and Archaea are prokaryotic, having a very simple cell structure with no internal organelles. Bacteria can be classified into gram positive/no outer membrane, gram negative/outer membrane present and ungrouped phyla. Archaea constitute a domain or kingdom of single-celled microorganisms. Although visually similar to bacteria, archaea possess genes and several metabolic pathways that are more closely related to those of eukaryotes, notably the enzymes involved in transcription and translation. Other aspects of archaeal biochemistry are unique, such as the presence of ether lipids in their cell membranes. The Archaea are divided into four recognized phyla: Thaumarchaeota, Aigarchaeota, Crenarchaeota and Korarchaeota.

DNA methylation profiles inferred from SMRT sequencing provide an informative orthogonal epigenomic feature that can improve contig clustering. The DNA methylation profile is analogous to the sequence composition profile and the differential coverage profile, where normalized k-mer frequencies across k-mers and normalized coverage values across samples provide features for discriminative binning, respectively.

In the case of contig methylation profiles, each contig has a feature set consisting of contig-level DNA methylation scores across sequence motifs (see Examples).

The methylation score for a given motif on a contig reflects the extent to which all instances of that motif on the contig are methylated. It is calculated using inter-pulse duration (IPD) values, which records the time it takes a DNA polymerase to translocate from one nucleotide to the next during real-time DNA synthesis, often referred to as the polymerase kinetics. The methylation score for a motif on a contig becomes more reliable for predicting DNA methylation with an increase in two values: (<NUM>) the number of motif sites on the contig, which is generally larger for shorter motifs, and (<NUM>) the number of reads aligning to the contig, as each read contributes independent IPD measurements of methylation likelihood at the motif site. Evaluation based on methylation data from a bacterium with a set of well-characterized N<NUM>-methyladenine (6mA) motifs suggests that the specificity and sensitivity of methylation scores for detecting methylated motifs improve dramatically with an increase in the number of individual IPD values used to calculate them (<FIG>; see Examples).

A critical first step in using methylation profiles for binning is to identify the methylated motifs in the metagenomic assembly, as only those motifs that are methylated on one or more contig will contribute to the discriminative power of the binning. Therefore, a motif filtering method was designed to identify the relatively small number of motifs with scores suggesting likely methylation, excluding from the downstream analysis the vast majority of motifs that lack evidence of methylation on any contigs in the assembly (see Examples). In the Examples presented below, motif filtering simplifies the motif feature space from over <NUM>,<NUM> to between <NUM>-<NUM> motifs in metagenomic assemblies. The precise number of motifs that remain after filtering is often not critically important as long as the set of remaining motifs jointly captures the most significant differences between contig methylation profiles. This property contrasts with existing methods for methylation motif discovery that attempt to rigorously identify the single most parsimonious version of a motif. The proposed motif filtering is more robust to noise and different threshold choices, making it more effective and flexible for leveraging SMRT sequencing polymerase kinetics in a metagenomic setting.

To evaluate the ability of this procedure to segregate contigs based solely on DNA methylation profiles, a synthetic metagenomic mixture was created consisting of SMRT sequencing reads from eight separately sequenced bacterial species (Table <NUM>, below), four of which belong to the genus Bacteroides (see Examples).

The reads were combined and de novo assembly was done using the hierarchical genome-assembly process (HGAP3). The motif filtering procedure of the description de novo identified <NUM> motifs from the metagenomic contigs, <NUM> (<NUM>%) of which are exact matches to the true methylated motifs (as determined by separate methylation analysis for each species independent from the creation or analyses of the synthetic mixture; (Table <NUM>, below). The remaining two motifs are closely related to and provide similar methylation signals to the true motifs. Hierarchical clustering of the largest contigs from each species and their motif methylations scores shows that among the <NUM> motifs selected by motif filtering, each species in the mixture has a unique methylation profile (<FIG>).

To ease visualization and interpretation of high-dimensional features of many metagenomic contigs, dimensionality reduction was used to reduce the feature space to two dimensions that are amenable to plotting. The dimensionality reduction algorithm primarily used in this study is the Barnes-Hut approximation of t-distributed stochastic neighbor embedding (t-SNE) (see Examples), which has already been demonstrated to be effective at segregating metagenomic contigs based on k-mer frequency. Because t-SNE is a non-linear dimensionality reduction algorithm that is designed to preserve local pairwise distances, it differs from linear methods, such as principal components analysis (PCA) that captures global variance, making t-SNE well suited for complex microbiome communities with subpopulation structures that are not effectively captured by PCA.

The 2D map generated by applying t-SNE to the matrix of methylation profiles (<NUM> motifs for each contig) reveals contigs that are generally well separated based on their known species (<FIG>). Specifically, the four species from the Bacteroides genus show remarkably clear separation from each other, despite the fact that the genomes share significant sequence similarity (Table <NUM>, below). This separation of the four Bacteroides species is clearer than is possible using composition methods alone (<FIG>) and cleaner than when the contig coverage values are included with composition (<FIG>). The methylation-based map results in a cluster silhouette coefficient, which ranges between -<NUM> (significant mixing) and <NUM> (complete separation), of <NUM>, while the composition-based clustering results in a <NUM> silhouette coefficient.

Interestingly, there is some mixing of small contigs that are likely too short to contain IPD values from the full set of methylated motifs for a species. This is supported by the observation that several contigs belonging to Clostridium bolteae, which are too small to contain the full diversity of C. bolteae methylated motifs (<FIG>), cluster more closely with Ruminococcus gnavus, a species without any detectable methylation motifs. While some organisms will be, like R. gnavus, absent any detectable methylation, these are relatively rare.

Having demonstrated how methylation profiles can be used for contig binning in a mock metagenomic community, next the approach was applied to examine a microbial community sampled from an adult mouse gut. <NUM> rRNA sequencing (see Examples) indicated that the sample was complex and dominated by an undefined number of organisms from the S24-<NUM> family of the order Bacteroidales (<FIG>). SMRT sequencing reads were assembled using the HGAP3 assembler (Table <NUM>).

<NUM> methylated motifs were detected from the assembled contigs and visualized the methylation landscape of the sample by using t-SNE to reduce the <NUM> dimensions to a 2D scatter plot (<FIG>). The resulting scatter plot reveals nine distinct bins of contigs with consistent methylation profiles. In eight of the nine bins, the uniform contig coverage values within each bin support that the contigs correspond to eight single organisms, while the split coverage values in bin7 suggest that it may contain contigs from two different genomes (<FIG>).

Next, CheckM was used to assess the genome completeness and contamination of each bin based on single-copy gene counts. Eight of the nine bins have ><NUM>% completeness and only bin7 has significant contamination, likely from the second genome in the bin (Table <NUM>, below).

Querying the contig sequences in each bin against a manually curated set of <NUM> publicly available mouse gut microbial references revealed significant reference hits with eight of the nine bins (<FIG>; Table <NUM>, below), providing further support that the bins identified using methylation profiles represent the genomes of distinct organisms.

Bin4 and bin5 have high-quality, nearly full-length matches with the finished genomes for Akkermansia mucinophilia YL-<NUM> (average nucleotide identity (ANI) = <NUM>%) and Parabacteroides sp. YL-<NUM> (ANI = <NUM>%), respectively. The remaining six bins have high-quality matches with genome assemblies of species that have been identified in the mouse gut in other studies but lack finished reference sequences. Three of these six bins have full-length matches with three draft assemblies of uncultured members of the Bacteroidales S24-<NUM> family: bin1 matches Bacteroidales bacterium M1 (ANI = <NUM>%), bin3 matches Bacteroidales bacterium M12 (ANI = <NUM>%), and bin8 matches Bacteroidales bacterium M2 (ANI = <NUM>%). The final three bins have high-quality matches with three unidentified metagenomic species (MGS) previously binned in a large study of mouse gut microbiomes: bin2 matches MGS:<NUM> (ANI = <NUM>%), bin8 matches MGS:<NUM> (ANI = <NUM>%), and bin9 matches MGS:<NUM> (ANI = <NUM>%). The seven Bacteroidales bins all share high ANI with each other (<NUM>-<NUM>% ANI), but at values suggesting inter- rather than intraspecies relationships (Table <NUM>).

Because the only other family of Bacteroidales identified in the sample by <NUM> sequencing was the family Rikenellaceae at <NUM>% abundance, it is likely that these seven highly contiguous genome bins all belong to the poorly characterized S24-<NUM> family of Bacteroidales that dominated the <NUM> abundance profile for the sample (<FIG>). Quality alignment of the bin5 contigs to the reference for Parabacteroides sp. YL-<NUM> was observed, which is classified as belonging to the closely related Bacteroidales family Tannerellaceae, but there is some apparent divergence in the alignment that raises doubts about it being an exact match (<FIG>). Collectively, these comprehensive evaluations demonstrate that the nine bins isolated using methylation profiles represent highly contiguous draft assemblies for organisms that were previously uncharacterized or only represented by fragmented WGS assemblies.

Next, the mouse gut microbiome community was explored by leveraging the complementarity of methylation-based binning with existing methods that utilize differential coverage and sequence composition, such as CONCOCT, GroopM, and MetaBAT, which have been demonstrated to be powerful methods for isolating genomes in complex metagenomic samples. Illumina WGS data from <NUM> publically available mouse gut samples was aligned to the assembled contigs in order to generate coverage values for each sample. CONCOCT was then applied, which combines contig <NUM>-mer frequency profiles with the coverage profiles to call genome bins. This analysis generated high-quality bins of near-complete genomes for several organisms, including members of the order Clostridiales (mapped to MGS:<NUM>), Verrucomicrobiales (mapped to A. mucinophilia YL-<NUM>), and two organisms that do not have methylation bins, Burkholderiales and Lactobacillales (<FIG>; Table <NUM>, below). However, CONCOCT assigned multiple Bacteroidales genomes to a single bin containing <NUM> Mbp of sequence. A further analysis showed that the co-binning of several Bacteroidales genomes by CONCOCT is due to the high similarity in their abundance profiles across microbiome samples, even after excluding genomic regions where sequence similarity might cause reads to map to multiple Bacteroidales genomes (<FIG> and Examples). Therefore, although differential coverage binning proved very effective for binning many organisms in the sample, it did not effectively handle organisms with similar coverage covariance profiles.

Collectively, the above analyses highlight the great discriminative power of methylation-based binning and its complementarity with existing methods for improving binning resolution in complex microbiome samples. In recognition of this, the present analysis pipeline was extended to assess methylation profiles at the level of reads, contigs and bins, where the binning assignments can come from various differential coverage binning software. This approach allowed to discover eight additional motifs at the bin level that were not detectable by focusing on individual contigs (Table <NUM>, above).

An analysis of an infant gut microbiome was also performed to illustrate additional ways in which methylation profiles can be integrated with sequence composition features (see Example <NUM>).

Bacterial communities often contain a significant extra-chromosomal genetic potential in the form of mobile genetic elements (MGEs). MGEs may include, without limitation, plasmids, transposons (including class I or retrotransposons, class II or DNA transposons, and insertion sequences), bacteriophages (including bacteriophage elements such as Mu), and introns (including group I introns and group II introns).

Transposons (transposable elements, or TEs) are DNA sequences that can change their position within a genome, sometimes creating or reversing mutations and altering the cell's genetic identity and genome size. It has been shown that transposons are important in genome function and evolution. Transposons are also useful to researchers as a means to alter DNA inside a living organism. There are at least two classes of TEs: Class I TEs or retrotransposons generally function via reverse transcription, while Class II TEs or DNA transposons encode the protein transposase, which they require for insertion and excision, and some of these TEs also encode other proteins.

Bacteriophages (phages) are viruses that infect and replicate within a bacterium. Bacteriophages are composed of proteins that encapsulate a DNA or RNA genome, and may have relatively simple or elaborate structures. Their genomes may encode as few as four genes, and as many as hundreds of genes. Phages replicate within the bacterium following the injection of their genome into its cytoplasm. Bacteriophages are ubiquitous viruses, found wherever bacteria exist. It's estimated there are more than <NUM><NUM> bacteriophages on the planet.

Plasmids are small (typically <NUM>-<NUM> kbp), circular, and highly mobile DNA elements can be transferred among host bacteria during conjugation events or through natural transformation of extracellular plasmids into competent cells, making them an important mediator of HGT in bacteria. The genes encoded by plasmids can confer antibiotic resistance, encode virulence factors or provide specific metabolic functions that allow the host cell to survive under conditions that would otherwise be hostile. If a plasmid has a broad range of acceptable host species, the genes encoded by that plasmid, for instance those conferring antibiotic resistance, can be added to the genetic repertoire of a large number of species. It is therefore critically important to determine the host species of plasmids in microbiomes, as this information not only reflect the full genetic catalog of the host, but can also be used to track the transmission of antibiotic resistance elements across different members of a bacterial community.

MGE replication can be independent of chromosomal replication, meaning that the sequenced coverage of, e.g., a plasmid will likely differ significantly from the sequenced coverage of the chromosomal contigs of its host. Furthermore, empirical evidence supports the hypothesis that sequence composition alone is often not capable of mapping a plasmid to its host in a metagenomic setting. By examining the WGS sequencing data from <NUM>,<NUM> plasmids and the chromosomes of their host species in the REBASE database, it was observed that the plasmid sequence composition profile (i.e. the vector of <NUM>-mer frequencies) can differ significantly from that of the host chromosome (<FIG>). While the majority of Euclidian distances (d) between plasmid composition profiles and those of their host chromosome fell between <NUM> and <NUM>, many greatly exceeded that distance and fell under the empirical distribution created by calculating sequence distance between plasmids and randomly sampled chromosomes (see Examples). Highly dissimilar composition between host and plasmid might suggest a recent HGT event where the plasmid was acquired from a distant donor species. However, even in cases of moderate dissimilarity between host chromosome and plasmid sequence composition, there is not a clear strategy for determining which organism might be host to a particular plasmid.

Due to the difficulty of resolving complex repeats and mobile genetic elements, assembling complete plasmid sequences using short-read technologies has proved challenging. While SMRT sequencing is capable of generating high-quality, closed plasmid assemblies from clinical isolates, little work has been done to generate whole plasmid sequences from a metagenomic sample and associate the plasmids to their host bacterial species in the community. To do this, the present invention takes advantage of the fact that plasmid DNA and the chromosomal DNA of the bacterial host are both methylated by the same set of MTases. The result is that the methylation profiles of the plasmids match the methylation profile of its host bacterium. This phenomenon is demonstrated by transforming the <NUM> kbp plasmid pHel3 from Escherichia coli DH5α into E. coli CFT073 and Helicobacter pylori JP26. In each case, SMRT sequencing was used (Table <NUM>, below) to show that the methylation profile of pHel3 inherits that of its new host strain (<FIG>).

In order to evaluate the general potential of using methylation profiles for mapping plasmids in a community, the wealth of publicly available SMRT sequenced bacteria in the REBASE database was next surveyed, which consists of the assembled sequences and the observed methylated motifs for <NUM> genomes and <NUM> plasmids. Because successful mapping of a plasmid to its host requires a sufficient diversity of methylated motifs within a specific community, communities of different sizes were simulated by randomly selecting entries in the REBASE database and assessed the methylome diversity in each mock community. As the number of organisms in a community increases, the number of organisms with unique methylomes, expressed as a fraction of the community size, decreases but still remains fairly high even in communities consisting of <NUM> species (<FIG>). As expected, the decrease is more pronounced when multiple strains of a species are added to a community. Similar values of methylome uniqueness are observed when surveying only organisms in REBASE that have at least one known plasmid (<FIG>).

Plasmid size is another consideration for methylation-based host mapping, as shorter plasmids are less likely to possess instances of the full suite of methylated motifs that can help conclusively demonstrate a matching methylation profile with that of a host genome. Sequences of different lengths were simulated from the REBASE genomes and assessed how frequently these sequences contained the full set of the methylated motifs from the source genome (<FIG>). It was found that, on average, <NUM>% of <NUM> kbp sequences will contain instances of at least three quarters of the 6mA motifs, while <NUM>% of <NUM> kbp sequences will capture instances of all the 6mA motifs. Therefore, the rich methylation profiles required for mapping to a host genome are more likely to occur with larger, rather than smaller, plasmids. However, a partially complete methylation profile (i.e. lacking one or more methylated motifs) might be sufficient to unambiguously map the plasmid to its host if the methylated motifs included in the plasmid sequence are uniquely methylated by the host bacteria in a specific microbiome sample. Additional analysis of methylation motifs in an outbreak strain of Klebsiella pneumoniae underscores how methylation profiling could help identify the host of a 362kb plasmid carrying thirteen antibiotic resistance genes in a metagenomic sample (see Example <NUM>).

Building on the important considerations learned from the above analysis, the methylation-based plasmid-host mapping procedure was first applied using the mock community of eight bacterial species, where the true mappings are known. Six closed circular sequences are identified from the SMRT contigs assembled from the mock community by HGAP3 (see Examples). A confident mapping of a plasmid to a host is defined if contigs accounting for ><NUM>% of the host genome contain (<NUM>) the same methylated motifs (i.e. motifs with methylation score ≥ <NUM> calculated from ≥ <NUM> IPD values) that are found on the plasmid, and (<NUM>) no additional methylated motifs. Using this approach, the correct host was recovered using methylation profiles in four of the six circular contigs (<NUM>%), including the only known plasmid in the group, B. thetaiotaomicron plasmid p5482 (GenBank accession AY171301. The remaining two circular contigs were not mapped to the wrong host, but were just too short (<<NUM> kbp) to contain sufficient motif sites for a conclusively mapping, consistent with the estimations from the above simulation analysis (<FIG>).

Next, the methylation-based plasmid-host mapping procedure was applied to the adult mouse gut microbiome sample. <NUM> contigs between <NUM>-<NUM> kbp were identified, of which eleven are fully circularized and nine are conjugative transposon elements (encoding at least five genes annotated as conjugative transposon-related). Thirteen of these mobile genetic elements (MGE) did not assemble using the original complex metagenomic reads, but were only discovered by isolating the reads that map to contigs in each methylation bin and reassembling them in a single genome setting (see Examples). Using the same methylation-based criteria defined above, eight of the <NUM> discovered MGEs were confidently mapped to distinct methylation bins containing genomes from the order Bacteroidales (Table <NUM>, above). These eight mapped MGEs include five highly likely plasmids (<50kb circular contigs containing origins of replication) and three conjugative transposons. Conjugative transposons are known to play an important role in HGT and the spread of antibiotics in Bacteroidales, and they have been implicated in sequence sharing between multiple Bacteroidales species in the human gut. Collectively, these analyses demonstrate that DNA methylation can be exploited as a novel discriminative feature for MGE-host (e.g., plasmid-host) mapping in complex microbiome samples.

Highly variable organism abundances in metagenomic samples often present significant challenges to de novo assembly tools, especially for the low abundance species. Because it can be expected that some community members will not be represented among the assembled contigs, a more complete representation of the community might be achieved by binning unassembled metagenomic sequencing reads alongside the assembled contigs. Multiple tools use unsupervised binning of metagenomic short reads, but the insufficient sequence information content in short reads limits their accuracy and practical applicability outside of very low-complexity metagenomic samples. While third generation sequencing platforms produce amplification-free reads with much longer read lengths, the raw reads are confounded by a high single-pass error rate (typically ~<NUM>% for SMRT sequencing). Although it has been shown that longer contig sequences result in greater segregation using <NUM>-mer frequency vectors and t-SNE, it remained a fundamental question whether this would also apply to high-error unaligned SMRT reads.

To evaluate the ability of <NUM>-mer frequency metrics to bin unassembled SMRT reads and assembled contigs together, a synthetic microbiome (mixed DNA from the <NUM>-member Mock Community B) created as part of the Human Microbiome Project (HMP) was first analyzed. The original mock community contained each member in roughly equal proportion, making it an unrealistic mixture. The reads were therefore downsampled (see Examples) to create a distribution of relative abundances that follows a log curve, where the most predominant species, Streptococcus mutans (294x coverage), is present at <NUM> times the abundance of the most minor species, Rhodobacter sphaeroides (2x coverage) (<FIG>; Table <NUM>).

<NUM>-mer frequency metrics for all HMP mock community sequences (unassembled SMRT reads and assembled contigs) were subjected to t-SNE. In the resulting 2D map, only the contigs were first visualized and annotated using Kraken, revealing a clean separation of contigs from species for which there is a significant number of assembled bases (<FIG>). To ensure that the contig separation in the 2D map was not biased due to poor assembly results of the low-abundance members of the downsampled community, these findings were confirmed using the even abundance community, finding consistent results (<FIG>). Next, binning quality of the unassembled SMRT long reads was assessed. Notably, the <NUM>-mer frequency profiles are resilient to the random error in the long reads; the clusters of unassembled reads are highly species-specific. sphaeroides, while poorly represented in the set of assembled contigs (<FIG>), is clearly present as a distinct cluster of unassembled reads (<FIG>), highlighting the benefit of including unassembled reads in the composition-based binning to reveal the presence of very low-abundance species that are not captured by metagenomic assembly. The additional sequence information in longer reads provides more stable <NUM>-mer frequency profiles and tighter clusters compared to clustering of shorter reads (<FIG>). Furthermore, a 2D histogram provides an overview of global community complexity even absent any sequence annotation (<FIG>), making it possible to identify a set of novel sequences from a particular taxon and investigate them further. This analysis highlights the feasibility of direct binning of single molecule long reads even though the raw error rate is high, and the promise of joint binning of unassembled long reads and assembled contigs for more complete representation of a microbiome with low abundance species.

Next, single molecule long reads from third generation sequencing were also binned using their read-level methylation profiles. This can help avoid or resolve chimeric contigs, which occur when multiple strains in a mixture are assembled into contigs built from reads originating from different strains. The significant challenges associated with chimeric contigs affect coverage- and k-mer-based binning methods, hinder strain-specific variant calling and, in the case of single-molecule long-read sequencing, confound the identification of strain-specific methylation on each contig. Importantly, because MTases often transmit across species and strains by HGT, closely related strains with high sequence similarity often encode different MTases that target unique combinations of methylation motifs and provide a novel opportunity to de-convolve co-existing strains in a microbiome sample. A measure of read-level methylation that was developed for the study of epigenetic heterogeneity in single organisms was built on and extended to assess read-level epigenetic heterogeneity in a metagenomic setting (see Examples).

To demonstrate how this can improve multi-strain assemblies, two synthetic mixtures of reads were constructed (see Examples) from (<NUM>) two strains of H. pylori (Table <NUM>) and (<NUM>) three strains of E. coli (Table <NUM>).

Despite the high sequence similarity of the strains in each mixture (Tables <NUM> and <NUM>), they encode different MTases that result in different sets of methylated motifs.

The first mixture contained reads from the H. pylori strains J99 and <NUM> that assembled together into one small contig from strain <NUM> and another large, highly chimeric contig (<FIG>). To reduce the chimerism in the assembly, a pre-assembly binning strategy was adopted analogous to that described by Cleary et al. , but instead of using k-mer co-abundance for binning reads, single-molecule long reads were separated into bins based on their methylation profiles and subsequently assembled each bin. A small set of four high-density motifs (GATC, GAGG, TGCA, CATG) are sufficient to differentiate these two H. pylori strains (Table <NUM>) and were selected to generate methylation profiles for individual single molecule reads.

Principal component analysis (PCA) was then used for the dimensionality reduction step to generate a 2D plot of each mixture, revealing a bimodal concentration of reads organized solely by their methylation profiles (<FIG>). Dimensionality reduction using t-SNE also revealed two strain-specific clusters, but the resulting clusters did not follow a Gaussian distribution, making delineation of them less straightforward than with PCA (<FIG>). Because the number of features is small (i.e. four motifs), PCA provides cleaner separation of Gaussian subpopulations in this application than does t-SNE; this also suggests that different dimensionality reduction methods may complement each other in different applications. Finally, the epigenetically binned reads were assembled separately using HGAP3 with the same parameters used for the mixed assembly, resulting in separately assembled contigs with improved contiguity, including chromosome-scale contigs for both strains, and minimal chimerism (<FIG>).

The read-level methylation binning procedure was next applied to another data set consisting of SMRT reads from three strains of E. coli from distinct serotypes: O26:H11, O103:H11, and O111 (see Examples). An assembly of these mixed reads results in many highly chimeric contigs and very few contigs that are specific to a strain (<FIG>). The motifs that differentiate these strains, AGCACY, CRARCAG, GGTNACC, and CTGCAG, are longer (Table <NUM>), which are more likely to be disrupted by the random nature of the sequencing errors in unaligned single molecule long reads, causing incorrect IPD values for each long motif.

Addressing this required an additional alignment step to error correct the reads prior to calculating the scores for the methylation profiles. Specifically, the reads from each strain were aligned to the standard E. coli K12 MG1655 reference sequence (RefSeq accession NC_000913. <NUM>) then calculated read-level methylation scores for each motif. Methylation profiles were again visualized using PCA and reads were binned based on visible subpopulations (<FIG>). Finally, isolated assembly of reads from each bin resulted in a substantial reduction of contig chimerism and an increase in contigs containing sequence specific to each E. coli strain (<FIG>).

Recent advances in library preparation protocols for Illumina sequencing have made it possible to generate synthetic long reads of several kilobases in length. The read lengths of synthetic long reads can approach those generated by SMRT sequencing, yet important differences between the technologies have implications for their specific applications in metagenomics and therefore warrant a detailed investigation. Because the capability to infer methylation events is a unique strength of SMRT sequencing as studied above, other aspects of the two techniques and their potential complementarity are emphasized here.

The read lengths and high accuracy of synthetic reads have enabled researchers to phase substrain-level bacterial haplotypes in metagenomic samples. By aligning synthetic long reads to contigs generated through de novo metagenomic assembly, the study revealed the presence of multiple genotypes within the same strain. A prerequisite for substrain haplotyping with synthetic long reads is a metagenomic assembly that serves as a reference for the read alignment. Kuleshov et al. acknowledge that SMRT reads are more likely to result in large draft assemblies, and indeed point out that contigs assembled from SMRT reads are significantly larger than those assembled using synthetic long reads, even when the latter was supplemented by traditional short reads.

Given the multi-kb read lengths and high accuracy of synthetic long reads, it was sought to understand why they resulted in more fragmented and less comprehensive assemblies than did SMRT reads. To this end, both the synthetic long reads sequenced from the <NUM>-member HMP Mock Community B (staggered abundance; HM-277D) and the SMRT reads from the same community were aligned to their reference genomes. Because the SMRT reads were sequenced from a different version of the HMP Mock Community B (even abundance; HM-276D), the aligned reads were downsampled so that total numbers of aligned bases for each organism were roughly equal for both sequencing technologies (see Examples; Table <NUM>, above).

Despite considering approximately the same number of aligned bases for each technology, SMRT reads covered a higher percentage of genome positions in <NUM> of the <NUM> species and matched the percentage of genome positions covered by synthetic long reads in the remaining three species (<FIG>; Table <NUM>).

In several cases, the increases in genome coverage over synthetic long reads were dramatic: SMRT sequencing of D. radiodurans, A. odontolyticus, E. faecalis, and S. pneumoniae covered an additional <NUM>%, <NUM>%, <NUM>%, and <NUM>% of their genomes, respectively. The genomes with the highest GC-content (R. sphaeroides, <NUM>% GC; D. radiodurans, <NUM>% GC; P. aeruginosa, <NUM>% GC; A. odontolyticus, <NUM>% GC) were among those that saw significant increases in genome coverage with SMRT reads compared to synthetic long reads (Table <NUM>). This observation is consistent with previous studies showing that the PCR amplification of DNA fragments required for synthetic long read sequencing is sensitive to genomic GC-content and can result in significant coverage biases (i.e. highly nonuniform sequence coverage).

SMRT sequencing, however, is an amplification-free protocol and is not subject to GC bias, resulting in more uniform coverage profiles across genomes (<FIG>). Further illustrating this phenomenon are three small regions from the genomes of S. agalactiae, S. aureus, and P. aeruginosa (<FIG>), which are representative of many of the genomes in the mock community (<FIG>). The synthetic long reads coverage profiles consist of peaks and valleys, representing over- and under-amplified DNA fragments, respectively. Some of the valleys result in complete coverage dropouts, across which genome assembly becomes impossible. The SMRT sequencing protocol, on the other hand, results in much more uniform coverage profiles and fewer coverage dropouts, making it more amenable to metagenomic assembly and more likely to result in chromosome-scale contigs.

Two additional sources of systematic error in the synthetic long reads, resulting from dilution and sub-assembly steps in the protocol, make it more difficult to assemble high abundance species and regions containing tandem repeats. These steps are unique to synthetic long reads and do not apply to SMRT sequencing, which might further contribute to the superiority of SMRT reads for generating large metagenomic assemblies. The strengths of synthetic long reads, however, lie in their ability to call (and phase) local genomic features, such as single nucleotide variants (SNVs) or short insertions and deletions. Overall, this suggests a complementary strategy for maximizing assembly quality with SMRT sequencing and leveraging synthetic long reads for variant calling and haplotyping.

Methylation binning of contigs alone may, in some instances, to be challenging for organisms that are present at low-abundance in high-complexity samples, as it is difficult to detect methylated motifs from the small contigs that are typically assembled from low-abundance organisms. However, this can be complemented by binning assignments from cross-coverage and composition-based binning tools, such as CONCOCT, because contigs can be phased together according to third-party binning assignments to aid the discovery of methylated motifs, as was demonstrated with the mouse gut microbiome analysis. De novo methylation motif detection is well powered at the levels of contigs or bins, but is challenging at the level of single reads due to the requirement for long read length, especially for large, sparsely distributed motifs. However, read-level binning by methylation profiles can build on a priori knowledge of the methylation motifs in a species of interest for the de-convolution of multiple co-existing strains, as illustrated in this study. Continued increases in read length of third-generation sequencing also raise the prospect of more reliable de novo detection of methylated motifs at the single read-level in the near future.

The choice of SMRT sequencing libraries of long insert size can improve contiguity in a metagenomic assembly, but the size selection procedure may filter out short MGEs like plasmids and phages. The choice of library size would depend on goals specific to the particular research study. When resource allows, combinations of long and short libraries can be integrated to achieve both good assembly contiguity and the good coverage of short MGEs, although challenges currently exist in assembling complex MGEs from shorter reads. Integrating additional sequence data from a rolling circle amplification library might help to highlight plasmids that are excluded from the standard SMRT library or do not fully circularize in the SMRT assembly.

Although the long reads and methylation profiles made possible by SMRT sequencing (and other third-generation sequencing technologies) hold great promise for studying microbial communities, they currently require more input DNA than second generation sequencing technologies. However, this requirement has decreased recently as the SMRT technology has matured and further reductions are anticipated in the future, given the active development and pace of technological improvement.

The following examples illustrate specific aspects of the instant description. The examples should not be construed as limiting, as the examples merely provide specific understanding and practice of the embodiments and their various aspects.

Using metagenomic sequencing data from several synthetic and real microbiome samples, comprehensive evaluations of the proposed approach were performed and it was demonstrated that DNA methylation is a novel and rich feature that provides significant discriminative power capable of complementing existing methods for high-resolution metagenomic binning.

Code availability. The software supporting all proposed methods is implemented in Python and is available with full documentation at the world wide web github. com/fanglab/mbin.

Epigenetic information was used to segregate contigs assembled from highly similar strains that would be otherwise indistinguishable using k-mer frequency-based methods. Two sets of infant gut microbiota obtained from stool samples of children who were selected for sequencing based on a high genetic risk for development of T1D were examined.

Interestingly, it has been observed that the particular species of Bacteroides that dominates the composition of both samples, Bacteroides dorei, often spikes in relative abundance prior to onset of T1D in children, making it an important species to understand and potentially monitor during early adolescence. <NUM> sequencing showed that the two samples contained two distinct strains of B. dorei: Sample A consisted of <NUM>% B. <NUM> (CP007619), while Sample B contained <NUM>% B. <NUM> (CP008741). Despite a high sequence similarity between the two B. dorei strains (Table <NUM>), each strain has a unique set of methylated sequence motifs and therefore a unique methylation profile.

SMRT sequencing data were collected for the two microbiome samples from a previous study (Table <NUM>) and performed a metagenomic de novo assembly using a combination of both gut samples to generate a mixture of contigs from both B. dorei strains in the output set of metagenomic contigs. Lacking any labeling for these contigs, the sequence annotation tool Kraken was applied for labeling of all non-B. dorei contigs and an alignment-based labeling approach for distinguishing the two B. dorei strains (See Examples).

Composition-based binning was first conducted using <NUM>-mer frequency profiles , followed by t-SNE dimensionality reduction (<FIG>). The map has five distinct clusters of contigs, four of which consist primarily of a combination of contigs from multiple species or strains. This suggests that composition-based binning is insufficient to segregate the two strains of B. dorei due to their high sequence similarity. Notably, composition-based binning also fails to segregate Bacteroides fragilis from Bacteroides thetaiotaomicron, Bifidobacterium breve from Bifidobacterium longum, and Alistipes finegoldii from Alistipes shahii.

Motif filtering identified seven motifs with significant methylation scores on at least one contig in the assembly: GGATCA, GATCA, TTCGAA, GATC, CTCAT, GAATC, and GGATC. The resulting t-SNE map constructed using methylation profiles alone (<FIG>) resolves the contigs into four clusters. In contrast to the k-mer frequency-based map and as a consequence of their unique methylation profiles, the two strains of B. dorei are very well segregated in the methylation-based binning analysis. However, methylation-based binning alone did not fully segregate all other species due to an insufficient diversity of methylated motifs among them. This suggests that both methylation-based and composition-based binning methods can complement each other to compensate for the shortcomings of each approach. By combining k-mer frequency and methylation profiles, both reduced separately by t-SNE to 2D, into a single matrix with four columns, t-SNE was again used to reduce the matrix and generate a 2D scatter plot (<FIG>). This approach succeeds in separating the two strains of B. dorei from each other, B. fragilis from B. thetaiotaomicron, and B. breve from B. Only the two species from the Alistipes genus remain convoluted in the combined map, due to high sequence similarity and likely identical methylomes. Again using a silhouette coefficient to assess the contig clustering, it was found that while composition-based binning alone results in a silhouette coefficient of <NUM>, the integration with methylation-based binning increases the coefficient to <NUM>, demonstrating that contig methylation profile can help deconvolute contigs with high sequence similarity.

To assess the methylome diversity across strains of a clinically relevant bacterial species, the <NUM> bacterial strains in the REBASE database for which methylated motifs have been identified through SMRT sequencing were analyzed. Among these was a virulent and antibiotic-resistant strain of Klebsiella pneumoniae (strain <NUM>-<NUM>) isolated from a patient during a <NUM> outbreak in Germany. A single 362kb plasmid (pKpn23412-<NUM>) hosted by this strain contained thirteen antibiotic-resistance genes, including the blaCTX-M-<NUM> (Kpn23412_5431) gene responsible for conferring the extended spectrum β-lactamase (ESBL) phenotype of the bacteria. The plasmid also contained multiple replicons, which helps to expand the range of organisms in which the plasmid can successfully replicate.

The sequence composition profiles of this plasmid and the K. pneumoniae chromosome differed to an extent (Euclidian distance, d = <NUM>) that would prohibit any sequence-based mapping of plasmid to host in a metagenomic sample. However, the methylated motifs, including GATC and CCAYNNNNNTCC (SEQ ID NO: <NUM>), present an opportunity for linking the plasmid and host epigenetically. To demonstrate this, the methylated motifs of nine other species contained in the REBASE database were examined, all of which had chromosome sequence composition profiles closer to K. pneumoniae plasmid pKpn23412-<NUM> (d < <NUM>) than did the true host chromosome. Although some of the composition profiles are relatively similar to the plasmid, the methylation profiles are diverse, making it possible to match plasmid pKpn23412-<NUM> to its K. pneumoniae host (<FIG>). Finally, all <NUM> strains of K. pneumoniae contained in the Rebase database were examined, and it was found that the sequence of plasmid pKpn23412-<NUM> was roughly the same Euclidian distance from the chromosomes of each strain (<FIG>). However, these <NUM> strains include <NUM> distinct methylation profiles (i.e. different combinations of methylation motifs), one of which is found only in strain <NUM>-<NUM>. This means that if multiple K. pneumoniae strains were present in the same metagenomic sample, DNA methylation profiles may be able to help map plasmid pKpn23412-<NUM> to its true host strain directly from metagenomic data. This epigenetic plasmid-host mapping approach highlights the broad range of applications in which epigenetic profiles can be exploited to address difficult challenges in a variety of clinically relevant situations.

Bacteroides caccae ATCC <NUM>, Bacteroides ovatus ATCC <NUM>, Bacteroides thetaiotaomicron VPI-<NUM>, Bacteroides vulgatus ATCC <NUM>, Collinsella aerofaciens ATCC <NUM>, Clostridium bolteae ATCC BAA-<NUM>, and Ruminococcus gnavus ATCC <NUM> were grown individually in <NUM> of supplemented Brain-heart infusion broth in an anaerobic chamber from Coy Laboratory Products. Escherichia coli MG1655 was grown aerobically in <NUM> of LB broth. Construction of the 10kb DNA libraries for SMRT sequencing was performed according to the manufacturer's instructions.

A male <NUM>-week-old NOD/shiltj mouse (no. <NUM>, Jackson Labs) was housed in a Specific Pathogen Free (SPF) room at New York University Langone Medical Center (NYUMC). At the week <NUM> of life, the mouse was placed into a clean plastic container in a fume hood, and its fresh fecal pellets were collected in sterilized microcentrifuge tubes and frozen at -<NUM>. Fecal DNA was extracted using PowerSoil DNA isolation kit (MoBio Labs, Carsbad, CA). 10kb library preparation for SMRT sequencing was performed according to the manufacturer's instructions. The bacterial <NUM> rRNA gene V4 regions were amplified and libraries constructed as previously described by Livanos et al.

pylori shuttle plasmid pHel3 was electroporated from E. coli strain DH5α to strain CFT073 using MicroPulser following procedures recommended by the manufacturer (Bio-Rad Lab. , Hercules, CA). The same plasmid was also introduced from E. coli strain DH5α into H. pylori strain JP26 by natural transformation as previously described. coli DH5α carrying pHel3 and CFT073 carrying pHel3 were grown in Luria-Bertani (LB) medium with kanamycin (Km; <NUM>µg/ml) at <NUM> for <NUM> hours. pylori JP26 carrying pHel3 were grown in Brucella broth (BB) medium supplemented with <NUM>% newborn calf serum (NBCS) and Km (<NUM>µg/ml) at <NUM> in microaerophilic condition for <NUM> hours. Bacterial cell pellets of E. pylori cultures were collected by centrifugation, genomic DNA of each culture was purified using Wizard Genomic DNA Purification Kit (Promega, Madison, WI), and plasmid DNA of each culture was purified using QIAprep Spin Miniprep Kit (QIAgen, Valencia, CA). 2kb library preparation for SMRT sequencing genomic and plasmid DNA for each culture was performed according to the manufacturer's instructions.

Genomic DNA for the three strains of E. coli, BAA-<NUM>, BAA-<NUM>, and BAA-<NUM>, were purchased from ATCC and construction of the 10kb DNA libraries for SMRT sequencing was performed according to the manufacturer's instructions.

Primer was annealed to the size-selected SMRTbell with the full-length libraries (<NUM> for <NUM> minute <NUM> seconds followed by decreasing the temperature by <NUM> to <NUM>). The polymerase-template complex was then bound to the P6 enzyme using a ratio of <NUM>:<NUM> polymerase to SMRTbell at <NUM> for <NUM> hours at <NUM> and then held at <NUM> until ready for magnetic bead loading, prior to sequencing. The magnetic bead-loading step was conducted at <NUM> for <NUM>-minutes per manufacturer's guidelines. The magnetic bead-loaded, polymerase-bound, SMRTbell libraries were placed onto the RSII machine at a sequencing concentration of <NUM> - <NUM> pM and configured for a <NUM>-minute continuous sequencing run.

Sequencing of the <NUM> V4 region was performed using the Illumina MiSeq platform as previously described by Livanos et al.

All k-mer frequency metrics in this study used a k-mer size of <NUM>. Counts of pairs of pentamers that are reverse complements of each other were combined, resulting in a set of <NUM><NUM>-mers as composition features for each sequence (contig or single-molecule read). Following the procedure described by Alneberg et al. , a small pseudo-count was added to each <NUM>-mer count to ensure all counts are non-zero then normalize by the total number of <NUM>-mers in the sequence and log<NUM>-transform the normalized values.

The contig- and read-level polymerase kinetics scores are calculated using the inter-pulse duration (IPD) values provided in the SMRT sequencing reads. Subread normalization, done by log-transforming the ratio of each subread IPD value to the mean of all IPD values in the subread, corrects for any potential slowing of polymerase kinetics over the course of an entire read (which can consists of multiple subreads). Each normalized IPD (nIPD) value in the subread is calculated as follows:
<MAT>
where the subread is N bases long and therefore contains N IPD values. To calculate the observed read-level methylation score (Ro) for motif i on read j, <MAT>, the mean of all nIPD values was taken from all sites of motif i across all subreads of read j:
<MAT>
where each of the S subreads in the read contains Ms motif sites. Longer subreads typically contain more distinct sites of a given motif and generate more reliable methylation scores.

Kinetic variation in the polymerase activity exists even in the absence of methylated bases and is highly correlated with the local nucleotide context surrounding the polymerase as it processes along the template. To account for this baseline variation and remove it from the final methylation score, a corresponding set of control kinetics scores, <MAT> was subtracted from the observed kinetics scores, <MAT>. These control kinetics scores are motif-matched and calculated similar to <MAT> using a sampling of SMRT sequencing unaligned reads (N=<NUM>,<NUM>) known to be free of any methylation:
<MAT>.

As no methylated motifs were detected after sequencing an isolate of Ruminococcus gnavus, this data served as the non-methylated control set for calculating values of <MAT>. These non-methylated control values are used for the motif filtering procedure, but not for the final calculation of methylation profiles. Because the dimensionality reduction with t-SNE calculates a Euclidian distance between two points (i.e. two methylation profiles), the subtraction of a constant (control) vector from both methylation profiles has no effect on their pairwise distances.

Contig-level methylation scores (C) for motif i on contig j, Cij, are calculated in a similar manner. The difference is that the scores take into account not just the subreads from a single read, but rather all subreads that align to the contig:
<MAT>
where each of the S* subreads that align to the contig contain Ms motif sites. Similar to the read-level methylation scores, matching control kinetics scores, <MAT>, are generated using a sample of aligned reads (N=<NUM>,<NUM>) known to be free of methylation and subtracted from the observed kinetics scores, <MAT>, in order to remove the baseline kinetics variation stemming from local sequence context:
<MAT>.

As with the read-level methylation scoring, non-methylated control values are used only during the motif filtering procedure but not in the final contig-level methylation scores. Much like the read-level methylation assessment, the reliability of the motif score on a contig increases with the number of motif sites on the contig. Typically, short motifs are present at higher density in the genome than longer, more complex motifs, although exceptions to this rule exist. Therefore, while even the shortest contigs in an assembly are able to return reliable methylation scores for short motifs, longer contigs are usually required to accurately assess the methylation status of more complex motifs. A default methylation score of zero is assigned if no instances of the motif occur on the read or contig.

The optional parameter --cross_cov_bins in the mBin program accepts a file containing contig assignments to bins (in the format contig_name, bin_id) identified from coverage- and composition-based binning tools. If this parameter is specified, the IPD values used to calculate each contig-level methylation score are aggregated based on binning assignment and bin-level methylation scores are calculated.

An initial motif-filtering step is necessary to reduce the space of motifs down to only those that have a significant methylation score in the metagenomic mixture. First, due to memory considerations and because a motif could theoretically describe any arbitrary string of bases, the maximum motif length and allowable base configuration of motifs was defined in the initial query space. All possible 4mers, 5mers, and 6mers were considered, for a total of <NUM>,<NUM> contiguous motifs. For bipartite motifs, where a string of non-specific Ns was bookended by sets of specific bases (e.g. CCANNNNNNCAT (SEQ ID NO: <NUM>)), several common configurations often found in prokaryotes were considered. All combinations of the following were considered: <NUM> or <NUM> specific bases (beginning), <NUM> or <NUM> non-specific Ns (middle), and <NUM> or <NUM> specific bases (end). This adds an additional <NUM>,<NUM> possible bipartite motifs to space of motifs to consider for the initial filtering step, for a total of <NUM>,<NUM> motifs. The exact same method can be used to further incorporate <NUM>-mer and <NUM>-mer motifs.

Next, the motif query space was dramatically reduce by randomly sampling a small number of reads (N=<NUM>,<NUM>) from the mixture and removing from further analysis all motifs that do not return a methylation score above a chosen threshold (<NUM>) on at least one contig in the assembly (or on at least twenty unaligned reads for read-level binning). Despite choosing a lenient threshold to include many variations of the truly modified motif, this typically reduces the number of motifs to be included in the further analysis by multiple orders of magnitude. A further step searches for multiple specifications representing a single degenerate motif that, if identified, replaces the individual specifications in the final set of motifs. The remaining motifs need not exactly match the most parsimonious versions of the methylated motifs, but they nonetheless will carry some methylation signature that is useful for binning the sequences through subsequent dimensionality reduction analysis. Put another way, the precise number of motifs that remain after filtering is not usually critically important as long as the set of remaining motifs captures the most significant differences between methylation profiles. This property contrasts with existing methods for methylation motif discovery that attempt to identify the single most parsimonious version of a motif.

The combination of k-mer frequency and methylation scores used to segregate similar species and strains in the combined infant gut microbiome samples A and B (<FIG>) was done by z-score transforming both feature matrices after each had been reduced to 2D using t-SNE. The two 2D matrices of z-scores were then combined and the resulting 4D matrix of z-scores was subjected to a second round of t-SNE to get a final 2D matrix.

To assess the sequence similarity between two reference genomes, average nucleotide identity (ANI) was calculated using the web-based portal at the world wide web enve-omics.

A database of <NUM> reference genomes isolated from the mouse gut was compiled from four recent studies. Blastn was first run to identify which of the reference sequences had significant matches with the contigs in the nine bins identified using methylation profiles. Significant hits were considered to be alignments >100bp in length with ><NUM>% identity. For each bin, the reference genomes were ranked based on the percentage of the total binned contig sequences that were covered by a significant hit with the reference. The mummer package was then used to align the highest ranked matching references to the contigs in each bin and visualized the alignments (<FIG>) with the mummer package.

After aligning reads from <NUM> publicly available mouse gut microbiome sequencing data sets to the largest contigs in each of the nine methylation bins, coverage values were normalized according to the standard normalization procedures employed by CONCOCT. To exclude regions where high sequence similarity with other contigs might result in ambiguous mapping and unreliable coverage values, each contig was divided into 10kb subsequences and excluded any subsequences that displayed any alignments using nucmer. Mean coverage values were calculated for the unique remaining subsequences and these were used to construct the coverage profiles across all <NUM> samples (<FIG>).

The long reads used in this study often result in a bacterial genome being represented by a small number of very large contigs. The t-SNE dimensionality reduction algorithm places data points in low-dimensional space based on the local similarities in the original high-dimensional space. Species with few large contigs that are represented by only a few points in the high-dimensional space do not contribute significantly to the objective function of the t-SNE algorithm. To adjust for this bias from different contig sizes, a length-weighted representation of all large contigs over <NUM> kbp in length was use so that each large contig is represented in the matrix of features not by one row, but by N rows, where N is the contig length divided by <NUM> kbp. The features (column values) for each <NUM> kbp sub-contig, either k-mer frequency or methylation scores, are the same values that were computed for the original large contig.

In order to assess the power of methylation scores to distinguish a contig methylated at a motif sites (case) from a contig that is not methylated at that motif (control), <NUM>,<NUM> normalized IPD (nIPD) values were sampled from GATC sites on each of two large assembled contigs from the mixture of eight bacterial species. The case was the <NUM>. 6Mb contig representing the E. coli chromosome, while the second <NUM>. 7Mb contig (control) represents a large assembled portion of the R. gnavus genome, which does not contain any methylated motifs based on SMRT sequencing data (see Table <NUM>). The two sets of <NUM>,<NUM> nIPD values were then used as pools from which to sample <NUM>, <NUM>, <NUM>, and <NUM> values for both the case and control. The nIPD values were used to construct methylation scores for GATC on both the case and control contigs, for each of the four specified nIPD sampling numbers (<NUM>, <NUM>, <NUM>, and <NUM>). This process was repeated <NUM>,<NUM> times to create a receiver operating characteristic (ROC) curve (<FIG>) showing the effect of the number of nIPD values on creating methylation scores that can distinguish a methylated contig/motif from a non-methylated contig/motif.

When calculating the Euclidian distance between a plasmid and the chromosome of its host bacterium, the largest chromosome was selected when a bacterium contained more than one chromosome. The empirical distribution of Euclidian distances between the plasmids and randomly selected bacteria was constructed by iterating over all plasmids in REBASE, randomly selecting a bacterium for each plasmid, and computing the distance between the plasmid <NUM>-mer frequency vector and that of the largest chromosome of the selected bacterium.

Methylation motifs were gathered for each of the <NUM> SMRT sequenced bacterial genomes stored in the REBASE database and mock communities of N species were constructed, where N = <NUM>, <NUM>, <NUM>,. , <NUM> and each community was created <NUM>,<NUM> times by randomly selecting from the <NUM> organisms. For each mock community, the methylation motifs for each constituent organism were analyzed and number of organisms with a unique methylome in the community was returned, reported as the fraction of total organisms in the community. Multiple curves in <FIG> represent the different results obtained by varying the multi-strain content of the mock communities. The same procedure was again used to analyze only those <NUM> organisms in REBASE that are known to host at least one plasmid sequence. Mock communities of N species were again constructed, where N = <NUM>, <NUM>, <NUM> and each community was created <NUM>,<NUM> times by randomly selecting from the <NUM> organisms. Multiple curves in <FIG> represent the different results obtained by varying the multi-strain content of the mock communities.

For each SMRT sequenced genome in the REBASE database, <NUM> random sequences of length L were simulated, where L = <NUM>, <NUM>, <NUM>,. Given the known methylation motifs for each genome, the number of sequences containing the motifs was returned, reported as the fraction of the <NUM> total simulated sequences. Multiple curves in <FIG> represent the different results obtained by varying the percentage of the genome's methylation motifs that are required to be present on each sequence. For instance, the <NUM>% curve represents the number of simulated sequences that contain at least one instance of at least three quarters of the genome's total set of methylation motifs.

In each methylation bin, the reads aligning to each binned contig were reassembled with the HGAP3 assembler using a genomeSize parameter modified to reflect the total number of contig bases in each bin.

A combination of two methods was used to identify circular contigs in metagenomic assemblies: (<NUM>) a custom script aligned the 20kb sequences at the beginning and end of contigs to look for evidence of circularization, and (<NUM>) the freely available program Circlator was used with default parameters. Contigs identified as circularized were then manually checked using Gepard to look for visual evidence of circularization, as opposed to signs of mis-assembly.

Small (<200kb) contigs were classified as conjugative transposons if they contained at least five genes encoding conjugative transposon-related genes. The contigs from each methylation bin (#<NUM>-<NUM>) were annotated by submission to the RAST server.

Eight species synthetic mixture. SMRT reads were obtained separately from eight individual bacterial species (Table <NUM>) and the reads were mixed, without any labeling, by combining one SMRT cell of sequencing from each species to create a synthetic metagenomic mixture at similar relative abundances. Read labels were applied for evaluation purposes only after all binning procedures were completed.

Human Microbiome Project Mock Community B. Equimolar amounts of genomic DNA were extracted from twenty different species (Table <NUM>) then combined and sequenced using a Pacific Biosciences RSII instrument. The <NUM> SMRT cells of reads are publicly available at this GitHub link on the world wide web at github. com/PacificBiosciences/DevNet/wiki/Human_Microbiome_Project_MockB_Shotgun. In order to simulate a more realistic mixture with widely varying relative abundances, the raw sequencing reads were downsampled to impose relative species abundances that follow a natural log decay curve (<FIG>). The species identity for all reads were first determined by aligning the reads to reference assemblies for each of the twenty species. After determining the species mappings for all reads (excluding those with ambiguous alignments), reads from each species were then selected to impose the desired relative abundances. The alignment and labeling procedures were used strictly for data downsampling and were not part of the read-level binning procedure.

Multi-strain mixture of Helicobacter pylori. Two strains of H. pylori, str. <NUM> and str. J99, were sequenced separately using a Pacific Biosciences RSII instrument as part of a previous study. In order to create a multi-strain mixture, reads from one SMRT cell per strain were combined. These strain-specific sets of reads were downsampled using their SMRT cell labels then combined to a mixture containing both strains at 150x coverage (Table <NUM>). Binning procedures did not use any information from the labels.

Multi-strain mixture of Escherichia coli. Three strains of E. coli, BAA-<NUM> O26:H11, BAA-<NUM> O103:H11, and BAA-<NUM> O111, were sequenced separately using a Pacific Biosciences RSII instrument (see See Examples section entitled Three E. coli strains for synthetic mixture). The synthetic, multi-strain mixture was created by combining a single SMRT cell from each of these separate sequencing runs (Table <NUM>). Binning procedures did not use any information from the labels.

The microbial DNA HM-277D was obtained from BEI Resources and was sequenced in a previous study by Kuleshov et al. using the Illumina TruSeq protocol. These sequencing results were downloaded for the current study using the SRA accession code SRR2822454.

Both synthetic long reads and SMRT reads were aligned to the <NUM> reference sequences of the genomes contained in the HMP Mock Community B. The synthetic long reads were aligned using the SMRT read aligner blasr with default parameters and "-bestn <NUM> -sam" options. The synthetic long reads were aligned using bwa-mem with default parameters.

bam files containing the aligned synthetic long reads and SMRT reads for the <NUM> species in the HMP Mock Community B were analyzed to count the total number of aligned bases in each. For each reference, the smaller number of aligned bases was chosen as the target number of aligned bases and the file with the larger number of aligned bases was selected for downsampling. The target fraction is calculated by dividing the target number of aligned bases by the original number of bases. The following samtools command was used to generate the downsampled file:
samtools view -s <NUM>. [target_frac] -h -b original. bam > downsampled.

The results of this downsampling are summarized in Table <NUM>.

DNA was isolated from stool samples taken from two Finnish children. The donor of Sample A (containing B. <NUM>) was <NUM> months of age, while Sample B (containing B. <NUM>) was obtained from child at <NUM> months of age. Full details on sample isolation and DNA extraction are provided by Leonard et al. A summary of the SMRT sequencing statistics can be found in Table <NUM>.

The high-dimensional matrix of features (e.g. k-mer frequencies, methylation scores, or a combination) for all sequences was subjected to the Barnes-Hut implementation of t-distributed stochastic neighbor embedding (t-SNE). The Barnes-Hut approximation of t-SNE reduces the computational complexity from <IMG>(N<NUM>) to <IMG>(N log N), making it feasible to generate 2D maps of hundreds of thousands of metagenomic sequences containing hundreds of features. All runs used the default parameters for perplexity (<NUM>) and theta (<NUM>).

All metagenomic assemblies in this study used the hierarchical genome-assembly process (HGAP3). With the exception of the parameter specifying the expected genome size to be assembled, all default parameters were used. The expected genome size parameter is used to determine the optimum number of long seed reads and was adjusted based on the expected complexity of the metagenome. Specifically, the genome size was set to 40Mb for the synthetic mixture of eight bacterial species assembly, 66Mb for the <NUM>-member HMP assembly, 20Mb for the combined infant gut microbiome samples A and B assembly, <NUM>. 6Mb for the combined and separate H. pylori strain assemblies, and 20Mb for the infant gut microbiome sample A assembly.

Kraken version <NUM>. <NUM>-beta was configured to use two databases. The database used to annotate sequences from the Human Microbiome Project (HMP) Mock Community B consisted of reference sequences for the twenty known species included in the mock community (Table <NUM>). All other Kraken annotations used a database consisting of the RefSeq complete set of bacterial/archaeal genomes (using "--download-library bacteria") and draft assemblies of five Bacteroides dorei strains. Database construction from these libraries and all Kraken annotations used default parameters.

In the infant gut microbiome t-SNE maps showing the combined assemblies of samples A and B (<FIG>), all contigs other than those labeled as belonging to B. dorei were annotated using Kraken. The contigs belonging to the two B. dorei strains, however, were manually labeled by first aligning the reads from the combined samples to the fully assembled references for each B. dorei strain (strain <NUM>: CP007619; strain <NUM>: CP008741). The contig-labeling assignments were determined by examining the reads aligning to the either of the B. dorei references and counting how many of these reads aligned to each of the assembled contigs. For example, if the majority of the reads aligning to a contig aligned to the strain <NUM> reference, the contig was labeled as belonging to strain <NUM>. However, if the majority aligned to the strain <NUM> reference, the contig was labeled as belonging to strain <NUM>.

Claim 1:
A method of deconvoluting genomes of prokaryotic organisms in a microbiome sample, said method comprising the steps of:
a) sequencing nucleic acids of the prokaryotic organisms comprised in a sample obtained from a subject using single-molecule long reads sequencing technology, wherein the sequencing comprises the step of identifying methylated nucleotides, and at least one of the steps of:
i. sequencing single molecule reads of nucleic acids;
ii. assembling contigs from single molecule reads of the nucleic acids; and
b) assigning a methylation score reflecting the extent of methylation for sequence motifs of the nucleic acids on the assembled contig and/or the single molecule read;
c) applying motif filtering to identify sequence motifs with methylation scores indicating methylation on the assembled contigs and/or the single molecule reads;
d) determining nucleic acid methylation profiles of the assembled contigs or the single molecule reads in the microbiome sample based on motifs identified in step (c);
e) separating the assembled contigs and/or the single molecule reads into bins corresponding to distinct prokaryotic organisms based on the methylation profiles of step (d);
f) assembling the bins of step (e), thereby obtaining assembled genomes of the distinct bacterial organisms in the microbiome sample,
thereby deconvoluting genomes of the prokaryotic organisms in a microbiome sample.