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.

Genetic variability and natural selection at the ligand domain of the Duffy binding protein in brazilian Plasmodium vivax populations

Background: Plasmodium vivax malaria is a major public health challenge in Latin America, Asia and Oceania, with 130-435 million clinical cases per year worldwide. Invasion of host blood cells by P. vivax mainly depends on a type I membrane protein called Duffy binding protein (PvDBP). The erythrocyte-binding motif of PvDBP is a 170 amino-acid stretch located in its cysteine-rich region II (PvDBP(II)), which is the most variable segment of the protein. Methods: To test whether diversifying natural selection has shaped the nucleotide diversity of PvDBP(II) in Brazilian populations, this region was sequenced in 122 isolates from six different geographic areas. A Bayesian method was applied to test for the action of natural selection under a population genetic model that incorporates recombination. The analysis was integrated with a structural model of PvDBP(II), and T-and B-cell epitopes were localized on the 3-D structure. Results: The results suggest that: (i) recombination plays an important role in determining the haplotype structure of PvDBP(II), and (ii) PvDBP(II) appears to contain neutrally evolving codons as well as codons evolving under natural selection. Diversifying selection preferentially acts on sites identified as epitopes, particularly on amino acid residues 417, 419, and 424, which show strong linkage disequilibrium. Conclusions: This study shows that some polymorphisms of PvDBP(II) are present near the erythrocyte-binding domain and might serve to elude antibodies that inhibit cell invasion. Therefore, these polymorphisms should be taken into account when designing vaccines aimed at eliciting antibodies to inhibit erythrocyte invasion.

SIZE AS A LINE OF LEAST RESISTANCE II: DIRECT SELECTION ON SIZE OR CORRELATED RESPONSE DUE TO CONSTRAINTS?

Evolutionary change in New World Monkey (NWM) skulls occurred primarily along the line of least resistance defined by size (including allometric) variation (g(max)). Although the direction of evolution was aligned with this axis, it was not clear whether this macroevolutionary pattern results from the conservation of within population genetic covariance patterns (long-term constraint) or long-term selection along a size dimension, or whether both, constraints and selection, were inextricably involved. Furthermore, G-matrix stability can also be a consequence of selection, which implies that both, constraints embodied in g(max) and evolutionary changes observed on the trait averages, would be influenced by selection Here, we describe a combination of approaches that allows one to test whether any particular instance of size evolution is a correlated by-product due to constraints (g(max)) or is due to direct selection on size and apply it to NWM lineages as a case study. The approach is based on comparing the direction and amount of evolutionary change produced by two different simulated sets of net-selection gradients (beta), a size (isometric and allometric size) and a nonsize set. Using this approach it is possible to distinguish between the two hypotheses (indirect size evolution due to constraints or direct selection on size), because although both may produce an evolutionary response aligned with g(max), the amount of change produced by random selection operating through the variance/covariance patterns (constraints hypothesis) will be much smaller than that produced by selection on size (selection hypothesis). Furthermore, the alignment of simulated evolutionary changes with g(max) when selection is not on size is not as tight as when selection is actually on size, allowing a statistical test of whether a particular observed case of evolution along the line of least resistance is the result of selection along it or not. Also, with matrix diagonalization (principal components [PC]) it is possible to calculate directly the net-selection gradient on size alone (first PC [PC1]) by dividing the amount of phenotypic difference between any two populations by the amount of variation in PC1, which allows one to benchmark whether selection was on size or not

Sob o signo de Darwin? Sobre o mau uso de uma quimera

No ano de 2003 Francisco de Oliveria publicou um artigo intitulado "O Ornitorrinco" no qual fez considerações críticas sobre a conjectura politico-social daquele momento histórico. Tal artigo é permeado por um paralelo entre o evolucionismo darwinista e a visão do autor sobre a sociedade brasileira contemporânea. Entretanto, ao fazer tal analogia ele incorre numa série de equívocos teóricos sobre a teoria evolucionista. Tais equívocos consistem, em grande parte, numa substuição indevida entre aquilo que ficou conhecido como Darwinismo Social e a teoria neodarwinista como entendida pelos seus atuais proponentes. O presente trabalho identifica estes equívocos e os contextualiza dentro da teoria neodarwiniana. Além disso, fazemos um recorte histórico do processo de formação do pensamento evolucionista para enfatizar que a associação entre biologia e darwinismo social é mais complexa do que geralmente se assume.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

Genetic Variation among Major Human Geographic Groups Supports a Peculiar Evolutionary Trend in PAX9

