Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

Evolution of mitochondrial DNA in Drosophila subobscura

The colonization of the New World by the Palearctic species Drosophila subobscura was first detected in 1978 in South America and around 1982 in western North America. The ensuing dramatic expansion of the species, in territory as well as numbers, provides an opportunity for studying evolution in a scale rarely possible. We have used 10 restriction endonucleases to analyze the mitochondrial DNA (mtDNA) of individuals from 23 widely dispersed localities. Only two mtDNA composite morphs have been detected in the Americas. None of the two morphs has been found in Africa, and only one in the Atlantic islands; but both are widespread in Europe, which provides no clue of the precise geographic origin of the colonizers. The amount of nucleotide-substitution polymorphism detected in D. subobscura is typical for animals, but it is greater in the Old than in the New World, presumably due to the recent colonization by a limited number of colonizers. Assuming standard evolutionary rates of mtDNA base substitution, the mtDNA morphs found in D. subobscura can be traced to a single one that existed no less than one million years ago. We argue against the inference that the D. subobscura flies now living descend from only one or a few females that lived at that time. This type of inference, which we call the "Mother Eve hypothesis," has been made to conclude that the human population went through a severe constriction about 200,000 years ago, so that all living humans descend from only one or a few women who lived at that time. The Mother Eve hypothesis is fallacious.

Trade-off between r-selection and K-selection in Drosophila populations

Density-dependent genetic evolution was tested in experimental populations of Drosophila melanogaster subject for eight generations to natural selection under high (K-selection) or low (r-selection) population density regimes. The test consisted of determining at high and at low densities the per capita rate of population growth of the selected populations. At high densities, the K-selected populations showed a higher per capita rate of population growth than did r-selected populations, but the reverse was true at low densities. These results corroborate the predictions derived from formal models of density-dependent selection. However, no evidence of a trade-off in per capita rate of growth was observed in 25 populations of D. melanogaster, each homozygous for a different second chromosome sampled from a natural population.

Alloprocoptic selection: A mode of natural selection promoting polymorphism

The fitness effect of genetic variation at three loci coding for enzymes is studied in Drosophila melanogaster. The fertility of a female is determined by the female genotype as well as by the genotype of the male with which she mates. Significant interactions exist between female and male genotype, so that the fertility of a given mating combination cannot be predicted from the average fertility of the two genotypes involved. Multiple stable equilibria are possible when such interactions exist. At two loci, the fertility is greater than expected when the two mating individuals are homozygous for different alleles and smaller than expected when they are homozygous for the same allele. This mode of selection in which association of opposites increases their fitness is herein named alloprocoptic selection. It will contribute to maintaining genetic polymorphism in nature.

Expression and diversity of Echinococcus multilocularis AgB genes in secondarily infected mice: evaluating the influence of T-cell immune selection on antigenic variation

The T-cell-mediated immune response exhibits a crucial function in the control of the intrahepatic proliferation of Echinococcus multilocularis larvae in mice and humans, both being natural intermediate hosts of the parasite. Antigen B (AgB), a metabolized Echinococcus spp. lipoprotein, contributes to the modulation of the T-cell immune response, and distinct sites of the corresponding AgB1, AgB3 and AgB4 genes were shown to be under positive selection pressure. Since several AgB gene variants are present in a single Echinococcus metacestode, we used secondary E. multilocularis infections in BALB/c and in athymic nude mice (devoid of T-cell responses) to analyze the effect of the cellular immune response on the expression and diversity of EmAgB1-EmAgB4 genes. We demonstrated hereby that EmAgB transcripts were less abundant in nude mice during the early phase of infection (at one month post-infection), and that EmAgB2 is simultaneously down-regulated when compared to the other three genes. A negative relationship exists between the level of transcription and diversity of EmAgB genes. Moreover, no excess of non-synonymous substitutions was found among the distinct EmAgB alleles from a single host. Together, these results pointed to the effect of purifying selection, which seemed to eliminate the detrimental AgB variants generated during the development of the metacestode within the peritoneal cavity of its intermediate host.

Redundancy and recombination in the Echinococcus AgB multigene family: is there any similarity with protozoan contingency genes?

Numerous genetic variants of the Echinococcus antigen B (AgB) are encountered within a single metacestode. This could be a reflection of gene redundancy or the result of a somatic hypermutation process. We evaluate the complexity of the AgB multigene family by characterizing the upstream promoter regions of the 4 already known genes (EgAgB1-EgAgB4) and evaluating their redundancy in the genome of 3 Echinococcus species (E. granulosus, E. ortleppi and E. multilocularis) using PCR-based approaches. We have ascertained that the number of AgB gene copies is quite variable, both within and between species. The most repetitive gene seems to be AgB3, of which there are more than 110 copies in E. ortleppi. For E. granulosus, we have cloned and characterized 10 distinct upstream promoter regions of AgB3 from a single metacestode. Our sequences suggest that AgB1 and AgB3 are involved in gene conversion. These results are discussed in light of the role of gene redundancy and recombination in parasite evasion mechanisms of host immunity, which at present are known for protozoan organisms, but virtually unknown for multicellular parasites.

