Horizontal gene transfer and nucleotide compositional anomaly in large DNA viruses

The GC% landscape of the Nucleocytoviricota

Genomic studies on sequence composition employ various approaches, such as calculating the proportion of guanine and cytosine within a given sequence (GC% content), which can shed light on various aspects of the organism's biology. In this context, GC% can provide insights into virus-host relationships and evolution. Here, we present a comprehensive gene-by-gene analysis of 61 representatives belonging to the phylum Nucleocytoviricota, which comprises viruses with the largest genomes known in the virosphere. Parameters were evaluated not only based on the average GC% of a given viral species compared to the entire phylum but also considering gene position and phylogenetic history. Our results reveal that while some families exhibit similar GC% among their representatives (e.g., Marseilleviridae), others such as Poxviridae, Phycodnaviridae, and Mimiviridae have members with discrepant GC% values, likely reflecting adaptation to specific biological cycles and hosts. Interestingly, certain genes located at terminal regions or within specific genomic clusters show GC% values distinct from the average, suggesting recent acquisition or unique evolutionary pressures. Horizontal gene transfer and the presence of potential paralogs were also assessed in genes with the most discrepant GC% values, indicating multiple evolutionary histories. Taken together, to the best of our knowledge, this study represents the first global and gene-by-gene analysis of GC% distribution and profiles within genomes of Nucleocytoviricota members, highlighting their diversity and identifying potential new targets for future studies.

LINE-1 retrotransposons facilitate horizontal gene transfer into poxviruses

There is ample phylogenetic evidence that many critical virus functions, like immune evasion, evolved by the acquisition of genes from their hosts through horizontal gene transfer (HGT). However, the lack of an experimental system has prevented a mechanistic understanding of this process. We developed a model to elucidate the mechanisms of HGT into vaccinia virus, the prototypic poxvirus. All identified gene capture events showed signatures of long interspersed nuclear element-1 (LINE-1)-mediated retrotransposition, including spliced-out introns, polyadenylated tails, and target site duplications. In one case, the acquired gene integrated together with a polyadenylated host U2 small nuclear RNA. Integrations occurred across the genome, in some cases knocking out essential viral genes. These essential gene knockouts were rescued through a process of complementation by the parent virus followed by nonhomologous recombination during serial passaging to generate a single, replication-competent virus. This work links multiple evolutionary mechanisms into one adaptive cascade and identifies host retrotransposons as major drivers for virus evolution.

Virus-host coexistence in phytoplankton through the genomic lens

Virus-microbe interactions in the ocean are commonly described by "boom and bust" dynamics, whereby a numerically dominant microorganism is lysed and replaced by a virus-resistant one. Here, we isolated a microalga strain and its infective dsDNA virus whose dynamics are characterized instead by parallel growth of both the microalga and the virus. Experimental evolution of clonal lines revealed that this viral production originates from the lysis of a minority of virus-susceptible cells, which are regenerated from resistant cells. Whole-genome sequencing demonstrated that this resistant-susceptible switch involved a large deletion on one chromosome. Mathematical modeling explained how the switch maintains stable microalga-virus population dynamics consistent with their observed growth pattern. Comparative genomics confirmed an ancient origin of this "accordion" chromosome despite a lack of sequence conservation. Together, our results show how dynamic genomic rearrangements may account for a previously overlooked coexistence mechanism in microalgae-virus interactions.

Diversity and evolution of the emerging Pandoraviridae family

With DNA genomes reaching 2.5 Mb packed in particles of bacterium-like shape and dimension, the first two Acanthamoeba-infecting pandoraviruses remained up to now the most complex viruses since their discovery in 2013. Our isolation of three new strains from distant locations and environments is now used to perform the first comparative genomics analysis of the emerging worldwide-distributed Pandoraviridae family. Thorough annotation of the genomes combining transcriptomic, proteomic, and bioinformatic analyses reveals many non-coding transcripts and significantly reduces the former set of predicted protein-coding genes. Here we show that the pandoraviruses exhibit an open pan-genome, the enormous size of which is not adequately explained by gene duplications or horizontal transfers. As most of the strain-specific genes have no extant homolog and exhibit statistical features comparable to intergenic regions, we suggest that de novo gene creation could contribute to the evolution of the giant pandoravirus genomes.

Giant viruses are visible by light microscopy and have unusually long genomes. Here, the authors report three new members of the Pandoraviridae family and investigate their evolution and diversity.

Re-Assembly and Analysis of an Ancient Variola Virus Genome

We report a major improvement to the assembly of published short read sequencing data from an ancient variola virus (VARV) genome by the removal of contig-capping sequencing tags and manual searches for gap-spanning reads. The new assembly, together with camelpox and taterapox genomes, permitted new dates to be calculated for the last common ancestor of all VARV genomes. The analysis of recently sequenced VARV-like cowpox virus genomes showed that single nucleotide polymorphisms (SNPs) and amino acid changes in the vaccinia virus (VACV)-Cop-O1L ortholog, predicted to be associated with VARV host specificity and virulence, were introduced into the lineage before the divergence of these viruses. A comparison of the ancient and modern VARV genome sequences also revealed a measurable drift towards adenine + thymine (A + T) richness.

