Epstein-Barr virus and human herpesvirus 6 detection in a non-Hodgkin's diffuse large B-cell lymphoma cohort by using RNA sequencing

EBV DNA methylation profiles and its application in distinguishing nasopharyngeal carcinoma and nasal NK/T-cell lymphoma

BACKGROUND

As an oncovirus, EBV is associated with multiple cancers, including solid tumors and hematological malignancies. EBV methylation plays an important role in regulating tumor occurrence. However, the EBV methylation profiles in EBV-associated tumor tissues are poorly understood.

RESULTS

In this study, EBV methylation capture sequencing was conducted in several different tumor tissue samples, including NPC, EBVaGC, lung LELC and parotid LELC. Besides, EBV capture sequencing and following qMSP were performed on nasopharyngeal brushing samples from NPC and nasal NKTCL patients. Our results showed that the EBV genome among different types of tumors displayed specific methylation patterns. Among the four types of tumors from epithelial origin (NPC, EBVaGC, lung LELC and parotid LELC), the most significant differences were found between EBVaGC and the others. For example, in EBVaGC, all CpG sites within 1,44,189-1,45,136 bp of the EBV genome sequence on gene RPMS1 were hyper-methylated compared to the others. Differently, significant differences of EBV CpG sites, particularly those located on gene BILF2, were observed between NPC and nasal NKTCL patients in nasopharyngeal brushing samples. Further, the methylated level of BILF2 was further detected using qMSP, and a diagnostic model distinguishing NPC and nasal NKTCL was established. The AUC of the model was 0.9801 (95% CI 0.9524-1.0000), with the sensitivity and specificity of 98.81% (95% CI 93.63-99.94%) and 76.92% (95% CI 49.74-91.82%), respectively.

CONCLUSIONS

Our study reveals more clues for further understanding the pathogenesis of EBV, and provides a possibility for distinguishing EBV-related tumor by detecting specific EBV CpG sites.

Functional Implications of Epstein-Barr Virus Lytic Genes in Carcinogenesis

Epstein-Barr virus (EBV) is associated with a diverse range of tumors of both lymphoid and epithelial origin. Similar to other herpesviruses, EBV displays a bipartite life cycle consisting of latent and lytic phases. Current dogma indicates that the latent genes are key drivers in the pathogenesis of EBV-associated cancers, while the lytic genes are primarily responsible for viral transmission. In recent years, evidence has emerged to show that the EBV lytic phase also plays an important role in EBV tumorigenesis, and the expression of EBV lytic genes is frequently detected in tumor tissues and cell lines. The advent of next generation sequencing has allowed the comprehensive profiling of EBV gene expression, and this has revealed the consistent expression of several lytic genes across various types of EBV-associated cancers. In this review, we provide an overview of the functional implications of EBV lytic gene expression to the oncogenic process and discuss possible avenues for future investigations.

The Roles of the Virome in Cancer

Viral infections as well as changes in the composition of the intestinal microbiota and virome have been linked to cancer. Moreover, the success of cancer immunotherapy with checkpoint inhibitors has been correlated with the intestinal microbial composition of patients. The transfer of feces-which contain mainly bacteria and their viruses (phages)-from immunotherapy responders to non-responders, known as fecal microbiota transplantation (FMT), has been shown to be able to convert some non-responders to responders. Since phages may also increase the response to immunotherapy, for example by inducing T cells cross-reacting with cancer antigens, modulating phage populations may provide a new avenue to improve immunotherapy responsiveness. In this review, we summarize the current knowledge on the human virome and its links to cancer, and discuss the potential utility of bacteriophages in increasing the responder rate for cancer immunotherapy.

Pathogenic Role of Epstein-Barr Virus in Lung Cancers

Human oncogenic viruses account for at least 12% of total cancer cases worldwide. Epstein–Barr virus (EBV) is the first identified human oncogenic virus and it alone causes ~200,000 cancer cases and ~1.8% of total cancer-related death annually. Over the past 40 years, increasing lines of evidence have supported a causal link between EBV infection and a subgroup of lung cancers (LCs). In this article, we review the current understanding of the EBV-LC association and the etiological role of EBV in lung carcinogenesis. We also discuss the clinical impact of the knowledge gained from previous research, challenges, and future directions in this field. Given the high clinical relevance of EBV-LC association, there is an urgent need for further investigation on this topic.

Role of Virally Encoded Circular RNAs in the Pathogenicity of Human Oncogenic Viruses

Human oncogenic viruses are a group of important pathogens that etiologically contribute to at least 12% of total cancer cases in the world. As an emerging class of non-linear regulatory RNA molecules, circular RNAs (circRNAs) have gained increasing attention as a crucial player in the regulation of signaling pathways involved in viral infection and oncogenesis. With the assistance of current circRNA enrichment and detection technologies, numerous novel virally-encoded circRNAs (vcircRNAs) have been identified in the human oncogenic viruses, initiating an exciting new era of vcircRNA research. In this review, we discuss the current understanding of the roles of vcircRNAs in the respective viral infection cycles and in virus-associated pathogenesis.

Nanopore Sequencing and Its Clinical Applications

Nanopore sequencing is a method for determining the order and modifications of DNA/RNA nucleotides by detecting the electric current variations when DNA/RNA oligonucleotides pass through the nanometer-sized hole (nanopore). Nanopore-based DNA analysis techniques have been commercialized by Oxford Nanopore Technologies, NabSys, and Sequenom, and widely used in scientific researches recently including human genomics, cancer, metagenomics, plant sciences, etc., moreover, it also has potential applications in the field of healthcare due to its fast turn-around time, portable and real-time data analysis. Those features make it a promising technology for the point-of-care testing (POCT) and its potential clinical applications are briefly discussed in this chapter.

Sequencing facility and DNA source associated patterns of virus-mappable reads in whole-genome sequencing data

