Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age.

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so-called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral-like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low-dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake-stream divergence in parapatric Swiss stickleback.

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecological speciation. Ecological speciation is well studied in parts of the native range of the three-spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving different stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajectories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known ecological relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.

Population genomic signatures of divergent adaptation, gene flow and hybrid speciation in the rapid radiation of Lake Victoria cichlid fishes.

Adaptive radiations are an important source of biodiversity and are often characterized by many speciation events in very short succession. It has been proposed that the high speciation rates in these radiations may be fuelled by novel genetic combinations produced in episodes of hybridization among the young species. The role of such hybridization events in the evolutionary history of a group can be investigated by comparing the genealogical relationships inferred from different subsets of loci, but such studies have thus far often been hampered by shallow genetic divergences, especially in young adaptive radiations, and the lack of genome-scale molecular data. Here, we use a genome-wide sampling of SNPs identified within restriction site-associated DNA (RAD) tags to investigate the genomic consistency of patterns of shared ancestry and adaptive divergence among five sympatric cichlid species of two genera, Pundamilia and Mbipia, which form part of the massive adaptive radiation of cichlids in the East African Lake Victoria. Species pairs differ along several axes: male nuptial colouration, feeding ecology, depth distribution, as well as the morphological traits that distinguish the two genera and more subtle morphological differences. Using outlier scan approaches, we identify signals of divergent selection between all species pairs with a number of loci showing parallel patterns in replicated contrasts either between genera or between male colour types. We then create SNP subsets that we expect to be characterized to different extents by selection history and neutral processes and describe phylogenetic and population genetic patterns across these subsets. These analyses reveal very different evolutionary histories for different regions of the genome. To explain these results, we propose at least two intergeneric hybridization events (between Mbipia spp. and Pundamilia spp.) in the evolutionary history of these five species that would have lead to the evolution of novel trait combinations and new species.

Population structure and biogeography of migratory freshwater fishes (Prochilodus: Characiformes) in major South American rivers.

Mitochondrial DNA (mtDNA) sequences were used to infer the phylogenetic relationships of Prochilodus species in the Paraná, Amazonas, Orinoco, and Magdalena basins. Sequences of ATPase subunits 6 and 8 (total 840 bp) were obtained for 21 Prochilodus specimens from the four river systems. Using Semaprochilodus as an outgroup, phylogenetic analyses showed that: (i) each river basin contains a monophyletic group of mtDNA lineages; and (ii) the branching order places Magdalena in a basal position with subsequent branching of Orinoco, Amazon and Paraná. The mitochondrial control region was sequenced for 26 P. lineatus (from the Paraná basin) and six other Prochilodus specimens from the Magdalena, Orinoco and Amazon. All 26 control region haplotypes were unique with sequence divergence ranging from 0.3 to 3.6%. The control region phylogeny is well resolved but phylogenetic structure is not associated with geography. For example, mtDNA haplotypes from the upper Paraná (Mogui Guassú) and the upper Bermejo, separated by at least 2600 km, have close genealogical ties. Phylogeographic analyses, including nested clade analysis, suggest high levels of gene flow within this basin.