Genome sequence of Perigonia lusca single nucleopolyhedrovirus: insights into the evolution of a nucleotide metabolism enzyme in the family Baculoviridae

The genome of a novel group II alphabaculovirus, Perigonia lusca single nucleopolyhedrovirus (PeluSNPV), was sequenced and shown to contain 132,831 bp with 145 putative ORFs (open reading frames) of at least 50 amino acids. An interesting feature of this novel genome was the presence of a putative nucleotide metabolism enzyme-encoding gene (pelu112). The pelu112 gene was predicted to encode a fusion of thymidylate kinase (tmk) and dUTP diphosphatase (dut). Phylogenetic analysis indicated that baculoviruses have independently acquired tmk and dut several times during their evolution. Two homologs of the tmk-dut fusion gene were separately introduced into the Autographa californica multiple nucleopolyhedrovirus (AcMNPV) genome, which lacks tmk and dut. The recombinant baculoviruses produced viral DNA, virus progeny, and some viral proteins earlier during in vitro infection and the yields of viral occlusion bodies were increased 2.5-fold when compared to the parental virus. Interestingly, both enzymes appear to retain their active sites, based on separate modeling using previously solved crystal structures. We suggest that the retention of these tmk-dut fusion genes by certain baculoviruses could be related to accelerating virus replication and to protecting the virus genome from deleterious mutation.

A betabaculovirus encoding a gp64 homolog

Background

A betabaculovirus (DisaGV) was isolated from Diatraea saccharalis (Lepidoptera: Crambidae), one of the most important insect pests of the sugarcane and other monocot cultures in Brazil.

Results

The complete genome sequence of DisaGV was determined using the 454-pyrosequencing method. The genome was 98,392 bp long, which makes it the smallest lepidopteran-infecting baculovirus sequenced to date. It had a G + C content of 29.7 % encoding 125 putative open reading frames (ORF). All the 37 baculovirus core genes and a set of 19 betabaculovirus-specific genes were found. A group of 13 putative genes was not found in any other baculovirus genome sequenced so far. A phylogenetic analysis indicated that DisaGV is a member of Betabaculovirus genus and that it is a sister group to a cluster formed by ChocGV, ErelGV, PiraGV isolates, ClanGV, CaLGV, CpGV, CrleGV, AdorGV, PhopGV and EpapGV. Surprisingly, we found in the DisaGV genome a G protein-coupled receptor related to lepidopteran and other insect virus genes and a gp64 homolog, which is likely a product of horizontal gene transfer from Group 1 alphabaculoviruses.

Conclusion

DisaGV represents a distinct lineage of the genus Betabaculovirus. It is closely related to the CpGV-related group and presents the smallest genome in size so far. Remarkably, we found a homolog of gp64, which was reported solely in group 1 alphabaculovirus genomes so far.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2408-9) contains supplementary material, which is available to authorized users.

Homology-independent metrics for comparative genomics

A mainstream procedure to analyze the wealth of genomic data available nowadays is the detection of homologous regions shared across genomes, followed by the extraction of biological information from the patterns of conservation and variation observed in such regions. Although of pivotal importance, comparative genomic procedures that rely on homology inference are obviously not applicable if no homologous regions are detectable. This fact excludes a considerable portion of “genomic dark matter” with no significant similarity — and, consequently, no inferred homology to any other known sequence — from several downstream comparative genomic methods. In this review we compile several sequence metrics that do not rely on homology inference and can be used to compare nucleotide sequences and extract biologically meaningful information from them. These metrics comprise several compositional parameters calculated from sequence data alone, such as GC content, dinucleotide odds ratio, and several codon bias metrics. They also share other interesting properties, such as pervasiveness (patterns persist on smaller scales) and phylogenetic signal. We also cite examples where these homology-independent metrics have been successfully applied to support several bioinformatics challenges, such as taxonomic classification of biological sequences without homology inference. They where also used to detect higher-order patterns of interactions in biological systems, ranging from detecting coevolutionary trends between the genomes of viruses and their hosts to characterization of gene pools of entire microbial communities. We argue that, if correctly understood and applied, homology-independent metrics can add important layers of biological information in comparative genomic studies without prior homology inference.

Dip in the gene pool: metagenomic survey of natural coccolithovirus communities

Despite the global oceanic distribution and recognised biogeochemical impact of coccolithoviruses (EhV), their diversity remains poorly understood. Here we employed a metagenomic approach to study the occurrence and progression of natural EhV community genomic variability. Analysis of EhV metagenomes from the early and late stages of an induced bloom led to three main discoveries. First, we observed resilient and specific genomic signatures in the EhV community associated with the Norwegian coast, which reinforce the existence of limitations to the capacity of dispersal and genomic exchange among EhV populations. Second, we identified a hyper-variable region (approximately 21kbp long) in the coccolithovirus genome. Third, we observed a clear trend for EhV relative amino-acid diversity to reduce from early to late stages of the bloom. This study validated two new methodological combinations, and proved very useful in the discovery of new genomic features associated with coccolithovirus natural communities.

