Genic incompatibilities in two hybrid bacteriophages

Interspecies Recombination Has Driven the Macroevolution of Cassava Mosaic Begomoviruses

Begomoviruses (family Geminiviridae, genus Begomovirus) significantly hamper crop production and threaten food security around the world. The frequent emergence of new begomovirus genotypes is facilitated by high mutation frequencies and the propensity to recombine and reassort. Homologous recombination has been especially implicated in the emergence of novel cassava mosaic begomovirus (CMB) genotypes, which cause cassava mosaic disease (CMD). Cassava (Manihot esculenta) is a staple food crop throughout Africa and an important industrial crop in Asia, two continents where production is severely constrained by CMD. The CMD species complex is comprised of 11 bipartite begomovirus species with ample distribution throughout Africa and the Indian subcontinent. While recombination is regarded as a frequent occurrence for CMBs, a revised, systematic assessment of recombination and its impact on CMB phylogeny is currently lacking. We assembled data sets of all publicly available, full-length DNA-A (n = 880) and DNA-B (n = 369) nucleotide sequences from the 11 recognized CMB species. Phylogenetic networks and complementary recombination detection methods revealed extensive recombination among the CMB sequences. Six out of the 11 species descended from unique interspecies recombination events. Estimates of recombination and mutation rates revealed that all species experience mutation more frequently than recombination, but measures of population divergence indicate that recombination is largely responsible for the genetic differences between species. Our results support that recombination has significantly impacted the CMB phylogeny and has driven speciation in the CMD species complex. IMPORTANCE Cassava mosaic disease (CMD) is a significant threat to cassava production throughout Africa and Asia. CMD is caused by a complex comprised of 11 recognized virus species exhibiting accelerated rates of evolution, driven by high frequencies of mutation and genetic exchange. Here, we present a systematic analysis of the contribution of genetic exchange to cassava mosaic virus species-level diversity. Most of these species emerged as a result of genetic exchange. This is the first study to report the significant impact of genetic exchange on speciation in a group of viruses.

Engineered Bacteriophage Therapeutics: Rationale, Challenges and Future

The current problems with increasing bacterial resistance to antibacterial therapies, resulting in a growing frequency of incurable bacterial infections, necessitates the acceleration of studies on antibacterials of a new generation that could offer an alternative to antibiotics or support their action. Bacteriophages (phages) can kill antibiotic-sensitive as well as antibiotic-resistant bacteria, and thus are a major subject of such studies. Their efficacy in curing bacterial infections has been demonstrated in in vivo experiments and in the clinic. Unlike antibiotics, phages have a narrow range of specificity, which makes them safe for commensal microbiota. However, targeting even only the most clinically relevant strains of pathogenic bacteria requires large collections of well characterized phages, whose specificity would cover all such strains. The environment is a rich source of diverse phages, but due to their complex relationships with bacteria and safety concerns, only some naturally occurring phages can be considered for therapeutic applications. Still, their number and diversity make a detailed characterization of all potentially promising phages virtually impossible. Moreover, no single phage combines all the features required of an ideal therapeutic agent. Additionally, the rapid acquisition of phage resistance by bacteria may make phages already approved for therapy ineffective and turn the search for environmental phages of better efficacy and new specificity into an endless race. An alternative strategy for acquiring phages with desired properties in a short time with minimal cost regarding their acquisition, characterization, and approval for therapy could be based on targeted genome modifications of phage isolates with known properties. The first example demonstrating the potential of this strategy in curing bacterial diseases resistant to traditional therapy is the recent successful treatment of a progressing disseminated Mycobacterium abscessus infection in a teenage patient with the use of an engineered phage. In this review, we briefly present current methods of phage genetic engineering, highlighting their advantages and disadvantages, and provide examples of genetically engineered phages with a modified host range, improved safety or antibacterial activity, and proven therapeutic efficacy. We also summarize novel uses of engineered phages not only for killing pathogenic bacteria, but also for in situ modification of human microbiota to attenuate symptoms of certain bacterial diseases and metabolic, immune, or mental disorders.

Recurrent evolution of high virulence in isolated populations of a DNA virus

