Population genetic analysis of large sequence polymorphisms in Plasmodium falciparum blood-stage antigens

Plasmodium falciparum, the causative agent of human malaria, invades host erythrocytes using several proteins on the surface of the invasive merozoite, which have been proposed as potential vaccine candidates. Members of the multi-gene PfRh family are surface antigens that have been shown to play a central role in directing merozoites to alternative erythrocyte receptors for invasion. Recently, we identified a large structural polymorphism, a 0.58 Kb deletion, in the C-terminal region of the PfRh2b gene, present at a high frequency in parasite populations from Senegal. We hypothesize that this region is a target of humoral immunity. Here, by analyzing 371 P. falciparum isolates we show that this major allele is present at varying frequencies in different populations within Senegal, Africa, and throughout the world. For allelic dimorphisms in the asexual stage antigens, Msp-2 and EBA-175, we find minimal geographic differentiation among parasite populations from Senegal and other African localities, suggesting extensive gene flow among these populations and/or immune-mediated frequency-dependent balancing selection. In contrast, we observe a higher level of inter-population divergence (as measured by F(st)) for the PfRh2b deletion, similar to that observed for SNPs from the sexual stage Pfs45/48 loci, which is postulated to be under directional selection. We confirm that the region containing the PfRh2b polymorphism is a target of humoral immune responses by demonstrating antibody reactivity of endemic sera. Our analysis of inter-population divergence suggests that in contrast to the large allelic dimorphisms in EBA-175 and Msp-2, the presence or absence of the large PfRh2b deletion may not elicit frequency-dependent immune selection, but may be under positive immune selection, having important implications for the development of these proteins as vaccine candidates. (C) 2009 Elsevier B.V. All rights reserved.

Comparative molecular field analysis of a series of inhibitors of HIV-1 protease

Several protease inhibitors have reached the world market in the last fifteen years, dramatically improving the quality of life and life expectancy of millions of HIV-infected patients. In spite of the tremendous research efforts in this area, resistant HIV-1 variants are constantly decreasing the ability of the drugs to efficiently inhibit the enzyme. As a consequence, inhibitors with novel frameworks are necessary to circumvent resistance to chemotherapy. In the present work, we have created 3D QSAR models for a series of 82 HIV-1 protease inhibitors employing the comparative molecular field analysis (CoMFA) method. Significant correlation coefficients were obtained (q(2) = 0.82 and r(2) = 0.97), indicating the internal consistency of the best model, which was then used to evaluate an external test set containing 17 compounds. The predicted values were in good agreement with the experimental results, showing the robustness of the model and its substantial predictive power for untested compounds. The final QSAR model and the information gathered from the CoMFA contour maps should be useful for the design of novel anti-HIV agents with improved potency.

A molecular phylogeny of the stingless bee genus Melipona (Hymenoptera: Apidae)

Stingless bees (Meliponini) constitute a diverse group of highly eusocial insects that occur throughout tropical regions around the world. The meliponine genus Melipona is restricted to the New World tropics and has over 50 described species. Melipona, like Apis, possesses the remarkable ability to use representational communication to indicate the location of foraging patches. Although Melipona has been the subject of numerous behavioral, ecological, and genetic studies, the evolutionary history of this genus remains largely unexplored. Here, we implement a multigene phylogenetic approach based on nuclear, mitochondrial, and ribosomal loci, coupled with molecular clock methods, to elucidate the phylogenetic relationships and antiquity of subgenera and species of Melipona. Our phylogenetic analysis resolves the relationship among subgenera and tends to agree with morphology-based classification hypotheses. Our molecular clock analysis indicates that the genus Melipona shared a most recent common ancestor at least similar to 14-17 million years (My) ago. These results provide the groundwork for future comparative analyses aimed at understanding the evolution of complex communication mechanisms in eusocial Apidae. (C) 2010 Elsevier Inc. All rights reserved.

Unexpected genetic heterogeneity in a large consanguineous Brazilian pedigree presenting deafness

Nonsyndromic autosomal recessive deafness accounts for 80% of hereditary deafness. To date, 52 loci responsible for autosomal recessive deafness have been mapped and 24 genes identified. Here, we report a large inbred Brazilian pedigree with 26 subjects affected by prelingual deafness. Given the extensive consanguinity found in this pedigree, the most probable pattern of inheritance is autosomal recessive. However, our linkage and mutational analysis revealed, instead of an expected homozygous mutation in a single gene, two different mutant alleles and a possible third undetected mutant allele in the MYO15A gene (DFNB3 locus), as well as evidence for other causes for deafness in the same pedigree. Among the 26 affected subjects, 15 were homozygous for the novel c.10573delA mutation in the MYO15A gene, 5 were compound heterozygous for the mutation c.10573delA and the novel deletion c.9957_9960delTGAC and one inherited only a single c.10573delA mutant allele, while the other one could not be identified. Given the extensive consanguinity of the pedigree, there might be at least one more deafness locus segregating to explain the condition in some of the subjects whose deafness is not clearly associated with MYO15A mutations, although overlooked environmental causes could not be ruled out. Our findings illustrate a high level of etiological heterogeneity for deafness in the family and highlight some of the pitfalls of genetic analysis of large genes in extended pedigrees, when homozygosity for a single mutant allele is expected.