Horizontal gene transfer in plants

Horizontal gene transfer (HGT) describes the transmission of genetic material across species boundaries. HGT often occurs in microbic and eukaryotic genomes. However, the pathways by which HGTs occur in multicellular eukaryotes, especially in plants, are not well understood. We systematically summarized more than ten possible pathways for HGT. The intimate contact which frequently occurs in parasitism, symbiosis, pathogen, epiphyte, entophyte, and grafting interactions could promote HGTs between two species. Besides these direct transfer methods, genes can be exchanged with a vector as a bridge: possible vectors include pollen, fungi, bacteria, viruses, viroids, plasmids, transposons, and insects. HGT, especially when involving horizontal transfer of transposable elements, is recognized as a significant force propelling genomic variation and biological innovation, playing an important functional and evolutionary role in both eukaryotic and prokaryotic genomes. We proposed possible mechanisms by which HGTs can occur, which is useful in understanding the genetic information exchange among distant species or distant cellular components.

Rapid quantification of sequence repeats to resolve the size, structure and contents of bacterial genomes

BACKGROUND

The numerous classes of repeats often impede the assembly of genome sequences from the short reads provided by new sequencing technologies. We demonstrate a simple and rapid means to ascertain the repeat structure and total size of a bacterial or archaeal genome without the need for assembly by directly analyzing the abundances of distinct k-mers among reads.

RESULTS

The sensitivity of this procedure to resolve variation within a bacterial species is demonstrated: genome sizes and repeat structure of five environmental strains of E. coli from short Illumina reads were estimated by this method, and total genome sizes corresponded well with those obtained for the same strains by pulsed-field gel electrophoresis. In addition, this approach was applied to read-sets for completed genomes and shown to be accurate over a wide range of microbial genome sizes.

CONCLUSIONS

Application of these procedures, based solely on k-mer abundances in short read data sets, allows aspects of genome structure to be resolved that are not apparent from conventional short read assemblies. This knowledge of the repetitive content of genomes provides insights into genome evolution and diversity.

Open questions about giant viruses

The recent discovery of giant viruses exhibiting double-stranded DNA genomes larger than a million base pairs, encoding more than a thousand proteins and packed in near micron-sized icosahedral particles, opened a new and unexpected chapter in virology. As of today, these giant viruses and their closest relatives of lesser dimensions infect unicellular eukaryotes found in aquatic environments, but belonging to a wide diversity of early branching phyla. This broad phylogenetic distribution of hosts is consistent with the hypothesis that giant viruses originated prior to the radiation of the eukaryotic domain and/or might have been involved in the partition of nuclear versus cytoplasmic functions in ancestral cells. The distinctive features of the known giant viruses, in particular the recurrent presence of components of the translation apparatus in their proteome, raise a number of fundamental questions about their origin, their mode of evolution, and the relationship they may entertain with other dsDNA viruses, the genome size of which exhibits the widest distribution among all biological entities, from less than 5 kb to more than 1.25 Mb (a ratio of 1:250). At a more conceptual level, the convergence between the discovery of increasingly reduced parasitic cellular organisms and that of giant viruses exhibiting a widening array of cellular-like functions may ultimately abolish the historical discontinuity between the viral and the cellular world.

Exploring nucleo-cytoplasmic large DNA viruses in Tara Oceans microbial metagenomes

Nucleo-cytoplasmic large DNA viruses (NCLDVs) constitute a group of eukaryotic viruses that can have crucial ecological roles in the sea by accelerating the turnover of their unicellular hosts or by causing diseases in animals. To better characterize the diversity, abundance and biogeography of marine NCLDVs, we analyzed 17 metagenomes derived from microbial samples (0.2–1.6 μm size range) collected during the Tara Oceans Expedition. The sample set includes ecosystems under-represented in previous studies, such as the Arabian Sea oxygen minimum zone (OMZ) and Indian Ocean lagoons. By combining computationally derived relative abundance and direct prokaryote cell counts, the abundance of NCLDVs was found to be in the order of 10⁴–10⁵ genomes ml⁻¹ for the samples from the photic zone and 10²–10³ genomes ml⁻¹ for the OMZ. The Megaviridae and Phycodnaviridae dominated the NCLDV populations in the metagenomes, although most of the reads classified in these families showed large divergence from known viral genomes. Our taxon co-occurrence analysis revealed a potential association between viruses of the Megaviridae family and eukaryotes related to oomycetes. In support of this predicted association, we identified six cases of lateral gene transfer between Megaviridae and oomycetes. Our results suggest that marine NCLDVs probably outnumber eukaryotic organisms in the photic layer (per given water mass) and that metagenomic sequence analyses promise to shed new light on the biodiversity of marine viruses and their interactions with potential hosts.