Hosts and viruses are constantly evolving in response to each other: as a host attempts to suppress a virus, the virus attempts to evade and suppress the host's immune system. Here, we describe the recurrent evolution of a virulent strain of a DNA virus, which infects multiple Drosophila species. Specifically, we identified two distinct viral types that differ 100-fold in viral titer in infected individuals, with similar differences observed in multiple species. Our analysis suggests that one of the viral types recurrently evolved at least four times in the past ~30,000 years, three times in Arizona and once in another geographically distinct species. This recurrent evolution may be facilitated by an effective mutation rate which increases as each prior mutation increases viral titer and effective population size. The higher titer viral type suppresses the host-immune system and an increased virulence compared to the low viral titer type.

ssDNA viruses: key players in global virome

Single-stranded (ss)DNA viruses are extremely widespread, infect diverse hosts from all three domains of life and include important pathogens. Most ssDNA viruses possess small genomes that replicate by the rolling-circle-like mechanism initiated by a distinct virus-encoded endonuclease. High throughput genome sequencing and improved bioinformatics tools have yielded vast information on presence of ssDNA viruses in diverse habitats. The simple genome of ssDNA viruses have high propensity to undergo mutation and recombination often emerging as threat to human civilization. Interestingly their genome is found embedded in fossils dating back to million years. The unusual evolutionary history of ssDNA viruses reveal evidences of horizontal gene transfer, sometimes between different species and genera.

Mutation-Driven Parallel Evolution during Viral Adaptation

Convergent evolution has been demonstrated across all levels of biological organization, from parallel nucleotide substitutions to convergent evolution of complex phenotypes, but whether instances of convergence are the result of selection repeatedly finding the same optimal solution to a recurring problem or are the product of mutational biases remains unsettled. We generated 20 replicate lineages allowed to fix a single mutation from each of four bacteriophage genotypes under identical selective regimes to test for parallel changes within and across genotypes at the levels of mutational effect distributions and gene, protein, amino acid, and nucleotide changes. All four genotypes shared a distribution of beneficial mutational effects best approximated by a distribution with a finite upper bound. Parallel adaptation was high at the protein, gene, amino acid, and nucleotide levels, both within and among phage genotypes, with the most common first-step mutation in each background fixing on an average in 7 of 20 replicates and half of the substitutions in two of the four genotypes occurring at shared sites. Remarkably, the mutation of largest beneficial effect that fixed for each genotype was never the most common, as would be expected if parallelism were driven by selection. In fact, the mutation of smallest benefit for each genotype fixed in a total of 7 of 80 lineages, equally as often as the mutation of largest benefit, leading us to conclude that adaptation was largely mutation-driven, such that mutational biases led to frequent parallel fixation of mutations of suboptimal effect.

Elevating fitness after a horizontal gene exchange in bacteriophage φX174

In an earlier study, protein-based barriers to horizontal gene transfer were investigated by placing the bacteriophage G4 G gene, encoding the major spike protein, into the φX174 genome. The foreign G protein promoted off-pathway assembly reactions, resulting in a lethal phenotype. After three targeted genetic selections, one of two foreign spike proteins was productively integrated into the φX174 system: the complete G4 or a recombinant G4/φX174 protein (94% G4:6% φX174). However, strain fitness was very low. In this study, the chimeras were characterized and experimentally evolved. Inefficient assembly was the primary contributor to low fitness: accordingly, mutations affecting assembly restored fitness. The spike protein preference of the ancestral and evolved strains was determined in competition experiments between the foreign and φX174G proteins. Before adaptation, both G proteins were incorporated into virions; afterwards, the foreign proteins were strongly preferred. Thus, a previously inhibitory protein became the preferred substrate during assembly.

The microviridae: Diversity, assembly, and experimental evolution

The Microviridae, comprised of ssDNA, icosahedral bacteriophages, are a model system for studying morphogenesis and the evolution of assembly. Historically limited to the φX174-like viruses, recent results demonstrate that this richly diverse family is broadly divided into two groups. The defining feature appears to be whether one or two scaffolding proteins are required for assembly. The single-scaffolding systems contain an internal scaffolding protein, similar to many dsDNA viruses, and have a more complex coat protein fold. The two-scaffolding protein systems (φX174-like) encode an internal and external species, as well as an additional structural protein: a spike on the icosahedral vertices. Here, we discuss recent in silico and in vivo evolutionary analyses conducted with chimeric viruses and/or chimeric proteins. The results suggest 1) how double scaffolding systems can evolve into single and triple scaffolding systems; and 2) how assembly is the critical factor governing adaptation and the maintenance of species boundaries.

