Mapping natural polymorphisms of hepatitis C virus NS3/4A protease and antiviral resistance to inhibitors in worldwide isolates

BACKGROUND

Several inhibitors for the hepatitis C virus (HCV) NS3/4A protease are under development. Although previous studies identified viral resistance mutations, there is little information on the natural variability of proteases from the different viral subtypes. Here, we aimed to determine both the natural variability and presence of resistance or compensatory mutations to new protease inhibitors (PI) in NS3/4A proteases from worldwide HCV isolates.

METHODS

A comprehensive analysis was performed in 380 HCV NS3 sequences (275 genotype 1; 105 other genotypes) from public HCV databases (EuHCVdb and Los Alamos). Amino acid polymorphism and signature patterns were deduced in the protease domain, including all sites associated with resistance to the PIs BILN-2061, Telaprevir (VX-950), Boceprevir (SCH-503034), SCH-6 and ITMN-191.

RESULTS

Few of the residues in the catalytic triad or in substrate/metal-binding sites were polymorphic, and were identified in only 4/380 isolates. However, a relevant polymorphism was found in sites associated either with resistance to PI (V36, 1170 and D168) or with compensatory mutations (171, T72, Q86 and 1153). Furthermore, some unique genotype-specific signature patterns associated with resistance to PI were also identified.

CONCLUSIONS

We describe for the first time the relevant natural polymorphisms of the HCV NS3/4A protease in worldwide isolates. Although the prevalence of major resistance mutations is very low, many compensatory sites are naturally polymorphic among proteases from several HCV subtypes. These data will help to determine whether HCV resistance is likely to be selected with new PIs and will aid the design of genotypic resistance testing.

Using evolutionary tools to refine the new hypervariable region 3 within the envelope 2 protein of hepatitis C virus

The envelope 2 protein of hepatitis C virus (HCV) presents three hypervariable regions, named HVR1, HVR2 and HVR3, in which the presence of antigenic sites has been described. Genetic variability in these regions may reflect the generation of escape mutants as a consequence of the immune response. Therefore, these regions would tend to accumulate amino acid changes along the infection process, an effect that could be accelerated by antiviral treatments. In this study, we have analyzed the E1-E2 region of 23 HCV patients non-responders to antiviral treatment, 7 of which were infected with subtype 1a, 15 with subtype 1b, and 1 with a new HCV-1 subtype, before and after 6 and/or 12 months of peg-interferon+ribavirin treatment. We have sequenced about 100 clones from each sample, analyzing a total of 4906 sequences. A detailed analysis of the evolutionary forces acting along the genome region studied confirmed the existence of the three hypervariable regions, characterized by significant changes in amino acid composition between samples taken at different times from the same patient and a high number of sites evolving under positive selection. Moreover, for the recently described HVR3, our results suggest that its location could be restricted to residues 434-450, instead of the originally postulated 431-466.

Structural analyses of a hypothetical minimal metabolism

By integrating data from comparative genomics and large-scale deletion studies, we previously proposed a minimal gene set comprising 206 protein-coding genes. To evaluate the consistency of the metabolism encoded by such a minimal genome, we have carried out a series of computational analyses. Firstly, the topology of the minimal metabolism was compared with that of the reconstructed networks from natural bacterial genomes. Secondly, the robustness of the metabolic network was evaluated by simulated mutagenesis and, finally, the stoichiometric consistency was assessed by automatically deriving the steady-state solutions from the reaction set. The results indicated that the proposed minimal metabolism presents stoichiometric consistency and that it is organized as a complex power-law network with topological parameters falling within the expected range for a natural metabolism of its size. The robustness analyses revealed that most random mutations do not alter the topology of the network significantly, but do cause significant damage by preventing the synthesis of several compounds or compromising the stoichiometric consistency of the metabolism. The implications that these results have on the origins of metabolic complexity and the theoretical design of an artificial minimal cell are discussed.

The frontier between cell and organelle: genome analysis of Candidatus Carsonella ruddii

BACKGROUND

Bacterial symbioses are widespread among insects. The early establishment of such symbiotic associations has probably been one of the key factors for the evolutionary success of insects, since it may have allowed access to novel ecological niches and to new imbalanced food resources, such as plant sap or blood. Several genomes of bacterial endosymbionts of different insect species have been recently sequenced, and their biology has been extensively studied. Recently, the complete genome sequence of Candidatus Carsonella ruddii, considered the primary endosymbiont of the psyllid Pachpsylla venusta, has been published. This genome consists of a circular chromosome of 159,662 bp and has been proposed as the smallest bacterial endosymbiont genome known to date.

RESULTS

The detailed analysis of the gene content of C. ruddii shows that the extensive degradation of the genome is not compatible with its consideration as a mutualistic endosymbiont and, even more, as a living organism. The ability to perform most essential functions for a cell to be considered alive is heavily impaired by the lack of genes involved in DNA replication, transcription and translation. Furthermore, the shortening of genes causes, in some cases, the loss of essential domains and functional residues needed to fulfill such vital functions. In addition, at least half of the pathways towards the biosynthesis of essential amino acids, its proposed symbiotic function, are completely or partially lost.

CONCLUSION

We propose that this strain of C. ruddii can be viewed as a further step towards the degeneration of the former primary endosymbiont and its transformation in a subcellular new entity between living cells and organelles. Although the transition of genes from C. ruddii to the host nucleus has been proposed, the amount of genes that should have been transferred to the germinal line of the insect would be so big that it would be more plausible to consider the implication of the mitochondrial machinery encoded in the insect nucleus. Furthermore, since most genes for the biosynthesis of essential amino acids have also been lost, it is likely that the host depends on another yet unidentified symbiont to complement its deficient diet.

Contribution of insertions and deletions to the variability of hepatitis C virus populations

Little is known about the potential effects of insertions and deletions (indels) on the evolutionary dynamics of hepatitis C virus (HCV). In fact, the consequences of indels on antiviral treatment response are a field of investigation completely unexplored. Here, an extensive sequencing project was undertaken by cloning and sequencing serum samples from 25 patients infected with HCV subtype 1a and 48 patients with subtype 1b. For 23 patients, samples obtained after treatment with alpha interferon plus ribavirin were also available. Two genome fragments containing the hypervariable regions in the envelope 2 glycoprotein and the PKR-BD domain in NS5A were sequenced, yielding almost 16 000 sequences. Our results show that insertions are quite rare, but they are often present in biologically relevant domains of the HCV genome. Moreover, their frequency distributions between different time samples reflect the quasispecies dynamics of HCV populations. Deletions seem to be subject to negative selection.