Towards defining the chloroviruses: a genomic journey through a genus of large DNA viruses

Background

Giant viruses in the genus Chlorovirus (family Phycodnaviridae) infect eukaryotic green microalgae. The prototype member of the genus, Paramecium bursaria chlorella virus 1, was sequenced more than 15 years ago, and to date there are only 6 fully sequenced chloroviruses in public databases. Presented here are the draft genome sequences of 35 additional chloroviruses (287 – 348 Kb/319 – 381 predicted protein encoding genes) collected across the globe; they infect one of three different green algal species. These new data allowed us to analyze the genomic landscape of 41 chloroviruses, which revealed some remarkable features about these viruses.

Results

Genome colinearity, nucleotide conservation and phylogenetic affinity were limited to chloroviruses infecting the same host, confirming the validity of the three previously known subgenera. Clues for the existence of a fourth new subgenus indicate that the boundaries of chlorovirus diversity are not completely determined. Comparison of the chlorovirus phylogeny with that of the algal hosts indicates that chloroviruses have changed hosts in their evolutionary history. Reconstruction of the ancestral genome suggests that the last common chlorovirus ancestor had a slightly more diverse protein repertoire than modern chloroviruses. However, more than half of the defined chlorovirus gene families have a potential recent origin (after Chlorovirus divergence), among which a portion shows compositional evidence for horizontal gene transfer. Only a few of the putative acquired proteins had close homologs in databases raising the question of the true donor organism(s). Phylogenomic analysis identified only seven proteins whose genes were potentially exchanged between the algal host and the chloroviruses.

Conclusion

The present evaluation of the genomic evolution pattern suggests that chloroviruses differ from that described in the related Poxviridae and Mimiviridae. Our study shows that the fixation of algal host genes has been anecdotal in the evolutionary history of chloroviruses. We finally discuss the incongruence between compositional evidence of horizontal gene transfer and lack of close relative sequences in the databases, which suggests that the recently acquired genes originate from a still largely un-sequenced reservoir of genomes, possibly other unknown viruses that infect the same hosts.

Genomics of Megavirus and the elusive fourth domain of Life

We recently described Megavirus chilensis, a giant virus isolated off the coast of Chile, also replicating in fresh water acanthamoeba. Its 1,259,197-bp genome encodes 1,120 proteins and is the largest known viral genome. Megavirus and its closest relative Mimivirus only share 594 orthologous genes, themselves sharing only 50% of identical residues in average. Despite this divergence, comparable to the maximal divergence exhibited by bacteria within the same division (e.g., gamma proteobacteria), Megavirus retained all of the genomic features unique to Mimivirus, in particular its genes encoding key-elements of the translation apparatus, a trademark of cellular organisms. Besides homologs to the four aminoacyl-tRNA synthetases (aaRS) encoded by Mimivirus, Megavirus added three additional ones, raising the total of known virus-encoded aaRS to seven: IleRS, TrpRS, AsnRS, ArgRS, CysRS, MetRS, TyrRS. This finding strongly suggests that large DNA viruses derived from an ancestral cellular genome by reductive evolution. The nature of this cellular ancestor remains hotly debated.

Lateral gene exchanges shape the genomes of amoeba-resisting microorganisms

Based on Darwin's concept of the tree of life, vertical inheritance was thought to be dominant, and mutations, deletions, and duplication were streaming the genomes of living organisms. In the current genomic era, increasing data indicated that both vertical and lateral gene inheritance interact in space and time to trigger genome evolution, particularly among microorganisms sharing a given ecological niche. As a paradigm to their diversity and their survival in a variety of cell types, intracellular microorganisms, and notably intracellular bacteria, were considered as less prone to lateral genetic exchanges. Such specialized microorganisms generally have a smaller gene repertoire because they do rely on their host's factors for some basic regulatory and metabolic functions. Here we review events of lateral gene transfer (LGT) that illustrate the genetic exchanges among intra-amoebal microorganisms or between the microorganism and its amoebal host. We tentatively investigate the functions of laterally transferred genes in the light of the interaction with their host as they should confer a selective advantage and success to the amoeba-resisting microorganisms (ARMs).

Horizontal gene transfer confers fermentative metabolism in the respiratory-deficient plant trypanosomatid Phytomonas serpens