Searching for antigen B genes and their adaptive sites in distinct strains and species of the helminth Echinococcus

Twenty-seven PCR-derived antigen B (AgB) nucleotide sequences from four Echinococcus species (Echinococcus granulosus, Echinococcus multilocularis, Echinococcus oligarthrus and Echinococcus vogeli) were aligned with 78 already published sequences, to generate a maximum likelihood phylogeny of the AgB multigene family. The phylogenetic analysis confirms that the family is constituted by four groups of genes present in each one of the four species (AgB1, AgB2, AgB3 and AgB4), and suggests that it originated by ancient duplication events preceding speciation within the genus. AgB5 sequences, which had been formerly suggested to correspond to a putatively new AgB subunit, cluster with AgB3. Likelihood tests suggest that AgB gene evolution may have been driven by heterogeneous selection pressures acting on particular AgB1, AgB3 and AgB4 codons. No selection is detected in AgB2. We discuss implications of our findings in terms of AgB biology and its use as a diagnostic tool.

"Clonal" population structure of the malaria agent Plasmodium falciparum in high-infection regions

The population genetic structure of Plasmodium falciparum, the agent of malignant malaria, has been shown to be predominantly "clonal" (i.e., highly inbred) in regions of low infectivity; in high-infectivity regions, it is often thought to be panmictic, or nearly so, although there is little supporting evidence for this. The matter can be settled by investigating the parasite's genetic makeup in the midgut oocysts of the mosquito vector, where the products of meiosis can directly be observed. The developmental stages of P. falciparum are haploid, except in the oocysts of infected mosquito vectors, where two gametes fuse, diploidy occurs, and meiosis ensues. We have investigated genetic polymorphisms at seven microsatellite loci located on five chromosomes by assaying 613 oocysts in 145 mosquitoes sampled from 11 localities of Kenya, where malignant malaria is perennial and intense. There is considerable allelic variation, 16.3 +/- 2.1 alleles per locus, and considerable inbreeding, approximately 50% on the average. The inbreeding is caused by selfing (approximately 25%) and nonrandom genotype distribution of oocysts among mosquito guts (35%). The observed frequency of heterozygotes is 0.43 +/- 0.03; the expected frequency, assuming random mating, is 0.80 +/- 0.05. Linkage disequilibrium is statistically significant for all 21 pairwise comparisons between loci, even though 19 comparisons are between loci in different chromosomes, which is consistent with strong deviation from panmixia and the consequent reproduction of genomes as clones, without recombination between gene loci. This is of considerable evolutionary significance and of epidemiological consequence, concerning the spread of multilocus drug and vaccine resistance.

Towards a population genetics of microorganisms: The clonal theory of parasitic protozoa

Over the past 15 years, molecular investigations, including the study of isozymes and DNA markers, have provided much information on the genetic variation, population structure, breeding system and other population characteristics of parasitic protozoa. For some parasitic protozoa, but not for others, the evidence indicates that their reproduction is prevailingly clonal. In this article, Michel Tibayrenc and Francisco Ayala propose that the issue of whether the predominant mode of reproduction of a given micro-organism is clonal or sexual can only be settled by population genetics information, and they summarize evidence favoring a clonal population structure for a number of parasitic protozoa.

Meager genetic variability of the human malaria agent Plasmodium vivax

Malaria is a major human parasitic disease caused by four species of Plasmodium protozoa. Plasmodium vivax, the most widespread, affects millions of people across Africa, Asia, the Middle East, and Central and South America. We have studied the genetic variability of 13 microsatellite loci in 108 samples from 8 localities in Asia, Africa, South America, and New Guinea. Only one locus is polymorphic; nine are completely monomorphic, and the remaining three are monomorphic in all but one or two populations, which have a rare second allele. In contrast, Plasmodium falciparum displays extensive microsatellite polymorphism within and among populations. We further have analyzed, in 96 samples from the same 8 localities, 8 tandem repeats (TRs) located on a 100-kb contiguous chromosome segment described as highly polymorphic. Each locus exhibits 2-10 alleles in the whole sample but little intrapopulation polymorphism (1-5 alleles with a prevailing allele in most cases). Eight microsatellite loci monomorphic in P. vivax are polymorphic in three of five Plasmodium species related to P. vivax (two to seven individuals sampled). Plasmodium simium, a parasite of New World monkeys, is genetically indistinguishable from P. vivax. At 13 microsatellite loci and at 7 of the 8 TRs, both species share the same (or most common) allele. Scarce microsatellite polymorphism may reflect selective sweeps or population bottlenecks in recent evolutionary history of P. vivax; the differential variability of the TRs may reflect selective processes acting on particular regions of the genome. We infer that the world expansion of P. vivax as a human parasite occurred recently, perhaps <10,000 years ago.

