Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

Host genome integration and giant virus-induced reactivation of the virophage mavirus

Abstract

Endogenous viral elements are increasingly found in eukaryotic genomes1, yet little is known about their origins, dynamics, or function. Here we provide a compelling example of a DNA virus that readily integrates into a eukaryotic genome where it acts as an inducible antiviral defence system. We found that the virophage mavirus2, a parasite of the giant Cafeteria roenbergensis virus (CroV)3, integrates at multiple sites within the nuclear genome of the marine protozoan Cafeteria roenbergensis4. The endogenous mavirus is structurally and genetically similar to eukaryotic DNA transposons and endogenous viruses of the Maverick/Polinton family5,6,7. Provirophage genes are not constitutively expressed, but are specifically activated by superinfection with CroV, which induces the production of infectious mavirus particles. Virophages can inhibit the replication of mimivirus-like giant viruses and an anti-viral protective effect of provirophages on their hosts has been hypothesized2,8. We find that provirophage-carrying cells are not directly protected from CroV; however, lysis of these cells releases infectious mavirus particles that are then able to suppress CroV replication and enhance host survival during subsequent rounds of infection. The microbial host–parasite interaction described here involves an altruistic aspect and suggests that giant-virus-induced activation of provirophages might be ecologically relevant in natural protist populations.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Generation and characterization of endogenous mavirus virophages in Cafeteria roenbergensis.
Figure 2: Gene expression analysis of the endogenous mavirus genome.
Figure 3: CroV infection induces replication and virion production of the endogenous mavirus.
Figure 4: Reactivated mavirus inhibits CroV and promotes host survival in subsequent infections.

Similar content being viewed by others