Selection for robustness in mutagenized RNA viruses

Mutational robustness is defined as the constancy of a phenotype in the face of deleterious mutations. Whether robustness can be directly favored by natural selection remains controversial. Theory and in silico experiments predict that, at high mutation rates, slow-replicating genotypes can potentially outcompete faster counterparts if they benefit from a higher robustness. Here, we experimentally validate this hypothesis, dubbed the "survival of the flattest," using two populations of the vesicular stomatitis RNA virus. Characterization of fitness distributions and genetic variability indicated that one population showed a higher replication rate, whereas the other was more robust to mutation. The faster replicator outgrew its robust counterpart in standard competition assays, but the outcome was reversed in the presence of chemical mutagens. These results show that selection can directly favor mutational robustness and reveal a novel viral resistance mechanism against treatment by lethal mutagenesis.

Phylogenetic signal and functional categories in Proteobacteria genomes

BACKGROUND

A comprehensive evolutionary analysis of bacterial genomes implies to identify the hallmark of vertical and non-vertical signals and to discriminate them from the presence of mere phylogenetic noise. In this report we have addressed the impact of factors like the universal distribution of the genes, their essentiality or their functional role in the cell on the inference of vertical signal through phylogenomic methods.

RESULTS

We have established that supermatrices derived from data sets composed mainly by genes suspected to be essential for bacterial cellular life perform better on the recovery of vertical signal than those composed by widely distributed genes. In addition, we show that the "Transcription" category of genes seems to harbor a better vertical signal than other functional categories. Moreover, the "Poorly characterized" category performs better than other categories related with metabolism or cellular processes.

CONCLUSION

From these results we conclude that different data sets allow addressing different questions in phylogenomic analyses. The vertical signal seems to be more present in essential genes although these also include a significant degree of incongruence. From a functional perspective, as expected, informational genes perform better than operational ones but we have also shown the surprising behavior of poorly annotated genes, which points to their importance in the genome evolution of bacteria.

From phylogenetics to phylogenomics: the evolutionary relationships of insect endosymbiotic gamma-Proteobacteria as a test case

The increasing availability of complete genome sequences and the development of new, faster methods for phylogenetic reconstruction allow the exploration of the set of evolutionary trees for each gene in the genome of any species. This has led to the development of new phylogenomic methods. Here, we have compared different phylogenetic and phylogenomic methods in the analysis of the monophyletic origin of insect endosymbionts from the gamma-Proteobacteria, a hotly debated issue with several recent, conflicting reports. We have obtained the phylogenetic tree for each of the 579 identified protein-coding genes in the genome of the primary endosymbiont of carpenter ants, Blochmannia floridanus, after determining their presumed orthologs in 20 additional Proteobacteria genomes. A reference phylogeny reflecting the monophyletic origin of insect endosymbionts was further confirmed with different approaches, which led us to consider it as the presumed species tree. Remarkably, only 43 individual genes produced exactly the same topology as this presumed species tree. Most discrepancies between this tree and those obtained from individual genes or by concatenation of different genes were due to the grouping of Xanthomonadales with beta-Proteobacteria and not to uncertainties over the monophyly of insect endosymbionts. As previously noted, operational genes were more prone to reject the presumed species tree than those included in information-processing categories, but caution should be exerted when selecting genes for phylogenetic inference on the basis of their functional category assignment. We have obtained strong evidence in support of the monophyletic origin of gamma-Proteobacteria insect endosymbionts by a combination of phylogenetic and phylogenomic methods. In our analysis, the use of concatenated genes has shown to be a valuable tool for analyzing primary phylogenetic signals coded in the genomes. Nevertheless, other phylogenomic methods such as supertree approaches were useful in revealing alternative phylogenetic signals and should be included in comprehensive phylogenomic studies.

The cost of replication fidelity in human immunodeficiency virus type 1

Mutation rates should be governed by at least three evolutionary factors: the need for beneficial mutations, the benefit of minimizing the mutational load and the cost of replication fidelity. RNA viruses show high mutation rates compared with DNA micro-organisms, and recent findings suggest that the cost of fidelity might play a role in the evolution of increased mutation rates. Here, by analysing previously published data from HIV-1 reverse transcriptase in vitro assays, we show a trade-off between enzymatic accuracy and the maximum rate of polymerization, thus providing a biochemical basis for the fitness cost of fidelity in HIV-1. This trade-off seems to be related to inefficient extension of mispairs, which increases fidelity at the expense of the polymerization rate. Since in RNA viruses fast replication is critical for survival, this could impose a high cost of fidelity and favour the evolution of high mutation rates.

Epidemic dynamics of two coexisting hepatitis C virus subtypes

Hepatitis C virus (HCV) infection affects about 3% of the human population. Phylogenetic analyses have grouped its variants into six major genotypes, which have a star-like distribution and several minor subtypes. The most abundant genotype in Europe is the so-called genotype 1, with two prevalent subtypes, 1a and 1b. In order to explain the higher prevalence of subtype 1b over 1a, a large-scale sequence analysis (100 virus clones) has been carried out over 25 patients of both subtypes in two regions of the HCV genome: one comprising hypervariable region 1 and another including the interferon sensitivity-determining region. Neither polymorphism analysis nor molecular variance analysis (attending to intra- and intersubtype differences, age, sex and previous history of antiviral treatment) was able to show any particular difference between subtypes that might account for their different prevalence. Only the demographic history of the populations carrying both subtypes and analysis of molecular variance (AMOVA) for risk practice suggested that the route of transmission may be the most important factor to explain the observed difference.

The effect of co- and superinfection on the adaptive dynamics of vesicular stomatitis virus