A total of 172 persons from nine South Amerindian, three African and one Eskimo populations were studied in relation to the Paired box gene 9 (PAX9) exon 3 (138 base pairs) as well as its 5' and 3' flanking intronic segments (232 bp and 220 bp, respectively) and integrated with the information available for the same genetic region from individuals of different geographical origins. Nine mutations were scored in exon 3 and six in its flanking regions; four of them are new South American tribe-specific singletons. Exon3 nucleotide diversity is several orders of magnitude higher than its intronic regions. Additionally, a set of variants in the PAX9 and 101 other genes related with dentition can define at least some dental morphological differences between Sub-Saharan Africans and non-Africans, probably associated with adaptations after the modern human exodus from Africa. Exon 3 of PAX9 could be a good molecular example of how evolvability works.

Testing Evolutionary and Dispersion Scenarios for the Settlement of the New World

Background: Discussion surrounding the settlement of the New World has recently gained momentum with advances in molecular biology, archaeology and bioanthropology. Recent evidence from these diverse fields is found to support different colonization scenarios. The currently available genetic evidence suggests a ""single migration'' model, in which both early and later Native American groups derive from one expansion event into the continent. In contrast, the pronounced anatomical differences between early and late Native American populations have led others to propose more complex scenarios, involving separate colonization events of the New World and a distinct origin for these groups. Methodology/Principal Findings: Using large samples of Early American crania, we: 1) calculated the rate of morphological differentiation between Early and Late American samples under three different time divergence assumptions, and compared our findings to the predicted morphological differentiation under neutral conditions in each case; and 2) further tested three dispersal scenarios for the colonization of the New World by comparing the morphological distances among early and late Amerindians, East Asians, Australo-Melanesians and early modern humans from Asia to geographical distances associated with each dispersion model. Results indicate that the assumption of a last shared common ancestor outside the continent better explains the observed morphological differences between early and late American groups. This result is corroborated by our finding that a model comprising two Asian waves of migration coming through Bering into the Americas fits the cranial anatomical evidence best, especially when the effects of diversifying selection to climate are taken into account. Conclusions: We conclude that the morphological diversity documented through time in the New World is best accounted for by a model postulating two waves of human expansion into the continent originating in East Asia and entering through Beringia.

Singing in the Rain Forest: How a Tropical Bird Song Transfers Information

How information transmission processes between individuals are shaped by natural selection is a key question for the understanding of the evolution of acoustic communication systems. Environmental acoustics predict that signal structure will differ depending on general features of the habitat. Social features, like individual spacing and mating behavior, may also be important for the design of communication. Here we present the first experimental study investigating how a tropical rainforest bird, the white-browed warbler Basileuterus leucoblepharus, extracts various information from a received song: species-specific identity, individual identity and location of the sender. Species-specific information is encoded in a resistant acoustic feature and is thus a public signal helping males to reach a wide audience. Conversely, individual identity is supported by song features susceptible to propagation: this private signal is reserved for neighbors. Finally, the receivers can locate the singers by using propagation-induced song modifications. Thus, this communication system is well matched to the acoustic constraints of the rain forest and to the ecological requirements of the species. Our results emphasize that, in a constraining acoustic environment, the efficiency of a sound communication system results from a coding/decoding process particularly well tuned to the acoustic properties of this environment.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background: The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized. Results: We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for similar to 40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance. Conclusion: These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.

Harvest-ironman: heavy armature, and not its defensive secretions, protects a harvestman against a spider

Natural selection has caused prey species to evolve distinct defensive mechanisms. One of such mechanisms was the evolution of noxious or distasteful chemicals, which have appeared independently in a number of vertebrates and invertebrates. In detailed analyses of arthropod behaviour, scent gland secretions have consistently been shown to be responsible for repelling specific predators. Because using such chemicals is costly, animals with alternative cheaper defences are expected not to release such secretions when alternative options exist. In this study, we sought to determine the defensive mechanisms of the harvestman Discocyrtus invalidus, a heavy bodied species that bears a pair of repugnatorial glands. The spider Enoploctenus cyclothorax was used as the predator, and the cricket Gryllus sp. was used as a control. In a first set of experiments, the harvestmen were preyed upon significantly less than the crickets. In two other experiments, we found that harvestmen did not use their scent gland secretions to deter the predator. Moreover, results of a fourth experiment revealed that these spiders are not repelled by defensive secretions. Discocyrtus invalidus has a thick cuticle on the entire body: scanning electron micrographs revealed that only the mouth, the articulations of appendages and the tips of the legs are not covered by a hard integument. In a fifth experiment, we found that these spiders had difficulty piercing the harvestmen body. This is the first experimental evidence that a chemically defended arachnid does not use its scent gland secretions to repel a much larger predator but instead relies on its heavily built body. (c) 2010 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

Parallel Dynamic Optimization of Steel Risers

This paper addresses the use of optimization techniques in the design of a steel riser. Two methods are used: the genetic algorithm, which imitates the process of natural selection, and the simulated annealing, which is based on the process of annealing of a metal. Both of them are capable of searching a given solution space for the best feasible riser configuration according to predefined criteria. Optimization issues are discussed, such as problem codification, parameter selection, definition of objective function, and restrictions. A comparison between the results obtained for economic and structural objective functions is made for a case study. Optimization method parallelization is also addressed. [DOI: 10.1115/1.4001955]

