Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Geographic variation and positive selection on M7 lysin, an acrosomal sperm protein in mussels (Mytilus spp.)

Successful fertilization in free-spawning marine organisms depends on the interactions between genes expressed on the surfaces of eggs and sperm. Positive selection frequently characterizes the molecular evolution of such genes, raising the possibility that some common deterministic process drives the evolution of gamete recognition genes and may even be important for understanding the evolution of prezygotic isolation and speciation in the marine realm. One hypothesis is that gamete recognition genes are subject to selection for prezygotic isolation, namely reinforcement. In a previous study, positive selection on the gene coding for the acrosomal sperm protein M7 lysin was demonstrated among allopatric populations of mussels in the Mytilus edulis species group (M. edulis, M. galloprovincialis, and M. trossulus). Here, we expand sampling to include M7 lysin haplotypes from populations where mussel species are sympatric and hybridize to determine whether there is a pattern of reproductive character displacement, which would be consistent with reinforcement driving selection on this gene. We do not detect a strong pattern of reproductive character displacement; there are no unique haplotypes in sympatry nor is there consistently greater population structure in comparisons involving sympatric populations. One distinct group of haplotypes, however, is strongly affected by natural selection and this group of haplotypes is found within M. galloprovincialis populations throughout the Northern Hemisphere concurrent with haplotypes common to M. galloprovincialis and M. edulis. We suggest that balancing selection, perhaps resulting from sexual conflicts between sperm and eggs, maintains old allelic diversity within M. galloprovincialis.

Positive selection on an acrosomal sperm protein, M7 lysin, in three species of the mussel genus Mytilus

Marine invertebrate sperm proteins are particularly interesting because they are characterized by positive selection and are likely to be involved in prezyogotic isolation and, thus, speciation. Here, we present the first survey of inter and intraspecific variation of a bivalve sperm protein among a group of species that regularly hybridize in nature. M7 lysin is found in sperm acrosomes of mussels and dissolves the egg vitelline coat, permitting fertilization. We sequenced multiple alleles of the mature protein-coding region of M7 lysin from allopatric populations of mussels in the Mytilus edulis species group (M. edulis, M. galloprovincialis, and M. trossulus). A significant McDonald-Kreitman test showed an excess of fixed amino acid replacing substitutions between species, consistent with positive selection. In addition, Kolmogorov-Smirnov tests showed significant heterogeneity in polymorphism to divergence ratios for both synonymous variation and combined synonymous and non-synonymous variation within M. galloprovincialis. These results indicate that there has been adaptive evolution at M7 lysin and, furthermore, shows that positive selection on sperm proteins can occur even when post-zygotic reproductive isolation is incomplete.

Life-history consequences of divergent selection on egg size in Drosophila melanogaster

Life histories are generally assumed to evolve via antagonistic pleiotropy (negative genetic correlations) among traits, and trade-offs between life-history traits are typically studied using either phenotypic manipulations or selection experiments. We investigated the trade-off between egg size and fecundity in Drosophila melanogaster by examining both the phenotypic and genetic relationships between these traits after artificial selection for large and small eggs, relative to female body size. Egg size responded strongly to selection in both directions, increasing in the large-egg selected lines and decreasing in the small-egg selected lines. Phenotypic correlations between egg size and fecundity in the large-egg selected lines were negative, but no relationship between these traits occurred in either the control or small-egg selected lines. There was no negative genetic correlation between egg size and fecundity. Total reproductive allocation decreased in the small-egg selected lines but did not increase in the large-egg lines. Our results have three implications. First, our selection procedure may have forced females selected for large eggs into a physiological trade-off not reflected in a negative genetic correlation between these traits. Second, the lack of a negative genetic correlation between egg size and number suggests that the phenotypic trade-off frequently observed between egg size and number in other organisms may not evolve over the short term via a direct genetic trade-off whereby increases in egg size are automatically accompanied by decreased fecundity. Finally, total reproductive allocation may not evolve independently of egg size as commonly assumed.

Microevolution in island forms: the roles of drift and directional selection in morphological divergence of a passerine bird

Theory predicts that in small isolated populations random genetic drift can lead to phenotypic divergence; however this prediction has rarely been tested quantitatively in natural populations. Here we utilize natural repeated island colonization events by members of the avian species complex, Zosterops lateralis, to assess whether or not genetic drift alone is an adequate explanation for the observed patterns of microevolutionary divergence in morphology. Morphological and molecular genetic characteristics of island and mainland populations are compared to test three predictions of drift theory: (1) that the pattern of morphological change is idiosyncratic to each island; (2) that there is concordance between morphological and neutral genetic shifts across island populations; and (3) for populations whose time of colonization is known, that the rate of morphological change is sufficiently slow to be accounted for solely by genetic drift. Our results are not consistent with these predictions. First, the direction of size shifts was consistently towards larger size, suggesting the action of a nonrandom process. Second, patterns of morphological divergence among recently colonized populations showed little concordance with divergence in neutral genetic characters. Third, rate tests of morphological change showed that effective population sizes were not small enough for random processes alone to account for the magnitude of microevolutionary change. Altogether, these three lines of evidence suggest that drift alone is not an adequate explanation of morphological differentiation in recently colonized island Zosterops and therefore we suggest that the observed microevolutionary changes are largely a result of directional natural selection.

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.

