Jaime Iranzo

The Double-Stranded DNA Virosphere as a Modular Hierarchical Network of Gene Sharing

ABSTRACT Virus genomes are prone to extensive gene loss, gain, and exchange and share no universal genes. Therefore, in a broad-scale study of virus evolution, gene and genome network analyses can complement traditional phylogenetics. We performed an exhaustive comparative analysis of the genomes of double-stranded DNA (dsDNA) viruses by using the bipartite network approach and found a robust hierarchical modularity in the dsDNA virosphere. Bipartite networks consist of two classes of nodes, with nodes in one class, in this case genomes, being connected via nodes of the second class, in this case genes. Such a network can be partitioned into modules that combine nodes from both classes. The bipartite network of dsDNA viruses includes 19 modules that form 5 major and 3 minor supermodules. Of these modules, 11 include tailed bacteriophages, reflecting the diversity of this largest group of viruses. The module analysis quantitatively validates and refines previously proposed nontrivial evolutionary relationships. An expansive supermodule combines the large and giant viruses of the putative order “Megavirales” with diverse moderate-sized viruses and related mobile elements. All viruses in this supermodule share a distinct morphogenetic tool kit with a double jelly roll major capsid protein. Herpesviruses and tailed bacteriophages comprise another supermodule, held together by a distinct set of morphogenetic proteins centered on the HK97-like major capsid protein. Together, these two supermodules cover the great majority of currently known dsDNA viruses. We formally identify a set of 14 viral hallmark genes that comprise the hubs of the network and account for most of the intermodule connections. IMPORTANCE Viruses and related mobile genetic elements are the dominant biological entities on earth, but their evolution is not sufficiently understood and their classification is not adequately developed. The key reason is the characteristic high rate of virus evolution that involves not only sequence change but also extensive gene loss, gain, and exchange. Therefore, in the study of virus evolution on a large scale, traditional phylogenetic approaches have limited applicability and have to be complemented by gene and genome network analyses. We applied state-of-the art methods of such analysis to reveal robust hierarchical modularity in the genomes of double-stranded DNA viruses. Some of the identified modules combine highly diverse viruses infecting bacteria, archaea, and eukaryotes, in support of previous hypotheses on direct evolutionary relationships between viruses from the three domains of cellular life. We formally identify a set of 14 viral hallmark genes that hold together the genomic network. Viruses and related mobile genetic elements are the dominant biological entities on earth, but their evolution is not sufficiently understood and their classification is not adequately developed. The key reason is the characteristic high rate of virus evolution that involves not only sequence change but also extensive gene loss, gain, and exchange. Therefore, in the study of virus evolution on a large scale, traditional phylogenetic approaches have limited applicability and have to be complemented by gene and genome network analyses. We applied state-of-the art methods of such analysis to reveal robust hierarchical modularity in the genomes of double-stranded DNA viruses. Some of the identified modules combine highly diverse viruses infecting bacteria, archaea, and eukaryotes, in support of previous hypotheses on direct evolutionary relationships between viruses from the three domains of cellular life. We formally identify a set of 14 viral hallmark genes that hold together the genomic network.

The impact of quasispecies dynamics on the use of therapeutics

The application of quasispecies theory to viral populations has boosted our understanding of how endogenous and exogenous features condition their adaptation. Mounting empirical evidence demonstrates that internal interactions within mutant spectra may cause unexpected responses to antiviral treatments. In this scenario, increased mutagenesis could be efficient at low mutagen doses due to the lethal action of defective genomes, whereas sequential administration of antiviral drugs might be superior to combination therapies. Our ability to predict the outcome of a particular therapy takes advantage of the complementary use of in vivo observations, in vitro experiments, and mathematical models.

Tempo and mode of inhibitor–mutagen antiviral therapies: A multidisciplinary approach

