Lake Cisco Research Articles

Historical and morphological evidence for a remnant population of Lake Erie cisco Coregonus artedi (albus) in Crystal Lake, Pennsylvania

The cisco (Coregonus artedi) population in Crystal Lake, Pennsylvania, is of great scientific interest as it either originated from Lake Erie or Lake Ontario. Cisco in Lake Erie once supported the largest freshwater fishery in the world, but populations were extirpated by 1960. We conducted a morphological analysis of Crystal Lake cisco to determine whether it was consistent with a distinctive Lake Erie form (albus), which was also historically documented, albeit rarely, in western Lake Ontario. Using principal component analysis, we compared eight morphometric ratios and one meristic from our Crystal Lake cisco collection with historical and contemporary collections of cisco from Lakes Erie and Ontario. Maximum likelihood ellipse overlaps between Crystal Lake cisco and presumed albus (the dominant Lake Erie form prior to extirpation) collections averaged 54%. For all groups, the greatest morphological overlap (73.9%) occurred between Crystal Lake and 1957 Lake Erie cisco, which only differed from Crystal Lake cisco in dorsal fin length. Alternatively, overlap between Crystal Lake cisco and all other Lake Ontario collections averaged 3.2%. Our results demonstrate that Crystal Lake cisco are likely an albus form; furthermore, historical documentation and our morphological results suggest a Lake Erie origin. Substantial overlap between Crystal Lake cisco and Lake Ontario albus collected in 1917 is likely explained by continuous entrainment of Lake Erie larvae into Lake Ontario. We suspect this created an albus metapopulation spanning Lakes Erie and Ontario, yet albus are no longer observed in either lake today.

Body size and condition, age and catch‐and‐release injury of ice‐angled cisco (Coregonus artedi) from four connected lakes in Saskatchewan, Canada

AbstractThe lake cisco (Coregonus artedi) is not often caught by recreational anglers in open‐water in western Canada but becomes a common component of catch‐and‐release ice fisheries in the winter. Little is known about the composition of recreational cisco fisheries or factors that affect the success of winter catch and release. We examined the size and body condition of 555 cisco caught by anglers in four interconnected lakes in southern Saskatchewan, Canada. We also examined the influence of length on the probability of hooking injuries and determined the age of a subset of angled fish. Cisco caught by anglers ranged in total length from 24 to 48 cm and 110 to 1370 g mass; however, 60.2% of cisco caught ranged from 37 to 40 cm and 470 to 560 g. The scale ages of cisco caught by anglers ranged from 4 to 10 years, with the most common ages caught being 6–8 years; however, scales likely underestimated age beginning around 8 years. Overall, cisco from our four study lakes had lower body conditions compared to Lake Superior populations based on standard weight comparisons. The length, mass, age and body condition of cisco caught by anglers varied significantly among our four study lakes, following a west‐to‐east productivity gradient. Of the 555 cisco assessed, 54 (9.7%) experienced gill bleeding and 21 (3.8%) had other visible hooking injuries. Longer cisco had a significantly increased probability of gill bleeding, but length did not affect the probability of other injuries. Cisco are relatively slow growing in our study area and warrant more consideration from a catch‐and‐release angling perspective.

Variability in the drivers of microplastic consumption by fish across four lake ecosystems

Microplastic (<5 mm) pollution has been documented globally throughout freshwater and marine ecosystems. Exposure to and ingestion of microplastics presents a threat to the health of aquatic and marine organisms and humans through the consumption of fish and crustaceans. Understanding the factors which influence microplastic ingestion by fish is a key step in predicting the potential health risks. Drivers of microplastic consumption have been studied in lab settings, but there has been limited ability to confirm in field studies. Here we examine the roles and contributions of feeding guild, pelagic microplastic concentrations, and fish length to microplastic consumption by three species of fish across four lake ecosystems in Minnesota, United States. Fish samples were collected in the summers of 2019 and 2020 and processed to determine variability in microplastic ingestion. Identifying particles between 0.18 and 5 mm, plastic ingestion ranged from 0.6 microplastics fish−1 in Elk Lake (low surface water microplastics) bluegill to 1.09 microplastics fish-1 in White Iron Lake (moderate surface water microplastics) cisco. Results indicate that microplastic consumption by filter feeding cisco is driven by surface water microplastic concentrations, while microplastic consumption by visual feeding bluegill and yellow perch is not. Additionally, the high variability of ingestion between lake ecosystems coupled with the complex behavior of some fish species presents difficulties in identifying primary drivers of microplastic consumption that would be broadly applicable across ecosystems and species.

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range.

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well‐oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site‐associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation‐by‐distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non‐neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

