Mutational and selective pressures on codon and amino acid usage in Buchnera, endosymbiotic bacteria of aphids

We have explored compositional variation at synonymous (codon usage) and nonsynonymous (amino acid usage) positions in three complete genomes of Buchnera, endosymbiotic bacteria of aphids, and also in their orthologs in Escherichia coli, a close free-living relative. We sought to discriminate genes of variable expression levels in order to weigh the relative contributions of mutational bias and selection in the genomic changes following symbiosis. We identified clear strand asymmetries, distribution biases (putative high-expression genes were found more often on the leading strand), and a residual slight codon bias within each strand. Amino acid usage was strongly biased in putative high-expression genes, characterized by avoidance of aromatic amino acids, but above all by greater conservation and resistance to AT enrichment. Despite the almost complete loss of codon bias and heavy mutational pressure, selective forces are still strong at nonsynonymous sites of a fraction of the genome. However, Buchnera from Baizongia pistaciae appears to have suffered a stronger symbiotic syndrome than the two other species.

Molecular systematics of aphids (Homoptera: Aphididae): new insights from the long-wavelength opsin gene

Viviparous aphids (Aphididae) constitute a monophyletic group within the Homoptera with more than 4000 extant species worldwide but higher diversity in temperate regions. Several aspects of their biology account for attention paid to this group of insects. Their plant-sap-sucking way of feeding with many species transmitting viruses to crop plants has important implications on crop management strategies. Cyclical parthenogenesis associated in many groups to host alternation and elaborate polyphenisms is of special interests for evolutionists. Finally, the ancient association of most aphid species with intracellular endosymbiotic bacteria (Buchnera sp.) has also received much attention from evolutionists interested in mechanisms involved in the symbiotic process. Knowing the phylogenetic relationships among major aphid taxa is of special interest to evolutionists interested in the above issues. However, until recently, molecular approaches to aphid phylogeny were absent and discussions on the evolution of aphid life-cycles and on evolutionary aspects of their symbiotic association with Buchnera were framed by morphology-based phylogenies. Recently, two reports using molecular approaches attempted to address the yet unresolved phylogeny of Aphididae with limited although somehow different conclusions. In the present report we study the utility of the long-wave opsin gene in resolving phylogenetic relationships among seven subfamilies within the Aphididae. Our results corroborate some previously proposed relationships and suggest a revision of some others. In particular, our data support grouping the analysed aphid species into three main clades, being the subfamily Lachninae one of them, which contradicts its generally accepted sistership relationship with the subfamily Aphidinae. Moreover, our data also suggest a basal position of Lachninae which has implications on current discussions about the ancestrality of conifer-feeding in modern aphids.

Selection on coding regions determined Hox7 genes evolution

The important role of Hox genes in determining the regionalization of the body plan of the vertebrates makes them invaluable candidates for evolutionary analyses regarding functional and morphological innovation. Gene duplication and gene loss led to a variable number of Hox genes in different vertebrate lineages. The evolutionary forces determining the conservation or loss of Hox genes are poorly understood. In this study, we show that variable selective pressures acted on Hox7 genes in different evolutionary lineages, with episodes of positive selection occurring after gene duplications. Tests for functional divergence in paralogs detected significant differentiation in a region known to modulate HOX7 protein activity. Our results show that both positive and negative selection on coding regions are influencing Hox7 genes evolution.

Assortative mating and fertility in two Drosophila subobscura strains with different mitochondrial DNA haplotypes

The mating pattern and female fertility on the two main mitochondrial DNA haplotypes (I and II) of Drosophila subobscura were studied, in an attempt to find possible differences between them in relation to sexual selection or isolation that could explain the populational dynamics and the co-existence of these two strains in nature. The mating pattern indicated an assortative mating in population cages, where couples of the same haplotype, mainly those of haplotype I, mated more often. However, the significations detected in laboratory conditions disappeared in wild populations, where random mating was the rule. The female fertility also showed differences in the laboratory compared to the wild, since couples with haplotype I males were more efficient in the laboratory populations. These results, together with others that we previously obtained, either point to selection acting directly on the mtDNA or to the presence of some kind of cytonuclear co-adaptation in these two haplotypes, although this must be modulated by other factors that change with the seasons and time. The end result could well be a balance of opposite forces acting on both haplotypes.