Numerous viral sequences have been reported in the whole-genome sequencing (WGS) data of human blood. However, it is not clear to what degree the virus-mappable reads represent true viral sequences rather than random-mapping or noise originating from sample preparation, sequencing processes, or other sources. Identification of patterns of virus-mappable reads may generate novel indicators for evaluating the origins of these viral sequences. We characterized paired-end unmapped reads and reads aligned to viral references in human WGS datasets, then compared patterns of the virus-mappable reads among DNA sources and sequencing facilities which produced these datasets. We then examined potential origins of the source- and facility-associated viral reads. The proportions of clean unmapped reads among the seven sequencing facilities were significantly different (P < 2 × 10-16). We identified 260,339 reads that were mappable to a total of 99 viral references in 2535 samples. The majority (86.7%) of these virus-mappable reads (corresponding to 47 viral references), which can be classified into four groups based on their distinct patterns, were strongly associated with sequencing facility or DNA source (adjusted P value <0.01). Possible origins of these reads include artificial sequences in library preparation, recombinant vectors in cell culture, and phages co-contaminated with their host bacteria. The sequencing facility-associated virus-mappable reads and patterns were repeatedly observed in other datasets produced in the same facilities. We have constructed an analytic framework and profiled the unmapped reads mappable to viral references. The results provide a new understanding of sequencing facility- and DNA source-associated batch effects in deep sequencing data and may facilitate improved bioinformatics filtering of reads.

Ubiquitin Modification of the Epstein-Barr Virus Immediate Early Transactivator Zta

The Epstein-Barr virus (EBV) immediate early transactivator Zta plays a key role in regulating the transition from latency to the lytic replication stages of EBV infection. Regulation of Zta is known to be controlled through a number of transcriptional and posttranscriptional events. Here, we show that Zta is targeted for ubiquitin modification and that this can occur in EBV-negative and in EBV-infected cells. Genetic studies show critical roles for both an amino-terminal region of Zta and the basic DNA binding domain of Zta in regulating Zta ubiquitination. Pulse-chase experiments demonstrate that the bulk population of Zta is relatively stable but that at least a subset of ubiquitinated Zta molecules are targeted for degradation in the cell. Mutation of four out of a total of nine lysine residues in Zta largely abrogates its ubiquitination, indicating that these are primary ubiquitination target sites. A Zta mutant carrying mutations at these four lysine residues (lysine 12, lysine 188, lysine 207, and lysine 219) cannot induce latently infected cells to produce and/or release infectious virions. Nevertheless, this mutant can induce early gene expression, suggesting a possible defect at the level of viral replication or later in the lytic cascade. As far as we know, this is the first study that has investigated the targeting of Zta by ubiquitination or its role in Zta function.IMPORTANCE Epstein-Barr virus (EBV) is a ubiquitous human pathogen and associated with various human diseases. EBV undergoes latency and lytic replication stages in its life cycle. The transition into the lytic replication stage, at which virus is produced, is mainly regulated by the viral gene product, Zta. Therefore, the regulation of Zta function becomes a central issue regarding viral biology and pathogenesis. Known modifications of Zta include phosphorylation and sumoylation. Here, we report the role of ubiquitination in regulating Zta function. We found that Zta is subjected to ubiquitination in both EBV-infected and EBV-negative cells. The ubiquitin modification targets 4 lysine residues on Zta, leading to both mono- and polyubiquitination of Zta. Ubiquitination of Zta affects the protein's stability and likely contributes to the progression of viral lytic replication. The function and fate of Zta may be determined by the specific lysine residue being modified.

Virome assembly and annotation in brain tissue based on next-generation sequencing

The glioblastoma multiforme (GBM) is one of the deadliest tumors. It has been speculated that virus plays a role in GBM but the evidences are controversy. Published researches are mainly limited to studies on the presence of human cytomegalovirus (HCMV) in GBM. No comprehensive assessment of the brain virome, the collection of viral material in the brain, based on recently sequenced data has been performed. Here, we characterized the virome from 111 GBM samples and 57 normal brain samples from eight projects in the SRA database by a tested and comprehensive assembly approach. The annotation of the assembled contigs showed that most viral sequences in the brain belong to the viral family Retroviridae. In some GBM samples, we also detected full genome sequence of a novel picornavirus recently discovered in invertebrates. Unlike previous reports, our study did not detect herpes virus such as HCMV in GBM from the data we used. However, some contigs that cannot be annotated with any known genes exhibited antibody epitopes in their sequences. These findings provide several avenues for potential cancer therapy: the newly discovered picornavirus could be a starting point to engineer novel oncolytic virus; and the exhibited antibody epitopes could be a source to explore potential drug targets for immune cancer therapy. By characterizing the virosphere in GBM and normal brain at a global level, the results from this study strengthen the link between GBM and viral infection which warrants the further investigation.

High-Throughput Sequence Analysis of Peripheral T-Cell Lymphomas Indicates Subtype-Specific Viral Gene Expression Patterns and Immune Cell Microenvironments

