Sequence analysis of the mitochondrial DNA control region of ciscoes (genus Coregonus): taxonomic implications for the Great Lakes species flock

Genotyping-by-sequencing illuminates high levels of divergence among sympatric forms of coregonines in the Laurentian Great Lakes

Effective resource management depends on our ability to partition diversity into biologically meaningful units. Recent evolutionary divergence, however, can often lead to ambiguity in morphological and genetic differentiation, complicating the delineation of valid conservation units. Such is the case with the "coregonine problem," where recent postglacial radiations of coregonines into lacustrine habitats resulted in the evolution of numerous species flocks, often with ambiguous taxonomy. The application of genomics methods is beginning to shed light on this problem and the evolutionary mechanisms underlying divergence in these ecologically and economically important fishes. Here, we used restriction site-associated DNA (RAD) sequencing to examine genetic diversity and differentiation among sympatric forms in the Coregonus artedi complex in the Apostle Islands of Lake Superior, the largest lake in the Laurentian Great Lakes. Using 29,068 SNPs, we were able to clearly distinguish among the three most common forms for the first time, as well as identify putative hybrids and potentially misidentified specimens. Population assignment rates for these forms using our RAD data were 93%-100% with the only mis-assignments arising from putative hybrids, an improvement from 62% to 77% using microsatellites. Estimates of pairwise differentiation (F ST: 0.045-0.056) were large given the detection of hybrids, suggesting that reduced fitness of hybrid individuals may be a potential mechanism for the maintenance of differentiation. We also used a newly built C. artedi linkage map to look for islands of genetic divergence among forms and found widespread differentiation across the genome, a pattern indicative of long-term drift, suggesting that these forms have been reproductively isolated for a substantial amount of time. The results of this study provide valuable information that can be applied to develop well-informed management strategies and stress the importance of re-evaluating conservation units with genomic tools to ensure they accurately reflect species diversity.

SPECIES FLOCK IN THE NORTH AMERICAN GREAT LAKES: MOLECULAR ECOLOGY OF LAKE NIPIGON CISCOES (TELEOSTEI: COREGONIDAE: COREGONUS)

Studies on north temperate fish species indicate that new habitat availability following the last ice sheet retreat has promoted ecological speciation in postglacial lakes. Extensive ecophenotypic polymorphisms observed among the North American Great Lakes ciscoes suggest that this fish group has radiated through trophic adaptation and reproductive isolation. This study aims at relating the ecomorphological and genetic polymorphisms expressed by the Lake Nipigon ciscoes to evaluate the likelihood of an intralacustrine divergence driven by the exploitation of alternative resources. Morphological variation and trophic and spatial niches are characterized and contrasted among 203 individuals. Genetic variation at six microsatellite loci is also analyzed to appraise the extent of genetic differentiation among these morphotypes. Ecomorphological data confirm the existence of four distinct morphotypes displaying various levels of trophic and depth niche overlap and specialization. However, ecological and morphological variations were not coupled as expected, suggesting that trophic morphology is not always predictive of ecology. Although extensive genetic variability was observed, little genetic differentiation was found among morphotypes, with only one morph being slightly but significantly differentiated. Contrasting patterns of morphological, ecological, and genetic polymorphisms did not support the hypothesis of ecological speciation: the most ecologically different forms were morphologically most similar, while the only genetically differentiated morph was the least ecologically specialized. The low levels of genetic differentiation and the congruence between θ and φ estimates altogether suggest a recent (most likely postglacial) process of divergence and/or high gene flow among morphs A, C, and D, whereas higher φ estimates for comparison involving morph B suggest that this morph may be derived from another colonizing lineage exchanging little genes with the other morphs. Patterns of ecophenotypic and genetic diversity are also compatible with a more complex evolutionary history involving hybridization and introgression.

Genetic structure of the cisco (Coregonus albula L.) from lakes of glacial origin in northern Poland