References

  1. Aiewsakun, P. & Katzourakis, A. Endogenous viruses: connecting recent and ancient viral evolution. Virology 479–480, 26–37 (2015)

    Article  Google Scholar 

  2. Fischer, M. G. & Suttle, C. A. A virophage at the origin of large DNA transposons. Science 332, 231–234 (2011)

    Article  ADS  CAS  Google Scholar 

  3. Fischer, M. G., Allen, M. J., Wilson, W. H. & Suttle, C. A. Giant virus with a remarkable complement of genes infects marine zooplankton. Proc. Natl Acad. Sci. USA 107, 19508–19513 (2010)

    Article  ADS  CAS  Google Scholar 

  4. Fenchel, T. & Patterson, D. J. Cafeteria roenbergensis nov. gen., nov. sp., a heterotrophic microflagellate from marine plankton. Mar. Microb. Food Webs 3, 9–19 (1988)

    Google Scholar 

  5. Pritham, E. J., Putliwala, T. & Feschotte, C. Mavericks, a novel class of giant transposable elements widespread in eukaryotes and related to DNA viruses. Gene 390, 3–17 (2007)

    Article  CAS  Google Scholar 

  6. Kapitonov, V. V. & Jurka, J. Self-synthesizing DNA transposons in eukaryotes. Proc. Natl Acad. Sci. USA 103, 4540–4545 (2006)

    Article  ADS  CAS  Google Scholar 

  7. Krupovic, M. & Koonin, E. V. Polintons: a hotbed of eukaryotic virus, transposon and plasmid evolution. Nature Rev. Microbiol. 13, 105–115 (2015)

    Article  CAS  Google Scholar 

  8. Blanc, G., Gallot-Lavallée, L. & Maumus, F. Provirophages in the Bigelowiella genome bear testimony to past encounters with giant viruses. Proc. Natl Acad. Sci. USA 112, E5318–E5326 (2015)

    Article  ADS  CAS  Google Scholar 

  9. Krupovic, M., Kuhn, J. H. & Fischer, M. G. A classification system for virophages and satellite viruses. Arch. Virol. 161, 233–247 (2016)

    Article  CAS  Google Scholar 

  10. La Scola, B. et al. The virophage as a unique parasite of the giant mimivirus. Nature 455, 100–104 (2008)

    Article  ADS  CAS  Google Scholar 

  11. Krupovic, M. & Koonin, E. V. Self-synthesizing transposons: unexpected key players in the evolution of viruses and defense systems. Curr. Opin. Microbiol. 31, 25–33 (2016)

    Article  CAS  Google Scholar 

  12. Yutin, N., Shevchenko, S., Kapitonov, V., Krupovic, M. & Koonin, E. V. A novel group of diverse Polinton-like viruses discovered by metagenome analysis. BMC Biol. 13, 95 (2015)

    Article  Google Scholar 

  13. Desnues, C. et al. Provirophages and transpovirons as the diverse mobilome of giant viruses. Proc. Natl Acad. Sci. USA 109, 18078–18083 (2012)

    Article  ADS  CAS  Google Scholar 

  14. Santini, S. et al. Genome of Phaeocystis globosa virus PgV-16T highlights the common ancestry of the largest known DNA viruses infecting eukaryotes. Proc. Natl Acad. Sci. USA 110, 10800–10805 (2013)

    Article  ADS  CAS  Google Scholar 

  15. Fischer, M. G., Kelly, I., Foster, L. J. & Suttle, C. A. The virion of Cafeteria roenbergensis virus (CroV) contains a complex suite of proteins for transcription and DNA repair. Virology 466-467, 82–94 (2014)

    Article  CAS  Google Scholar 

  16. Resch, W., Hixson, K. K., Moore, R. J., Lipton, M. S. & Moss, B. Protein composition of the vaccinia virus mature virion. Virology 358, 233–247 (2007)

    Article  CAS  Google Scholar 

  17. Renesto, P. et al. Mimivirus giant particles incorporate a large fraction of anonymous and unique gene products. J. Virol. 80, 11678–11685 (2006)

    Article  CAS  Google Scholar 

  18. Legendre, M. et al. In-depth study of Mollivirus sibericum, a new 30,000-y-old giant virus infecting Acanthamoeba. Proc. Natl Acad. Sci. USA 112, E5327–E5335 (2015)

    Article  CAS  Google Scholar 

  19. González, J. M. & Suttle, C. A. Grazing by marine nanoflagellates on viruses and virus-sized particles - ingestion and digestion. Mar. Ecol. Prog. Ser. 94, 1–10 (1993)

    Article  ADS  Google Scholar 

  20. Garza, D. R. & Suttle, C. A. Large double-stranded DNA viruses which cause the lysis of a marine heterotrophic nanoflagellate (Bodo sp) occur in natural marine viral communities. Aquat. Microb. Ecol. 9, 203–210 (1995)

    Article  Google Scholar 

  21. Reed, L. & Muench, H. A simple method of estimating fifty per cent end points. Am. J. Hyg. 27, 493–497 (1938)

    Google Scholar 

  22. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014)

    Article  CAS  Google Scholar 

  23. Hackl, T., Hedrich, R., Schultz, J. & Förster, F. proovread: large-scale high-accuracy PacBio correction through iterative short read consensus. Bioinformatics 30, 3004–3011 (2014)

    Article  CAS  Google Scholar 

  24. Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012)

    Article  MathSciNet  CAS  Google Scholar 

  25. Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013)

    Article  CAS  Google Scholar 

  26. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410 (1990)

    Article  CAS  Google Scholar 

  27. Wick, R. R., Schultz, M. B., Zobel, J. & Holt, K. E. Bandage: interactive visualization of de novo genome assemblies. Bioinformatics 31, 3350–3352 (2015)

    Article  CAS  Google Scholar 

  28. Magoč, T. & Salzberg, S. L. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963 (2011)

    Article  Google Scholar 

  29. Li, H. Toward better understanding of artifacts in variant calling from high-coverage samples. Bioinformatics 30, 2843–2851 (2014)

    Article  CAS  Google Scholar 

  30. Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009)

    Article  Google Scholar 

  31. Skinner, M. E., Uzilov, A. V., Stein, L. D., Mungall, C. J. & Holmes, I. H. JBrowse: a next-generation genome browser. Genome Res. 19, 1630–1638 (2009)

    Article  CAS  Google Scholar 

  32. Campbell, M. S., Holt, C., Moore, B. & Yandell, M. in Current Protocols in Bioinformatics vol. 48, 4.11.1–4.11.39 (John Wiley, 2014)

    Article  Google Scholar 

  33. Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014)

    Article  CAS  Google Scholar 

  34. Kent, W. J. et al. The human genome browser at UCSC. Genome Res. 12, 996–1006 (2002)

    Article  CAS  Google Scholar 

  35. Marçais, G. & Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27, 764–770 (2011)

    Article  Google Scholar 