Species richness in agamid lizards: chance, body size, sexual selection or ecology?

Why does species richness vary so greatly across lineages? Traditionally, variation in species richness has been attributed to deterministic processes, although it is equally plausible that it may result from purely stochastic processes. We show that, based on the best available phylogenetic hypothesis, the pattern of cladogenesis among agamid lizards is not consistent with a random model, with some lineages having more species, and others fewer species, than expected by chance. We then use phylogenetic comparative methods to test six types of deterministic explanation for variation in species richness: body size, life history, sexual selection, ecological generalism, range size and latitude. Of eight variables we tested, only sexual size dimorphism and sexual dichromatism predicted species richness. Increases in species richness are associated with increases in sexual dichromatism but reductions in sexual size dimorphism. Consistent with recent comparative studies, we find no evidence that species richness is associated with small body size or high fecundity. Equally, we find no evidence that species richness covaries with ecological generalism, latitude or range size.

Signal trait sexual dimorphism and mutual sexual selection in Drosophila serrata

The evolution of sexual dimorphism may occur when natural and sexual selection result in different optimum trait values for males and females. Perhaps the most prominent examples of sexual dimorphism occur in sexually selected traits, for which males usually display exaggerated trait levels, while females may show reduced expression of the trait. In some species, females also exhibit secondary sexual traits that may either be a consequence of a correlated response to sexual selection on males or direct sexual selection for female secondary sexual traits. In this experiment, we simultaneously measure the intersex genetic correlations and the relative strength of sexual selection on males and females for a set of cuticular hydrocarbons in Drosophila serrata. There was significant directional sexual selection on both male and female cuticular hydrocarbons: the strength of sexual selection did not differ among the sexes but males and females preferred different cuticular hydrocarbons. In contrast with many previous studies of sexual dimorphism, intersex genetic correlations were low. The evolution of sexual dimorphism in D. serrata appears to have been achieved by sex-limited expression of traits controlled by genes on the X chromosome and is likely to be in its final stages.

Agronomic performance of 27 cocoa progenies and plant selection based on productivity, self-compatibility and disease resistance

The objective of the present work was to evaluate 27 progenies of cocoa crosses considering the agronomic traits and select F1 plants within superior crosses. The experiment was installed in March 2005, in the Experimental Station Joaquim Bahiana (ESJOB), in Itajuipe, Bahia. The area of the experiment is of approximately 3 ha, with a total of 3240 plants. Thirteen evaluations of vegetative brooms, five of cushion brooms and 15 of number of pods per plant were accomplished. Thirty pollinations were made for each selected plant to test for self-compatibility. The production, based on the number of pods per plant, and resistance to witches´ broom indicated CEPEC 94 x CCN 10, RB 39 x CCN 51 and CCN 10 x VB 1151 as superior progenies. All selections tested were self-compatible. The analyses of progenies and individual tree data, associated to visual field observations, allowed the selection of 17 plants which were included in a network of regional tests to determine the phenotypic stability.

Studies on the genetics of artemisinin resistance in malaria

A existência de estirpes de Plasmodium falciparum resistentes a multiplos fármacos é um dos problemas mais graves no controlo da malária. Novos fármacos, como a artemisinina (ART) e seus derivados são cada vez mais utilizados no tratamento da malaria e muito embora até ao momento não haja registos de fármaco-resistência estável à ART o seu surgimento seria desastroso devido á falta de alternativas. A investigação apresentada nesta tese descreve a selecção de resistência estável à ART e ao artesunato (ATN) utilizando um modelo roedor de malária, o parasita Plasmodium chabaudi chabaudi (Plasmodium chabaudi). Dois clones de Plasmodium chabaudi diferentes, AS-15CQ e AS-30CQ, foram inoculados em murganhos que por sua vez foram tratados na presença de concentrações sucessivamente crescentes de ATN e ART, sendo que no final do processo de seleção de resistência, os parasitas obtidos apresentavam uma resistência de 6 e 15 vezes superior ao ATN e à ART, respectivamente, em relação aos parasitas iniciais. Os clones obtidos foram nomeados respectivamente AS-ATN (obtido a partir de AS-15CQ por seleção com pressão de ATN) e AS-ART (obtido a partir de AS-30CQ por seleção com pressão de ART). A resistência obtida durante o processo de seleção é estável após clonagem, congelamento/descongelamento, passagem sanguínea na ausência de pressão de fármaco e transmissão natural através do mosquito vector. A sequência nucleotídica e o número de cópias dos genes previamente descritos na literatura como moduladores putativos de resistência à ART e seus derivados: mdr1, cg10, tctp e atp6; foi comparada entre parasitas resistentes e sensíveis, não tendo sido encontradas nenhumas alterações, quer na sequência quer no número de cópias destes genes. Posteriormente, numa tentativa de identificar os genes envolvidos na resistância à ART e ao ATN a técnica de Linkage Group Selection (LGS) foi utilizada. Para tal dois cruzamentos genéticos foram realizados. Estes cruzamentos foram realizados entre os clones fármaco-resistentes; AS-ART e AS-ATN e um clone geneticamente distinto dos anteriores e sensível aos fármacos em estudos, AJ. Após realização do LGS quatro loci genéticos; nos cromossomas de P. chabaudi 1, 2, 6 e 8 foram encontrados associados à resistência. Atendendo a que, a selecção no cromossoma 2 era a mais forte, este locus foi submetido a subsequentes análises genéticas, tendo sido encontradas duas mutações diferentes (V739F e V770F) num gene que codifica para um enzima de desubiquitinação (gene ubp-1).

