Temperature-independent genome-wide DNA methylation profile in turbot post-embryonic development.

The morphological and biological characteristics of ectothermic vertebrates are known to be strongly influenced by environmental conditions, particularly temperature. Epigenetic mechanisms such as DNA methylation have been reported to contribute to the phenotypic plasticity observed in vertebrates in response to environmental changes. Additionally, DNA methylation is a dynamic process that occurs throughout vertebrate ontogeny and it has been associated with the activation and silencing of gene expression during post-embryonic development and metamorphosis. In this study, we investigated genome-wide DNA methylation profiles during turbot metamorphosis, as well as the epigenetic effects of temperature on turbot post-embryonic development. Fish growth and rates of development were greatly affected by rearing temperature. Thus, turbot raised at ambient temperature (18 °C) achieved greater body weights and progressed through development more quickly than those reared at a colder temperature (14 °C). Genome-wide DNA methylation dynamics analyzed via a methylation-sensitive amplified polymorphism (MSAP) technique were not significantly different between animals reared within the two different thermal environments. Furthermore, comparisons between phenotypically similar fish revealed that genome-wide DNA methylation profiles do not necessarily correlate with specific developmental stages in turbot.

Population structure and effective/census population size ratio in threatened three-spined stickleback populations from an isolated river basin in northwest Spain.

Variability at 20 microsatellite loci was examined to assess the population genetic structure, gene flow, and effective population size (N(e)) in three populations of three-spined stickleback (Gasterosteus aculeatus) from the upper basin of the Miño River in Galicia, NW Spain, where this species is threatened. The three populations showed similar levels of genetic diversity. There is a significant genetic differentiation between the three populations, but also significant gene flow. N(e) estimates based on linkage disequilibrium yielded values of 355 for the Miño River population and 241 and 311 for the Rato and Guisande Rivers, respectively, although we expect that these are overestimates. N(e) estimates based on temporal methods, considering gene flow or not, for the tributaries yielded values of 30-56 and 47-56 for the Rato and Guisande Rivers, respectively. Estimated census size (N(c)) for the Rato River was 880 individuals. This yielded a N(e)/N(c) estimate of 3-6 % for temporal estimation of N(e), which is within the empirical range observed in freshwater fishes. We suggest that the three populations analyzed have a sufficient level of genetic diversity with some genetic structure. Additionally, the absence of physical barriers suggests that conservation efforts and monitoring should focus in the whole basin as a unit.

Genome-wide methylation study of diploid and triploid brown trout (Salmo trutta L.).

The induction of triploidization in fish is a very common practice in aquaculture. Although triploidization has been applied successfully in many salmonid species, little is known about the epigenetic mechanisms implicated in the maintenance of the normal functions of the new polyploid genome. By means of methylation-sensitive amplified polymorphism (MSAP) techniques, genome-wide methylation changes associated with triploidization were assessed in DNA samples obtained from diploid and triploid siblings of brown trout (Salmo trutta). Simple comparative body measurements showed that the triploid trout used in the study were statistically bigger, however, not heavier than their diploid counterparts. The statistical analysis of the MSAP data showed no significant differences between diploid and triploid brown trout in respect to brain, gill, heart, liver, kidney or muscle samples. Nonetheless, local analysis pointed to the possibility of differences in connection with concrete loci. This is the first study that has investigated DNA methylation alterations associated with triploidization in brown trout. Our results set the basis for new studies to be undertaken and provide a new approach concerning triploidization effects of the salmonid genome while also contributing to the better understanding of the genome-wide methylation processes.

A simulation study on the performance of differentiation-based methods to detect selected loci using linked neutral markers.

We investigated the performance of two of the most popular differentiation-based methods to detect loci under selection (dfdist/fdist and bayescan) in order to ascertain the average chromosome map distance between the detected outlier markers and the nearest loci under selection. We used a model of neutral markers genetically linked to selected loci (QTL) controlling a quantitative trait subject to divergent selection in two subpopulations connected by migration. The results are not particularly encouraging because for chromosome lengths above 0.5 morgan, at least 30% of outliers detected were positioned in chromosomes where QTL were absent, clearly denoting false positives. Outliers linked to QTL were on average closer to the nearest QTL than randomly chosen markers, but the methods showed a substantial uncertainty about the genetic association between markers and selected loci, as this association could be shown significantly only in a moderate number of replicates for most scenarios. At equal conditions, bayescan seemed to perform somewhat more efficiently than dfdist/fdist, with little difference between results for dominant and codominant markers.