In many infectious diseases, hosts are often simultaneously infected with several genotypes of the same pathogen. Much theoretical work has been done on modelling multiple infection dynamics, but empirical evidences are relatively scarce. Previous studies have demonstrated that coinfection allows faster adaptation than single infection in RNA viruses. Here, we use experimental populations of the vesicular stomatitis Indiana virus derived from an infectious cDNA, to show that superinfection dynamics promotes faster adaptation than single infection. In addition, we have analysed two different periodicities of multiple infection, daily and separated 5 days in time. Daily multiple infections allow higher fitness increases than multiple infections taking place every 5 days. We propose that the effect of superinfection on fitness is mainly influenced by the time elapsed between the first and the second infection, since shorter time intervals offer more opportunities to competition between resident and invading populations.

Modular organization in the reductive evolution of protein-protein interaction networks

BACKGROUND

The variation in the sizes of the genomes of distinct life forms remains somewhat puzzling. The organization of proteins into domains and the different mechanisms that regulate gene expression are two factors that potentially increase the capacity of genomes to create more complex systems. High-throughput protein interaction data now make it possible to examine the additional complexity generated by the way that protein interactions are organized.

RESULTS

We have studied the reduction in genome size of Buchnera compared to its close relative Escherichia coli. In this well defined evolutionary scenario, we found that among all the properties of the protein interaction networks, it is the organization of networks into modules that seems to be directly related to the evolutionary process of genome reduction.

CONCLUSION

In Buchnera, the apparently non-random reduction of the modular structure of the networks and the retention of essential characteristics of the interaction network indicate that the roles of proteins within the interaction network are important in the reductive process.

The evolutionary origin of Xanthomonadales genomes and the nature of the horizontal gene transfer process

Determining the influence of horizontal gene transfer (HGT) on phylogenomic analyses and the retrieval of a tree of life is relevant for our understanding of microbial genome evolution. It is particularly difficult to differentiate between phylogenetic incongruence due to noise and that resulting from HGT. We have performed a large-scale, detailed evolutionary analysis of the different phylogenetic signals present in the genomes of Xanthomonadales, a group of Proteobacteria. We show that the presence of phylogenetic noise is not an obstacle to infer past and present HGTs during their evolution. The scenario derived from this analysis and other recently published reports reflect the confounding effects on bacterial phylogenomics of past and present HGT. Although transfers between closely related species are difficult to detect in genome-scale phylogenetic analyses, past transfers to the ancestor of extant groups appear as conflicting signals that occasionally might make impossible to determine the evolutionary origin of the whole genome.

A small microbial genome: the end of a long symbiotic relationship?

Intracellular bacteria are characterized by genome reduction. The 422,434-base pair genome of Buchnera aphidicola BCc, primary endosymbiont of the aphid Cinara cedri, is approximately 200 kilobases smaller than the previously sequenced B. aphidicola genomes. B. aphidicola BCc has lost most metabolic functions, including the ability to synthesize the essential amino acid tryptophan and riboflavin. In addition, most retained genes are evolving rapidly. Possibly, B. aphidicola BCc is losing its symbiotic capacity and is being complemented (and might be replaced) by the highly abundant coexisting secondary symbiont.

Tempo and mode of early gene loss in endosymbiotic bacteria from insects

BACKGROUND

Understanding evolutionary processes that drive genome reduction requires determining the tempo (rate) and the mode (size and types of deletions) of gene losses. In this study, we analysed five endosymbiotic genome sequences of the gamma-proteobacteria (three different Buchnera aphidicola strains, Wigglesworthia glossinidia, Blochmannia floridanus) to test if gene loss could be driven by the selective importance of genes. We used a parsimony method to reconstruct a minimal ancestral genome of insect endosymbionts and quantified gene loss along the branches of the phylogenetic tree. To evaluate the selective or functional importance of genes, we used a parameter that measures the level of adaptive codon bias in E. coli (i.e. codon adaptive index, or CAI), and also estimates of evolutionary rates (Ka) between pairs of orthologs either in free-living bacteria or in pairs of symbionts.

RESULTS

Our results demonstrate that genes lost in the early stages of symbiosis were on average less selectively constrained than genes conserved in any of the extant symbiotic strains studied. These results also extend to more recent events of gene losses (i.e. among Buchnera strains) that still tend to concentrate on genes with low adaptive bias in E. coli and high evolutionary rates both in free-living and in symbiotic lineages. In addition, we analyzed the physical organization of gene losses for early steps of symbiosis acquisition under the hypothesis of a common origin of different symbioses. In contrast with previous findings we show that gene losses mostly occurred through loss of rather small blocks and mostly in syntenic regions between at least one of the symbionts and present-day E. coli.

CONCLUSION

At both ancient and recent stages of symbiosis evolution, gene loss was at least partially influenced by selection, highly conserved genes being retained more readily than lowly conserved genes: although losses might result from drift due to the bottlenecking of endosymbiontic populations, we demonstrated that purifying selection also acted by retaining genes of greater selective importance.

Analysis of the overdispersed clock in the short-term evolution of hepatitis C virus: Using the E1/E2 gene sequences to infer infection dates in a single source outbreak

The assumption of a molecular clock for dating events from sequence information is often frustrated by the presence of heterogeneity among evolutionary rates due, among other factors, to positively selected sites. In this work, our goal is to explore methods to estimate infection dates from sequence analysis. One such method, based on site stripping for clock detection, was proposed to unravel the clocklike molecular evolution in sequences showing high variability of evolutionary rates and in the presence of positive selection. Other alternatives imply accommodating heterogeneity in evolutionary rates at various levels, without eliminating any information from the data. Here we present the analysis of a data set of hepatitis C virus (HCV) sequences from 24 patients infected by a single individual with known dates of infection. We first used a simple criterion of relative substitution rate for site removal prior to a regression analysis. Time was regressed on maximum likelihood pairwise evolutionary distances between the sequences sampled from the source individual and infected patients. We show that it is indeed the fastest evolving sites that disturb the molecular clock and that these sites correspond to positively selected codons. The high computational efficiency of the regression analysis allowed us to compare the site-stripping scheme with random removal of sites. We demonstrate that removing the fast-evolving sites significantly increases the accuracy of estimation of infection times based on a single substitution rate. However, the time-of-infection estimations improved substantially when a more sophisticated and computationally demanding Bayesian method was used. This method was used with the same data set but keeping all the sequence positions in the analysis. Consequently, despite the distortion introduced by positive selection on evolutionary rates, it is possible to obtain quite accurate estimates of infection dates, a result of especial relevance for molecular epidemiology studies.

