Humberto Quesada

PHYLOGENETIC EVIDENCE FOR MULTIPLE SYMPATRIC ECOLOGICAL DIVERSIFICATION IN A MARINE SNAIL

Abstract Parallel speciation can occur when traits determining reproductive isolation evolve independently in different populations that experience a similar range of environments. However, a common problem in studies of parallel evolution is to distinguish this hypothesis from an alternative one in which different ecotypes arose only once in allopatry and now share a sympatric scenario with substantial gene flow between them. Here we show that the combination of a phylogenetic approach with life-history data is able to disentangle both hypotheses in the case of the intertidal marine snail Littorina saxatilis on the rocky shores of Galicia in northwestern Spain. In this system, numerous phenotypic and genetic differences have evolved between two sympatric ecotypes spanning a sharp ecological gradient, and as aside effect of the former have produced partial reproductive isolation. A mitochondrial phylogeny of these populations strongly suggests that the two sympatric ecotypes have originated independently several times. Building upon earlier work demonstrating size-based assortative mating as the main contributor to reproductive isolation among ecotypes, our analysis provides strong evidence that divergent selection across a sharp ecological gradient promoted the parallel divergence of body size and shape between two sympatric ecotypes. Thus, divergent selection occurring independently in different populations has produced the marine equivalent of host races, which may represent the first step in speciation.

Parallel evolution of local adaptation and reproductive isolation in the face of gene flow.

Parallel evolution of similar phenotypes provides strong evidence for the operation of natural selection. Where these phenotypes contribute to reproductive isolation, they further support a role for divergent, habitat‐associated selection in speciation. However, the observation of pairs of divergent ecotypes currently occupying contrasting habitats in distinct geographical regions is not sufficient to infer parallel origins. Here we show striking parallel phenotypic divergence between populations of the rocky‐shore gastropod, Littorina saxatilis, occupying contrasting habitats exposed to either wave action or crab predation. This divergence is associated with barriers to gene exchange but, nevertheless, genetic variation is more strongly structured by geography than by ecotype. Using approximate Bayesian analysis of sequence data and amplified fragment length polymorphism markers, we show that the ecotypes are likely to have arisen in the face of continuous gene flow and that the demographic separation of ecotypes has occurred in parallel at both regional and local scales. Parameter estimates suggest a long delay between colonization of a locality and ecotype formation, perhaps because the postglacial spread of crab populations was slower than the spread of snails. Adaptive differentiation may not be fully genetically independent despite being demographically parallel. These results provide new insight into a major model of ecologically driven speciation.

Evaluating the relationship between evolutionary divergence and phylogenetic accuracy in AFLP data sets

Using in silico amplified fragment length polymorphism (AFLP) fingerprints, we explore the relationship between sequence similarity and phylogeny accuracy to test when, in terms of genetic divergence, the quality of AFLP data becomes too low to be informative for a reliable phylogenetic reconstruction. We generated DNA sequences with known phylogenies using balanced and unbalanced trees with recent, uniform and ancient radiations, and average branch lengths (from the most internal node to the tip) ranging from 0.02 to 0.4 substitutions per site. The resulting sequences were used to emulate the AFLP procedure. Trees were estimated by maximum parsimony (MP), neighbor-joining (NJ), and minimum evolution (ME) methods from both DNA sequences and virtual AFLP fingerprints. The estimated trees were compared with the reference trees using a score that measures overall differences in both topology and relative branch length. As expected, the accuracy of AFLP-based phylogenies decreased dramatically in the more divergent data sets. Above a divergence of approximately 0.05, AFLP-based phylogenies were largely inaccurate irrespective of the distinct topology, radiation model, or phylogenetic method used. This value represents an upper bound of expected tree accuracy for data sets with a simple divergence history; AFLP data sets with a similar divergence but with unbalanced topologies and short ancestral branches produced much less accurate trees. The lack of homology of AFLP bands quickly increases with divergence and reaches its maximum value (100%) at a divergence of only 0.4. Low guanine-cytosine (GC) contents increase the number of nonhomologous bands in AFLP data sets and lead to less reliable trees. However, the effect of the lack of band homology on tree accuracy is surprisingly small relative to the negative impact due to the low information content of AFLP characters. Tree-building methods based on genetic distance displayed similar trends and outperformed parsimony at low but not at high divergences. However, the impact of using alternative phylogenetic methods on tree accuracy was generally small relative to the uncertainty arising from factors such as divergence, nonhomology of bands, or the low information content of AFLP characters. Nevertheless, our data suggest that under certain circumstances, AFLPs may be suitable to reconstruct deeper phylogenies than usually accepted.

Homoplasy and Distribution of AFLP Fragments: An Analysis In Silico of the Genome of Different Species