Certain peripheral T-cell lymphomas (PTCLs) have been associated with viral infection, particularly infection with Epstein-Barr virus (EBV). However, a comprehensive virome analysis across PTCLs has not previously been reported. Here we utilized published whole-transcriptome RNA sequencing (RNA-seq) data sets from seven different PTCL studies and new RNA-seq data from our laboratory to screen for virus association, to analyze viral gene expression, and to assess B- and T-cell receptor diversity paradigms across PTCL subtypes. In addition to identifying EBV in angioimmunoblastic T-cell lymphoma (AITL) and extranodal NK/T-cell lymphoma (ENKTL), two PTCL subtypes with well-established EBV associations, we also detected EBV in several cases of anaplastic large-cell lymphoma (ALCL), and we found evidence of infection by the oncogenic viruses Kaposi's sarcoma-associated herpesvirus and human T-cell leukemia virus type 1 in isolated PTCL cases. In AITLs, EBV gene expression analysis showed expression of immediate early, early, and late lytic genes, suggesting either low-level lytic gene expression or productive infection in a subset of EBV-infected B-lymphocyte stromal cells. Deconvolution of immune cell subpopulations demonstrated a greater B-cell signal in AITLs than in other PTCL subtypes, consistent with a larger role for B-cell support in the pathogenesis of AITL. Reconstructed T-cell receptor (TCR) and B-cell receptor (BCR) repertoires demonstrated increased BCR diversity in AITLs, consistent with a possible EBV-driven polyclonal response. These findings indicate potential alternative roles for EBV in PTCLs, in addition to the canonical oncogenic mechanisms associated with EBV latent infection. Our findings also suggest the involvement of other viruses in PTCL pathogenesis and demonstrate immunological alterations associated with these cancers.IMPORTANCE In this study, we utilized next-generation sequencing data from 7 different studies of peripheral T-cell lymphoma (PTCL) patient samples to globally assess viral associations, provide insights into the contributions of EBV gene expression to the tumor phenotype, and assess the unique roles of EBV in modulating the immune cell tumor microenvironment. These studies revealed potential roles for EBV replication genes in some PTCL subtypes, the possible role of additional human tumor viruses in rare cases of PTCLs, and a role for EBV in providing a unique immune microenvironmental niche in one subtype of PTCLs. Together, these studies provide new insights into the understudied role of tumor viruses in PTCLs.

Detection of Epstein-Barr Virus Infection in Non-Small Cell Lung Cancer

Previous investigations proposed a link between the Epstein-Barr virus (EBV) and lung cancer (LC), but the results are highly controversial largely due to the insufficient sample size and the inherent limitation of the traditional viral screening methods such as PCR. Unlike PCR, current next-generation sequencing (NGS) utilizes an unbiased method for the global assessment of all exogenous agents within a cancer sample with high sensitivity and specificity. In our current study, we aim to resolve this long-standing controversy by utilizing our unbiased NGS-based informatics approaches in conjunction with traditional molecular methods to investigate the role of EBV in a total of 1127 LC. In situ hybridization analysis of 110 LC and 10 normal lung samples detected EBV transcripts in 3 LC samples. Comprehensive virome analyses of RNA sequencing (RNA-seq) data sets from 1017 LC and 110 paired adjacent normal lung specimens revealed EBV transcripts in three lung squamous cell carcinoma and one lung adenocarcinoma samples. In the sample with the highest EBV coverage, transcripts from the BamHI A region accounted for the majority of EBV reads. Expression of EBNA-1, LMP-1 and LMP-2 was observed. A number of viral circular RNA candidates were also detected. Thus, we for the first time revealed a type II latency-like viral transcriptome in the setting of LC in vivo. The high-level expression of viral BamHI A transcripts in LC suggests a functional role of these transcripts, likely as long non-coding RNA. Analyses of cellular gene expression and stained tissue sections indicated an increased immune cell infiltration in the sample expressing high levels of EBV transcripts compared to samples expressing low EBV transcripts. Increased level of immune checkpoint blockade factors was also detected in the sample with higher levels of EBV transcripts, indicating an induced immune tolerance. Lastly, inhibition of immune pathways and activation of oncogenic pathways were detected in the sample with high EBV transcripts compared to the EBV-low LC indicating the direct regulation of cancer pathways by EBV. Taken together, our data support the notion that EBV likely plays a pathological role in a subset of LC.

The 2017 Network Tools and Applications in Biology (NETTAB) workshop: aims, topics and outcomes

The 17th International NETTAB workshop was held in Palermo, Italy, on October 16-18, 2017. The special topic for the meeting was "Methods, tools and platforms for Personalised Medicine in the Big Data Era", but the traditional topics of the meeting series were also included in the event. About 40 scientific contributions were presented, including four keynote lectures, five guest lectures, and many oral communications and posters. Also, three tutorials were organised before and after the workshop. Full papers from some of the best works presented in Palermo were submitted for this Supplement of BMC Bioinformatics. Here, we provide an overview of meeting aims and scope. We also shortly introduce selected papers that have been accepted for publication in this Supplement, for a complete presentation of the outcomes of the meeting.

From trash to treasure: detecting unexpected contamination in unmapped NGS data

Background

Next Generation Sequencing (NGS) experiments produce millions of short sequences that, mapped to a reference genome, provide biological insights at genomic, transcriptomic and epigenomic level. Typically the amount of reads that correctly maps to the reference genome ranges between 70% and 90%, leaving in some cases a consistent fraction of unmapped sequences. This ’misalignment’ can be ascribed to low quality bases or sequence differences between the sample reads and the reference genome. Investigating the source of the unmapped reads is definitely important to better assess the quality of the whole experiment and to check for possible downstream or upstream ’contamination’ from exogenous nucleic acids.

Results

Here we propose DecontaMiner, a tool to unravel the presence of contaminating sequences among the unmapped reads. It uses a subtraction approach to identify bacteria, fungi and viruses genome contamination. DecontaMiner generates several output files to track all the processed reads, and to provide a complete report of their characteristics. The good quality matches on microorganism genomes are counted and compared among samples. DecontaMiner builds an offline HTML page containing summary statistics and plots. The latter are obtained using the state-of-the-art D3 javascript libraries. DecontaMiner has been mainly used to detect contamination in human RNA-Seq data. The software is freely available at http://www-labgtp.na.icar.cnr.it/decontaminer.

Conclusions

DecontaMiner is a tool designed and developed to investigate the presence of contaminating sequences in unmapped NGS data. It can suggest the presence of contaminating organisms in sequenced samples, that might derive either from laboratory contamination or from their biological source, and in both cases can be considered as worthy of further investigation and experimental validation. The novelty of DecontaMiner is mainly represented by its easy integration with the standard procedures of NGS data analysis, while providing a complete, reliable, and automatic pipeline.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2684-x) contains supplementary material, which is available to authorized users.

