Abstract
Nowadays, we are experiencing the big data era with the emerging challenge of single data interpretation. Although the advent of high-throughput technologies as well as chemo- and bio-informatics tools presents pan-omics data as the way forward to precision medicine, personalized health care and tailored-made therapeutics can be only envisaged when interindividual variability in response to/toxicity of xenobiotics can be interpreted and thus, predicted. We know that such variability is the net outcome of genetics (host and microbiota) and environmental factors (diet, lifestyle, polypharmacy, and microbiota) and for this, tremendous efforts have been made to clarify key-molecules from correlation to causality to clinical significance. Herein, we focus on the host–microbiome interplay and its direct and indirect impact on efficacy and toxicity of xenobiotics and we inevitably wonder about the role of viruses, as the least acknowledged ones. We present the emerging discipline of pharmacometabolomics-informed viromics, in which pre-dose metabotypes can assist modeling and prediction of interindividual response to/toxicity of xenobiotics. Such features, either alone or in combination with host genetics, can power biomarker discovery so long as the features are variable among patients, stable enough to be of predictive value, and better than pre-existing tools for predicting therapeutic efficacy/toxicity.
Introduction
Complex interactions between the host immune system and microbiota are of dynamic nature and hence, of fundamental importance when homeostasis is considered, as the host is exposed to trillions of indigenous microorganisms; bacteria, archaea, fungi, and viruses (Hooper et al., 2012; Virgin, 2014) (Figure 1). Metagenomics data suggest that the microbiome of healthy humans includes several viral genes and the intestines and skin of healthy individuals are associated with viruses that replicate in eukaryotic and prokaryotic cells (Reyes et al., 2010). Temperate bacteriophages, which are assembled and propagated stably within the bacterial host chromosome, are considered as the predominant viral group (Reyes et al., 2010; Minot et al., 2013). Plant-infecting viruses, herpesviruses, poxviruses and picornaviruses are less abundant and often derived from diet (Kim et al., 2011). Yet, the great majority of viruses remains unidentified (Minot et al., 2012a,b; Reyes et al., 2012) and much less is known about the role of resident viruses.
FIGURE 1
We believe in a systems-level understanding of the host–microbiome interplay and its impact on the efficacy/toxicity of xenobiotics with an emphasis on resident viruses, as their diversity, abundance and role is still poorly understood. Herein, we wonder about the role of viruses, as the least acknowledged ones and we present the emerging discipline of pharmacometabolomics-informed viromics, in which pre-dose metabotypes (metabolic phenotypes) can assist modeling and prediction of interindividual variability in response to/toxicity of xenobiotics. Such features, either alone or in combination with host genetics, can power biomarker discovery so long as the features are variable among patients, stable enough to be of predictive value, and better than pre-existing tools for predicting therapeutic efficacy/toxicity.
Horizontal Genetic Transfer Results in Evolutionary Changes and Confers Interindividual Variability
Viral–microbial interactions in the human gut seem to be rather different from the predator-prey relationship that is known as “kill the winner” and dominates in other ecosystems (Reyes et al., 2010). Horizontal genetic transfer occurs frequently among gut microbes (Roberts et al., 2008; Tamames and Moya, 2008). Viruses, particularly bacteriophages, are one of the main drivers of the evolutionary change seen in microorganisms through horizontal genetic transfer. The observation that gut microbial metagenomes consist of a large number of phage-related genes implies a viral role in gut homeostasis (Qin et al., 2010). In healthy gut, viruses may elicit innate immune responses. Pattern-recognition receptors can detect viral components (Takeuchi and Akira, 2010) and initiate the crosstalk between resident viruses and the innate immune system. Commensal bacteria-depleted mice have been reported to recover from severe gut inflammation upon administration of lipoteichoic acid or lipopolysaccharide, implying that toll-like receptors have a protective role that extends beyond the recognition of commensal bacteria to other microbes (Rakoff-Nahoum et al., 2004). Resident gut viruses have reduced inflammation via toll-like receptor 3- and 7-mediated interferon-beta production, shedding light for the first time on the interplay of host innate immune system and viruses (Yang et al., 2016). Horizontal genetic transfer confers increased pathogenicity, antibiotic resistance and new metabolic activity (Reyes et al., 2010; Minot et al., 2011). Therefore, bacteriophages are of biomedical importance.
Within-Individual and Inter-individual Viral Evolution
Metagenomic data emphasize that bacterial strains are significantly different among individuals, although their gut typically contains bacteria from only a few phyla (Turnbaugh et al., 2009). Different body sites are also inhabited by different bacterial strains (Peterson et al., 2009; Blaser, 2010; Chen et al., 2010; Qin et al., 2010). This variation may account for their highly variable phage predators, as phages can be highly selective (Rodriguez-Valera et al., 2009). Notably, phage sensitivity is used in the clinic to distinguish bacterial strains (Sell et al., 1983; Mahony et al., 1991). Furthermore, the great variability reported in phage populations among individuals may come from within-individual viral evolution, when long-term viral residents are considered. Rapid within-host viral evolution has been reported, suggesting that multiple new viral species arise in the gut of a typical human over the course of human life (Minot et al., 2013). Notwithstanding, the forces diversifying bacteriophage genomes in human hosts have not been studied in detail.
We strongly feel that changes in abundance and composition of microbiota (dysbiosis) underline interindividual variability. That said, gut microbiota and particularly, resident viruses have not been a focus for the drug metabolism and toxicology communities, despite several early studies showing their importance in xenobiotic biotransformation (Boxenbaum et al., 1979; Illing, 1981). Notwithstanding, re-evaluation and awareness are critical, as microbiota represent a source of physiological variability among individuals and populations that can readily affect the disposition and toxicity of xenobiotics and their metabolites. These effects can be direct or indirect (Nayak and Turnbaugh, 2016). Indirect effects include the metabolic exchange and the co-metabolism and processing of endogenous and dietary substrates (Nicholson and Wilson, 2003). Gut microbiome can alter the expression of host’s pharmacogenes (Björkholm et al., 2009). Microbiome-derived metabolites have been reported to modulate the drug metabolizing systems of the host (Wilson and Nicholson, 2016). Orally administered xenobiotics are exposed to gut microbes prior to their absorption and hence, their bioavailability and half-life is altered (Al-Hilal et al., 2013). The compositional and/or functional alterations in gut microbiota brought about by polypharmacy or the administration of antibiotics, probiotics, or prebiotics greatly increase interindividual variability in response to xenobiotics (Li and Jia, 2013). Notably, the role of resident viruses is still poorly understood.
Challenges and Opportunities
If microbiome in total has escaped our full attention so far, resident viruses in humans and animals are poorly acknowledged today. Even though viruses are extremely important to the ecology (Roossinck, 2015; Puxty et al., 2016), the discipline of viromics (or viral metagenomics, as it is used interchangeably) can be currently considered as the most elusive of the -omics fields.
“Which Entities are Present?”
The human body is a host for a complex living microbial community consisting of bacteria, archaea, fungi, protozoa and viruses, which together constitute the human microbiota. The total genome of microorganisms is also the human metagenome and hence, metagenomic studies aim to address the question of “which entities are present?” via high-throughput sequencing or microarrays. Technical challenges are still profound, when seeking for an answer to that question (Figure 1).
Viruses outnumber microbial cells 10:1 in most environments, yet, viral DNA only represents 0.1% of the total DNA in a microbial community, such as the human gut (Qin et al., 2010). If a deep sequence coverage of the human viruses is desired, viral particles should be isolated (Reyes et al., 2010; Hurwitz et al., 2013). As the amount of the nucleic acids extracted from purified viral particles is often below the required threshold for sequencing, several amplification methods have been developed; random amplified shotgun library (Rohwer et al., 2001), linker-amplified shotgun library (Breitbart et al., 2003) and multiple displacement amplification (Hutchison et al., 2005) are a few examples. The latter method is greatly advantageous, as it allows the amplification of complete viral genomes, although recent reports have implied that critical biases and contamination is an issue (Zhang et al., 2006; Yilmaz et al., 2010).
Today, the majority of the viral sequences are novel or enriched in regions of low-complexity repeats. Thus, sequencing technologies that prioritize long-read lengths over those of short-read lengths are preferred (Wommack et al., 2008). However, sequencing technologies of long-read lengths, such as 454/Roche pyrosequencing, are about to be discontinued (Bikel et al., 2015). Bioinformatics tools that have been developed to analyze viruses from short sequence reads have accepted this challenge (Liu et al., 2007). Moreover, the majority (usually, 60–99%) of the viral sequences studied so far have no significant similarity to other sequences in databases or have higher homology to prokaryotic or eukaryotic genes (Breitbart et al., 2002, 2003; Blomström et al., 2010). A crucial step in viromics is the filtering of bad quality sequences or the decontamination of 16S rRNA, 18S rRNA and human sequences by mapping. Mapping algorithms or tools, such as ViroBLAST (Deng et al., 2007), BLASTX (Segata et al., 2012) or USEARCH (Edgar, 2010) are employed to obtain the taxonomic composition of a viral community. Taxonomic and functional assignments are crucial for the viral community profile to be created, since it reflects sample diversity. Today, the number of the deposited genomes in databases is far less than the expected number of virotypes (Rohwer and Thurber, 2009), while most of the new sequences are poorly annotated (Brister et al., 2014; Martınez et al., 2014). The percentage of sequence reads with similarity to known viral sequences depends on the database used and on how well sequences have been filtered. This percentage is less than 0.01% (Yang et al., 2011; Bikel et al., 2015). MetaVir, VIROME, and iVirus are few publicly available databases (Sullivan, 2015).
To assist viromics data analysis similarity-independent methods have been also developed. PHACCS assesses the biodiversity of uncultured viral communities and quantifies virotypes (Breitbart et al., 2002; Angly et al., 2005) and CRASS allows the simultaneous cross-assembly of all the samples in a dataset to identify shared viruses (Dutilh et al., 2012). MaxiPhi uses pairwise assemblies from pooled viromes (Angly et al., 2005). To tackle chimeras, the overlap–layout–consensus algorithms have achieved the assembly of viral genomes. Newbler has been extensively used in viral and bacterial shotgun metagenomic projects (Qin et al., 2010; Reyes et al., 2010; Bikel et al., 2015). Yet, it remains to be determined, if Newbler will be discontinued with the 454/Roche in 2016 (Bikel et al., 2015). Minimo is designed for the assembly of small datasets (Treangen et al., 2011) and VICUNA is an assembler that is specialized in de novo assembly of data from heterogeneous viral populations (Yang et al., 2012). MetaVelvet (Namiki et al., 2012) and other de Brujin graph assemblers are an alternative to the overlap–layout–consensus assemblers and have also been used on the assembly of viral genomes (Hurwitz and Sullivan, 2013). Sequencing of total RNA viruses has been proven impractical (Robertson et al., 2010).
“Which Entities are Active?”
One might have already thought that this question is even more difficult to address. We share an optimistic view on the basis of the metatranscriptomic, metaproteomic and metabonomic studies being reported (Lim et al., 2013; Aguiar-Pulido et al., 2016; Yu et al., 2016). The NIH-sponsored Human Microbiome Project has been established to comprehensively characterize human microbiota from multiple body sites and analyze their impact in human health and disease (Peterson et al., 2009; Proctor, 2016). The Vaginal Human Microbiome Project has validated a protocol that achieves species-level classification of V1–V3 16S rRNA sequence from the vaginal microbiome (Fettweis et al., 2012). Culturomics, coupled to matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry, has been also introduced as complementary to metagenomics to study complex microbial ecosystems, such as the human gut. Lagier et al. (2012) provided a proof-of-concept applying 212 different culture conditions and successfully culturing 340 different bacterial species, five fungi and for the first time, Senegalvirus (Greub, 2012; Lagier et al., 2012). Dubourg et al. (2014) coupled culturomics with pyrosequencing to address the effect of antibiotics in gut microbiota diversity (Dubourg et al., 2014). Yet, resident viruses and their interplay with the host (human or other microbes) are poorly studied.
“Which Entities are Predictive of Response to/Toxicity of Xenobiotics?”
Every human being is unique and microbiota that colonize the human body shape individuality to a great extent. Interindividual variability is a prerequisite to tailored-made therapeutics and personalized health care. A systems-level understanding is required to delineate the host–microbiota interplay. The advent of high-throughput technologies and bio- and chemo-informatics tools presents pan-omics approaches as the way forward. Can we predict xenobiotics efficacy/toxicity by pre-dose profiling?
Traditionally, the host has been our focus, embracing the potential of pharmacogenomics and pharmacometabolomics. Pharmacogenomics has pioneered the prediction of the outcome of a xenobiotic intervention in an individual based on an analysis of that individual’s genetic profile (Smalley and Sondak, 2010; Ritchie, 2012; Hershfield et al., 2013; Everett, 2015). Today, FDA-approved drug labeling may contain information on pharmacogenomic biomarkers. The U.S. National Human Genome Research Institute and the National Academy of Medicine recently brought together 25 innovative genomic medicine programs to coalesce innovative genomic medicine programs around concrete and at the same time, compelling signature projects accelerating the responsible implementation of genomic medicine in efforts to improve clinical care worldwide (Manolio et al., 2015). Pharmacometagenomics are anticipated to perform an analogous task to that of pharmacogenomics focusing on human metagenome (the total genome of microorganisms). Indeed, a more mechanistic understanding of which microbes and genes contribute to xenobiotics efficacy/toxicity may enable prediction of which patients will derive greatest benefit from a therapeutic intervention (Haiser et al., 2013; Nayak and Turnbaugh, 2016). Pharmacometabolomics is predictive of the outcome of a xenobiotic intervention in an individual based on a mathematical model of pre-intervention metabolite signatures (Nicholson et al., 1999; Clayton et al., 2006). Pharmacometabolomics is based on metabotypes, as they are the net result of genetic, physiological, chemical, and environmental influences (Holmes et al., 2008; Everett, 2015). Recently, a new concept has arisen, namely “pharmacometabolomics-aided pharmacogenomics,” to reinforce the identification and validation of clinically relevant associations (Ji et al., 2011; Suhre et al., 2011; Abo et al., 2012). We have recently proposed that pharmacometabolomics-aided pharmacogenomics may have an even greater impact if coupled to information technologies to facilitate data analysis and sense- and decision-making on the basis of a synergy between artificial and human intelligence (Katsila et al., 2016).
If we want to turn information growth into knowledge growth and better informed decisions, the implementation of new working practices is more than necessary. A paradigm may well be pharmacometabolomics-informed viromics. In an analogy to pharmacometabolomics-aided pharmacogenomics, pharmacometabolomics-informed viromics benefits from pre-dose metabotypes that can assist modeling and prediction of interindividual response to and/or toxicity of xenobiotics, this time not focusing on the host, but on the resident viruses and their interactions. Although bacteria have been directly associated with dysbiosis and interindividual variability upon xenobiotics administration, the role of the virome in the microbial community should be further explored. So far, scientists have been focusing on the “which entities are present?” and “which entities are active?” questions with the aim to address the “which entities are predictive of response to/toxicity of xenobiotics?” challenge. Instead, we suggest a global focus on metabolites as the biochemical end points of complex interactions between the host and its environment (including the gut microbiome). Thus, a global readout will be revealed mapping microbial, viral and mammalian interindividual variability upon xenobiotics administration. Our next step is to integrate phenotypic and genotyping approaches with information technologies.
Instead of any single omics approach, we propose an integrated (transomic) analysis that is anticipated to provide more insights into the emergence of the phenotypes in question (response to/toxicity of xenobiotics) than any layer can by itself, highlighting the complementarity of a multilayered strategy. In this context, we have two layers of systems-scale molecular measurements; the pharmacometabolome (layer 1) and the virome (layer 2; Figure 2). Layer 1 characterization includes sample acquisition and preparation, analysis (NMR or mass spectrometry technologies), data processing and data analysis (targeted and untargeted). Particularly, untargeted analysis is of great benefit as a tool to shape hypothesis; multiple analytes are quantified simultaneously and pharmacometabolomic modeling is not limited by prior understanding or hypotheses. Such metabotype-based findings may be patient and/or xenobiotic profiling. MALDI-TOF mass spectrometry has recently enabled spatial information. Dorrestein’s methodologies aim to scan microbial communities, capture their exact image and identify what microorganisms do in their complex communities by identifying and locating the metabolites produced (Tullis, 2016). Pirhaji et al. (2016) reported PIUMet, a network-based approach, prize-collecting Steiner forest algorithm for integrative analysis of untargeted metabolomics (Pirhaji et al., 2016). PIUMet infers molecular pathways and components via integrative analysis of metabolite features, without requiring their identification. If a hypothesis is already in place, targeted pharmacometabolomics is advantageous (Williams et al., 2016). Then, viromics (layer 2) is performed as of today. VirCapSeq-VERT is ideally suited for the analysis of virome composition and dynamics, as its highly multiplexed nature allows the simultaneous identification and comprehensive genetic characterization of all known vertebrate viruses and their genetic variants as well as the novel ones (Briese et al., 2015). Similarly, quantitative temporal viromics are applied to a wide range of viruses (Weekes et al., 2014). At this point, transomic data sets are generated that consist of pharmacometabolomics (layer 1) and viromics (layer 2) data. We propose in-depth data mining, analysis and argumentation, according to which information technologies provide the means for filtering and systems-level dynamic parameters from fewer samples across broad molecular interaction networks. Data mining, analysis, collaboration and decision-making in such diverse data-intensive and cognitively complex settings is performed via the Dicode approach, supporting artificial and human intelligence (the Dicode platform and services provide a remedy to the information and cognitive overload as users can customize the Dicode workbench via a proper assembly of tools that suit their needs and properly structured data lead to more informed decisions). The envisioned architecture combines batch and stream processing (Karacapilidis, 2013; Tsiliki et al., 2014). Current limitations, such as lack of annotation, lack of a conserved region in the virome and/or “non-cultivable” viral entities are by-passed. Viral–microbiota–host–xenobiotics relationships are revealed and modeled following a multilayered approach.
FIGURE 2
Conclusion and Future Perspectives
Microbiota in humans is a collection of microscopic organisms that inhabit the body and contains representatives from all the domains of life: archaea, bacteria, eukarya, viruses. If interindividual variability in response to and/or toxicity of xenobiotics is to be acknowledged, viral–microbial–host–xenobiotics dynamics need to be further clarified. Pharmacometabolomics-informed viromics coupled to information technologies highlight the complementarity of a multilayered approach to turn information growth into knowledge growth.
Statements
Author contributions
Conception and design: AB, TK, and GP. Writing, review, and/or revision of the manuscript: AB, TK, and GP.
Funding
This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 668353.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AboR.HebbringS.JiY.ZhuH.ZengZ.-B.BatzlerA.et al (2012). Merging pharmacometabolomics with pharmacogenomics using’1000 Genomes’ single-nucleotide polymorphism imputation: selective serotonin reuptake inhibitor response pharmacogenomics.Pharmacogenet. Genomics22247–253. 10.1097/FPC.0b013e32835001c9
2
Aguiar-PulidoV.HuangW.Suarez-UlloaV.CickovskiT.MatheeK.NarasimhanG. (2016). Metagenomics, metatranscriptomics, and metabolomics approaches for microbiome analysis.Evol. Bioinform. Online125–16.
3
Al-HilalT. A.AlamF.ByunY. (2013). Oral drug delivery systems using chemical conjugates or physical complexes.Adv. Drug Deliv. Rev.65845–864. 10.1016/j.addr.2012.11.002
4
AnglyF.Rodriguez-BritoB.BangorD.McNairnieP.BreitbartM.SalamonP.et al (2005). PHACCS, an online tool for estimating the structure and diversity of uncultured viral communities using metagenomic information.BMC Bioinformatics6:41. 10.1186/1471-2105-6-41
5
BikelS.Valdez-LaraA.Cornejo-GranadosF.RicoK.Canizales-QuinterosS.SoberónX.et al (2015). Combining metagenomics, metatranscriptomics and viromics to explore novel microbial interactions: towards a systems-level understanding of human microbiome.Comput. Struct. Biotechnol. J.13390–401. 10.1016/j.csbj.2015.06.001
6
BjörkholmB.BokC. M.LundinA.RafterJ.HibberdM. L.PetterssonS. (2009). Intestinal microbiota regulate xenobiotic metabolism in the liver.PLoS ONE4:e6958. 10.1371/journal.pone.0006958
7
BlaserM. J. (2010). Harnessing the power of the human microbiome.Proc. Natl. Acad. Sci. U.S.A.1076125–6126. 10.1073/pnas.1002112107
8
BlomströmA.-L.WidénF.HammerA.-S.BelákS.BergM. (2010). Detection of a novel astrovirus in brain tissue of mink suffering from shaking mink syndrome by use of viral metagenomics.J. Clin. Microbiol.484392–4396. 10.1128/JCM.01040-10
9
BoxenbaumH.BekerskyI.JackM.KaplantS. (1979). Influence of gut microflora on bioavailability.Drug Metab. Rev.9259–279. 10.3109/03602537908993894
10
BreitbartM.HewsonI.FeltsB.MahaffyJ. M.NultonJ.SalamonP.et al (2003). Metagenomic analyses of an uncultured viral community from human feces.J. Bacteriol.1856220–6223. 10.1128/JB.185.20.6220-6223.2003
11
BreitbartM.SalamonP.AndresenB.MahaffyJ. M.SegallA. M.MeadD.et al (2002). Genomic analysis of uncultured marine viral communities.Proc. Natl. Acad. Sci. U.S.A.9914250–14255. 10.1073/pnas.202488399
12
BrieseT.KapoorA.MishraN.JainK.KumarA.JabadoO. J.et al (2015). Virome capture sequencing enables sensitive viral diagnosis and comprehensive virome analysis.MBio6:e1491-15. 10.1128/mBio.01491-15
13
BristerJ. R.Ako-AdjeiD.BaoY.BlinkovaO. (2014). NCBI viral genomes resource.Nucleic Acids Res.43D571–D577. 10.1093/nar/gku1207
14
ChenT.YuW.-H.IzardJ.BaranovaO. V.LakshmananA.DewhirstF. E. (2010). The Human Oral Microbiome Database: a web accessible resource for investigating oral microbe taxonomic and genomic information.Database2010:baq013. 10.1093/database/baq013
15
ClaytonT. A.LindonJ. C.CloarecO.AnttiH.CharuelC.HantonG.et al (2006). Pharmaco-metabonomic phenotyping and personalized drug treatment.Nature4401073–1077. 10.1038/nature04648
16
DengW.NickleD. C.LearnG. H.MaustB.MullinsJ. I. (2007). ViroBLAST: a stand-alone BLAST web server for flexible queries of multiple databases and user’s datasets.Bioinformatics232334–2336. 10.1093/bioinformatics/btm331
17
DubourgG.LagierJ. C.RobertC.ArmougomF.HugonP.MetidjiS.et al (2014). Culturomics and pyrosequencing evidence of the reduction in gut microbiota diversity in patients with broad-spectrum antibiotics.Int. J. Antimicrob. Agents44117–124. 10.1016/j.ijantimicag.2014.04.020
18
DutilhB. E.SchmiederR.NultonJ.FeltsB.SalamonP.EdwardsR. A.et al (2012). Reference-independent comparative metagenomics using cross-assembly: crAss.Bioinformatics283225–3231. 10.1093/bioinformatics/bts613
19
EdgarR. C. (2010). Search and clustering orders of magnitude faster than BLAST.Bioinformatics262460–2461. 10.1093/bioinformatics/btq461
20
EverettJ. R. (2015). Pharmacometabonomics in humans: a new tool for personalized medicine.Pharmacogenomics16737–754. 10.2217/pgs.15.20
21
FettweisJ. M.SerranoM. G.ShethN. U.MayerC. M.GlascockA. L.BrooksJ. P.et al (2012). Species-level classification of the vaginal microbiome.BMC Genomics13:S17. 10.1186/1471-2164-13-S8-S17
22
GreubG. (2012). Culturomics: a new approach to study the human microbiome.Clin. Microbiol. Infect.181157–1159. 10.1111/1469-0691.12032
23
HaiserH. J.GootenbergD. B.ChatmanK.SirasaniG.BalskusE. P.TurnbaughP. J. (2013). Predicting and manipulating cardiac drug inactivation by the human gut bacterium Eggerthella lenta.Science341295–298. 10.1126/science.1235872
24
HershfieldM.CallaghanJ. T.TassaneeyakulW.MushirodaT.ThornC.KleinT.et al (2013). Clinical pharmacogenetics implementation consortium guidelines for human leukocyte Antigen-B genotype and allopurinol dosing.Clin. Pharmacol. Ther.93153–158. 10.1038/clpt.2012.209
25
HolmesE.WilsonI. D.NicholsonJ. K. (2008). Metabolic phenotyping in health and disease.Cell134714–717. 10.1016/j.cell.2008.08.026
26
HooperL. V.LittmanD. R.MacphersonA. J. (2012). Interactions between the microbiota and the immune system.Science3361268–1273. 10.1126/science.1223490
27
HurwitzB. L.DengL.PoulosB. T.SullivanM. B. (2013). Evaluation of methods to concentrate and purify ocean virus communities through comparative, replicated metagenomics.Environ. Microbiol.151428–1440. 10.1111/j.1462-2920.2012.02836.x
28
HurwitzB. L.SullivanM. B. (2013). The Pacific Ocean Virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology.PLoS ONE8:e57355. 10.1371/journal.pone.0057355
29
HutchisonC. A.SmithH. O.PfannkochC.VenterJ. C. (2005). Cell-free cloning using φ29 DNA polymerase.Proc. Natl. Acad. Sci. U.S.A.10217332–17336. 10.1073/pnas.0508809102
30
IllingH. (1981). Techniques for microfloral and associated metabolic studies in relation to the absorption and enterohepatic circulation of drugs.Xenobiotica11815–830. 10.3109/00498258109045319
31
JiY.HebbringS.ZhuH.JenkinsG. D.BiernackaJ.SnyderK.et al (2011). Glycine and a glycine dehydrogenase (GLDC) SNP as citalopram/escitalopram response biomarkers in depression: pharmacometabolomics-informed pharmacogenomics.Clin. Pharmacol. Ther.8997–104. 10.1038/clpt.2010.250
32
KaracapilidisN. I. (2013). Mastering Data-Intensive Collaboration and Decision Making. KDIR/KMIS.Berlin: Springer.
33
KatsilaT.KonstantinouE.LavdaI.MalakisH.PapantoniI.SkondraL.et al (2016). Pharmacometabolomics-aided pharmacogenomics in autoimmune disease.EBioMedicine540–45. 10.1016/j.ebiom.2016.02.001
34
KimM.-S.ParkE.-J.RohS. W.BaeJ.-W. (2011). Diversity and abundance of single-stranded DNA viruses in human feces.Appl. Environ. Microbiol.778062–8070. 10.1128/AEM.06331-11
35
LagierJ. C.ArmougomF.MillionM.HugonP.PagnierI.RobertC.et al (2012). Microbial culturomics: paradigm shift in the human gut microbiome study.Clin. Microbiol. Infect.181185–1193. 10.1111/1469-0691.12023
36
LiH.JiaW. (2013). Cometabolism of microbes and host: implications for drug metabolism and drug-induced toxicity.Clin. Pharmacol. Ther.94574–581. 10.1038/clpt.2013.157
37
LimY. W.SchmiederR.HaynesM.WillnerD.FurlanM.YouleM.et al (2013). Metagenomics and metatranscriptomics: windows on CF-associated viral and microbial communities.J. Cyst. Fibros.12154–164. 10.1016/j.jcf.2012.07.009
38
LiuZ.LozuponeC.HamadyM.BushmanF. D.KnightR. (2007). Short pyrosequencing reads suffice for accurate microbial community analysis.Nucleic Acids Res.35:e120. 10.1093/nar/gkm541
39
MahonyD.ClowJ.AtkinsonL.VakhariaN.SchlechW. (1991). Development and application of a multiple typing system for Clostridium difficile.Appl. Environ. Microbiol.571873–1879.
40
ManolioT. A.AbramowiczM.Al-MullaF.AndersonW.BallingR.BergerA. C.et al (2015). Global implementation of genomic medicine: we are not alone.Sci. Transl. Med.7:290s13. 10.1126/scitranslmed.aab0194
41
MartınezJ. M.SwanB. K.WilsonW. H. (2014). Marine viruses, a genetic reservoir revealed by targeted viromics.ISME J.81079–1088. 10.1038/ismej.2013.214
42
MinotS.BrysonA.ChehoudC.WuG. D.LewisJ. D.BushmanF. D. (2013). Rapid evolution of the human gut virome.Proc. Natl. Acad. Sci. U.S.A.11012450–12455. 10.1073/pnas.1300833110
43
MinotS.GrunbergS.WuG. D.LewisJ. D.BushmanF. D. (2012a). Hypervariable loci in the human gut virome.Proc. Natl. Acad. Sci. U.S.A.1093962–3966. 10.1073/pnas.1119061109
44
MinotS.WuG. D.LewisJ. D.BushmanF. D. (2012b). Conservation of gene cassettes among diverse viruses of the human gut.PLoS ONE7:e42342. 10.1371/journal.pone.0042342
45
MinotS.SinhaR.ChenJ.LiH.KeilbaughS. A.WuG. D.et al (2011). The human gut virome: inter-individual variation and dynamic response to diet.Genome Res.211616–1625. 10.1101/gr.122705.111
46
NamikiT.HachiyaT.TanakaH.SakakibaraY. (2012). MetaVelvet: an extension of Velvet assembler to de novo metagenome assembly from short sequence reads.Nucleic Acids Res.40:e155. 10.1093/nar/gks678
47
NayakR. R.TurnbaughP. J. (2016). Mirror, mirror on the wall: which microbiomes will help heal them all?BMC Med.14:72. 10.1186/s12916-016-0622-6
48
NicholsonJ. K.LindonJ. C.HolmesE. (1999). ’Metabonomics’: understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data.Xenobiotica291181–1189. 10.1080/004982599238047
49
NicholsonJ. K.WilsonI. D. (2003). Understanding’global’systems biology: metabonomics and the continuum of metabolism.Nat. Rev. Drug Discov.2668–676. 10.1038/nrd1157
50
PetersonJ.GargesS.GiovanniM.McInnesP.WangL.SchlossJ. A.et al (2009). The NIH human microbiome project.Genome Res.192317–2323. 10.1101/gr.096651.109
51
PirhajiL.MilaniP.LeidlM.CurranT.Avila-PachecoJ.ClishC. B.et al (2016). Revealing disease-associated pathways by network integration of untargeted metabolomics.Nat. Methods13770–776. 10.1038/nmeth.3940
52
ProctorL. M. (2016). The National Institutes of Health Human Microbiome Project.Semin. Fetal Neonatal Med.10.1016/j.siny.2016.05.002 [Epub ahead of print].
53
PuxtyR. J.MillardA. D.EvansD. J.ScanlanD. J. (2016). Viruses inhibit CO2 fixation in the most abundant phototrophs on earth.Curr. Biol.261585–1589. 10.1016/j.cub.2016.04.036
54
QinJ.LiR.RaesJ.ArumugamM.BurgdorfK. S.ManichanhC.et al (2010). A human gut microbial gene catalogue established by metagenomic sequencing.Nature46459–65. 10.1038/nature08821
55
Rakoff-NahoumS.PaglinoJ.Eslami-VarzanehF.EdbergS.MedzhitovR. (2004). Recognition of commensal microflora by toll-like receptors is required for intestinal homeostasis.Cell118229–241. 10.1016/j.cell.2004.07.002
56
ReyesA.HaynesM.HansonN.AnglyF. E.HeathA. C.RohwerF.et al (2010). Viruses in the faecal microbiota of monozygotic twins and their mothers.Nature466334–338. 10.1038/nature09199
57
ReyesA.SemenkovichN. P.WhitesonK.RohwerF.GordonJ. I. (2012). Going viral: next-generation sequencing applied to phage populations in the human gut.Nat. Rev. Microbiol.10607–617. 10.1038/nrmicro2853
58
RitchieM. D. (2012). The success of pharmacogenomics in moving genetic association studies from bench to bedside: study design and implementation of precision medicine in the post-GWAS era.Hum. Genet.1311615–1626. 10.1007/s00439-012-1221-z
59
RobertsA. P.ChandlerM.CourvalinP.GuédonG.MullanyP.PembrokeT.et al (2008). Revised nomenclature for transposable genetic elements.Plasmid60167–173. 10.1016/j.plasmid.2008.08.001
60
RobertsonG.ScheinJ.ChiuR.CorbettR.FieldM.JackmanS. D.et al (2010). De novo assembly and analysis of RNA-seq data.Nat. Methods7909–912. 10.1038/nmeth.1517
61
Rodriguez-ValeraF.Martin-CuadradoA.-B.Rodriguez-BritoB.PašiæL.ThingstadT. F.RohwerF.et al (2009). Explaining microbial population genomics through phage predation.Nat. Rev. Microbiol.7828–836. 10.1038/nrmicro2235
62
RohwerF.SeguritanV.ChoiD.SegallA.AzamF. (2001). Production of shotgun libraries using random amplification.Biotechniques31108–119.
63
RohwerF.ThurberR. V. (2009). Viruses manipulate the marine environment.Nature459207–212. 10.1038/nature08060
64
RoossinckM. J. (2015). Metagenomics of plant and fungal viruses reveals an abundance of persistent lifestyles.Front. Microbiol.5:767. 10.3389/fmicb.2014.00767
65
SegataN.WaldronL.BallariniA.NarasimhanV.JoussonO.HuttenhowerC. (2012). Metagenomic microbial community profiling using unique clade-specific marker genes.Nat. Methods9811–814. 10.1038/nmeth.2066
66
SellT. L.SchabergD. R.FeketyF. R. (1983). Bacteriophage and bacteriocin typing scheme for Clostridium difficile.J. Clin. Microbiol.171148–1152.
67
SmalleyK. S.SondakV. K. (2010). Melanoma—an unlikely poster child for personalized cancer therapy.N. Engl. J. Med.363876–878. 10.1056/NEJMe1005370
68
SuhreK.ShinS.-Y.PetersenA.-K.MohneyR. P.MeredithD.WägeleB.et al (2011). Human metabolic individuality in biomedical and pharmaceutical research.Nature47754–60. 10.1038/nature10354
69
SullivanM. B. (2015). Viromes, not gene markers, for studying double-stranded DNA virus communities.J. Virol.892459–2461. 10.1128/JVI.03289-14
70
TakeuchiO.AkiraS. (2010). Pattern recognition receptors and inflammation.Cell140805–820. 10.1016/j.cell.2010.01.022
71
TamamesJ.MoyaA. (2008). Estimating the extent of horizontal gene transfer in metagenomic sequences.BMC Genomics9:136. 10.1186/1471-2164-9-136
72
TreangenT. J.SommerD. D.AnglyF. E.KorenS.PopM. (2011). Next generation sequence assembly with AMOS.Curr. Protoc. BioinformaticsChapter 11:Unit11.8.
73
TsilikiG.KaracapilidisN.ChristodoulouS.TzagarakisM. (2014). Collaborative mining and interpretation of large-scale data for biomedical research insights.PLoS ONE9:e108600. 10.1371/journal.pone.0108600
74
TullisP. (2016). The man who can map the chemicals all over your body.Nature534170–172. 10.1038/534170a
75
TurnbaughP. J.HamadyM.YatsunenkoT.CantarelB. L.DuncanA.LeyR. E.et al (2009). A core gut microbiome in obese and lean twins.Nature457480–484. 10.1038/nature07540
76
VirginH. W. (2014). The virome in mammalian physiology and disease.Cell157142–150. 10.1016/j.cell.2014.02.032
77
WeekesM. P.TomasecP.HuttlinE. L.FieldingC. A.NusinowD.StantonR. J.et al (2014). Quantitative temporal viromics: an approach to investigate host-pathogen interaction.Cell1571460–1472. 10.1016/j.cell.2014.04.028
78
WilliamsE. G.WuY.JhaP.DubuisS.BlattmannP.ArgmannC. A.et al (2016). Systems proteomics of liver mitochondria function.Science352:aad0189. 10.1126/science.aad0189
79
WilsonI. D.NicholsonJ. K. (2016). Gut Microbiome Interactions with Drug Metabolism, Efficacy and Toxicity.Transl. Res.10.1016/j.trsl.2016.08.002 [Epub ahead of print].
80
WommackK. E.BhavsarJ.RavelJ. (2008). Metagenomics: read length matters.Appl. Environ. Microbiol.741453–1463. 10.1128/AEM.02181-07
81
YangJ.YangF.RenL.XiongZ.WuZ.DongJ.et al (2011). Unbiased parallel detection of viral pathogens in clinical samples by use of a metagenomic approach.J. Clin. Microbiol.493463–3469. 10.1128/JCM.00273-11
82
YangJ.-Y.KimM.-S.KimE.CheonJ. H.LeeY.-S.KimY.et al (2016). Enteric viruses ameliorate gut inflammation via Toll-like Receptor 3 and Toll-like Receptor 7-Mediated Interferon-β Production.Immunity44889–900. 10.1016/j.immuni.2016.03.009
83
YangX.CharleboisP.GnerreS.CooleM. G.LennonN. J.LevinJ. Z.et al (2012). De novo assembly of highly diverse viral populations.BMC Genomics13:475. 10.1186/1471-2164-13-475
84
YilmazS.AllgaierM.HugenholtzP. (2010). Multiple displacement amplification compromises quantitative analysis of metagenomes.Nat. Methods7943–944. 10.1038/nmeth1210-943
85
YuJ.MarshS.HuJ.FengW.WuC. (2016). Gut microbiota and metagenomic advancement in digestive disease.Gastroenterol. Res. Pract.2016:4703406. 10.1155/2016/4703406
86
ZhangK.MartinyA. C.ReppasN. B.BarryK. W.MalekJ.ChisholmS. W.et al (2006). Sequencing genomes from single cells by polymerase cloning.Nat. Biotechnol.24680–686. 10.1038/nbt1214
Summary
Keywords
pharmacometagenomics, pharmacometabolomics, pharmacogenomics, precision medicine, viromics, metagenomics
Citation
Balasopoulou A, Patrinos GP and Katsila T (2016) Pharmacometabolomics Informs Viromics toward Precision Medicine. Front. Pharmacol. 7:411. doi: 10.3389/fphar.2016.00411
Received
30 July 2016
Accepted
17 October 2016
Published
27 October 2016
Volume
7 - 2016
Edited by
Daniel Rotroff, North Carolina State University, USA
Reviewed by
John Frederick Pearson, University of Otago, New Zealand; Mohamed H. Shahin, University of Florida, USA
Updates
Copyright
© 2016 Balasopoulou, Patrinos and Katsila.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Theodora Katsila, thkatsila@upatras.gr
This article was submitted to Pharmacogenetics and Pharmacogenomics, a section of the journal Frontiers in Pharmacology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.