Michael J. Strong

Whole-Genome Sequencing of the Akata and Mutu Epstein-Barr Virus Strains

ABSTRACT Using a simple viral genome enrichment approach, we report the de novo assembly of the Akata and Mutu Epstein-Barr virus (EBV) genomes from a single lane of next-generation sequencing (NGS) reads. The Akata and Mutu viral genomes are type I EBV strains of approximately 171 kb in length. Evidence for genome heterogeneity was found for the Akata but not for the Mutu strain. A comparative analysis of Akata with another four completely sequenced EBV strains, B95-8/Raji, AG876, Mutu, and GD1, demonstrated that the Akata strain is most closely related to the GD1 strain and exhibits the greatest divergence from the type II strain, AG876. A global comparison of latent and lytic gene sequences showed that the four latency genes, EBNA2, EBNA3A, EBNA3B, and EBNA3C, are uniquely defining of type I and type II strain differences. Within type I strains, LMP1, the latency gene, is among the most divergent of all EBV genes, with three insertion or deletion loci in its CTAR2 and CTAR3 signaling domains. Analysis of the BHLF1 and LF3 genes showed that the reading frames identified in the B95-8/Raji genome are not conserved in Akata (or Mutu, for BHLF1), suggesting a primarily non-protein-coding function in EBVs life cycle. The Akata and Mutu viral-genome sequences should be a useful resource for homology-based functional prediction and for molecular studies, such as PCR, RNA-seq, recombineering, and transcriptome studies. As an illustration, we identified novel RNA-editing events in ebv-miR-BART6 antisense transcripts using the Akata and Mutu reference genomes.

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.

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.

Galectin-9 Functionally Impairs Natural Killer Cells in Humans and Mice

ABSTRACT Galectin-9 is a pleiotropic immune modulator affecting numerous cell types of innate and adaptive immunity. Patients with chronic infection with either hepatitis C virus (HCV) or HIV have elevated circulating levels. Limited data exist on the regulation of natural killer (NK) cell function through interaction with galectin-9. We found that galectin-9 ligation downregulates multiple immune-activating genes, including eight involved in the NK cell-mediated cytotoxicity pathway, impairs lymphokine-activated killing, and decreases the proportion of gamma interferon (IFN-γ)-producing NK cells that had been stimulated with interleukin-12 (IL-12)/IL-15. We demonstrate that the transcriptional and functional changes induced by galectin-9 are independent of Tim-3. Consistent with these results for humans, we find that the genetic absence of galectin-9 in mice is associated with greater IFN-γ production by NK cells and enhanced degranulation. We also show that in the setting of a short-term (4-day) murine cytomegalovirus infection, terminally differentiated NKs accumulate in the livers of galectin-9 knockout mice, and that hepatic NKs spontaneously produce significantly more IFN-γ in this setting. Taken together, our results indicate that galectin-9 engagement impairs the function of NK cells, including cytotoxicity and cytokine production.

Comprehensive High-Throughput RNA Sequencing Analysis Reveals Contamination of Multiple Nasopharyngeal Carcinoma Cell Lines with HeLa Cell Genomes

ABSTRACT 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.

High-Throughput RNA Sequencing-Based Virome Analysis of 50 Lymphoma Cell Lines from the Cancer Cell Line Encyclopedia Project

ABSTRACT 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), Kaposis 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 Hodgkins 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), Kaposis 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 Hodgkins 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.

Global Bidirectional Transcription of the Epstein-Barr Virus Genome during Reactivation

ABSTRACT Epstein-Barr virus (EBV) reactivation involves the ordered induction of approximately 90 viral genes that participate in the generation of infectious virions. Using strand-specific RNA-seq to assess the EBV transcriptome during reactivation, we found extensive bidirectional transcription extending across nearly the entire genome. In contrast, only 4% of the EBV genome is currently bidirectionally annotated. Most of the newly identified transcribed regions show little evidence of coding potential, supporting noncoding roles for most of these RNAs. Based on previous cellular long noncoding RNA size calculations, we estimate that there are likely hundreds more EBV genes expressed during reactivation than was previously known. Limited 5′ and 3′ rapid amplification of cDNA ends (RACE) experiments and findings of novel splicing events by RNA-seq suggest that the complexity of the viral genome during reactivation may be even greater. Further analysis of antisense transcripts at some of the EBV latency gene loci showed that they are “late” genes, they are nuclear, and they tend to localize in areas of the nucleus where others find newly synthesized viral genomes. This raises the possibility that these transcripts perform functions such as new genome processing, stabilization, organization, etc. The finding of a significantly more complex EBV transcriptome during reactivation changes our view of the viral production process from one that is facilitated and regulated almost entirely by previously identified viral proteins to a process that also involves the contribution of a wide array of virus encoded noncoding RNAs. IMPORTANCE Epstein-Barr virus (EBV) is a herpesvirus that infects the majority of the worlds population, in rare cases causing serious disease such as lymphoma and gastric carcinoma. Using strand-specific RNA-seq, we have studied viral gene expression during EBV reactivation and have discovered hundreds more viral transcripts than were previously known. The finding of alternative splicing and the prevalence of overlapping transcripts indicate additional complexity. Most newly identified transcribed regions do not encode proteins but instead likely function as noncoding RNA molecules which could participate in regulating gene expression, gene splicing or even activities such as viral genome processing. These findings broaden the scope of what we need to consider to understand the viral manufacturing process. As more detailed studies are undertaken they will likely change the way we view this process as a whole.

