Targeted whole-genome recovery of single viral species in a complex environmental sample.

Characterizing unknown viruses is essential for understanding viral ecology and preparing against viral outbreaks. Recovering complete genome sequences from environmental samples remains computationally challenging using metagenomics, especially for low-abundance species with uneven coverage. We present an experimental method for reliably recovering complete viral genomes from complex environmental samples. Individual genomes are encapsulated into droplets and amplified using multiple displacement amplification. A unique gene detection assay, which employs an RNA-based probe and an exonuclease, selectively identifies droplets containing the target viral genome. Labeled droplets are sorted using a microfluidic sorter, and genomes are extracted for sequencing. We demonstrate this method's efficacy by spiking two known viral genomes, Simian virus 40 (SV40, 5,243 bp) and Human Adenovirus 5 (HAd5, 35,938 bp), into a sewage sample with a final abundance in the droplets of around 0.1% and 0.015%, respectively. We achieve 100% recovery of the complete sequence of the spiked-in SV40 genome with uniform coverage distribution. For the larger HAd5 genome, we cover approximately 99.4% of its sequence. Notably, genome recovery is achieved with as few as one sorted droplet, which enables the recovery of any desired genomes in complex environmental samples, regardless of their abundance. This method enables single-genome whole-genome amplification and targeting characterizations of rare viral species and will facilitate our ability to access the mutational profile in single-virus genomes and contribute to an improved understanding of viral ecology.

Cultured Renal Proximal Tubular Epithelial Cells Resemble a Stressed/Damaged Kidney While Supporting BK Virus Infection.

BK virus (BKV; human polyomavirus 1) infections are asymptomatic in most individuals, and the virus persists throughout life without harm. However, BKV is a threat to transplant patients and those with immunosuppressive disorders. Under these circumstances, the virus can replicate robustly in proximal tubule epithelial cells (PT). Cultured renal proximal tubule epithelial cells (RPTE) are permissive to BKV and have been used extensively to characterize different aspects of BKV infection. Recently, lines of hTERT-immortalized RPTE have become available, and preliminary studies indicate they support BKV infection as well. Our results indicate that BKV infection leads to a similar response in primary and immortalized RPTE. In addition, we examined the patterns of global gene expression of primary and immortalized RPTE and compared them with uncultured PT freshly dissociated from human kidney. As expected, PT isolated from the healthy kidney express a number of differentiation-specific genes that are associated with kidney function. However, the expression of most of these genes is absent or repressed in cultured RPTE. Rather, cultured RPTE exhibit a gene expression profile indicative of a stressed or injured kidney. Inoculation of cultured RPTE with BKV results in the suppression of many genes associated with kidney stress. In summary, this study demonstrated similar global gene expression patterns and responses to BKV infection between primary and immortalized RPTE. Moreover, results from bulk transcriptome sequencing (RNA-seq) and SCT experiments revealed distinct transcriptomic signatures representing cell injury and stress in primary RPTE in contrast to the uncultured, freshly dissociated PT from human kidney. IMPORTANCE Cultured primary human cells provide powerful tools for the study of viral infectious cycles and host virus interactions. In the case of BKV-associated nephropathy, viral replication occurs primarily in the proximal tubule epithelia in the kidney. Consequently, cultured primary and immortalized renal proximal tubule epithelial cells (RPTE) are widely used to study BKV infection. In this work, using bulk and single-cell transcriptomics, we found that primary and immortalized RPTE responded similarly to BKV infection. However, both uninfected primary and immortalized RPTE have gene expression profiles that are markedly different from healthy proximal tubule epithelia isolated directly from human kidney without culture. Cultured RPTE are in a gene expression state indicative of an injured or stressed kidney. These results raise the possibility that BKV replicates preferentially in injured or stressed kidney epithelial cells during nephropathy.

Single-Cell Transcriptomics Reveals a Heterogeneous Cellular Response to BK Virus Infection.