A new subtype of hepatitis C virus genotype 1: complete genome and phylogenetic relationships of an Equatorial Guinea isolate

Hepatitis C virus (HCV) is the leading cause of chronic liver disease and is associated with hepatocellular carcinoma. However, there have been few studies on the distribution and genetic diversity of HCV isolates in non-developed countries. Here, the complete genome sequence of an HCV genotype 1 isolate from Equatorial Guinea is reported, the first complete HCV-1 genome of African origin. Phylogenetic analysis revealed that this sequence always grouped with sequences of genotype 1, but did not group clearly with any subtype described so far. An analysis of partial NS5B gene sequences with additional sequences of African origin also failed to find close similarities between the new sequence and any previously known isolate. Genetic divergence of the coding region of this new sequence with respect to the recognized subtypes of HCV-1 ranged from 20 to 22%. It is proposed that this isolate is a representative of a new, distinct variant of HCV subtype 1.

Evolutionary analysis of genetic variation observed in citrus tristeza virus (CTV) after host passage

We have studied the genetic variability in two genes (p18 and p20) from two groups of Citrus tristeza virus (CTV) isolates. One group (isolates T385, T317, T318, and T305) was derived from a Spanish source by successive host passages while the other (isolates T388 and T390) was obtained after aphid transmission from a Japanese source. A total of 274 sequences were obtained for gene p18 and 451 for p20. In the corresponding phylogenetic trees, sequences derived from the severe isolates (T318, T305, and T388) clustered together and separately from those derived from mild or moderate isolates (T385, T317, and T390), regardless of their geographic origin. Hierarchical analyses of molecular variance showed that up to 53% of the total genetic variability in p18 and up to 87% of the variation in p20 could be explained by differences in the pathogenicity features of the isolates. Neutrality tests revealed that different selection forces had been acting between isolates and between genes, with purifying selection being suggested for p18 from isolates T385 and T390 and for p20 from isolates T385, T317, and T388, and balancing selection for p18 from isolates T318, T305, and T388 and for p20 from isolates T318 and T390. Furthermore, several models of codon selection were observed, with purifying selection being the most notable one, compatible with low effective population size of the virus populations resulting from transmission bottlenecks. We found no evidence of recombination playing a significant role during p18 and p20 evolution in these isolates. These results suggest that hosts can be an important evolutionary factor for CTV isolates.

Variations in two gene sequences of Citrus tristeza virus after host passage

We estimated genetic variation in two groups of Citrus tristeza virus (CTV) isolates: one of them (isolates T385, T317, T318 and T305) derived from a Spanish source by successive host passages, and the other (isolates T388 and T390), obtained after aphid transmission of a Japanese source. The population structure of these isolates had been characterized by single-strand conformation polymorphism analysis of genes p18 and p20. The nucleotide sequences of representative haplotypes of each isolate and gene were used to estimate genetic diversity within and between isolates and to evaluate genetic differentiation between populations. Phylogenetic analysis of p18 and p20 sequence variants showed two main groups: one them included variants predominant in the severe isolates (T318, T305 and T388), and the other comprised variants present in both mild (T385, T317) and severe isolates. Most sequence variants of isolate T390 were not associated to these groups. In some isolates, within-isolate diversity was higher than diversity with other isolates because their population contained distantly related sequence variants, some of which were genetically close to variants predominant in the second isolate. Isolates T388 and T390 were genetically different for the two genes, as estimated by the F statistic. Furthermore, genetic differentiation between T385 and T317, T318 and T305 increased after each host passage. Our results suggest that aphid transmission and host passage may significantly alter the composition of CTV populations and thus be an important factor in their evolution.

Molecular epidemiology and evolution in an outbreak of fulminant hepatitis B virus

In order to establish the transmission pathway for two outbreak patients affected by fulminant hepatitis B (FHB) following a shared period of hospitalization, we sequenced the complete genomes of the hepatitis B viruses (HBV) isolated from them as well as from the suspected common source and 11 additional controls. Phylogenetic and statistical analyses of these sequences revealed that the two FHB patients were indeed infected by a common source and that the fatal development of the disease did not appear to be associated with any mutation previously reported to be related to FHB. These data have also allowed us to estimate the extent and distribution of genetic variability along the genomes of HBV genotype D samples from the same source population. As a result of these analyses, we provide an improved statistical method to individualize the assignment of each suspected patient and the source of an outbreak and information on which genome region to analyze in the molecular epidemiological assessment of hepatitis B virus transmission cases.

Genetic similarity of hepatitis C virus and fibrosis progression in chronic and recurrent infection after liver transplantation

The effect of hepatitis C virus (HCV) genetic heterogeneity on clinical features of post-transplantation hepatitis C is controversial. Different regions of the HCV genome have been associated with apoptosis, fibrosis, and other pathways leading to liver damage in chronic HCV infection. Besides, differences in immunodominant regions, such as NS3, may influence HCV-specific immune responses and disease outcome. In the liver transplant setting, a recent study has reported a positive association between HCV-1b Core region genetic relatedness 5-year post-transplantation and histological severity of recurrent hepatitis C. We have compared nucleotide sequences of HCV Core, NS3 and NS5b regions in HCV-1b-infected patients 3 years post-transplantation (n = 22). A cohort of nontransplanted patients (n = 22) was used as control of natural chronic HCV-1b infection. Histological evaluation was used to define the rate of fibrosis progression. Molecular variance analysis did not show significant differences in HCV sequences between transplanted and nontransplanted patients, or between those with fast or slow fibrosis progression. The same results were obtained when analysing phylogenetic trees for Core, NS3 and NS5b regions. A more appropriate clustering method (using minimum spanning networks) revealed a significant positive relationship between HCV genetic similarity in Core (r = 0.550, P < 0.01) and NS5b regions (r = 0.847, P < 0.01) and the yearly rate of fibrosis progression in nontransplanted patients which, in contrast, was not observed in transplanted patients. Our results indicate that some strains of HCV-1b might be more pathogenic in the natural course of chronic infection by this virus subtype. In the liver transplant setting, when the immune response is severely compromised, other mechanisms are probably more important in determining hepatitis C progression.