Among trypanosomatids, the genus Phytomonas is the only one specifically adapted to infect plants. These hosts provide a particular habitat with a plentiful supply of carbohydrates. Phytomonas sp. lacks a cytochrome-mediated respiratory chain and Krebs cycle, and ATP production relies predominantly on glycolysis. We have characterised the complete gene encoding a putative pyruvate/indolepyruvate decarboxylase (PDC/IPDC) (548 amino acids) of P. serpens, that displays high amino acid sequence similarity with phytobacteria and Leishmania enzymes. No orthologous PDC/IPDC genes were found in Trypanosoma cruzi or T. brucei. Conservation of the PDC/IPDC gene sequence was verified in 14 Phytomonas isolates. A phylogenetic analysis shows that Phytomonas protein is robustly monophyletic with Leishmania spp. and C. fasciculata enzymes. In the trees this clade appears as a sister group of indolepyruvate decarboxylases of γ-proteobacteria. This supports the proposition that a horizontal gene transfer event from a donor phytobacteria to a recipient ancestral trypanosome has occurred prior to the separation between Phytomonas, Leishmania and Crithidia. We have measured the PDC activity in P. serpens cell extracts. The enzyme has a Km value for pyruvate of 1.4mM. The acquisition of a PDC, a key enzyme in alcoholic fermentation, explains earlier observations that ethanol is one of the major end-products of glucose catabolism under aerobic and anaerobic conditions. This represents an alternative and necessary route to reoxidise part of the NADH produced in the highly demanding glycolytic pathway and highlights the importance of this type of event in metabolic adaptation.

Phosphate transporters in marine phytoplankton and their viruses: cross-domain commonalities in viral-host gene exchanges

Phosphate (PO(4)) is an important limiting nutrient in marine environments. Marine cyanobacteria scavenge PO(4) using the high-affinity periplasmic phosphate binding protein PstS. The pstS gene has recently been identified in genomes of cyanobacterial viruses as well. Here, we analyse genes encoding transporters in genomes from viruses that infect eukaryotic phytoplankton. We identified inorganic PO(4) transporter-encoding genes from the PHO4 superfamily in several virus genomes, along with other transporter-encoding genes. Homologues of the viral pho4 genes were also identified in genome sequences from the genera that these viruses infect. Genome sequences were available from host genera of all the phytoplankton viruses analysed except the host genus Bathycoccus. Pho4 was recovered from Bathycoccus by sequencing a targeted metagenome from an uncultured Atlantic Ocean population. Phylogenetic reconstruction showed that pho4 genes from pelagophytes, haptophytes and infecting viruses were more closely related to homologues in prasinophytes than to those in what, at the species level, are considered to be closer relatives (e.g. diatoms). We also identified PHO4 superfamily members in ocean metagenomes, including new metagenomes from the Pacific Ocean. The environmental sequences grouped with pelagophytes, haptophytes, prasinophytes and viruses as well as bacteria. The analyses suggest that multiple independent pho4 gene transfer events have occurred between marine viruses and both eukaryotic and bacterial hosts. Additionally, pho4 genes were identified in available genomes from viruses that infect marine eukaryotes but not those that infect terrestrial hosts. Commonalities in marine host-virus gene exchanges indicate that manipulation of host-PO(4) uptake is an important adaptation for viral proliferation in marine systems. Our findings suggest that PO(4) -availability may not serve as a simple bottom-up control of marine phytoplankton.

An integrative method for identifying the over-annotated protein-coding genes in microbial genomes

The falsely annotated protein-coding genes have been deemed one of the major causes accounting for the annotating errors in public databases. Although many filtering approaches have been designed for the over-annotated protein-coding genes, some are questionable due to the resultant increase in false negative. Furthermore, there is no webserver or software specifically devised for the problem of over-annotation. In this study, we propose an integrative algorithm for detecting the over-annotated protein-coding genes in microorganisms. Overall, an average accuracy of 99.94% is achieved over 61 microbial genomes. The extremely high accuracy indicates that the presented algorithm is efficient to differentiate the protein-coding genes from the non-coding open reading frames. Abundant analyses show that the predicting results are reliable and the integrative algorithm is robust and convenient. Our analysis also indicates that the over-annotated protein-coding genes can cause the false positive of horizontal gene transfers detection. The webserver of the proposed algorithm can be freely accessible from www.cbi.seu.edu.cn/RPGM.

The distribution of a phage-related insertion sequence element in the cyanobacterium, Microcystis aeruginosa

The cyanophage Ma-LMM01, specifically-infecting Microcystis aeruginosa, has an insertion sequence (IS) element that we named IS607-cp showing high nucleotide similarity to a counterpart in the genome of the cyanobacterium Cyanothece sp. We tested 21 strains of M. aeruginosa for the presence of IS607-cp using PCR and detected the element in strains NIES90, NIES112, NIES604, and RM6. Thermal asymmetric interlaced PCR (TAIL-PCR) revealed each of these strains has multiple copies of IS607-cp. Some of the ISs were classified into three types based on their inserted positions; IS607-cp-1 is common in strains NIES90, NIES112 and NIES604, whereas IS607-cp-2 and IS607-cp-3 are specific to strains NIES90 and RM6, respectively. This multiplicity may reflect the replicative transposition of IS607-cp. The sequence of IS607-cp in Ma-LMM01 showed robust affinity to those found in M. aeruginosa and Cyanothece spp. in a phylogenetic tree inferred from counterparts of various bacteria. This suggests the transfer of IS607-cp between the cyanobacterium and its cyanophage. We discuss the potential role of Ma-LMM01-related phages as donors of IS elements that may mediate the transfer of IS607-cp; and thereby partially contribute to the genome plasticity of M. aeruginosa.