RNA CoMPASS: A Dual Approach for Pathogen and Host Transcriptome Analysis of RNA-Seq Datasets

High-throughput RNA sequencing (RNA-seq) has become an instrumental assay for the analysis of multiple aspects of an organisms transcriptome. Further, the analysis of a biological specimens associated microbiome can also be performed using RNA-seq data and this application is gaining interest in the scientific community. There are many existing bioinformatics tools designed for analysis and visualization of transcriptome data. Despite the availability of an array of next generation sequencing (NGS) analysis tools, the analysis of RNA-seq data sets poses a challenge for many biomedical researchers who are not familiar with command-line tools. Here we present RNA CoMPASS, a comprehensive RNA-seq analysis pipeline for the simultaneous analysis of transcriptomes and metatranscriptomes from diverse biological specimens. RNA CoMPASS leverages existing tools and parallel computing technology to facilitate the analysis of even very large datasets. RNA CoMPASS has a web-based graphical user interface with intrinsic queuing to control a distributed computational pipeline. RNA CoMPASS was evaluated by analyzing RNA-seq data sets from 45 B-cell samples. Twenty-two of these samples were derived from lymphoblastoid cell lines (LCLs) generated by the infection of naïve B-cells with the Epstein Barr virus (EBV), while another 23 samples were derived from Burkitts lymphomas (BL), some of which arose in part through infection with EBV. Appropriately, RNA CoMPASS identified EBV in all LCLs and in a fraction of the BLs. Cluster analysis of the human transcriptome component of the RNA CoMPASS output clearly separated the BLs (which have a germinal center-like phenotype) from the LCLs (which have a blast-like phenotype) with evidence of activated MYC signaling and lower interferon and NF-kB signaling in the BLs. Together, this analysis illustrates the utility of RNA CoMPASS in the simultaneous analysis of transcriptome and metatranscriptome data. RNA CoMPASS is freely available at http://rnacompass.sourceforge.net/.

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

ABSTRACT Many cell lines commonly used for biological studies have been found to harbor exogenous agents such as the human tumor viruses Epstein-Barr virus (EBV) and human papillomavirus. Nevertheless, broad-based, unbiased approaches to globally assess the presence of ectopic organisms within cell model systems have not previously been available. We reasoned that high-throughput sequencing should provide unparalleled insights into the microbiomes of tissue culture cell systems. Here we have used our RNA-seq analysis pipeline, PARSES (Pipeline for Analysis of RNA-Seq Exogenous Sequences), to investigate the presence of ectopic organisms within two EBV-positive B-cell lines commonly used by EBV researchers. Sequencing data sets from both the Akata and JY B-cell lines were found to contain reads for EBV, and the JY data set was found to also contain reads from the murine leukemia virus (MuLV). Further investigation revealed that MuLV transcription in JY cells is highly active. We also identified a number of MuLV alternative splicing events, and we uncovered evidence of APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G)-dependent DNA editing. Finally, reverse transcription-PCR analysis showed the presence of MuLV in three other human B-cell lines (DG75, Ramos, and P3HR1 Cl.13) commonly used by investigators in the Epstein-Barr virus field. We believe that a thorough examination of tissue culture microbiomes using RNA-seq/PARSES-like approaches is critical for the appropriate utilization of these systems in biological studies.

Effects of the endocrine-disrupting chemical DDT on self-renewal and differentiation of human mesenchymal stem cells.

Background: Although the global use of the endocrine-disrupting chemical DDT has decreased, its persistence in the environment has resulted in continued human exposure. Accumulating evidence suggests that DDT exposure has long-term adverse effects on development, yet the impact on growth and differentiation of adult stem cells remains unclear. Objectives: Human mesenchymal stem cells (MSCs) exposed to DDT were used to evaluate the impact on stem cell biology. Methods: We assessed DDT-treated MSCs for self-renewal, proliferation, and differentiation potential. Whole genome RNA sequencing was performed to assess gene expression in DDT-treated MSCs. Results: MSCs exposed to DDT formed fewer colonies, suggesting a reduction in self-renewal potential. DDT enhanced both adipogenic and osteogenic differentiation, which was confirmed by increased mRNA expression of glucose transporter type 4 (GLUT4), lipoprotein lipase (LpL), peroxisome proliferator-activated receptor gamma (PPARγ), leptin, osteonectin, core binding factor 1 (CBFA1), and FBJ murine osteosarcoma viral oncogene homolog (c-Fos). Expression of factors in DDT-treated cells was similar to that in estrogen-treated MSCs, suggesting that DDT may function via the estrogen receptor (ER)-mediated pathway. The coadministration of ICI 182,780 blocked the effects of DDT. RNA sequencing revealed 121 genes and noncoding RNAs to be differentially expressed in DDT-treated MSCs compared with controls cells. Conclusion: Human MSCs provide a powerful biological system to investigate and identify the molecular mechanisms underlying the effects of environmental agents on stem cells and human health. MSCs exposed to DDT demonstrated profound alterations in self-renewal, proliferation, differentiation, and gene expression, which may partially explain the homeostatic imbalance and increased cancer incidence among those exposed to long-term EDCs. Citation: Strong AL, Shi Z, Strong MJ, Miller DF, Rusch DB, Buechlein AM, Flemington EK, McLachlan JA, Nephew KP, Burow ME, Bunnell BA. 2015. Effects of the endocrine-disrupting chemical DDT on self-renewal and differentiation of human mesenchymal stem cells. Environ Health Perspect 123:42–48; http://dx.doi.org/10.1289/ehp.1408188