Contingent, non-neutral evolution in a multicellular parasite: natural selection and gene conversion in the Echinococcus granulosus antigen B gene family

Recent studies have demonstrated that the Echinococcus granulosus antigen B (AgB) interferes with the intermediate hosts' immune response and is encoded by a multigene family. The number of members within the family is still uncertain, but there are several evidences of a large genetic variability. The E. granulosus AgB genomic sequences available in nucleotide databases can be grouped into four clades, corresponding to genes EgAgB1, EgAgB2, EgAgB3 and EgAgB4. In the present study, we use PCR amplifications followed by cloning and sequencing to evaluate the genetic variability for AgB isoforms. Two pairs of primers were independently used for PCR amplification. Both PCR reactions from each of three isolated protoscolex (larvae) were cloned in a plasmid vector and the plasmid inserts of 30 colonies from each cloning experiment were sequenced. Using phylogenetic tools, the 113 EgAgB clones are classified as follows: 25 are related to EgAgB1, 24 to EgAgB2, 9 to EgAgB3 and 39 to EgAgB4. The remaining 16 clones form a separate cluster, which we name EgAgB5, more closely related to EgAgB3 than to any of the other genes. Within each gene group, a number of variant sequences occur, which differ from one another by one or few nucleotides. One EgAgB3 clone has a premature stop codon (pseudogene) and an EgAgB2 clone lacks the region corresponding to the intron. The overall variation cannot be explained by differences among the asexual protoscoleces, or by experimental artifacts. Using Echinococcuss AgB genes from other species/strains as outgroups, neutrality is rejected for EgAgB2, and balancing selection is detected for EgAgB5, which also seems to be involved in gene conversion. We suggest that EgAgB1-EgAgB5 represent a family of contingency genes, that is, genes that are variably expressed, so that some but not others are expressed in each individual parasite. Contingency genes are common in parasitic protozoa and other microparasites, but the EgAgB family is the first set identified in a multicellular parasite.

The Echinococcus granulosus antigen B shows a high degree of genetic variability

Echinococcus granulosus larvae secret a polymeric lipoprotein known as antigen B (AgB) into the metacestode hydatid fluid. Three similar AgB subunits have been previously identified (AgB1, AgB2, and AgB3), and their respective genes isolated, but the actual number of genes encoding AgB subunits remains uncertain. In this study, we characterize the variability of genes encoding the AgB2 subunit, using PCR and RT-PCR followed by cloning and sequencing. We have analyzed 32 cDNA and 34 genomic sequences from a single metacestode, showing a high degree of sequence polymorphism. In addition, we have identified a possibly new AgB subunit, which we call AgB4. Additionally, we describe an AgB2 genomic clone lacking (i) a segment corresponding to the intron and (ii) a short, 45 bp sequence within exon II. The 45 bp segment encompasses the conserved splicing signals and corresponds to a highly conserved insect promoter motif.

[Pseudogenes: structure conservation, expression, and functions]

Pseudogenes have been defined as nonfunctional sequences of genomic DNA (junk DNA) originally derived from functional genes. It is therefore assumed that pseudogenes are not subject to natural selection and consequently pseudogene mutations are selectively neutral and have equal probability to become fixed in the population. We describe some unexpected features of pseudogenes in diverse organisms that are inconsistent with this widely accepted point of view. Pseudogenes are often evolutionary conserved and transcriptionally active. Moreover, pseudogenes that have been suitably investigated often exhibit functional roles, such as gene regulation, generation of genetic diversity, and other features that are expected in genes or DNA sequences that have functional roles. A review of the evidence leads to the conclusion that pseudogenes are important components of genomes, representing a repertoire of sequences available for functional evolution and subject to non-neutral evolutionary changes. Pseudogenes might be considered as potogenes, i.e. DNA sequences with a potentiality for becoming new genes or acquire new functions. Furthermore we conjecture that some pseudogenes along with their parental sequences may constitute sets of indivisible functionally interacting entities (intergenic complexes or "intergenes"), in which all the component elements are required in order to fulfill a collective functional role.

Livestock trade history, geography, and parasite strains: the mitochondrial genetic structure of Echinococcus granulosus in Argentina

A sample of 114 isolates of Echinococcus granulosus (Cestoda: Taeniidae) collected from different host species and sites in Argentina has been sequenced for 391 bp from the mitochondrial cytochrome c oxidase subunit I gene to analyze genetic variability and population structure. Nine different haplotypes were identified, 5 of which correspond to already characterized strains. Analysis of molecular variance and nested clade analysis of the distribution of haplotypes among localities within 3 main geographic regions indicate that geographic differentiation accounts for the overall pattern of genetic variability in E. granulosus populations. Significant geographic differentiation is also present when the sheep strain alone is considered. Our results suggest that geographic patterns are not due to actual restricted gene flow between regions but are rather a consequence of past history, probably related to the time and origin of livestock introduction in Argentina.