We carried out an in silico analysis of the complete genome sequences of 14 species, including eukaryotes, prokaryotes, and archaea, to investigate the proportion of amplified fragment length polymorphism bands that are homoplasious for the different species, as well as the distribution of fragment lengths. We investigated several possible reasons for the disagreement, previously observed in Arabidopsis thaliana, between the observed fragment length distribution and the null random sequence distribution, which occurs in the direction of a deficit of fragments of small length and an excess of those of large length with respect to the null distribution. We made the following findings: 1) The positive relationship previously found between the percentage of homoplasy and genome size is a direct consequence of the number of observed bands and the GC content. For the same number of observed bands, the percentage of homoplasy is independent of the genome size of the species. 2) The disagreement between the observed fragment length distribution and the null random sequence distribution observed in A. thaliana is a phenomenon that also occurs in other species. 3) This disagreement is due neither to the structure of the genomes in isochores nor the possible impact of indels in reducing the number of restriction sites, two hypotheses discussed in the literature. 4) Nor is the disagreement eliminated by using restriction enzymes with balanced motifs. 5) The discrepancy seems to be caused, rather, by the nonrandom distribution of restriction enzyme motifs.

Allelic diversity for neutral markers retains a higher adaptive potential for quantitative traits than expected heterozygosity

The adaptive potential of a population depends on the amount of additive genetic variance for quantitative traits of evolutionary importance. This variance is a direct function of the expected frequency of heterozygotes for the loci which affect the trait (QTL). It has been argued, but not demonstrated experimentally, that long‐term response to selection is more dependent on QTL allelic diversity than on QTL heterozygosity. Conservation programmes, aimed at preserving this variation, usually rely on neutral markers rather than on quantitative traits for making decisions on management. Here, we address, both through simulation analyses and experimental studies with Drosophila melanogaster, the question of whether allelic diversity for neutral markers is a better indicator of a high adaptive potential than expected heterozygosity. In both experimental and simulation studies, we established synthetic populations for which either heterozygosity or allelic diversity was maximized using information from QTL (simulations) or unlinked neutral markers (simulations and experiment). The synthetic populations were selected for the quantitative trait to evaluate the evolutionary potential provided by the two optimization methods. Our results show that maximizing the number of alleles of a low number of markers implies higher responses to selection than maximizing their heterozygosity.

Gene-Expression Changes Caused by Inbreeding Protect Against Inbreeding Depression in Drosophila

We present a transcriptomic analysis aimed at investigating whether the changes in gene expression that occur under inbreeding generally reduce or enhance inbreeding depression. Discerning between these two alternatives can be addressed only when both changes in expression due to inbreeding and to inbreeding depression are estimated simultaneously. We used Affymetrix 2.0 arrays to study the changes in gene expression associated with both inbreeding and inbreeding depression for fitness in four sets of inbred sublines of Drosophila melanogaster. We found that for most genes showing changes in expression associated with inbreeding, the least depressed sublines were those showing the largest departures in expression from that of the outbred control. This suggests a pattern consistent with a protective role of expression changes against inbreeding effects, and would reveal a new dimension of the transcriptomics of inbreeding. The variation in depression observed could then be due not only to the genetic damages primarily originating that depression, but also possibly to differences in the ability to carry out the appropriate adjustments in gene expression to cope with the inbreeding. We also found that these expression changes with a putative protective role against inbreeding effects show a clear specificity on RNA synthesis and splicing and energy derivation functions.

Insights into the role of differential gene expression on the ecological adaptation of the snail Littorina saxatilis

BackgroundIn the past 40 years, there has been increasing acceptance that variation in levels of gene expression represents a major source of evolutionary novelty. Gene expression divergence is therefore likely to be involved in the emergence of incipient species, namely, in a context of adaptive radiation. In this study, a genome-wide expression profiling approach (cDNA-AFLP), validated by quantitative real-time polymerase chain reaction (qPCR) were used to get insights into the role of differential gene expression on the ecological adaptation of the marine snail Littorina saxatilis. This gastropod displays two sympatric ecotypes (RB and SU) which are becoming one of the best studied systems for ecological speciation.ResultsAmong the 99 transcripts shared between ecotypes, 12.12% showed significant differential expression. At least 4% of these transcripts still displayed significant differences after correction for multiple tests, highlighting that gene expression can differ considerably between subpopulations adapted to alternative habitats in the face of gene flow. One of the transcripts identified was Cytochrome c Oxidase subunit I (COI). In addition, 6 possible reference genes were validated to normalize and confirm this result using qPCR. α-Tubulin and histone H3.3 showed the more stable expression levels, being therefore chosen as the best option for normalization. The qPCR analysis confirmed a higher COI expression in SU individuals.ConclusionsAt least 4% of the transcriptome studied is being differentially expressed between ecotypes living in alternative habitats, even when gene flow is still substantial between ecotypes. We could identify a candidate transcript of such ecotype differentiation: Cytochrome c Oxidase Subunit I (COI), a mitochondrial gene involved in energy metabolism. Quantitative PCR was used to confirm the differences found in COI and its over-expression in the SU ecotype. Interestingly, COI is involved in the oxidative phosphorylation, suggesting an enhanced mitochondrial gene expression (or increased number of mitochondria) to improve energy supply in the ecotype subjected to the strongest wave action.