Comparative Analysis of Gammaherpesvirus Circular RNA Repertoires: Conserved and Unique Viral Circular RNAs

Recent studies have identified circular RNAs (circRNAs) expressed from the Epstein-Barr virus (EBV) and Kaposi's sarcoma herpesvirus (KSHV) human DNA tumor viruses. To gain initial insights into the potential relevance of EBV circRNAs in virus biology and disease, we assessed the circRNAome of the interspecies homologue rhesus macaque lymphocryptovirus (rLCV) in a naturally occurring lymphoma from a simian immunodeficiency virus (SIV)-infected rhesus macaque. This analysis revealed rLCV orthologues of the latency-associated EBV circular RNAs circRPMS1_E4_E3a and circEBNA_U. Also identified in two samples displaying unusually high lytic gene expression was a novel rLCV circRNA that contains both conserved and rLCV-specific RPMS1 exons and whose backsplice junctions flank an rLCV lytic origin of replication (OriLyt). Analysis of a lytic infection model for the murid herpesvirus 68 (MHV68) rhadinovirus identified a cluster of circRNAs near an MHV68 lytic origin of replication, with the most abundant of these, circM11_ORF69, spanning the OriLyt. Lastly, analysis of KSHV latency and reactivation models revealed the latency associated circRNA originating from the vIRF4 gene as the predominant viral circRNA. Together, the results of this study broaden our appreciation for circRNA repertoires in the Lymphocryptovirus and Rhadinovirus genera of gammaherpesviruses and provide evolutionary support for viral circRNA functions in latency and viral replication.IMPORTANCE Infection with oncogenic gammaherpesviruses leads to long-term viral persistence through a dynamic interplay between the virus and the host immune system. Critical for remodeling of the host cell environment after the immune responses are viral noncoding RNAs that modulate host signaling pathways without attracting adaptive immune recognition. Despite the importance of noncoding RNAs in persistent infection, the circRNA class of noncoding RNAs has only recently been identified in gammaherpesviruses. Accordingly, their roles in virus infection and associated oncogenesis are unknown. Here we report evolutionary conservation of EBV-encoded circRNAs determined by assessing the circRNAome in rLCV-infected lymphomas from an SIV-infected rhesus macaque, and we report latent and lytic circRNAs from KSHV and MHV68. These experiments demonstrate utilization of the circular RNA class of RNAs across 4 members of the gammaherpesvirus subfamily, and they identify orthologues and potential homoplastic circRNAs, implying conserved circRNA functions in virus biology and associated malignancies.

RNA Sequencing of the In Vivo Human Herpesvirus 6B Transcriptome To Identify Targets for Clinical Assays Distinguishing between Latent and Active Infections

Human herpesvirus 6B (HHV-6B) DNA is frequently detected in human samples. Diagnostic assays distinguishing HHV-6B reactivation from latency are limited. This has impaired strategies to diagnose and treat HHV-6B-associated diseases. We used RNA sequencing to characterize and compare the HHV-6B transcriptome in multiple sample types, including (i) whole blood from hematopoietic cell transplant (HCT) recipients with and without HHV-6B plasma viremia, (ii) tumor tissue samples from subjects with large B cell lymphoma infected with HHV-6B, (iii) lymphoblastoid cell lines (LCLs) from subjects with inherited chromosomally integrated HHV-6B or latent infection with HHV-6B, and (iv) HHV-6B Z29 infected SupT1 CD4⁺ T cells. We demonstrated substantial overlap in the HHV-6B transcriptome observed in in vivo and in vitro samples, although there was variability in the breadth and quantity of gene expression across samples. The HHV-6B viral polymerase gene U38 was the only HHV-6B transcript detected in all next-generation RNA sequencing (RNA-seq) data sets and was one of the most highly expressed genes. We developed a novel reverse transcription-PCR assay targeting HHV-6B U38, which identified U38 mRNA in all tested whole-blood samples from patients with concurrent HHV-6B viremia. No HHV-6B U38 transcripts were detected by RNA-seq or reverse transcription-real-time quantitative PCR (RT-qPCR) in whole-blood samples from subjects without HHV-6B plasma detection or from latently infected LCLs. A RT-qPCR assay for HHV-6B U38 may be useful to identify lytic HHV-6B infection in nonplasma samples and samples from individuals with inherited chromosomally integrated HHV-6B. This study also demonstrates the feasibility of transcriptomic analyses for HCT recipients.IMPORTANCE Human herpesvirus 6B (HHV-6B) is a DNA virus that infects most children within the first few years of life. After primary infection, HHV-6B persists as a chronic, latent infection in many cell types. Additionally, HHV-6B can integrate into germ line chromosomes, resulting in individuals with viral DNA in every nucleated cell. Given that PCR to detect viral DNA is the mainstay for diagnosing HHV-6B infection, the characteristics of HHV-6B infection complicate efforts to distinguish between latent and active viral infection, particularly in immunocompromised patients who have frequent HHV-6B reactivation. In this study, we used RNA sequencing to characterize the HHV-6B gene expression profile in multiple sample types, and our findings identified evidence-based targets for diagnostic tests that distinguish between latent and active viral infection.

Profiling of lymphoma from formalin-fixed paraffin-embedded tissue

Molecular profiling of lymphoma samples has contributed enormously to our understanding of disease biology leading to detailed descriptions of diagnostic categories. These studies have also helped the field to recognize different subtypes of disease, different diseases that share similar cellular pathway perturbations, different immune responses, and different prognostic groups. While nearly all of these discoveries were made using unfixed, snap-frozen materials, with few exceptions, clinical biopsy materials are comprised of formalin-fixed and paraffin-embedded (FFPE) tissues. Here, we describe the impact of molecular profiling on the field of lymphoma, the challenges associated with using FFPE tissues for downstream molecular diagnostic testing, the various molecular profiling techniques, and also provide an example of the clinical application of a molecular profiling test of lymphoma using FFPE tissues.

