A note on the accuracy of PAC-likelihood inference with microsatellite data

Stephens and Donnelly have introduced a simple yet powerful importance sampling scheme for computing the likelihood in population genetic models. Fundamental to the method is an approximation to the conditional probability of the allelic type of an additional gene, given those currently in the sample. As noted by Li and Stephens, the product of these conditional probabilities for a sequence of draws that gives the frequency of allelic types in a sample is an approximation to the likelihood, and can be used directly in inference. The aim of this note is to demonstrate the high level of accuracy of "product of approximate conditionals" (PAC) likelihood when used with microsatellite data. Results obtained on simulated microsatellite data show that this strategy leads to a negligible bias over a wide range of the scaled mutation parameter theta. Furthermore, the sampling variance of likelihood estimates as well as the computation time are lower than that obtained with importance sampling on the whole range of theta. It follows that this approach represents an efficient substitute to IS algorithms in computer intensive (e.g. MCMC) inference methods in population genetics. (c) 2006 Elsevier Inc. All rights reserved.

The evolutionary ecology of resistance to parasitoids by Drosophila

Parasitoids are the most important natural enemies of many insect species. Larvae of many Drosophila species can defend themselves against attack by parasitoids through a cellular immune response called encapsulation. The paper reviews recent studies of the evolutionary biology and ecological genetics of resistance in Drosophila, concentrating on D. melanogaster. The physiological basis of encapsulation, and the genes known to interfere with resistance are briefly summarized. Evidence for within- and between-population genetic variation in resistance from isofemale line, artificial selection and classical genetic studies are reviewed. There is now firm evidence that resistance is costly to Drosophila, and the nature of this cost is discussed, and the possibility that it may involve a reduction in metabolic rate considered. Comparative data on encapsulation and metabolic rates across seven Drosophila species provides support for this hypothesis. Finally, the possible population and community ecological consequences of evolution in the levels of host resistance are examined.

Adaptive divergence and speciation among sexual and pseudoviviparous populations of Festuca

Pseudovivipary is an environmentally induced flowering abnormality in which vegetative shoots replace seminiferous (sexual) inflorescences. Pseudovivipary is usually retained in transplantation experiments, indicating that the trait is not solely induced by the growing environment. Pseudovivipary is the defining characteristic of Festuca vivipara, and arguably the only feature separating this species from its closest seminiferous relative, Festuca ovina. We performed phylogenetic and population genetic analysis on sympatric F. ovina and F. vivipara samples to establish whether pseudovivipary is an adaptive trait that accurately defines the separation of genetically distinct Festuca species. Chloroplast and nuclear marker-based analyses revealed that variation at a geographical level can exceed that between F. vivipara and F. ovina. We deduced that F. vivipara is a recent species that frequently arises independently within F. ovina populations and has not accumulated significant genetic differentiation from its progenitor. We inferred local gene flow between the species. We identified one amplified fragment length polymorphism marker that may be linked to a pseudovivipary-related region of the genome, and several other markers provide evidence of regional local adaptation in Festuca populations. We conclude that F. vivipara can only be appropriately recognized as a morphologically and ecologically distinct species; it lacks genetic differentiation from its relatives. This is the first report of a ‘failure in normal flowering development’ that repeatedly appears to be adaptive, such that the trait responsible for species recognition constantly reappears on a local basis.

Evolution of discrimination in populations at equilibrium between selfishness and altruism

Where there is genetically based variation in selfishness and altruism, as in man, altruists with an innate ability to recognise and thereby only help their altruistic relatives may evolve. Here we use diploid population genetic models to chart the evolution of genetically-based discrimination in populations initially in stable equilibrium between altruism and selfishness. The initial stable equilibria occur because help is assumed subject to diminishing returns. Similar results were obtained whether we used a model with two independently inherited loci, one controlling altruism the other discrimination, or a one locus model with three alleles. The latter is the opposite extreme to the first model, and can be thought of as involving complete linkage between two loci on the same chromosome. The introduction of discrimination reduced the benefits obtained by selfish individuals, more so as the number of discriminators increased, and selfishness was eventually eliminated in some cases. In others selfishness persisted and the evolutionary outcome was a stable equilibrium involving selfish individuals and both discriminating and non-discriminating altruists. Heritable variation in selfishness, altruism and discrimination is predicted to be particularly evident among full sibs. The suggested coexistence of these three genetic dispositions could explain widespread interest within human social groups as to who will and who will not help others. These predictions merit experimental and observational investigation by primatologists, anthropologists and psychologists. Keywords: Population genetics, Diploid, Heritability, Prosocial, Behaviour genetics

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