BK virus (BKV) is a human polyomavirus that is generally harmless but can cause devastating disease in immunosuppressed individuals. BKV infection of renal cells is a common problem for kidney transplant patients undergoing immunosuppressive therapy. In cultured primary human renal proximal tubule epithelial (RPTE) cells, BKV undergoes a productive infection. The BKV-encoded large T antigen (LT) induces cell cycle entry, resulting in the upregulation of numerous genes associated with cell proliferation. Consistently, microarray and transcriptome sequencing (RNA-seq) experiments performed on bulk infected cell populations identified several proliferation-related pathways that are upregulated by BKV. These studies revealed few genes that are downregulated. In this study, we analyzed viral and cellular transcripts in single mock- or BKV-infected cells. We found that the levels of viral mRNAs vary widely among infected cells, resulting in different levels of LT and viral capsid protein expression. Cells expressing the highest levels of viral transcripts account for approximately 20% of the culture and have a gene expression pattern that is distinct from that of cells expressing lower levels of viral mRNAs. Surprisingly, cells expressing low levels of viral mRNA do not progress with time to high expression, suggesting that the two cellular responses are determined prior to or shortly following infection. Finally, comparison of cellular gene expression patterns of cells expressing high levels of viral mRNA with those of mock-infected cells or cells expressing low levels of viral mRNA revealed previously unidentified pathways that are downregulated by BKV. Among these are pathways associated with drug metabolism and detoxification, tumor necrosis factor (TNF) signaling, energy metabolism, and translation.IMPORTANCE The outcome of viral infection is determined by the ability of the virus to redirect cellular systems toward progeny production countered by the ability of the cell to block these viral actions. Thus, an infected culture consists of thousands of cells, each fighting its own individual battle. Bulk measurements, such as PCR or RNA-seq, measure the average of these individual responses to infection. Single-cell transcriptomics provides a window to the one-on-one battle between BKV and each cell. Our studies reveal that only a minority of infected cells are overwhelmed by the virus and produce large amounts of BKV mRNAs and proteins, while the infection appears to be restricted in the remaining cells. Correlation of viral transcript levels with cellular gene expression patterns reveals pathways manipulated by BKV that may play a role in limiting infection.

Human polyomavirus BKV infection of endothelial cells results in interferon pathway induction and persistence.

Polyomavirus BKV is highly prevalent among humans. The virus establishes an asymptomatic persistent infection in the urinary system in healthy people, but uncontrolled productive infection of the virus in immunocompromised patients can lead to serious diseases. In spite of its high prevalence, our knowledge regarding key aspects of BKV polyomavirus infection remains incomplete. To determine tissue and cell type tropism of the virus, primary human epithelial cells, endothelial cells and fibroblasts isolated from the respiratory and urinary systems were tested. Results from this study demonstrated that all 9 different types of human cells were infectable by BKV polyomavirus but showed differential cellular responses. In microvascular endothelial cells from the lung and the bladder, BKV persistent infection led to prolonged viral protein expression, low yield of infectious progeny and delayed cell death, in contrast with infection in renal proximal tubular epithelial cells, a widely used cell culture model for studying productive infection of this virus. Transcriptomic profiling revealed the activation of interferon signaling and induction of multiple interferon stimulated genes in infected microvascular endothelial cells. Further investigation demonstrated production of IFNß and secretion of chemokine CXCL10 by infected endothelial cells. Activation of IRF3 and STAT1 in infected endothelial cells was also confirmed. In contrast, renal proximal tubular epithelial cells failed to mount an interferon response and underwent progressive cell death. These results demonstrated that microvascular endothelial cells are able to activate interferon signaling in response to polyomavirus BKV infection. This raises the possibility that endothelial cells might provide initial immune defense against BKV infection. Our results shed light on the persistence of and immunity against infection by BKV polyomavirus.

DNA Tumor Viruses and Their Contributions to Molecular Biology.

This summer marks the 51st anniversary of the DNA tumor virus meetings. Scientists from around the world will gather in Trieste, Italy, to report their latest results and to agree or disagree on the current concepts that define our understanding of this diverse class of viruses. This article offers a brief history of the impact the study of these viruses has had on molecular and cancer biology and discusses obstacles and opportunities for future progress.

Analysis of viruses present in urine from patients with interstitial cystitis.

The question of whether some cases of interstitial cystitis may have an infectious etiology has been debated for some time. Previous studies have looked for the presence of certain specific viruses, but generally did not use the types of sensitive and unbiased approaches that are currently available. As part of the MAPP (Multidisciplinary Approach to the Study of Chronic Pelvic Pain) Research Network, we examined urine specimens from interstitial cystitis patients who provided specimens over time and also reported various symptoms at the time of urine collection. We first performed next-generation sequencing to look for the presence of viruses in urines, and detected two human polyomaviruses that are known to be excreted into urine, BKPyV and JCPyV. We were especially interested in BKPyV because it is a known cause of another bladder disease, hemorrhagic cystitis, in bone marrow transplant recipients. Further analysis of individual samples indicates a trend toward higher excretion of polyomaviruses in patients experiencing increased symptoms.