Analysis of and function predictions for previously conserved hypothetical or putative proteins in Blochmannia floridanus

BACKGROUND

There is an increasing interest to better understand endosymbiont capabilities in insects both from an ecological point of view and for pest control. Blochmannia floridanus provides important nutrients for its host, the ant Camponotus, while the bacterium in return is provided with a niche to proliferate. Blochmannia floridanus proteins and metabolites are difficult to study due to its endosymbiontic life style; however, its complete genome sequence became recently available.

RESULTS

Improved sequence analysis algorithms, databanks and gene and pathway context methods allowed us to reveal new information on various enzyme and pathways from the Blochmannia floridanus genome sequence [EMBL-ID BX248583]. Furthermore, these predictions are supported and linked to experimental data for instance from structural genomics projects (e.g. Bfl341, Bfl 499) or available biochemical data on proteins from other species which we show here to be related. We were able to assign a confirmed or at least a putative molecular function for 21 from 27 previously conserved hypothetical proteins. For 48 proteins of 66 with a previous putative assignment the function was further clarified. Several of these proteins occur in many proteobacteria and are found to be conserved even in the compact genome of this endosymbiont. To extend and re-test predictions and links to experimentally verified protein functions, functional clusters and interactions were assembled. These included septum initiation and cell division (Bfl165, Bfl303, Bfl248 et al.); translation; transport; the ubiquinone (Bfl547 et al.), the inositol and nitrogen pathways.

CONCLUSION

Taken together, our data allow a better and more complete description of the pathway capabilities and life style of this typical endosymbiont.

Population dynamics of the 2 major mitochondrial DNA haplotypes in experimental populations of Drosophila subobscura

The evolution of Drosophila subobscura mitochondrial DNA has been studied in experimental populations, founded with flies from a natural population from Calvià (Majorca, Balearic Islands, Spain). This population, like others founded in Europe, is characterized by the presence of 2 very common (>95%) mitochondrial haplotypes (named I and II) and rare and endemic haplotypes that appear at very low frequencies. Four experimental populations were established with flies having a heterogeneous nuclear genetic background, which was representative of the composition of the natural population. The populations were started with haplotypes I and II at an initial frequency of 50% each. After 33 generations, the 2 haplotypes coexisted. Random drift could be rejected as the only force responsible for the observed changes in haplotype frequencies. A slight but significant linear trend favouring a mtDNA (haploid) fitness effect has been detected, with a nonlinear deviation that could be due to a nuclear component. An analysis of chromosomal arrangements was made before the foundations of the cages and at generation 23. Our results indicated that the hypothesis that the maintenance of the frequencies of haplotypes I and II in natural populations could be due to their association with chromosomal arrangements remains controversial.

Chromosomal stasis versus plasmid plasticity in aphid endosymbiont Buchnera aphidicola

The study of three genomes of the aphid endosymbiont Buchnera aphidicola has revealed an extraordinary stasis: conservation of gene order and genetic composition of the chromosome, while the chromosome size and number of genes has reduced. The reduction in genome size appears to be ongoing since some lineages we now know to have even smaller chromosomes than the first B. aphidicola analysed. The current sequencing by our group of one of these smaller genomes with an estimated size of 450 kb, and its comparison with the other three available genomes provide insights into the nature of processes involved in shrinkage. We discuss whether B. aphidicola might be driven to extinction and be replaced by secondary aphid endosymbionts. In some lineages, genes encoding key enzymes in the pathways leading to tryptophan and leucine biosynthesis (trpEG and leuABCD, respectively) are located on plasmids, rather than the main chromosome. In contrast to the stasis of the main chromosome, plasmid genes have frequently been transferred to the main chromosome and undergone other gene rearrangements. We propose a two-step scenario to explain these contrasting modes of evolution. Essential genes may have escaped regulation by moving to plasmids in a moving B. aphidicola ancestor. B. aphidicola became polyploidy at a given stage of its evolution and plasmid genes have been transferred to the main chromosome through several independent events.

The cost of replication fidelity in an RNA virus

It is often argued that high mutation rates are advantageous for RNA viruses, because they confer elevated rates of adaptation. However, there is no direct evidence showing a positive correlation between mutation and adaptation rates among RNA viruses. Moreover, theoretical work does not argue in favor of this prediction. We used a series of vesicular stomatitis virus clones harboring single amino acid substitutions in the RNA polymerase to demonstrate that changes inducing enhanced fidelity paid a fitness cost, but that there was no positive correlation between mutation an adaptation rates. We demonstrate that the observed mutation rate in vesicular stomatitis virus can be explained by a trade-off between replication rate and replication fidelity.

Epistasis and the adaptability of an RNA virus

We have explored the patterns of fitness recovery in the vesicular stomatitis RNA virus. We show that, in our experimental setting, reversions to the wild-type genotype were rare and fitness recovery was at least partially driven by compensatory mutations. We compared compensatory adaptation for genotypes carrying (1) mutations with varying deleterious fitness effects, (2) one or two deleterious mutations, and (3) pairs of mutations showing differences in the strength and sign of epistasis. In all cases, we found that the rate of fitness recovery and the proportion of reversions were positively affected by population size. Additionally, we observed that mutations with large fitness effect were always compensated faster than mutations with small fitness effect. Similarly, compensatory evolution was faster for genotypes carrying a single deleterious mutation than for those carrying pairs of mutations. Finally, for genotypes carrying two deleterious mutations, we found evidence of a negative correlation between the epistastic effect and the rate of compensatory evolution.

Gene expression levels influence amino acid usage and evolutionary rates in endosymbiotic bacteria

Most endosymbiotic bacteria have extremely reduced genomes, accelerated evolutionary rates, and strong AT base compositional bias thought to reflect reduced efficacy of selection and increased mutational pressure. Here, we present a comparative study of evolutionary forces shaping five fully sequenced bacterial endosymbionts of insects. The results of this study were three-fold: (i) Stronger conservation of high expression genes at not just nonsynonymous, but also synonymous, sites. (ii) Variation in amino acid usage strongly correlates with GC content and expression level of genes. This pattern is largely explained by greater conservation of high expression genes, leading to their higher GC content. However, we also found indication of selection favoring GC-rich amino acids that contrasts with former studies. (iii) Although the specific nutritional requirements of the insect host are known to affect gene content of endosymbionts, we found no detectable influence on substitution rates, amino acid usage, or codon usage of bacterial genes involved in host nutrition.