Sampling and repeatability in the evaluation of hepatitis C virus genetic variability

Among the experimental techniques available to study the genetic variability of RNA virus populations, the most informative involve reverse transcription (RT), amplification, cloning and sequencing. The effects of several aspects of these techniques on the estimation of genetic variability in a virus population were analysed. Hepatitis C virus populations from four patients were examined. For each patient, ten series of data derived from independent PCR amplifications of a single RT reaction were obtained. The sample size of each data set was 10 sequences (in nine series) and 100 sequences (in one series). An additional data set derived from an independent RT reaction (about 10 sequences) performed on RNA extracted from the same serum sample was also analysed. The availability of data sets of different sample sizes allowed the effect of sample size on the amount and nature of the genetic variability recovered to be examined. The repeatability of the data obtained in different amplification experiments as well as from different RT reactions was also determined, together with the best strategy to obtain a given number of sequences by comparing the set of 100 sequences obtained from a single amplification with those obtained by pooling the nine sets of 10 sequences. In all cases, these results confirm the high repeatability of the conclusions and parameters derived from the sets of 10 sequences. These results validate the use of relatively small sample sets for the evaluation of genetic variability and for the estimation of phylogenetic relationships of RNA viruses in population and epidemiological studies.

The genome sequence of Blochmannia floridanus: comparative analysis of reduced genomes

Bacterial symbioses are widespread among insects, probably being one of the key factors of their evolutionary success. We present the complete genome sequence of Blochmannia floridanus, the primary endosymbiont of carpenter ants. Although these ants feed on a complex diet, this symbiosis very likely has a nutritional basis: Blochmannia is able to supply nitrogen and sulfur compounds to the host while it takes advantage of the host metabolic machinery. Remarkably, these bacteria lack all known genes involved in replication initiation (dnaA, priA, and recA). The phylogenetic analysis of a set of conserved protein-coding genes shows that Bl. floridanus is phylogenetically related to Buchnera aphidicola and Wigglesworthia glossinidia, the other endosymbiotic bacteria whose complete genomes have been sequenced so far. Comparative analysis of the five known genomes from insect endosymbiotic bacteria reveals they share only 313 genes, a number that may be close to the minimum gene set necessary to sustain endosymbiotic life.

Evolution of RNA virus in spatially structured heterogeneous environments

A hallmark of the infectious cycle for many RNA viruses parasitizing multicellular hosts is the need to invade and successfully replicate in tissues that comprise a variety of cell types. Thus, multicellular hosts represent a heterogeneous environment to evolving viral populations. To understand viral adaptation to multicellular hosts, we took a double approach. First, we developed a mathematical model that served to make predictions concerning the dynamics of viral populations evolving in heterogeneous environments. Second, the predictions were tested by evolving vesicular stomatitis virus in vitro on a spatially structured environment formed by three different cell types. In the absence of gene flow, adaptation was tissue-specific, but fitness in all tissues decreased with migration rate. The performance in a given tissue was negatively correlated with its distance to the tissue hosting the population. This correlation decreased with migration rate.

Why are the genomes of endosymbiotic bacteria so stable?

The comparative analysis of three strains of the endosymbiotic bacterium Buchnera aphidicola has revealed high genome stability associated with an almost complete absence of chromosomal rearrangements and horizontal gene transfer events during the past 150 million years. The loss of genes involved in DNA uptake and recombination in the initial stages of endosymbiosis probably underlies this stability. Gene loss, which was extensive during the initial steps of Buchnera evolution, has continued in the different Buchnera lineages since their divergence.

Side-stepping secondary symbionts: widespread horizontal transfer across and beyond the Aphidoidea

To elucidate the co-evolutionary relationships between phloem-feeding insects and their secondary, or facultative, bacterial symbionts, we explore the distributions of three such microbes--provisionally named the R-type (or PASS, or S-sym), T-type (or PABS), and U-type--across a number of aphid and psyllid hosts through the use of diagnostic molecular screening techniques and DNA sequencing. Although typically maternally transmitted, phylogenetic and pairwise divergence analyses reveal that these bacteria have been independently acquired by a variety of unrelated insect hosts, indicating that horizontal transfer has helped to shape their distributions. Based on the high genetic similarity between symbionts in different hosts, we argue that transfer events have occurred recently on an evolutionary timescale. In several instances, however, closely related symbionts associate with related hosts, suggesting that horizontal transfer between distant relatives may be rarer than transmission between close relatives. Our findings on the prevalence of these symbionts within many aphid taxa, along with published observations concerning their effects on host fitness, imply a significant role of facultative symbiosis in aphid ecology and evolution.

