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.

High-sensitivity whole-genome recovery of single viral species in environmental samples.

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. This work presents a method for reliably recovering complete viral genomes from complex environmental samples. Individual genomes are encapsulated into droplets and amplified using multiple displacement amplification. A novel 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. Validation experiments using a sewage sample spiked with two known viruses demonstrate the method's efficacy. We achieve 100% recovery of the spiked-in SV40 (Simian virus 40, 5243bp) genome sequence with uniform coverage distribution, and approximately 99.4% for the larger HAd5 genome (Human Adenovirus 5, 35938bp). 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 targeted characterizations of rare viral species and whole-genome amplification of single genomes for accessing the mutational profile in single virus genomes, contributing 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.

Recurrent integration of human papillomavirus genomes at transcriptional regulatory hubs.

Oncogenic human papillomavirus (HPV) genomes are often integrated into host chromosomes in HPV-associated cancers. HPV genomes are integrated either as a single copy or as tandem repeats of viral DNA interspersed with, or without, host DNA. Integration occurs frequently in common fragile sites susceptible to tandem repeat formation and the flanking or interspersed host DNA often contains transcriptional enhancer elements. When co-amplified with the viral genome, these enhancers can form super-enhancer-like elements that drive high viral oncogene expression. Here we compiled highly curated datasets of HPV integration sites in cervical (CESC) and head and neck squamous cell carcinoma (HNSCC) cancers, and assessed the number of breakpoints, viral transcriptional activity, and host genome copy number at each insertion site. Tumors frequently contained multiple distinct HPV integration sites but often only one "driver" site that expressed viral RNA. As common fragile sites and active enhancer elements are cell-type-specific, we mapped these regions in cervical cell lines using FANCD2 and Brd4/H3K27ac ChIP-seq, respectively. Large enhancer clusters, or super-enhancers, were also defined using the Brd4/H3K27ac ChIP-seq dataset. HPV integration breakpoints were enriched at both FANCD2-associated fragile sites and enhancer-rich regions, and frequently showed adjacent focal DNA amplification in CESC samples. We identified recurrent integration "hotspots" that were enriched for super-enhancers, some of which function as regulatory hubs for cell-identity genes. We propose that during persistent infection, extrachromosomal HPV minichromosomes associate with these transcriptional epicenters and accidental integration could promote viral oncogene expression and carcinogenesis.

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.

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.

Detecting viral sequences in NGS data.

Next generation sequencing (NGS) technologies provide an increasingly important avenue for detecting known viruses, and for discovering novel viruses present in clinical or environmental samples. Several computational pipelines capable of identifying and classifying viral sequences in NGS data have been developed and used to search for viruses in human or animal samples, microbiomes, and in various environments. In this review we summarize the different approaches used to determine viral presence in sequence data. Strategies for avoiding confounding factors such as physical contamination and computational artifacts that lead to false virus identification are discussed. The application of these methodologies to cancer data sets has led to important insights on viruses both as drivers of and biomarkers for specific tumors.

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.

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.

Identification of Head and Neck Cancer Subtypes Based on Human Papillomavirus Presence and E2F-Regulated Gene Expression.

Human papillomavirus (HPV) is present in a subset of head and neck squamous cell carcinomas (HNSCCs). The cell cycle regulatory Rb-E2F pathway is a major target of HPV and is perturbed by these viruses in cell culture and animal models, as well as in human tumors. In this study, we examined differences in the Rb-E2F pathway displayed by HPV-positive (HPV+) and HPV-negative (HPV-) HNSCC tumors. We created a computational approach that effectively categorizes gene expression as unchanged, downregulated, or upregulated by comparing the gene's mRNA levels in the tumor to the corresponding mRNA levels across normal tissue samples. Our findings suggest that there are three major HNSCC subtypes, defined by differences in the presence of HPV and in E2F-regulated gene expression. Most HPV+ HNSCC tumors show upregulation of E2F-regulated genes, which is consistent with inactivation of Rb by the virus-encoded E7 protein. In contrast, many HPV- HNSCCs show little or no change in the Rb-E2F pathway. However, we also identified a set of tumors that show alterations in the Rb-E2F pathway in the absence of HPV. Thus, one class of HPV- HNSCCs arise without significant alterations of the Rb-E2F pathway, while a second class of tumors appear to deregulate this pathway independently of the presence of HPV. IMPORTANCE Cancer is a complex disease that can be caused by a multitude of factors. HNSCC is complicated because some of these cancers are clearly associated with HPV, while others have no viral involvement. Determining the pathways that are commonly altered in both types of HNSCC, as well as those that are unique to viral and nonviral tumors, is important for a basic understanding of how these cancers arise and progress and critical to the development of targeted therapies. In this work, we show that all HPV-associated tumors have increased expression of E2F target genes, indicating that the tumor suppressor function of Rb is blocked. Importantly, Rb is also inhibited in a subset of nonviral tumors, suggesting that mutations present in these cancers mimic the action of the HPV E6 and E7 oncogenes.