Complex evolution of orthologous and paralogous decarboxylase genes

The decarboxylases are involved in neurotransmitter synthesis in animals, and in pathways of secondary metabolism in plants. Different decarboxylase proteins are characterized for their different substrate specificities, but are encoded by homologous genes. We study, within a maximum-likelihood framework, the evolutionary relationships among dopa decarboxylase (Ddc), histidine decarboxylase (Hdc) and alpha-methyldopa hypersensitive (amd) in animals, and tryptophan decarboxylase (Wdc) and tyrosine decarboxylase (Ydc) in plants. The evolutionary rates are heterogeneous. There are differences between paralogous genes in the same lineages: 4.13 x 10(-10) nucleotide substitutions per site per year in mammalian Ddc vs. 1.95 in Hdc; between orthologous genes in different lineages, 7.62 in dipteran Ddc vs. 4.13 in mammalian Ddc; and very large temporal variations in some lineages, from 3.7 up to 54.9 in the Drosophila Ddc lineage. Our results are inconsistent with the molecular clock hypothesis.

Xanthine dehydrogenase (XDH): episodic evolution of a "neutral" protein

We investigated the evolution of xanthine dehydrogenase (Xdh) in 34 species from the three multicellular kingdoms, including one plant, two fungi, and three animal phyla, two classes of vertebrates, four orders of mammals, and two orders of insects. We adopted a model-based maximum-likelihood framework of inference. After accounting for among-site rate variation and heterogeneous nucleotide composition of the sequences using the discrete gamma distribution, and using nonhomogeneous nonstationary representations of the substitution process, the rate of amino acid replacement is 30.4 x 10(-10)/site/year when Drosophila species are compared but only approximately 18 x 10(-10)/site/year when comparisons are made between mammal orders, between insect orders, or between different animal phyla and approximately 11 x 10(-10)/site/year when comparisons are made between birds and mammals, between fungi, or between the three multicellular kingdoms. To account for these observations, the rate of amino acid replacement must have been eight or more times higher in some lineages and at some times than in others. Spastic evolution of Xdh appears to be related to the particularities of the genomes in which the locus is embedded.

Erratic overdispersion of three molecular clocks: GPDH, SOD, and XDH

The neutrality theory predicts that the rate of neutral molecular evolution is constant over time, and thus that there is a molecular clock for timing evolutionary events. It has been observed that the variance of the rate of evolution is generally larger than expected according to the neutrality theory, which has raised the question of how reliable the molecular clock is or, indeed, whether there is a molecular clock at all. We have carried out an extensive investigation of three proteins, glycerol-3-phosphate dehydrogenase (GPDH), superoxide dismutase (SOD), and xanthine dehydrogenase (XDH). We have observed that (i) the three proteins evolve erratically through time and across lineages and (ii) the erratic patterns of acceleration and deceleration differ from locus to locus, so that one locus may evolve faster in one than another lineage, whereas the opposite may be the case for another locus. The observations are inconsistent with the predictions made by various subsidiary hypotheses proposed to account for the overdispersion of the molecular clock.

Shared nucleotide composition biases among species and their impact on phylogenetic reconstructions of the Drosophilidae

Compositional changes are a major feature of genome evolution. Overlooking nucleotide composition differences among sequences can seriously mislead phylogenetic reconstructions. Large compositional variation exists among the members of the family Drosophilidae. Until now, however, base composition differences have been largely neglected in the formulations of the nucleotide substitution process used to reconstruct the phylogeny of this important group of species. The present study adopts a maximum-likelihood framework of phylogenetic inference in order to analyze five nuclear gene regions and shows that (1) the pattern of compositional variation in the Drosophilidae does not match the phylogeny of the species; (2) accounting for the heterogeneous GC content with Galtier and Gouy's nucleotide substitution model leads to a tree that differs in significant aspects from the tree inferred when the nucleotide composition differences are ignored, even though both phylogenetic hypotheses attain strong nodal support in the bootstrap analyses; and (3) the LogDet distance correction cannot completely overcome the distorting effects of the compositional variation that exists among the species of the Drosophilidae. Our analyses confidently place the Chymomyza genus as an outgroup closer than the genus Scaptodrosophila to the Drosophila genus and conclusively support the monophyly of the Sophophora subgenus.

Nucleotide sequences provide evidence of genetic exchange among distantly related lineages of Trypanosoma cruzi

