Detection of murine leukemia virus in the Epstein-Barr virus-positive human B-cell line JY, using a computational RNA-Seq-based exogenous agent detection pipeline, PARSES

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.

Characterization of a new case of XMLV (Bxv1) contamination in the human cell line Hep2 (clone 2B)

The use of misidentified cell lines contaminated by other cell lines and/or microorganisms has generated much confusion in the scientific literature. Detailed characterization of such contaminations is therefore crucial to avoid misinterpretation and ensure robustness and reproducibility of research. Here we use DNA-seq data produced in our lab to first confirm that the Hep2 (clone 2B) cell line (Sigma-Aldrich catalog number: 85011412-1VL) is indistinguishable from the HeLa cell line by mapping integrations of the human papillomavirus 18 (HPV18) at their expected loci on chromosome 8. We then show that the cell line is also contaminated by a xenotropic murine leukemia virus (XMLV) that is nearly identical to the mouse Bxv1 provirus and we characterize one Bxv1 provirus, located in the second intron of the pseudouridylate synthase 1 (PUS1) gene. Using an RNA-seq dataset, we confirm the high expression of the E6 and E7 HPV18 oncogenes, show that the entire Bxv1 genome is moderately expressed, and retrieve a Bxv1 splicing event favouring expression of the env gene. Hep2 (clone 2B) is the fourth human cell line so far known to be contaminated by the Bxv1 XMLV. This contamination has to be taken into account when using the cell line in future experiments.

Murine xenograft bioreactors for human immunopeptidome discovery

The study of peptides presented by MHC class I and class II molecules is limited by the need for relatively large cell numbers, especially when studying post-translationally modified or otherwise rare peptide species. To overcome this problem, we pose the hypothesis that human cells grown as xenografts in immunodeficient mice should produce equivalent immunopeptidomes as cultured cells. Comparing human cell lines grown either in vitro or as murine xenografts, we show that the immunopeptidome is substantially preserved. Numerous features are shared across both sample types, including peptides and proteins featured, length distributions, and HLA-binding motifs. Peptides well-represented in both groups were from more abundant proteins, or those with stronger predicted HLA binding affinities. Samples grown in vivo also recapitulated a similar phospho-immunopeptidome, with common sequences being those found at high copy number on the cell surface. These data indicate that xenografts are indeed a viable methodology for the production of cells for immunopeptidomic discovery.

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.

Use of a new RNA next generation sequencing approach for the specific detection of virus infection in cells

The utilization of the current combination of in vitro, in vivo and PCR assays for the identification of adventitious viruses in production cells has a limited range of detection. While Next Generation Sequencing (NGS) has a broader breadth of detection, it is unable to differentiate sequences from replicating viruses versus background inert sequences. In order to improve NGS specificity, we have designed a new NGS approach which targets subsets of viral RNAs only synthesized during cell infection. In order to evaluate the performance of this approach for detecting low levels of adventitious viruses, we selected two difficult virus/cell systems. This included B95-8 cells persistently infected by Human herpesvirus 4 (HHV-4) and serially diluted into HHV-4 negative Ramos cells and Madin-Darby bovine kidney cells with an early infection produced via a low dose of Bovine viral diarrhea virus. We demonstrated that the sensitivity of our RNA NGS approach was equivalent to targeted PCR with an increased specificity for the detection of viral infection. We were also able to identify a previously undetected Murine Leukemia Virus contaminant in Ramos cells. Based on these results, we conclude that this new RNA NGS approach is suitable for conducting viral safety evaluations of cells.

The Protein/Peptide Direct Virus Inactivation During Chromatographic Process: Developing Approaches