Studies on the genetics of artemisinin resistance in malaria

A existência de estirpes de Plasmodium falciparum resistentes a multiplos fármacos é um dos problemas mais graves no controlo da malária. Novos fármacos, como a artemisinina (ART) e seus derivados são cada vez mais utilizados no tratamento da malaria e muito embora até ao momento não haja registos de fármaco-resistência estável à ART o seu surgimento seria desastroso devido á falta de alternativas. A investigação apresentada nesta tese descreve a selecção de resistência estável à ART e ao artesunato (ATN) utilizando um modelo roedor de malária, o parasita Plasmodium chabaudi chabaudi (Plasmodium chabaudi). Dois clones de Plasmodium chabaudi diferentes, AS-15CQ e AS-30CQ, foram inoculados em murganhos que por sua vez foram tratados na presença de concentrações sucessivamente crescentes de ATN e ART, sendo que no final do processo de seleção de resistência, os parasitas obtidos apresentavam uma resistência de 6 e 15 vezes superior ao ATN e à ART, respectivamente, em relação aos parasitas iniciais. Os clones obtidos foram nomeados respectivamente AS-ATN (obtido a partir de AS-15CQ por seleção com pressão de ATN) e AS-ART (obtido a partir de AS-30CQ por seleção com pressão de ART). A resistência obtida durante o processo de seleção é estável após clonagem, congelamento/descongelamento, passagem sanguínea na ausência de pressão de fármaco e transmissão natural através do mosquito vector. A sequência nucleotídica e o número de cópias dos genes previamente descritos na literatura como moduladores putativos de resistência à ART e seus derivados: mdr1, cg10, tctp e atp6; foi comparada entre parasitas resistentes e sensíveis, não tendo sido encontradas nenhumas alterações, quer na sequência quer no número de cópias destes genes. Posteriormente, numa tentativa de identificar os genes envolvidos na resistância à ART e ao ATN a técnica de Linkage Group Selection (LGS) foi utilizada. Para tal dois cruzamentos genéticos foram realizados. Estes cruzamentos foram realizados entre os clones fármaco-resistentes; AS-ART e AS-ATN e um clone geneticamente distinto dos anteriores e sensível aos fármacos em estudos, AJ. Após realização do LGS quatro loci genéticos; nos cromossomas de P. chabaudi 1, 2, 6 e 8 foram encontrados associados à resistência. Atendendo a que, a selecção no cromossoma 2 era a mais forte, este locus foi submetido a subsequentes análises genéticas, tendo sido encontradas duas mutações diferentes (V739F e V770F) num gene que codifica para um enzima de desubiquitinação (gene ubp-1).

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

Mitochondrial targeting adaptation of the hominoid-specific glutamate dehydrogenase driven by positive Darwinian selection.

Many new gene copies emerged by gene duplication in hominoids, but little is known with respect to their functional evolution. Glutamate dehydrogenase (GLUD) is an enzyme central to the glutamate and energy metabolism of the cell. In addition to the single, GLUD-encoding gene present in all mammals (GLUD1), humans and apes acquired a second GLUD gene (GLUD2) through retroduplication of GLUD1, which codes for an enzyme with unique, potentially brain-adapted properties. Here we show that whereas the GLUD1 parental protein localizes to mitochondria and the cytoplasm, GLUD2 is specifically targeted to mitochondria. Using evolutionary analysis and resurrected ancestral protein variants, we demonstrate that the enhanced mitochondrial targeting specificity of GLUD2 is due to a single positively selected glutamic acid-to-lysine substitution, which was fixed in the N-terminal mitochondrial targeting sequence (MTS) of GLUD2 soon after the duplication event in the hominoid ancestor approximately 18-25 million years ago. This MTS substitution arose in parallel with two crucial adaptive amino acid changes in the enzyme and likely contributed to the functional adaptation of GLUD2 to the glutamate metabolism of the hominoid brain and other tissues. We suggest that rapid, selectively driven subcellular adaptation, as exemplified by GLUD2, represents a common route underlying the emergence of new gene functions.