Simple phylogenetic tests were applied to a large data set of nucleotide sequences from two nuclear genes and a region of the mitochondrial genome of Trypanosoma cruzi, the agent of Chagas' disease. Incongruent gene genealogies manifest genetic exchange among distantly related lineages of T. cruzi. Two widely distributed isoenzyme types of T. cruzi are hybrids, their genetic composition being the likely result of genetic exchange between two distantly related lineages. The data show that the reference strain for the T. cruzi genome project (CL Brener) is a hybrid. Well-supported gene genealogies show that mitochondrial and nuclear gene sequences from T. cruzi cluster, respectively, in three or four distinct clades that do not fully correspond to the two previously defined major lineages of T. cruzi. There is clear genetic differentiation among the major groups of sequences, but genetic diversity within each major group is low. We estimate that the major extant lineages of T. cruzi have diverged during the Miocene or early Pliocene (3-16 million years ago).

Evidence for clonal propagation in natural isolates of Plasmodium falciparum from Venezuela

We have analyzed 75 isolates of Plasmodium falciparum, collected in Venezuela during both the dry (November) and rainy (May-July) seasons, with a range of genetic markers including antigen genes and 14 random amplified polymorphic DNA (RAPD) primers. Thirteen P. falciparum stocks from Kenya and four other Plasmodium species are included in the analysis for comparison. Cross-hybridization shows that the 14 RAPD primers reveal 14 separate regions of the parasite's genome. The P. falciparum isolates are a monophyletic clade, significantly different from the other Plasmodium species. We identify three RAPD characters that could be useful as "tags" for rapid species identification. The Venezuelan genotypes fall into two discrete genetic subdivisions associated with either the dry or the rainy season; the isolates collected in the rainy season exhibit greater genetic diversity. There is significant linkage disequilibrium in each seasonal subsample and in the full sample. In contrast, no linkage disequilibrium is detected in the African sample. These results support the hypothesis that the population structure of P. falciparum in Venezuela, but not in Africa, is predominantly clonal. However, the impact of genetic recombination on Venezuelan P. falciparum seems higher than in parasitic species with long-term clonal evolution like Trypanosoma cruzi, the agent of Chagas' disease. The genetic structure of the Venezuelan samples is similar to that of Escherichia coli, a bacterium that propagates clonally, with occasional genetic recombination.

Genetic variation and the recent worldwide expansion of Plasmodium falciparum

Plasmodium falciparum, the agent of human malignant malaria, diverged from Plasmodium reichenowi, the chimpanzee parasite, about the time the human and chimpanzee lineages diverged from each other. The absence of synonymous nucleotide variation at ten loci indicates that the world populations of P. falciparum derive most recently from one single strain, or 'cenancestor,' which lived a few thousand years ago. Antigenic genes of P. falciparum (such as Csp, Msp-1, and Msp-2) exhibit numerous polymorphisms that have been estimated to be millions of years old. We have discovered in these antigenic genes short repetitive sequences that distort the alignment of alleles and account for the apparent old age of the polymorphisms. The processes of intragenic recombination that generate the repeats occur at rates about 10(-3) to 10(-2), several orders of magnitude greater than the typical mutational process of nucleotide substitutions. We conclude that the antigenic polymorphisms of P. falciparum are consistent with a recent expansion of the world populations of the parasite from a cenancestor that lived in tropical Africa a few thousand years ago.

Neutralism and selectionism: the molecular clock

The neutrality theory predicts that the rate of molecular evolution will be constant over time, and thus that there is a molecular clock for timing evolutionary events. It has been observed that the variance of the rate of evolution is generally larger than expected according to the neutrality theory. Several modifications of the theory have been proposed to account for the 'overdispersion' of the molecular clock, by postulating effects attributed to generation-time, population size, slightly deleterious mutations, repair mechanisms, and the like. An extensive investigation of two proteins, glycerol-3-phosphate dehydrogenase (GPDH) and superoxide dismutase (SOD), manifests that none of these modifications can simultaneously account for the disparate patterns observed in both proteins. GPDH evolves very slowly in Drosophila species, but several times faster in mammals, other animals, plants, and fungi. SOD evolves very fast in Drosophila species and also in mammals, but much more slowly in other animals and still slower when plants and fungi are compared to one another, or to animals.

Evidence for a high ancestral GC content in Drosophila

Study of the nucleotide composition in Drosophila, focusing on the saltans and willistoni groups, has revealed unanticipated differences in nucleotide composition among lineages. Compositional differences are associated with an accelerated rate of nucleotide substitution in functionally less constrained regions. These observations have been set forth against the extended opinion that the pattern of point mutation has remained constant during the evolution of the genus. A crucial assumption has been that the most recent common ancestor of the subgenus Sophophora had an elevated GC content. Until now, this assumption has been supported by indirect arguments, consisting of extrapolations from closely related outgroups and limited by the robustness of mathematical descriptions concerning the extensive nucleotide composition differences among sequences. The present study seeks to test the assumption of a high ancestral GC content using realistic representations of the nucleotide substitution process to account for potential biases induced by the heterogeneous GC content of the taxa. The analysis of eight nuclear genes unambiguously corroborates that the common ancestor of Sophophora had an elevated GC content.