Associations between omega-3 fatty acids, selenium content, and mercury levels in wild-harvested fish from the Dehcho Region, Northwest Territories, Canada

ABSTRACTTo better understand the risks and benefits of eating wild-harvested fish from the Northwest Territories, Canada, levels of total mercury (HgT) and selenium (Se) and composition of omega-3 fatty acid (n-3 FA) were measured in muscle tissue of fish harvested from lakes in the Dehcho Region, Northwest Territories, Canada. Average HgT levels ranged from 0.057 mg/kg (cisco) to 0.551 mg/kg (northern pike), while average n-3 FA concentrations ranged from 101 mg/100 g (burbot) to 1,689 mg/100 g (lake trout). In contrast to HgT and n-3 FA, mean Se concentrations were relatively similar among species. Consequently, species such as lake whitefish, cisco, and longnose sucker displayed the highest nutrient levels relative to HgT content. Levels of HgT tended to increase with fish size, while Se and n-3 FA levels were typically not associated with fork length or fish weight. Interestingly, HgT concentration was occasionally inversely related to tissue nutrient content. Significant negative correlations were observed between Hg and n-3 FA for lake trout, northern pike, and walleye. There were also significant negative correlations between Hg and Se noted for lake whitefish, cisco, and northern pike. Samples with the highest nutritional content displayed, on occasion, lower levels of HgT. This study provides valuable information for the design of probabilistic models capable of refining public health messaging related to minimizing Hg risks and maximizing nutrient levels in wild-harvested fish in the Canadian subarctic.

Spatial synchrony in cisco recruitment

We examined the spatial scale of recruitment variability for disparate cisco (Coregonus artedi) populations in the Great Lakes (n=8) and Minnesota inland lakes (n=4). We found that the scale of synchrony was approximately 400km when all available data were utilized; much greater than the 50-km scale suggested for freshwater fish populations in an earlier global analysis. The presence of recruitment synchrony between Great Lakes and inland lake cisco populations supports the hypothesis that synchronicity is driven by climate and not dispersal. We also found synchrony in larval densities among three Lake Superior populations separated by 25–275km, which further supports the hypothesis that broad-scale climatic factors are the cause of spatial synchrony. Among several candidate climate variables measured during the period of larval cisco emergence, maximum wind speeds exhibited the most similar spatial scale of synchrony to that observed for cisco. Other factors, such as average water temperatures, exhibited synchrony on broader spatial scales, which suggests they could also be contributing to recruitment synchrony. Our results provide evidence that abiotic factors can induce synchronous patterns of recruitment for populations of cisco inhabiting waters across a broad geographic range, and show that broad-scale synchrony of recruitment can occur in freshwater fish populations as well as those from marine systems.

Morphology and life history of the Great Slave Lake ciscoes (Salmoniformes: Coregonidae)

AbstractThe taxonomy of the North American ciscoes (Salmoniformes: Coregonidae) remains unresolved. We provide the first comprehensive description of the Great Slave Lake ciscoes. Our analysis supports the hypothesis that the Great Slave Lake cisco complex includes at least two nominate species (Coregonus artedi and Coregonus sardinella) and an adfluvial C. artedi morph that is distinct from its lacustrine conspecific in terms of life history, morphology, age, growth and mortality. Coregonus sardinella has previously been identified from Great Slave Lake, but we provide the first comprehensive description of this species in the lake and confirm a significant range extension for the species. The lacustrine C. artedi differs little from descriptions throughout its range. In addition to these three ciscoes, linear phenotypic traits, gillraker number and morphology, and growth data support the possible occurrence of two other, less‐distinct morphs, the big‐eye cisco and a shortjaw‐like morph Coregonus zenithicus. Although the big‐eye morph was not identified by body shape and linear phenotypic measures, it was visually identified on the basis of differences in traditional phenotypic proportions, such as orbital length, paired fin lengths, head and gillraker morphology expressed as thousands of standard length and showed different age and growth structure compared with the other lacustrine cisco morphs. Coregonus zenithicus was distinguished visually and by a statistical model of linear phenotypic traits as well as by gillraker number and morphology. Identifying, characterising and managing locally adapted cisco morphs that reflect important ecological and bioenergetic linkages are critical to conserving the ecological integrity of northern ecosystems.