Human Herpesvirus 6 and Malignancy: A Review

In order to determine the role of human herpesvirus 6 (HHV-6) in human disease, several confounding factors, including methods of detection, types of controls, and the ubiquitous nature of the virus, must be considered. This is particularly problematic in the case of cancer, in which rates of detection vary greatly among studies. To determine what part, if any, HHV-6 plays in oncogenesis, a review of the literature was performed. There is evidence that HHV-6 is present in certain types of cancer; however, detection of the virus within tumor cells is insufficient for assigning a direct role of HHV-6 in tumorigenesis. Findings supportive of a causal role for a virus in cancer include presence of the virus in a large proportion of cases, presence of the virus in most tumor cells, and virus-induced in-vitro cell transformation. HHV-6, if not directly oncogenic, may act as a contributory factor that indirectly enhances tumor cell growth, in some cases by cooperation with other viruses. Another possibility is that HHV-6 may merely be an opportunistic virus that thrives in the immunodeficient tumor microenvironment. Although many studies have been carried out, it is still premature to definitively implicate HHV-6 in several human cancers. In some instances, evidence suggests that HHV-6 may cooperate with other viruses, including EBV, HPV, and HHV-8, in the development of cancer, and HHV-6 may have a role in such conditions as nodular sclerosis Hodgkin lymphoma, gastrointestinal cancer, glial tumors, and oral cancers. However, further studies will be required to determine the exact contributions of HHV-6 to tumorigenesis.

Lymphoproliferative Syndromes Associated with Human Herpesvirus-6A and Human Herpesvirus-6B

Human herpesvirus 6A and 6B (HHV-6A and HHV-6B) have been noted since their discovery for their T-lymphotropism. Although it has proven difficult to determine the extent to which HHV-6A and HHV-6B are involved in the pathogenesis of many diseases, evidence suggests that primary infection and reactivation of both viruses may induce or contribute to the progression of several lymphoproliferative disorders, ranging from benign to malignant and including infectious mononucleosis-like illness, drug induced hypersensitivity syndrome/drug reaction with eosinophilia and systemic symptoms (DIHS/DRESS), and nodular sclerosis Hodgkin's lymphoma. Herein, we discuss the conditions associated with the lymphoproliferative capacity of HHV-6, as well as the potential mechanisms behind them. Continued exploration on this topic may add to our understanding of the interactions between HHV-6 and the immune system and may open the doors to more accurate diagnosis and treatment of certain lymphoproliferative disorders.

Tumour virology in the era of high-throughput genomics

With the advent of massively parallel sequencing, oncogenic viruses in tumours can now be detected in an unbiased and comprehensive manner. Additionally, new viruses or strains can be discovered based on sequence similarity with known viruses. Using this approach, the causative agent for Merkel cell carcinoma was identified. Subsequent studies using data from large collections of tumours have confirmed models built during decades of hypothesis-driven and low-throughput research, and a more detailed and comprehensive description of virus-tumour associations have emerged. Notably, large cohorts and high sequencing depth, in combination with newly developed bioinformatical techniques, have made it possible to rule out several suggested virus-tumour associations with a high degree of confidence. In this review we discuss possibilities, limitations and insights gained from using massively parallel sequencing to characterize tumours with viral content, with emphasis on detection of viral sequences and genomic integration events.This article is part of the themed issue 'Human oncogenic viruses'.

A case report of Rosai-Dorfman disease with kidney involvement

Rosai-Dorfman disease (RDD) is a rare histiocytic disorder of unclear etiology, which commonly presented with the enlargement of lymph nodes of the neck and the head. Here, we report an unusual case of 77-year-old male patient presenting with left kidney lesion with several small enlarged lymph nodes around the abdominal aorta. The diagnosis of left kidney cancer was suspected and the patient underwent left laparoscopic exploration and lymph node biopsy. Only saponification of the renal surrounding fat and enlargement of the left renal pedicle and 5 abdominal aortic lymph nodes were found; no kidney cancer was found. Surrenalectomy and lymphadenectomy dissection were then performed and the left kidney was retained. Intraoperative frozen and postoperative pathology indicates Rosai-Dorfman disease. RDD with kidney involvement is uncommon, and its x-ray imaging appearances are atypical, and often resemble kidney cancer leading to kidney loss. A systematic literature review was also performed to investigate the x-ray imaging and treatment features of this disease.

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.

Epstein-Barr virus lytic cycle involvement in diffuse large B cell lymphoma

Epstein-Barr virus (EBV)-mediated B cell transformation is achieved predominantly through the action of latent proteins, but recent evidence suggests that lytic EBV replication has also a certain pathogenic role in lymphomagenesis, at least in the early phases of cell transformation. Particularly, in diffuse large B cell lymphoma (DLBCL), the EBV lytic cycle is by and large unexplored, so to disclose lytic cell contribution to lymphomagenesis, our aim was to evaluate viral early and late lytic gene expression in relation to several immune response markers in a series of EBV+ DLBCL from Argentina. An unexpected number of cells expressed lytic transcripts, being transcribed at the BZLF1, BHRF1, and BLLF1 locus, by real-time quantitative polymerase chain reaction. This lytic antigen expression was confirmed by immunohistochemical staining for BMRF1 early lytic protein, and a positive correlation between lytic and latent genes was confirmed, revealing a close link between their expressions in EBV+ DLBCL pathogenesis. Remarkably, BZLF1 displayed a negative correlation with CD4 cell counts, and this could be in part justified by the restriction of antigen presentation previously reported. The direct correlation for the late lytic gene BLLF1 and IFNγ in this series could represent a specific response directed towards this antigen. Interleukin 10 transcripts also displayed a positive correlation with lytic expression, indicating that regulatory mechanisms could be also involved on EBV-associated DLBCL pathogenesis in our series. Complete lytic reactivation in EBV-positive tumours could potentially kill EBV-positive malignant cells, providing a tool to promote tumour cell killing mediated by EBV as a complementary treatment strategy.