Molecular epidemiology of a hepatitis C virus outbreak in a hemodialysis unit

We analyzed a hepatitis C virus (HCV) transmission case in the hemodialysis unit of a private clinic by sequencing two genome regions of virus isolates from a number of patients attending this unit and some external controls. The analysis of 337 nucleotides (nt) in the NS5B region did not provide enough resolution to ascertain which patients were actually involved in the outbreak and the potential source. Nevertheless, this region allowed the exclusion of several patients as putative sources of the transmission case based on their genotypes and phylogenetic relationships. On the other hand, the analysis of several 472-nt-long clone sequences per sample in a more rapidly evolving region of the HCV genome, coding for the envelope proteins and encompassing hypervariable region 1, allowed us to establish the existence of at least two independent transmission events involving two different source patients and three recipients. The direction of the transmissions was further corroborated by different measures of genetic variability within and among samples.

Adaptive evolution in GroEL from distantly related endosymbiotic bacteria of insects

Many symbioses between bacteria and insects resulted from ancient infections followed by strict vertical transmission within host lineages. The strong bottlenecks under which this transmission occurs promote the neutral fixation of slightly deleterious mutations by genetic drift. As predicted by Muller's ratchet, this fixation will drive endosymbiotic bacteria through an irreversible dynamics of fitness loss. The chaperonin GroEL has been proposed as a compensatory mechanism whereby endosymbiotic bacteria of aphids persist. Here, we show that endosymbiotic bacteria of insects from two phylogenetically very distant bacterial phyla have fixed amino acid substitutions by positive selection in functionally important GroEL regions involved in either GroES/peptide binding or in the en bloc movement of the GroEL apical domain. These results, together with the high levels of constitutive expression of GroEL in these endosymbionts, provide valuable insights into the evolution of a molecular mechanism responsible for the maintenance of the symbiotic lifestyle.

Following the very initial growth of biological RNA viral clones

Due to their extremely high genetic diversity, which is a direct consequence of high mutation rates, RNA viruses are often described as molecular quasispecies. According to this theory, RNA virus populations cannot be understood in terms of individual viral clones, as they are clouds of interconnected mutants, but this prediction has not yet been demonstrated experimentally. The goal of this study was to determine the fitness of individual clones sampled from a given RNA virus population, a necessary previous step to test the above prediction. To do so, limiting dilutions of a vesicular stomatitis virus population were employed to isolate single viral clones and their initial growth dynamics were followed, corresponding to the release of the first few hundred viral particles. This technique is useful for estimating basic fitness parameters, such as intracellular growth rate, viral yield per cell, rate at which cells are infected and time spent in cell-to-cell transmission. A combination of these parameters allows estimation of the fitness of individual clones, which seems to be determined mainly by their ability to complete infection cycles more quickly. Interestingly, fitness was systematically higher for initial clones than for their derived populations. In addition to environmental changes, such as cellular defence mechanisms, these differences are attributable to high RNA virus mutation rates.

Comparative analysis of two genomic regions among four strains of Buchnera aphidicola, primary endosymbiont of aphids

Preliminary analysis of two selected genomic regions of Buchnera aphidicola BCc, the primary endosymbiont of the cedar aphid Cinara cedri, has revealed a number of interesting features when compared with the corresponding homologous regions of the three B. aphidicola genomes previously sequenced, that are associated with different aphid species. Both regions exhibit a significant reduction in length and gene number in B. aphidicola BCc, as it could be expected since it possess the smallest bacterial genome. However, the observed genome reduction is not even in both regions, as it appears to be dependent on the nature of their gene content. The region fpr-trxA, that contains mainly metabolic genes, has lost almost half of its genes (45.6%) and has reduced 52.9% its length. The reductive process in the region rrl-aroK, that contains mainly ribosomal protein genes, is less dramatic, since it has lost 9.3% of genes and has reduced 15.5% of its length. Length reduction is mainly due to the loss of protein-coding genes, not to the shortening of ORFs or intergenic regions. In both regions, G+C content is about 4% lower in BCc than in the other B. aphidicola strains. However, when only conserved genes and intergenic regions of the four B. aphidicola strains are compared, the G+C reduction is higher in the fpr-trxA region.

Validating viral quasispecies with digital organisms: a re-examination of the critical mutation rate

Background

In this report we re-examine some recent experiments with digital organisms to test some predictions of quasispecies theory. These experiments revealed that under high mutation rates populations of less fit organisms previously adapted to such high mutation rates were able to outcompete organisms with higher average fitness but adapted to low mutation rates.

Results

We have verified that these results do hold in the original conditions and, by extending the set of initial parameters, we have also detected that the critical mutation rate was independent of population size, a result that we have found to be dependent on a different, contingent factor, the initial fitness vector. Furthermore, in all but one case, the critical mutation rate is higher than the error threshold, a key parameter in quasispecies theory, which prevents its extrapolation to natural viral populations.

Conclusion

From these results we conclude that digital organisms are useful tools for investigating evolutionary patterns and processes including some predictions from the quasispecies theory.

Mode of selection and experimental evolution of antiviral drugs resistance in vesicular stomatitis virus

The possession of an antiviral resistance mutation benefits a virus when the corresponding antiviral is present. But does the resistant virus pay a fitness cost when the antiviral is absent? Would an evolutionary history of association between a genotype and a resistance mutation overcome this cost by changes compensating the harmful side-effect of resistance mutations? Are combined therapies more effective against the rise of resistant viruses or against evolutionary compensations? To explore all these questions, we took an experimental evolution approach. After selecting vesicular stomatitis virus (VSV) populations able to replicate under increasing concentrations of ribavirin and/or alpha-interferon, we evolved them for more than 100 generations under sub-inhibitory concentrations of the corresponding antivirals, allowing for evolutionary compensation. Our results show that resistance for these two antivirals was not easily achieved, being the selected populations generally less fit than the ancestrals both in presence and absence of the antivirals. Evolution in presence of sub-inhibitory concentrations of antivirals compensated for the reduction in fitness in presence of antiviral therapy.

Effect of oligonucleotide primers in determining viral variability within hosts

BACKGROUND