Two new subfamilies of DNA mismatch repair proteins (MutS) specifically abundant in the marine environment

MutS proteins are ubiquitous in cellular organisms and have important roles in DNA mismatch repair or recombination. In the virus world, the amoeba-infecting Mimivirus, as well as the recently sequenced Cafeteria roenbergensis virus are known to encode a MutS related to the homologs found in octocorals and ɛ-proteobacteria. To explore the presence of MutS proteins in other viral genomes, we performed a genomic survey of four giant viruses ('giruses') (Pyramimonas orientalis virus (PoV), Phaeocystis pouchetii virus (PpV), Chrysochromulina ericina virus (CeV) and Heterocapsa circularisquama DNA virus (HcDNAV)) that infect unicellular marine algae. Our analysis revealed the presence of a close homolog of Mimivirus MutS in all the analyzed giruses. These viral homologs possess a specific domain structure, including a C-terminal HNH-endonuclease domain, defining the new MutS7 subfamily. We confirmed the presence of conserved mismatch recognition residues in all members of the MutS7 subfamily, suggesting their role in DNA mismatch repair rather than DNA recombination. PoV and PpV were found to contain an additional type of MutS, which we propose to call MutS8. The MutS8 proteins in PoV and PpV were found to be closely related to homologs from 'Candidatus Amoebophilus asiaticus', an obligate intracellular amoeba-symbiont belonging to the Bacteroidetes. Furthermore, our analysis revealed that MutS7 and MutS8 are abundant in marine microbial metagenomes and that a vast majority of these environmental sequences are likely of girus origin. Giruses thus seem to represent a major source of the underexplored diversity of the MutS family in the microbial world.

Mimivirus: the emerging paradox of quasi-autonomous viruses

What is a virus? Are viruses alive? Should they be classified among microorganisms? One would expect these simple questions to have been settled a century after the discovery of the first viral disease. For years, modern virology successfully unravelled the huge diversity of viruses in terms of genetic material, replication mechanism, pathogenicity, host infection, and more recently particle structure, planet-wide distribution and ecological significance. Yet, little progress was made in understanding their evolutionary origin(s), as well as the fundamental nature of their relationship with the cellular world. Thanks to the recent studies on Mimivirus and other large DNA viruses, we are now entering a new era where the most basic concepts about viruses are revisited, including their true nature, how fundamentally different they are from cellular microorganisms, and how essential they might have been in the major innovations that punctuated the evolution of life.

Identification of horizontal gene transfer and recombination of PsaA gene in streptococcus mitis group

Pneumococcal surface adhesin A (psaA) gene is universally confirmed as one of the Streptococcus pneumoniae adhesion genes, but it is disputed whether the psaA gene is a Streptococcus pneumoniae species-specific gene. In the present study, the presence of the psaA gene in 34 streptococcus mitis group isolates was identified by the PCR approach and a comparison of sequencing PCR products (Streptococcus pneumoniae R6 as the control strain). Also, the evolutionary scenarios of these psaA genes in these streptococcus mitis group isolates were analyzed by a phylogenetic tree based on the housekeeping genes (sodA and rnpB) and psaA genes. As a result, a high degree of conservation of open reading frame sequences in all six Streptococcus pneumoniae strains (100% similarity) and in the other species of the streptococcus mitis group (92.6-100% similarity) was revealed. Further genetics research based on housekeeping genes and psaA gene phylogenies showed that the psaA gene was of vertical inheritance only in Streptococcus pneumoniae; however, high-frequency horizontal psaA gene transfer and recombination occurred in the other species of the streptococcus mitis group. These findings confirmed that the psaA gene was not a Streptococcus pneumoniae species-specific gene, and high-frequency HGT and recombination events may explain the presence of the psaA gene in the other species of the streptococcus mitis group.

The source of self: genetic parasites and the origin of adaptive immunity

Stable colonization of the host by viruses (genetic parasites) can alter the systems of host identity and provide immunity against related viruses. To attain the needed stability, some viruses of prokaryotes (P1 phage) use a strategy called an addiction module. The linked protective and destructive gene functions of an addiction module insures both virus persistence but will also destroy cells that interrupt this module and thereby prevent infection by competitors. Previously, I have generalized this concept to also include persistent and lytic states of virus infection, which can be considered as a virus addiction module. Such states often involve defective viruses. In this report, I examine the origin of the adaptive immune system from the perspective of a virus addiction module. The likely role of both endogenous and exogenous retroviruses, DNA viruses, and their defective elements is considered in the origin of all the basal components of adaptive immunity (T-cell receptor, RAG-mediated gene rearrangement, clonal lymphocyte proliferation, antigen surface presentation, apoptosis, and education of immune cells). It is concluded that colonization by viruses and their defectives provides a more coherent explanation for the origin of adaptive immunity.

Mimivirus and its virophage