Complete Unique Genome Sequence, Expression Profile, and Salivary Gland Tissue Tropism of the Herpesvirus 7 Homolog in Pigtailed Macaques

Human herpesvirus 6A (HHV-6A), HHV-6B, and HHV-7 are classified as roseoloviruses and are highly prevalent in the human population. Roseolovirus reactivation in an immunocompromised host can cause severe pathologies. While the pathogenic potential of HHV-7 is unclear, it can reactivate HHV-6 from latency and thus contributes to severe pathological conditions associated with HHV-6. Because of the ubiquitous nature of roseoloviruses, their roles in such interactions and the resulting pathological consequences have been difficult to study. Furthermore, the lack of a relevant animal model for HHV-7 infection has hindered a better understanding of its contribution to roseolovirus-associated diseases. Using next-generation sequencing analysis, we characterized the unique genome of an uncultured novel pigtailed macaque roseolovirus. Detailed genomic analysis revealed the presence of gene homologs to all 84 known HHV-7 open reading frames. Phylogenetic analysis confirmed that the virus is a macaque homolog of HHV-7, which we have provisionally named Macaca nemestrina herpesvirus 7 (MneHV7). Using high-throughput RNA sequencing, we observed that the salivary gland tissue samples from nine different macaques had distinct MneHV7 gene expression patterns and that the overall number of viral transcripts correlated with viral loads in parotid gland tissue and saliva. Immunohistochemistry staining confirmed that, like HHV-7, MneHV7 exhibits a natural tropism for salivary gland ductal cells. We also observed staining for MneHV7 in peripheral nerve ganglia present in salivary gland tissues, suggesting that HHV-7 may also have a tropism for the peripheral nervous system. Our data demonstrate that MneHV7-infected macaques represent a relevant animal model that may help clarify the causality between roseolovirus reactivation and diseases.

IMPORTANCE Human herpesvirus 6A (HHV-6A), HHV-6B, and HHV-7 are classified as roseoloviruses. We have recently discovered that pigtailed macaques are naturally infected with viral homologs of HHV-6 and HHV-7, which we provisionally named MneHV6 and MneHV7, respectively. In this study, we confirm that MneHV7 is genetically and biologically similar to its human counterpart, HHV-7. We determined the complete unique MneHV7 genome sequence and provide a comprehensive annotation of all genes. We also characterized viral transcription profiles in salivary glands from naturally infected macaques. We show that broad transcriptional activity across most of the viral genome is associated with high viral loads in infected parotid glands and that late viral protein expression is detected in salivary duct cells and peripheral nerve ganglia. Our study provides new insights into the natural behavior of an extremely prevalent virus and establishes a basis for subsequent investigations of the mechanisms that cause HHV-7 reactivation and associated disease.

A comprehensive next generation sequencing-based virome assessment in brain tissue suggests no major virus - tumor association

Next generation sequencing (NGS) can globally interrogate the genetic composition of biological samples in an unbiased yet sensitive manner. The objective of this study was to utilize the capabilities of NGS to investigate the reported association between glioblastoma multiforme (GBM) and human cytomegalovirus (HCMV). A large-scale comprehensive virome assessment was performed on publicly available sequencing datasets from the Cancer Genome Atlas (TCGA), including RNA-seq datasets from primary GBM (n = 157), recurrent GBM (n = 13), low-grade gliomas (n = 514), recurrent low-grade gliomas (n = 17), and normal brain (n = 5), and whole genome sequencing (WGS) datasets from primary GBM (n = 51), recurrent GBM (n = 10), and normal matched blood samples (n = 20). In addition, RNA-seq datasets from MRI-guided biopsies (n = 92) and glioma stem-like cell cultures (n = 9) were analyzed. Sixty-four DNA-seq datasets from 11 meningiomas and their corresponding blood control samples were also analyzed. Finally, three primary GBM tissue samples were obtained, sequenced using RNA-seq, and analyzed. After in-depth analysis, the most robust virus findings were the detection of papillomavirus (HPV) and hepatitis B reads in the occasional LGG sample (4 samples and 1 sample, respectively). In addition, low numbers of virus reads were detected in several datasets but detailed investigation of these reads suggest that these findings likely represent artifacts or non-pathological infections. For example, all of the sporadic low level HCMV reads were found to map to the immediate early promoter intimating that they likely originated from laboratory expression vector contamination. Despite the detection of low numbers of Epstein-Barr virus reads in some samples, these likely originated from infiltrating B-cells. Finally, human herpesvirus 6 and 7 aligned viral reads were identified in all DNA-seq and a few RNA-seq datasets but detailed analysis demonstrated that these were likely derived from the homologous human telomeric-like repeats. Other low abundance viral reads were detected in some samples but for most viruses, the reads likely represent artifacts or incidental infections. This analysis argues against associations between most known viruses and GBM or mengiomas. Nevertheless, there may be a low percentage association between HPV and/or hepatitis B and LGGs.

Age and cellular composition influence overall survival in a collective of non-immunocompromised patients with EBV-positive diffuse large B-cell lymphoma from a German lymphoma center

We investigated 41 diffuse large B-cell lymphomas (DLBCL) diagnosed at one center harboring ≥50% of latently Epstein-Barr virus (EBV)-infected neoplastic cells occurring in 34 patients aged ≥50 years and in 7 patients younger than 50 years in the absence of any known immunodeficiency for the expression patterns of EBV latent and immediate-early proteins, for the differentiation stage of the neoplastic cells, the presence of cytogenetic alterations and a possible co-infection with the human herpes virus (HHV)-8. Here, we show that EBV-positive DLBCLs rarely arise from naïve and more frequently from post-germinal center B-cells that often contain crippling immunoglobulin gene mutations. Most of the lymphomas did not exhibit breaks in the BCL2, BCL6, and MYC genes and none of the cases investigated contained HHV-8 sequences. Patients aged <50 years performed better than older ones while in patients aged ≥50 years only the cellular composition had an impact on overall survival.