Virus clearance is required for pharmaceutical preparations derived from animal or human sources such as blood products, vaccines, recombinant proteins produced in mammalian cell lines, etc. High cost and substantial protein losses during virus inactivation are significant problems for protein/peptide manufacturing. The goal of this project was to develop a method to perform virus inactivation in a course of protein chromatographic purification. Another goal was to show that the chromatographic adsorbent can serve as reliable "sieva" for mechanical washing away of infecting viruses. Using chromatographic, photometric, IFA, and RT-PCR approaches, it was discovered that high temperature-depending dynamic capacity of adsorbent allowed to perform a virus inactivation directly in a chromatographic column by solvent/detergent treatment. The peptide/protein biological activity was completely preserved. Using this new approach enveloped and nonenveloped viruses were effectively removed protein preparation. In addition, it was shown that RT-PCR method demonstrates more precise and reproducible results and robust properties for assessment of virus reduction than virus titer followed by infectivity studies. Presented method allowed to obtain the factor of virus concentration decrease (FVD) values that were higher than those provided by known technologies and was sufficient for a full inactivation of viruses. The method is recommended to use in pharmaceutical industry.

Global transcript structure resolution of high gene density genomes through multi-platform data integration

Annotation of herpesvirus genomes has traditionally been undertaken through the detection of open reading frames and other genomic motifs, supplemented with sequencing of individual cDNAs. Second generation sequencing and high-density microarray studies have revealed vastly greater herpesvirus transcriptome complexity than is captured by existing annotation. The pervasive nature of overlapping transcription throughout herpesvirus genomes, however, poses substantial problems in resolving transcript structures using these methods alone. We present an approach that combines the unique attributes of Pacific Biosciences Iso-Seq long-read, Illumina short-read and deepCAGE (Cap Analysis of Gene Expression) sequencing to globally resolve polyadenylated isoform structures in replicating Epstein-Barr virus (EBV). Our method, Transcriptome Resolution through Integration of Multi-platform Data (TRIMD), identifies nearly 300 novel EBV transcripts, quadrupling the size of the annotated viral transcriptome. These findings illustrate an array of mechanisms through which EBV achieves functional diversity in its relatively small, compact genome including programmed alternative splicing (e.g. across the IR1 repeats), alternative promoter usage by LMP2 and other latency-associated transcripts, intergenic splicing at the BZLF2 locus, and antisense transcription and pervasive readthrough transcription throughout the genome.

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.

Latent Expression of the Epstein-Barr Virus (EBV)-Encoded Major Histocompatibility Complex Class I TAP Inhibitor, BNLF2a, in EBV-Positive Gastric Carcinomas

The Epstein-Barr virus (EBV) BNLF2a gene product provides immune evasion properties to infected cells through inhibition of transporter associated with antigen processing (TAP)-mediated transport of antigen peptides. Although BNLF2a is considered to be a lytic gene, we demonstrate that it is expressed in nearly half of the EBV-associated gastric carcinomas analyzed. Further, we show that BNLF2a expression is dissociated from lytic gene expression. BNLF2a is therefore expressed in this latency setting, potentially helping protect the infected tumor cells from immunosurveillance.

Frequent infection of human cancer xenografts with murine endogenous retroviruses in vivo

Infection of human cancer xenografts in mice with murine leukemia viruses (MLVs) is a long-standing observation, but the likelihood of infection in vivo and its biological consequences are poorly understood. We therefore conducted a prospective study in commonly used xenograft recipient strains. From BALB/c nude mice engrafted with MCF7 human mammary carcinoma cells, we isolated a virus that was virtually identical to Bxv1, a locus encoding replication-competent xenotropic MLV (XMLV). XMLV was detected in 9/17 (53%) independently isolated explants. XMLV was not found in primary leukemias or in THP1 leukemia cells grown in Bxv1-negative NSG (NOD/SCID/γCnull) mice, although MCF7 explants harbored replication-defective MLV proviruses. To assess the significance of infection for xenograft behavior in vivo, we examined changes in growth and global transcription in MCF7 and the highly susceptible Raji Burkitt lymphoma cell line chronically infected with XMLV. Raji cells showed a stronger transcriptional response that included up-regulation of chemokines and effectors of innate antiviral immunity. In conclusion, the risk of de novo XMLV infection of xenografts is high in Bxv1 positive mice, while infection can have positive or negative effects on xenograft growth potential with significant consequences for interpretation of many xenograft studies.