Mimivirus, a virus infecting amoebae of the acanthamoeba genus, is the prototype member of the Mimiviridae, the latest addition to the family of the nucleocytoplasmic large DNA viruses, already including the Poxviridae, the Iridoviridae, the Asfarviridae, and the Phycodnaviridae. Because of the size of its particle-a fiber-covered icosahedral protein capsid 0.75 microm in diameter-Mimivirus was initially mistaken for a parasitic bacterium. Its 1.2-Mb genome sequence encodes more than 900 proteins, many of them associated with functions never before encountered in a virus, such as four aminoacyl-tRNA synthetases. These findings revived the debate about the origin of DNA viruses and their possible role in the emergence of the eukaryotic nucleus. The recent isolation of a new type of satellite virus, called a virophage, associated with a second strain of Mimivirus, confirmed its unique position within the virus world. Post-genomic studies are now in progress, slowly shedding some light on the physiology of the most complex virus isolated to date.

Viruses are essential agents within the roots and stem of the tree of life

In contrast with former definitions of life limited to membrane-bound cellular life forms which feed, grow, metabolise and replicate (i) a role of viruses as genetic symbionts, (ii) along with peripheral phenomena such as cryptobiosis and (iii) the horizontal nature of genetic information acquisition and processing broaden our view of the tree of life. Some researchers insist on the traditional textbook conviction of what is part of the community of life. In a recent review [Moreira, D., Lopez-Garcia, P., 2009. Ten reasons to exclude viruses from the tree of life. Nat. Rev. Microbiol. 7, 306-311.] they assemble four main arguments which should exclude viruses from the tree of life because of their inability to self-sustain and self-replicate, their polyphyly, the cellular origin of their cell-like genes and the volatility of their genomes. In this article we will show that these features are not coherent with current knowledge about viruses but that viral agents play key roles within the roots and stem of the tree of life.

Viruses manipulate the marine environment

Marine viruses affect Bacteria, Archaea and eukaryotic organisms and are major components of the marine food web. Most studies have focused on their role as predators and parasites, but many of the interactions between marine viruses and their hosts are much more complicated. A series of recent studies has shown that viruses have the ability to manipulate the life histories and evolution of their hosts in remarkable ways, challenging our understanding of this almost invisible world.

Whole-proteome phylogeny of large dsDNA virus families by an alignment-free method

The vast sequence divergence among different virus groups has presented a great challenge to alignment-based sequence comparison among different virus families. Using an alignment-free comparison method, we construct the whole-proteome phylogeny for a population of viruses from 11 viral families comprising 142 large dsDNA eukaryote viruses. The method is based on the feature frequency profiles (FFP), where the length of the feature (l-mer) is selected to be optimal for phylogenomic inference. We observe that (i) the FFP phylogeny segregates the population into clades, the membership of each has remarkable agreement with current classification by the International Committee on the Taxonomy of Viruses, with one exception that the mimivirus joins the phycodnavirus family; (ii) the FFP tree detects potential evolutionary relationships among some viral families; (iii) the relative position of the 3 herpesvirus subfamilies in the FFP tree differs from gene alignment-based analysis; (iv) the FFP tree suggests the taxonomic positions of certain "unclassified" viruses; and (v) the FFP method identifies candidates for horizontal gene transfer between virus families.

A hidden reservoir of integrative elements is the major source of recently acquired foreign genes and ORFans in archaeal and bacterial genomes

A large-scale survey of potential recently acquired integrative elements in 119 archaeal and bacterial genomes reveals that many recently acquired genes have originated from integrative elements

Background

Archaeal and bacterial genomes contain a number of genes of foreign origin that arose from recent horizontal gene transfer, but the role of integrative elements (IEs), such as viruses, plasmids, and transposable elements, in this process has not been extensively quantified. Moreover, it is not known whether IEs play an important role in the origin of ORFans (open reading frames without matches in current sequence databases), whose proportion remains stable despite the growing number of complete sequenced genomes.

Results

We have performed a large-scale survey of potential recently acquired IEs in 119 archaeal and bacterial genomes. We developed an accurate in silico Markov model-based strategy to identify clusters of genes that show atypical sequence composition (clusters of atypical genes or CAGs) and are thus likely to be recently integrated foreign elements, including IEs. Our method identified a high number of new CAGs. Probabilistic analysis of gene content indicates that 56% of these new CAGs are likely IEs, whereas only 7% likely originated via horizontal gene transfer from distant cellular sources. Thirty-four percent of CAGs remain unassigned, what may reflect a still poor sampling of IEs associated with bacterial and archaeal diversity. Moreover, our study contributes to the issue of the origin of ORFans, because 39% of these are found inside CAGs, many of which likely represent recently acquired IEs.

Conclusions

Our results strongly indicate that archaeal and bacterial genomes contain an impressive proportion of recently acquired foreign genes (including ORFans) coming from a still largely unexplored reservoir of IEs.

The evolution of guanylyl cyclases as multidomain proteins: conserved features of kinase-cyclase domain fusions

