The pea aphid phylome: a complete catalogue of evolutionary histories and arthropod orthology and paralogy relationships for Acyrthosiphon pisum genes

Phylogenetic analyses serve many purposes, including the establishment of orthology relationships, the prediction of protein function and the detection of important evolutionary events. Within the context of the sequencing of the genome of the pea aphid, Acyrthosiphon pisum, we undertook a phylogenetic analysis for every protein of this species. The resulting phylome includes the evolutionary relationships of all predicted aphid proteins and their homologues among 13 other fully-sequenced arthropods and three out-group species. Subsequent analyses have revealed multiple gene expansions that are specific to aphids and have served to transfer functional annotations to 4058 pea aphid genes that display one-to-one orthology relationships with Drosophila melanogaster annotated genes. All phylogenies and alignments are accessible through the PhylomeDB database. Here we provide a description of this dataset and provide some examples on how can it be exploited.

A genetic background with low mutational robustness is associated with increased adaptability to a novel host in an RNA virus

Although mutational robustness is central to many evolutionary processes, its relationship to evolvability remains poorly understood and has been very rarely tested experimentally. Here, we measure the evolvability of Vesicular stomatitis virus in two genetic backgrounds with different levels of mutational robustness. We passaged the viruses into a novel cell type to model a host-jump episode, quantified changes in infectivity and fitness in the new host, evaluated the cost of adaptation in the original host and analyzed the genetic basis of this adaptation. Lineages evolved from the less robust genetic background demonstrated increased adaptability, paid similar costs of adaptation to the new host and fixed approximately the same number of mutations as their more robust counterparts. Theory predicts that robustness can promote evolvability only in systems where large sets of genotypes are connected by effectively neutral mutations. We argue that this condition might not be fulfilled generally in RNA viruses.

Indications for and results of liver retransplantation

Liver retransplantation (LRT) is the only therapeutic option for the irreversible failure of a hepatic graft. The aim of this study was to evaluate our rate, indications, postoperative morbidity and mortality and patient survival at one and 5 years after LRT. 1,260 liver transplants (LT) were performed between 1991 and 2006, 79 were LRT (6.3%). During the first LT, there were no apparent differences between patients who did or did not require LRT. The most common reasons for LRT were hepatic artery thrombosis (31.6%), recurrence of hepatitis C virus cirrhosis (30.4%) and primary graft non function (21.5%). The actuarial survival rates at one and 5 years were 83% and 69% among those without LRT versus 71% and 61% among those with early LRT, and 64% and 34% among those with late LRT (p < 0.001). Although high morbidity and mortality were associated with LRT, it seems that this therapeutic option is valid for patients with early hepatic loss, although not when the graft loss is late. It becomes necessary to define the minimal acceptable results so that patients can benefit from LRT.

The evolutionary history of symbiotic associations among bacteria and their animal hosts: a model

A model to explain the evolutionary history of animal-bacteria obligatory mutualistic symbiosis is presented. Dispensability of genes and genetic isolation are key factors in the reduction process of these bacterial genomes. Major steps in such genome reductive evolution, leading towards primary endosimbiosis, and the possibility of complementation or replacement by a secondary symbiont are also indicated. Yet, we need to understand what happens at the beginning of the adaptative process towards an obligate mutualistic relationship. For this purpose, we propose to sequence the complete genome of SOPE, the primary endosymbiont of the rice weevil.

Refined analysis of genetic variability parameters in hepatitis C virus and the ability to predict antiviral treatment response

Hepatitis C virus (HCV) infects approximately 3% of the world population. The chronicity of hepatitis C seems to depend on the level of genetic variability. We have recently (Torres-Puente et al., J Viral Hepat, 2008; 15: 188) reported genetic variability estimates from a large-scale sequence analysis of 67 patients infected with HCV subtypes 1a (23 patients) and 1b (44 patients) and related them to response, or lack of, to alpha-interferon plus ribavirin treatment.. Two HCV genome regions were analysed in samples prior to antiviral therapy, one compressing the three hypervariable regions of the E2 glycoprotein and another one including the interferon sensitive determining region and the V3 domain of the NS5A protein. Haplotype and nucleotide diversity measures showed a clear tendency to higher genetic variability levels in nonresponder than in responder patients. Here, we have refined the analysis of genetic variability (haplotype and nucleotide diversity, number of haplotypes and mutations) by considering their distribution in each of the biologically meaningful subregions mentioned above, as well as in their surrounding and intervening regions. Variability levels are very heterogeneous among the different subregions, being higher for nonresponder patients. Interestingly, significant differences were detected in the biologically relevant regions, but also in the surrounding regions, suggesting that the level of variability of the whole HCV genome, rather than exclusively that from the hypervariable regions, is the main indicator of the treatment response. Finally, the number of haplotypes and mutations seem to be better discriminators than haplotype and nucleotide diversity, especially in the NS5A region.

Genetic variability in hepatitis C virus and its role in antiviral treatment response