A comprehensive transcriptional portrait of human cancer cell lines

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

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.

Targeting xCT, a cystine-glutamate transporter induces apoptosis and tumor regression for KSHV/HIV-associated lymphoma

Kaposi’s sarcoma-associated herpesvirus (KSHV) is the etiological agent of primary effusion lymphoma (PEL), which represents a rapidly progressing malignancy arising in HIV-infected patients. Conventional chemotherapy for PEL treatment induces unwanted toxicity and is ineffective — PEL continues to portend nearly 100% mortality within a period of months, which requires novel therapeutic strategies. The amino acid transporter, xCT, is essential for the uptake of cystine required for intracellular glutathione (GSH) synthesis and for maintaining the intracellular redox balance. Inhibition of xCT induces growth arrest in a variety of cancer cells, although its role in virus-associated malignancies including PEL remains unclear. In the current study, we identify that xCT is expressed on the surface of patient-derived KSHV⁺ PEL cells, and targeting xCT induces caspase-dependent cell apoptosis. Further experiments demonstrate the underlying mechanisms including host and viral factors: reducing intracellular GSH while increasing reactive oxygen species (ROS), repressing cell-proliferation-related signaling, and inducing viral lytic genes. Using an immune-deficient xenograft model, we demonstrate that an xCT selective inhibitor, Sulfasalazine (SASP), prevents PEL tumor progression in vivo. Together, our data provide innovative and mechanistic insights into the role of xCT in PEL pathogenesis, and the framework for xCT-focused therapies for AIDS-related lymphoma in future.

Partial molecular cloning of the JHK retrovirus using gammaretrovirus consensus PCR primers

The JHK virus (JHKV) was previously described as a type C retrovirus that has some distinctive ultrastructural features and replicates constitutively in a human B-lymphoblastoid cell line, JHK-3. In order to facilitate the cloning of sequences from JHKV, a series of partially degenerate consensus retroviral PCR primers were created by a data-driven design approach based on an alignment of 14 diverse gammaretroviral genomes. These primers were used in the PCR amplification of purified JHK virion cDNA, and ana lysis of the resulting amplified sequence indicates that the JHKV is in the murine leukemia virus (MLV) family. The JHK sequence is nearly identical to the corresponding region of the Bxv-1 endogenous mouse retrovirus (GenBank accession AC115959) and distinct from XMRV. JHKV gag-specific amplification was demonstrated with nucleic acids from uncultivated, frozen, peripheral blood mononuclear cells (PBMCs) of the index patient, but not in PBMCs from nine healthy blood donors. Unlike earlier reports, in which MLV-like sequences were identified in human source material, which may have been due to murine contamination, budding retrovirions were demonstrated repeatedly by electron microscopy in uncultivated lymphocytes of the index patient that were morphologically identical in their development to the virions in the JHK-3 cells, and immunological evidence was obtained that the index patient produced IgG antibodies that bound to the budding viral particles in patient PBMCs and in the JHK-3 cells. These data indicate that the patient had been infected by JHKV, lending significance to the demonstration of JHKV amplicons in nucleic acids of the patient's PBMCs. In future studies, the PCR primer sets described herein may expand the detection of an amplifiable subset of viruses related to MLV.

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

Comprehensive virome analysis of RNA sequence (RNA-seq) data sets from 118 non-Hodgkin's B-cell lymphomas revealed a small subset that is positive for Epstein-Barr virus (EBV) or human herpesvirus 6B (HHV-6B), with one coinfection. EBV transcriptome analysis revealed expression of the latency genes RPMS1, LMP1, and LMP2, with one sample additionally showing a high level of early lytic expression and another sample showing a high level of EBNA2 expression. HHV-6B transcriptome analysis revealed that the majority of genes were transcribed.

GobyWeb: simplified management and analysis of gene expression and DNA methylation sequencing data