The relative contribution of band number to phylogenetic accuracy in AFLP data sets

We examined the effect of increasing the number of sampled amplified fragment length polymorphism (AFLP) bands to reconstruct an accurate and well‐supported AFLP‐based phylogeny. In silico AFLP was performed using simulated DNA sequences evolving along balanced and unbalanced model trees with recent, uniform and ancient radiations and average branch lengths (from the most internal node to the tip) ranging from 0.02 to 0.05 substitutions per site. Trees were estimated by minimum evolution (ME) and maximum parsimony (MP) methods from both DNA sequences and virtual AFLP fingerprints. The comparison of the true tree with the estimated AFLP trees suggests that moderate numbers of AFLP bands are necessary to recover the correct topology with high bootstrap support values (i.e. >70%). Fewer numbers of bands are necessary for shorter tree lengths and for balanced than for unbalanced tree topologies. However, branch length estimation was rather unreliable and did not improve substantially after a certain number of bands were sampled. These results hold for different levels of genome coverage and number of taxa analysed. In silico AFLP using bacterial genomic DNA sequences recovered a well‐supported tree topology that mirrored an empirical phylogeny based on a set of 31 orthologous gene sequences when as few as 263 AFLP bands were scored. These results suggest that AFLPs may be an efficient alternative to traditional DNA sequencing for accurate topology reconstruction of shallow trees when not very short ancestral branches exist.

Historical Biogeography of the Marine Snail Littorina saxatilis Inferred from Haplotype and Shell Morphology Evolution in NW Spain.

The marine snail Littorina saxatilis exhibits extreme morphological variation between and within geographical regions and represents an excellent model for assessing local adaptation. Previous studies support the hypothesis of parallel evolution in sympatry of two morphologically different ecotypes (named as RB and SU) that co-inhabit different habitats from Galician rocky shores (NW Spain), and which are interrupted by sheltered areas inhabited by a different morph never studied before (named as SRB). Here, we use morphological and mitochondrial DNA (mtDNA) sequence data to test hypotheses on the origin and diversification of SRB snails and to assess their evolutionary relationships with RB and SU ecotypes. Our results show that the SRB morph displays the largest size and shell elongation and the smallest relative shell aperture, representing an extreme type of the RB vs. SU polymorphism, which has been linked to adaptation to sheltered ecological factors. Phylogenetic analysis shows that the SRB morph shares ancestry with RB and SU ecotypes, rejecting the hypothesis that the SRB morph marks relict populations from which these ecotypes evolved in Galician coasts. Our data support that genetic differentiation among SRB, RB and SU morphs results from a general pattern of restricted gene flow and isolation by distance linked to the colonization of Galician coasts by two independent mtDNA lineages, rather than from a random fragmentation of the initial distributional range. Therefore, the confinement of distinct lineages to specific geographical areas denote evident limits to the distances these snails can disperse. Morphological analysis indicates no association between mtDNA lineage and a specific morphotype, and suggests the independent gain of convergent morphological patterns within each mtDNA lineage in populations occupying contrasting habitats following the colonization of Galician coasts.

Impact of deep coalescence and recombination on the estimation of phylogenetic relationships among species using AFLP markers

Deep coalescence and the nongenealogical pattern of descent caused by recombination have emerged as a common problem for phylogenetic inference at the species level. Here we use computer simulations to assess whether AFLP-based phylogenies are robust to the uncertainties introduced by these factors. Our results indicate that phylogenetic signal can prevail even in the face of extensive deep coalescence allowing recovering the correct species tree topology. The impact of recombination on tree accuracy was related to total tree depth and species effective population size. The correct tree topology could be recovered upon many simulation settings due to a trade-off between the conflicting signals resulting from intra-locus recombination and the benefits of the joint consideration of unlinked loci that better matched overall the true species tree. Errors in tree topology were not only determined by deep coalescence, but also by the timing of divergence and the tree-building errors arising from an insufficient number of characters. DNA sequences generally outperformed AFLPs upon any simulated scenario, but this difference in performance was nearly negligible when a sufficient number of AFLP characters were sampled. Our simulations suggest that the impact of deep coalescence and intra-locus recombination on the reliability of AFLP trees could be minimal for effective population sizes equal to or lower than 10,000 (typical of many vertebrates and tree plants) given tree depths above 0.02 substitutions per site.