Evolution on the Biophysical Fitness Landscape of an RNA Virus.

Viral evolutionary pathways are determined by the fitness landscape, which maps viral genotype to fitness. However, a quantitative description of the landscape and the evolutionary forces on it remain elusive. Here, we apply a biophysical fitness model based on capsid folding stability and antibody binding affinity to predict the evolutionary pathway of norovirus escaping a neutralizing antibody. The model is validated by experimental evolution in bulk culture and in a drop-based microfluidics that propagates millions of independent small viral subpopulations. We demonstrate that along the axis of binding affinity, selection for escape variants and drift due to random mutations have the same direction, an atypical case in evolution. However, along folding stability, selection and drift are opposing forces whose balance is tuned by viral population size. Our results demonstrate that predictable epistatic tradeoffs between molecular traits of viral proteins shape viral evolution.

The Ancient Evolutionary History of Polyomaviruses.

Polyomaviruses are a family of DNA tumor viruses that are known to infect mammals and birds. To investigate the deeper evolutionary history of the family, we used a combination of viral metagenomics, bioinformatics, and structural modeling approaches to identify and characterize polyomavirus sequences associated with fish and arthropods. Analyses drawing upon the divergent new sequences indicate that polyomaviruses have been gradually co-evolving with their animal hosts for at least half a billion years. Phylogenetic analyses of individual polyomavirus genes suggest that some modern polyomavirus species arose after ancient recombination events involving distantly related polyomavirus lineages. The improved evolutionary model provides a useful platform for developing a more accurate taxonomic classification system for the viral family Polyomaviridae.

Large T antigens of polyomaviruses: amazing molecular machines.

The large tumor antigen (T antigen) encoded by simian virus 40 is an amazing molecular machine because it orchestrates viral infection by modulating multiple fundamental viral and cellular processes. T antigen is required for viral DNA replication, transcription, and virion assembly. In addition, T antigen targets multiple cellular pathways, including those that regulate cell proliferation, cell death, and the inflammatory response. Ectopic T antigen expression results in the immortalization and transformation of many cell types in culture and T antigen induces neoplasia when expressed in rodents. The analysis of the mechanisms by which T antigen carries out its many functions has proved to be a powerful way of gaining insights into cell biology. The accelerating pace at which new polyomaviruses are being discovered provides a collection of novel T antigens that, like simian virus 40, can be used to discover and study key cellular regulatory systems.

Cellular transformation of mouse embryo fibroblasts in the absence of activator E2Fs.

UNLABELLED: The E2F family of transcription factors, broadly divided into activator and repressor E2Fs, regulates cell cycle genes. Current models indicate that activator E2Fs are necessary for cell cycle progression and tumorigenesis and are also required to mediate transformation induced by DNA tumor viruses. E2Fs are negatively regulated by the retinoblastoma (RB) family of tumor suppressor proteins, and virus-encoded oncogenes disrupt the RB-E2F repressor complexes. This results in the release of activator E2Fs and induction of E2F-dependent genes. In agreement, expression of large tumor T antigens (TAg) encoded by polyomaviruses in mammalian cells results in increased transcriptional levels of E2F target genes. In addition, tumorigenesis induced by transgenic expression of simian virus 40 (SV40) TAg in choroid plexus or intestinal villi requires at least one activator E2F. In contrast, we show that SV40 TAg-induced transformation in mouse embryonic fibroblasts is independent of activator E2Fs. This work, coupled with recent studies showing that proliferation in stem and progenitor cells is independent of activator E2Fs, suggests the presence of parallel pathways governing cell proliferation and tumorigenesis. IMPORTANCE: The RB-E2F pathway is altered in many cancers and is also targeted by DNA tumor viruses. Viral oncoprotein action on RBs results in the release of activator E2Fs and upregulation of E2F target genes; thus, activator E2Fs are considered essential for normal and tumorigenic cell proliferation. However, we have observed that SV40 large T antigen can induce cell proliferation and transformation in the absence of activator E2Fs. Our results also suggest that TAg action on pRBs regulates both E2F-dependent and -independent pathways that govern proliferation. Thus, specific cell proliferation pathways affected by RB alterations in cancer may be a factor in tumor behavior and response to therapy.

HeLa nucleic acid contamination in the cancer genome atlas leads to the misidentification of human papillomavirus 18.

UNLABELLED: We searched The Cancer Genome Atlas (TCGA) database for viruses by comparing non-human reads present in transcriptome sequencing (RNA-Seq) and whole-exome sequencing (WXS) data to viral sequence databases. Human papillomavirus 18 (HPV18) is an etiologic agent of cervical cancer, and as expected, we found robust expression of HPV18 genes in cervical cancer samples. In agreement with previous studies, we also found HPV18 transcripts in non-cervical cancer samples, including those from the colon, rectum, and normal kidney. However, in each of these cases, HPV18 gene expression was low, and single-nucleotide variants and positions of genomic alignments matched the integrated portion of HPV18 present in HeLa cells. Chimeric reads that match a known virus-cell junction of HPV18 integrated in HeLa cells were also present in some samples. We hypothesize that HPV18 sequences in these non-cervical samples are due to nucleic acid contamination from HeLa cells. This finding highlights the problems that contamination presents in computational virus detection pipelines. IMPORTANCE: Viruses associated with cancer can be detected by searching tumor sequence databases. Several studies involving searches of the TCGA database have reported the presence of HPV18, a known cause of cervical cancer, in a small number of additional cancers, including those of the rectum, kidney, and colon. We have determined that the sequences related to HPV18 in non-cervical samples are due to nucleic acid contamination from HeLa cells. To our knowledge, this is the first report of the misidentification of viruses in next-generation sequencing data of tumors due to contamination with a cancer cell line. These results raise awareness of the difficulty of accurately identifying viruses in human sequence databases.

Isolation and Analysis of Rare Norovirus Recombinants from Coinfected Mice Using Drop-Based Microfluidics.

UNLABELLED: Human noroviruses (HuNoVs) are positive-sense RNA viruses that can cause severe, highly infectious gastroenteritis. HuNoV outbreaks are frequently associated with recombination between circulating strains. Strain genotyping and phylogenetic analyses show that noroviruses often recombine in a highly conserved region near the junction of the viral polyprotein (open reading frame 1 [ORF1]) and capsid (ORF2) genes and occasionally within the RNA-dependent RNA polymerase (RdRP) gene. Although genotyping methods are useful for tracking changes in circulating viral populations, they report only the dominant recombinant strains and do not elucidate the frequency or range of recombination events. Furthermore, the relatively low frequency of recombination in RNA viruses has limited studies to cell culture or in vitro systems, which do not reflect the complexities and selective pressures present in an infected organism. Using two murine norovirus (MNV) strains to model coinfection, we developed a microfluidic platform to amplify, detect, and recover individual recombinants following in vitro and in vivo coinfection. One-step reverse transcriptase PCR (RT-PCR) was performed in picoliter drops with primers that identified the wild-type and recombinant progenies and scanned for recombination breakpoints at ∼1-kb intervals. We detected recombination between MNV strains at multiple loci spanning the viral protease, RdRP, and capsid ORFs and isolated individual recombinant RNA genomes that were present at a frequency of 1/300,000 or higher. This study is the first to examine norovirus recombination following coinfection of an animal and suggests that the exchange of RNA among viral genomes in an infected host occurs in multiple locations and is an important driver of genetic diversity. IMPORTANCE: RNA viruses increase diversity and escape host immune barriers by genomic recombination. Studies using a number of viral systems indicate that recombination occurs via template switching by the virus-encoded RNA-dependent RNA polymerase (RdRP). However, factors that govern the frequency and positions of recombination in an infected organism remain largely unknown. This work leverages advances in the applied physics of drop-based microfluidics to isolate and sequence rare recombinants arising from the coinfection of mice with two distinct strains of murine norovirus. This study is the first to detect and analyze norovirus recombination in an animal model.

Simian virus 40 large T antigen induces IFN-stimulated genes through ATR kinase.

Polyomaviruses encode a large T Ag (LT), a multifunctional protein essential for the regulation of both viral and host cell gene expression and productive viral infection. Previously, we have shown that stable expression of LT protein results in upregulation of genes involved in the IFN induction and signaling pathway. In this study, we focus on the cellular signaling mechanism that leads to the induction of IFN responses by LT. Our results show that ectopic expression of SV40 LT results in the induction of IFN-stimulated genes (ISGs) in human fibroblasts and confers an antiviral state. We describe a LT-initiated DNA damage response (DDR) that activates IFN regulatory factor 1, causing IFN-ß production and consequent ISG expression in human cells. This IFN-ß and ISG induction is dependent on ataxia-telangiectasia mutated and Rad3-related (ATR) kinase, but independent of ATM. ATR kinase inhibition using a selective kinase inhibitor (ETP-46464) caused a decrease in IFN regulatory factor 1 stabilization and ISG expression. Furthermore, expression of a mutant LT that does not induce DDR also does not induce IFN-ß and ISGs. These results show that, in the absence of viral infection, LT-initiated activation of ATR-dependent DDR is sufficient for the induction of an IFN-ß-mediated innate immune response in human cells. Thus, we have uncovered a novel and critical role for ATR as a mediator of antiviral responses utilizing LT.

Artifact-Free Quantification and Sequencing of Rare Recombinant Viruses by Using Drop-Based Microfluidics.

Recombination is an important driver in the evolution of viruses and thus is key to understanding viral epidemics and improving strategies to prevent future outbreaks. Characterization of rare recombinant subpopulations remains technically challenging because of artifacts such as artificial recombinants, known as chimeras, and amplification bias. To overcome this, we have developed a high-throughput microfluidic technique with a second verification step in order to amplify and sequence single recombinant viruses with high fidelity in picoliter drops. We obtained the first artifact-free estimate of in vitro recombination rate between murine norovirus strains MNV-1 and WU20 co-infecting a cell (P(rec) = 3.3 × 10(-4) ± 2 × 10(-5) ) for a 1205 nt region. Our approach represents a time- and cost-effective improvement over current methods, and can be adapted for genomic studies requiring artifact- and bias-free selective amplification, such as microbial pathogens, or rare cancer cells.

Flexibility in surface-exposed loops in a virus capsid mediates escape from antibody neutralization.

UNLABELLED: New human norovirus strains emerge every 2 to 3 years, partly due to mutations in the viral capsid that allow escape from antibody neutralization and herd immunity. To understand how noroviruses evolve antibody resistance, we investigated the structural basis for the escape of murine norovirus (MNV) from antibody neutralization. To identify specific residues in the MNV-1 protruding (P) domain of the capsid that play a role in escape from the neutralizing monoclonal antibody (MAb) A6.2, 22 recombinant MNVs were generated with amino acid substitutions in the A'B' and E'F' loops. Six mutations in the E'F' loop (V378F, A382K, A382P, A382R, D385G, and L386F) mediated escape from MAb A6.2 neutralization. To elucidate underlying structural mechanisms for these results, the atomic structure of the A6.2 Fab was determined and fitted into the previously generated pseudoatomic model of the A6.2 Fab/MNV-1 virion complex. Previously, two distinct conformations, A and B, of the atomic structures of the MNV-1 P domain were identified due to flexibility in the two P domain loops. A superior stereochemical fit of the A6.2 Fab to the A conformation of the MNV P domain was observed. Structural analysis of our observed escape mutants indicates changes toward the less-preferred B conformation of the P domain. The shift in the structural equilibrium of the P domain toward the conformation with poor structural complementarity to the antibody strongly supports a unique mechanism for antibody escape that occurs via antigen flexibility instead of direct antibody-antigen binding. IMPORTANCE: Human noroviruses cause the majority of all nonbacterial gastroenteritis worldwide. New epidemic strains arise in part by mutations in the viral capsid leading to escape from antibody neutralization. Herein, we identify a series of point mutations in a norovirus capsid that mediate escape from antibody neutralization and determine the structure of a neutralizing antibody. Fitting of the antibody structure into the virion/antibody complex identifies two conformations of the antibody binding domain of the viral capsid: one with a superior fit and the other with an inferior fit to the antibody. These data suggest a unique mode of antibody neutralization. In contrast to other viruses that largely escape antibody neutralization through direct disruption of the antibody-virus interface, we identify mutations that acted indirectly by limiting the conformation of the antibody binding loop in the viral capsid and drive the antibody binding domain into the conformation unable to be bound by the antibody.

The conserved core enzymatic activities and the distinct dynamics of polyomavirus large T antigens.

Several human polyomaviruses including JCV, BKV and TSV are associated with diseases, particularly in immunosuppressed patients. While the large T antigen (LT) encoded by the monkey polyomavirus SV40 is well studied, and possesses intrinsic ATPase and DNA helicase activities, the LTs of the human polyomaviruses are relatively uncharacterized. In order to evaluate whether these enzymatic activities, which are required for viral DNA replication, are conserved between polyomaviruses, we performed a comparative study using the LTs from JCV, TSV and SV40. The ATPase and DNA helicase activities and the interaction with the cellular tumor suppressor p53 were assayed for the purified Zn-ATPase domains of the three LTs. We found that all Zn-ATPases were active ATPases. The Zn-ATPase domains also functioned as DNA helicases, although the measured kinetic constants differed among the three proteins. In addition, when tested against four small molecule ATPase inhibitors, the Zn-ATPase domains of TSV was more resistant than that of SV40 and JCV. Our results show that, while LTs from JCV and TSV share the core ATPase and DNA helicase activities, they possess important functional differences that might translate into their respective abilities to infect and replicate in hosts.

E2f1-3 switch from activators in progenitor cells to repressors in differentiating cells.

In the established model of mammalian cell cycle control, the retinoblastoma protein (Rb) functions to restrict cells from entering S phase by binding and sequestering E2f activators (E2f1, E2f2 and E2f3), which are invariably portrayed as the ultimate effectors of a transcriptional program that commit cells to enter and progress through S phase. Using a panel of tissue-specific cre-transgenic mice and conditional E2f alleles we examined the effects of E2f1, E2f2 and E2f3 triple deficiency in murine embryonic stem cells, embryos and small intestines. We show that in normal dividing progenitor cells E2f1-3 function as transcriptional activators, but contrary to the current view, are dispensable for cell division and instead are necessary for cell survival. In differentiating cells E2f1-3 function in a complex with Rb as repressors to silence E2f targets and facilitate exit from the cell cycle. The inactivation of Rb in differentiating cells resulted in a switch of E2f1-3 from repressors to activators, leading to the superactivation of E2f responsive targets and ectopic cell divisions. Loss of E2f1-3 completely suppressed these phenotypes caused by Rb deficiency. This work contextualizes the activator versus repressor functions of E2f1-3 in vivo, revealing distinct roles in dividing versus differentiating cells and in normal versus cancer-like cell cycles.

Identification of a novel polyomavirus in a pancreatic transplant recipient with retinal blindness and vasculitic myopathy.

BACKGROUND: A 33 year-old pancreatic transplant recipient developed weakness, retinal blindness, and necrotic plaques on her face, scalp, and hands. METHODS: A muscle biopsy was analyzed by light and electron microscopy and high-throughput nucleic acid sequencing. RESULTS: The biopsy revealed microthrombosis and viral particles in swollen endothelial cell nuclei. High-throughput sequencing of nucleic acid revealed a novel polyomavirus. In situ hybridization confirmed the presence of the polyomavirus in endothelial cells at sites of myositis and cutaneous necrosis. CONCLUSIONS: New Jersey polyomavirus (NJPyV-2013) is a novel polyomavirus that may have tropism for vascular endothelial cells.

Whole-Genome Sequencing of a Single Viral Species from a Highly Heterogeneous Sample.

Metagenomic studies suggest that only a small fraction of the viruses that exist in nature have been identified and studied. Characterization of unknown viral genomes is hindered by the many genomes populating any virus sample. A new method is reported that integrates drop-based microfluidics and computational analysis to enable the purification of any single viral species from a complex mixed virus sample and the retrieval of complete genome sequences. By using this platform, the genome sequence of a 5243 bp dsDNA virus that was spiked into wastewater was retrieved with greater than 96% sequence coverage and more than 99.8% sequence identity. This method holds great potential for virus discovery since it allows enrichment and sequencing of previously undescribed viruses as well as known viruses.

SummonChimera infers integrated viral genomes with nucleotide precision from NGS data.

BACKGROUND: Viral integration into a host genome is defined by two chimeric junctions that join viral and host DNA. Recently, computational tools have been developed that utilize NGS data to detect chimeric junctions. These methods identify individual viral-host junctions but do not associate chimeric pairs as an integration event. Without knowing the chimeric boundaries of an integration, its genetic content cannot be determined. RESULTS: Summonchimera is a Perl program that associates chimera pairs to infer the complete viral genomic integration event to the nucleotide level within single or paired-end NGS data. SummonChimera integration prediction was verified on a set of single-end IonTorrent reads from a purified Salmonella bacterium with an integrated bacteriophage. Furthermore, SummonChimera predicted integrations from experimentally verified Hepatitis B Virus chimeras within a paired-end Whole Genome Sequencing hepatocellular carcinoma tumor database. CONCLUSIONS: SummonChimera identified all experimentally verified chimeras detected by current computational methods. Further, SummonChimera integration inference precisely predicted bacteriophage integration. The application of SummonChimera to cancer NGS accurately identifies deletion of host and viral sequence during integration. The precise nucleotide determination of an integration allows prediction of viral and cellular gene transcription patterns.