Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2010 Dec;84(24):13053-8.
doi: 10.1128/JVI.01521-10. Epub 2010 Oct 13.

Quantitative and qualitative RNA-Seq-based evaluation of Epstein-Barr virus transcription in type I latency Burkitt's lymphoma cells

Affiliations

Quantitative and qualitative RNA-Seq-based evaluation of Epstein-Barr virus transcription in type I latency Burkitt's lymphoma cells

Zhen Lin et al. J Virol. 2010 Dec.

Abstract

RNA-seq provides a rich source of transcriptome information with high qualitative and quantitative value. Here, we provide a pipeline for Epstein-Barr virus (EBV) transcriptome analysis using RNA-seq and we apply it to two type I latency cell lines, Mutu I and Akata. This analysis revealed substantial average expression levels of many lytic genes in predominantly latent cell populations. The lytic transcripts BHLF1 and LF3 were expressed at levels greater than those for 98% of all cellular polyadenylated transcripts. Exon junction mapping accurately identified the Qp-derived EBNA1 splicing pattern, lytic gene splicing, and a complex splicing pattern within the BamHI A region.

PubMed Disclaimer

Figures

FIG. 1.
FIG. 1.
Visualization of RNA-seq coverage across the EBV genome. Coverage (Wiggle) files generated from SAMMate and the EBV annotation file were loaded onto the Integrated Genome Viewer (IGV [http://www.broadinstitute.org/igv/], developed at the Broad Institute). The y axis shows the number of reads mapping to each location of the genome. (A) Whole-genome view; (B) zoomed view of the intergenic region between the BMRF2 and BSLF2 genes; (C) lack of reads corresponding to the EBNA2 locus. The data range for coverage data was set to 20 (for Mutu I cells) or 30 (for Akata cells), meaning that maximal peaks represent genomic positions where there were at least 20 or 30 reads that crossed that position.
FIG. 2.
FIG. 2.
RPKM values for EBV genes in Mutu I (A) and Akata (B) cells. Mutu I cell results are the averages from two technical replicates (TR) from each of two separate RNA preparations. Error bars indicate the standard deviation for each gene. (C) The number and percentage of genes showing higher RPKM values than those for BHLF1 and LF3 in Mutu I and Akata cells out of a total of 22,803 annotated cellular and viral genes.
FIG. 3.
FIG. 3.
Illustration of specificity for RNA-seq in assessing EBV transcriptomes. The total number of reads that mapped to the EBV genome per 10 million mapped reads from the EBV-positive cell lines, Akata and Mutu I, and the EBV-negative cell lines, A549 and MCF7. No EBV-specific reads in either of the EBV-negative cell lines were identified. RNA1 and RNA2 refer to biological replicate RNA samples from Mutu I cells. TR1 and TR2 refer to technical sequencing replicates.
FIG. 4.
FIG. 4.
Visualization of junction evidence for EBNA1 (A), BZLF1 (B), and BLLF1/BLLF2 (C) genes. Junction (browser extensible data [BED]) files were generated by the junction mapper TopHat as outlined in file S1 in the supplemental data.

Similar articles

Cited by

References

    1. Bernasconi, M., C. Berger, J. A. Sigrist, A. Bonanomi, J. Sobek, F. K. Niggli, and D. Nadal. 2006. Quantitative profiling of housekeeping and Epstein-Barr virus gene transcription in Burkitt lymphoma cell lines using an oligonucleotide microarray. Virol. J. 3:43. - PMC - PubMed
    1. Costa, V., C. Angelini, I. De Feis, and A. Ciccodicola. 2010. Uncovering the complexity of transcriptomes with RNA-Seq. J. Biomed Biotechnol. 2010:853916. - PMC - PubMed
    1. Dolan, A., C. Addison, D. Gatherer, A. J. Davison, and D. J. McGeoch. 2006. The genome of Epstein-Barr virus type 2 strain AG876. Virology 350:164-170. - PubMed
    1. Edwards, R. H., A. R. Marquitz, and N. Raab-Traub. 2008. Epstein-Barr virus BART microRNAs are produced from a large intron prior to splicing. J. Virol. 82:9094-9106. - PMC - PubMed
    1. Gao, Y., P. R. Smith, L. Karran, Q. L. Lu, and B. E. Griffin. 1997. Induction of an exceptionally high-level, nontranslated, Epstein-Barr virus-encoded polyadenylated transcript in the Burkitt's lymphoma line Daudi. J. Virol. 71:84-94. - PMC - PubMed

Publication types

MeSH terms

LinkOut - more resources