Computational tools for viral metagenomics and their application in clinical research

AbstractThere are 100 times more virions than eukaryotic cells in a healthy human body. The characterization of human-associated viral communities in a non-pathological state and the detection of viral pathogens in cases of infection are essential for medical care and epidemic surveillance. Viral metagenomics, the sequenced-based analysis of the complete collection of viral genomes directly isolated from an organism or an ecosystem, bypasses the “single-organism-level” point of view of clinical diagnostics and thus the need to isolate and culture the targeted organism. The first part of this review is dedicated to a presentation of past research in viral metagenomics with an emphasis on human-associated viral communities (eukaryotic viruses and bacteriophages). In the second part, we review more precisely the computational challenges posed by the analysis of viral metagenomes, and we illustrate the problem of sequences that do not have homologs in public databases and the possible approaches to characterize them

Tropheryma whipplei, the Whipple's disease bacillus, induces macrophage apoptosis through the extrinsic pathway

Tropheryma whipplei, the etiological agent of Whipple's disease, is an intracellular bacterium that infects macrophages. We previously showed that infection of macrophages results in M2 polarization associated with induction of apoptosis and interleukin (IL)-16 secretion. In patients with Whipple's disease, circulating levels of apoptotic markers and IL-16 are increased and correlate with the activity of the disease. To gain insight into the understanding of the pathophysiology of this rare disease, we examined the molecular pathways involved in T. whipplei-induced apoptosis of human macrophages. Our data showed that apoptosis induction depended on bacterial viability and inhibition of bacterial protein synthesis reduced the apoptotic program elicited by T. whipplei. Induction of apoptosis was also associated with a massive degradation of both pro- and anti-apoptotic mediators. Caspase-specific inhibition experiments revealed that initiator caspases 8 and 10 were required for apoptosis, in contrast to caspases 2 and 9, in spite of cytochrome-c release from mitochondria. Finally, the effector caspases 3 and 6 were mandatory for apoptosis induction. Collectively, these data suggest that T. whipplei induces apoptosis through the extrinsic pathway and that, beside M2 polarization of macrophages, apoptosis induction contributes to bacterial replication and represents a virulence trait of this intracellular pathogen

Expanding networks of RNA virus evolution

In a recent BMC Evolutionary Biology article, Huiquan Liu and colleagues report two new genomes of double-stranded RNA (dsRNA) viruses from fungi and use these as a springboard to perform an extensive phylogenomic analysis of dsRNA viruses. The results support the old scenario of polyphyletic origin of dsRNA viruses from different groups of positive-strand RNA viruses and additionally reveal extensive horizontal gene transfer between diverse viruses consistent with the network-like rather than tree-like mode of viral evolution. Together with the unexpected discoveries of the first putative archaeal RNA virus and a RNA-DNA virus hybrid, this work shows that RNA viral genomics has major surprises to deliver

Type I Interferon Induction Is Detrimental during Infection with the Whipple's Disease Bacterium, Tropheryma whipplei

Macrophages are the first line of defense against pathogens. Upon infection macrophages usually produce high levels of proinflammatory mediators. However, macrophages can undergo an alternate polarization leading to a permissive state. In assessing global macrophage responses to the bacterial agent of Whipple's disease, Tropheryma whipplei, we found that T. whipplei induced M2 macrophage polarization which was compatible with bacterial replication. Surprisingly, this M2 polarization of infected macrophages was associated with apoptosis induction and a functional type I interferon (IFN) response, through IRF3 activation and STAT1 phosphorylation. Using macrophages from mice deficient for the type I IFN receptor, we found that this type I IFN response was required for T. whipplei-induced macrophage apoptosis in a JNK-dependent manner and was associated with the intracellular replication of T. whipplei independently of JNK. This study underscores the role of macrophage polarization in host responses and highlights the detrimental role of type I IFN during T. whipplei infection

Global Analysis of Circulating Immune Cells by Matrix-Assisted Laser Desorption Ionization Time-of-Flight Mass Spectrometry

Background: MALDI-TOF mass spectrometry is currently used in microbiological diagnosis to characterize bacterial populations. Our aim was to determine whether this technique could be applied to intact eukaryotic cells, and in particular, to cells involved in the immune response. Methodology/Principal Findings: A comparison of frozen monocytes, T lymphocytes and polymorphonuclear leukocytes revealed specific peak profiles. We also found that twenty cell types had specific profiles, permitting the establishment of a cell database. The circulating immune cells, namely monocytes, T lymphocytes and polymorphonuclear cells, were distinct from tissue immune cells such as monocyte-derived macrophages and dendritic cells. In addition, MALDI-TOF mass spectrometry was valuable to easily identify the signatures of monocytes and T lymphocytes in peripheral mononuclear cells. Conclusions/Significance: This method was rapid and easy to perform, and unlike flow cytometry, it did not require any additional components such as specific antibodies. The MALDI-TOF mass spectrometry approach could be extended t

Metagenomic Analysis of Lysogeny in Tampa Bay: Implications for Prophage Gene Expression