European cisco (Coregonus albula L.) is one of the most precious species of the European lake ichtyofauna, however, due to progressive eutrophication of water, the range of its occurrence has decreased. Deteriorating ecological conditions are the main cause of this decline in population, and most of the existing populations in lakes are maintained thanks to reintroduction. Thus, it is important to determine the genetic structure of the European cisco. The study involved PCR-RFLP-based genetic analysis of C. albula caught in 15 lakes in northern Poland, including four lakes located in national parks. The analysis covered 3 genes located in the mitochondrial DNA: ND1, ND3/4 and cytochrome b, as well as a control region (D-loop). The PCR product was digested with 4 endonucleases (RsaI, MspI, BsuI and HhaI) and the resulting haplotypes were grouped into combinations. Statistical analysis were then performed on these groups. Based on the genetic distance, a phenogram was constructed in which two groups could be distinguished. One group was represented by C. albula populations from most lakes in north-western Poland, including the three protected lakes. The other group consisted of the European cisco population from Lake Wigry (north-eastern Poland) and commercially exploited lakes from the southern part of the investigated north-western Poland area. The results of the study of C. albula from northern Poland present a valuable molecular characterization of the populations and can be a starting point for further genetic monitoring.

Mitochondrial DNA identification of game and harvested freshwater fish species

The use of DNA in forensics has grown rapidly for human applications along with the concomitant development of bioinformatics and demographic databases to help fully realize the potential of this molecular information. Similar techniques are also used routinely in many wildlife cases, such as species identification in food products, poaching and the illegal trade of endangered species. The use of molecular techniques in forensic cases related to wildlife and the development of associated databases has, however, mainly focused on large mammals with the exception of a few high-profile species. There is a need to develop similar databases for aquatic species for fisheries enforcement, given the large number of exploited and endangered fish species, the intensity of exploitation, and challenges in identifying species and their derived products. We sequenced a 500bp fragment of the mitochondrial cytochrome b gene from representative individuals from 26 harvested fish taxa from Ontario, Canada, focusing on species that support major commercial and recreational fisheries. Ontario provides a unique model system for the development of a fish species database, as the province contains an evolutionarily diverse array of freshwater fish families representing more than one third of all freshwater fish in Canada. Inter- and intraspecific sequence comparisons using phylogenetic analysis and a BLAST search algorithm provided rigorous statistical metrics for species identification. This methodology and these data will aid in fisheries enforcement, providing a tool to easily and accurately identify fish species in enforcement investigations that would have otherwise been difficult or impossible to pursue.

Mitochondrial DNA phylogeography of lake cisco (Coregonus artedi): evidence supporting extensive secondary contacts between two glacial races

The comparative molecular phylogeography of regional fish fauna has revealed the wide distribution of young clades in freshwater fishes of formerly glaciated areas as well as interspecific incongruences in their refugial origins and recolonization routes. In this study, we employed single-strand conformation polymorphism (SSCP) and sequence analyses to describe mitochondrial DNA (mtDNA) polymorphism among 27 populations of the lake cisco (Coregonus artedi) from its entire range of distribution in order to evaluate the hypothesis of dual glacial refuges proposed by Bernatchez & Dodson against the traditional view that this species is solely of Mississippian origin. Results indicate that this taxon is composed of two closely related groups that are widely distributed and intermixed over most of the sampled range. The estimated level of divergence (0.48%), the contrast in the geographical distribution of each group, as well as the general distribution of C. artedi in North America together support the hypothesis that one group dispersed from a Mississippian refuge via the proglacial lakes, while the other is of Atlantic origin and also took advantages of earlier dispersal routes towards eastern Hudson Bay drainages. However, the signal of past range fragmentation revealed by a nested clade analysis was weak, and did not allow to formally exclude the hypothesis of a single Mississippian origin for both lineages. Comparisons with the phylogeographic patterns of other Nearctic freshwater fishes suggest that the salinity tolerance and thermal sensitivity of lake cisco may have been determinant for its extensive postglacial dispersal. The presence or co-occurrence of sympatric or allopatric eco/morphotypes were not found to be necessarily associated with the presence of both haplotype groups.