The origin of antigenic diversity in Plasmodium falciparum

Most studies of genetic variability of Plasmodium falciparum have focused on protein antigens and the genes that encode them. The consensus is that populations exhibit high levels of genetic polymorphism, most notably the genes encoding surface proteins of the merozoite (Msp1, Msp2) and the sporozoite (Csp). The age and derivation of this variation is a subject that warrants further careful consideration, as discussed here by Stephen Rich, Marcelo Ferreira and Francisco Ayala.

Tree rooting with outgroups when they differ in their nucleotide composition from the ingroup: the Drosophila saltans and willistoni groups, a case study

Rooting is frequently the most precarious step in any phylogenetic analysis. Outgroups can become useless for rooting if they are too distantly related to the ingroup. Specifically, little attention has been paid to scenarios where outgroups have evolved different nucleotide frequencies from the ingroup. We investigate one empirical example that arose seeking to determine the phylogenetic relationship between the saltans and the willistoni groups of Drosophila (subgenus Sophophora). We have analyzed 2085 coding nucleotides from the xanthine dehydrogenase (Xdh) gene in 14 species, 6 from the saltans group and 8 from the willistoni group. We adopt a two-step strategy: (1) we investigate the phylogeny without outgroups, rooting the network by the midpoint method; (2) we reinvestigate the rooting of this phylogeny using predefined outgroups in both a parsimony- and a model-based maximum-likelihood framework. A satisfactory description of the substitution process along the Xdh region calls for six substitution types and substitution rate variation among codon positions. When the ingroup sequences are considered alone, the phylogeny obtained using this description corroborates the known relationships derived from anatomical criteria. Inclusion of the outgroups makes the root unstable, apparently because of differences between ingroups and outgroups in the substitution processes; these differences are better accounted for by a simplified model of evolution than by more complex, realistic descriptions of the substitution process.

Population structure and recent evolution of Plasmodium falciparum

Plasmodium falciparum is the agent of malignant malaria, one of mankind's most severe maladies. The parasite exhibits antigenic polymorphisms that have been postulated to be ancient. We have proposed that the extant world populations of P. falciparum have derived from one single parasite, a cenancestor, within the last 5, 000-50,000 years. This inference derives from the virtual or complete absence of synonymous nucleotide polymorphisms at genes not involved in immune or drug responses. Seeking to conciliate this claim with extensive antigenic polymorphism, we first note that allele substitutions or polymorphisms can arise very rapidly, even in a single generation, in large populations subject to strong natural selection. Second, new alleles can arise not only by single-nucleotide mutations, but also by duplication/deletion of short simple-repeat DNA sequences, a process several orders of magnitude faster than single-nucleotide mutation. We analyze three antigenic genes known to be extremely polymorphic: Csp, Msp-1, and Msp-2. We identify regions consisting of tandem or proximally repetitive short DNA sequences, including some previously unnoticed. We conclude that the antigenic polymorphisms are consistent with the recent origin of the world populations of P. falciparum inferred from the analysis of nonantigenic genes.

Disparate evolution of paralogous introns in the Xdh gene of Drosophila

Drosophila nuclear introns are commonly assumed to change according to a single rate of substitution, yet little is known about the evolution of these non-coding sequences. The hypothesis of a uniform substitution rate for introns seems to be at odds with recent findings that the nucleotide composition of introns varies at a scale unknown before, and that their base content variation is correlated with that of the adjacent exons. However, no direct attempt at comparing substitution rates in introns seems to have been addressed so far. We have studied the rate of nucleotide substitution over a region of the Xdh gene containing two adjacent short, constitutively spliced introns, in several species of Drosophila and related genera. The two introns differ significantly in base composition and substitution rate, with one intron evolving at least twice as fast as the other. In addition, the substitution pattern of the introns is positively associated with that of the surrounding coding regions, evidencing that the molecular evolution of these introns is impacted by the region in which they are embedded. The observed differences cannot be attributed to selection acting differently at the level of the secondary structure of the pre-mRNA. Rather, they are better accounted for by locally heterogeneous patterns of mutation.

A new method for characterizing replacement rate variation in molecular sequences. Application of the Fourier and wavelet models to Drosophila and mammalian proteins