Successional Change in the Lake Superior Fish Community: Population Trends in Ciscoes, Rainbow Smelt, and Lake Trout, 1958--2008

The Lake Superior fi sh community underwent massive changes in the second half of the 20 century. Those changes are largely refl ected in changes in abundance of the adults of principal prey species, the ciscoes (Coregonus spp.), the invasive rainbow smelt (Osmerus mordax), and the principal predator, lake trout (Salvelinus namaycush). To better understand changes in species abundances, a comprehensive series of gillnet and bottom trawl data collected from 1958 to 2008 were examined. In the late 1950s/early 1960s, smelt abundance was at its maximum, wild lake trout was at its minimum, and an abundance of hatchery lake trout was increasing rapidly. The bloater (Coregonus hoyi) was the prevalent cisco in the lake; abundance was more than 300% greater than the next most abundant cisco, shortjaw cisco (C. zenithicus), followed by kiyi (C. kiyi) and lake cisco (C. artedi). By the mid-1960s, abundance of hatchery lake trout was nearing maximum, smelt abundance was beginning to decline, and abundances of all ciscoes declined, but especially that of shortjaw cisco and kiyi. By the late 1970s, recovery of wild lake trout stocks was well underway and abundances of hatchery lake trout and smelt were declining and the ciscoes were reaching their nadir. During 1980–1990, the fi sh community underwent a dramatic shift in organization and structure. The rapid increase in abundance of wild lake trout, concurrent with a rapid decline in hatchery lake trout, signaled the impending recovery. Rainbow smelt abundance dropped precipitously and within four years, lake cisco and bloater populations rebounded on the heels of a series of strong recruitment events. Kiyi populations showed signs of recovery by 1989, and shortjaw by 2000, though well below historic maximum abundances. High abundance of adult smelt prior to 1980 appears to be the only factor linked to recruitment failure in the ciscoes. Life history traits of the cisco species were examined to better understand their different responses to conditions of low and high predator levels, i.e., late 1950s–early 1960s vs. post 1980. Bloaters are most likely to become the predominant cisco in the absence of strong predation and the least abundant under prolonged predation; smelt share this pattern. Conversely, the lake cisco and shortjaw cisco fare better when predator abundance is high. The recovery of lake trout in Lake Superior reestablished a strong top-down infl uence on the fi sh community and its present structure and organization appears to be approaching an equilibrium that refl ects a more natural state. If lake trout recovery is sustained, shortjaw cisco abundance is expected to increase and join lake cisco and kiyi as dominant cisco species, and bloater and smelt will oscillate at lower abundances.

Differential expression of DTSsa4 Tc1-like transposons in closely related populations of Baikal ciscoes

Two representatives of Baikal ciscoes - lake cisco and omul - diverged from a common ancestor as recently as 10-20 thousand years ago. We have found an increasing expression level of DTSsa4 Tc1-like DNA transposons in cisco and omul brains. The mapping of the sequences of these transposons from Salmo salar and Danio rerio genomes has shown that in some cases, these transposons are located in the 5' and 3' regions, as well as in the promoter regions of various genes. Probably, Tc1-like transposons affect the activity of neighboring genes, providing the adaptive divergence of the cisco population.

Candidate Genes and Adaptive Radiation: Insights from Transcriptional Adaptation to the Limnetic Niche among Coregonine Fishes (Coregonus spp., Salmonidae)