Identification of a gene overexpressed in aphids reared under short photoperiod

Most aphids develop a cyclic parthenogenesis life-cycle. After several generations of viviparously produced parthenogenetic females, 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 (together with temperature) is a key factor inducing the sexual response. Currently no genes involved in the cascade of events that lead to the appearance of sexual forms have been reported. After a Differential Display RT-PCR survey performed on Acyrthosiphon pisum aphids, we identified a gene that is overexpressed in aphids reared under short photoperiod conditions that induce sexuality in this species. This cDNA (called ApSDI-1) shows similarities with a protein involved in amino acid transport in GABAergic neurons. Since several studies implicate GABAergic transmission in the generation and modulation of circadian rhythmicity, we propose that ApSDI-1 could be involved in the transduction of the photoperiodic message and therefore be a candidate to participate at some point in processes that trigger the sexual response in aphids. This is the first gene identified in aphids whose expression is governed by the photoperiod.

Variation of haplotype distributions of two genomic regions of Citrus tristeza virus populations from eastern Spain

Genetic variation in natural populations of Citrus tristeza virus (CTV) was studied using haplotypes detected by single-strand conformation polymorphism (SSCP) analysis of two genomic regions (p20 gene and segment A, located in ORF1a). Analysis of 254 samples from 125 trees, collected at 12 different sites, yielded 8 different haplotypes for p20 and 5 for segment A. The most frequent haplotype of p20 was predominant at all sites, but several sites differed in the predominance of segment A haplotypes. At most sites, the homozygosity observed for the p20 gene tended to be higher than expected in a neutral evolution, whereas the opposite was true for segment A. Comparison of the populations at different sites showed that 44 of the 66 possible population pairs were genetically distinct for segment A, but only six pairs differed for the p20 gene. Analysis of molecular variance grouping trees by site, scion variety, rootstock or age, showed that variation in segment A was significantly affected by site, tree age and rootstock, and that variation between trees in each group and within trees was even more important. In contrast, variation in p20 was affected only by site and rootstock, each factor contributing to < 2% of the variation. The data suggest that sequence variations in segment A must be functionally less important and that it has less evolutionary constraints than p20. Detection of different haplotypes in neighbour trees or in samples from the same tree may help explain part of the variability observed in CTV symptom expression.

Molecular epidemiology and forensic genetics: application to a hepatitis C virus transmission event at a hemodialysis unit

Molecular phylogenetic analyses are frequently used in epidemiologic testing, although only occasionally in forensics. Their acceptability is hampered by a lack of statistical confidence in the conclusions. However, maximum likelihood testing provides a sound statistical framework for the testing of phylogenetic hypotheses relevant for forensic analysis. We present the results of applying this method to a small hepatitis C outbreak produced in a hospital hemodialysis unit that involved 6 patients. Polymerase chain reaction products from a 472-nt fragment of the E1-E2 region, including the hypervariable region, HVR-1, of the hepatitis C virus genome were cloned, and an average of 10 clones/patient and from 11 additional control patients were sequenced. The method allows a statistical evaluation that the likelihood of each sample belonging or not to a given group, a question of relevance in many forensic and epidemiological analyses of molecular sequences.

Reductive genome evolution in Buchnera aphidicola

We have sequenced the genome of the intracellular symbiont Buchnera aphidicola from the aphid Baizongia pistacea. This strain diverged 80-150 million years ago from the common ancestor of two previously sequenced Buchnera strains. Here, a field-collected, nonclonal sample of insects was used as source material for laboratory procedures. As a consequence, the genome assembly unveiled intrapopulational variation, consisting of approximately 1,200 polymorphic sites. Comparison of the 618-kb (kbp) genome with the two other Buchnera genomes revealed a nearly perfect gene-order conservation, indicating that the onset of genomic stasis coincided closely with establishment of the symbiosis with aphids, approximately 200 million years ago. Extensive genome reduction also predates the synchronous diversification of Buchnera and its host; but, at a slower rate, gene loss continues among the extant lineages. A computational study of protein folding predicts that proteins in Buchnera, as well as proteins of other intracellular bacteria, are generally characterized by smaller folding efficiency compared with proteins of free living bacteria. These and other degenerative genomic features are discussed in light of compensatory processes and theoretical predictions on the long-term evolutionary fate of symbionts like Buchnera.