Viral sequences in human cancer.

We have developed a virus detection and discovery computational pipeline, Pickaxe, and applied it to NGS databases provided by The Cancer Genome Atlas (TCGA). We analyzed a collection of whole genome (WGS), exome (WXS), and RNA (RNA-Seq) sequencing libraries from 3052 participants across 22 different cancers. NGS data from nearly all tumor and normal tissues examined contained contaminating viral sequences. Intensive computational and manual efforts are required to remove these artifacts. We found that several different types of cancers harbored Herpesviruses including EBV, CMV, HHV1, HHV2, HHV6 and HHV7. In addition to the reported associations of Hepatitis B and C virus (HBV & HCV) with liver cancer, and Human papillomaviruses (HPV) with cervical cancer and a subset of head and neck cancers, we found additional cases of HPV integrated in a small number of bladder cancers. Gene expression and mutational profiles suggest that HPV drives tumorigenesis in these cases.

Complete Genome Sequence of a Polyomavirus Recovered from a Pomona Leaf-Nosed Bat (Hipposideros pomona) Metagenome Data Set.

We report here the complete genome sequence of a polyomavirus found in a nasal/rectal metagenome of Hipposideros pomona (Pomona leaf-nosed bat). Interestingly, the genetic organization and phylogenetic relationships of the new virus suggest greater similarity to recently discovered fish-associated polyomaviruses rather than to polyomavirus species previously observed in bats.

Merkel Cell Polyomavirus Exhibits Dominant Control of the Tumor Genome and Transcriptome in Virus-Associated Merkel Cell Carcinoma.

Merkel cell polyomavirus is the primary etiological agent of the aggressive skin cancer Merkel cell carcinoma (MCC). Recent studies have revealed that UV radiation is the primary mechanism for somatic mutagenesis in nonviral forms of MCC. Here, we analyze the whole transcriptomes and genomes of primary MCC tumors. Our study reveals that virus-associated tumors have minimally altered genomes compared to non-virus-associated tumors, which are dominated by UV-mediated mutations. Although virus-associated tumors contain relatively small mutation burdens, they exhibit a distinct mutation signature with observable transcriptionally biased kataegic events. In addition, viral integration sites overlap focal genome amplifications in virus-associated tumors, suggesting a potential mechanism for these events. Collectively, our studies indicate that Merkel cell polyomavirus is capable of hijacking cellular processes and driving tumorigenesis to the same severity as tens of thousands of somatic genome alterations. IMPORTANCE: A variety of mutagenic processes that shape the evolution of tumors are critical determinants of disease outcome. Here, we sequenced the entire genome of virus-positive and virus-negative primary Merkel cell carcinomas (MCCs), revealing distinct mutation spectra and corresponding expression profiles. Our studies highlight the strong effect that Merkel cell polyomavirus has on the divergent development of viral MCC compared to the somatic alterations that typically drive nonviral tumorigenesis. A more comprehensive understanding of the distinct mutagenic processes operative in viral and nonviral MCCs has implications for the effective treatment of these tumors.

Dosage-dependent copy number gains in E2f1 and E2f3 drive hepatocellular carcinoma.