In 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 the lake whitefish species complex (Coregonus clupeaformis), previous microarray experiments have led to the identification of candidate genes potentially implicated in the parallel evolution of the limnetic dwarf lake whitefish, which is highly distinct from the benthic normal lake whitefish in life history, morphology, metabolism, and behavior, and yet diverged from it only approximately 15,000 years before present. The aim of the present study was to address transcriptional divergence for six candidate genes among lake whitefish and European whitefish (Coregonus lavaretus) species pairs, as well as lake cisco (Coregonus artedi) and vendace (Coregonus albula). The main goal was to test the hypothesis that parallel phenotypic adaptation toward the use of the limnetic niche in coregonine fishes is accompanied by parallelism in candidate gene transcription as measured by quantitative real-time polymerase chain reaction. Results obtained for three candidate genes, whereby parallelism in expression was observed across all whitefish species pairs, provide strong support for the hypothesis that divergent natural selection plays an important role in the adaptive radiation of whitefish species. However, this parallelism in expression did not extend to cisco and vendace, thereby infirming transcriptional convergence between limnetic whitefish species and their limnetic congeners for these genes. As recently proposed (Lynch 2007a. The evolution of genetic networks by non-adaptive processes. Nat Rev Genet. 8:803-813), these results may suggest that convergent phenotypic evolution can result from nonadaptive shaping of genome architecture in independently evolved coregonine lineages.

Clinal variation at microsatellite loci reveals historical secondary intergradation between glacial races of Coregonus artedi (Teleostei: Coregoninae).

Classical models of the spatial structure of population genetics rely on the assumption of migration-drift equilibrium, which is seldom met in natural populations having only recently colonized their current range (e.g., postglacial). Population structure then depicts historical events, and counfounding effects due to recent secondary contact between recently differentiated lineages can further counfound analyses of association between geographic and genetic distances. Mitochondrial polymorphisms have revealed the existence of two closely related lineages of the lake cisco, Coregonus artedi, whose significantly different but overlaping geographical distributions provided a weak signal of past range fragmentation blurred by putative subsequent extensive secondary contacts. In this study, we analyzed geographical patterns of genetic variation at seven microsatellite loci among 22 populations of lake cisco located along the axis of an area covered by proglacial lakes 12,000-8,000 years ago in North America. The results clearly confirmed the existence of two genetically distinct races characterized by different sets of microsatellite alleles whose frequencies varied clinally across some 3000 km. Equilibrium and nonequilibrium analyses of isolation by distance revealed historical signal of gene flow resulting from the nearly complete admixture of these races following neutral secondary contacts in their historical habitat and indicated that the colonization process occurred by a stepwise expansion of an eastern (Atlantic) race into a previously established Mississippian race. This historical signal of equilibrium contrasted with the current migration-drift disequilibrium within major extant watersheds and was apparently maintained by high effective population sizes and low migration regimes.

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.

CLINAL VARIATION AT MICROSATELLITE LOCI REVEALS HISTORICAL SECONDARY INTERGRADATION BETWEEN GLACIAL RACES OF COREGONUS ARTEDI (TELEOSTEI: COREGONINAE)

Classical models of the spatial structure of population genetics rely on the assumption of migration-drift equilibrium, which is seldom met in natural populations having only recently colonized their current range (e.g., postglacial). Population structure then depicts historical events, and counfounding effects due to recent secondary contact between recently differentiated lineages can further counfound analyses of association between geographic and genetic distances. Mitochondrial polymorphisms have revealed the existence of two closely related lineages of the lake cisco, Coregonus artedi, whose significantly different but overlaping geographical distributions provided a weak signal of past range fragmentation blurred by putative subsequent extensive secondary contacts. In this study, we analyzed geographical patterns of genetic variation at seven microsatellite loci among 22 populations of lake cisco located along the axis of an area covered by proglacial lakes 12,000–8000 years ago in North America. The results clearly confirmed the existence of two genetically distinct races characterized by different sets of microsatellite alleles whose frequencies varied clinally across some 3000 km. Equilibrium and nonequilibrium analyses of isolation by distance revealed historical signal of gene flow resulting from the nearly complete admixture of these races following neutral secondary contacts in their historical habitat and indicated that the colonization process occurred by a stepwise expansion of an eastern (Atlantic) race into a previously established Mississippian race. This historical signal of equilibrium contrasted with the current migration-drift disequilibrium within major extant watersheds and was apparently maintained by high effective population sizes and low migration regimes.Corresponding Editor: T. Smith

Multiple Trophic Level Structuring in Salvelinus–Coregonus Assemblages in Boreal Forest Lakes