We present GobyWeb, a web-based system that facilitates the management and analysis of high-throughput sequencing (HTS) projects. The software provides integrated support for a broad set of HTS analyses and offers a simple plugin extension mechanism. Analyses currently supported include quantification of gene expression for messenger and small RNA sequencing, estimation of DNA methylation (i.e., reduced bisulfite sequencing and whole genome methyl-seq), or the detection of pathogens in sequenced data. In contrast to previous analysis pipelines developed for analysis of HTS data, GobyWeb requires significantly less storage space, runs analyses efficiently on a parallel grid, scales gracefully to process tens or hundreds of multi-gigabyte samples, yet can be used effectively by researchers who are comfortable using a web browser. We conducted performance evaluations of the software and found it to either outperform or have similar performance to analysis programs developed for specialized analyses of HTS data. We found that most biologists who took a one-hour GobyWeb training session were readily able to analyze RNA-Seq data with state of the art analysis tools. GobyWeb can be obtained at http://gobyweb.campagnelab.org and is freely available for non-commercial use. GobyWeb plugins are distributed in source code and licensed under the open source LGPL3 license to facilitate code inspection, reuse and independent extensions http://github.com/CampagneLaboratory/gobyweb2-plugins.

IMSA: integrated metagenomic sequence analysis for identification of exogenous reads in a host genomic background

Metagenomics, the study of microbial genomes within diverse environments, is a rapidly developing field. The identification of microbial sequences within a host organism enables the study of human intestinal, respiratory, and skin microbiota, and has allowed the identification of novel viruses in diseases such as Merkel cell carcinoma. There are few publicly available tools for metagenomic high throughput sequence analysis. We present Integrated Metagenomic Sequence Analysis (IMSA), a flexible, fast, and robust computational analysis pipeline that is available for public use. IMSA takes input sequence from high throughput datasets and uses a user-defined host database to filter out host sequence. IMSA then aligns the filtered reads to a user-defined universal database to characterize exogenous reads within the host background. IMSA assigns a score to each node of the taxonomy based on read frequency, and can output this as a taxonomy report suitable for cluster analysis or as a taxonomy map (TaxMap). IMSA also outputs the specific sequence reads assigned to a taxon of interest for downstream analysis. We demonstrate the use of IMSA to detect pathogens and normal flora within sequence data from a primary human cervical cancer carrying HPV16, a primary human cutaneous squamous cell carcinoma carrying HPV 16, the CaSki cell line carrying HPV16, and the HeLa cell line carrying HPV18.

Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant therapy

Epstein-Barr virus (EBV) is associated with roughly 10% of gastric carcinomas worldwide (EBVaGC). Although previous investigations provide a strong link between EBV and gastric carcinomas, these studies were performed using selected EBV gene probes. Using a cohort of gastric carcinoma RNA-seq data sets from The Cancer Genome Atlas (TCGA), we performed a quantitative and global assessment of EBV gene expression in gastric carcinomas and assessed EBV associated cellular pathway alterations. EBV transcripts were detected in 17% of samples but these samples varied significantly in EBV coverage depth. In four samples with the highest EBV coverage (hiEBVaGC – high EBV associated gastric carcinoma), transcripts from the BamHI A region comprised the majority of EBV reads. Expression of LMP2, and to a lesser extent, LMP1 were also observed as was evidence of abortive lytic replication. Analysis of cellular gene expression indicated significant immune cell infiltration and a predominant IFNG response in samples expressing high levels of EBV transcripts relative to samples expressing low or no EBV transcripts. Despite the apparent immune cell infiltration, high levels of the cytotoxic T-cell (CTL) and natural killer (NK) cell inhibitor, IDO1, was observed in the hiEBVaGCs samples suggesting an active tolerance inducing pathway in this subgroup. These results were confirmed in a separate cohort of 21 Vietnamese gastric carcinoma samples using qRT-PCR and on tissue samples using in situ hybridization and immunohistochemistry. Lastly, a panel of tumor suppressors and candidate oncogenes were expressed at lower levels in hiEBVaGC versus EBV-low and EBV-negative gastric cancers suggesting the direct regulation of tumor pathways by EBV.