Synergistic Pleiotropy Overrides the Costs of Complexity in Viral Adaptation

Adaptive evolution progresses as a series of steps toward a multidimensional phenotypic optimum, and organismal or environmental complexity determines the number of phenotypic dimensions, or traits, under selection. Populations evolving in complex environments may experience costs of complexity such that improvement in one or more traits is impeded by selection on others. We compared the fitness effects of the first fixed mutations for populations of single-stranded DNA bacteriophage evolving under simple selection for growth rate to those of populations evolving under more complex selection for growth rate as well as capsid stability. We detected a cost of complexity manifested as a smaller growth rate improvement for mutations fixed under complex conditions. We found that, despite imposing a cost for growth rate improvement, strong complex selection resulted in the greatest overall fitness improvement, even for single mutations. Under weaker secondary selective pressures, tradeoffs between growth rate and stability were pervasive, but strong selection on the secondary trait resulted largely in mutations beneficial to both traits. Strength of selection therefore determined the nature of pleiotropy governing observed trait evolution, and strong positive selection forced populations to find mutations that improved multiple traits, thereby overriding costs incurred as a result of a more complex selective environment. The costs of complexity, however, remained substantial when considering the effects on a single trait in the context of selection on multiple traits.

Intergenic incompatibilities reduce fitness in hybrids of extremely closely related bacteriophages

Horizontal gene transfer and recombination occur across many groups of viruses and play key roles in important viral processes such as host-range expansion and immune-system avoidance. To have any predictive power regarding the ability of viruses to readily recombine, we must determine the extent to which epistasis restricts the success of recombinants, particularly as it relates to the genetic divergence between parental strains. In any hybridization event, the evolutionary success or failure of hybrids is largely determined by the pervasiveness of epistasis in the parental genomes. Recombination has previously been shown to incur steep fitness costs in highly divergent viruses as a result of disrupted epistatic interactions. We used a pair of bacteriophages of the family Microviridae to demonstrate that epistasis may evidence itself in the form of fitness costs even in the case of the exchange of alleles at a locus with amino acid divergence as low as 1%. We explored a possible biophysical source of epistasis in the interaction of viral coat and scaffolding proteins and examined a recovery mutation that likely repairs interactions disrupted by recombination.

The Kinetic and Thermodynamic Aftermath of Horizontal Gene Transfer Governs Evolutionary Recovery

Shared host cells can serve as melting pots for viral genomes, giving many phylogenies a web-like appearance due to horizontal gene transfer. However, not all virus families exhibit web-like phylogenies. Microviruses form three distinct clades, represented by φX174, G4, and α3. Here, we investigate protein-based barriers to horizontal gene transfer between clades. We transferred gene G, which encodes a structural protein, between φX174 and G4, and monitored the evolutionary recovery of the resulting chimeras. In both cases, particle assembly was the major barrier after gene transfer. The G4φXG chimera displayed a temperature-sensitive assembly defect that could easily be corrected through single mutations that promote productive assembly. Gene transfer in the other direction was more problematic. The initial φXG4G chimera required an exogenous supply of both the φX174 major spike G and DNA pilot H proteins. Elevated DNA pilot protein levels may be required to compensate for off-pathway reactions that may have become thermodynamically and/or kinetically favored when the foreign spike protein was present. After three targeted genetic selections, the foreign spike protein was productively integrated into the φX174 background. The first adaption involved a global decrease in gene expression. This was followed by modifications affecting key protein-protein interactions that govern assembly. Finally, gene expression was re-elevated. Although the first selection suppresses nonproductive reactions, subsequent selections promote productive assembly and ultimately viability. However, viable chimeric strains exhibited reduced fitness compared with wild-type. This chimera's path to recovery may partially explain how unusual recombinant viruses could persist long enough to naturally emerge.

Contrasting genetic structure between two begomoviruses infecting the same leguminous hosts