Download references

Acknowledgements

This research was supported by the Max Planck Society. We are grateful to C. Suttle for access to host and virus strains, and to the Roscoff team for maintaining and distributing protist strains. We thank K.-A. Seifert, K. Fenzl and K. Barenhoff for technical assistance, U. Mersdorf for electron microscopy, C. Roome for IT support, L. Czaja and the Max Planck Genome Centre in Cologne for bioinformatic assistance, S. Higgins for suggestions, K. Haslinger and J. Reinstein for comments on the manuscript, and I. Schlichting for mentoring and support.

Author information

Authors and Affiliations

Authors

Contributions

M.G.F. conceived the study, designed and performed experiments, collected, analysed, and interpreted data, and wrote the manuscript. T.H. corrected and assembled sequence data, and analysed, interpreted, and visualized data. Both authors revised and approved the manuscript.

Corresponding author

Correspondence to Matthias G. Fischer.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Additional information

Reviewer Information

Nature thanks J.-M. Claverie, E. Koonin and the other anonymous reviewer(s) for their contribution to the peer review of this work.

Extended data figures and tables

Extended Data Figure 1 CroV-mavirus co-infection experiment of C. roenbergensis strain E4-10P to generate the mavirus-carrying strain E4-10M1.

C. roenbergensis strain E4-10P was either mock-infected, infected with CroV, or co-infected with CroV and mavirus. Cell densities are based on microscopy counts and were monitored for 8 days, viral numbers were monitored for 6 days and are derived from qPCR data assaying short amplicons of the mavirus MV18 gene (MCP) and the crov283 gene (D11-like transcription factor), respectively. The detection limit for both methods was about 103 per millilitre. These experiments were performed in single copies. The mavirus-positive host strain E4-10M1 was isolated from the pool of surviving cells in the +mavirus, +CroV infection. See also Supplementary Spreadsheet 4.

Source data

Extended Data Figure 2 The k-mer frequency distribution of the diploid genome of Cafeteria roenbergensis strain E4-10P.

Frequency distribution of in silico generated random 19-mers in the genomic read set of E4-10P. The distribution exhibits a major peak at 120× coverage (solid red line) corresponding to the majority of homozygous k-mers of the underlying diploid genome, a smaller peak at half the diploid coverage (60×, dotted green line) comprising haplotype-specific k-mers, and a weak third peak at three times the haploid coverage (180×, dashed blue line) indicating a primarly diploid, partly triploid genome structure. Low-coverage k-mers derive from sequencing errors and bacterial contamination.

Extended Data Figure 3 Gene expression of the endogenous mavirus genome is inhibited by cycloheximide (CHX) and aphidicolin (APH).

Selected cellular and viral transcripts isolated at 0 and 24 h p.i. from mock-infected or CroV-infected E4-10P and E4-10M1 cultures in the presence of 5 μg/ml aphidicolin or 50 μg/ml CHX were quantified by qRT-PCR. Shown are the average quantification cycle (Cq) values of three independent experiments; error bars, s.d. The following genes were assayed: host AspRS, C. roenbergensis E4-10 aspartyl-tRNA synthetase; crov342, CroV major capsid protein; crov497, CroV DNA polymerase B; crov505, CroV isoleucyl-tRNA synthetase; MV03, mavirus DNA polymerase B; MV15, mavirus genome-packaging ATPase; MV16, mavirus maturation protease; MV17, mavirus minor capsid protein; MV18, mavirus major capsid protein. Cq values of the reverse transcriptase-negative (−RT) reactions are shown directly to the right of the respective +RT results. Accession numbers are listed in Extended Data Table 3. See also Supplementary Spreadsheet 1.

