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.

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.

Coalescent analysis of mtDNA indicates Pleistocene divergence among three species of howler monkey (Alouatta spp.) and population subdivision within the Atlantic Coastal Forest species, A. guariba

We have used coalescent analysis of mtDNA cytochrome b (cyt b) sequences to estimate times of divergence of three species of Alouatta-A. caraya, A. belzebul, and A. guariba-which are in close geographic proximity. A. caraya is inferred to have diverged from the A. guariba/A. belzebul clade approximately 3.83 million years ago (MYA), with the later pair diverging approximately 1.55 MYA. These dates are much more recent than previous dates based on molecular-clock methods. In addition, analyses of new sequences from the Atlantic Coastal Forest species A. guariba indicate the presence of two distinct haplogroups corresponding to northern and southern populations with both haplogroups occurring in sympatry within Sao Paulo state. The time of divergence of these two haplogroups is estimated to be 1.2 MYA and so follows quite closely after the divergence of A. guariba and A. belzebul. These more recent dates point to the importance of Pleistocene environmental events as important factors in the diversification of A. belzebul and A. guariba. We discuss the diversification of the three Alouatta species in the context of recent models of climatic change and with regard to recent molecular phylogeographic analyses of other animal groups distributed in Brazil.

Divergence Involving Global Regulatory Gene Mutations in an Escherichia coli Population Evolving under Phosphate Limitation

Many of the important changes in evolution are regulatory in nature. Sequenced bacterial genomes point to flexibility in regulatory circuits but we do not know how regulation is remodeled in evolving bacteria. Here, we study the regulatory changes that emerge in populations evolving under controlled conditions during experimental evolution of Escherichia coli in a phosphate-limited chemostat culture. Genomes were sequenced from five clones with different combinations of phenotypic properties that coexisted in a population after 37 days. Each of the distinct isolates contained a different mutation in 1 of 3 highly pleiotropic regulatory genes (hfq, spoT, or rpoS). The mutations resulted in dissimilar proteomic changes, consistent with the documented effects of hfq, spoT, and rpoS mutations. The different mutations do share a common benefit, however, in that the mutations each redirect cellular resources away from stress responses that are redundant in a constant selection environment. The hfq mutation lowers several individual stress responses as well the small RNA-dependent activation of rpoS translation and hence general stress resistance. The spoT mutation reduces ppGpp levels, decreasing the stringent response as well as rpoS expression. The mutations in and upstream of rpoS resulted in partial or complete loss of general stress resistance. Our observations suggest that the degeneracy at the core of bacterial stress regulation provides alternative solutions to a common evolutionary challenge. These results can explain phenotypic divergence in a constant environment and also how evolutionary jumps and adaptive radiations involve altered gene regulation.

Population genetic structuring in pacu (Piaractus mesopotamicus) across the Paraná-Paraguay basin: evidence from microsatellites

The Paraná-Paraguay basin encompasses central western Brazil, northeastern Paraguay, eastern Bolivia and northern Argentina. The Pantanal is a flooded plain with marked dry and rainy seasons that, due to its soil characteristics and low declivity, has a great water holding capacity supporting abundant fish fauna. Piaractus mesopotamicus, or pacu, endemic of the Paraná-Paraguay basin, is a migratory species economically important in fisheries and ecologically as a potential seed disperser. In this paper we employ eight microsatellite loci to assess the population structure of 120 pacu sampled inside and outside the Pantanal of Mato Grosso. Our main objective was to test the null hypothesis of panmixia and to verify if there was a different structuring pattern between the Pantanal were there were no physical barriers to fish movement and the heavily impounded Paraná and Paranapanema rivers. All loci had moderate to high levels of polymorphism, the number of alleles varied from three to 18. The average observed heterozygosity varied from 0.068 to 0.911. After the Bonferroni correction three loci remained significant for deviations from Hardy-Weinberg, and for those the frequency of null alleles was estimated. F ST and R ST pairwise comparisons detected low divergence among sampling sites, and differentiation was significant only between Paranapanema and Cuiabá and Paranapanema and Taquari. No correlation between genetic distance and the natural logarithm of the geographic distance was detected. Results indicate that for conservation purposes and for restoration programs small genetic differences detected in the Cuiabá and Paranapanema rivers should be taken in consideration.

Lineage divergence detected in the malaria vector Anopheles marajoara (Diptera: Culicidae) in Amazonian Brazil