Begomoviruses are whitefly-transmitted, ssDNA plant viruses and are among the most damaging pathogens causing epidemics in economically important crops worldwide. Wild/non-cultivated plants play a crucial epidemiological role, acting as begomovirus reservoirs and as ‘mixing vessels' where recombination can occur. Previous work suggests a higher degree of genetic variability in begomovirus populations from non-cultivated hosts compared with cultivated hosts. To assess this supposed host effect on the genetic variability of begomovirus populations, cultivated (common bean, Phaseolus vulgaris, and lima bean, Phaseolus lunatus) and non-cultivated (Macroptilium lathyroides) legume hosts were sampled from two regions of Brazil. A total of 212 full-length DNA-A genome segments were sequenced from samples collected between 2005 and 2012, and populations of the begomoviruses Bean golden mosaic virus (BGMV) and Macroptilium yellow spot virus (MaYSV) were obtained. We found, for each begomovirus species, similar genetic variation between populations infecting cultivated and non-cultivated hosts, indicating that the presumed genetic variability of the host did not a priori affect viral variability. We observed a higher degree of genetic variation in isolates from MaYSV populations than BGMV populations, which was explained by numerous recombination events in MaYSV. MaYSV and BGMV showed distinct distributions of genetic variation, with the BGMV population (but not MaYSV) being structured by both host and geography.

Extensive recombination-induced disruption of genetic interactions is highly deleterious but can be partially reversed by small numbers of secondary recombination events

Although homologous recombination can potentially provide viruses with vastly more evolutionary options than are available through mutation alone, there are considerable limits on the adaptive potential of this important evolutionary process. Primary among these is the disruption of favorable coevolved genetic interactions that can occur following the transfer of foreign genetic material into a genome. Although the fitness costs of such disruptions can be severe, in some cases they can be rapidly recouped by either compensatory mutations or secondary recombination events. Here, we used a maize streak virus (MSV) experimental model to explore both the extremes of recombination-induced genetic disruption and the capacity of secondary recombination to adaptively reverse almost lethal recombination events. Starting with two naturally occurring parental viruses, we synthesized two of the most extreme conceivable MSV chimeras, each effectively carrying 182 recombination breakpoints and containing thorough reciprocal mixtures of parental polymorphisms. Although both chimeras were severely defective and apparently noninfectious, neither had individual movement-, encapsidation-, or replication-associated genome regions that were on their own "lethally recombinant." Surprisingly, mixed inoculations of the chimeras yielded symptomatic infections with viruses with secondary recombination events. These recombinants had only 2 to 6 breakpoints, had predominantly inherited the least defective of the chimeric parental genome fragments, and were obviously far more fit than their synthetic parents. It is clearly evident, therefore, that even when recombinationally disrupted virus genomes have extremely low fitness and there are no easily accessible routes to full recovery, small numbers of secondary recombination events can still yield tremendous fitness gains. Importance: Recombination between viruses can generate strains with enhanced pathological properties but also runs the risk of producing hybrid genomes with decreased fitness due to the disruption of favorable genetic interactions. Using two synthetic maize streak virus genome chimeras containing alternating genome segments derived from two natural viral strains, we examined both the fitness costs of extreme degrees of recombination (both chimeras had 182 recombination breakpoints) and the capacity of secondary recombination events to recoup these costs. After the severely defective chimeras were introduced together into a suitable host, viruses with between 1 and 3 secondary recombination events arose, which had greatly increased replication and infective capacities. This indicates that even in extreme cases where recombination-induced genetic disruptions are almost lethal, and 91 consecutive secondary recombination events would be required to reconstitute either one of the parental viruses, moderate degrees of fitness recovery can be achieved through relatively small numbers of secondary recombination events.

The adaptive potential of hybridization demonstrated with bacteriophages