Genetic variability in viral populations is usually estimated by means of polymerase chain reaction (PCR) based methods in which the relative abundance of each amplicon is assumed to be proportional to the frequency of the corresponding template in the initial sample. Although bias in template-to-product ratios has been described before, its relevance in describing viral genetic variability at the intrapatient level has not been fully assessed yet.

RESULTS

To investigate the role of oligonucleotide design in estimating viral variability within hosts, genetic diversity in hepatitis C virus (HCV) populations from eight infected patients was characterised by two parallel PCR amplifications performed with two slightly different sets of primers, followed by cloning and sequencing (mean = 89 cloned sequences per patient). Population genetics analyses of viral populations recovered by pairs of amplifications revealed that in seven patients statistically significant differences were detected between populations sampled with different set of primers.

CONCLUSIONS

Genetic variability analyses demonstrates that PCR selection due to the choice of primers, differing in their degeneracy degree at some nucleotide positions, can eclipse totally or partially viral variants, hence yielding significant different estimates of viral variability within a single patient and therefore eventually producing quite different qualitative and quantitative descriptions of viral populations within each host.

Genomic determinants of protein folding thermodynamics in prokaryotic organisms

Here we investigate how thermodynamic properties of orthologous proteins are influenced by the genomic environment in which they evolve. We performed a comparative computational study of 21 protein families in 73 prokaryotic species and obtained the following main results. (i) Protein stability with respect to the unfolded state and with respect to misfolding are anticorrelated. There appears to be a trade-off between these two properties, which cannot be optimized simultaneously. (ii) Folding thermodynamic parameters are strongly correlated with two genomic features, genome size and G+C composition. In particular, the normalized energy gap, an indicator of folding efficiency in statistical mechanical models of protein folding, is smaller in proteins of organisms with a small genome size and a compositional bias towards A+T. Such genomic features are characteristic for bacteria with an intracellular lifestyle. We interpret these correlations in light of mutation pressure and natural selection. A mutational bias toward A+T at the DNA level translates into a mutational bias toward more hydrophobic (and in general more interactive) proteins, a consequence of the structure of the genetic code. Increased hydrophobicity renders proteins more stable against unfolding but less stable against misfolding. Proteins with high hydrophobicity and low stability against misfolding occur in organisms with reduced genomes, like obligate intracellular bacteria. We argue that they are fixed because these organisms experience weaker purifying selection due to their small effective population sizes. This interpretation is supported by the observation of a high expression level of chaperones in these bacteria. Our results indicate that the mutational spectrum of a genome and the strength of selection significantly influence protein folding thermodynamics.

Bacterial endosymbionts of insects: insights from comparative genomics

The development of molecular techniques for the study of uncultured bacteria allowed the extensive study of the widespread association between insects and intracellular symbiotic bacteria. Most of the bacterial endosymbionts involved in such associations are gamma-proteobacteria, closely related to Escherichia coli. In recent years, five genomes from insect endosymbionts have been sequenced, allowing the performance of extensive genome comparative analysis that, as a complement of phylogenetic studies, and analysis on individual genes, can help to understand the different traits of this particular association, including how the symbiotic process is established, the explanation of the special features of these microbial genomes, the bases of this intimate association and the possible future that awaits the endosymbionts with extremely reduced genomes.

The contribution of epistasis to the architecture of fitness in an RNA virus

The tendency for genetic architectures to exhibit epistasis among mutations plays a central role in the modern synthesis of evolutionary biology and in theoretical descriptions of many evolutionary processes. Nevertheless, few studies unquestionably show whether, and how, mutations typically interact. Beneficial mutations are especially difficult to identify because of their scarcity. Consequently, epistasis among pairs of this important class of mutations has, to our knowledge, never before been explored. Interactions among genome components should be of special relevance in compacted genomes such as those of RNA viruses. To tackle these issues, we first generated 47 genotypes of vesicular stomatitis virus carrying pairs of nucleotide substitution mutations whose separated and combined deleterious effects on fitness were determined. Several pairs exhibited significant interactions for fitness, including antagonistic and synergistic epistasis. Synthetic lethals represented 50% of the latter. In a second set of experiments, 15 genotypes carrying pairs of beneficial mutations were also created. In this case, all significant interactions were antagonistic. Our results show that the architecture of the fitness depends on complex interactions among genome components.

Coexistence of Wolbachia with Buchnera aphidicola and a secondary symbiont in the aphid Cinara cedri

Intracellular symbiosis is very common in the insect world. For the aphid Cinara cedri, we have identified by electron microscopy three symbiotic bacteria that can be characterized by their different sizes, morphologies, and electrodensities. PCR amplification and sequencing of the 16S ribosomal DNA (rDNA) genes showed that, in addition to harboring Buchnera aphidicola, the primary endosymbiont of aphids, C. cedri harbors a secondary symbiont (S symbiont) that was previously found to be associated with aphids (PASS, or R type) and an alpha-proteobacterium that belongs to the Wolbachia genus. Using in situ hybridization with specific bacterial probes designed for symbiont 16S rDNA sequences, we have shown that Wolbachia was represented by only a few minute bacteria surrounding the S symbionts. Moreover, the observed B. aphidicola and the S symbionts had similar sizes and were housed in separate specific bacterial cells, the bacteriocytes. Interestingly, in contrast to the case for all aphids examined thus far, the S symbionts were shown to occupy a similarly sized or even larger bacteriocyte space than B. aphidicola. These findings, along with the facts that C. cedri harbors the B. aphidicola strain with the smallest bacterial genome and that the S symbionts infect all Cinara spp. analyzed so far, suggest the possibility of bacterial replacement in these species.

Fitness and life-history traits of the two major mitochondrial DNA haplotypes of Drosophila subobscura

Mitochondrial DNA restriction site analyses on natural populations of Drosophila subobscura have proved the existence of two common, coexisting haplotypes (I and II), as well as a set of less frequent ones derived from them. To explain this distribution, experiments to date point practically to all possible genetic mechanisms being involved in the changes of gene frequencies (cytonuclear coadaptation, direct natural selection on mtDNA and genetic drift). In an attempt to find differences that help to understand the dynamics of these haplotypes and to detect the effect of selection, we measured certain fitness components and life-history traits (egg-larva and larva-adult viabilities and developmental times, longevity, resistance to desiccation and optimal density) of the two main haplotypes I and II when maintained in laboratory population cages. As a general trend, haplotype II showed a higher net fitness than haplotype I, which explains the superiority of haplotype II over haplotype I in experimental populations but not their coexistence in nature, where additional factors must be considered.