A sliding window-based method to detect selective constraints in protein-coding genes and its application to RNA viruses

Here we present a new sliding window-based method specially designed to detect selective constraints in specific regions of a multiple protein-coding sequence alignment. In contrast to previous window-based procedures, our method is based on a nonarbitrary statistical approach to find the appropriate codon-window size to test deviations of synonymous (d(S)) and nonsynonymous (d(N)) nucleotide substitutions from the expectation. The probabilities of d(N) and d(S) are obtained from simulated data and used to detect significant deviations of d(N) and d(S) in a specific window region of the real sequence alignment. The nonsynonymous-to-synonymous rate ratio (w = d(N)/d(S)) was used to highlight selective constraints in any window wherein d(S) or d(N) was significantly different from the expectation. In these significant windows, w and its variance [V(w)] were calculated and used to test the neutral hypothesis. Computer simulations showed that the method is accurate even for highly divergent sequences. The main advantages of the new method are that it (i) uses a statistically appropriate window size to detect different selective patterns, (ii) is computationally less intensive than maximum likelihood methods, and (iii) detects saturation of synonymous sites, which can give deviations from neutrality. Hence, it allows the analysis of highly divergent sequences and the test of different alternative hypothesis as well. The application of the method to different human immunodeficiency virus type 1 and to foot-and-mouth disease virus genes confirms the action of positive selection on previously described regions as well as on new regions.

Molecular basis of adaptive convergence in experimental populations of RNA viruses

Characterizing the molecular basis of adaptation is one of the most important goals in modern evolutionary genetics. Here, we report a full-genome sequence analysis of 21 independent populations of vesicular stomatitis ribovirus evolved on the same cell type but under different demographic regimes. Each demographic regime differed in the effective viral population size. Evolutionary convergences are widespread both at synonymous and nonsynonymous replacements as well as in an intergenic region. We also found evidence for epistasis among sites of the same and different loci. We explain convergences as the consequence of four factors: (1) environmental homogeneity that supposes an identical challenge for each population, (2) structural constraints within the genome, (3) epistatic interactions among sites that create the observed pattern of covariation, and (4) the phenomenon of clonal interference among competing genotypes carrying different beneficial mutations. Using these convergences, we have been able to estimate the fitness contribution of the identified mutations and epistatic groups. Keeping in mind statistical uncertainties, these estimates suggest that along with several beneficial mutations of major effect, many other mutations got fixed as part of a group of epistatic mutations.

Comparative molecular evolution of primary (Buchnera) and secondary symbionts of aphids based on two protein-coding genes

A+T content, phylogenetic relationships, codon usage, evolutionary rates, and ratio of synonymous versus non-synonymous substitutions have been studied in partial sequences of the atpD and aroQ/pheA genes of primary ( Buchnera) and secondary symbionts of aphids and a set of selected non-symbiotic bacteria, belonging to the five subdivisions of the Proteobacteria. Compared to the homologous genes of the last group, both genes belonging to Buchnera behave in a similar way, showing a higher A+T content, forming a monophyletic group, a loss in codon bias, especially in third base position, an evolutionary acceleration and an increase in the number of non-synonymous substitutions, confirming previous results reported elsewhere for other genes. When available, these properties have been partly observed with the secondary symbionts, but with values that are intermediate between Buchnera and free living Proteobacteria. They show high A+T content, but not as high as Buchnera, a non-solved phylogenetic position between Buchnera, and the other gamma-Proteobacteria, a loss in codon bias, again not as high as in Buchnera and a significant evolutionary acceleration in the case of the three atpD genes, but not when considering aroQ/pheA genes. These results give support to the hypothesis that they are symbionts at different stages of the symbiotic accommodation to the host.

The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection

The heat-shock protein GroEL is a double-ring-structured chaperonin that assists the folding of many newly synthesized proteins in Escherichia coli and the refolding in vitro, with the cochaperonin GroES, of conformationally damaged proteins. This protein is constitutively overexpressed in the primary symbiotic bacteria of many insects, constituting approximately 10% of the total protein in Buchnera, the primary endosymbiont of aphids. In the present study, we perform a maximum likelihood (ML) analysis to unveil the selective constraints in GroEL. In addition, we apply a new statistical approach to determine the patterns of evolution in this highly interesting protein. The main conclusion derived from our analysis is that GroEL has suffered an accelerated rate of amino acid substitution upon the symbiotic integration of Buchnera into the aphids. It is most interesting that the ML analysis of codon substitutions in the different branches of the phylogenetic tree strongly supports the action of positive selection in the different lineages of BUCHNERA: Additionally, the new sliding window analysis of the complete groEL sequence reveals different regions of the molecule under the action of positive selection, mainly located in the apical domain, that are important for both peptide and GroES binding.

Extreme genome reduction in Buchnera spp.: toward the minimal genome needed for symbiotic life

Buchnera is a mutualistic intracellular symbiont of aphids. Their association began about 200 million years ago, with host and symbiont lineages evolving in parallel since that time. During this coevolutionary process, Buchnera has experienced a dramatic decrease of genome size, retaining only essential genes for its specialized lifestyle. Previous studies reported that genome size in Buchnera spp. is very uniform, suggesting that genome shrinkage occurred early in evolution, and that modern lineages retain the genome size of a common ancestor. Our physical mapping of Buchnera genomes obtained from five aphid lineages shows that the genome size is not conserved among them, but has been reduced down to 450 kb in some species. Here we show evidence of six species with a genome size smaller than Mycoplasma genitalium, the smallest bacterial genome reported thus far (580 kb). Our findings strongly suggest that the Buchnera genome is still experiencing a reductive process toward a minimum set of genes necessary for its symbiotic lifestyle.

Frequency-dependent selection in human immunodeficiency virus type 1

Genetic variation is the main evolutionary strategy adopted by RNA viruses and retroviruses. Evolution operates through competition between different individuals in the same environment, resulting in the imposition of the fittest variant. The process of competition could be affected by various factors, including the frequency of the different competing individuals. In order to investigate this aspect, individual virus populations derived from a human immunodeficiency virus type 1 isolate were studied at different competing proportions. The dynamics of variant imposition in each competition experiment permitted the detection of frequency-dependent selection (FDS); i.e. the imposition of variants is related to their biological fitness, which is also affected by the proportions at which they compete. The existence of FDS in different viruses with RNA genomes would indicate a general mechanism favouring genetic heterogeneity.

Genome size reduction through multiple events of gene disintegration in Buchnera APS

The evolution of the endosymbiont Buchnera during its adaptation to intracellular life involved a massive reduction in its genome. By comparing the orthologous genes of Buchnera, Escherichia coli and Vibrio cholerae, we show that the minimal genome size of Buchnera arose from multiple events of gene disintegration dispersed over the whole genome. The elimination of the genes was a continuous process that began with gene inactivation and progressed until the DNA corresponding to the pseudogenes were completely deleted.

Phylogenetic analysis of viroid and viroid-like satellite RNAs from plants: a reassessment

The proposed monophyletic origin of a group of subviral plant pathogens (viroids and viroid-like satellite RNAs), as well as the phylogenetic relationships and the resulting taxonomy of these entities, has been recently questioned. The criticism comes from the (apparent) lack of sequence similarity among these RNAs necessary to reliably infer a phylogeny. Here we show that, despite their low overall sequence similarity, a sequence alignment manually adjusted to take into account all the local similarities and the insertions/deletions and duplications/rearrangements described in the literature for viroids and viroid-like satellite RNA, along with the use of an appropriate estimator of genetic distances, constitutes a data set suitable for a phylogenetic reconstruction. When the likelihood-mapping method was applied to this data set, the tree-likeness obtained was higher than that corresponding to a sequence alignment that does not take into consideration the local similarities. In addition, bootstrap analysis also supports the major groups previously proposed and the reconstruction is consistent with the biological properties of this RNAs.