Possible Human Papillomavirus 38 Contamination of Endometrial Cancer RNA Sequencing Samples in The Cancer Genome Atlas Database

Viruses are causally associated with a number of human malignancies. In this study, we sought to identify new virus-cancer associations by searching RNA sequencing data sets from >2,000 patients, encompassing 21 cancers from The Cancer Genome Atlas (TCGA), for the presence of viral sequences. In agreement with previous studies, we found human papillomavirus 16 (HPV16) and HPV18 in oropharyngeal cancer and hepatitis B and C viruses in liver cancer. Unexpectedly, however, we found HPV38, a cutaneous form of HPV associated with skin cancer, in 32 of 168 samples from endometrial cancer. In 12 of the HPV38-positive (HPV38(+)) samples, we observed at least one paired read that mapped to both human and HPV38 genomes, indicative of viral integration into the host DNA, something not previously demonstrated for HPV38. The expression levels of HPV38 transcripts were relatively low, and all 32 HPV38(+) samples belonged to the same experimental batch of 40 samples, whereas none of the other 128 endometrial carcinoma samples were HPV38(+), raising doubts about the significance of the HPV38 association. Moreover, the HPV38(+) samples contained the same 10 novel single nucleotide variations (SNVs), leading us to hypothesize that one patient was infected with this new isolate of HPV38, which was integrated into his/her genome and may have cross-contaminated other TCGA samples within batch 228. Based on our analysis, we propose guidelines to examine the batch effect, virus expression level, and SNVs as part of next-generation sequencing (NGS) data analysis for evaluating the significance of viral/pathogen sequences in clinical samples.

IMPORTANCE

High-throughput RNA sequencing (RNA-Seq), followed by computational analysis, has vastly accelerated the identification of viral and other pathogenic sequences in clinical samples, but cross-contamination during the processing of the samples remain a major problem that can lead to erroneous conclusions. We found HPV38 sequences specifically present in RNA-Seq samples from endometrial cancer patients from TCGA, a virus not previously associated with this type of cancer. However, multiple lines of evidence suggest possible cross-contamination in these samples, which were processed together in the same batch. Despite this potential cross-contamination, our data indicate that we have detected a new isolate of HPV38 that appears to be integrated into the human genome. We also provide general guidelines for computational detection and interpretation of pathogen-disease associations.

No significant viral transcription detected in whole breast cancer transcriptomes

Background

Studies evaluating the presence of viral sequences in breast cancer (BC), including various strains of human papillomavirus and human herpes virus, have yielded conflicting results. Most were based on RT-PCR and in situ hybridization.

Methods

In this report we searched for expressed viral sequences in 58 BC transcriptomes using five distinct in silico methods. In addition, we complemented our RNA sequencing results with exome sequencing, PCR and immunohistochemistry (IHC) analyses. A control sample was used to test our in silico methods.

Results

All of the computational methods correctly detected viral sequences in the control sample. We identified a small number of viral sequences belonging to human herpesvirus 4 and 6 and Merkel cell polyomavirus. The extremely low expression levels—two orders of magnitude lower than in a typical hepatitis B virus infection in hepatocellular carcinoma—did not suggest active infections. The presence of viral elements was confirmed in sample-matched exome sequences, but could not be confirmed by PCR or IHC.

Conclusions

Our results show that no viral sequences are expressed in significant amounts in the BC investigated. The presence of non-transcribed viral DNA cannot be excluded.

Epstein-Barr virus-positive diffuse large B-cell lymphoma association is not only restricted to elderly patients

Diffuse large B-cell lymphoma (DLBCL), the most common group of malignant lymphomas, account for 30% of adult non-Hodgkin lymphomas. The 2008 World Health Organization (WHO) classification included a new entity, Epstein-Barr virus (EBV)+ DLBCL of the elderly, affecting patients aged 50 years or older. However, some reports of younger EBV+ DLBCL cases, without evidence of underlying immunosuppression, can be found. The role of EBV in tumor microenvironment composition in DLBCL is still not well understood. Our aim was to assess EBV presence and latency pattern as well as tumor T-cell population in an adult DLBCL series of Argentina. The study was conducted on biopsies from 75 DLBCL patients. EBERs expression was performed by in situ hybridization, while EBV gene expression was analyzed using real-time polymerase chain reaction. LMP1, LMP2A, EBNA2, EBNA3A, CD4, CD8 and Foxp3 expression was assessed by immunohistochemistry. Nine percent of cases showed EBV expression, with similar frequency among patients younger than 50 years and 50 years or older (13% and 8%, respectively). T-cell subsets were not altered by EBV presence. Latency type II was the most frequently observed, together with lytic gene expression in EBV+ DLBCL, with ≥20% of EBERs+ cells. These findings suggest that EBV+ DLBCL in our series was similar to the previously described in Asia and Latin-America, displaying latency II or III expression profile and no age-specific characteristics. Finally, EBV+ DLBCL may be an entity that is not only restricted to patients who are older than 50 years of age, in consequence the age cutoff revision may be a current goal.

Roseomics: a blank slate

Recent technological advances have led to an explosion in the system-wide profiling of biological processes in the study of herpesvirus biology, herein referred to as '-omics'. In many cases these approaches have revealed novel virus-induced changes to host cell biology that can be targeted with new antiviral therapeutics. Despite these successes, -omics approaches are not widely applied in the study of roseoloviruses. Here we describe examples of how -omics studies have shaped our understanding of herpesvirus biology, and discuss how these approaches might be used to identify host and viral factors that mediate roseolovirus pathogenesis.

Microbial contamination in next generation sequencing: implications for sequence-based analysis of clinical samples