Development of new nuclear markers and characterization of single nucleotide polymorphisms in kelp gull (Larus dominicanus)

The present study seeks to develop nuclear markers for the kelp gull (Larus dominicanus). We hereby report the characterization of 12 independent nuclear introns, where 104 single nucleotide polymorphisms (SNPs) in 8138 sequenced base pairs were observed. These SNP markers are the first to be designed for genotyping a gull species. The markers will provide useful tools for understanding which processes act or acted upon kelp gulls to cause their low genetic variability in mitochondrial DNA. In addition, these markers open a new opportunity for population genetic and evolutionary studies in the Laridae group.

Evolution of Dendrocolaptes platyrostris (Aves: Furnariidae) between the South American open vegetation corridor and the Atlantic forest

The open vegetation corridor of South America is a region dominated by savanna biomes. It contains forests (i.e. riverine forests) that may act as corridors for rainforest specialists between the open vegetation corridor and its neighbouring biomes (i.e. the Amazonian and Atlantic forests). A prediction for this scenario is that populations of rainforest specialists in the open vegetation corridor and in the forested biomes show no significant genetic divergence. We addressed this hypothesis by studying plumage and genetic variation of the Planalto woodcreeper Dendrocolaptes platyrostris Spix (1824) (Aves: Furnariidae), a forest specialist that occurs in both open habitat and in the Atlantic forest. The study questions were: (1) is there any evidence of genetic continuity between populations of the open habitat and the Atlantic forest and (2) is plumage variation congruent with patterns of neutral genetic structure or with ecological factors related to habitat type? We used cytochrome b and mitochondrial DNA control region sequences to show that D. platyrostris is monophyletic and presents substantial intraspecific differentiation. We found two areas of plumage stability: one associated with Cerrado and the other associated with southern Atlantic Forest. Multiple Mantel tests showed that most of the plumage variation followed the transition of habitats but not phylogeographical gaps, suggesting that selection may be related to the evolution of the plumage of the species. The results were not compatible with the idea that forest specialists in the open vegetation corridor and in the Atlantic forest are linked at the population level because birds from each region were not part of the same genetic unit. Divergence in the presence of gene flow across the ecotone between both regions might explain our results. Also, our findings indicate that the southern Atlantic forest may have been significantly affected by Pleistocene climatic alteration, although such events did not cause local extinction of most taxa, as occurred in other regions of the globe where forests were significantly affected by global glaciations. Finally, our results neither support plumage stability areas, nor subspecies as full species. (C) 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103, 801-820.

Plasmodium falciparum Accompanied the Human Expansion out of Africa

Plasmodium falciparum is distributed throughout the tropics and is responsible for an estimated 230 million cases of malaria every year, with a further 1.4 billion people at risk of infection [1-3]. Little is known about the genetic makeup of P. falciparum populations, despite variation in genetic diversity being a key factor in morbidity, mortality, and the success of malaria control initiatives. Here we analyze a worldwide sample of 519 P. falciparum isolates sequenced for two housekeeping genes (63 single nucleotide polymorphisms from around 5000 nucleotides per isolate). We observe a strong negative correlation between within-population genetic diversity and geographic distance from sub-Saharan Africa (R(2) = 0.95) over Africa, Asia, and Oceania. In contrast, regional variation in transmission intensity seems to have had a negligible impact on the distribution of genetic diversity. The striking geographic patterns of isolation by distance observed in P. falciparum mirror the ones previously documented in humans [4-7] and point to a joint sub-Saharan African origin between the parasite and its host. Age estimates for the expansion of P. falciparum further support that anatomically modern humans were infected prior to their exit out of Africa and carried the parasite along during their colonization of the world.

Worldwide sequence conservation of transmission-blocking vaccine candidate Pvs230 in Plasmodium vivax