We propose models for describing replacement rate variation in genes and proteins, in which the profile of relative replacement rates along the length of a given sequence is defined as a function of the site number. We consider here two types of functions, one derived from the cosine Fourier series, and the other from discrete wavelet transforms. The number of parameters used for characterizing the substitution rates along the sequences can be flexibly changed and in their most parameter-rich versions, both Fourier and wavelet models become equivalent to the unrestricted-rates model, in which each site of a sequence alignment evolves at a unique rate. When applied to a few real data sets, the new models appeared to fit data better than the discrete gamma model when compared with the Akaike information criterion and the likelihood-ratio test, although the parametric bootstrap version of the Cox test performed for one of the data sets indicated that the difference in likelihoods between the two models is not significant. The new models are applicable to testing biological hypotheses such as the statistical identity of rate variation profiles among homologous protein families. These models are also useful for determining regions in genes and proteins that evolve significantly faster or slower than the sequence average. We illustrate the application of the new method by analyzing human immunoglobulin and Drosophilid alcohol dehydrogenase sequences.

Fluctuating mutation bias and the evolution of base composition in Drosophila

The idea that the pattern of point mutation in Drosophila has remained constant during the evolution of the genus has recently been challenged. A study of the nucleotide composition focused on the Drosophila saltans group has evidenced unsuspected nucleotide composition differences among lineages. Compositional differences are associated with an accelerated rate of amino acid replacement in functionally less constrained regions. Here we reassess this issue from a different perspective. Adopting a maximum-likelihood estimation approach, we focus on the different predictions that mutation and selection make about the nonsynonymous-to-synonymous rate ratio. We investigate two gene regions, alcohol dehydrogenase (Adh) and xanthine dehydrogenase (Xdh), using a balanced data set that comprises representatives from the melangaster, obscura, saltans, and willistoni groups. We also consider representatives of the Hawaiian picture-winged group. These Hawaiian species are known to have experienced repeated bottlenecks and are included as a reference for comparison. Our results confirm patterns previously detected. The branch ancestral to the fast-evolving willistoni/saltans lineage, where most of the change in GC content has occurred, exhibits an excess of synonymous substitutions. The shift in mutation bias has affected the extent of the rate variation among sites in Xdh.

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.

Phylogeny of Drosophila and related genera: conflict between molecular and anatomical analyses

Drosophila species are extensively used in biological research; yet, important phylogenetic relationships within the genus and with related genera remain unresolved. The combined data for three genes (Adh, Sod, and Gpdh) statistically resolves outstanding issues. We define the genus Drosophila inclusively so as to include Scaptomyza and Zaprionus (considered distinct genera in the taxonomy of Wheeler, 1981) but excluding Scaptodrosophila. The genus Drosophila so defined is monophyletic. The subgenus Sophophora (including the melanogaster, obscura, and willistoni groups) is monophyletic and the sister clade to all other Drosophila subgenera. The Hawaiian Drosophila (including Scaptomyza) is a monophyletic group, but the subgenus Drosophila is not monophyletic, because the immigrans group is more closely related to the subgenus Hirtodrosophila than to other species of the subgenus Drosophila, such as the virilis and repleta groups.

Molecular evolution and phylogeny of the Drosophila saltans species group inferred from the Xdh gene

The Drosophila saltans group of the subgenus Sophophora consists of five species subgroups whose phylogenetic relationships are poorly known. We have analyzed 2085 coding nucleotides from the xanthine dehydrogenase (Xdh) gene in six species, at least one from each subgroup. We follow a model-based maximum likelihood framework. We first model the substitution process using a tree topology that is approximately accurate. Then we evaluate several candidate tree topologies using a working model of nucleotide substitution. We found that a minimally realistic description of the substitution process along the Xdh region should allow two transition and four transversion rate parameters and different fixed rates for codon positions, which are distributed statistically according to different gamma distributions. The phylogeny obtained using this description differs in significant respects from a phylogeny based on anatomical criteria. We have also analyzed data from five additional (three nuclear and two mitochondrial) gene regions. In our analysis, these relatively short DNA sequences, either separately or jointly, fail to discriminate statistically among alternative phylogenies. When the data for these five gene regions are combined with the Xdh sequences, the strong phylogenetic signal emerging from Xdh becomes somewhat diluted rather than reinforced. The phylogeny of the species and biogeographical considerations suggest that the D. saltans group originated in the tropics of the New World, similarly as the closely related D. willistoni group.

Fitness of wild-caught Drosophila melanogaster females: allozyme variants of GPDH, ADH, PGM, and EST

We have collected several hundred Drosophila melanogaster flies (near Davis, California), isolated them individually, without anesthesia, at the collecting site, and estimated the fitness components of the wild-caught females under different environmental conditions. The fitness parameters measured are fecundity, oviposition rate, and productivity (egg-to-adult viability, development rate, and number of progeny). The environmental variables are two temperatures (22 degrees C and 28 degrees C) and two densities ('scant' and 'crowded'). After the fitness measurements are completed for each individual female, its genotype is determined at four loci encoding enzymes: GPDH and ADH, located on chromosome II; and PGM and EST-C, located on chromosome III. Density has a large significant effect on productivity; temperature has significant effects on fecundity, oviposition rate, and development rate. The experiments show that allozyme polymorphisms are associated with selection effects. Fitness differences between allozyme genotypes occur for all fitness components, except oviposition rate. But which genotype is superior depends on the environmental conditions; heterozygotes exhibit higher fitness than homozygotes in a number of cases, but inferior in others. A unique feature of the present experiments is that the experimental flies are wild-caught females rather than laboratory-bred individuals.

