Patent Application
20200102599
The present invention relates to the determination of growth rates and viability from genome coverage, and more specifically, to assessing the growth and hazard potential of one or more species of microorganisms in a multi-species sample.
Metagenomic or metatranscriptomic sequencing datasets allow the determination of the presence or absence of a given organism but do not provide a means of assessing the viability, rate of growth, and/or rate of replication of an organism at the time of sampling. This deficiency can preclude the determination of a hazard potential associated with the organism.
Many methods exist for inferring the different life states of bacteria (Davis, C., “Enumeration of probiotic strains: Review of culture-dependent and alternative techniques to quantify viable bacteria,” Journal of Microbiological Methods, 103, 9-17 (2014); Singer et al., “Capturing the genetic makeup of the active microbiome in situ,” The ISME Journal (2017), 11, 1949-1963). However, no currently implemented method exists for utilizing sequencing data for this purpose.
Most methods of viability determination are based on quantitative PCR (polymerase chain reaction) of pre-rRNA sequences before and after stimulation of the sample with nutrients (Cangelosi et al., “Dead or Alive: Molecular Assessment of Microbial Viability,” Applied and Environmental Microbiology 80, 5884-5891 (2014)). These methods cannot be applied to passively collected data and are tailored to assay only a select set of species, as PCR probes need to be synthesized ahead of time (U.S. Pat. No. 9,115,407 B2). The same probe-related bias persists with methods that are based on fluorescence (e.g., fluorescence activated cell sorting (FACS)).
The advent of high-throughput next-generation sequencing has allowed for sampling the entire microbial community present in a food, environmental, or patient sample. Central to this approach is broad fragmentation of nucleic acids extracted from bacterial cells, referred to herein as “shotgun sequencing.” For each sequencing experiment, researchers have the choice of extracting either RNA or DNA as the primary template of interest. Both can provide insight into the taxonomic and functional characteristics of the microbes, but they skew differently in terms of their coverage of various parts of the investigated genomes and in terms of answering different biological questions.
Computational methods in the bioinformatics field have been developed to align the short reads generated by shotgun sequencing to high quality assembled reference genomes. This is done by minimizing the disagreement in the nucleotides of the two biological sequences being compared.
Alignments of short reads to reference genomes can be seen as a histogram of counts of respective nucleotides along the entire length of the genome template. Commonly, the differential peak heights of this histogram are informative of the sample either having an alternate copy number of a segment compared to the reference or having differential expression for DNA and RNA, respectively (Korem et al., “Growth dynamics of gut microbiota in health and disease inferred from single metagenomic samples,” Science, 2015, 349, 1101-1106).
A general comparison of metagenomic and metatranscriptomic sequence data of a small set of samples validates that some genes are overexpressed compared to their genomic copy number (e.g., ribosomal proteins), suggesting more active replication, while others are underexpressed (e.g., sporulation genes) suggesting senescence (Franzosa et al., “Relating the metatranscriptome and metagenome of the human gut,” Proc. Natl. Acad. Sci. U.S.A., 2014, 111, E2329-38). However, this work has not been used to further strengthen the link between differences in expression and microbial activity, nor has it been used to develop a machine learning method that is capable of predicting replication rates.
In addition to studies in bacteria, some yeast genes are predictive of yeast growth rates measured experimentally (Airoldi et al., “Predicting cellular growth from gene expression signatures,” PLoS Computational Biology, 2009, 5, e1000257). This validates that transcription is altered by the underlying changes in the number of DNA copies. More complex dynamics arise than just a linear amplification of transcript copies, especially in eukaryotes (Bar-Ziv et al., “Dealing with Gene-Dosage Imbalance during S Phase,” Trends in Genetics, 2016, 32(11), 717-723), necessitating a more careful tuning of predictive sequences.
Methods of identifying high growth species in a community sample are needed that do not rely on culturing the samples and are adaptable to machine learning.
Accordingly, a method is disclosed, comprising:
sequencing nucleic acids of a non-cultured sample, the sample containing one or more prokaryotic species, thereby producing a plurality of sequences of the nucleic acids, the sequences comprising i) metagenomic sequences corresponding to DNA and ii) metatranscriptomic sequences corresponding to RNA;
mapping the sequences to reference genomes of the one or more prokaryotic species, thereby obtaining i) a list of taxonomically identified prokaryotic species of the sample and ii) coverages of the sequences with respect to two or more regions of the reference genomes active during growth of the prokaryotic species, the regions designated active regions; and
calculating a growth rate and a viability of the one or more prokaryotic species of the sample based on the coverages.
Also disclosed is a system, comprising one or more computer processor circuits configured and arranged to:
sequence nucleic acids of a non-cultured sample, the sample containing one or more prokaryotic species, thereby producing a plurality of sequences of the nucleic acids, the sequences comprising i) metagenomic sequences corresponding to DNA and ii) metatranscriptomic sequences corresponding to RNA;
map the sequences to reference genomes of the one or more prokaryotic species, thereby obtaining i) a list of taxonomically identified prokaryotic species of the sample and ii) coverages of the sequences with respect to two or more regions of the reference genomes active during growth of the prokaryotic species, the regions designated active regions; and
calculate a growth rate and a viability of the one or more prokaryotic species of the sample based on the coverages.
Further disclosed is a computer program product for calculating growth rates and viability of prokaryotic species of a non-cultured sample, the sample comprising one or more prokaryotic species, the computer program product comprising a non-transitory computer readable storage medium having program code embodied therewith, the program code executable by a processing device to perform a method comprising:
sequencing nucleic acids of a non-cultured sample, the sample containing one or more prokaryotic species, thereby producing a plurality of sequences of the nucleic acids, the sequences comprising i) metagenomic sequences corresponding to DNA and ii) metatranscriptomic sequences corresponding to RNA;
mapping the sequences to reference genomes of one or more prokaryotic species, thereby obtaining i) a list of taxonomically identified prokaryotic species of the sample and ii) coverages of the sequences with respect to two or more regions of the reference genomes active during growth of the prokaryotic species, the regions designated active regions; and
calculating the growth rate and the viability of the one or more prokaryotic species of the sample based on the coverages.
The above-described and other features and advantages of the present invention will be appreciated and understood by those skilled in the art from the following detailed description, drawings, and appended claims.
A method is disclosed for predicting viability and growth rates of at least one prokaryotic species of a working sample (e.g., environmental, food, patient) utilizing nucleic acid (i.e., DNA and/or RNA) sequencing data. The method is not dependent on amplification of the nucleic acid, hybridization of nucleic acid, or physical properties of the nucleic acid once the nucleic acid has been sequenced. The method is general to any prokaryotic organism (i.e., archaea, bacteria), as opposed to a well-characterized single prokaryotic species. The method does not require culturing of the sample or any implicit annotation of the genome of interest. The method is also applicable to samples previously frozen/shipped from a remote location. The method is adaptable to machine learning.
As bacteria reproduce asexually, the increase in their population is accomplished by binary fission, which features among many things the duplication of the circular genome. In general, this occurs via a symmetric process that starts at a known locus, referred to in the art as the origin of replication, and proceeds in both directions along the template genome. This locus has certain identifiable features (e.g., high adenine-thymine content) that promote the replication “machinery” (e.g., DNA polymerase, circular sliding clamp, pentameric clamp loader, helicase, primase, and single-strand binding protein) binding to it.
During the replication process, a transient difference can occur in the copy number of DNA regions close to the origin of replication compared to more distant areas of the genome. The copy number imbalance gathered from a working sample containing multiple bacteria has been found to be useful to infer statistically which species of the working sample are actively replicating and to assign viability scores based on the relative activity of each of the species. The analysis to predict viability is therefore linked to present replication activity rather than to natural selection to regions surrounding the loci (i.e., population evolution).
Other loci of the circular genome can also serve as indicators of active replication. For example, in order to replicate the entire genome and successfully reproduce, bacteria produce more RNA transcripts (and subsequently proteins) related to making the replication machinery and the substrates and energy required therefor. These loci share homology among multiple bacterial genomes, as the need to reproduce and hence replicate one's genome is a universal feature of living organisms. Additionally, bacterial genomes prioritize genes vital for replication and transcription to be clustered closer to the origin of replication (Couturier et al., “Replication-associated gene dosage effects shape the genomes of fast-growing bacteria but only for transcription and translation genes,” Mol. Microbiol., 2006, 59, 1506-1518). The disclosed method utilizes these loci, selecting the most predictive loci in samples containing multiple species.
The following definitions are applicable.
A “contig” is a set of overlapping DNA sequences that together represent a consensus sequence of DNA or a region thereof.
A “consensus sequence” is the calculated order of the most frequent residues found at each position in a sequence alignment.
“Copy number” means the number of copies of a gene or plasmid within a genome. The copy number can vary from individual to individual.
“Coverage” or “depth of coverage” is the number of times a given sequence from a genome is represented in the set of sequences derived from that genome.
“DNA” is deoxyribonucleic acid. DNA can be double-stranded or single-stranded.
A “gene” is the basic unit of heredity, a linear sequence of nucleotides along a segment of DNA that provides the coded instructions for synthesis of RNA, which, when translated into protein, leads to the expression of a hereditary trait.
A “genome” is the total genetic content of a microorganism. In the case of bacteria, the genome is the DNA.
A “ground truth dataset” is a dataset formed by direct observation (measured data) as opposed to data obtained by inference or assumption.
Herein “high-throughput sequencing” (HTS) is any method of sequencing a nucleic acid that is highly parallel and does not involve cloning the nucleic acid. A genome or metagenome is cut into a large number of fragments, and the fragments are sequenced in parallel.
“Homology” refers to the similarity of sequences (e.g., DNA, RNA, Protein, etc.) arising from a common ancestry.
“Hybridization” is the formation of a double-stranded helix from single-stranded complimentary pairs of DNA and/or RNA by annealing.
The term “k-mer” means a sub-sequence of a read obtained through DNA sequencing having k number of nucleotide base units, where k is a positive whole number greater than 1.
Herein, a “database” is an electronic file for storing and retrieving data. Databases are also referred to herein as data tables. Data tables comprise rows and columns (i.e., fields) of data. The rows are formally called tuples or records. A data table comprises one or more records, each record comprising one or more defined fields having respective defined data types (e.g., text, numeric, date, time, memo, and so on) and defined field lengths where applicable. A working data table comprises at least one record containing data in one or more fields of the record. The data tables are located on data storage devices, which can be remote or local relative to the user input/output devices. A “database system” comprises at least one data table and a database management software program for managing the storage and retrieval of data to and from the data tables. The database management programs can be remote or local relative to the data tables and/or the end user. A Relational Database Management System (RDBMS) is a database management system (DBMS) that uses relational techniques for storing and retrieving data using data tables. A relational database system can have many data tables, and each data table can have multiple records and multiple fields within each record. A data table in a relational database system can be accessed using an index. An index is an ordered set of references (e.g., pointers) to the records or rows in a data table. The index is used to access each record in the file using a key (e.g., one or more of the fields of the record or attributes of the row). Without an index, finding information in a large data table would require a resource-intensive time-consuming scan (e.g., linearly) of each record of a table. Indexes provide a faster alternate technique of accessing data contained in one or more data tables that are linked by a common key. Users can create indexes on a table after the table is built. An index is based on one or more columns (fields) of a given table.
A “k-mer database” is a database in which a given record comprises a field for storing a k-mer of a nucleic acid sequence of one or more organisms. Another field of the record stores a taxonomic ID that associates the k-mer to a lowest common ancestor node (LCA) of a taxonomic tree. As will be described below in more detail, other fields of the record can store reference IDs to a reference taxonomy. Still other fields of the record can store metadata associated with the k-mer and/or the nucleic acid sequence from which the k-mer originated.
A “locus” (plural loci) is a position on a genome of a region of interest (e.g., gene, regulatory element, origin of replication).
A “metagenome” is all the genetic information from a microbiome (see microbiome below).
“Metagenomics” is the analysis or study of metagenomes.
“Metatranscriptome” is the collection of all RNA transcripts from a microbiome (see microbiome below)
“Metatranscriptomics” is the analysis or study of metatranscriptomes.
A “microbiome” is a community of microorganisms that inhabit a particular environment (e.g., microbes of the human gut), or a sample taken therefrom.
“Origin of replication” is the locus at which DNA replication begins.
A “plasmid” is a self-replicating extrachromosomal circular DNA that replicates independently of the bacterial chromosome and carries genes for functions not essential for growth.
“RNA” is ribonucleic acid.
“mRNA” refers to messenger RNA. The mRNA codes for amino acid sequences of proteins.
“rRNA” refers to ribosomal RNA.
“tRNA” refers to transfer RNA. A tRNA transports a specific amino acid to a ribosome for synthesis of a protein.
An “RNA transcript” is an RNA produced through the process of transcription of DNA.
“Sample” means any sample capable of undergoing analysis using the disclosed methods. The sample can contain one or more microorganisms. Non-limiting samples include water samples obtained from tap water, lakes, streams, field runoff, and sewage; swabbed samples from contact surfaces (e.g., building surfaces, countertops, furniture, utensils, clinical instruments, computer hardware, cell phones, door handles, doors, windows, screens, cabinets, cabinet doors, sinks, faucet); animal samples (e.g., blood, blood plasma, serum, cells, a cellular extract, a cellular aspirate, expectorant, sputum, saliva, mucous, urine, sweat, tears); and samples obtained from food, food-handling equipment, and surfaces contacted by food. The samples can be a solid or liquid containing water or no water.
“Sequencing” refers to a process of determining the precise order of residues (e.g., nucleotides) in a nucleic acid (e.g., DNA, RNA).
“Sequence” refers to a fragment of a nucleic acid (RNA or DNA) that has been sequenced (i.e., the order of base pairs is known).
A “sequence read” or “read” is a finite length or fragment of a nucleic acid that is output by a sequencing instrument. For example, a read from an Illumina sequencer is 100-150 base pairs in length today. Sequencing may also be done on “paired end” reads where two reads are connected by a spacer (that is not read), increasing the effective read length to 300 or more and covering a larger region of the genome.
A “sequence alignment” is a way of arranging sequences to identify regions of similarity, which may be a consequence of functional, structural, or evolutionary relationships between the sequences.
“Shotgun sequencing” is a quasi-random process in which a nucleic acid is broken up into many random smaller fragments that are individually sequenced. The sequences are ordered based on overlapping regions of genetic code, and reassembled into the complete sequence of the nucleic acid.
“Taxonomy” is a biological scheme of classification of organisms. Herein, for prokaryotes, the heirarchy is domain, kingdom, division, phylum, class, order, family, genus, species, sub-species, and strain.
A “taxonomic tree” herein is a data structure for classifying organisms. The taxonomic tree comprises nodes (i.e., taxa, singular taxon) that are grouped into “parent nodes” linked to “child nodes”. Parent nodes are depicted above child nodes in the tree diagram. Child nodes are taxonomic descendants of parent nodes. For example, a genus (parent node) can be linked to two or more species (child nodes). The taxonomic tree can be rooted (i.e., known ancestral root) or unrooted (i.e., unknown ancestral root), bifurcating (i.e., two child nodes per parent node) or multi-furcating (i.e., more than two child nodes per parent node). Typically, the taxonomic tree is in the form of a “binary tree” (i.e., each parent node has two child nodes). A “leaf node” is a child node having no descendants (e.g., the species of a genus). In the self-consistent taxonomy, each leaf node has one genome. “Internal nodes” are all nodes other than the leaf nodes.
“Transcription” is the process of forming an RNA from a DNA template.
The abbreviation “bp” means “base-pair” (e.g., a read of 100-bp means that one DNA read has 100 nucleotides in the polymer chain.
“Miscalling” refers to a sequencing error where a nucleotide in a sequence read is different from the true nucleotide.
A quality value is an assigned value given to each nucleotide in a sequence read that reflects the likelihood of miscalling the nucleotide. The higher the quality value is, the lower the likelihood of miscalling.
A “reference genome” is a genome from the same species or close species that has already been sequenced.
“Mapping” a sequence read is a process of finding the position or coordinate of a sequence read on the reference genome.
A “perfect match prefix” is a k-mer of a sequence read that is identical to, or a perfect match to some equal-length k-mer(s) of the reference genome. The k-mer of the sequence read is used to initially anchor the sequence read on the reference genome.
Base substitution: After a sequence read is mapped to the reference genome, certain bases are different from the corresponding bases on the reference genome.
Insertion: Compared with the reference genome, some continuous bases are inserted between two adjacent bases on the sequence read.
Deletion: Compared with the reference genome, the sequence read loses some continuous bases.
INDEL: an insertion or deletion in a read when trying to find the best alignment of a read to a reference genome.
The disclosed method is composed of two parts. Preferably, each part is performed using a computer system. In Part 1, depicted in the process flow diagram of FIG. 1, which includes
In Part 2, depicted in the process flow diagram of
More specifically, the reference database is used as a component of computer system to perform a statistical analysis of the metagenomic coverage data of the working sample, from which is generated a viability and/or growth score of one or more microorganisms therein. Both viability and growth rate predictions are based on the alignment of sequenced nucleic acid fragments over the active regions defined in the reference database. The computer system includes a computer program and end user interface, described in more detail further below.
Parts 1 and 2 are described in more detail in the following sections.
To construct a reference database, the population growth rate of one of more bacterial control samples are measured with time (
To create a control sample, the user selects a bacterial isolate and measures an initial quantity or concentration of bacterial cells. Non-limiting methods of measuring bacterial growth rates include i) optical density (OD) measurements, ii) CFU (colony forming units) counting assays, and iii) quantitative PCR. These methods can be used singularly or in combination. The initial quantity of bacterial cells is then added to a liquid growth medium and incubated at a user selected temperature that supports growth (25° C. or 37° C. are most common). A standard bacterial growth assay is performed as follows. Samples of the liquid culture are taken at user-defined times that should include the 3 phases of bacterial growth: the lag phase, log phase, and stationary phase. It should be understood that different bacteria can have different rates of growth that affect the choice of time interval between samples. The bacterial growth of each of the samples is then characterized using one or a combination of the above-described methods of determining the concentration of bacterial cells. The samples are also sequenced (below). As an example, Escherichia coli completes all three growth phases in approximately 20 hours at 37° C. The user may elect to sample the liquid culture once every two hours for preparation for sequencing. After the 20 hours and 11 time points (0-20), the user has samples from each phase of growth.
At each time interval, the nucleic acids (i.e., RNA and/or DNA) are extracted from the control samples. The extracted nucleic acids are then subjected to high-throughput sequencing (
Optionally, quality control of the raw sequences can be performed by purging sequences of poor quality, removing contaminating sequences introduced by the sequence methodology, and/or removing (trimming) sequences of low complexity (
Optionally, the raw sequences (or the clean sequences) can be assembled by software to form what is referred to herein as raw contigs or clean contigs, respectively (not shown). If using assembled contigs, the assembly must represent complete coverage of the genome or substantially complete coverage of the genome.
Each of the time-dependent sets of raw sequences, clean sequences, raw contigs, and/or clean contigs of the control sample, referred to herein collectively as control reads, are then mapped (aligned) to a known reference genome assembly of the microorganism of the control sample (
The per-position-coverage and windows of the set of mapped control reads are then calculated (
Optionally, the mapped control reads can be filtered based on a coverage threshold (not shown), effectively removing them.
Optionally, the coverage data of the mapped control reads can be smoothed using LOESS (locally weighted polynomial regression) or another genome feature-aware method to remove unlikely spikes of coverage (not shown). More specifically, LOESS is a method for capturing diverse score-length relationships in sequence data and correcting sequence scores for length-dependent artifacts.
The time-dependent coverages of the mapped control reads are then analyzed in order to identify all “active regions” of the genomes that are associated with growth and viability of the microorganism of the control sample (
To evaluate the performance of different methods at predicting replication rate, standard evaluation criteria can be used such as, for example, mean squared error or mean absolute error. To assess performance at predicting viability using the ground truth data generated from the control experiment and comparing to the predicted features from a model trained on that data, standard metrics can be used such as misclassification rate, sensitivity, specificity, and other techniques known in the art.
As used herein, the term “specificity” refers to a measure of correctly predicting growth rate and viability by identifying genomic regions (based on coverages) that should not be included in the estimation of growth and viability (i.e., a measure of true negatives). The term “sensitivity” refers to a measure of correctly predicting growth rate and viability by identifying genomic regions (based on coverages) that should be included in the estimation of growth and viability (i.e., a measure of true positives).
The reference database is then constructed as a set of sequences (
The reference database is used in Part 2 to predict growth rates of one or more microorganism of a working sample.
In Part 2 of the method, a working sample is collected (
Nucleic acids (i.e., RNA and/or DNA) are extracted from the working sample without prior nutrient stimulation of the working sample or use of other DNA amplification techniques. The extracted nucleic acids are then subjected to high throughput sequencing (
Optionally, quality control of the raw working sequences can be performed as described above by removing sequence data of poor quality introduced by the sequencing methodology, removing contaminating sequences, and/or removing low complexity sequences (
Optionally, the raw working sequences or clean working sequences are assembled to form raw working contigs or clean working contigs, respectively (not shown).
The raw working sequences, clean working sequences, raw working contigs, or clean working contigs, collectively referred to herein as working reads, can then be mapped to the collection of known taxonomically classified whole reference genomes of prokaryotes (
The mapping also determines the per-position-coverage (
In the next step (
The reference database formed in Part 1 can be further refined using the positive and/or negative outcomes of the predicted growth rates and viability scores in order to improve the likelihood and accuracy of future successful predictions.
Microorganisms include bacteria, fungi, viruses, protozoans, and parasites.
Exemplary non-limiting bacterial species include Acetobacter aurantius, Acinetobacter baumannii, Actinomyces israelii, Agrobacterium radiobacter, Agrobacterium tumefaciens, Akkermansia muciniphila, Alistipes_onderdonkii, Alistipes putredinis, Anaplasma phagocytophilum, Azorhizobium caulinodans, Azotobacter vinelandii, Bacillus anthracis, Bacillus brevis, Bacillus cereus, Bacillus fusiformis, Bacillus licheniformis, Bacillus megaterium, Bacillus mycoides, Bacillus stearothermophilus, Bacillus subtilis, Bacillus Thuringiensis, Bacteroides fragilis, Bacteroides gingivalis, Bacteroides melaninogenicus (also known as Prevotella melaninogenica), Bartonella henselae, Bartonella quintana, Bordetella, Bordetella bronchiseptica, Bordetella pertussis, Borrelia afzelii, Borrelia burgdorferi, Borrelia garinii, Borrelia recurrentis, Brucella abortus, Brucella canis, Brucella melitensis, Brucella suis, Burkholderia mallei, Burkholderia pseudomallei, Burkholderia cepacia, Calymmatobacterium granulomatis, Campylobacter, Campylobacter coli, Campylobacter fetus, Campylobacter jejuni, Campylobacter pylori, Chlamydophila pneumoniae (previously called Chlamydia pneumoniae), Chlamydophila psittaci (previously called Chlamydia psittaci), Chlamydia trachomatis, Clostridium botulinum, Clostridium difficile, Clostridium perfringens (previously called Clostridium welchii), Clostridium tetani, Corynebacterium diphtheriae, Corynebacterium fusiforme, Coxiella burnetii, Ehrlichia canis, Ehrlichia chaffeensis, Enterobacter cloacae, Enterococcus avium, Enterococcus durans, Enterococcus faecalis, Enterococcus faecium, Enterococcus galllinarum, Enterococcus maloratus, Escherichia coli, Francisella tularensis, Fusobacterium nucleatum, Gardnerella vaginalis, Haemophilus ducreyi, Haemophilus influenzae, Haemophilus parainfluenzae, Haemophilus pertussis, Haemophilus vaginalis, Helicobacter pylori, Klebsiella pneumoniae, Lactobacillus acidophilus, Lactobacillus bulgaricus, Lactobacillus casei, Lactococcus lactis, Legionella pneumophila, Leptospira interrogans, Leptospira santarosai, Leptospira weilii, Leptospira noguchii, Listeria monocytogenes, Methanobacterium extroquens, Microbacterium multiforme, Micrococcus luteus, Moraxella catarrhalis, Mycobacterium avium, Mycobacterium bovis, Mycobacterium diphtheriae, Mycobacterium intracellulare, Mycobacterium leprae, Mycobacterium lepraemurium, Mycobacterium phlei, Mycobacterium smegmatis, Mycobacterium tuberculosis, Mycobacterium ulcerans, Mycoplasma fermentans, Mycoplasma genitalium, Mycoplasma hominis, Mycoplasma penetrans, Mycoplasma pneumoniae, Neisseria gonorrhoeae, Neisseria meningitidis, Pasteurella multocida, Pasteurella tularensis, Pepto streptococcus, Porphyromonas gingivalis, Prevotella melaninogenica (previously called Bacteroides melaninogenicus), Pseudomonas aeruginosa, Rhizobium radiobacter, Rickettsia prowazekii, Rickettsia psittaci, Rickettsia quintana, Rickettsia rickettsii, Rickettsia trachomae, Rochalimaea henselae, Rochalimaea quintana, Rothia dentocariosa, Salmonella enteritidis, Salmonella typhi, Salmonella typhimurium, Serratia marcescens, Shigella dysenteriae, Shigella sonnei, Spirillum volutans, Streptococcus agalactiae, Staphylococcus aureus, Staphylococcus epidermidis, Staphylococcus saprophyticus, Stenotrophomonas maltophilia, Streptococcus agalactiae, Streptococcus avium, Streptococcus bovis, Streptococcus cricetus, Streptococcus faceium, Streptococcus faecalis, Streptococcus ferus, Streptococcus gallinarum, Streptococcus lactis, Streptococcus mitior, Streptococcus mitis, Streptococcus mutans, Streptococcus oralis, Streptococcus pneumoniae, Streptococcus pyogenes, Streptococcus rattus, Streptococcus salivarius, Streptococcus sanguis, Streptococcus sobrinus, Streptococcus viridans, Treponema pallidum, Treponema denticola, Ureaplasma urealyticum, Vibrio cholerae, Vibrio comma, Vibrio parahaemolyticus, Vibrio vulnificus, Yersinia enterocolitica, Yersinia pestis, Yersinia pseudotuberculosis,
Non-limiting exemplary viruses include the family Retroviridae, such as human deficiency viruses, such as HIV-I (also referred to as HTLV-III), HIV-II, LAC, IDLV-III/LAV, HIV-III or other isolates such as HIV-LP, the family Picornaviridae, such as poliovirus, hepatitis A, enteroviruses, human Coxsackie viruses, rhinoviruses, echoviruses, the family Calciviridae, such as viruses that cause gastroenteritis, the family Togaviridae, such as equine encephalitis viruses and rubella viruses, the family Flaviviridae, such as dengue viruses, encephalitis viruses and yellow fever viruses, the family Coronaviridae, such as coronaviruses, the family Rhabdoviridae, such as vesicular stomata viruses and rabies viruses, the family Filoviridae, such as Ebola viruses, the family Paramyxoviridae, such as parainfluenza viruses, mumps viruses, measles virus and respiratory syncytial virus, the family Orthomyxoviridae, such as influenza viruses, the family Bungaviridae, such as Hataan viruses, bunga viruses, phleoboviruses and Nairo viruses, the family Arena viridae, such as hemorrhagic fever viruses, the family Reoviridae, such as reoviruses, orbiviruses and rotaviruses, the family Bimaviridae, the family Hepadnaviridae, such as hepatitis B virus, the family Parvoviridae, such as parvoviruses, the Papovaviridae, such as papilloma viruses and polyoma viruses, the family Adenoviridae, such as adenoviruses, the family Herpesviridae, such as herpes simplex virus (HSV) I and II, varicella zoster virus and pox viruses, or the family Iridoviridae, such as African swine fever virus). The virus can be an unclassified virus, such as the etiologic agents of Spongiform encephalopathies, the agent of delta hepatitis, the agents of non-A, non-B hepatitis (class 1 enterally transmitted; class 2 parenterally transmitted such as Hepatitis C); Norwalk and related viruses and astroviruses.
Other non-limiting exemplary viruses include Varicella-zoster virus, Epstein-barr virus, Human cytomegalovirus, Human herpesvirus type 8, Human papillomavirus, BK virus, JC virus, Smallpox, Parvovirus B19, poliovirus, yellow fever virus, West Nile virus, TBE virus, Rubella virus, Hepatitis E virus, Influenza virus, Lassa virus, Crimean-Congo hemorrhagic fever virus, Hantaan virus, Marburg virus, Coltivirus, Banna virus, and zika virus.
Non-limiting exemplary fungi include Candida albicans, Aspergillus fumigatus, Aspergillus flavus, Aspergillus clavatus, Cryptococcus neoformans, Cryptococcus laurentii, Cryptococcus albidus, Cryptococcus gattii, Histoplasma capsulatum, Pneumocystis jirovecii, Pneumocystis carinii, and Stachybotrys chartarurn.
Non-limiting exemplary protozoa include Entamoeba histolytica, Entamoeba coli, Entamoeba dispar, Entamoeba moshkovskii, Entamoeba Bangladeshi, Entamoeba hartmanni, Dientamoeba fragilis, Endolimax nana, Lodarnoeba butschlii, Plasmodium malariae, Plasmodium falciparum, Plasmodium vivax, Plasmodium ovale, Naegleria fowleri, Acanthamoeba species, Balamuthia mandrillaris, Sappinia diploidea, Giardia larnblia, Giardia intestinalis, Giardia duodenalis, Toxoplasma gondii, Nippostrongylus brasiliensis, Cryptosporidium parvum, Cryptosporidium hominis, Cryptosporidium cams, Cryptosporidium felis, Cryptosporidium meleagridis, Cryptosporidium muris, Trichomonas vaginalis, Trypanosoma cruzi, Leishmania major, Leishmania tropica, Leishmania barziliensis, Leishmania mexicana, Leishmania guyanesis, Leishmania panamensis, and Trypanosoma brucei.
The computer system for implementing the present invention can take the form of an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, microcode, etc.), or a combination of software and hardware that may all generally be referred to herein as a “circuit,” “module,” or “system.”
The present invention may be a system, a method, and/or a computer program product at any possible technical detail level of integration. The computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.
The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. A non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or raised structures in a groove having instructions recorded thereon, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.
Computer readable program instructions described herein can be downloaded to respective computing/processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network. The network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers. A network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.
Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, configuration data for integrated circuitry, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++, or the like, and procedural programming languages, such as the “C” programming language or similar programming languages. The computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.
Aspects of the present invention are described herein with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer readable program instructions.
These computer readable program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks. These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function/act specified in the flowchart and/or block diagram block or blocks.
The computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.
The flowchart and block diagrams in the Figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s). In some alternative implementations, the functions noted in the blocks may occur out of the order noted in the Figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems that perform the specified functions or acts or carry out combinations of special purpose hardware and computer instructions.
In
Hardware data storage devices 111 may include, but are not limited to, magnetic tape drives, fixed or removable hard disks, optical discs, storage-equipped mobile devices, and solid-state random-access or read-only storage devices. I/O devices may comprise, but are not limited to: input devices 113, such as keyboards, scanners, handheld telecommunications devices, touch-sensitive displays, tablets, biometric readers, joysticks, trackballs, or computer mice; and output devices 115, which may comprise, but are not limited to printers, plotters, tablets, mobile telephones, displays, or sound-producing devices. Data storage devices 111, input devices 113, and output devices 115 may be located either locally or at remote sites from which they are connected to I/O Interface 109 through a network interface.
Processor 103 may also be connected to one or more memory devices 105, which may include, but are not limited to, Dynamic RAM (DRAM), Static RAM (SRAM), Programmable Read-Only Memory (PROM), Field-Programmable Gate Arrays (FPGA), Secure Digital memory cards, SIM cards, or other types of memory devices.
At least one memory device 105 contains stored computer program code 107, which is a computer program that comprises computer-executable instructions. The stored computer program code can include a program for natural-language processing that implements the disclosed methods. The data storage devices 111 may store the computer program code 107. Computer program code 107 stored in the storage devices 111 can be configured to be executed by processor 103 via the memory devices 105. Processor 103 can execute the stored computer program code 107.
Thus the present invention discloses a process for supporting computer infrastructure, integrating, hosting, maintaining, and deploying computer-readable code into the computer system 101, wherein the code in combination with the computer system 101 is capable of performing the analysis of sequence data pertinent to the formation of the self-consistent k-mer database from a reference k-mer database, and generating reports therefrom. The computer system 101 is capable of performing the analysis of sequence data of a sample pertinent to the determination of identifying species using the self-consistent k-mer database.
Any of the components of the present invention could be created, integrated, hosted, maintained, deployed, managed, serviced, supported, etc. by a service provider. Thus, the present invention discloses a process for deploying or integrating computing infrastructure, comprising integrating computer-readable code into the computer system 101, wherein the code in combination with the computer system 101 is capable of performing the analysis of sequence data pertinent to the determination of identifying the viable species of the sample.
One or more data storage units 111 (or one or more additional memory devices not shown in
While it is understood that program code 107 may be deployed by manually loading the program code 107 directly into client, server, and proxy computers (not shown) by loading the program code 107 into a computer-readable storage medium (e.g., computer data storage device 111), program code 107 may also be automatically or semi-automatically deployed into computer system 101 by sending program code 107 to a central server (e.g., computer system 101) or to a group of central servers. Program code 107 may then be downloaded into client computers (not shown) that will execute program code 107.
Alternatively, program code 107 may be sent directly to the client computer via e-mail. Program code 107 may then either be detached to a directory on the client computer or loaded into a directory on the client computer by an e-mail option that selects a program that detaches program code 107 into the directory.
Another alternative is to send program code 107 directly to a directory on the client computer hard drive. If proxy servers are configured, the process selects the proxy server code, determines on which computers the proxy servers' code should be placed, transmits the proxy server code, and then installs the proxy server code on the proxy computer. Program code 107 is then transmitted to the proxy server and stored on the proxy server.
In one embodiment, program code 107 is integrated into a client, server and network environment by providing for program code 107 to coexist with software applications (not shown), operating systems (not shown) and network operating systems software (not shown) and then installing program code 107 on the clients and servers in the environment where program code 107 will function.
The first step of the aforementioned integration of code included in program code 107 is to identify any software, including the network operating system (not shown), on the clients and servers where program code 107 will be deployed, that is required by program code 107 or that works in conjunction with program code 107. This identified software includes the network operating system, where the network operating system comprises software that enhances a basic operating system by adding networking features. Next, the software applications and version numbers are identified and compared to a list of software applications and correct version numbers that have been tested to work with program code 107. A software application that is missing or that does not match a correct version number is upgraded to the correct version.
A program instruction that passes parameters from program code 107 to a software application is checked to ensure that the instruction's parameter list matches a parameter list required by the program code 107. Conversely, a parameter passed by the software application to program code 107 is checked to ensure that the parameter matches a parameter required by program code 107. The client and server operating systems, including the network operating systems, are identified and compared to a list of operating systems, version numbers, and network software programs that have been tested to work with program code 107. An operating system, version number, or network software program that does not match an entry of the list of tested operating systems and version numbers is upgraded to the listed level on the client computers and upgraded to the listed level on the server computers.
After ensuring that the software, where program code 107 is to be deployed, is at a correct version level that has been tested to work with program code 107, the integration is completed by installing program code 107 on the clients and servers.
Embodiments of the present invention may be implemented as a method performed by a processor of a computer system, as a computer program product, as a computer system, or as a processor-performed process or service for supporting computer infrastructure.
The following example illustrates the practice of the invention and is depicted in the flow chart of
The raw nucleotide sequences (reads) were downloaded via NCBI's prefetch and Fastq-dump utilities. Quality control produced minimally processed paired-end files with technical reads (i.e., identical reads) removed and left- and right-clipping performed.
The metagenomic subset of the working reads was input to a taxonomic profiling tool (Kraken), which generated a list of prevalent bacterial species most likely to be in the sample. Using this list, the metagenomic reads were aligned to representative reference genomes of each of the bacterial species. The resulting alignments were stored in sequence alignment map (SAM) files.
The SAM file was then assessed with the iRep software program found in the public Github repository [www.github.com/christophertbrown/iRep], which generated an estimated growth rate for each species of the working sample. iRep implements the method of Brown et al., “Measurement of bacterial replication rates in microbial communities”, Nature Biotechnology 2016, pp 1256-1263). The growth rate estimates were collated into a single data file with the following columns: sample, species, and growth. Table 1 illustrates a small portion of this data file for three bacteria, Akkermansia_muciniphila, Alistipes_onderdonki, and Alistipes putredinis.
The metatranscriptomic data was aligned to annotated reference genomes of the PATRIC database using the same list of prevalent bacterial species as before, which generated a list of expression levels of the genes for each of the identified bacterial species. PATRIC is an all-bacterial bioinformatics database. Wattam et al., “Improvements to PATRIC, the all-bacterial Bioinformatics Database and Analysis Resource Center,” Nucleic Acids Res., 2017, 45(D1):D535-D542.
The relationship between a real-valued phenotype (growth rate from metagenomic data) and the expression levels of a large selection of genes (median coverages from metatranscriptomic data) was evaluated using a differential expression testing framework called limma with the voom extension, developed by Ritchie et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, 2015. This generated a collection of species-specific genes that had non-random linear association with estimated growth rate, which was used to predict growth rate of each of the species. Table 2 illustrates a small portion of this data for some of the genes of the species Akkermansia muciniphila. “Log FC” is the estimate of the log2-fold-change corresponding to the effect or contrast between two or more experimental conditions, “Ave_expr” is the average log2-expression level for that gene across all the arrays and channels in the experiment, “t” is the moderated t-statistic, “P Value” is the raw p-value. “Adj.P.Value” is the p-value adjusted for multiple testing. “B” is the log-odds that the gene is differentially expressed. The probability of a gene being differentially expressed increases with increasing B value.
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
Akkermansia_muciniphila_55290
With a larger training dataset, and more elaborate machine learning algorithms, predictive performance could by further refined, with the ultimate goal of using only metatranscriptomic data to accurately predict the growth rate that previously would need metagenomic data to estimate.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the singular forms “a”, “an” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms “comprises” and/or “comprising,” when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. When a range is used to express a possible value using two numerical limits X and Y (e.g., a concentration of X ppm to Y ppm), unless otherwise stated the value can be X, Y, or any number between X and Y.
The corresponding structures, materials, acts, and equivalents of all means or step plus function elements in the claims below are intended to include any structure, material, or act for performing the function in combination with other claimed elements as specifically claimed. The description of the present invention has been presented for purposes of illustration and description, but is not intended to be exhaustive or limited to the invention in the form disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the invention. The embodiments were chosen and described in order to best explain the principles of the invention and their practical application, and to enable others of ordinary skill in the art to understand the invention.