Background: Cryptic species complexes are common among anophelines. Previous phylogenetic analysis based on the complete mtDNA COI gene sequences detected paraphyly in the Neotropical malaria vector Anopheles marajoara. The ""Folmer region"" detects a single taxon using a 3% divergence threshold. Methods: To test the paraphyletic hypothesis and examine the utility of the Folmer region, genealogical trees based on a concatenated (white + 3' COI sequences) dataset and pairwise differentiation of COI fragments were examined. The population structure and demographic history were based on partial COI sequences for 294 individuals from 14 localities in Amazonian Brazil. 109 individuals from 12 localities were sequenced for the nDNA white gene, and 57 individuals from 11 localities were sequenced for the ribosomal DNA (rDNA) internal transcribed spacer 2 (ITS2). Results: Distinct A. marajoara lineages were detected by combined genealogical analysis and were also supported among COI haplotypes using a median joining network and AMOVA, with time since divergence during the Pleistocene (< 100,000 ya). COI sequences at the 3' end were more variable, demonstrating significant pairwise differentiation (3.82%) compared to the more moderate 2.92% detected by the Folmer region. Lineage 1 was present in all localities, whereas lineage 2 was restricted mainly to the west. Mismatch distributions for both lineages were bimodal, likely due to multiple colonization events and spatial expansion (similar to 798 - 81,045 ya). There appears to be gene flow within, not between lineages, and a partial barrier was detected near Rio Jari in Amapa state, separating western and eastern populations. In contrast, both nDNA data sets (white gene sequences with or without the retention of the 4th intron, and ITS2 sequences and length) detected a single A. marajoara lineage. Conclusions: Strong support for combined data with significant differentiation detected in the COI and absent in the nDNA suggest that the divergence is recent, and detectable only by the faster evolving mtDNA. A within subgenus threshold of >2% may be more appropriate among sister taxa in cryptic anopheline complexes than the standard 3%. Differences in demographic history and climatic changes may have contributed to mtDNA lineage divergence in A. marajoara.

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

Allozymatic divergence between border populations of two cryptic species of the Drosophila buzzatii cluster species (Diptera: Drosophilidae)

Drosophila antonietae and Drosophila gouveai are allopatric, cactophilic, cryptic and endemic of South America species, which aedeagus morphology is considered the main diagnostic character. In this work, single close populations from the edge distributions of each species, located in an ""introgressive corridor"", were analyzed regarding temporal isozenzymatic genetic variability. Isocitrate dehydrogenase (Idh) appeared as a diagnostic locus between D. antonieate and D. gouveai because each population was fixed for different alleles. Moreover, several polymorphic loci showed accentuated divergence in the allele frequency, as evidenced by Nei`s l(0.3188) and D (1.1432), and also by Reynolds` genetic distance and identity (1.3207 and 0.7331, respectively). Our results showed that, in spite of the very similar external morphology, related evolutionary histories, close distributions, and events of introgression in the studied area, these cryptic species have high allozymatic differentiation, and this is discussed here. (C) 2010 Elsevier Ltd. All rights reserved.

Phylogeography of the cactophilic species Drosophila gouveai: demographic events and divergence timing in dry vegetation enclaves in eastern Brazil

Aim The aim of this study was to assess the causal mechanisms underlying populational subdivision in Drosophila gouveai, a cactophilic species associated with xeric vegetation enclaves in eastern Brazil. A secondary aim was to investigate the genetic effects of Pleistocene climatic fluctuations on these environments. Location Dry vegetation enclaves within the limits of the Cerrado domain in eastern Brazil. Methods We determined the mitochondrial DNA haplotypes of 55 individuals (representing 12 populations) based on sequence data of a 483-bp fragment from the cytochrome c oxidase subunit II (COII) gene. Phylogenetic and coalescent analyses were used to test for the occurrence of demographic events and to infer the time of divergence amongst genetically independent groups. Results Our analyses revealed the existence of two divergent subclades (G1 and G2) plus an introgressed clade restricted to the southernmost range of D. gouveai. Subclades G1 and G2 displayed genetic footprints of range expansion and segregated geographical distributions in south-eastern and some central highland regions, east and west of the Parana River valley. Molecular dating indicated that the main demographic and diversification events occurred in the late to middle Pleistocene. Main conclusions The phylogeographical and genetic patterns observed for D. gouveai in this study are consistent with changes in the distribution of dry vegetation in eastern Brazil. All of the estimates obtained by molecular dating indicate that range expansion and isolation pre-dated the Last Glacial Maximum, occurring during the late to middle Pleistocene, and were probably triggered by climatic changes during the Pleistocene. The current patchy geographical distribution and population subdivision in D. gouveai is apparently closely linked to these past events.

Population structure of the migratory fish Prochilodus lineatus (Characiformes) from rio Grande basin (Brazil), an area fragmented by dams

1. Prochilodus lineatus (Prochilodontidae, Characiformes) is a migratory species of great economic importance both in fisheries and aquaculture that is found throughout the Jacui, Paraiba do Sul, Parana, Paraguay and Uruguay river basins in South America. Earlier population studies of P. lineatus in the rio Grande basin (Parana basin) indicated the existence of a single population; however, the range of this species has been fragmented by the construction of several dams. Such dams modified the environmental conditions and could have constrained the reproductive migration of P. lineatus, possibly leading to changes in the population genetic structure. 2. In order to evaluate how genetic diversity is allocated in the rio Grande basin, 141 specimens of P. lineatus from eight collection sites were analysed using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) with 15 restriction enzymes. 3. Forty-six haplotypes were detected, and 70% of them are restricted. The mean genetic variability indexes (h = 0.7721 and pi = 1.6%) were similar to those found in natural populations with a large effective size. Fst and Exact Test values indicated a lack of structuring among the samples, and the model of isolation by distance was tested and rejected. 4. The haplotype network indicated that this population of P. lineatus has been maintained as a single variable stock with some differences in the genetic composition (haplotypes) between samples. Indications of population expansion were detected, and this finding was supported by neutrality tests and mismatch distribution analyses. 5. The present study focused on regions between dams to serve as a parameter for further evaluations of genetic variability and the putative impact of dams and repopulation programmes in natural populations of P. lineatus. Copyright (C) 2011 John Wiley & Sons, Ltd.