The continuous emergence of drug-resistant viruses is a major obstacle for the successful treatment of viral infections, thus representing a persistent spur to the search for new therapeutic strategies. Among them, multidrug treatments are currently at the forefront of pharmaceutical, clinical, and computational investigation. Still, there are many unknowns in the way that different drugs interact among themselves and with the pathogen that they aim to control. Inspired by experimental studies with picornavirus, here, we discuss the performance of sequential vs. combination therapies involving two dissimilar drugs: the mutagen ribavirin and an inhibitor of viral replication, guanidine. Because a systematic analysis of viral response to drug doses demands a precious amount of time and resources, we present and analyze an in silico model describing the dynamics of the viral population under the action of the two drugs. The model predicts the response of the viral population to any dose combination, the optimal therapy to be used in each case, and the way to minimize the probability of appearance of resistant mutants. In agreement with the theoretical predictions, in vitro experiments with foot-and-mouth disease virus confirm that the suitability of simultaneous or sequential administration depends on the drug doses. In addition, intrinsic replicative characteristics of the virus (e.g., replication through RNA only or a DNA intermediate) play a key role to determine the appropriateness of a sequential or combination therapy. Knowledge of several model parameters can be derived by means of few, simple experiments, such that the model and its predictions can be extended to other viral systems.

Evolutionary dynamics of genome segmentation in multipartite viruses

Multipartite viruses are formed by a variable number of genomic fragments packed in independent viral capsids. This fact poses stringent conditions on their transmission mode, demanding, in particular, a high multiplicity of infection (MOI) for successful propagation. The actual advantages of the multipartite viral strategy are as yet unclear. The origin of multipartite viruses represents an evolutionary puzzle. While classical theories suggested that a faster replication rate or higher replication fidelity would favour shorter segments, recent experimental results seem to point to an increased stability of virions with incomplete genomes as a factor able to compensate for the disadvantage of mandatory complementation. Using as main parameters differential stability as a function of genome length and MOI, we calculate the conditions under which a set of complementary segments of a viral genome would outcompete the non-segmented variant. Further, we examine the likeliness that multipartite viral forms could be the evolutionary outcome of the competition among the defective genomes of different lengths that spontaneously arise under replication of a complete, wild-type genome. We conclude that only multipartite viruses with a small number of segments could be produced in our scenario, and discuss alternative hypotheses for the origin of multipartite viruses with more than four segments.

Bipartite network analysis of the archaeal virosphere: evolutionary connections between viruses and capsid-less mobile elements

ABSTRACT Archaea and particularly hyperthermophilic crenarchaea are hosts to many unusual viruses with diverse virion shapes and distinct gene compositions. As is typical of viruses in general, there are no universal genes in the archaeal virosphere. Therefore, to obtain a comprehensive picture of the evolutionary relationships between viruses, network analysis methods are more productive than traditional phylogenetic approaches. Here we present a comprehensive comparative analysis of genomes and proteomes from all currently known taxonomically classified and unclassified, cultivated and uncultivated archaeal viruses. We constructed a bipartite network of archaeal viruses that includes two classes of nodes, the genomes and gene families that connect them. Dissection of this network using formal community detection methods reveals strong modularity, with 10 distinct modules and 3 putative supermodules. However, compared to similar previously analyzed networks of eukaryotic and bacterial viruses, the archaeal virus network is sparsely connected. With the exception of the tailed viruses related to bacteriophages of the order Caudovirales and the families Turriviridae and Sphaerolipoviridae that are linked to a distinct supermodule of eukaryotic and bacterial viruses, there are few connector genes shared by different archaeal virus modules. In contrast, most of these modules include, in addition to viruses, capsidless mobile elements, emphasizing tight evolutionary connections between the two types of entities in archaea. The relative contributions of distinct evolutionary origins, in particular from nonviral elements, and insufficient sampling to the sparsity of the archaeal virus network remain to be determined by further exploration of the archaeal virosphere. IMPORTANCE Viruses infecting archaea are among the most mysterious denizens of the virosphere. Many of these viruses display no genetic or even morphological relationship to viruses of bacteria and eukaryotes, raising questions regarding their origins and position in the global virosphere. Analysis of 5,740 protein sequences from 116 genomes allowed dissection of the archaeal virus network and showed that most groups of archaeal viruses are evolutionarily connected to capsidless mobile genetic elements, including various plasmids and transposons. This finding could reflect actual independent origins of the distinct groups of archaeal viruses from different nonviral elements, providing important insights into the emergence and evolution of the archaeal virome.