Phage integrase genes often play a role in the establishment of lysogeny in temperate phage by catalyzing the integration of the phage into one of the host's replicons. To investigate temperate phage gene expression, an induced viral metagenome from Tampa Bay was sequenced by 454/Pyrosequencing. The sequencing yielded 294,068 reads with 6.6% identifiable. One hundred-three sequences had significant similarity to integrases by BLASTX analysis (e≤0.001). Four sequences with strongest amino-acid level similarity to integrases were selected and real-time PCR primers and probes were designed. Initial testing with microbial fraction DNA from Tampa Bay revealed 1.9×107, and 1300 gene copies of Vibrio-like integrase and Oceanicola-like integrase L−1 respectively. The other two integrases were not detected. The integrase assay was then tested on microbial fraction RNA extracted from 200 ml of Tampa Bay water sampled biweekly over a 12 month time series. Vibrio-like integrase gene expression was detected in three samples, with estimated copy numbers of 2.4-1280 L−1. Clostridium-like integrase gene expression was detected in 6 samples, with estimated copy numbers of 37 to 265 L−1. In all cases, detection of integrase gene expression corresponded to the occurrence of lysogeny as detected by prophage induction. Investigation of the environmental distribution of the two expressed integrases in the Global Ocean Survey Database found the Vibrio-like integrase was present in genome equivalents of 3.14% of microbial libraries and all four viral metagenomes. There were two similar genes in the library from British Columbia and one similar gene was detected in both the Gulf of Mexico and Sargasso Sea libraries. In contrast, in the Arctic library eleven similar genes were observed. The Clostridium-like integrase was less prevalent, being found in 0.58% of the microbial and none of the viral libraries. These results underscore the value of metagenomic data in discovering signature genes that play important roles in the environment through their expression, as demonstrated by integrases in lysogeny

Minimum information about an uncultivated virus genome (MIUVIG)

This is the final version. Available on open access from Nature Research via the DOI in this recordNOTE: the full list of funders and grants is in the acknowledgements section of the articleWe present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.Simons Foundation InternationalNatural Environment Research Council (NERC

The complete genome sequence and genetic analysis of ΦCA82 a novel uncultured microphage from the turkey gastrointestinal system

The genomic DNA sequence of a novel enteric uncultured microphage, ΦCA82 from a turkey gastrointestinal system was determined utilizing metagenomics techniques. The entire circular, single-stranded nucleotide sequence of the genome was 5,514 nucleotides. The ΦCA82 genome is quite different from other microviruses as indicated by comparisons of nucleotide similarity, predicted protein similarity, and functional classifications. Only three genes showed significant similarity to microviral proteins as determined by local alignments using BLAST analysis. ORF1 encoded a predicted phage F capsid protein that was phylogenetically most similar to the Microviridae ΦMH2K member's major coat protein. The ΦCA82 genome also encoded a predicted minor capsid protein (ORF2) and putative replication initiation protein (ORF3) most similar to the microviral bacteriophage SpV4. The distant evolutionary relationship of ΦCA82 suggests that the divergence of this novel turkey microvirus from other microviruses may reflect unique evolutionary pressures encountered within the turkey gastrointestinal system

Analysis of a viral metagenomic library from 200 m depth in Monterey Bay, California constructed by direct shotgun cloning

<p>Abstract</p> <p>Background</p> <p>Viruses have a profound influence on both the ecology and evolution of marine plankton, but the genetic diversity of viral assemblages, particularly those in deeper ocean waters, remains poorly described. Here we report on the construction and analysis of a viral metagenome prepared from below the euphotic zone in a temperate, eutrophic bay of coastal California.</p> <p>Methods</p> <p>We purified viruses from approximately one cubic meter of seawater collected from 200m depth in Monterey Bay, CA. DNA was extracted from the virus fraction, sheared, and cloned with no prior amplification into a plasmid vector and propagated in <it>E. coli </it>to produce the MBv200m library. Random clones were sequenced by the Sanger method. Sequences were assembled then compared to sequences in GenBank and to other viral metagenomic libraries using BLAST analyses.</p> <p>Results</p> <p>Only 26% of the 881 sequences remaining after assembly had significant (E ≤ 0.001) BLAST hits to sequences in the GenBank nr database, with most being matches to bacteria (15%) and viruses (8%). When BLAST analysis included environmental sequences, 74% of sequences in the MBv200m library had a significant match. Most of these hits (70%) were to microbial metagenome sequences and only 0.7% were to sequences from viral metagenomes. Of the 121 sequences with a significant hit to a known virus, 94% matched bacteriophages (Families <it>Podo</it>-, <it>Sipho</it>-, and <it>Myoviridae</it>) and 6% matched viruses of eukaryotes in the Family <it>Phycodnaviridae </it>(5 sequences) or the Mimivirus (2 sequences). The largest percentages of hits to viral genes of known function were to those involved in DNA modification (25%) or structural genes (17%). Based on reciprocal BLAST analyses, the MBv200m library appeared to be most similar to viral metagenomes from two other bays and least similar to a viral metagenome from the Arctic Ocean.</p> <p>Conclusions</p> <p>Direct cloning of DNA from diverse marine viruses was feasible and resulted in a distribution of virus types and functional genes at depth that differed in detail, but were broadly similar to those found in surface marine waters. Targeted viral analyses are useful for identifying those components of the greater marine metagenome that circulate in the subcellular size fraction.</p