The high level of accuracy and sensitivity of next generation sequencing for quantifying genetic material across organismal boundaries gives it tremendous potential for pathogen discovery and diagnosis in human disease. Despite this promise, substantial bacterial contamination is routinely found in existing human-derived RNA-seq datasets that likely arises from environmental sources. This raises the need for stringent sequencing and analysis protocols for studies investigating sequence-based microbial signatures in clinical samples.

High-throughput RNA sequencing-based virome analysis of 50 lymphoma cell lines from the Cancer Cell Line Encyclopedia project

Using high-throughput RNA sequencing data from 50 common lymphoma cell culture models from the Cancer Cell Line Encyclopedia project, we performed an unbiased global interrogation for the presence of a panel of 740 viruses and strains known to infect human and other mammalian cells. This led to the findings of previously identified infections by Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV), and human T-lymphotropic virus type 1 (HTLV-1). In addition, we also found a previously unreported infection of one cell line (DEL) with a murine leukemia virus. High expression of murine leukemia virus (MuLV) transcripts was observed in DEL cells, and we identified four transcriptionally active integration sites, one being in the TNFRSF6B gene. We also found low levels of MuLV reads in a number of other cell lines and provided evidence suggesting cross-contamination during sequencing. Analysis of HTLV-1 integrations in two cell lines, HuT 102 and MJ, identified 14 and 66 transcriptionally active integration sites with potentially activating integrations in immune regulatory genes, including interleukin-15 (IL-15), IL-6ST, STAT5B, HIVEP1, and IL-9R. Although KSHV and EBV do not typically integrate into the genome, we investigated a previously identified integration of EBV into the BACH2 locus in Raji cells. This analysis identified a BACH2 disruption mechanism involving splice donor sequestration. Through viral gene expression analysis, we detected expression of stable intronic RNAs from the EBV BamHI W repeats that may be part of long transcripts spanning the repeat region. We also observed transcripts at the EBV vIL-10 locus exclusively in the Hodgkin's lymphoma cell line, Hs 611.T, the expression of which were uncoupled from other lytic genes. Assessment of the KSHV viral transcriptome in BCP-1 cells showed expression of the viral immune regulators, K2/vIL-6, K4/vIL-8-like vCCL1, and K5/E2-ubiquitin ligase 1 that was significantly higher than expression of the latency-associated nuclear antigen. Together, this investigation sheds light into the virus composition across these lymphoma model systems and provides insights into common viral mechanistic principles.

IMPORTANCE

Viruses cause cancer in humans. In lymphomas the Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV) and human T-lymphotropic virus type 1 are major contributors to oncogenesis. We assessed virus-host interactions using a high throughput sequencing method that facilitates the discovery of new virus-host associations and the investigation into how the viruses alter their host environment. We found a previously unknown murine leukemia virus infection in one cell line. We identified cellular genes, including cytokine regulators, that are disrupted by virus integration, and we determined mechanisms through which virus integration causes deregulation of cellular gene expression. Investigation into the KSHV transcriptome in the BCP-1 cell line revealed high-level expression of immune signaling genes. EBV transcriptome analysis showed expression of vIL-10 transcripts in a Hodgkin's lymphoma that was uncoupled from lytic genes. These findings illustrate unique mechanisms of viral gene regulation and to the importance of virus-mediated host immune signaling in lymphomas.

Comprehensive high-throughput RNA sequencing analysis reveals contamination of multiple nasopharyngeal carcinoma cell lines with HeLa cell genomes

In an attempt to explore infectious agents associated with nasopharyngeal carcinomas (NPCs), we employed our high-throughput RNA sequencing (RNA-seq) analysis pipeline, RNA CoMPASS, to investigate the presence of ectopic organisms within a number of NPC cell lines commonly used by NPC and Epstein-Barr virus (EBV) researchers. Sequencing data sets from both CNE1 and HONE1 were found to contain reads for human papillomavirus 18 (HPV-18). Subsequent real-time reverse transcription-PCR (RT-PCR) analysis on a panel of NPC cell lines identified HPV-18 in CNE1 and HONE1 as well as three additional NPC cell lines (CNE2, AdAH, and NPC-KT). Further analysis of the chromosomal integration arrangement of HPV-18 in NPCs revealed patterns identical to those observed in HeLa cells. Clustering based on human single nucleotide variation (SNV) analysis of two separate HeLa cell lines and several NPC cell lines demonstrated two distinct clusters with CNE1, as well as HONE1 clustering with the two HeLa cell lines. In addition, duplex-PCR-based genotyping showed that CNE1, CNE2, and HONE1 do not have a HeLa cell-specific L1 retrotransposon insertion, suggesting that these three HPV-18(+) NPC lines are likely products of a somatic hybridization with HeLa cells, which is also consistent with our RNA-seq-based gene level SNV analysis. Taking all of these findings together, we conclude that a widespread HeLa contamination may exist in many NPC cell lines, and authentication of these cell lines is recommended. Finally, we provide a proof of concept for the utility of an RNA-seq-based approach for cell authentication.

IMPORTANCE

Nasopharyngeal carcinoma (NPC) cell lines are important model systems for analyzing the complex life cycle and pathogenesis of Epstein-Barr virus (EBV). Using an RNA-seq-based approach, we found HeLa cell contamination in several NPC cell lines that are commonly used in the EBV and related fields. Our data support the notion that contamination resulted from somatic hybridization with HeLa cells, likely occurring at the point of cell line establishment. Given the rarity of NPCs, the long history of NPC cell lines, and the lack of rigorous cell line authentication, it is likely that the actual prevalence and impact of HeLa cell contamination on the EBV field might be greater. We therefore recommend cell line authentication prior to performing experiments using NPC cell lines to avoid inaccurate conclusions. The novel RNA-seq-based cell authentication approach reported here can serve as a comprehensive method for validating cell lines.