The success or failure of hybrids and the factors that determine their fitness have ecological, evolutionary, medical, and economic implications. Hybrid fitness is a major determinant of the size of hybrid zones and the maintenance of related species with overlapping ranges. It also influences the evolution of emerging pathogens and the success of economically important crop species experimentally hybridized in search of strains with increased yields or disease resistance. Hybrid fitness may largely be determined by the pervasiveness of epistasis in the genome, as epistasis is known to debilitate hybrids through disrupted inter- and intragenic interactions. We identified two bacteriophages isolated from their natural environment, one the result of a past hybridization event involving an ancestor of the other phage and a third, unknown phage. By performing a reciprocal cross of the affected region of the genome, consisting of a single complete gene, we both approximately recreated and reversed this original hybridization event in two chimeric bacteriophage genomes. Subsequent adaptation of the hybrid phages allowed for the recovery of fitness losses incurred by the hybrid genotypes. Furthermore, adaptation led to the ascension of a substantially higher and previously inaccessible adaptive peak. We show that by allowing genotypes to take large leaps across the adaptive landscape rather than single mutational steps, hybridization can lead to huge long-term fitness gains in spite of short-term costs resulting from disrupted epistatic interactions, demonstrating that the success or failure of hybrids may be determined not by their initial fitness, but rather by their adaptive potential.

Patterns of recombination in HIV-1M are influenced by selection disfavouring the survival of recombinants with disrupted genomic RNA and protein structures

Genetic recombination is a major contributor to the ongoing diversification of HIV. It is clearly apparent that across the HIV-genome there are defined recombination hot and cold spots which tend to co-localise both with genomic secondary structures and with either inter-gene boundaries or intra-gene domain boundaries. There is also good evidence that most recombination breakpoints that are detectable within the genes of natural HIV recombinants are likely to be minimally disruptive of intra-protein amino acid contacts and that these breakpoints should therefore have little impact on protein folding. Here we further investigate the impact on patterns of genetic recombination in HIV of selection favouring the maintenance of functional RNA and protein structures. We confirm that chimaeric Gag p24, reverse transcriptase, integrase, gp120 and Nef proteins that are expressed by natural HIV-1 recombinants have significantly lower degrees of predicted folding disruption than randomly generated recombinants. Similarly, we use a novel single-stranded RNA folding disruption test to show that there is significant, albeit weak, evidence that natural HIV recombinants tend to have genomic secondary structures that more closely resemble parental structures than do randomly generated recombinants. These results are consistent with the hypothesis that natural selection has acted both in the short term to purge recombinants with disrupted RNA and protein folds, and in the longer term to modify the genome architecture of HIV to ensure that recombination prone sites correspond with those where recombination will be minimally deleterious.

Identifying geographic hot spots of reassortment in a multipartite plant virus

Reassortment between different species or strains plays a key role in the evolution of multipartite plant viruses and can have important epidemiological implications. Identifying geographic locations where reassortant lineages are most likely to emerge could be a valuable strategy for informing disease management and surveillance efforts. We developed a predictive framework to identify potential geographic hot spots of reassortment based upon spatially explicit analyses of genome constellation diversity. To demonstrate the utility of this approach, we examined spatial variation in the potential for reassortment among Cardamom bushy dwarf virus (CBDV; Nanoviridae, Babuvirus) isolates in Northeast India. Using sequence data corresponding to six discrete genome components for 163 CBDV isolates, a quantitative measure of genome constellation diversity was obtained for locations across the sampling region. Two key areas were identified where viruses with highly distinct genome constellations cocirculate, and these locations were designated as possible geographic hot spots of reassortment, where novel reassortant lineages could emerge. Our study demonstrates that the potential for reassortment can be spatially dependent in multipartite plant viruses and highlights the use of evolutionary analyses to identify locations which could be actively managed to facilitate the prevention of outbreaks involving novel reassortant strains.

Environment determines epistatic patterns for a ssDNA virus

Despite the accumulation of substantial quantities of information about epistatic interactions among both deleterious and beneficial mutations in a wide array of experimental systems, neither consistent patterns nor causal explanations for these interactions have yet emerged. Furthermore, the effects of mutations depend on the environment in which they are characterized, implying that the environment may also influence epistatic interactions. Recent work with beneficial mutations for the single-stranded DNA bacteriophage ID11 demonstrated that interactions between pairs of mutations could be understood by means of a simple model that assumes that mutations have additive phenotypic effects and that epistasis arises through a nonlinear phenotype-fitness map with a single intermediate optimum. To determine whether such a model could also explain changes in epistatic patterns associated with changes in environment, we measured epistatic interactions for these same mutations under conditions for which we expected to find the wild-type ID11 at different distances from its phenotypic optimum by assaying fitnesses at three different temperatures: 33°, 37°, and 41°. Epistasis was present and negative under all conditions, but became more pronounced as temperature increased. We found that the additive-phenotypes model explained these patterns as changes in the parameters of the phenotype-fitness map, but that a model that additionally allows the phenotypes to vary across temperatures performed significantly better. Our results show that ostensibly complex patterns of fitness effects and epistasis across environments can be explained by assuming a simple structure for the genotype-phenotype relationship.