Comparing three different methods to detect selective loci using dominant markers.

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.

Lack of nonadditive genetic effects on early fecundity in Drosophila melanogaster.

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.

msap: a tool for the statistical analysis of methylation-sensitive amplified polymorphism data.

In this study msap, an R package which analyses methylation-sensitive amplified polymorphism (MSAP or MS-AFLP) data is presented. The program provides a deep analysis of epigenetic variation starting from a binary data matrix indicating the banding pattern between the isoesquizomeric endonucleases HpaII and MspI, with differential sensitivity to cytosine methylation. After comparing the restriction fragments, the program determines if each fragment is susceptible to methylation (representative of epigenetic variation) or if there is no evidence of methylation (representative of genetic variation). The package provides, in a user-friendly command line interface, a pipeline of different analyses of the variation (genetic and epigenetic) among user-defined groups of samples, as well as the classification of the methylation occurrences in those groups. Statistical testing provides support to the analyses. A comprehensive report of the analyses and several useful plots could help researchers to assess the epigenetic and genetic variation in their MSAP experiments. msap is downloadable from CRAN (http://cran.r-project.org/) and its own webpage (http://msap.r-forge.R-project.org/). The package is intended to be easy to use even for those people unfamiliar with the R command line environment. Advanced users may take advantage of the available source code to adapt msap to more complex analyses.

The evolutionary forces maintaining a wild polymorphism of Littorina saxatilis: model selection by computer simulations.

Two rocky shore ecotypes of Littorina saxatilis from north-west Spain live at different shore levels and habitats and have developed an incomplete reproductive isolation through size assortative mating. The system is regarded as an example of sympatric ecological speciation. Several experiments have indicated that different evolutionary forces (migration, assortative mating and habitat-dependent selection) play a role in maintaining the polymorphism. However, an assessment of the combined contributions of these forces supporting the observed pattern in the wild is absent. A model selection procedure using computer simulations was used to investigate the contribution of the different evolutionary forces towards the maintenance of the polymorphism. The agreement between alternative models and experimental estimates for a number of parameters was quantified by a least square method. The results of the analysis show that the fittest evolutionary model for the observed polymorphism is characterized by a high gene flow, intermediate-high reproductive isolation between ecotypes, and a moderate to strong selection against the nonresident ecotypes on each shore level. In addition, a substantial number of additive loci contributing to the selected trait and a narrow hybrid definition with respect to the phenotype are scenarios that better explain the polymorphism, whereas the ecotype fitnesses at the mid-shore, the level of phenotypic plasticity, and environmental effects are not key parameters.

Genetic differentiation and estimation of effective population size and migration rates in two sympatric ecotypes of the marine snail Littorina saxatilis.

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 (theta(ST) = 0.067) than in SU ones (theta(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.

Mutation-selection balance accounting for genetic variation for viability in Drosophila melanogaster as deduced from an inbreeding and artificial selection experiment.

We carried out an experiment of inbreeding and upward artificial selection for egg-to-adult viability in a recently captured population of Drosophila melanogaster, as well as computer simulations of the experimental design, in order to obtain information on the nature of genetic variation for this important fitness component. The inbreeding depression was linear with a rate of 0.70 +/- 0.11% of the initial mean per 1% increase in inbreeding coefficient, and the realized heritability was 0.06 +/- 0.07. We compared the empirical observations of inbreeding depression and selection response with computer simulations assuming a balance between the occurrence of partially recessive deleterious mutations and their elimination by selection. Our results suggest that a model assuming mutation-selection balance with realistic mutational parameters can explain the genetic variation for viability in the natural population studied. Several mutational models are incompatible with some observations and can be discarded. Mutational models assuming a low rate of mutations of large average effect and highly recessive gene action, and others assuming a high rate of mutations of small average effect and close to additive gene action, are compatible with all the observations.