Epstein-Barr virus (EBV) is detected in roughly 10% of gastric carcinoma (GC) cases worldwide. Despite a strong link between EBV and gastric carcinoma, the contribution of EBV to the tumor environment in EBV associated gastric carcinoma is unclear. We performed a global assessment of EBV and host cell gene expression in gastric carcinoma tumors from 71 patients to link EBV genes (and expression intensities) to cell and microenvironmental changes. In addition to the finding that EBV is associated with down-regulated tumor regulatory genes, this study revealed that samples with high levels of EBV gene expression (hiEBVaGCs) displayed elevated immune cell infiltration with high interferon-gamma (IFNG) expression compared to samples with low or no EBV gene expression. Despite this evidence of increased immune posturing, hiEBVaGC samples also showed elevated expression of the potent immune cell inhibitor, IDO1. This finding may partly explain the persistence of these virus associated tumors in the face of local immune cell concentration. Importantly, the small molecule IDO inhibitor, 1MT (1-methyl Tryptophan), has been shown to reverse the tolerance inducing effects of IDO1 in other tumors. We propose that stratification of gastric carcinomas into EBV-negative, EBV-low and EBV-high may provide indicator value for the use of IDO1 inhibitors as adjuvant therapies against hiEBVaGCs.

Limited RNA editing in exons of mouse liver and adipose

Several studies have investigated RNA-DNA differences (RDD), presumably due to RNA editing, with conflicting results. We report a rigorous analysis of RDD in exonic regions in mice, taking into account critical biases in RNA-Seq analysis. Using deep-sequenced F1 reciprocal inbred mice, we mapped 40 million RNA-Seq reads per liver sample and 180 million reads per adipose sample. We found 7300 apparent hepatic RDDs using a multiple-site mapping procedure, compared with 293 RDD found using a unique-site mapping procedure. After filtering for repeat sequence, splice junction proximity, undirectional strand, and extremity read bias, 63 RDD remained. In adipose tissue unique-site mapping identified 1667 RDD, and after applying the same four filters, 188 RDDs remained. In both tissues, the filtering procedure increased the proportion of canonical (A-to-I and C-to-U) editing events. The genomic DNA of 12 RDD sites among the potential 63 hepatic RDD was tested by Sanger sequencing, three of which proved to be due to unreferenced SNPs. We validated seven liver RDD with Sequenom technology, including two noncanonical, Gm5424 C-to-I(G) and Pisd I(G)-to-A RDD. Differences in diet, sex, or genetic background had very modest effects on RDD occurrence. Only a small number of apparent RDD sites overlapped between liver and adipose, indicating a high degree of tissue specificity. Our findings underscore the importance of properly filtering for bias in RNA-Seq investigations, including the necessity of confirming the DNA sequence to eliminate unreferenced SNPs. Based on our results, we conclude that RNA editing is likely limited to hundreds of events in exonic RNA in liver and adipose.

Vironome of Kaposi sarcoma associated herpesvirus-inflammatory cytokine syndrome in an AIDS patient reveals co-infection of human herpesvirus 8 and human herpesvirus 6A

KSHV inflammatory cytokine syndrome (KICS) is a newly described condition characterized by systemic illness as a result of systemic, lytic KSHV infection. We used Illumina sequencing to establish the DNA vironome of blood from such a patient. It identified concurrent high-level viremia of human herpesvirus (HHV) 8 and 6a. The HHV8 plasma viral load was 5,300,000 copies/ml, which is the highest reported to date; this despite less than five skin lesions and no HHV8 associated lymphoma. This is the first report of systemic HHV6a/KSHV co-infection in a patient. It is the first whole genome KSHV sequence to be determined directly from patient plasma rather than cultured or biopsied tumor material. This case supports KICS as a new clinical entity associated with KSHV.