Whole-genome synthesis and characterization of viable S13-like bacteriophages

Background

Unprecedented progresses in high-throughput DNA sequencing and de novo gene synthesis technologies have allowed us to create living organisms in the absence of natural template.

Methodology/Principal Findings

The sequence of wild-type S13 phage genome was downloaded from GenBank. Two synonymous mutations were introduced into wt-S13 genome to generate m1-S13 genome. Another mutant, m2-S13 genome, was obtained by engineering two nonsynonymous mutations in the capsid protein coding region of wt-S13 genome. A chimeric phage genome was designed by replacing the F capsid protein open reading frame (ORF) from phage S13 with the F capsid protein ORF from phage G4. The whole genomes of all four phages were assembled from a series of chemically synthesized short overlapping oligonucleotides. The linear synthesized genomes were circularized and electroporated into E.coli C, the standard laboratory host of S13 phage. All four phages were recovered and plaques were visualized. The results of sequencing showed the accuracy of these synthetic genomes. The synthetic phages were capable of lysing their bacterial host and tolerating general environmental conditions. While no phenotypic differences among the variant strains were observed when grown in LB medium with CaCl₂, the S13/G4 chimera was found to be much more sensitive to the absence of calcium and to have a lower adsorption rate under calcium free condition.

Conclusions/Significance

The bacteriophage S13 and its variants can be chemically synthesized. The major capsid gene of phage G4 is functional in the phage S13 life cycle. These results support an evolutional hypothesis which has been proposed that a homologous recombination event involving gene F of quite divergent ancestral lineages should be included in the history of the microvirid family.

Evolutionary recovery of a recombinant viral genome

It is well appreciated that the evolutionary divergence of genes and genomes from a common ancestor ultimately leads to incompatibilities if those genomes are hybridized. Far less is known about the ability and nature of compensatory evolution to yield the recovery of function in hybrid genomes. Here the major capsid gene of the bacteriophage T7 (40-kb dsDNA) was replaced with the homologous gene of either T3 or K11, each 22% different at the protein level from the T7 homolog. Initial fitness was moderately impaired for the T3 exchange, but the K11 exchange was not viable without a compensatory change in the T7 scaffolding protein. Subsequent adaptation of the transgenic phages led to nearly complete fitness recoveries. Compensatory changes were few, mostly in the transgene and its main interacting partner, the scaffolding protein gene. The large magnitude of fitness recovery with relatively few mutations suggests that the fitness costs of hybridizations and horizontal gene exchanges between moderately diverged genomes can potentially be short-lived through compensatory evolution.

Complex recombination patterns arising during geminivirus coinfections preserve and demarcate biologically important intra-genome interaction networks

Genetic recombination is an important process during the evolution of many virus species and occurs particularly frequently amongst begomoviruses in the single stranded DNA virus family, Geminiviridae. As in many other recombining viruses it is apparent that non-random recombination breakpoint distributions observable within begomovirus genomes sampled from nature are the product of variations both in basal recombination rates across genomes and in the over-all viability of different recombinant genomes. Whereas factors influencing basal recombination rates might include local degrees of sequence similarity between recombining genomes, nucleic acid secondary structures and genomic sensitivity to nuclease attack or breakage, the viability of recombinant genomes could be influenced by the degree to which their co-evolved protein-protein and protein-nucleotide and nucleotide-nucleotide interactions are disreputable by recombination. Here we investigate patterns of recombination that occur over 120 day long experimental infections of tomato plants with the begomoviruses Tomato yellow leaf curl virus and Tomato leaf curl Comoros virus. We show that patterns of sequence exchange between these viruses can be extraordinarily complex and present clear evidence that factors such as local degrees of sequence similarity but not genomic secondary structure strongly influence where recombination breakpoints occur. It is also apparent from our experiment that over-all patterns of recombination are strongly influenced by selection against individual recombinants displaying disrupted intra-genomic interactions such as those required for proper protein and nucleic acid folding. Crucially, we find that selection favoring the preservation of co-evolved longer-range protein-protein and protein DNA interactions is so strong that its imprint can even be used to identify the exact sequence tracts involved in these interactions.