Viruses of archaea: Structural, functional, environmental and evolutionary genomics

Viruses of archaea represent one of the most enigmatic parts of the virosphere. Most of the characterized archaeal viruses infect extremophilic hosts and display remarkable diversity of virion morphotypes, many of which have never been observed among viruses of bacteria or eukaryotes. The uniqueness of the virion morphologies is matched by the distinctiveness of the genomes of these viruses, with ∼75% of genes encoding unique proteins, refractory to functional annotation based on sequence analyses. In this review, we summarize the state-of-the-art knowledge on various aspects of archaeal virus genomics. First, we outline how structural and functional genomics efforts provided valuable insights into the functions of viral proteins and revealed intricate details of the archaeal virus-host interactions. We then highlight recent metagenomics studies, which provided a glimpse at the diversity of uncultivated viruses associated with the ubiquitous archaea in the oceans, including Thaumarchaeota, Marine Group II Euryarchaeota, and others. These findings, combined with the recent discovery that archaeal viruses mediate a rapid turnover of thaumarchaea in the deep sea ecosystems, illuminate the prominent role of these viruses in the biosphere. Finally, we discuss the origins and evolution of archaeal viruses and emphasize the evolutionary relationships between viruses and non-viral mobile genetic elements. Further exploration of the archaeal virus diversity as well as functional studies on diverse virus-host systems are bound to uncover novel, unexpected facets of the archaeal virome.

Large-scale genomic analysis suggests a neutral punctuated dynamics of transposable elements in bacterial genomes.

Insertion sequences (IS) are the simplest and most abundant form of transposable DNA found in bacterial genomes. When present in multiple copies, it is thought that they can promote genomic plasticity and genetic exchange, thus being a major force of evolutionary change. The main processes that determine IS content in genomes are, though, a matter of debate. In this work, we take advantage of the large amount of genomic data currently available and study the abundance distributions of 33 IS families in 1811 bacterial chromosomes. This allows us to test simple models of IS dynamics and estimate their key parameters by means of a maximum likelihood approach. We evaluate the roles played by duplication, lateral gene transfer, deletion and purifying selection. We find that the observed IS abundances are compatible with a neutral scenario where IS proliferation is controlled by deletions instead of purifying selection. Even if there may be some cases driven by selection, neutral behavior dominates over large evolutionary scales. According to this view, IS and hosts tend to coexist in a dynamic equilibrium state for most of the time. Our approach also allows for a detection of recent IS expansions, and supports the hypothesis that rapid expansions constitute transient events—punctuations—during which the state of coexistence of IS and host becomes perturbated.

Inevitability of Genetic Parasites

Abstract Almost all cellular life forms are hosts to diverse genetic parasites with various levels of autonomy including plasmids, transposons and viruses. Theoretical modeling of the evolution of primordial replicators indicates that parasites (cheaters) necessarily evolve in such systems and can be kept at bay primarily via compartmentalization. Given the (near) ubiquity, abundance and diversity of genetic parasites, the question becomes pertinent: are such parasites intrinsic to life? At least in prokaryotes, the persistence of parasites is linked to the rate of horizontal gene transfer (HGT). We mathematically derive the threshold value of the minimal transfer rate required for selfish element persistence, depending on the element duplication and loss rates as well as the cost to the host. Estimation of the characteristic gene duplication, loss and transfer rates for transposons, plasmids and virus-related elements in multiple groups of diverse bacteria and archaea indicates that most of these rates are compatible with the long term persistence of parasites. Notably, a small but non-zero rate of HGT is also required for the persistence of non-parasitic genes. We hypothesize that cells cannot tune their horizontal transfer rates to be below the threshold required for parasite persistence without experiencing highly detrimental side-effects. As a lower boundary to the minimum DNA transfer rate that a cell can withstand, we consider the process of genome degradation and mutational meltdown of populations through Muller’s ratchet. A numerical assessment of this hypothesis suggests that microbial populations cannot purge parasites while escaping Muller’s ratchet. Thus, genetic parasites appear to be virtually inevitable in cellular organisms.

