Widespread Horizontal Gene Transfer from Circular ... - BioMedSearch

7 downloads 0 Views 3MB Size Report
Sep 26, 2011 - Some of the transferred viral genes were conserved and expressed, suggesting that these ... Viruses with circular single stranded DNA (ssDNA).
Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

RESEARCH ARTICLE

Open Access

Widespread Horizontal Gene Transfer from Circular Single-stranded DNA Viruses to Eukaryotic Genomes Huiquan Liu1,2,3, Yanping Fu2, Bo Li1,2, Xiao Yu1,2, Jiatao Xie2, Jiasen Cheng2, Said A Ghabrial4, Guoqing Li1,2, Xianhong Yi2 and Daohong Jiang1,2*

Abstract Background: In addition to vertical transmission, organisms can also acquire genes from other distantly related species or from their extra-chromosomal elements (plasmids and viruses) via horizontal gene transfer (HGT). It has been suggested that phages represent substantial forces in prokaryotic evolution. In eukaryotes, retroviruses, which can integrate into host genome as an obligate step in their replication strategy, comprise approximately 8% of the human genome. Unlike retroviruses, few members of other virus families are known to transfer genes to host genomes. Results: Here we performed a systematic search for sequences related to circular single-stranded DNA (ssDNA) viruses in publicly available eukaryotic genome databases followed by comprehensive phylogenetic analysis. We conclude that the replication initiation protein (Rep)-related sequences of geminiviruses, nanoviruses and circoviruses have been frequently transferred to a broad range of eukaryotic species, including plants, fungi, animals and protists. Some of the transferred viral genes were conserved and expressed, suggesting that these genes have been coopted to assume cellular functions in the host genomes. We also identified geminivirus-like and parvovirus-like transposable elements in genomes of fungi and lower animals, respectively, and thereby provide direct evidence that eukaryotic transposons could derive from ssDNA viruses. Conclusions: Our discovery extends the host range of circular ssDNA viruses and sheds light on the origin and evolution of these viruses. It also suggests that ssDNA viruses act as an unforeseen source of genetic innovation in their hosts.

Background In addition to vertical transmission and gene acquisition from other distantly related species via horizontal gene transfer (HGT), organisms can also capture genetic material from extra-chromosomal elements (plasmids and viruses) during evolution. It is widely accepted that phages represent substantial forces in prokaryotic evolution, with the integrated phages (prophages) accounting for as much as 10-20% of some bacterial genomes [1,2]. In eukaryotes, animal retroviruses, which can integrate into host genome as an obligate step in their replication strategy, comprise approximately 8% of the human * Correspondence: [email protected] 1 State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan 430070, Hubei Province, P R China Full list of author information is available at the end of the article

genome in the form of inherited endogenous retroviruses [3]. Moreover, the integrated retroviral genes have been demonstrated to play critical role in mammalian reproduction [4,5]. Recent data reveal that several non-retroviral viruses have also contributed to the genetic makeup of many eukaryotic organisms [6-15]. Especially, genes derived from ancestral nudiviruses have been co-opted to facilitate a parasitic lifestyle in parasitoid wasps [7]; and a gene derived from partitiviruses was exapted to regulate the activities of the phytohormone auxin, indole-3-acetic acid (IAA) in Arabidopsis thaliana [12]. Still, this type of transfer is thought to be rare in eukaryotes. Viruses with circular single stranded DNA (ssDNA) genomes are the smallest viruses known to infect eukaryotes and are currently grouped into four families:

© 2011 Liu et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Anelloviridae, Circoviridae, Geminiviridae and Nanoviridae (Virus Taxonomy: 2009, ICTV, http://www.ictvonline.org/virusTaxonomy.asp?version = 2009). The members of the first two families infect vertebrates and of the last two families infect plants. Recently a virus distantly related to circoviruses carrying a covalently closed circular, partially double-stranded ssDNA genome has been found to infect the marine diatom Chaetoceros salsugineum [16]. A similar virus was also discovered in C. debilis [17]. Moreover, recent viral metagenomic studies have shown that small circular ssDNA viruses are more prevalent and diverse in the environment than previously recognized [18-22]. Small circular ssDNA viruses commonly replicate their genomes in the nuclei of infected cells via a rolling circle replication (RCR) mechanism initiated by virus-encoded replication initiation protein (Rep), and there are clear similarities among the sequences of these proteins [23,24]. So far, no associated integrase activity has been identified for these viruses. However, Bejarano et al [25] reported multiple repeats of geminiviral Rep DNA that have been integrated into the nuclear genome of tobacco. In addition, Rep-like genes were also found in genomes of the parasitic protozoan Entamoeba histolytica and Giardia intestinalis [26]. These discoveries suggest that the small circular ssDNA viruses could also contribute to the genetic heritage of eukaryotic organisms. Considering that the circular ssDNA viruses are widespread in nature, the role played by these viruses in eukaryotic evolution needs to be evaluated. Accordingly, we performed a systematic search for sequences related to known small circular ssDNA viruses in the publicly available eukaryotic genome databases. As our study was being prepared for publication, an independent group of investigators reported that sequences related to circoviruses were detected in the genomes of six vertebrate species [27]. Here we report our more comprehensive and convincing results based on sufficiently critical data analysis, bench research and phylogenetic analysis. Our studies have not only corroborated the integration of circovirus-related sequences in these six species, but they have also revealed that numerous sequences related to circoviruses, geminiviruses and nanoviruses have been integrated into the germlines of diverse eukaryotes including plants, fungi, animals and protists. Furthermore, we have demonstrated some of these integrated genes were conserved and expressed in eukaryotic organisms. In addition, we also identified geminivirus-like and parvovirus-like transposable elements in the genomes of fungi and lower animals, respectively. The origin and evolution of small circular ssDNA viruses were also discussed.

Page 2 of 15

Results and discussion Identification of circular ssDNA viral Rep-related proteins in eukaryotic systems

Rep proteins are commonly encoded by mobile elements (most phages and eukaryotic ssDNA viruses, some plasmids of Gram-positive bacteria, eukaryotic helitron transposons, etc.) but without cellular homologs and therefore have been recognized as virus/plasmid-specific proteins (hallmark proteins) [28,29]. The Rep proteins of eukaryotic ssDNA viruses contain RCR catalytic domain and a C-terminal NTPase/helicase domain [30,31]. With such structure, the sequence of the Rep protein of geminiviruses readily detected those of the geminivirus Rep catalytic domain (Gemini_AL1), central domain (Gemini_AL1_M) and the RNA helicase domain (RNA_helicase) by using the NCBI Conserved Domain Database searches (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) (see Additional file 1: Figure S1). Likewise, the Rep of nanoviruses and circoviruses detected the putative viral replication protein domain (Viral_Rep) and the RNA helicase domain. On the other hand, the Rep of pLS1 family of prokaryotic plasmids [32] comprises only plasmid replication protein domain (Rep_2). Recently, the plasmid Rep containing an additional helicase domain has also been reported in a phytoplasma [33]. To investigate sequences closely related to Rep proteins of eukaryotic circular ssDNA viruses in other systems, we used the representative Rep proteins of geminiviruses and circoviruses to search against NCBI non-redundant (NR) protein database by PSI-BLAST [34]. After removing the known eukaryotic circular ssDNA viruses and subsequent reverse BLAST comparisons, we obtained a final dataset containing 113 Replike protein sequences from plasmids, other eukaryotic viruses and cellular genomes respectively (see Additional file 2: Tabular data S1). Most of these belonged to bacterial plasmids (typically phytoplasmal plasmids) and bacterial genomes. Conserved Domain searches showed that, apart from the known Rep_2 plus RNA_helicase domains, Rep proteins from bacterial plasmids had two other types of domain arrangements: Gemini_AL1 plus RNA_helicase and Viral_Rep plus RNA_helicase (see Additional file 1: Figure S1 and Additional file 2: Tabular data S1). The domain organization of Rep-like proteins from bacterial genomes showed high similarity to those of plasmids suggesting that these cellular homologs may originate from their plasmid counterparts. One apparently truncated Rep-like protein from the mitochondrion of oomycete Phytophthora sojae [35] contained Gemini_AL1 domain. Interestingly, we also found its coordinates in mitochondrion type II haplotypes of P. infestans [36] (Figure 1A), which has not

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 3 of 15

Phytophthora infestans haplotype IIa mitochondrion AY898627.1 4490 44 49

ORF99

ORF39

ORF45

ORF80

6800 80000 bpp

ORF198

ORF170

ORF51

A

~70% nt ORF98 30 30500 05

ORF101

28500 000 bpp

ORF112a ORF112b

ORF116

Phytophthora sojae mitochondrion NC_009385.1

Giardia intestinalis ATCC 50581 strain GS/M H7 ACGJ01002178.1 1

2155 bp

EET01111 35% aa

Human gut viral metagenome FI579842.1

1

1128 bp

ψ

G. intestinalis isolate BRIS/92/HEPU/1541 AF059664.1

60% aa

B

22000 0

ψ

AAF28770

5583 bp

AAF28772

AAF28771 65% aa

Human gut viral metagenome FI587624.1

1

ψ

564 5

ψ

2678 bp 267

ψ 98% nt

2128

Figure 1 Integrated plasmid or virus-like genes in Phytophthora sp. (A) and Giardia intestinalis (B). Arrowhead boxes indicate ORFs (orange, Rep-like genes; other colors, unknown genes). Gray sectors connect corresponding homologous regions and the % nucleotide (nt) or amino acid (aa) identity are indicated. The annotated ORF names or accession numbers are indicated. ψ, interrupted ORF.

been identified in previous reports. The mitochondrial regions containing the Rep-like sequences were absent in mitochondrion of P. ramorum and the type I haplotypes of P. infestans [36]. In addition, they are most closely related to the Reps of plasmids from the red algae Porphyra tenera [37]. These findings suggest that the Rep-like sequences are most likely derived from an integrated plasmid. Two circovirus Rep-related sequences have been found in the genome of canarypox virus [38] but not in other poxviruses, thus suggesting that they were acquired horizontally. Our dataset also included the previously reported Replike genes in the genome of G. intestinalis isolate BRIS/ 92/HEPU/1541 and E. histolytica HM-1:IMSS [26]. These genes have been shown to be present in G. intestinalis ATCC 50581 strain GS/M H7 but not present in the G. intestinalis ATCC 50803 WB genome [39]. We

found that two sequences of human gut viral metagenome from Genomic Survey Sequence (GSS) database share sequence similarities not only with the Rep-like genes of G. intestinalis but also with their neighboring genes (Figure 1B). This could provide evidence that these genes were originated from integrated proviruses. In addition to the known Rep-like genes in genomes of E. histolytica and G. intestinalis, we also identified 20 new Rep-like genes in genomes of other protozoan species as well as fungi, placozoans, and roundworms (see Additional file 2: Tabular data S1). Among these, the fungal Rep-like proteins contained the geminivirus-like domain, while the rest have domain similarities with Reps of nanoviruses or circoviruses. These findings suggest that the eukaryotic circular ssDNA viral genes may be of widespread occurrence in their host genomes but have yet to be discovered.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 4 of 15

Identification and validation of eukaryotic endogenous circular ssDNA virus-like sequences in germline genomes

Characteristics and phylogenies of endogenous circular ssDNA virus-like sequences

To conduct this survey, we performed a comprehensive BLAST searches using as queries the viral Rep-like proteins in eukaryotic genomes and the protein sequences of representative eukaryotic circular ssDNA viruses against the genomic assemblies of 209 eukaryotes plus other uncompleted eukaryotic genomes in HTGS, WGS and GSS databases. This process identified 305 significant matches to Rep or CP proteins of geminiviruses, nanoviruses and circoviruses (Table 1 and Additional file 2: Tabular data S2). The most abundant of these virus-like sequences were related to the Rep genes, while only three sequences (one geminivirus-like, two circovirus-like) were related to CP genes, which were detected in tobacco and sloth genomes respectively. This finding is consistent with previous metagenomic research; possibly due to the fact that the Rep gene is more conserved during evolution. However, the possibility that possession of Rep sequences may offer some selective advantage to the host species cannot be ruled out either. Because circular ssDNA viruses replicate their genomes in the nuclei of their host cells, sequence similarities between these viruses and host genomes could be attributed to trivial contamination of eukaryotic genomic DNA with viral sequences during cloning or sequence assembly. To rule out this possibility, we did a closer inspection of the raw sequence reads used for WGS assembly and the results indicated deep sequencing coverage across the junctions between the endogenous virus-like sequences and adjacent cellular sequences (see Additional file 2: Tabular data S2). These results suggest that the endogenous viral sequences were not artifacts of cloning or sequence assembly. For some eukaryotic genomes, the trace archives were not available, but their endogenous viral sequences underwent various degrees of degradation (see Additional file 2: Tabular data S2), suggesting that the viral sequences had invaded host genomes millions of years ago and therefore represent established germline infections. To validate these observations, we amplified and sequenced the endogenous circovirus-like sequences and their flanking host sequences from dog and cat tissues (Figure 2). The results revealed that the PCR products were of the expected sizes and the experimental sequences were identical to relevant regions of sequenced animal genomes. Altogether, we discovered endogenous virus-like sequences in at least 35 species broadly distributed among nuclear genomes of plants, fungi, animals and protists. Remarkably, no anellovirus-like sequence was detected in any eukaryotic genome, although these viruses have been noted in various animal species [40].

Compared to related exogenous viral genes, some endogenous virus-like sequences are full-length or near fulllength genes while many others comprise only gene fragments. Despite pronounced sequence divergence, the conserved motifs of Rep protein can still be easily found in putative protein sequences of endogenous virus-like genes (see Additional file 1: Figure S2 and S3). The endogenous viral sequences are generally interspersed within non-coding regions of host genomes but several were found to be inserted into the coding regions of host genes or transposons (Figure 3). This finding suggests that these viral sequences have influenced host genome evolution through gene disruption. We next performed a comprehensive phylogenetic analysis to determine the relationship between the endogenous virus-like sequences and known circular ssDNA viruses as well as among endogenous sequences (Figure 4 and Additional file 1: Figure S4). As shown in Additional file 1: Figure S4, all the sequences formed three large clades: geminivirus-like, nanovirus-like and circovirus-like. In each clade, the endogenous virus-like sequences generally clustered distinctly with the known viruses but did not fall into the established viral families, suggesting that these virus-like sequences may have originated from previously undescribed circular ssDNA viral lineages. An exception was the only virus-like sequence detected in opossum (Monodelphis domestica) genome, which clustered within the Circoviridae clade and was most closely related to pig circoviruses. Unlike Rep-like sequence in tobacco that were acquired more recently from members of begomovirus [41], one genera in the family Geminiviridae, the geminiviral Rep-like sequence in Populus, was located at the base of the Geminiviridae clade in the phylogenetic tree (Figure 4), suggesting that it was derived from integration of a Geminiviridae ancestor. Indeed, this sequence was degenerate, containing three inframe stop codons and one frameshift, an indication that it has been inserted a million years ago. Alternatively, it represents a distantly related geminiviral lineage infecting Populus. All the virus-like sequences from fungi clustered together and were most closely related to the Sclerotinia sclerotiorum hypovirulence associated DNA virus 1 (SsHADV-1) (Figure 4), a mycovirus recently reported by us [42], suggesting that these endogenous viral sequences originated from SsHADV-1 like mycoviruses. Moreover, the SsHADV-1-like Reps were prevalent in viral metagenomes of different samples, including freshwater, human gut, rice paddy soil, marine environments and mosquito (see Additional file 1: Figure S5). Our phylogenetic analysis also suggests that the circular ssDNA viruses were likely to co-evolve with their

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 5 of 15

Table 1 Numbers of endogenous circular ssDNA virus-like sequences in eukaryotic genomes Organism group

Organism

No. of virus-related genes Rep

Capsid

Plants land plants

Populus trichocarpa (black cottonwood)

1

Nicotiana tabacum (common tobacco) green algae

1

Micromonas pusilla (green algae) CCMP1545

1

Aspergillus nidulans FGSC A4

1

Aspergillus fumigatus A1163 Aspergillus niger CBS 513.88

1 1

Fungi ascomycetes

basidiomycetes

Trichoderma atroviride IMI 206040

1

Magnaporthe oryzae 70-15 (rice blast fungus)

1

Nectria haematococca mpVI 77-13-4

4

Tuber melanosporum Mel28 (Perigord truffle)

42

Laccaria bicolor S238N-H82 (Bicoloured deceiver)

5

Entamoeba invadens IP1 Entamoeba terrapinae

10 3

Entamoeba histolytica HM-1:IMSS

14

Entamoeba dispar SAW760

7

Protists protozoans

diatoms Animals mammals

Blastocystis hominis Singapore isolate B (sub-type 7)

7

Giardia intestinalis ATCC 50581 strain GS/M H7

13

Giardia intestinalis isolate BRIS/92/HEPU/1541

2

Phaeodactylum tricornutum (diatom)

1

Canis lupus familiaris (dog) *

4

Monodelphis domestica (gray short-tailed opossum) *

1

Felis catus(domestic cat) *

6

Ailuropoda melanoleuca (giant panda) *

12

Choloepus hoffmanni (Hoffmann’s two-fingered sloth) *

2†

gastropods

Aplysia californica (California sea hare)

1

amphibians

Xenopus (Silurana) tropicalis(western clawed frog) *

2

lancelets roundworms

Branchiostoma floridae (Florida lancelet) strain S238N-H82 Brugia malayi (agent of lymphatic filariasis)

7 1

Loa loa (African eyeworm)

10

Wuchereria bancrofti (agent of lymphatic filariasis)

3

crustaceans

Onchocerca volvulus (agent of onchocerciasis)

5

Lepeophtheirus salmonis (salmon louse) strain Pacific

59 56

mites & ticks

Varroa destructor (honeybee mite)strain Korean

placozoans

Trichoplax adhaerens (placozoan) strain Grell-BS-1999

2

hydrozoans

Hydra magnipapillata (hydrozoan) strain 105

18

Total

35

302

3

* Some endogenous virus-like sequences from the indicated species have been reported by Belyi et al.[27] † The two capsid-like sequences were incorrectly reported as Rep-like sequences by Belyi et al. [27]

hosts over a long evolutionary time frame. For example, the virus-like sequences from lower eukaryotes (such as protozoans) were generally present at the base of each clade while those in relatively higher eukaryotes were more closely related to the known circoviruses, geminiviruses and nanoviruses that were infecting higher eukaryotes (see Additional file 1: Figure S4). There were,

however, several exceptions, possibly due to horizontal viral transfers over short periods of time. In most cases, the endogenous virus-like sequences from one species clustered together (such as those in salmon louse, honeybee mite, Hydra and roundworm species) (Additional file 1: Figure S4). Sequence comparison showed that, in each genome, some endogenous

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 6 of 15

&DW

'RJ ES

0DUNHU

&KU 

&KU 

&KU 

&KU 









  

Figure 2 Circovirus-like sequences in dog and cat genome were validated by PCR amplification and sequencing. PCR products were fractionated by gel electrophoresis on 1% agarose gels and stained with ethidium bromide. Marker, DNA marker DL 2000. The sequences of bands of the expected sizes from lanes: Chr 5, Chr 22, 071, and 274 were deposited under Genbank accession numbers: JF414126-JF414131.

viral copies may have resulted from segmental duplication within host genomes after a single original integration, as similar levels of identity are observed between them as well as between their flanking genomic regions. While others may have been derived from multiple independent integration events, as no similarity was observed among their flanking genomic sequences (see Additional file 1: Figure S6). Generally, the copy numbers of integrated viral sequences are less than 10 copies per species; whereas near sixty copies were identified in genomes of salmon louse (Lepeophtheirus salmonis) and honeybee mite (Varroa destructor) (see Additional file 1: Figure S4).

Comparison of the viral copies and their adjacent host sequences in these two genomes showed that for most viral copies, no similarity was observed among their flanking genomic sequences, suggesting that most were derived from multiple invasions of same or very similar viruses. However, considering that the Rep protein of eukaryotic ssDNA viruses has DNA binding, endonuclease and NTPase activity required for viral DNA replication [30,31], the integrated genes encoding Rep-like proteins may catalyze their own single-strand excision and invasion, and therefore act as selfish genetic elements capable of parasite-like proliferation in the host genome. This scenario could be supported by the fact

glucan (1,4-alpha-), branching enzyme 1 like gene 2800

NADH dehydrogenase subunit 2 gene

2552 bp

4004

NW_002106000.1

4243 4296

4463 bp

ADND01076753.1 2721 2716

88% nt

1

91% nt

1257 bp

74% nt

76% nt

ADND01038826.1 380

1648 bp

NW_002173830.1

Hydra magnipapillata (hydrozoan) strain 105 139611

140159

Lepeophtheirus salmonis strain Pacific 141860 bp

Pao retrotransposon

Transposable element Tcb2

1

NW_002157217.1

1193 bp

ADND01213500.1 2729

96% nt

95% nt

6000 bp

92% nt

1

NW_002166165.1

89% nt

90% nt

ADND01306175.1 99% nt

10968

1055

9370 bp

1284

1687

2165 bp

Lepeophtheirus salmonis strain Pacific

NW_002161614.1

Hydra magnipapillata (hydrozoan) strain 105 1

Helicase-like gene

97% nt 13900 bp

12400

1748 bp

ADBW01002987.1 1

90% nt

519

93% nt

2412 bp

1800

Pao retrotransposon

76% nt

ADBU01000376.1 3300 bp

ADBW01001434.1

Onchocerca volvulus (agent of onchocerciasis)

Loa loa (African eyeworm)

ADBU01002969.1

Figure 3 Genomic comparisons showing the endogenous viral sequences inserted into coding regions of host genes. Rectangular boxes with arrowheads indicate genes (Red, viral Rep-like genes; blue, host genes). Gray sectors connect corresponding homologous regions and the % nucleotide (nt) identity scores are indicated.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Cyclovirus PK5006 (ADD62451.1) Porcine circovirus 1 (NP_065678.1) Faba bean necrotic yellows virus DNA 1 (NP_619565.1) Faba bean necrotic stunt virus DNA R (YP_003104737.1) 96 Entamoeba terrapinae (AM688208.1:76..438) Entamoeba terrapinae (AM680040.1:12..239) 98 Entamoeba invadens IP1 (AANW02001352.1:3997..4345) 74 Entamoeba invadens IP1 (AANW02000727.1:5113..5313) Entamoeba invadens IP1 (AANW02001627.1:3316..4164) 72 Entamoeba invadens IP1 (AANW02001352.1:429..731) Protozoan Entamoeba invadens IP1 (AANW02003860.1:1242..1544) 74 97 Entamoeba invadens IP1 (AANW02004798.1:108..431) Entamoeba invadens IP1 (AANW02003978.1:1291..1473) Entamoeba invadens IP1 (AANW02000727.1:4414..4848) 91 Entamoeba invadens IP1 (AANW02001627.1:561..1343) Entamoeba invadens IP1 (AANW02002709.1:1267..1488) 84 95 Porphyra tenera plasmid pPT4-NU (YP_003058229.1) 52 Porphyra tenera plasmid pPT4 RAP (PT4) (AAL50514.1) Phytophthora sojae mitochondrion (YP_001165402.1) 72 88 66 Porphyra pulchra plasmid (AAF36423.1) Porphyra pulchra plasmid (AAF36422.1) 99 Porphyra pulchra plasmid (AAF36424.1) Onion yellows phytoplasma plasmid EcOYW1 (YP_002600745.1) 100 Onion yellows phytoplasma plasmid EcOYW2 (BAD36752.1) Onion yellows phytoplasma plasmid EcOYNIM (BAI39443.1) 76 83 Onion yellows phytoplasma plasmid EcOYM(P60470.1) 94 Onion yellows phytoplasma plasmid unknown (YP_214979.1) Peanut witches'-broom phytoplasma plasmid pPNWB (NP_852617.1) 100 Phytoplasma sp. plasmid pPSE (YP_001687728.1) 100 Phytoplasma sp. extrachromosomal DNA (CAB38248.1) 89 100 Aster yellows witches'-broom phytoplasma AYWB (YP_456464.1) 82 Aster yellows witches'-broom phytoplasma AYWB (YP_456347.1) Candidatus Phytoplasma australiense plasmid pPAPh2 (YP_001965305.1) 78 Candidatus Phytoplasma australiense plasmid pPASb11 (YP_001965310.1) 65 Aster yellows witches'-broom phytoplasma AYWB plasmid pAYWB-I (YP_456868.1) 91 Aster yellows witches'-broom phytoplasma AYWB plasmid pAYWB-III (YP_456877.1) 98 Paulownia witches'-broom phytoplasma plasmid pPaWBNy-2 (YP_001708790.1) Paulownia witches'-broom phytoplasma plasmid pPaWBNy-1 (YP_001708784.1) 96 Populus trichocarpa (NC_008478.1:11410621..11411228) Plants 98 Beet mild curly top virus-[Worland] (AAC54875.1) 89 71 Turnip curly top virus (YP_003778178.1) 86 Beet curly top virus (NP_040557.1) 100 Euphorbia mosaic virus (YP_717930.1) Tomato pseudo-curly top virus (NP_620735.1) Spilanthes yellow vein virus (YP_001285764.1) Kudzu mosaic virus (YP_001333680.1) 86 100 Emilia yellow vein virus-[Fz1] (YP_001661461.1) 98 Tobacco yellow dwarf virus (NP_620726.1) Wheat dwarf virus (CAJ13704.1) 100 Chloris striate mosaic virus (P18919.2) 76 Miscanthus streak virus (isolate 91) (Q67590.2) 99 Aspergillus niger CBS 513.88 (NT_166527.1:794829..795464) Trichoderma atroviride IMI 206040(ABDG01000062.1:195366..195950) 56 Laccaria bicolor S238N-H82 (ABFE01000943.1:16201..16343) Laccaria bicolor S238N-H82 (ABFE01000173.1:7779..8450) 77 Laccaria bicolor S238N-H82(ABFE01002205.1:68683..69683) Laccaria bicolor S238N-H82 (ABFE01001411.1:2070..2168) 53 Laccaria bicolor S238N-H82 (ABFE01001115.1:410..484) 96 97 Sclerotinia sclerotiorum hypovirulence associated DNA virus 1 (YP_003104796.1) Aspergillus nidulans FGSC A4 (NT_107012.1 chr II:1300422..1301042) 86 Nectria haematococca mpVI 77-13-4 (ACJF01000114.1:24025..24699) 98 100 Nectria haematococca mpVI 77-13-4 (ACJF01000185.1:14474..15148) 80 Nectria haematococca mpVI 77-13-4 (ACJF01000064.1 chr14: 36904..37572) Nectria haematococca mpVI 77-13-4 (ACJF01000097.1:8768..9433) Magnaporthe oryzae 70-15 (AACU02000431.1:90509..91234) 93 Aspergillus fumigatus A1163(ABDB01000026.1:88877..89998) 61 Tuber melanosporum (CABJ01002799.1:51549..52532)

Tuber melanosporum (CABJ01003429.1:11547..12725) Tuber melanosporum (CABJ01002197.1:7831..8979) Tuber melanosporum (CABJ01003737.1:11320..12522) Tuber melanosporum (CABJ01000406.1:11432..12644) Tuber melanosporum (CABJ01000392.1:95208..95691) 82 Tuber melanosporum (CABJ01000392.1:81806..82291) Tuber melanosporum (CABJ01001447.1:116854..117984) Tuber melanosporum (CABJ01002068.1:20593..21724) Tuber melanosporum (CABJ01000817.1:114207..115349) Tuber melanosporum (CABJ01001134.1:7905..9035) Tuber melanosporum (CABJ01001345.1:33204..34334) Tuber melanosporum (CABJ01002081.1:5684..6820) Tuber melanosporum (CABJ01002898.1:26046..27182) Tuber melanosporum (CABJ01000406.1:22484..22894) Tuber melanosporum (CABJ01000723.1:49567..50703) Tuber melanosporum (CABJ01002081.1:9116..10254) Tuber melanosporum (CABJ01002854.1:20374..21496) Tuber melanosporum (CABJ01001208.1:9066..10188) Tuber melanosporum (CABJ01002901.1:59..1238) Tuber melanosporum (CABJ01000040.1:17816..19087) Tuber melanosporum (CABJ01002167.1:63885..65020) Tuber melanosporum (CABJ01000204.1:8064..9077) Tuber melanosporum (CABJ01002567.1:4533..5669) Tuber melanosporum (CABJ01003913.1:41855..42991) Tuber melanosporum (CABJ01002683.1:2250..3390) Tuber melanosporum (CABJ01002337.1:61365..62501) Tuber melanosporum (CABJ01002800.1:7107..8243) Tuber melanosporum (CABJ01003136.1:44683..45729) Tuber melanosporum (CABJ01003747.1:130358..131494) Tuber melanosporum (CABJ01003956.1:1..1144) Tuber melanosporum (CABJ01001364.1:71248..71733) Tuber melanosporum (CABJ01000605.1:10854..11990) Tuber melanosporum (CABJ01001521.1:102536..103673) Tuber melanosporum (CABJ01000406.1:14806..15459) Tuber melanosporum (CABJ01000762.1:46366..47477) Tuber melanosporum (CABJ01000883.1:27559..28697) Tuber melanosporum (CABJ01001536.1:23875..24997) Tuber melanosporum (CABJ01003050.1:93469..95390) Tuber melanosporum (CABJ01000593.1:44072..45208) Tuber melanosporum (CABJ01000301.1:30477..31614) Tuber melanosporum (CABJ01004147.1:17387..18208)

Fungi

100

Page 7 of 15

84

98

93 99

*

90 86 69 85 76 77

100 81

Algal plasmid 97

Mitochondrion Phytoplasmal plasmid Phytoplasma Geminivirus Mycovirus Eukaryote

84 83

94 96 75

91 83 87 52 92

0.5

67 78 75

Fungi

97

Figure 4 Phylogeny of geminiviral Rep-like sequences from eukaryotes, known viruses, plasmids and phytoplasma. The phylogenetic tree was built using PhyML-mixtures based on a multiple sequence alignment generated using COBALT with the Constraint E-value parameter setting to 0.1. This tree was rooted with circoviruses and nanoviruses. The topology of blue asterisk marked clade was evaluated independently. Only p-values of the approximate likelihood ratios (SH-test) > 0.5 (50%) are indicated. scale bars correspond to 0.5 amino acid substitutions per site. Sequence accession numbers are given for each sequence.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

that the putative products of Helitrons [43], a eukaryotic rolling-circle transposon, shares motifs with the Reps of RCR plasmids and ssDNA viruses. Based on this fact, it has been suggested that ssDNA viruses might have evolved from RC transposons [43]. Our finding of endogenous viral Rep-like genes, however, favors the hypothesis that Helitrons may have arisen from ssDNA viruses which were integrated into the genome of an early eukaryotic ancestor [44]. Identification of ssDNA virus-like transposable elements in eukaryotic genomes

There are 42 geminivirus-like Rep genes or remnants interspersed in the genome of Perigord black truffle (Tuber melanosporum), an ectomycorrhizal fungus. All but one are most closely related to each other and formed a distinct clade (Figure 4). They share high (> 95%) nucleotide sequence identities with each other and thus allow us to reconstruct a consensus sequence. The reconstructed copy contains one interrupted Replike open reading frame (ORF), two transposase ORFs (one is interrupted and the other is truncated), and one microsatellite sequence (Figure 5A). It also contains 37-bp terminal inverted repeats (TIRs) but no obvious target site duplications (TSDs). It is most likely that this copy represents a novel transposon related to geminivirus identified in a eukaryotic genome. The genetic distances among these transposable repeats are very short suggesting that the transposons have undergone recent large-scale amplification in the host genome. In addition to the geminivirus-like transposon, we have also identified a parvovirus (linear ssDNA)-like repetitive element in the acorn worm (Saccoglossus kowalevskii) genome (see Additional file 2: Tabular data S3). Like parvoviruses, this repetitive element contains two large ORFs: one putative ORF encodes a protein containing parvovirus non-structural protein NS1 domain (Parvo_NS, pfam01057) and the other putative ORF encodes a protein containing parvovirus coat protein VP1 domain (Parvo_coat_N, pfam08398) (Figure 5B). It also possesses a palindromic hairpin structure at its 5’ terminal sequence, which is commonly found in parvoviruses. There are over 50 copies of this repeat interspersed in the genome. Some of these contained degraded ORFs; and some contained only a single ORF or a fragment. Furthermore, we also identified numerous parvovirus non-structural protein-like sequences in genomes of the hydrozoan Hydra magnipapillata and the planarian Schmidtea mediterranea (Figure 5B). We noted that these have been annotated as integrated virus-like element: DENSOV_HM and PIVE in Repbase Update, respectively [Jurka J, Repbase Reports 8(3), 182182 (2008), http://www.girinst.org/2008/vol8/issue3/

Page 8 of 15

DENSOV_HM.html; Rebrikov DV et al. Repbase Reports 8(2), 166-166 (2008), http://www.girinst.org/2008/vol8/ issue2/PIVE.html]. Phylogenetic analysis revealed that PIVE is more closely related to planarian (Girardia tigrina) virus, Planaria asexual strain-specific virus-like element (PEVE) (see Additional file 1: Figure S7), which has not yet been found to integrate in host genome [45]. The DENSOV_HM was located at the base of the papillomavirus clade and did not cluster within family Papillomaviridae. Furthermore, their genome structure is different. Therefore it may have originated from integration of the virus in new family infecting hydrozoan. The parvovirus-like sequences from acorn worm grouped together with PIVE_1p and PEVE small segment and these further clustered with parvoviruses. Considering that the genome structure of acorn worm repeated element is also similar to parvoviruses, it is most likely that it derived from parvovirus lineage infecting acorn worm distantly related to known parvoviruses. Consequently, these findings provide direct evidence that eukaryotic transposons could originated from ssDNA viruses. Preservation and expression of endogenous viral genes in host genomes

Examination of the potential coding capacity of endogenous viral sequences indicates that most were truncated and degraded, containing numerous premature stop codons, frameshift mutations, and insertions or deletions (see Additional file 2: Tabular data S2), suggesting that these are unlikely to have functions. Some endogenous viral sequences however, were found to encode uninterrupted ORFs. The conservation of relatively long ORFs suggests that they may have evolved under functional constraints since integration. While intact ORFs could also reflect recent insertion rather than functional maintenance in a long-standing history within the host genome. Analysis of transcription products can provide stronger evidence of functional maintenance. Accordingly, we used putative amino acid sequences of endogenous viral sequences to search NCBI Expressed Sequence Tags (ESTs) database for the corresponding mRNAs with the tBLASTn program. Through subsequent sequence comparisons, it was clear that some endogenous viral sequences are expressed as mRNAs in host genomes because they share high sequence identity with most of their ESTs and adjacent host sequences over full-length sequences (Figure 6 and Additional file 2: Tabular data S4). Interestingly, ESTs related to endogenous viral sequences of the rice blast fungus M. oryzae and two roundworm species were found in related species whose genome sequences are not currently available, suggesting that there were

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 9 of 15

GCCGTTGACAGCTAGTTAATTTTGTCCCCGGACAATA

transposable element tc3 transposase

microsatellite

Rep-like

1

A

4809 bp TIR

654

palindromic hairpin structure

Tuber melanosporum geminivirus-like repeated element consensus sequence

T T C A C T G G A G T G A A C T A A A G T T

GA T G A C C T C A C T T G A T T T C A A

1070106

Saccoglossus kowalevskii NW_003137822.1

1070593

contig gap (813bp)

Parvo coat_N

1075686

Parvo_NS1

1070188

B 5699 bp

ORF4L

ORF1L

34% aa

1

PIVE_1p

TIR

Parvo_NS1

PIVE_3p

TIR

ORF2S

22% 5953 bp

PIVE_2p

33% aa

ORF1S TIR

ORF3L

Planaria asexual strain-specific virus-like element type 1 large DNA segment NC_003220.1 TIR 1 ORF2L

TIR

1

TIR

4392

99% QW

T A C G G C C G G C G C G C C G T A A T A T A T G C

G T GC C A T C GT G C A A

4077

969

1911 bp

Schmidtea mediterranea integrated virus-like element (PIVE)

TIR

Planaria asexual strain-specific virus-like element type 1 small DNA segment NC_003219.1 DENSOV_HM_2p

1

DENSOV_HM_1p

5133 bp

PPV_E1_C

Hydra magnipapillata integrated densovirus-like element (DENSOV_HM) Figure 5 Genomic organization of ssDNA virus-like transposons in fungi (A) and lower eukaryotes (B). (A) The genomic organization of geminivirus-like transposon in Tuber melanosporum. Arrowhead boxes indicate ORFs (orange, Rep-like gene; blue, transposase gene). The black vertical lines in the arrowhead boxes indicate stop codons. Green rectangular box indicates microsatellite sequence. The sequence of terminal inverted repeat (TIR) is shown at the top to the right. (B) The genomic organization and comparison of parvovirus-like transposon with related exogenous planaria virus. Yellow arrowhead boxes indicate Rep-like ORFs. Swallow tails indicate terminal inverted repeats (TIRs). The annotated ORF names are indicated. Purple rectangular boxes indicate protein domains and the domain family names are shown: Parvo_NS1, Parvovirus non-structural protein NS1 (pfam01057); Parvo_coat_N, Parvovirus coat protein VP1 (pfam08398); PPV_E1_C, Papillomavirus helicase (pfam00519). Gray sectors connect corresponding homologous regions and the % nucleotide (nt) or amino acid (aa) identity are indicated. The Planaria asexual strain-specific virus-like element has not been found to integrate in the host genome.

similarly expressed insertions in these species. In fact, we have also amplified these virus-like sequences from different strains of the rice blast fungus and detected their expression (data not shown). The fact that endogenous viral sequences are conserved and expressed in host organisms suggests that these viral genes have been coopted to assume cellular functions in eukaryotic genomes. It should be noted that expression of mRNA from endogenous viral sequences was also detected in the parasitic protozoan

E. histolytica, although their long ORFs were defective (Figure 6). Perhaps selection to maintain these viral sequences has recently been lost. We also detected endogenous viral sequence-related ESTs in some plants and animal species (see Additional file 2: Tabular data S5). Because the genome sequences of relevant host species are not available or available but not well matched with their ESTs, it remains to be established whether they represent authentic expressed endogenous viral sequences.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 10 of 15

Phaeodactylum tricornutum (diatom) ABQD01000138.1

Entamoeba histolytica HM-1:IMSS AAFB02000406.1 6730

7617

8149

98%

CX094342

E.histolytica cDNA

P. tricornutum cDNA

100%

Entamoeba histolytica HM-1:IMSS AAFB02000331.1 7333

6825

7801

27102 bp

26185

8835 bp

8958

9229

Entamoeba histolytica HM-1:IMSS AAFB02000468.1

10075 bp

5480 bp

4594

XP_652804 99%

100%

E.histolytica cDNA

Onchocerca volvulus (agent of onchocerciasis) ADBW01007856.1

E.histolytica cDNA

1049 1077 bp

667

1 Ancylostoma caninum EX559048 (dog hookworm) cDNA

78%

Wuchereria bancrofti (agent of lymphatic filariasis) ADBV01018817.1 1

1043 1091 bp

153

Magnaporthe oryzae 70-15 AACU02000431.1 90506

BQ457338 BQ456255

>66%

Dirofilaria immitis (dog heartworm nematode) cDNAs

FF140597 FF141572 Onchocerca flexuosa (deer filarial worm ) cDNAs

91234 bp

DC966977 DC962898

>97%

DC962897 CD030328 DC966976

M. grisea cDNA

Figure 6 Organization and transcription maps of endogenous viral Rep-like genes. Blue arrowhead boxes indicate Rep-like ORFs. Similar regions of expressed sequences are identified and the % nt identity with endogenous viral sequences are indicated. Note that the actual endogenous viral sequence extended beyond the ORF in two contigs (AAFB02000406.1 and AAFB02000468.1) of Entamoeba histolytica.

The host range of circular ssDNA viruses

Circoviruses are previously known to infect only birds and pigs [46]. These viruses have been detected in dragonflies, fish and human most recently [47-49]. Geminiviruses and nanoviruses only have been found to infect higher plants [23,50]. Recent metagenomic studies uncovered that these circular ssDNA viruses were commonly found in various environmental samples, but it is difficult to provide information on the host range and ecology for these viruses. The endogenous circovirus-like sequences in honeybee mite were most closely related to cycloviruses (Additional file 1: Figure S4), members of a recently proposed genus in the family Circoviridae, which were commonly found in faecal samples of human and chimpanzee by viral metagenomics [20]. In addition, the endogenous virus-like sequences from some species of various organisms (such as amphibians, algae, diatom, gastropod, etc.) were clustered with viral metagenomic sequences or circovirus-like genomes identified from environmental samples [18]. These findings suggest that these various species are the definitive hosts of relevant environmental viruses.

Altogether, we discovered endogenous ssDNA viruslike sequences in host species broadly distributed in four of the five supergroups of eukaryotes [51], namely Unikonta (including animals, fungi and Entamoeba), Plantae (including land plants and green algae), Chromalveolata (including diatoms and Blastocystis) and Excavata (including Giardia) (Figure 7). The endogenous viral sequences identified here as molecular fossil records of past viral invasions provide unambiguous definitive hosts for these viruses and extend the host range of circular ssDNA viruses. Interestingly, although the endogenous circovirus/ nanovirus-like sequences occurred widely in the genomes of eukaryotic species, ranging from unicellular organisms to mammals, we did not detect any of these sequences in plant, bird and pig genomes sequenced todate. In contrast, the nematodes (roundworms) were not known to be infected by ssDNA virus, but endogenous circovirus-like sequences occurred in some nematode species. In addition, geminivirus-like sequences were found in some fungal genomes. So far, however no genetically related exogenous counterparts were found in these fungi, even though some of these, such as the

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 11 of 15

Exogenous virus Endogenous virus

Carnivores Ungulates Primates

Mammals

Placentals Marsupials

Vertebrates

Birds Amphibians Bony fishes Lancelets Dragonflies

Arthropods

Unikonts

Crustaceans Mites

Animals

Nematodes Gastropods Hydrozoans Placozoans Fungi

Ascomycetes Basidiomycetes Entamoeba

Green plants

Green algae Land plants Diatoms

Plantae

Blastocystis

Chromalveolates

Giardia

Excavates

Figure 7 A tree of eukaryotes showing the known distribution of endogenous viral-like sequences and exogenous circo-, nano- and geminivirus like viruses. This tree was drawn base on The Tree of Life Web Project (http://tolweb.org/).

rice blast fungus Magnaporthe oryzae, were widely studied. Likewise, Populus is not known to be infected by geminiviruses but harbored one endogenous geminiviral sequence. These observations suggest that some of the

endogenous viral sequences could provide immune protection in the host similar to the endogenous retroviral capsid proteins in mice and sheep, which offer protection against exogenous retroviral infections [52,53].

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 12 of 15

The origin and evolution of circoviruses and nanoviruses

The origin and evolution of geminiviruses

Based on the different phylogenies between the N-terminal and the C-terminal regions of circovirus Rep, Gibbs and Weiller [54] suggested that circovirus Rep proteins may have evolved by a recombination event between the Rep protein of nanoviruses and an RNA binding protein encoded by picorna-like viruses after the nanoviruses switched hosts to infect a vertebrate. However, it seems unlikely that the virus recombination event took place in a vertebrate considering the fact that endogenous circovirus-like sequences were widely found in nonvertebrate species. To examine more thoroughly the origin and evolution of circoviruses and nanoviruses, we selected representative Rep-like proteins from viruses, plasmids and bacterial genomes and used sufficient samples to construct phylogenetic trees. In consideration of a possible recombination event, we aligned and performed phylogenetic analysis corresponding to fulllength Rep genes, the N-terminal and C-terminal regions respectively. As shown in Additional file 2: Figure S8, circovirus-like sequences and viral Rep-like sequences from bacterial plasmid and bacterial genomes clustered together in all trees. However, while nanovirus-like sequences clustered with circovirus-like sequences in the N-terminal tree they were grouped with geminivirus-like sequences in the C-terminal tree. In the full-length Rep tree, nanovirus-like sequences were placed between the geminivirus-like and circovirus-like sequences, possibly due to the compromise of different phylogenetic signals from the two parts of nanovirus-like Reps. Therefore, if a recombination event had occurred, it is likely to have taken place in the nanovirus-like Reps rather than in the circovirus-like Reps. It has been proposed that eukaryotic ssDNA viruses may have evolved from prokaryotic plasmids or phages [31]. In our phylogenetic trees, the virus-like sequences from bacterial plasmid and bacterial genomes were generally located at the base of circovirus-like sequences, suggesting that circoviruses might have originated from relevant bacterial plasmids. Considering that the nanovirus-like sequences clustered with circovirus-like sequences in the N-terminal tree, it is most likely that the nanoviruses shared the most recent common ancestor with circovirus-like viruses and subsequently the Cterminal sequences of ancestor nanoviral Reps may have recombined with those of geminivirus-like viruses or plasmids. But the possibility that nanovirus-like Reps were the result of convergent evolution cannot be ruled out. The Canarypox virus and the ancestor of picornalike viruses may have captured the helicase domain sequences from circovirus-like viruses by recombination.

Based on the observations that geminiviruses occupied a common ecological niche with phytoplasmas and their Reps shared a most recent common ancestor with phytoplasmal plasmids in phylogenetic analysis, Krupovic et al [55] proposed that the geminiviruses may have originated from phytoplasmal plasmid followed by acquisition of the capsid gene from an ssRNA plant virus. However, in view of the recent reports on the geminivirus-like mycovirus and numerous related sequences in fungal genomes, the evolutionary relationships among these geminivirus-like elements need to be revaluated. To address this question, we constructed phylogenetic trees for the representative Rep-like proteins from plants, fungi, phytoplasma and algae using the fulllength Rep genes, the N-terminal RCR catalytic domain and C-terminal helicase domain sequences respectively (Additional file 2: Figure S9). In all trees, the plant geminiviral Reps clustered together with fungal Reps, suggesting that they shared a more recent common ancestor with those from fungi rather than from phytoplasmal plasmids. Furthermore, although the Rep protein of SsHADV-1 is related to geminiviruses, the genome organization of SsHADV-1 and particle morphology is distinct from those of geminiviruses [42]. Although the capsid protein of SsHADV-1 lacks sequence similarity with those of any geminiviruses, its most similar sequences are commonly found in environmental samples. In addition, sequences related to SsHADV-1 were widely found in fungal genomes and diverse metagenomic samples. These results suggest that SsHADV-1 and related viruses from fungi and environment may have evolved independently rather than being descendent from geminiviruses or vice versa. Therefore, it is possible that the ancestor of geminiviruses and related fungal viruses may have occurred prior to the separation of plants and fungi, and subsequently perhaps they had a unique path to evolution in their hosts.

Conclusions Our study provided comprehensive and convincing evidence that the genes of small circular ssDNA viruses have been transferred into a broad range of eukaryotic genomes, and some of the transferred genes were also conserved and functional in host genomes. This discovery extends the host range of circular ssDNA viruses and offers insight into the origin and evolution of relevant viruses. Furthermore, our finding also revealed that the capture and functional assimilation of exogenous viral genes may represent an important force in the evolution of eukaryotes.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Page 13 of 15

Methods

Sequence comparison and phylogenetic analysis

Genome screening

The putative peptides of endogenous viral sequences were obtained according to BLASTx hits and manual editing. The in-frame stop codons were indicated as X. Multiple alignments of protein sequences were constructed either using MCOFFEE (when the number of sequences < 50) [56] or using COBALT [57] (http:// www.ncbi.nlm.nih.gov/tools/cobalt/cobalt.cgi?link_loc=BlastHomeAd) and manually edited. To give the best alignment, the alignment parameter Constraint E-value and Word Size were adjusted for different datasets when using COBALT. Although many of the endogenous viral sequences are of different lengths in alignments, it is now well known that sequences of very different lengths can be accurately placed on phylogenies [58]. Hence, all the putative peptides of endogenous viral sequences were used for the phylogenetic analysis with proteins of representative exogenous viruses to determine the full-scale evolutionary relationships among them. Maximum likelihood (ML) phylogenies were estimated using amino acid sequence alignments with PhyML-mixtures [59,60], assuming the EX2 mixture model [60] and SPR tree topologies search strategy [61]. Gaps in alignment are systematically treated as unknown characters. The reliability of internal branches was evaluated based on approximate likelihood ratio test (aLRT) statistics [62].

In order to screen for the circular ssDNA virus-related sequences in eukaryotic genomes, we performed tBLASTn searches against different NCBI sequence databases (http://blast.ncbi.nlm.nih.gov/Blast.cgi) using as queries the representative peptide sequences derived from viruses in families Anelloviridae, Circoviridae, Geminiviridae and Nanoviridae. NCBI databases used for sequence searches included nr (all GenBank + RefSeq Nucleotides + EMBL + DDBJ + PDB sequences + HTGS phase 3 but excluding HTGS phase 0,1,2, EST, GSS, STS, PAT, WGS), refseq_genomic (genomic entries from NCBI’s Reference Sequence project), NCBI Genomes/chromosome (a database with complete genomes and chromosomes from the NCBI Reference Sequence project), wgs (a database for whole genome shotgun sequence entries), gss (Genome Survey Sequence, includes single-pass genomic data, exontrapped sequences, and Alu PCR sequences), htgs (unfinished High Throughput Genomic Sequences: phases 0, 1 and 2), and the eukaryotes genomic BLAST database. All non-redundant matches from these searches with E-values ≤1e-5 were extracted along with 1 kb of flanking regions, and then were used to screen the non-redundant (NR) protein database using BLASTx. All genomic sequences from host genomes that unambiguously matched viral proteins were considered as candidate endogenous viral sequences. These candidate endogenous viral sequences were used to research the databases for other homologous sequences that would have been missed during initial searches using the known extant viruses. All database searches were performed online and were completed in June 2010. Examining possible chimeras or errors in assembling of endogenous viral sequences

To rule out the possibility that these endogenous viral sequences were chimeric clones or misassembled from contaminated sequences of exogenous incidental viral sequences, we searched against archival data of the eukaryotic genome sequencing using their endogenous viral sequences and flanking cellular sequences as megaBLAST queries on the NCBI Trace Archive (http://www.ncbi.nlm.nih.gov/blast/mmtrace. shtml) with the cut-off value: > 95% nt identity, respectively; and carefully examined the junctions between endogenous viral sequences and cellular sequences. The statistics of junction coverages that show the number of trace records containing the junctions between endogenous viral sequences and cellular sequences are listed in Additional file 2: Tabular data S2.

Detection of expression of endogenous viral sequences from host genomes

To investigate whether endogenous viral sequences could be expressed in host genomes, we first, used the endogenous viral sequences to screen the NCBI EST database using the method described in Genome screening. Subsequently, we used the identified virus-related ESTs to compare with host genomes and virus genomes by megaBLAST to determine whether they were expressed sequences from host genomes or the result of laboratory contamination. PCR amplification and DNA sequencing

Genomic DNA samples of dog (Canis lupus familiaris) and cat (Felis catus) were obtained from Zyagen Laboratories (USA). To PCR amplify the candidate DNA fragments from these DNA samples, primers pairs were designed based on the virus-like sequences and their flanking cellular sequences, see Additional file 3: Table S1 for the primers pairs used. PCR products were fractionated by gel electrophoresis on 1% agarose gels and stained with ethidium bromide. DNA was sequenced by Sanger methods at the Beijing Genomics Institute (BGI). New sequences generated in this study were deposited in the GenBank under accession numbers: JF414126JF414131.

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

Additional material

Page 14 of 15

2. 3.

Additional file 1: supplementary figures. This file includes 9 supplementary figures. Figure S1 illustrates the domain organization of different Rep-like proteins. Figure S2 and S3 show multiple alignments of circovirus, nanovirus or geminivirus Rep-like sequences, respectively. Figure S4 shows the phylogeny of viral Rep-like sequences from eukaryotes, known viruses and viral metagenomes. Figure S5 shows the phylogeny of geminiviral Rep-like sequences in viral metagenomes. Figure S6 shows the alignment of viral insertion loci in a genome. Figure S7 shows the phylogeny of parvovirus-like transposons. Figure S8 and S9 show the phylogenies of full-length, N-terminal and C-terminal regions of circoviral or geminiviral Rep-like proteins, respectively. Additional file 2: supplementary Tabular data. This file contains 5 supplementary Tabular data. Tabular data S1 shows the characterization of eukaryotic circular ssDNA viral Rep-like proteins in other systems. Tabular data S1 gives results of all detected endogenous circular ssDNA virus-like sequences in this study. Tabular data S3 shows the endogenous parvovirus-like sequences in acorn worm genome. Tabular data S4 summarizes the expressed endogenous viral sequences. Tabular data S5 lists the circular ssDNA virus-related ESTs which have not been determined if they were expressed endogenous viral sequences. Additional file 3: supplementary table S1. This file lists primers used for PCR of endogenous virus-like regions of dog and cat genomes.

4.

5.

6. 7.

8. 9.

10. 11.

List of abbreviations aLRT: approximate likelihood ratio test; BLAST: Basic Local Alignment Search Tool; ESTs: Expressed Sequence Tags; CP: capsid protein; GSS: Genomic Survey Sequence; HTGS: High Throughput Genomic Sequence; HGT: horizontal gene transfer; ICTV: International Committee on Taxonomy of Viruses; NCBI: National Center for Biotechnology Information; ORF: open reading frame; PCR: Polymerase chain reaction; RCR: rolling circle replication; ssDNA: single-stranded DNA; SPR: subtree prune and regraft; TIRs: terminal inverted repeats; TSDs: target site duplications; Rep: replication initiator protein; TEs: transposable elements; WGS: Whole Genome Shotgun. Acknowledgements We thank the anonymous reviewers for their constructive and helpful comments. This research was supported in part by the Program for New Century Excellent Talents in University (NCET-06-0665), the Commonweal Specialized Research Fund of China Agriculture (3-21) and the Huazhong Agricultural University Scientific & Technological Self-innovation Foundation. Author details State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan 430070, Hubei Province, P R China. 2The Provincial Key Lab of Plant Pathology of Hubei Province, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei Province, P R China. 3Current Address: Purdue-NWAFU Joint Research Center, Northwest A&F University, Yangling, 712100, Shaanxi Province, P R China. 4 Department of Plant Pathology, University of Kentucky, 201F Plant Science Building, 1405 Veterans Drive, University of Kentucky, Lexington, KY 405460312, USA. 1

12.

13. 14. 15. 16.

17.

18.

19.

20.

21. 22.

Authors’ contributions HL, YF and DJ conceived and designed the study; HL performed the computational analyses and lab experiments; HL, YF, DJ, BL, XY, SAG, GL, JX, JC and XY analyzed data; and HL, YF, DJ and SAG wrote the paper. All authors read and approved the final manuscript.

24.

Received: 28 March 2011 Accepted: 26 September 2011 Published: 26 September 2011

25.

References 1. Canchaya C, Proux C, Fournous G, Bruttin A, Brussow H: Prophage genomics. Microbiol Mol Biol Rev 2003, 67(2):238-276.

23.

26.

Casjens S: Prophages and bacterial genomics: what have we learned so far? Mol Microbiol 2003, 49(2):277-300. Belshaw R, Pereira V, Katzourakis A, Talbot G, Paces J, Burt A, Tristem M: Long-term reinfection of the human genome by endogenous retroviruses. Proc Natl Acad Sci USA 2004, 101(14):4894-4899. Dunlap KA, Palmarini M, Varela M, Burghardt RC, Hayashi K, Farmer JL, Spencer TE: Endogenous retroviruses regulate periimplantation placental growth and differentiation. Proc Natl Acad Sci USA 2006, 103(39):14390-14395. Dupressoir A, Vernochet C, Bawa O, Harper F, Pierron G, Opolon P, Heidmann T: Syncytin-A knockout mice demonstrate the critical role in placentation of a fusogenic, endogenous retrovirus-derived, envelope gene. Proc Natl Acad Sci USA 2009, 106(29):12127-12132. Staginnus C, Richert-Poggeler KR: Endogenous pararetroviruses: two-faced travelers in the plant genome. Trends Plant Sci 2006, 11(10):485-491. Bezier A, Annaheim M, Herbiniere J, Wetterwald C, Gyapay G, BernardSamain S, Wincker P, Roditi I, Heller M, Belghazi M, et al: Polydnaviruses of braconid wasps derive from an ancestral nudivirus. Science 2009, 323(5916):926-930. Taylor DJ, Bruenn J: The evolution of novel fungal genes from nonretroviral RNA viruses. BMC Biol 2009, 7:88. Horie M, Honda T, Suzuki Y, Kobayashi Y, Daito T, Oshida T, Ikuta K, Jern P, Gojobori T, Coffin JM, et al: Endogenous non-retroviral RNA virus elements in mammalian genomes. Nature 2010, 463(7277):84-87. Taylor DJ, Leach RW, Bruenn J: Filoviruses are ancient and integrated into mammalian genomes. BMC Evol Biol 2010, 10:193. Belyi VA, Levine AJ, Skalka AM: Unexpected inheritance: multiple integrations of ancient bornavirus and ebolavirus/marburgvirus sequences in vertebrate genomes. PLoS Pathog 2010, 6(7):e1001030.. Liu H, Fu Y, Jiang D, Li G, Xie J, Cheng J, Peng Y, Ghabrial SA, Yi X: Widespread horizontal gene transfer from double-stranded RNA viruses to eukaryotic nuclear genomes. J Virol 2010, 84(22):11876-11887. Kapoor A, Simmonds P, Lipkin WI: Discovery and characterization of mammalian endogenous parvoviruses. J Virol 2010, 84(24):12628-12635. Gilbert C, Feschotte C: Genomic fossils calibrate the long-term evolution of hepadnaviruses. PLoS Biol 2010, 8:(9). Katzourakis A, Gifford RJ: Endogenous viral elements in animal genomes. PLoS Genet 2010, 6(11):e1001191.. Nagasaki K, Tomaru Y, Takao Y, Nishida K, Shirai Y, Suzuki H, Nagumo T: Previously unknown virus infects marine diatom. Appl Environ Microbiol 2005, 71(7):3528-3535. Tomaru Y, Shirai Y, Suzuki H, Nagumo T, Nagasaki K: Isolation and characterization of a new single-stranded DNA virus infecting the cosmopolitan marine diatom Chaetoceros dehilis. Aquat Microb Ecol 2008, 50(2):103-112. Rosario K, Duffy S, Breitbart M: Diverse circovirus-like genome architectures revealed by environmental metagenomics. J Gen Virol 2009, 90:(Pt 10):2418-2424. Kim KH, Chang HW, Nam YD, Roh SW, Kim MS, Sung Y, Jeon CO, Oh HM, Bae JW: Amplification of uncultured single-stranded DNA viruses from rice paddy soil. Appl Environ Microbiol 2008, 74(19):5975-5985. Li L, Kapoor A, Slikas B, Bamidele OS, Wang C, Shaukat S, Masroor MA, Wilson ML, Ndjango JB, Peeters M, et al: Multiple diverse circoviruses infect farm animals and are commonly found in human and chimpanzee feces. J Virol 2010, 84(4):1674-1682. Rosario K, Nilsson C, Lim YW, Ruan Y, Breitbart M: Metagenomic analysis of viruses in reclaimed water. Environ Microbiol 2009, 11(11):2806-2820. Lopez-Bueno A, Tamames J, Velazquez D, Moya A, Quesada A, Alcami A: High diversity of the viral community from an Antarctic lake. Science 2009, 326(5954):858-861. Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA: Virus Taxonomy: Eighth Report of the International Committee on Taxonomy of Viruses. San Diego: Elsevier Academic Press; 2005. Campos-Olivas R, Louis JM, Clerot D, Gronenborn B, Gronenborn AM: The structure of a replication initiator unites diverse aspects of nucleic acid metabolism. Proc Natl Acad Sci USA 2002, 99(16):10310-10315. Bejarano ER, Khashoggi A, Witty M, Lichtenstein C: Integration of multiple repeats of geminiviral DNA into the nuclear genome of tobacco during evolution. Proc Natl Acad Sci USA 1996, 93(2):759-764. Gibbs MJ, Smeianov VV, Steele JL, Upcroft P, Efimov BA: Two families of rep-like genes that probably originated by interspecies recombination

Liu et al. BMC Evolutionary Biology 2011, 11:276 http://www.biomedcentral.com/1471-2148/11/276

27.

28. 29. 30.

31.

32. 33.

34.

35.

36.

37.

38. 39.

40. 41.

42.

43. 44.

45.

46.

47.

48.

are represented in viral, plasmid, bacterial, and parasitic protozoan genomes. Mol Biol Evol 2006, 23(6):1097-1100. Belyi VA, Levine AJ, Skalka AM: Sequences from ancestral single-stranded DNA viruses in vertebrate genomes: the parvoviridae and circoviridae are more than 40 to 50 million years old. J Virol 2010, 84(23):12458-12462. Koonin EV, Senkevich TG, Dolja VV: The ancient Virus World and evolution of cells. Biol Direct 2006, 1:29. Forterre P: Evolution, viral. In Encyclopedia of Microbiology.. 3 edition. Edited by: Schaechter M. Oxford: Elsevier; 2009:370-389. Gorbalenya AE, Koonin EV, Wolf YI: A new superfamily of putative NTPbinding domains encoded by genomes of small DNA and RNA viruses. FEBS Lett 1990, 262(1):145-148. Ilyina TV, Koonin EV: Conserved sequence motifs in the initiator proteins for rolling circle DNA replication encoded by diverse replicons from eubacteria, eucaryotes and archaebacteria. Nucleic Acids Res 1992, 20(13):3279-3285. Khan SA: Rolling-circle replication of bacterial plasmids. Microbiol Mol Biol Rev 1997, 61(4):442-455. Oshima K, Kakizawa S, Nishigawa H, Kuboyama T, Miyata S, Ugaki M, Namba S: A plasmid of phytoplasma encodes a unique replication protein having both plasmid- and virus-like domains: clue to viral ancestry or result of virus/plasmid recombination? Virology 2001, 285(2):270-277. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389-3402. Martin FN, Bensasson D, Tyler BM, Boore JL: Mitochondrial genome sequences and comparative genomics of Phytophthora ramorum and P. sojae. Curr Genet 2007, 51(5):285-296. Avila-Adame C, Gomez-Alpizar L, Zismann V, Jones KM, Buell CR, Ristaino JB: Mitochondrial genome sequences and molecular evolution of the Irish potato famine pathogen, Phytophthora infestans. Curr Genet 2006, 49(1):39-46. Na S, Shen T, Jia P, Men D, Chen Q: Characterization of the natural deletion mutant of plasmid pXZ10145 in Corynebacterium glutamicum and construction of a recombinant plasmid. Chin J Biotechnol 1991, 7(4):271-277. Tulman ER, Afonso CL, Lu Z, Zsak L, Kutish GF, Rock DL: The genome of canarypox virus. J Virol 2004, 78(1):353-366. Franzen O, Jerlstrom-Hultqvist J, Castro E, Sherwood E, Ankarklev J, Reiner DS, Palm D, Andersson JO, Andersson B, Svard SG: Draft genome sequencing of giardia intestinalis assemblage B isolate GS: is human giardiasis caused by two different species? PLoS Pathog 2009, 5(8): e1000560.. Okamoto H: TT viruses in animals. Curr Top Microbiol Immunol 2009, 331:35-52. Murad L, Bielawski JP, Matyasek R, Kovarik A, Nichols RA, Leitch AR, Lichtenstein CP: The origin and evolution of geminivirus-related DNA sequences in Nicotiana. Heredity 2004, 92(4):352-358. Yu X, Li B, Fu Y, Jiang D, Ghabrial SA, Li G, Peng Y, Xie J, Cheng J, Huang J, et al: A geminivirus-related DNA mycovirus that confers hypovirulence to a plant pathogenic fungus. Proc Natl Acad Sci USA 2010, 107(18):8387-8392. Kapitonov VV, Jurka J: Rolling-circle transposons in eukaryotes. Proc Natl Acad Sci USA 2001, 98(15):8714-8719. Feschotte C, Wessler SR: Treasures in the attic: rolling circle transposons discovered in eukaryotic genomes. Proc Natl Acad Sci USA 2001, 98(16):8923-8924. Rebrikov DV, Bulina ME, Bogdanova EA, Vagner LL, Lukyanov SA: Complete genome sequence of a novel extrachromosomal virus-like element identified in planarian Girardia tigrina. BMC Genomics 2002, 3(1):15.. Mankertz A: Circoviruses. In Encyclopedia of Virology. Volume 1.. 3 edition. Edited by: Mahy BWJ, van Regenmortel MHV. Oxford: Elsevier; 2008:513-519. Rosario K, Marinov M, Stainton D, Kraberger S, Wiltshire EJ, Collings DA, Walters M, Martin DP, Breitbart M, Varsani A: Dragonfly cyclovirus, a novel single-stranded DNA virus discovered in dragonflies (Odonata: Anisoptera). J Gen Virol 2011, 92(Pt 6):1302-1308. Lorincz M, Csagola A, Farkas SL, Szekely C, Tuboly T: First detection and analysis of a fish circovirus. J Gen Virol 2011, 92(Pt 8):1817-1821.

Page 15 of 15

49. Sauvage V, Cheval J, Foulongne V, Gouilh MA, Pariente K, Manuguerra JC, Richardson J, Dereure O, Lecuit M, Burguiere A, et al: Identification of the first human gyrovirus, a virus related to chicken anemia virus. J Virol 2011, 85(15):7948-7950. 50. Vetten HJ: Nanoviruses. In Encyclopedia of Virology. Volume 2.. 3 edition. Edited by: Mahy BWJ, van Regenmortel MHV. Oxford: Elsevier; 2008:385-391. 51. Keeling PJ, Burger G, Durnford DG, Lang BF, Lee RW, Pearlman RE, Roger AJ, Gray MW: The tree of eukaryotes. Trends Ecol Evol 2005, 20(12):670-676. 52. Best S, Le Tissier P, Towers G, Stoye JP: Positional cloning of the mouse retrovirus restriction gene Fv1. Nature 1996, 382(6594):826-829. 53. Arnaud F, Murcia PR, Palmarini M: Mechanisms of late restriction induced by an endogenous retrovirus. J Virol 2007, 81(20):11441-11451. 54. Gibbs MJ, Weiller GF: Evidence that a plant virus switched hosts to infect a vertebrate and then recombined with a vertebrate-infecting virus. Proc Natl Acad Sci USA 1999, 96(14):8022-8027. 55. Krupovic M, Ravantti JJ, Bamford DH: Geminiviruses: a tale of a plasmid becoming a virus. BMC Evol Biol 2009, 9:112. 56. Moretti S, Armougom F, Wallace IM, Higgins DG, Jongeneel CV, Notredame C: The M-Coffee web server: a meta-method for computing multiple sequence alignments by combining alternative alignment methods. Nucleic Acids Res 2007, , 35 Web Server: W645-648. 57. Papadopoulos JS, Agarwala R: COBALT: constraint-based alignment tool for multiple protein sequences. Bioinformatics 2007, 23(9):1073-1079. 58. Wiens JJ: Missing data and the design of phylogenetic analyses. J Biomed Inform 2006, 39(1):34-42. 59. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003, 52(5):696-704. 60. Le SQ, Lartillot N, Gascuel O: Phylogenetic mixture models for proteins. Philos Trans R Soc Lond B Biol Sci 2008, 363(1512):3965-3976. 61. Hordijk W, Gascuel O: Improving the efficiency of SPR moves in phylogenetic tree search methods based on maximum likelihood. Bioinformatics 2005, 21(24):4338-4347. 62. Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol 2006, 55(4):539-552. doi:10.1186/1471-2148-11-276 Cite this article as: Liu et al.: Widespread Horizontal Gene Transfer from Circular Single-stranded DNA Viruses to Eukaryotic Genomes. BMC Evolutionary Biology 2011 11:276.

Submit your next manuscript to BioMed Central and take full advantage of: • Convenient online submission • Thorough peer review • No space constraints or color figure charges • Immediate publication on acceptance • Inclusion in PubMed, CAS, Scopus and Google Scholar • Research which is freely available for redistribution Submit your manuscript at www.biomedcentral.com/submit