Hepatitis C virus (HCV) is a major health problem worldwide, infecting an estimated 170 million people. The high genetic variability of HCV contributes to the chronicity of hepatitis C. Here, we report results from a large-scale sequence analysis of 67 patients infected with HCV genotype 1, 23 with subtype 1a and 44 with subtype 1b. Two regions of the HCV genome were analysed in samples prior to combined therapy with alpha interferon plus ribavirin, one compressing the hypervariable regions (HVR1, HVR2 and HVR3) of the E2 glycoprotein and another one including the interferon-sensitive determining region (ISDR) and the V3 domain of the NS5A protein. Genetic diversity measures showed a clear tendency to higher genetic variability levels in nonresponder patients to antiviral treatment than in responder patients, although highly disperse values were present within each response group for both subtypes. A more detailed analysis of amino acid composition revealed the presence of several subtype-specific variants in a few positions, but no discriminating positions between responder and nonresponder patients were detected. Our results also revealed that most amino acid positions were highly conserved, especially for subtype 1a. We conclude that the outcome of the antiviral treatment might depend not only on the nature of one or a few independent positions, but more likely on the combination of several positions along the HCV genome. Moreover, the own host's ability to generate an appropriate systemic response, in combination with the action of antivirals, is also likely to be essential for treatment outcome.

Sex versus parthenogenesis: a transcriptomic approach of photoperiod response in the model aphid Acyrthosiphon pisum (Hemiptera: Aphididae)

Most aphids develop a cyclic parthenogenesis life-cycle. After several generations of viviparous parthenogenetic females, it follows a single annual generation of sexual individuals, usually in autumn, that mate and lay the sexual eggs. Shortening of photoperiod at the end of the summer is a key factor inducing the sexual response. With the survey here reported we aimed at identifying a collection of candidate genes to participate at some point in the cascade of events that lead to the sexual phenotypes. Following a suppression subtractive hybridization methodology (SSH) on the model aphid Acyrthosiphon pisum, we built and characterised two reciprocal cDNA libraries (SDU and SDD) enriched respectively in genes up-regulated or down-regulated by short photoperiod conditions that lead to the sexual response in this aphid species. A total of 557 ESTs were obtained altogether representing 223 non-overlapping contigs. 29% of these were new sequences not present in previous aphid EST libraries. BLAST searches allowed putative identification of about 54% of the contigs present in both libraries. Relative quantification of expression through real-time quantitative PCR demonstrated the differential expression in relation with the photoperiod of 6 genes (3 up-regulated and 3 down-regulated by shortening the day length). Among these, expression of a tubulin gene, two cuticular proteins and a yet unidentified sequence along the day-night cycle was further investigated. Implications for current studies on gene regulation of the dichotomy sex vs. parthenogenesis in aphids are discussed.

The Gypsy Database (GyDB) of mobile genetic elements

In this article, we introduce the Gypsy Database (GyDB) of mobile genetic elements, an in-progress database devoted to the non-redundant analysis and evolutionary-based classification of mobile genetic elements. In this first version, we contemplate eukaryotic Ty3/Gypsy and Retroviridae long terminal repeats (LTR) retroelements. Phylogenetic analyses based on the gag-pro-pol internal region commonly presented by these two groups strongly support a certain number of previously described Ty3/Gypsy lineages originally reported from reverse-transcriptase (RT) analyses. Vertebrate retroviruses (Retroviridae) are also constituted in several monophyletic groups consistent with genera proposed by the ICTV nomenclature, as well as with the current tendency to classify both endogenous and exogenous retroviruses by three major classes (I, II and III). Our inference indicates that all protein domains codified by the gag-pro-pol internal region of these two groups agree in a collective presentation of a particular evolutionary history, which may be used as a main criterion to differentiate their molecular diversity in a comprehensive collection of phylogenies and non-redundant molecular profiles useful in the identification of new Ty3/Gypsy and Retroviridae species. The GyDB project is available at http://gydb.uv.es.

A halocin acting on Na+/H+ exchanger of haloarchaea as a new type of inhibitor in NHE of mammals

The capability of halocin H6 (a bacteriocin-like protein produced by haloarchaea Haloferax gibbonsii) to inhibit Na+/H+ exchanger (NHE) in mammalian cells and its cardio-protective efficacy on the ischemic and reperfused myocardium were evaluated in the present study. H6 inhibits NHE activity (measured by a flow cytometry method) in a dose-dependent form of cell lines of mammalian origin (HEK293, NIH3T3, Jurkat and HL-1) as well as in primary cell culture from human skeletal muscle (myocytes and fibroblasts). In vivo, an ischemia-reperfusion model in dogs by coronary arterial occlusion was used (two hours of regional ischemia and three hours of reperfusion). In animals treated with halocin H6 there was a significant reduction of premature ventricular ectopic beats and infarct size, whereas blood pressure and heart rate remained unchanged. Up to date, halocin H6 is the only described biological molecule that exerts a specific inhibitory activity in NHE of eukaryotic cells.

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.

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.

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.

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.