Biochemistry, Genetics and Inmunology

We introduce a computer program for the dynamic and flexible management of conserved subdivided populations. Using molecular marker data or pedigree information, the software determines the optimal contributions (i.e., number of offspring) of each individual, the number of migrants, and the particular subpopulations involved in the exchange of individuals in order to maintain the largest level of gene diversity in the whole population with a desired control in the rate of inbreeding. Restrictions can be introduced for the total number of migrants, and the mating of particular pairs and their contribution. A full genetic diversity analysis of the population is carried out. The optimal contribution from each subpopulation to a pool of maximal gene diversity is also provided by the program.

Sexual isolation is a key component of reproductive isolation, involving mate choice among mature adults. While there are various statistics for estimating sexual isolation from mating frequencies, their ability to produce unbiased estimates varies considerably, depending on the particular situation. We investigated, under different biological scenarios, the estimation properties (statistical bias, efficiency, root mean square, statistical test) of 12 statistics commonly used in the literature for measuring sexual isolation. Yule's Q, V, YA (and related indices) and IPSI are revealed to be the most efficient, with the smallest biases and root mean square deviations. Yule's Q, YA and IPSI show better estimation properties when using infinite sample sizes, while IPSI is preferable using smaller sample sizes. Other statistics investigated should be avoided, at least within the range of conditions considered. Regarding the parametric test of hypothesis, the best alternative is YA. We discuss the advantages and drawbacks of the various estimators, and propose IPSI as the safest for biological sample sizes. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 85, 307–318.

Optimization of contributions of parents to progeny by minimizing the average coancestry of the progeny is an effective strategy for maintaining genetic diversity in ex situ conservation programs, but its application on the basis of molecular markers has the negative collateral effect of homogenizing the allelic frequencies at each locus. Because one of the objectives of a conservation program is to preserve the genetic composition of the original endangered population, we devised a method in which markers are used to maintain the allele frequency distribution at each locus as closely as possible to that of the native population. Contributions of parents were obtained so as to minimize changes in allele frequency for a set of molecular markers in a population of reduced size. We used computer simulations, under a range of scenarios, to assess the effectiveness of the method in comparison with methods in which contributions of minimum coancestry are sought, either making use of molecular markers or genealogical information. Our simulations indicated that the proposed method effectively maintained the original distribution of allele frequencies, particularly under strong linkage, and maintained acceptable levels of genetic diversity in the population. Nevertheless, contributions of minimum coancestry determined from pedigree information but ignoring the genealogy previous to the conservation program, was the most effective method for maintaining allelic frequencies in realistic situations.

On exposed rocky shores in Galicia (northwest Spain), a striking polymorphism exists between two ecotypes (RB and SU) of Littorina saxatilis that occupy different levels of the intertidal zone and exhibit an incomplete reproductive isolation. The setting has been suggested to represent ongoing sympatric speciation by ecological adaptation of the two ecotypes to their respective habitats. In this article we address whether or not the ecotypes have developed their own population structures in response to the rigors of their corresponding environments and life histories. We analyzed four to five allozymic loci from three surveys of the same sites, spanning a 14-year period. An experimental design including three localities with two transects per locality and three shore levels allowed studying temporal and spatial population structure and estimation of effective population sizes (N e ), neighborhood sizes (N n ), and migration rates (m). Genetic differentiation was significantly lower in RB populations (h ST ¼ 0.067) than in SU ones (h ST ¼ 0.124). Mean estimates of N e , N n , and m did not differ significantly between ecotypes, but local ecotype differences in migration between the two closest localities (larger migration rates in RB than in SU populations) could explain the pattern in population differentiation.

We carried out a simulation study to compare the efficiency of three alternative programs (dfdist, detseld and bayescan) to detect loci under directional selection from genome-wide scans using dominant markers. We also evaluated the efficiency of correcting for multiple testing those methods that use a classical probability approach. Under a wide range of scenarios, we conclude that bayescan appears to be more efficient than the other methods, detecting a usually high percentage of true selective loci as well as less than 1% of outliers (false positives) under a fully neutral model. In addition, the percentage of outliers detected by this software is always correlated with the true percentage of selective loci in the genome. Our results show, nevertheless, that false positives are common even with a combination of methods and multitest correction, suggesting that conclusions obtained from this approach should be taken with extreme caution.