Guanylyl cyclases (GCs) are enzymes that generate cyclic GMP and regulate different physiologic and developmental processes in a number of organisms. GCs possess sequence similarity to class III adenylyl cyclases (ACs) and are present as either membrane-bound receptor GCs or cytosolic soluble GCs. We sought to determine the evolution of GCs using a large-scale bioinformatic analysis and found multiple lineage-specific expansions of GC genes in the genomes of many eukaryotes. Moreover, a few GC-like proteins were identified in prokaryotes, which come fused to a number of different domains, suggesting allosteric regulation of nucleotide cyclase activity. Eukaryotic receptor GCs are associated with a kinase homology domain (KHD), and phylogenetic analysis of these proteins suggest coevolution of the KHD and the associated cyclase domain as well as a conservation of the sequence and the size of the linker region between the KHD and the associated cyclase domain. Finally, we also report the existence of mimiviral proteins that contain putative active kinase domains associated with a cyclase domain, which could suggest early evolution of the fusion of these two important domains involved in signal transduction.

Poxvirus protein evolution: family wide assessment of possible horizontal gene transfer events

To investigate the evolutionary origins of proteins encoded by the Poxviridae family of viruses, we examined all poxvirus protein coding genes using a method of characterizing and visualizing the similarity between these proteins and taxonomic subsets of proteins in GenBank. Our analysis divides poxvirus proteins into categories based on their relative degree of similarity to two different taxonomic subsets of proteins such as all eukaryote vs. all virus (except poxvirus) proteins. As an example, this allows us to identify, based on high similarity to only eukaryote proteins, poxvirus proteins that may have been obtained by horizontal transfer from their hosts. Although this method alone does not definitively prove horizontal gene transfer, it allows us to provide an assessment of the possibility of horizontal gene transfer for every poxvirus protein. Potential candidates can then be individually studied in more detail during subsequent investigation.

Results of our analysis demonstrate that in general, proteins encoded by members of the subfamily Chordopoxvirinae exhibit greater similarity to eukaryote proteins than to proteins of other virus families. In addition, our results reiterate the important role played by host gene capture in poxvirus evolution; highlight the functions of many genes poxviruses share with their hosts; and illustrate which host-like genes are present uniquely in poxviruses and which are also present in other virus families.

Bioinformatic detection of horizontally transferred DNA in bacterial genomes

We highlight a selection of recent research on computational methods and associated challenges surrounding the prediction of bacterial horizontal gene transfer. This research area continues to face controversy, but is becoming more critical as the importance of horizontal gene transfer in medically and ecologically important prokaryotic evolution is further appreciated.

Complete genome viral phylogenies suggests the concerted evolution of regulatory cores and accessory satellites

We consider the concerted evolution of viral genomes in four families of DNA viruses. Given the high rate of horizontal gene transfer among viruses and their hosts, it is an open question as to how representative particular genes are of the evolutionary history of the complete genome. To address the concerted evolution of viral genes, we compared genomic evolution across four distinct, extant viral families. For all four viral families we constructed DNA-dependent DNA polymerase-based (DdDp) phylogenies and in addition, whole genome sequence, as quantitative descriptions of inter-genome relationships. We found that the history of the polymerase gene was highly predictive of the history of the genome as a whole, which we explain in terms of repeated, co-divergence events of the core DdDp gene accompanied by a number of satellite, accessory genetic loci. We also found that the rate of gene gain in baculovirus and poxviruses proceeds significantly more quickly than the rate of gene loss and that there is convergent acquisition of satellite functions promoting contextual adaptation when distinct viral families infect related hosts. The congruence of the genome and polymerase trees suggests that a large set of viral genes, including polymerase, derive from a phylogenetically conserved core of genes of host origin, secondarily reinforced by gene acquisition from common hosts or co-infecting viruses within the host. A single viral genome can be thought of as a mutualistic network, with the core genes acting as an effective host and the satellite genes as effective symbionts. Larger virus genomes show a greater departure from linkage equilibrium between core and satellites functions.

Marine mimivirus relatives are probably large algal viruses

Background

Acanthamoeba polyphaga mimivirus is the largest known ds-DNA virus and its 1.2 Mb-genome sequence has revealed many unique features. Mimivirus occupies an independent lineage among eukaryotic viruses and its known hosts include only species from the Acanthamoeba genus. The existence of mimivirus relatives was first suggested by the analysis of the Sargasso Sea metagenomic data.

Results

We now further demonstrate the presence of numerous "mimivirus-like" sequences using a larger marine metagenomic data set. We also show that the DNA polymerase sequences from three algal viruses (CeV01, PpV01, PoV01) infecting different marine algal species (Chrysochromulina ericina, Phaeocystis pouchetii, Pyramimonas orientalis) are very closely related to their homolog in mimivirus.

Conclusion

Our results suggest that the numerous mimivirus-related sequences identified in marine environments are likely to originate from diverse large DNA viruses infecting phytoplankton. Micro-algae thus constitute a new category of potential hosts in which to look for new species of Mimiviridae.