Switch in codon bias and increased rates of amino acid substitution in the Drosophila saltans species group

We investigated the nucleotide composition of five genes, Xdh, Adh, Sod, Per, and 28SrRNA, in nine species of Drosophila (subgenus Sophophora) and one of Scaptodrosophila. The six species of the Drosophila saltans group markedly differ from the others in GC content and codon use bias. The GC content in the third codon position, and to a lesser extent in the first position and the introns, is higher in the D. melanogaster and D. obscura groups than in the D. saltans group (in Scaptodrosophila it is intermediate but closer to the melanogaster and obscura species). Differences are greater for Xdh than for Adh, Sod, Per, and 28SrRNA, which are functionally more constrained. We infer that rapid evolution of GC content in the saltans lineage is largely due to a shift in mutation pressure, which may have been associated with diminished natural selection due to smaller effective population numbers rather than reduced recombination rates. The rate of GC content evolution impacts the rate of protein evolution and may distort phylogenetic inferences. Previous observations suggesting that GC content evolution is very limited in Drosophila may have been distorted due to the restricted number of genes and species (mostly D. melanogaster) investigated.

Evolution of Plasmodium and the recent origin of the world populations of Plasmodium falciparum

We have investigated the evolution of Plasmodium parasites by analyzing DNA sequences of several genes. We reach the following conclusions: (1) The four human parasites, P. falciparum, P. malariae, P. ovale, and P. vivax are very remotely related to each other, so that their evolutionary divergence predates the origin of the hominids; several of these parasites became associated with the human lineage by lateral transfer from other hosts. (2) P. falciparum diverged from P. reichenowi about 8 million years ago, consistently with the time of divergence of the human lineage from the apes; a parsimonious inference is that falciparum has been associated with humans since the origin of the hominids. (3) P. malariae is genetically indistinguishable from P. brasilianum, a parasite of New World monkeys; and, similarly. (4) P. vivax is genetically indistinguishable from the New World monkey parasite P. simium. We infer in each of these two cases a very recent lateral transfer between the human and monkey hosts, and explore alternative hypotheses about the direction of the transfer. We have also investigated the population structure of P. falciparum by analyzing 10 genes and conclude that the extant world populations of this parasite have evolved from a single strain within the last several thousand years. The extensive polymorphisms observed in the highly repetitive central region of the Csp gene, as well as the apparently very divergent two classes of alleles at the Msa-1 gene, are consistent with this conclusion.

Plasmodium falciparum: population genetic analysis by multilocus enzyme electrophoresis and other molecular markers

Abderrazak, S. B., Oury, B, Lal, A. A., Bosseno, M.-F., Force-Barge, P., Dujardin, J.-P., Fandeur, T., Molez, J.-F., Kjellberg, F., Ayala, F. J., and Tibayrenc, M. 1999. Plasmodium falciparum: Population genetic analysis by multilocus enzyme electrophoresis and other molecular markers. Experimental Parasitology 92, 232-238. The population structure of Plasmodium falciparum, the agent of malignant malaria, is uncertain. We have analyzed multilocus enzyme electrophoresis (MLEE) polymorphisms at 7-12 gene loci in each of four populations (two populations in Burkina Faso, one in Sudan, one in Congo), plus one "cosmopolitan" sample consisting of parasite cultures from 15 distant localities in four different continents. We have also performed random amplified polymorphic DNA analysis (RAPD) and restriction fragment length polymorphism (RFLP) and characterized gene varia tion at four antigen genes in the Congo population. All genetic assays show abundant genetic variability in all populations analyzed. With the isoenzyme assays, strong linkage disequilibrium is apparent in at least two local populations, the Congo population and one population from Burkina Faso, as well as in the cosmopolitan sample, and less definitely in the other Burkina Faso population. However, no linkage disequilibrium is detected in the Congo population with the molecular assays. We failed to detect any nonrandom association between the different kinds of genetic markers; that is, MLEE with RAPD or RFLP, RAPD with RFLP, and so on. Although isoenzyme data show statistical departures from panmictic expectations, these results suggest that in the areas under survey, P. falciparum populations do not undergo predominant clonal evolution and show no clear-cut subdivisions, un like Trypanosoma cruzi, Leishmania sp., and other major parasitic species. We discuss the epidemiological and taxonomical significance of these results.