The CATH classification revisited-architectures reviewed and new ways to characterize structural divergence in superfamilies

The latest version of CATH (class, architecture, topology, homology) (version 3.2), released in July 2008 (http://www.cathdb.info), contains 1 14215 domains, 2178 Homologous superfamilies and 1110 fold groups. We have assigned 20 330 new domains, 87 new homologous superfamilies and 26 new folds since CATH release version 3.1. A total of 28 064 new domains have been assigned since our NAR 2007 database publication (CATH version 3.0). The CATH website has been completely redesigned and includes more comprehensive documentation. We have revisited the CATH architecture level as part of the development of a `Protein Chart` and present information on the population of each architecture. The CATHEDRAL structure comparison algorithm has been improved and used to characterize structural diversity in CATH superfamilies and structural overlaps between superfamilies. Although the majority of superfamilies in CATH are not structurally diverse and do not overlap significantly with other superfamilies, similar to 4% of superfamilies are very diverse and these are the superfamilies that are most highly populated in both the PDB and in the genomes. Information on the degree of structural diversity in each superfamily and structural overlaps between superfamilies can now be downloaded from the CATH website.

Prevalence and risk factors of apical periodontitis in endodontically treated teeth in a selected population of Brazilian adults

The aim of this study was to assess the prevalence and risk factors of apical periodontitis in endodontically treated teeth in a selected population of Brazilian adults. A total of 1,372 periapical radiographs of endodontically treated teeth were analyzed based on the quality of root filling, status of coronal restoration and presence of posts associated with apical periodontitis (AP). Data were analyzed statistically using odds ratio, confidence intervals and chi-square test. The prevalence of AP with adequate endodontic treatment was low (16.5%). This percentage dropped to 12.1% in cases with adequate root filling and adequate coronal restoration. Teeth with adequate endodontic treatment and poor coronal restoration had an AP prevalence of 27.9%. AP increased to 71.7% in teeth with poor endodontic treatment associated with poor coronal restoration. When poor endodontic treatment was combined with adequate coronal restoration, AP prevalence was 61.8%. The prevalence of AP was low when associated with high technical quality of root canal treatment. Poor coronal restoration increased the risk of AP even when endodontic treatment was adequate (OR=2.80; 95%CI=1.87-4.22). The presence of intracanal posts had no influence on AP prevalence.

Risk indicators for aggressive periodontitis in an untreated isolated young population from Brazil

This study aimed to assess the prevalence of aggressive periodontitis (AgP), and to investigate the association between demographic, socioeconomic and behavioral risk indicators with AgP in an untreated and isolated young population in Southeastern Brazil. For this cross-sectional survey, 134 subjects aged 12-29 years were selected by a census. Of those eligible, 101 subjects received a full-mouth clinical examination, and were interviewed using a structured written questionnaire. Cases were defined as individuals with 4 or more teeth with attachment loss > 4 mm or > 5 mm in the age groups 12-19 and 20-29, respectively. Overall, 9.9% of the subjects presented AgP (10.3% of the 12-19-year-olds and 9.7% of the 20-29-year-olds). The only risk indicator significantly associated with AgP in this isolated population was a high proportion of sites (> 30%) presenting supragingival calculus [OR = 23.2]. Having experienced an urgency dental treatment was a protective factor for AgP [OR = 0.1]. The authors concluded that this isolated and untreated population from Brazil presented a high prevalence of AgP. Local plaque-retaining factors played a major role in the prevalence of AgP in this isolated population, and should be included in further studies evaluating this destructive periodontal disease form.

Effects of dietary energy and nitrogen supplements on rumen fermentation and protozoa population in buffalo and zebu cattle

The effects were assessed of two energy sources in concentrate (ground grain corn vs. citrus pulp) and two nitrogen sources (soybean meal vs. urea) on rumen metabolism in four buffaloes and four zebu cattle (Nellore) with rumen cannula and fed in a 4 × 4 Latin square design with feeds containing 60% sugar cane. Energy supplements had no effect on the rumen ammonia concentration in cattle, but ground grain corn promoted higher ammonia level than citrus pulp in buffalo. Urea produced higher ammonia level than soybean meal in both animal species. On average, the buffaloes maintained a lower rumen ammonia concentration (11.7 mg/dL) than the cattle (14.5 mg/dL). Buffaloes had lower production of acetic acid than cattle (58.7 vs. 61.6 mol/100 mol) and higher of propionic acid (27.4 vs. 23.6 mol/100 mol). There was no difference in the butyric acid production between the buffaloes (13.6 mol/100 mol) and cattle (14.8 mol/100 mol) and neither in the total volatile fatty acids concentration (82.5 vs. 83.6 mM, respectively). The energy or nitrogen sources had no effect on rumen protozoa count in either animal species. The zebu cattle had higher rumen protozoa population (8.8 × 10(5)/mL) than the buffaloes (6.1 × 10(5)/mL). The rumen protozoa population differed between the animal species, except for Dasytricha and Charonina. The buffaloes had a lower Entodinium population than the cattle (61.0 vs 84.9%, respectively) and a greater percentage of species belonging to the Diplodiniinae subfamily than the cattle (28.6 vs. 1.4%, respectively). In cattle, ground corn is a better energy source than citrus pulp for use by Entodinium and Diplodiniinae. In the buffaloes, the Entodinium are favored by urea and Diplodiniinae species by soybean meal.

Population dynamics: seasonal variation of phytoplankton functional groups in brazilian reservoirs (Billings and Guarapiranga, São Paulo)

Phytoplankton may function as a "sensor" of changes in aquatic environment and responds rapidly to such changes. In freshwaters, coexistence of species that have similar ecological requirements and show the same environmental requirements frequently occurs; such species groups are named functional groups. The use of phytoplankton functional groups to evaluate these changes has proven to be very useful and effective. Thus, the aim of this study was to evaluate the occurrence of functional groups of phytoplankton in two reservoirs (Billings and Guarapiranga) that supply water to millions of people in São Paulo city Metropolitan Area, southeastern Brazil. Surface water samples were collected monthly and physical, chemical and biological (quantitative and qualitative analyses of the phytoplankton) were performed. The highest biovolume (mm³.L-1) of the descriptor species and functional groups were represented respectively by Anabaena circinalis Rabenh. (H1), Microcystis aeruginosa (Kützing) Kützing (L M/M) and Mougeotia sp. (T) in the Guarapiranga reservoir and Cylindrospermopsis raciborskii (Wolosz.) Seen. and Subba Raju (S N), Microcystis aeruginosa and M. panniformis Komárek et al. (L M/M), Planktothrix agardhii (Gom.) Anagn. and Komárek and P. cf. clathrata (Skuja) Anagn. and Komárek (S1) in the Billings reservoir. The environmental factors that most influenced the phytoplankton dynamics were water temperature, euphotic zone, turbidity, conductivity, pH, dissolved oxygen, nitrate and total phosphorous.