Fecundity is usually considered as a trait closely connected to fitness and is expected to exhibit substantial nonadditive genetic variation and inbreeding depression. However, two independent experiments, using populations of different geographical origin, indicate that early fecundity in Drosophila melanogaster behaves as a typical additive trait of low heritability. The first experiment involved artificial selection in inbred and non-inbred lines, all of them started from a common base population previously maintained in the laboratory for about 35 generations. The realized heritability estimate was 0.151 Ϯ 0.075 and the inbreeding depression was very small and nonsignificant (0.09 Ϯ 0.09% of the non-inbred mean per 1% increase in inbreeding coefficient). With inbreeding, the observed decrease in the within-line additive genetic variance and the corresponding increase of the between-line variance were very close to their expected values for pure additive gene action. This result is at odds with previous studies showing inbreeding depression and, therefore, directional dominance for the same trait and species. All experiments, however, used laboratory populations, and it is possible that the original genetic architecture of the trait in nature was subsequently altered by the joint action of random drift and adaptation to captivity. Thus, we carried out a second experiment, involving inbreeding without artificial selection in a population recently collected from the wild. In this case we obtained, again, a maximum-likelihood heritability estimate of 0.210 Ϯ 0.027 and very little nonsignificant inbreeding depression (0.06 Ϯ 0.12%). The results suggest that, for fitness-component traits, low levels of additive genetic variance are not necessarily associated with large inbreeding depression or high levels of nonadditive genetic variance.

For different fitness mutational models, with epistasis introduced, we simulated the consequences of drift (D scenario) or mutation, selection, and drift (MSD scenario) in populations at the MSD balance subsequently subjected to bottlenecks of size N ¼ 2, 10, 50 during 100 generations. No ''conversion'' of nonadditive into additive variance was observed, all components of the fitness genetic variance initially increasing with the inbreeding coefficient F and subsequently decreasing to zero (D) or to an equilibrium value (MSD). In the D scenario, epistasis had no appreciable effect on inbreeding depression and that on the temporal change of variance components was relevant only for high rates of strong epistatic mutation. In parallel, between-line differentiation in mean fitness accelerated with F and that in additive variance reached a maximum at F $ 0.6-0.7, both processes being intensified by strong epistasis. In the MSD scenario, however, the increase in additive variance was smaller, as it was used by selection to purge inbreeding depression (N $ 10), and selection prevented between-line differentiation. Epistasis, either synergistic or antagonistic (this leading to multiple adaptive peaks), had no appreciable effect on MSD results nor, therefore, on the evolutionary rate of fitness change.

Early detection of population declines is essential to prevent extinctions and to ensure sustainable harvest. We evaluated the performance of two Ne estimators to detect population declines: the two-sample temporal method and a one-sample method based on linkage disequilibrium (LD). We used simulated data representing a wide range of population sizes, sample sizes and number of loci. Both methods usually detect a population decline only one generation after it occurs if Ne drops to less than approximately 100, and 40 microsatellite loci and 50 individuals are sampled. However, the LD method often out performed the temporal method by allowing earlier detection of less severe population declines (Ne approximately 200). Power for early detection increased more rapidly with the number of individuals sampled than with the number of loci genotyped, primarily for the LD method. The number of samples available is therefore an important criterion when choosing between the LD and temporal methods. We provide guidelines regarding design of studies targeted at monitoring for population declines. We also report that 40 single nucleotide polymorphism (SNP) markers give slightly lower precision than 10 microsatellite markers. Our results suggest that conservation management and monitoring strategies can reliably use genetic based methods for early detection of population declines.

The study of hybrid zones resulting from Pleistocene vicariance is central in examining the potential of genetically diverged evolutionary units either to introgress and merge or to proceed with further isolation. The hybrid zone between two mitochondrial lineages of Chioglossa lusitanica is located near the Mondego River in Central Portugal. We used mitochondrial and nuclear diagnostic markers to conduct a formal statistical analysis of the Chioglossa hybrid zone in the context of tension zone theory. Key results are: (i) cline centres are not coincident for all markers, with average widths of ca . 2-15 km; (ii) heterozygote deficit was not observed across loci near the transect centre; (iii) associations of parental allele combinations ('linkage disequilibrium' R ) were not detected either across loci or across the transect. These observations suggest that the Chioglossa hybrid zone is not a tension zone with strong selection against hybrids but instead one shaped mostly by neutral mixing. The patterns uncovered suggest a complex history of populations over a small scale that may be common in southern Pleistocene refugia.

Western Indian Ocean islands (excluding Madagascar) are inhabited by at least five species of the lizard genus Mabuya, for which diversity and inter-island phylogenetic relationships are still unclear, particularly concerning island populations of M. maculilabris, M. comorensis and M. striata. With a comprehensive sampling across the Comoros archipelago and the islands of Zanzibar and Pemba (Tanzania) we use molecular data and previously published sequences of these species to describe genetic variation across their insular range and explore the possible distinction of insular populations of the different species, contributing to the clarification of their status and knowledge of their colonization patterns. We describe patterns of cryptic diversity within M. comorensis, revealing a West to East colonization of the Comoros archipelago and clarify the status of Malagasy populations of this species. Cryptic variation is also observed within M. maculilabris and M. striata, and the status of several continental and island populations discussed. M. sechellensis is shown to be paraphyletic relative to M. wrightii. Available data for all Afro-Malagasy Mabuya reveals additional possible instances of paraphyly within some species but does not allow for further resolution of basal relationships within this group, pointing to rapid and simultaneous divergences within this group and the need for multiple independent markers to recover the history of this group.