- Short report
- Open Access
Extensive retroviral diversity in shark
Retrovirology volume 12, Article number: 34 (2015)
Retroviruses infect a wide range of vertebrates. However, little is known about the diversity of retroviruses in basal vertebrates. Endogenous retrovirus (ERV) provides a valuable resource to study the ecology and evolution of retrovirus.
I performed a genome-scale screening for ERVs in the elephant shark (Callorhinchus milii) and identified three complete or nearly complete ERVs and many short ERV fragments. I designate these retroviral elements “C. milli ERVs” (CmiERVs). Phylogenetic analysis shows that the CmiERVs form three distinct lineages. The genome invasions by these retroviruses are estimated to take place more than 50 million years ago.
My results reveal the extensive retroviral diversity in the elephant shark. Diverse retroviruses appear to have been associated with cartilaginous fishes for millions of years. These findings have important implications in understanding the diversity and evolution of retroviruses.
Retroviruses infect a wide range of vertebrates and cause many notorious diseases, such as AIDS and cancers. However, much remains unknown about the diversity of retroviruses in basal vertebrate species. In particular, only several retroviruses have been identified in fishes, including Snakehead retrovirus, walleye dermal sarcoma virus, walleye epidermal hyperplasia virus, and Atlantic salmon swim bladder sarcoma virus [1-4]. Retrovirus employs a unique replication strategy, which requires reverse transcription of its RNA genome into DNA and integration of viral DNA into the host chromosomes. Occasionally, retroviruses infect germ line cells, and the resulting integrated retrovirus, known as endogenous retrovirus (ERV), becomes vertically inherited as a host genomic locus. Over time, some retroviral insertions are fixed in the host population. ERVs provide important insights into the ecology and evolutionary history of retroviruses.
Cartilaginous fishes (Chondrichthyes) are the most basal class of vertebrates from which retrovirus has been reported . Here, I analyzed the recently available genome sequence of the elephant shark (Callorhinchus milii), a high-quality genome assembly covering approximately 94% of the C. milii genome, for retroviral insertions . The tBLASTn algorithm with various representative retroviral Pol protein sequences was employed to screen the elephant shark genome for candidate ERV sequences. To distinguish ERVs from other LTR-retrotransposons, I used a strict criterion: only the retroviral Pol protein homolog sequence with a downstream Env protein homolog is defined as an ERV element. After initial identification of ERVs, the BLASTn algorithm was used to identify short ERV fragments. My genome-scale screening procedure identified three complete or nearly complete ERV insertions (within the C. milii genome scaffolds 2, 324, and 2324, respectively; Additional file 1: Dataset 1) and many short ERV fragments in the elephant shark genome. I designate these retroviral elements “C. milli ERVs” (CmiERVs).
To assess the relationship between CmiERVs and other retroviruses, CmiERV and representative retroviral Pol protein sequences (Additional file 2: Table S1) were aligned using MUSCLE . The ambiguous regions in the Pol protein alignments were removed using Gblocks 0.91b and then manually edited . Phylogenetic analyses were performed using MrBayes 3.1.2 . My phylogenetic analysis shows that these CmiERVs form three distinct lineages (Figure 1). Lineage I CmiERVs cluster with the retroviruses isolated from the snakehead fish (Ophicephalus striatus), while lineage II and III CmiERVs cluster with the epsilonretroviruses isolated from the walleye (Sander vitreus) and amphibians. CmiERV lineage I and lineages II/III are only distantly related to each other. These CmiERV lineages are likely to result from three independent retroviral invasion events. My analysis provides clear evidence there is extensive retroviral diversity in the elephant shark. It is possible that exogenous retroviruses related to CmiERVs identified here are still circulating in the elephant shark and possibly other Chondichthyes.
On endogenization, the 5′LTR and 3′LTR of a nascent ERV are identical and will accumulate mutations independently. Thus, the 5′LTR and 3′LTR genetic divergence could be used to estimate ERV invasion time . The invasion time of an ERV can be estimated through:
where t indicates the invasion time, u indicates the neutral evolutionary rate of host, and d indicates the genetic divergence between 5′ LTR and 3′ LTR. In this study, two complete CmiERV insertions were identified (Table 1). The genetic divergence between 5′-LTR and 3′-LTR was calculated with the Kimura two-parameter substitution model . The neutral evolutionary rate for the elephant shark is not available but is approximately an order of magnitude lower than those for mammals . The average neutral rate for mammals is estimated to be 2.2 x 10−9 substitutions per site per year . I thus assume 2.2 x 10−10 substitutions per site per year as the shark neutral rate. Then the insertion times for the two complete ERVs are approximately 75 and 54.5 million years ago, respectively. However, these estimates should be taken with cautions , given that I am not sure whether LTRs evolve at a neutral manner and what the actual neutral rate for the elephant shark is. Nevertheless, these results suggest these retroviruses were infecting the elephant shark millions of years ago.
Previously, a single ERV sequence was identified in the lemon shark (Negaprion brevirostris), which is closely related to human ERV; this ERV was thought to have a cross-transmission origin . However, I find that CmiERVs cluster together with retroviruses of fish origin. The phylogenetic pattern is compatible with the hypothesis of an ancient marine origin of retroviruses . Chondrichthyes are the most basal class of vertebrates from which retrovirus has been identified; no retrovirus is identified in earlier-diverging vertebrate lineages, the lampreys (Cephalaspidomorphi) and the hagfish (Myxini) . It follows that these CmiERV elements are likely to represent “primitive” retroviruses. However, the possibility that these elephant shark retroviruses originated from cross-transmission from other fishes cannot be formally excluded.
To date, only a limited number of exogenous/endogenous retroviruses have been identified in fishes [1-5,14]. My results reveal the unexpectedly extensive retroviral diversity of the elephant shark. The initial candidate ERVs were identified based on a strict criterion – whether there is a downstream Env protein homolog following the Pol protein homolog. This approach is conservative, given that the retroviral Env protein evolves rapidly and its similarity to other retroviral Env proteins will erode over a long time. On the other hand, the ERVs identified using this approach are authentic retroviruses. It is likely that there are additional ERV insertions that were not detected. Also, it should be noted that only a small proportion of retroviruses could leave endogenous copies in their host genomes . Therefore, I believe the actual diversity of retroviruses is more extensive in the elephant shark. Further analysis of ERV in basal vertebrates would improve our understanding of the diversity and evolution of retroviruses.
Hart D, Frerichs GN, Rambaut A, Onions DE. Complete nucleotide sequence and transcriptional analysis of snakehead fish retrovirus. J Virol. 1996;70:3606–16.
LaPierre LA, Holzschu DL, Bowser PR, Casey JW. Sequence and transcriptional analyses of the fish retroviruses walleye epidermal hyperplasia virus types 1 and 2: evidence for a gene duplication. J Virol. 1999;73:9393–403.
Holzschu DL, Martineau D, Fodor SK, Vogt VM, Bowser PR, Casey JW. Nucleotide sequence and protein analysis of a complex piscine retrovirus, walleye dermal sarcoma virus. J Virol. 1995;69:5320–31.
Paul TA, Quackenbush SL, Sutton C, Casey RN, Bowser PR, Casey JW. Identification and characterization of an exogenous retrovirus from atlantic salmon swim bladder sarcomas. J Virol. 2006;80:2941–8.
Herniou E, Martin J, Miller K, Cook J, Wilkinson M, Tristem M. Retroviral diversity and distribution in vertebrates. J Virol. 1998;72:5955–66.
Venkatesh B, Lee AP, Ravi V, Maurya AK, Lian MM, Swann JB, et al. Elephant shark genome provides unique insights into gnathostome evolution. Nature. 2014;505:174–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56:564–77.
Ronquist F, Huelsenbeck JP. Mrbayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Johnson WE, Coffin JM. Constructing primate phylogenies from ancient retrovirus sequences. Proc Natl Acad Sci USA. 1999;96:10254–60.
Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.
Martin AP. Substitution rates of organelle and nuclear genes in sharks: implicating metabolic rate (again). Mol Biol Evol. 1999;16:996–1002.
Kumar S, Subramanian S. Mutation rates in mammalian genomes. Proc Proc Natl Acad Sci USA. 2002;99:803–8.
Han GZ, Worobey M. A primitive endogenous lentivirus in a colugo: insights into the early evolution of lentiviruses. Mol Biol Evol. 2015;32:211–5.
Han GZ, Worobey M. An endogenous foamy-like viral element in the coelacanth genome. PLoS Pathog. 2012;8:e1002790.
Emerman M, Malik HS. Paleovirology–modern consequences of ancient viruses. PLoS Biol. 2010;8:e1000301.
This research was supported by the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions. I thank Chunyang Wang for helping prepare Figure 1.
The authors declare that they have no competing interests.
GZH designed and performed the research, analyzed the data, and wrote the manuscript.
About this article
Cite this article
Han, G. Extensive retroviral diversity in shark. Retrovirology 12, 34 (2015) doi:10.1186/s12977-015-0158-4
- Endogenous retroviruses