Disruption of the retinoblastoma (RB) tumor suppressor pathway, either through genetic mutation of upstream regulatory components or mutation of RB1 itself, is believed to be a required event in cancer. However, genetic alterations in the RB-regulated E2F family of transcription factors are infrequent, casting doubt on a direct role for E2Fs in driving cancer. In this work, a mutation analysis of human cancer revealed subtle but impactful copy number gains in E2F1 and E2F3 in hepatocellular carcinoma (HCC). Using a series of loss- and gain-of-function alleles to dial E2F transcriptional output, we have shown that copy number gains in E2f1 or E2f3b resulted in dosage-dependent spontaneous HCC in mice without the involvement of additional organs. Conversely, germ-line loss of E2f1 or E2f3b, but not E2f3a, protected mice against HCC. Combinatorial mapping of chromatin occupancy and transcriptome profiling identified an E2F1- and E2F3B-driven transcriptional program that was associated with development and progression of HCC. These findings demonstrate a direct and cell-autonomous role for E2F activators in human cancer.

E2f8 mediates tumor suppression in postnatal liver development.

E2F-mediated transcriptional repression of cell cycle-dependent gene expression is critical for the control of cellular proliferation, survival, and development. E2F signaling also interacts with transcriptional programs that are downstream of genetic predictors for cancer development, including hepatocellular carcinoma (HCC). Here, we evaluated the function of the atypical repressor genes E2f7 and E2f8 in adult liver physiology. Using several loss-of-function alleles in mice, we determined that combined deletion of E2f7 and E2f8 in hepatocytes leads to HCC. Temporal-specific ablation strategies revealed that E2f8's tumor suppressor role is critical during the first 2 weeks of life, which correspond to a highly proliferative stage of postnatal liver development. Disruption of E2F8's DNA binding activity phenocopied the effects of an E2f8 null allele and led to HCC. Finally, a profile of chromatin occupancy and gene expression in young and tumor-bearing mice identified a set of shared targets for E2F7 and E2F8 whose increased expression during early postnatal liver development is associated with HCC progression in mice. Increased expression of E2F8-specific target genes was also observed in human liver biopsies from HCC patients compared to healthy patients. In summary, these studies suggest that E2F8-mediated transcriptional repression is a critical tumor suppressor mechanism during postnatal liver development.

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.

Expression of the small T antigen of Lymphotropic Papovavirus is sufficient to transform primary mouse embryo fibroblasts.

Polyomaviruses induce cell proliferation and transformation through different oncoproteins encoded within the early region (ER): large T antigen (LT), small T antigen (sT) and, in some cases, additional components. Each virus utilizes different mechanisms to achieve transformation. For instance, the LTs of Simian virus 40 (SV40), BK and/or JC virus can induce transformation; but Merkel Cell Polyomavirus (MCPyV) requires expression of sT. Lymphotropic Papovavirus (LPV) is closely related to Human Polyomavirus 9 (HuPyV9) and, under similar conditions, mice expressing LPV.ER exhibit higher rates of tumor formation than mice expressing SV40.ER. We have investigated the contributions of individual LPV.ER components to cell transformation. In contrast to SV40, LPV.ER transforms mouse embryonic fibroblasts (MEFs), but expression of LPV LT is insufficient to transform MEFs. Furthermore, LPV sT induces immortalization and transformation of MEFs. Thus, in the case of LPV, sT is the main mediator of oncogenesis.

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.

Redeployment of Myc and E2f1-3 drives Rb-deficient cell cycles.

Robust mechanisms to control cell proliferation have evolved to maintain the integrity of organ architecture. Here, we investigated how two critical proliferative pathways, Myc and E2f, are integrated to control cell cycles in normal and Rb-deficient cells using a murine intestinal model. We show that Myc and E2f1-3 have little impact on normal G1-S transitions. Instead, they synergistically control an S-G2 transcriptional program required for normal cell divisions and maintaining crypt-villus integrity. Surprisingly, Rb deficiency results in the Myc-dependent accumulation of E2f3 protein and chromatin repositioning of both Myc and E2f3, leading to the 'super activation' of a G1-S transcriptional program, ectopic S phase entry and rampant cell proliferation. These findings reveal that Rb-deficient cells hijack and redeploy Myc and E2f3 from an S-G2 program essential for normal cell cycles to a G1-S program that re-engages ectopic cell cycles, exposing an unanticipated addiction of Rb-null cells on Myc.

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.