Pfs230, surface protein of gametocyte/gamete of the human malaria parasite, Plasmodium falciparum, is a prime candidate of malaria transmission-blocking vaccine. Plasmodium vivax has an ortholog of Pfs230 (Pvs230), however, there has been no study in any aspects on Pvs230 to date. To investigate whether Pvs230 can be a vivax malaria transmission-blocking vaccine, we performed evolutionary and population genetic analysis of the Pvs230 gene (pvs230: PVX_003905). Our analysis of Pvs230 and its orthologs in eight Plasmodium species revealed two distinctive parts: an interspecies variable part (IVP) containing species-specific oligopeptide repeats at the N-terminus and a 7.5 kb interspecies conserved part (ICP) containing 14 cysteine-rich domains. Pvs230 was closely related to its orthologs, Pks230 and Pcys230, in monkey malaria parasites. Analysis of 113 pvs230 sequences obtained from worldwide, showed that nucleotide diversity is remarkably low in the non-repeat 8-kb region of pvs230 (theta pi = 0.00118) with 77 polymorphic nucleotide sites, 40 of which results in amino acid replacements. A signature of purifying selection but not of balancing selection was seen on pvs230. Functional and/or structural constraints may limit the level of polymorphism in pvs230. The observed limited polymorphism in pvs230 should ground for utilization of Pvs230 as an effective transmission-blocking vaccine. (C) 2011 Elsevier Ltd. All rights reserved.

Cathepsin L-like genes of Trypanosoma vivax from Africa and South America - characterization, relationships and diagnostic implications

We characterized sequences from genes encoding cathepsin L-like (CatL-like) cysteine proteases from African and South American isolates of Trypanosoma vivax and T. vivax-like organisms, and evaluated their suitability as genetic markers for population structure analysis and diagnosis. Phylogenetic analysis of sequences corresponding to CatL-like catalytic domains revealed substantial polymorphism, and clades of sequences (TviCatL1-9) were separated by large genetic distances. TviCatL1-4 sequences were from cattle isolates from West Africa (Nigeria and Burkina Faso) and South America (Brazil and Venezuela), which belonged to the same T. vivax genotype. T. vivax-like genotypes from East Africa showed divergent sequences, including TviCatL5-7 for isolates from Mozambique and TviCatL8-9 for an isolate from Kenya. Phylogenetic analysis of CatL-like gene data supported the relationships among trypanosome species reflected in the phylogenies based on the analysis of small subunit (SSU) of ribosomal RNA gene sequence data. The discovery of different CatL-like sequences for each genotype, defined previously by ribosomal DNA data, indicate that these sequences provide useful targets for epidemiological and population genetic studies. Regions in CatL-like sequences shared by all T. vivax genotypes but not by other trypanosomes allowed the establishment of a specific and sensitive diagnostic PCR for epidemiological studies in South America and Africa. (C) 2008 Elsevier Ltd. All rights reserved.

Extensive microsatellite diversity in the human malaria parasite Plasmodium vivax

The population structure of Plasmodium vivax remains elusive. The markers of choice for large-scale population genetic studies of eukaryotes, short tandem repeats known as microsatellites, have been recently reported to be less polymorphic in R vivax. Here we investigate the microsatellite diversity and geographic structure in P vivax, at both local and global levels, using 14 new markers consisting of tri- or tetranucleotide repeats. The local-level analysis, which involved 50 field isolates from Sri Lanka, revealed unexpectedly high diversity (average virtual heterozygosity [H-E], 0.807) and significant multilocus linkage disequilibrium in this region of low malaria endemicity. Multiple-clone infections occurred in 60% of isolates sampled in 2005. The global-level analysis of field isolates or monkey-adapted strains identified 150 unique haplotypes among 164 parasites from four continents. Individual P. vivax isolates could not be unambiguously assigned to geographic populations. For example, we found relatively low divergence among parasites from Central America, Africa, Southeast Asia and Oceania, but substantial differentiation between parasites from the same continent (South Asia and Southeast Asia) or even from the same country (Brazil). Parasite relapses, which may extend the duration of P. vivax carriage in humans, are suggested to facilitate the spread of strains across continents, breaking down any pre-existing geographic structure. (C) 2008 Elsevier B.V. All rights reserved.

Restriction fragment length polymorphism (RFLP) in Exon 2 of the BoLA-DRB3 gene in South American cattle

