Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel PLOS GENETICS

No abstract

Sequence variation in the mtDNA control region of Arctic charr Salvelinus alpinus and Dolly Varden Salvelinus malma from 56 Siberian and North American populations was analysed to assess their phylogeographic relationships and the origins of sympatric forms. Phylogenetic trees confirm the integrity of phylogroups reported in previous mtDNA studies except that the Siberian group does not separate as a single cluster. Haplotype network analysis indicates the proximity of Siberian and Atlantic haplotypes. These are considered as one Eurasian group represented by the Atlantic, east Siberian (interior Siberia including Transbaikalia, Taimyr) and Eurosiberian (Finland, Spitsbergen, Taimyr) sub-groups. Salvelinus alpinus with presumably introgressed Bering group (malma) haplotypes were found along eastern Siberian coasts up to the Olenek Bay and the Lena Delta region, where they overlap with the Eurasian group and in the easternmost interior region. It is proposed that Siberia was colonized by S. alpinus in two stages: from the west by the Eurasian group and later from the east by the Bering group. The high diversity of Eurasian group haplotypes in Siberia indicates its earlier colonization by S. alpinus as compared with the European Alps. This colonization was rapid, proceeded from a diverse gene pool, and was followed by differential survival of ancestral mtDNA lineages in different basins and regions, and local mutational events in isolated populations. The results presented here support a northern origin of Transbaikalian S. alpinus, the dispersion of S. alpinus to the Lake Baikal Basin from the Lena Basin, segregation of S. alpinus between Lena tributaries and their restricted migration over the divides between sub-basins. These results also support sympatric origin of intralacustrine forms of S. alpinus.

Abstract. We use a comparative phylogeographical framework to evaluate the hypothesis of hydrological exchange during the Pleistocene among the now disjunct Lena, Amur, and Enisei basins in Siberia, and to provide evidence on the causal mechanism of their present day faunal dissimilarities. Approximately 600 bases of the mitochondrial control region were sequenced in five distinct lineages among three genera of salmonid fishes, Hucho, Brachymystax and Thymallus. All three basins were fixed for divergent (2–5.4%) lineages of Thymallus whereas a single shared haplotype was present in all three basins for Hucho taimen (Pallas, 1773) and one shared haplotype between the Lena and Amur basins out of a total of five for blunt‐snouted and one out of five for sharp‐snouted Brachymystax lenok (Pallas, 1773). For both blunt‐ and sharp‐snouted lenok the haplotypes found within each basin did not form clades, so no relationship between genotypes and geographical occurrence was found. Our data support relatively recent hydrological mixing of the major river drainage systems in eastern and far‐eastern Siberia, congruent with the hypothesis of large‐scale palaeo‐hydrological exchange stemming from glacial advance, retreat and melting during Pleistocene climate fluctuations. Furthermore, these results in conjunction with a comparison of overall faunal composition suggest that environmental differences rather than historical contingency may be responsible for the faunal dissimilarities of the Amur, Lena, and Enisei river basins.

Hucho taimen from eight populations spanning four drainage basins (Amur, Lena, Enisei and Khatanga) were analysed for nucleotide sequence variation across three mitochondrial genes (ATP6, NADH-1 and control region). Samples of H. hucho, Brachymystax lenok (sharp-snouted and blunt-snouted forms) and Parahucho perryi were also included for comparison. Nucleotide variation across a total of 1826 base pairs in H. taimen revealed shared haplotypes between the Amur and Lena basins, further supporting a previous hypothesis of late to post-Pleistocene hydrological exchange between these now disjunct basins. In contrast to an earlier study using the control region alone, clear phylogeographic structure was seen at a large geographic scale, reflected by two phylogroups, one corresponding to the Amur and Lena basins, and the other to the Enisei and Khatanga basins. Comparative rates of divergence revealed considerably faster and less heterogeneous substitution rates for the two coding genes, especially at interspecific levels compared to the mtDNA control region. # 2005 The Fisheries Society of the British Isles