Hypolimnetic trophic interactions were examined among lake trout (Salvelinus namaycush), cisco (Coregonus artedi), zooplankton, and macroinvertebrates in six northwestern Ontario lakes varying sevenfold in conductivity. Faster growth and larger body size of lake trout in high-conductivity lakes were associated with larger meal sizes, earlier entry of cisco into their diet, and consumption of relatively few, large-bodied prey. Cisco preyed on zooplankton and macroinvertebrates. Cisco were abundant, large bodied, and fast growing in high-conductivity lakes. In one low-conductivity lake (Greenwich Lake), lake trout were 12 times as abundant as cisco whereas in all other lakes, lake trout to cisco ratios were ~ 1:1 or less regardless of conductivity. In Greenwich Lake, lake trout consumed large quantities of age-0 cisco and Mysis relicta, which was associated with "top-down" trophic structuring. Food web analyses strongly suggest that in Greenwich Lake the high abundance of piscivores resulted directly in a low abundance of planktivores and indirectly in a large-bodied zooplankton community and planktivores with rapid growth rates.

Mitochondrial DNA Variation in Four Minnesota Populations of Lake Whitefish: Utility as Species and Population Markers

Abstract Restriction fragment length polymorphism (RFLP) analysis was performed on the mitochondria) DNA (mtDNA) of 33 lake whitefish Coregonus clupeaformis from four spawning populations in north-central Minnesota. Diagnostic differences between the RFLP patterns of mtDNA from lake whitefish and cisco C. artedi were useful in the identification of fresh fillets and fresh or frozen carcasses. However, mtDNA analysis was of little use for identification of lake whitefish populations because of the low diversity documented among study populations.

Genetic comparisons of New and Old World coregonid fishes*

The genetic relationships among New and Old World coregonid fishes were studied by electrophoresis. The genetic composition of 60 populations, representing perhaps nine commonly recognized species of Coregonus and Stenodus from Europe and North America was determined for 37 genetic loci. Six distinct genetic groups were evident. The first contained only populations of the inconnu, Stenodus leucichthys (Güldenstadt). The Nei genetic distance between Stenodus and Coregonus was 0.305, a relatively small value as compared to other salmonid inter‐generic comparisons. The second genetic grouping contained the Arctic cisco, C. autumnalis (Pallas), the N. American lake cisco, C. artedii Lesueur, and the Irish pollan, C. autumnalis pollan Thompson. These three taxa appear to be conspecific on the basis of genetic distances. The third genetic grouping contained the European whitefish, C. lavaretus (L.), and the N. American lake whitefish, C. clupeaformis (Mitchill). European and lake whitefish may be conspecific. Lake whitefish from northwestern N. America were more closely related to European whitefish (genetic distance 0.038) than to lake whitefish from central N. America (genetic distance 0.098). The fourth group contained the broad whitefish, C. nasus (Pallas), which is perhaps more closely related to the European/lake whitefish groups than other coregonids. The fifth genetic grouping contained only the Asian endemic, C. peled (Gmelin), and the sixth contained the least cisco, C. sardinella Valenciennes, and vendace, C. albula (L.), which also appear to be conspecific. The widespread genetic groupings obtained for ciscoes indicate that they do not constitute a single closely related group within the genus Coregonus.

Lake Trout (Salvelinus namaycush) Growth Potential Predicted from Cisco (Coregonus artedii) Population Structure and Conductivity

Trophic dynamics between lake trout (Salvelinus namaycush) and cisco (Coregonus artedii) were examined in a set of 10 northwestern Ontario lakes that ranged in conductivity by nearly an order of magnitude. Diets revealed the possibility of interspecific competition for invertebrate prey prior to initiation of piscivory by lake trout at approximately 200 mm fork length (FL). Von Bertalanffy growth equations (based on otolith ages) were used to estimate time to attain this length (1.4–3.0 yr). Contrary to anticipated results, lake trout reached this size more rapidly in lakes with high cisco abundance. Presumably, this relationship was caused by improved food conditions for young trout in some of the more productive waterbodies. However, after removing the collinearity between lake conductivity and nonpiscivorous growth by regression analysis, there was no relationship between the time for lake trout to reach 200 mm and cisco abundance. Piscivorous growth varied (L∞ = 508 to 740 mm) extensively among the 10 populations and was likely a function of growth efficiency from different food sources. Using stepwise multiple regression, 81% of the L∞ variation among lakes was explained by cisco abundance and mean size. Linear regression between L∞ and mean size of trout (but not trout abundance) on conductivity generated significant positive relationships.