GroEL and the maintenance of bacterial endosymbiosis

Many eukaryotic organisms have symbiotic associations with obligate intracellular bacteria. The clonal transmission of endosymbionts between host generations should lead to the irreversible fixation of slightly deleterious mutations in their non-recombinant genome by genetic drift. However, the stability of endosymbiosis indicates that some mechanism is involved in the amelioration of the effects of these mutations. We propose that the chaperone GroEL was involved in the acquisition of an endosymbiotic lifestyle not only by means of its over-production, as proposed by Moran, but also by its adaptive evolution mediated by positive selection to improve the interaction with the unstable endosymbiont proteome.

Determination of the core of a minimal bacterial gene set

The availability of a large number of complete genome sequences raises the question of how many genes are essential for cellular life. Trying to reconstruct the core of the protein-coding gene set for a hypothetical minimal bacterial cell, we have performed a computational comparative analysis of eight bacterial genomes. Six of the analyzed genomes are very small due to a dramatic genome size reduction process, while the other two, corresponding to free-living relatives, are larger. The available data from several systematic experimental approaches to define all the essential genes in some completely sequenced bacterial genomes were also considered, and a reconstruction of a minimal metabolic machinery necessary to sustain life was carried out. The proposed minimal genome contains 206 protein-coding genes with all the genetic information necessary for self-maintenance and reproduction in the presence of a full complement of essential nutrients and in the absence of environmental stress. The main features of such a minimal gene set, as well as the metabolic functions that must be present in the hypothetical minimal cell, are discussed.

The distribution of fitness effects caused by single-nucleotide substitutions in an RNA virus

Little is known about the mutational fitness effects associated with single-nucleotide substitutions on RNA viral genomes. Here, we used site-directed mutagenesis to create 91 single mutant clones of vesicular stomatitis virus derived from a common ancestral cDNA and performed competition experiments to measure the relative fitness of each mutant. The distribution of nonlethal deleterious effects was highly skewed and had a long, flat tail. As expected, fitness effects depended on whether mutations were chosen at random or reproduced previously described ones. The effect of random deleterious mutations was well described by a log-normal distribution, with -19% reduction of average fitness; the effects distribution of preobserved deleterious mutations was better explained by a beta model. The fit of both models was improved when combined with a uniform distribution. Up to 40% of random mutations were lethal. The proportion of beneficial mutations was unexpectedly high. Beneficial effects followed a gamma distribution, with expected fitness increases of 1% for random mutations and 5% for preobserved mutations.

Natural selection and the organ-specific differentiation of HIV-1 V3 hypervariable region

The existence of organ-specific HIV-1 populations within infected hosts has been studied for many years; nonetheless results reported by different authors are somewhat discrepant. To tackle this problem, we used a population genetics approach to analyze previously published data from the V3 hypervariable region of the envelope env gene. Our results are compatible with a population subdivision by organs in 95% of individuals analyzed at autopsy. In addition, populations infecting the nervous system and testicles clearly appear as differentiated subsets of the so-called macrophage-tropic variants. Liver and kidney may harbor differentiated populations as well. Although it is widely accepted that organ compartmentalization arises as a consequence of different selective pressures imposed by different organs, a definitive demonstration has not yet been provided. Our analysis of the pattern of synonymous and nonsynonymous nucleotide substitutions provides evidence supporting this hypothesis, without discarding the role of other evolutionary processes. In contrast, positive selection does not seem to be the mechanism responsible for the evolution of patient-specific sequences.

Evolution of the leucine gene cluster in Buchnera aphidicola: insights from chromosomal versions of the cluster

In Buchnera aphidicola strains associated with the aphid subfamilies Thelaxinae, Lachninae, Pterocommatinae, and Aphidinae, the four leucine genes (leuA, -B, -C, and -D) are located on a plasmid. However, these genes are located on the main chromosome in B. aphidicola strains associated with the subfamilies Pemphiginae and Chaitophorinae. The sequence of the chromosomal fragment containing the leucine cluster and flanking genes has different positions in the chromosome in B. aphidicola strains associated with three tribes of the subfamily Pemphiginae and one tribe of the subfamily Chaitophorinae. Due to the extreme gene order conservation of the B. aphidicola genomes, the variability in the position of the leucine cluster in the chromosome may be interpreted as resulting from independent insertions from an ancestral plasmid-borne leucine gene. These findings do not support a chromosomal origin for the leucine genes in the ancestral B. aphidicola and do support a back transfer evolutionary scenario from a plasmid to the main chromosome.

The population genetics and evolutionary epidemiology of RNA viruses

The authors discuss the main mechanisms of RNA virus evolution — mutation, recombination, natural selection, genetic drift and migration, and how these interact to shape the genetic structure of populations. The quasispecies model of RNA virus evolution is explained and the question of whether this model provides an accurate description of RNA virus evolution is discussed. Experiments that can be carried out to test the basic principles of evolutionary theory are briefly described. The authors review what such experiments have told us about virus evolution and, more widely, what these experiments have revealed in terms of general evolutionary principles. RNA viruses evolve quickly, so a detailed reconstruction of their epidemiological history can be undertaken. The authors show how epidemiological patterns of viruses result from their evolution at two different levels: within individual hosts (and vectors) and among hosts at the population level. Using several examples, including HIV and SARS, the authors describe how studying RNA virus evolution could be used to understand virus emergence. Finally, the important topics of the evolution of virulence and resistance to drugs are discussed.

RNA viruses are ubiquitous intracellular parasites that are responsible for many emerging diseases, including AIDS and SARS. Here, we discuss the principal mechanisms of RNA virus evolution and highlight areas where future research is required. The rapidity of sequence change in RNA viruses means that they are useful experimental models for the study of evolution in general and it enables us to watch them change in 'real time', and retrace the spread through populations with molecular phylogenies. An understanding of the mechanisms of RNA virus sequence change is also crucial to predicting important aspects of their emergence and long-term evolution.