PODOCARPUS SP.PL. REFUGIA IN BRAZIL`S ATLANTIC FOREST: ANALYSING THE PAST TO IMPROVE FUTURE PROTECTION

The Atlantic rainforest has the second highest biodiversity in Brazil. It has been shrinking rapidly in area as a result of intensive deforestation, and only 7% of the original cover now remains, as isolated patches or in ecological reserves. In order to obtain new information on the distribution of the Atlantic rainforest during the Quaternary, we examined herbarium data to locate relevant populations and extracted DNA from fresh leaves from 26 populations. The present-day distribution of endemic Podocarpus populations shows that they are widely dispersed across eastern Brazil, and that the expansion of Podocarpus recorded in single Amazonian pollen records may have originated from either western or eastern populations. Genetic analysis enabled us to determine the boundaries of their regional expansion: northern and central populations of P. sellowii appeared between 5 degrees and 15 degrees S some 16,000 years ago; populations of P lambertii or sellowii have appeared between 15 degrees and 23 degrees S at different times since the last glaciation at least; and P lambertii appeared between 23 degrees and 30 degrees S during the recent expansion of Araucaria forests. The combination of botanical, pollen, and molecular analyses proved to be a rapid means of inferring distribution boundaries for sparse populations and their regional evolution within tropical ecosystems. Today the rainforest refugia we identified have become hotspots that are crucial to the survival of the Atlantic forest under unfavourable climatic conditions and, as such, offer the only possible opportunity for this type of forest to expand in the event of future climate change.

Successful carnivore identification with faecal DNA across a fragmented Amazonian landscape

The use of scat surveys to obtain DNA has been well documented in temperate areas, where DNA preservation may be more effective than in tropical forests. Samples obtained in the tropics are often exposed to high humidity, warm temperatures, frequent rain and intense sunlight, all of which can rapidly degrade DNA. Despite these potential problems, we demonstrate successful mtDNA amplification and sequencing for faeces of carnivores collected in tropical conditions and quantify how sample condition and environmental variables influence the success of PCR amplification and species identification. Additionally, the feasibility of genotyping nuclear microsatellites from jaguar (Panthera onca) faeces was investigated. From October 2007 to December 2008, 93 faecal samples were collected in the southern Brazilian Amazon. A total of eight carnivore species was successfully identified from 71% of all samples obtained. Information theoretic analysis revealed that the number of PCR attempts before a successful sequence was an important negative predictor across all three responses (success of species identification, success of species identification from the first sequence and PCR amplification success), whereas the relative importance of the other three predictors (sample condition, season and distance from forest edge) varied between the three responses. Nuclear microsatellite amplification from jaguar faeces had lower success rates (15-44%) compared with those of the mtDNA marker. Our results show that DNA obtained from faecal samples works efficiently for carnivore species identification in the Amazon forest and also shows potential for nuclear DNA analysis, thus providing a valuable tool for genetic, ecological and conservation studies.

Skull Modularity in Neotropical Marsupials and Monkeys: Size Variation and Evolutionary Constraint and Flexibility

An organism is built through a series of contingent factors, yet it is determined by historical, physical, and developmental constraints. A constraint should not be understood as an absolute obstacle to evolution, as it may also generate new possibilities for evolutionary change. Modularity is, in this context, an important way of organizing biological information and has been recognized as a central concept in evolutionary biology bridging on developmental, genetics, morphological, biochemical, and physiological studies. In this article, we explore how modularity affects the evolution of a complex system in two mammalian lineages by analyzing correlation, variance/covariance, and residual matrices (without size variation). We use the multivariate response to selection equation to simulate the behavior of Eutheria and Metharia skulls in terms of their evolutionary flexibility and constraints. We relate these results to classical approaches based on morphological integration tests based on functional/developmental hypotheses. Eutherians (Neotropical primates) showed smaller magnitudes of integration compared with Metatheria (didelphids) and also skull modules more clearly delimited. Didelphids showed higher magnitudes of integration and their modularity is strongly influenced by within-groups size variation to a degree that evolutionary responses are basically aligned with size variation. Primates still have a good portion of the total variation based on size; however, their enhanced modularization allows a broader spectrum of responses, more similar to the selection gradients applied (enhanced flexibility). Without size variation, both groups become much more similar in terms of modularity patterns and magnitudes and, consequently, in their evolutionary flexibility. J. Exp. Zool. (Mol. Dev. Evol.) 314B:663-683, 2010. (C) 2010 Wiley-Liss, Inc.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Structural basis for selective inhibition of trypanosomatid glyceraldehyde-3-phosphate dehydrogenase: Molecular docking and 3D QSAR studies

The glycolytic enzyme glyceraldehyde-3 -phosphate dehydrogenase (GAPDH) is as an attractive target for the development of novel antitrypanosomatid agents. In the present work, comparative molecular field analysis and comparative molecular similarity index analysis were conducted on a large series of selective inhibitors of trypanosomatid GAPDH. Four statistically significant models were obtained (r(2) > 0.90 and q(2) > 0.70), indicating their predictive ability for untested compounds. The models were then used to predict the potency of an external test set, and the predicted values were in good agreement with the experimental results. Molecular modeling studies provided further insight into the structural basis for selective inhibition of trypanosomatid GAPDH.

Insights into the Molecular Requirements for the Anti-obesity Activity of a Series of CB1 Ligands

Two-dimensional and 3D quantitative structure-activity relationships studies were performed on a series of diarylpyridines that acts as cannabinoid receptor ligands by means of hologram quantitative structure-activity relationships and comparative molecular field analysis methods. The quantitative structure-activity relationships models were built using a data set of 52 CB1 ligands that can be used as anti-obesity agents. Significant correlation coefficients (hologram quantitative structure-activity relationships: r 2 = 0.91, q 2 = 0.78; comparative molecular field analysis: r 2 = 0.98, q 2 = 0.77) were obtained, indicating the potential of these 2D and 3D models for untested compounds. The models were then used to predict the potency of an external test set, and the predicted (calculated) values are in good agreement with the experimental results. The final quantitative structure-activity relationships models, along with the information obtained from 2D contribution maps and 3D contour maps, obtained in this study are useful tools for the design of novel CB1 ligands with improved anti-obesity potency.

Resurrection of Two Species From Synonymy of Anopheles (Nyssorhynchus) strodei Root, and Characterization of a Distinct Morphological Form From the Strodei Complex (Diptera: Culicidae)

Anopheles albertoi Unti and Anopheles arthuri Unti are revived from the synonymy with Anopheles strodei Root, and a distinct morphological form (designated in this study as Anopheles CP Form) from the Strodei Complex of Anopheles (Nyssorhynchus) is characterized. The male genitalia of An. arthuri and An. albertoi are described and illustrated for the first time. An. strodei, An. arthuri, and An. albertoi were first distinguished based on scanning electron microphotos of the eggs, and then each egg type was associated with diagnostic characters of the male genitalia. Identification of Anopheles CP Form was based on morphological characters of the male genitalia, characterized and illustrated in this study. Molecular phylogenetic analysis was most clear when an outgroup was not included, in which case using the nuclear white gene, or the white gene in combination with the mitochondrial cytochrome c oxidase subunit I (COI) gene, clearly separated these four taxa. When Anopheles quadrimaculatus Say and Anopheles stephensi Liston were included as an outgroup, combined white and COI data resolved An. strodei and An. albertoi, whereas An. arthuri was not well resolved. The single sequence of Anopheles CP Form was recovered well separated from other groups in all analyses.

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

Phylogenetic relationships, diversification and biogeography in Neotropical Brotogeris parakeets

Aim We present a molecular phylogenetic analysis of Brotogeris (Psittacidae) using several distinct and complementary approaches: we test the monophyly of the genus, delineate the basal taxa within it, uncover their phylogenetic relationships, and finally, based on these results, we perform temporal and spatial comparative analyses to help elucidate the historical biogeography of the Neotropical region. Location Neotropical lowlands, including dry and humid forests. Methods Phylogenetic relationships within Brotogeris were investigated using the complete sequences of the mitochondrial genes cyt b and ND2, and partial sequences of the nuclear intron 7 of the gene for Beta Fibrinogen for all eight species and 12 of the 17 taxa recognized within the genus (total of 63 individuals). In order to delinetae the basal taxa within the genus we used both molecular and plumage variation, the latter being based on the examination of 597 skin specimens. Dates of divergence and confidence intervals were estimated using penalized likelihood. Spatial and temporal comparative analyses were performed including several closely related parrot genera. Results Brotogeris was found to be a monophyletic genus, sister to Myiopsitta. The phylogenetic analyses recovered eight well-supported clades representing the recognized biological species. Although some described subspecies are diagnosably distinct based on morphology, there was generally little intraspecific mtDNA variation. The Amazonian species had different phylogenetic affinities and did not group in a monophyletic clade. Brotogeris diversification took place during the last 6 Myr, the same time-frame as previously found for Pionus and Pyrilia. Main conclusions The biogeographical history of Brotogeris implies a dynamic history for South American biomes since the Pliocene. It corroborates the idea that the geological evolution of Amazonia has been important in shaping its biodiversity, argues against the idea that the region has been environmentally stable during the Quaternary, and suggests dynamic interactions between wet and dry forest habitats in South America, with representatives of the Amazonian biota having several independent close relationships with taxa endemic to other biomes.

Covariance structure in the skull of Catarrhini: a case of pattern stasis and magnitude evolution

The study of the genetic variance/covariance matrix (G-matrix) is a recent and fruitful approach in evolutionary biology, providing a window of investigating for the evolution of complex characters. Although G-matrix studies were originally conducted for microevolutionary timescales, they could be extrapolated to macroevolution as long as the G-matrix remains relatively constant, or proportional, along the period of interest. A promising approach to investigating the constancy of G-matrices is to compare their phenotypic counterparts (P-matrices) in a large group of related species; if significant similarity is found among several taxa, it is very likely that the underlying G-matrices are also equivalent. Here we study the similarity of covariance and correlation structure in a broad sample of Old World monkeys and apes (Catarrhini). We made phylogenetically structured comparisons of correlation and covariance matrices derived from 39 skull traits, ranging from between species to the superfamily level. We also compared the overall magnitude of integration between skull traits (r(2)) for all Catarrhim genera. Our results show that P-matrices were not strictly constant among catarrhines, but the amount of divergence observed among taxa was generally low. There was significant and positive correlation between the amount of divergence in correlation and covariance patterns among the 30 genera and their phylogenetic distances derived from a recently proposed phylogenetic hypothesis. Our data demonstrate that the P-matrices remained relatively similar along the evolutionary history of catarrhines, and comparisons with the G-matrix available for a New World monkey genus (Saguinus) suggests that the same holds for all anthropoids. The magnitude of integration, in contrast, varied considerably among genera, indicating that evolution of the magnitude, rather than the pattern of inter-trait correlations, might have played an important role in the diversification of the catarrhine skull. (C) 2009 Elsevier Ltd. All rights reserved.

Morphology and phylogenetic relationships of a remarkable new genus and two new species of Neotropical freshwater stingrays from the Amazon basin (Chondrichthyes: Potamotrygonidae)

The morphology and phylogenetic relationships of a new genus and two new species of Neotropical freshwater stingrays, family Potamotrygonidae, are investigated and described in detail. The new genus, Heliotrygon, n. gen., and its two new species, Heliotrygon gomesi, n. sp. (type-species) and Heliotrygon rosai, n. sp., are compared to all genera and species of potamotrygonids, based on revisions in progress. Some of the derived features of Heliotrygon include its unique disc proportions (disc highly circular, convex anteriorly at snout region, its width and length very similar), extreme subdivision of suborbital canal (forming a complex honeycomb-like pattern anterolaterally on disc), stout and triangular pelvic girdle, extremely reduced caudal sting, basibranchial copula with very slender and acute anterior extension, and precerebral and frontoparietal fontanellae of about equal width, tapering very little posteriorly. Both new species can be distinguished by their unique color patterns: Heliotrygon gomesi is uniform gray to light tan or brownish dorsally, without distinct patterns, whereas Heliotrygon rosai is characterized by numerous white to creamy-white vermiculate markings over a light brown, tan or gray background color. Additional proportional characters that may further distinguish both species are also discussed. Morphological descriptions are provided for dermal denticles, ventral lateral-line canals, skeleton, and cranial, hyoid and mandibular muscles of Heliotrygon, which clearly corroborate it as the sister group of Paratrygon. Both genera share numerous derived features of the ventral lateral-line canals, neurocranium, scapulocoracoid, pectoral basals, clasper morphology, and specific patterns of the adductor mandibulae and spiracularis medialis muscles. Potamotrygon and Plesiotrygon are demonstrated to share derived characters of their ventral lateral-line canals, in addition to the presence of angular cartilages. Our morphological phylogeny is further corroborated by a molecular phylogenetic analysis of cytochrome b based on four sequences (637 base pairs in length), representing two distinct haplotypes for Heliotrygon gomesi. Parsimony analysis produced a single most parsimonious tree revealing Heliotrygon and Paratrygon as sister taxa (boot-strap proportion of 70%), which together are the sister group to a clade including Plesiotrygon and species of Potamotrygon. These unusual stingrays highlight that potamotrygonid diversity, both in terms of species composition and undetected morphological and molecular patterns, is still poorly known.