Absence of the HLA-G*0113N allele in Amerindian populations from the Brazilian Amazon region

The HLA-G gene is predominantly expressed at the maternal-fetal interface and has been associated with maternal-fetal tolerance. The HLA-G*0113N is a null allele defined by the insertion of a premature stop codon at exon 2, observed in a single Ghanaian individual. Likewise the G*0105N allele, the occurrence of the HL4-G*0113N in a population from an area with high pathogen load suggests that the reduced HLA-G expression in G*0113N heterozygous placentas could improve the intrauterine defense against infections. The presence of the G*0113N allele here was investigated in 150 Amerindians from five isolated tribes that inhabit the Central Amazon and in 295 admixed individuals from the State of Sao Paulo, Southeastern Brazil, previously genotyped for HLA-G. No copy of the G*0113N null allele was found in both population samples by exon 2 sequence-based analysis, reinforcing its restricted occurrence in Africa. (C) 2010 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.

Contrasting phylogenetic signals and evolutionary rates in floral traits of Neotropical lianas

The diversity of floral forms has long been considered a prime example of radiation through natural selection. However, little is still known about the evolution of floral traits, a critical piece of evidence for the understanding of the processes that may have driven flower evolution. We studied the pattern of evolution of quantitative floral traits in a group of Neotropical lianas (Bignonieae, Bignoniaceae) and used a time-calibrated phylogeny as basis to: (1) test for phylogenetic signal in 16 continuous floral traits; (2) evaluate the rate of evolution in those traits; and (3) reconstruct the ancestral state of the individual traits. Variation in floral traits among extant species of Bignonieae was highly explained by their phylogenetic history. However, opposite signals were found in floral traits associated with the attraction of pollinators (calyx and corolla) and pollen transfer (androecium and gynoecium), suggesting a differential role of selection in different floral whorls. Phylogenetic independent contrasts indicate that traits evolved at different rates, whereas ancestral character state reconstructions indicate that the ancestral size of most flower traits was larger than the mean observed sizes of the same traits in extant species. The implications of these patterns for the reproductive biology of Bignonieae are discussed. (C) 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102, 378-390.

Mating system and exclusive postzygotic paternal care in a Neotropical harvestman (Arachnida: Opiliones)

This study tests predictions of the hypothesis of evolution of paternal care via sexual selection by using the Neotropical harvestman Pseudopucrolia sp. as the model organism. Females use natural cavities in roadside banks as nesting sites, which are defended by males against other males. Females leave the nests after oviposition, and all postzygotic parental care is accomplished by males, which protect the eggs and nymphs from predators. We provided artificial mud nests to individuals in the laboratory and conducted observations on the reproduction of the species. Male reproductive success was directly related to nest ownership time: the longer a male held a nest, the higher his chances of obtaining copulations. All males that succeeded in mating and obtaining one clutch eventually mated with additional females that added eggs to the clutch. Thus, desirable males were not limited to monogamy by paternal care. Experimental manipulations demonstrated that guarding males were more attractive to females than were nonguarding males and also that males guarded unrelated eggs. Finally, we found that females and nonguarding males spent more time foraging than guarding males. We use our data to contrast hypotheses on the origin and maintenance of paternal care and to provide a critical assessment of the hypothesis of the evolution of paternal care via sexual selection. (C) 2009 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

Relation of pretreatment sequence diversity in NS5A region of HCV genotype 1 with immune response between pegylated-INF/ribavirin therapy outcomes

Hepatitis C virus (HCV) infection frequently persists despite substantial virus-specific immune responses and the combination of pegylated interferon (INF)-alpha and ribavirin therapy. Major histocompatibility complex class I restricted CD8+ T cells are responsible for the control of viraemia in HCV infection, and several studies suggest protection against viral infection associated with specific HLAs. The reason for low rates of sustained viral response (SVR) in HCV patients remains unknown. Escape mutations in response to cytotoxic T lymphocyte are widely described; however, its influence in the treatment outcome is ill understood. Here, we investigate the differences in CD8 epitopes frequencies from the Los Alamos database between groups of patients that showed distinct response to pegylated alpha-INF with ribavirin therapy and test evidence of natural selection on the virus in those who failed treatment, using five maximum likelihood evolutionary models from PAML package. The group of sustained virological responders showed three epitopes with frequencies higher than Non-responders group, all had statistical support, and we observed evidence of selection pressure in the last group. No escape mutation was observed. Interestingly, the epitope VLSDFKTWL was 100% conserved in SVR group. These results suggest that the response to treatment can be explained by the increase in immune pressure, induced by interferon therapy, and the presence of those epitopes may represent an important factor in determining the outcome of therapy.