The Bola-DRB3 gene participates in the development of the immune response and is highly polymorphic. For these reasons, it has been a candidate gene in studies of the genetic basis of disease resistance and in population genetic analysis. South American native cattle breeds have been widely replaced by improved exotic breeds leading to a loss of genetic resources. In particular South American native breeds have high levels of fertility and disease resistance. This work describes genetic variability in the BoLA-DRB3 gene in native (Caracu, Pantaneiro, Argentinean Creole) and exotic (Holstein, Jersey, Nelore, Gir) cattle breeds in Brazil and Argentina. PCR-RFLP alleles were identified by combining the restriction patterns for the BoLA-DRB3.2 locus obtained with RsaI, BstY, and HaeIII restriction enzymes. Allelic frequencies and deviations from the Hardy-Weinberg equilibrium were also calculated. Analysis of the 24 BoLA-DRB3 PCR-RFLP alleles identified showed differences in the allele distributions among breeds.

Estimativa de ganho genético a partir de diferenciais de seleção e parâmetros populacionais em um rebanho Caracu

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Diversidade genética em populações-núcleo de Eucalyptus grandis

No melhoramento genético de espécies florestais, uma população base ou indivíduos superiores pré-selecionados tem importância fundamental para a manutenção do programa. Indivíduos de melhores procedências e de ampla base genética propiciam a obtenção de ganhos de forma contínua. O objetivo deste trabalho foi avaliar a diversidade genética em duas populações-núcleo de Eucalyptus grandis. Foram avaliados 39 indivíduos, sendo 19 pertencentes à população 1 e 20, à população 2, utilizando-se 14 primers microssatélite. Os fragmentos foram identificados e analisados a partir dos programas GeneScan e Genotyper, utilizando-se um sequenciador automático ABI Prism 3100. O número de alelos encontrados para cada primer variou de cinco a 15 para a população 1 e, de 8 a 18 para a população 2. A heterozigosidade estimada foi maior na população 2, 0,869, contra 0,843 na população 1. A média da distância genética entre os indivíduos da população 1 foi 0,6220 e na população 2 foi 0,6112. Com a caracterização molecular dos indivíduos destas populações foi construído um banco de dados que permitirá, a partir dos parâmetros de genética de populações, monitorar esses programas de melhoramento em diferentes ciclos de seleção.

Estrutura genética de Chromis multilineata (Guichenot, 1853) (Perciformes Pomacentridae) na costa NE do Brasil e ilhas oceânicas: evidências de um passado evolutivo instável.

The fauna of Brazilian reef fishes comprises approximately 320 species distributed along the coast of the mainland and islands ocean. Little is known about the levels of connectivity between their populations, but has been given the interest in the relations between the offshore and the islands of the Brazil, in a biogeographical perspective. The oceanic islands Brazilian hosting a considerable number of endemic species, which are locally abundant, and divide a substantial portion of its reef fish fauna with the Western Atlantic. Among the richest families of reef fish in species are Pomacentridae. This study analyzed through analysis of sequences of the mitochondrial DNA control region (D-loop), the standards-breeding population of C. Multilineata in different areas of the NE coast of Brazil, involving both oceanic islands (Fernando de Noronha Archipelago and of St. Peter and St. Paul) and continental shelf (RN and BA). To this aim, partial sequences were used in the region HVR1 of mtDNA (312pb). The population structure and parameters for the estimates of genetic variability, molecular variance (AMOVA), estimation of the index for fixing (FST) and number of migrants were determined. The phylogenetic relationships between the populations were estimated using neighbor-joining (NJ) method. A group of Bayesian analysis was used to verify population structure, according to haplotype frequency of each individual. The genetic variability of populations was extremely high. The populations sampled show moderate genetic structure, with a higher degree of genetic divergence being observed for the sample of the Archipelago of St. Peter and St. Paul. At smaller geographical scale, the sample of Rio Grande do Norte and the Archipelago of Fernando de Noronha do not have genetic differentiation. Three moderately differentiated population groups were identified: a population group (I), formed by the Rio Grande do Norte (I') and the archipelago of Fernando de Noronha (I''), and two other different groups formed by the island population of the archipelago of Saint Peter and St. Paul (II) and Bahia (III). The genetic patterns found suggest that the species has suffered a relatively recent radiation favoring the absence of shared haplotypes. C. multilineata seems to constitute a relatively homogenous population along the West Atlantic coast, with evidence of a moderate population genetic structure in relation to the Archipelago of St. Peter and St. Paul. These data supports the importance of the dispersal larvae by marine current and the interpopulation similarity this species.