Genetic recombination between viruses is a form of parasexual reproduction during which two parental viruses each contribute genetic information to an offspring, or recombinant, virus. Unlike with sexual reproduction, however, recombination in viruses can even involve the transfer of sequences between the members of distantly related species. When parental genomes are very distantly related, it is anticipated that recombination between them runs the risk of producing defective offspring. The reason for this is that the interactions between different parts of genomes and the proteins they encode (such as between different viral proteins or between viral proteins and the virus genomic DNA or RNA) often depend on particular co-evolved binding sites that recognize one another. When in a recombinant genome the partners in a binding site pair are each inherited from different parents there is a possibility that they will not interact with one another properly. Here we examine recombinant genomes arising during experimental mixed infections of two distantly related viruses to detect evidence that intra-genome interaction networks are broadly preserved in these genomes. We show this preservation is so strict that patterns of recombination in these viruses can even be used to identify the interacting regions within their genomes.

Recombination in eukaryotic single stranded DNA viruses

Although single stranded (ss) DNA viruses that infect humans and their domesticated animals do not generally cause major diseases, the arthropod borne ssDNA viruses of plants do, and as a result seriously constrain food production in most temperate regions of the world. Besides the well known plant and animal-infecting ssDNA viruses, it has recently become apparent through metagenomic surveys of ssDNA molecules that there also exist large numbers of other diverse ssDNA viruses within almost all terrestrial and aquatic environments. The host ranges of these viruses probably span the tree of life and they are likely to be important components of global ecosystems. Various lines of evidence suggest that a pivotal evolutionary process during the generation of this global ssDNA virus diversity has probably been genetic recombination. High rates of homologous recombination, non-homologous recombination and genome component reassortment are known to occur within and between various different ssDNA virus species and we look here at the various roles that these different types of recombination may play, both in the day-to-day biology, and in the longer term evolution, of these viruses. We specifically focus on the ecological, biochemical and selective factors underlying patterns of genetic exchange detectable amongst the ssDNA viruses and discuss how these should all be considered when assessing the adaptive value of recombination during ssDNA virus evolution.

A high-throughput venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) and evidence for pervasive positive selection across toxin classes

Despite causing considerable human mortality and morbidity, animal toxins represent a valuable source of pharmacologically active macromolecules, a unique system for studying molecular adaptation, and a powerful framework for examining structure-function relationships in proteins. Snake venoms are particularly useful in the latter regard as they consist primarily of a moderate number of proteins and peptides that have been found to belong to just a handful of protein families. As these proteins and peptides are produced in dedicated glands, transcriptome sequencing has proven to be an effective approach to identifying the expressed toxin genes. We generated a venom-gland transcriptome for the Eastern Diamondback Rattlesnake (Crotalus adamanteus) using Roche 454 sequencing technology. In the current work, we focus on transcripts encoding toxins. We identified 40 unique toxin transcripts, 30 of which have full-length coding sequences, and 10 have only partial coding sequences. These toxins account for 24% of the total sequencing reads. We found toxins from 11 previously described families of snake-venom toxins and have discovered two putative, previously undescribed toxin classes. The most diverse and highly expressed toxin classes in the C. adamanteus venom-gland transcriptome are the serine proteinases, metalloproteinases, and C-type lectins. The serine proteinases are the most abundant class, accounting for 35% of the toxin sequencing reads. Metalloproteinases are the most diverse; 11 different forms have been identified. Using our sequences and those available in public databases, we detected positive selection in seven of the eight toxin families for which sufficient sequences were available for the analysis. We find that the vast majority of the genes that contribute directly to this vertebrate trait show evidence for a role for positive selection in their evolutionary history.