Myosin-driven transport network in plants

Significance A key feature of plant cells is rapid cytoplasmic streaming that is driven by myosin motors. However, specific mechanisms of myosin-dependent streaming are poorly understood. Here, we characterize a dense network of plant myosins and their receptors and adaptors that, jointly with myosins, appear to mediate cytoplasmic streaming through distinct endomembrane compartments. We additionally present data suggestive of a myosin-dependent nucleocytoplasmic trafficking pathway. The myosin network is an ancient functional module that was already present in the common ancestor of green algae and land plants but underwent a major expansion in the latter, probably contributing to land colonization by plants. We investigate the myosin XI-driven transport network in Arabidopsis using protein–protein interaction, subcellular localization, gene knockout, and bioinformatics analyses. The two major groups of nodes in this network are myosins XI and their membrane-anchored receptors (MyoB) that, together, drive endomembrane trafficking and cytoplasmic streaming in the plant cells. The network shows high node connectivity and is dominated by generalists, with a smaller fraction of more specialized myosins and receptors. We show that interaction with myosins and association with motile vesicles are common properties of the MyoB family receptors. We identify previously uncharacterized myosin-binding proteins, putative myosin adaptors that belong to two unrelated families, with four members each (MadA and MadB). Surprisingly, MadA1 localizes to the nucleus and is rapidly transported to the cytoplasm, suggesting the existence of myosin XI-driven nucleocytoplasmic trafficking. In contrast, MadA2 and MadA3, as well as MadB1, partition between the cytosolic pools of motile endomembrane vesicles that colocalize with myosin XI-K and diffuse material that does not. Gene knockout analysis shows that MadB1–4 contribute to polarized root hair growth, phenocopying myosins, whereas MadA1–4 are redundant for this process. Phylogenetic analysis reveals congruent evolutionary histories of the myosin XI, MyoB, MadA, and MadB families. All these gene families emerged in green algae and show concurrent expansions via serial duplication in flowering plants. Thus, the myosin XI transport network increased in complexity and robustness concomitantly with the land colonization by flowering plants and, by inference, could have been a major contributor to this process.

Immunity, suicide or both? Ecological determinants for the combined evolution of anti-pathogen defense systems

BackgroundParasite-host arms race is one of the key factors in the evolution of life. Most cellular life forms, in particular prokaryotes, possess diverse forms of defense against pathogens including innate immunity, adaptive immunity and programmed cell death (altruistic suicide). Coevolution of these different but interacting defense strategies yields complex evolutionary regimes.ResultsWe develop and extensively analyze a computational model of coevolution of different defense strategies to show that suicide as a defense mechanism can evolve only in structured populations and when the attainable degree of immunity against pathogens is limited. The general principle of defense evolution seems to be that hosts do not evolve two costly defense mechanisms when one is sufficient. Thus, the evolutionary interplay of innate immunity, adaptive immunity and suicide, leads to an equilibrium state where the combination of all three defense strategies is limited to a distinct, small region of the parameter space. The three strategies can stably coexist only if none of them are highly effective. Coupled adaptive immunity-suicide systems, the existence of which is implied by the colocalization of genes for the two types of defense in prokaryotic genomes, can evolve either when immunity-associated suicide is more efficacious than other suicide systems or when adaptive immunity functionally depends on the associated suicide system.ConclusionsComputational modeling reveals a broad range of outcomes of coevolution of anti-pathogen defense strategies depending on the relative efficacy of different mechanisms and population structure. Some of the predictions of the model appear compatible with recent experimental evolution results and call for additional experiments.