Source data

Extended Data Figure 4 Hypothesis for CroV-induced reactivation of endogenous mavirus.

Shown is a schematic C. roenbergensis cell displaying selected events of a CroV infection cycle in strains E4-10P (left) and E4-10M1 (right). Following CroV entry (1), the virion factory forms in the cytoplasm. At the onset of late phase, a CroV-encoded transcription factor recognizing the late CroV promoter motif is synthesized (2). We hypothesize that in E4-10M1 cells, the late transcription factor could enter the nucleus (3), bind the mavirus promoter sequences, and activate gene expression of the provirophage (4). Mavirus-specific transcripts would then be exported and translated (5) and some of the mavirus proteins would return to the nucleus to excise or replicate the provirophage genome (6). The mavirus genome could then translocate to the CroV factory (7), where genome replication, particle assembly, and genome packaging would occur (8). Cell lysis releases the newly synthesized CroV and mavirus particles (9) and the reactivated virophages inhibit CroV propagation during subsequent co-infections (10). By contrast, CroV infection of an E4-10P cell does not induce a virophage response and CroV continues to infect other host populations (11).

Extended Data Figure 5 Purification and characterization of reactivated mavirus particles.

a, Three-litre cultures of CroV- or mock-infected E4-10P and E4-10M1 cultures were concentrated 200-fold and stained with uranyl acetate for electron microscopy. Representative particles from each filtrate are boxed and shown at higher magnification. Mavirus-like particles are marked by arrows. Filtrate numbers refer to the infections shown in Fig. 3a. b, The concentrated samples were analysed on 1.1–1.5 g/ml linear CsCl density gradients. A concentrated sample of reference mavirus was run in parallel and yielded a band at approximately 1.29 g/ml CsCl (arrow). Only the CroV-infected E4-10M1 culture produced a band at a similar density. c, PCR analysis of band material extracted from the CsCl gradients shown in b at a density of about 1.29 g/ml CsCl. Primers MaV21F and MaV21R were used to generate a 956-bp-long product of the MV19 gene. d, Material from the 1.29 g/ml CsCl band or from equivalent positions was extracted from the gradients and visualized by negative-stain electron microscopy. Only the CroV-infected E4-10M1 culture contained mavirus-like particles.

Extended Data Figure 6 Ultraviolet light treatment abolishes infectivity of reactivated mavirus.

Reactivated mavirus contained in filtrate 4 from the infection experiments shown in Fig. 3a was irradiated with 500 J m−2 of ultraviolet-C light (λ = 254 nm). Ultraviolet light-treated and untreated mavirus suspensions were tested for infectivity by co-infection of host strain E4-10P with CroV. As shown in the lower right panel, reactivated mavirus treated with ultraviolet light was no longer able to replicate in the presence of CroV. Data are shown as the mean of biological triplicates ± s.d. See also Supplementary Spreadsheet 5.

Source data

Extended Data Table 1 Cell and virus concentrations in different size fractions of mock-infected and CroV-infected E4-10P and E4-10M1 populations
Extended Data Table 2 Details on the 11 bioinformatically well-supported mavirus integration sites in C. roenbergensis strain E4-10M1
Extended Data Table 3 PCR oligonucleotide primers used in this study

Supplementary information

Supplementary Spreadsheet 1

This file contains source data for figure 2 and Extended Data Figure 3. (XLSX 53 kb)

Supplementary Spreadsheet 2

This file contains source data for figure 3. (XLSX 54 kb)

Supplementary Spreadsheet 3

This file contains source data for figure 4. (XLSX 58 kb)

Supplementary Spreadsheet 4

This file contains source data for Extended Data Figure 1. (XLSX 12 kb)

Supplementary Spreadsheet 5

This file contains source data for Extended Data Figure 6. (XLSX 23 kb)

PowerPoint slides

Source data

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fischer, M., Hackl, T. Host genome integration and giant virus-induced reactivation of the virophage mavirus. Nature 540, 288–291 (2016). https://doi.org/10.1038/nature20593

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nature20593

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing