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.

A metagenomic analysis of the phase 2 Anopheles gambiae 1000 genomes dataset reveals a wide diversity of cobionts associated with field collected mosquitoes.

The Anopheles gambiae 1000 Genomes (Ag1000G) Consortium previously utilized deep sequencing methods to catalogue genetic diversity across African An. gambiae populations. We analyzed the complete datasets of 1142 individually sequenced mosquitoes through Microsoft Premonition's Bayesian mixture model based (BMM) metagenomics pipeline. All specimens were confirmed as either An. gambiae sensu stricto (s.s.) or An. coluzzii with a high degree of confidence ( > 98% identity to reference). Homo sapiens DNA was identified in all specimens indicating contamination may have occurred either at the time of specimen collection, preparation and/or sequencing. We found evidence of vertebrate hosts in 162 specimens. 59 specimens contained validated Plasmodium falciparum reads. Human hepatitis B and primate erythroparvovirus-1 viral sequences were identified in fifteen and three mosquito specimens, respectively. 478 of the 1,142 specimens were found to contain bacterial reads and bacteriophage-related contigs were detected in 27 specimens. This analysis demonstrates the capacity of metagenomic approaches to elucidate important vector-host-pathogen interactions of epidemiological significance.

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.

The pollen virome of wild plants and its association with variation in floral traits and land use.

Pollen is a unique vehicle for viral spread. Pollen-associated viruses hitchhike on or within pollen grains and are transported to other plants by pollinators. They are deposited on flowers and have a direct pathway into the plant and next generation via seeds. To discover the diversity of pollen-associated viruses and identify contributing landscape and floral features, we perform a species-level metagenomic survey of pollen from wild, visually asymptomatic plants, located in one of four regions in the United States of America varying in land use. We identify many known and novel pollen-associated viruses, half belonging to the Bromoviridae, Partitiviridae, and Secoviridae viral families, but many families are represented. Across the regions, species harbor more viruses when surrounded by less natural and more human-modified environments than the reverse, but we note that other region-level differences may also covary with this. When examining the novel connection between virus richness and floral traits, we find that species with multiple, bilaterally symmetric flowers and smaller, spikier pollen harbored more viruses than those with opposite traits. The association of viral diversity with floral traits highlights the need to incorporate plant-pollinator interactions as a driver of pollen-associated virus transport into the study of plant-viral interactions.

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.

Stability and detection of nucleic acid from viruses and hosts in controlled mosquito blood feeds.

Monitoring the presence and spread of pathogens in the environment is of critical importance. Rapid detection of infectious disease outbreaks and prediction of their spread can facilitate early responses of health agencies and reduce the severity of outbreaks. Current sampling methods are sorely limited by available personnel and throughput. For instance, xenosurveillance utilizes captured arthropod vectors, such as mosquitoes, as sampling tools to access blood from a wide variety of vertebrate hosts. Next generation sequencing (NGS) of nucleic acid from individual blooded mosquitoes can be used to identify mosquito and host species, and microorganisms including pathogens circulating within either host. However, there are practical challenges to collecting and processing mosquitoes for xenosurveillance, such as the rapid metabolization or decay of microorganisms within the mosquito midgut. This particularly affects pathogens that do not replicate in mosquitoes, preventing their detection by NGS or other methods. Accordingly, we performed a series of experiments to establish the windows of detection for DNA or RNA from human blood and/or viruses present in mosquito blood meals. Our results will contribute to the development of xenosurveillance techniques with respect to optimal timing of sample collection and NGS processing and will also aid trap design by demonstrating the stabilizing effect of temperature control on viral genome recovery from blood-fed mosquitoes.

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.

Coding-Complete Genome Sequence of a Pollen-Associated Virus Belonging to the Secoviridae Family Recovered from a Japanese Apricot (Prunus mume) Metagenome Data Set.

We report the coding-complete genome sequence of Japanese apricot pollen-associated secovirus 1 (JAPSV1), a virus belonging to the Secoviridae family, recovered from Japanese apricot (Prunus mume) pollen that is closely related to Peach leaf pitting-associated virus (PLPAV). This discovery adds to the number of known pollen-associated viruses.

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.

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.

Draft Genome Sequence of a Novel Rhabdovirus Isolated from Deinocerites Mosquitoes.

We present a draft genome of a novel rhabdovirus, called Grenada mosquito rhabdovirus 1 (GMRV1), with homology to Wuhan mosquito virus 9 (WMV9) (NCBI reference sequence NC_031303), isolated from Deinocerites mosquitoes. The genome has a length of 14,420 nucleotides and encodes five open reading frames.

Complete Genome Sequence of Pittsburgh Sewage-Associated Virus 1.

We present the complete genome sequence of a virus found in raw sewage collected in Pittsburgh, PA, USA. Pittsburgh sewage-associated virus 1 (PSAV1) encodes one large open reading frame with conserved domains typically found in the Picornavirales order of viruses. PSAV1 is closely related to Biomphalaria virus 2 (BV2).

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.

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.

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.