Geology topics

Daniel L. Yule

Publications and source records attributed to Daniel L. Yule.

At least 19 recordsLinked to original sources

USGS2026

Status and trends in the Lake Superior fish community, 2025

The U.S. Geological Survey has conducted annual trawl surveys across Lake Superior since 1978 that describe trends in fish species occurrence and relative abundance to support fisheries science and management. In 2025, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 72 nearshore stations in May and June and 36 offshore locations in July. Nearshore bottom trawls collected 63,005 fish represented by 30 species or morphotypes. The number of species collected at each location ranged from 1 to 13, with a median of 7.0 species. Estimated fish biomass density at individual stations ranged from <0.1 to 125.7 kg per ha with a lakewide mean of 7.6 kg per ha. Offshore bottom trawls collected 30,342 fish represented by 13 species or morphotypes. Estimated fish biomass density at individual stations ranged from 1.6 to 34.7 kg per ha with a lakewide mean of 9.1 kg per ha, which was the highest since the offshore survey began in 2011. Lakewide average numerical densities (fish per ha) of age-1 fish were 0.67 per ha for Bloater, 0.01 per ha for Cisco, 1.15 per ha for Kiyi, 0.95 per ha for Lake Whitefish, and 336.33 per ha for Rainbow Smelt. Surface trawling collected 1,562 larval Coregonus individuals which was the fewest Coregonus larvae collected in a whole lake survey since the larval fish survey began in 2014. Nearshore mean larval Coregonus numerical densities were 156 fish per ha in May and June 2025 and 21 fish per ha in July 2025. May mean surface water temperatures (6.1°C) were near the warmest for the period-of-record, while June (6.8°C) and July (10.0°C) were near average or below.

Lake Superior↗

USGS2026

Evaluating the use of uncrewed surface vessels to enhance Lake Erie acoustic prey-fish surveys

Incorporation of autonomous uncrewed surface vessels (USVs) into large-scale acoustic surveys may enhance spatiotemporal extent and quality of fish density estimates. Lake Erie is currently surveyed by three motorized research vessels (RVs), which annually collect acoustic data and estimate prey-fish abundances. To evaluate the feasibility of incorporating a USV into the existing survey to increase spatial coverage, we compared paired acoustic measurements of target strength (TS), area backscattering coefficient (ABC), and areal fish density from four motorized RVs (three currently and one historically used) and one USV (Saildrone Explorer). Acoustic data were collected along 2-km transects during five sampling events distributed throughout Lake Erie. Despite observing event-specific differences that were likely impacted by sampling conditions, survey-level estimates of areal density between RVs and the USV were comparable (i.e., RV 6% lower than USV). RVs and the USV appear to react differently (speed and stability) to variable environmental conditions (wind and waves), which impacted data quality and must be mitigated in future applications. However, similar to previous comparisons in marine and freshwater environments, incorporating a USV into an acoustic survey could enhance annual fish density estimates by increasing effort, spatial coverage, logistical flexibility, and biological data collection opportunities from RVs.

Michigan, New York, Ohio, Pennsylvania↗

USGS2026

A genomic tool to tackle cryptic diversity demonstrates the potential for off-target use of GT-seq panels

A comprehensive understanding of life history is vital to successful species conservation and management. When different life history stages are accompanied by considerable morphological or cryptic variation, such as the egg and larval phases exhibited by most fishes, genomic tools are essential for identifying species so that early-life ecology questions can be studied. Genotyping-in-thousands by sequencing (GT-seq) has recently emerged as a targeted and efficient approach for species identification. We leveraged existing genomic and transcriptomic data to develop a GT-seq panel capable of differentiating the members of the Coregonus artedi complex, a radiation of salmonids in the Laurentian Great Lakes whose members are indistinguishable with mitochondrial DNA barcoding loci and are the focus of bi-national conservation initiatives. Our panel of 494 loci was able to assign fishes in the C. artedi complex to species and lake. We examined cross-amplification in other coregonines with overlapping distributions and found that congeneric Lake Whitefish ( C. clupeaformis ) cross-amplified at 94% of loci and confamilial Round and Pygmy Whitefish ( Prosopium spp.) cross-amplified at 42% and 38% of loci, respectively. We adapted bioinformatic probes to account for Prosopium -specific variants including 22 new SNPs and developed a whitelist of 428 SNPs capable of distinguishing these whitefishes. Finally, we demonstrated performance by identifying 3,066 coregonine larvae and juveniles collected in spring 2019-2021 from Lake Superior. These results hold promise for future insights into the species-specific ecology of early life coregonines and demonstrate the flexibility of GT-seq panels, which may cross-amplify hundreds of informative genome-wide loci in related taxa.

BioRxiv↗

USGS2025

Food web structure of the Lake Superior fish community in 2021–2022

The trophic linkages and ecological requirements of the Lake Superior fish community have not been assessed on a whole lake scale in over a decade. Here, we investigated the trophic dynamics across multiple species and habitat zones of Lake Superior. From April to October of 2021 and 2022, a total of five piscivore, four planktivore, and six benthivore species were collected by region and length class during bottom-trawl, standard gillnet, commercial cisco gillnet, and recreational angler surveys. To assess trophic linkages, stomach contents were measured to estimate biomass consumed and a multivariate analysis was used to assess diet composition by species, length, and region. We found a high degree of interconnectedness in the Lake Superior food web, with Mysis as a critical diet item for most fishes. Native piscivore diets varied by region and with ontogeny. Lake charr were important habitat couplers in the lake, exhibiting a diverse diet and opportunistic foraging strategy. Conversely, Pacific salmon were more restricted in their foraging. Planktivores and benthivores primarily consumed Mysis , with less reliance on Diporeia compared to previous studies. Lake whitefish and cisco were the exception to this pattern, with broader bathymetric depth distributions represented in their diets. We found the food web to be supported by a predominantly native species assemblage, with redundancies at all trophic levels.

Lake Superior↗

USGS2025

Partial diel vertical migration and niche partitioning in Mysis revealed by stable isotopes

Diel vertical migration (DVM) is critical for moving energy and nutrients between surface and deep waters. Mysis sp. (Crustacea: Mysidae) facilitates this process by serving as predator and prey in both benthic and pelagic habitats. Mysis can also exhibit partial DVM (pDVM), where some individuals do not migrate into the pelagia at night or to the benthos during the day. However, whether Mysis pDVM is a fixed (i.e., same individuals migrate) or random (i.e., random individuals migrate) behavior remains unclear. To evaluate that, we tested whether Mysis exhibit pDVM and niche partitioning in Lake Superior by collecting Mysis from benthic and pelagic habitats day and night across depths ranging from 50 to 250 m and estimating their isotopic niche size (δ 13 C, δ 15 N; Corrected Standard Ellipse Area − SEA c ) and overlap among three life stages. At the population level, Mysis exhibited fixed pDVM structured by life stage. Benthic Mysis (12.9 ± 3.2 mm, mean ± SD) were larger than pelagic Mysis during night (9.6 ± 3.6 mm) and day (8.5 ± 3.6 mm). Adult Mysis (> 15 mm) had larger SEA c (1.8 ± 0.4 ‰ 2 ) compared to juveniles (< 10 mm; 0.3 ± 0.1 ‰ 2 ) and sub-adults (10–15 mm; 0.7 ± 0.2 ‰ 2 ), and their isotopic niche did not overlap with smaller life stages. Adults exhibited random pDVM (i.e., high isotopic niche overlap), whereas juveniles exhibited fixed pDVM (i.e., low overlap). Our observations indicate complex behaviors across and within Mysis life stages, likely due to varying pressures including size-selective predation, light and temperature thresholds, and nutritional requirements. Consequently, the benthic behavior of adult Mysis needs to be considered in monitoring programs.

Minnesota, Wisconsin↗

USGS2025

Lake depth and light conditions alter Mysis vertical distributions

Light regulates the vertical migration of many aquatic organisms. Mysis species couple pelagic and benthic habitats in lakes by diel vertical migrations (DVM), transporting energy and nutrients through the water column and food web. Although Mysis are generally assumed to remain on the bottom during the day, some have been observed in the pelagic zone during the day, indicating incomplete benthic-pelagic coupling in some systems. The degree to which light attenuation and lake depth interact to affect occurrence of mysids within the water column during the day is understudied. We used standardized Mysis net sampling in summers 2020 and 2021 across nine north-temperate lakes to test the hypotheses that 1) Mysis remain pelagic during the day at depths with sufficiently low light levels, and 2) pelagic-caught individuals during the day are, on average, smaller than those caught at night. To test these hypotheses, we assessed light, dissolved oxygen (DO), Mysis densities, and size distribution between night and day across bathymetric depths. In deep lakes and darkly colored shallow lakes, Mysis suspended in the water column during the day where light levels decreased to their light avoidance threshold (∼10 -5 to 10 -6 lx ). Mysis suspended in the water column during the day were smaller than those collected at night. Further, Mysis were not captured when DO reached levels < 3 mg/L, regardless of light conditions. Our results suggest that benthic-pelagic coupling by Mysis is mediated through light conditions, lake morphometry, and DO conditions, and may include some degree of size-dependent behavior.

Journal of Great Lakes Research↗

USGS2025

Archive of morphological data for the Coregonus artedi species complex of the Great Lakes, Lake Nipigon and Great Slave Lake

This publication is a user guide for an archive of morphological data recorded by various authors from North American ciscoes of the Coregonus artedi species complex (subfamily Coregoninae ). The archive is accessible from the Great Lakes Fishery Commission’s (GLFC) server, is open access, and contains data for the Laurentian Great Lakes; Lake Nipigon, Ontario; and Great Slave Lake, Northwest Territories. The archive comprises morphometrics and meristics (together metrics) for 6,700 individual Cisco of which 1,400 are accompanied by images. In addition, the archive contains metrics presented as arrays by W. N. Koelz, Coregonid fishes of the Great Lakes, Bulletin of the U.S. Bureau of Fisheries 43(2):297-643, which were based on 10,000 individuals. Spreadsheets in the Metrics folder of the archive are divided broadly into Contemporary and Historical subfolders and the Contemporary subfolder is further divided into Cisco Monograph and Extra Monograph subfolders to encourage statistical assessment of findings in GLFC Miscellaneous Publication 2023. The Images folder is organized into subfolders by lake. Tables in this user guide allow for quick determination of the availability of data by lake, subspecies, author, and year.

Laurentian↗

USGS2025

Status and trends in the Lake Superior fish community, 2024

The U.S. Geological Survey has conducted annual fishery surveys across Lake Superior since 1978 that describe trends in fish species occurrence and relative abundance to inform fisheries management and ecosystem health. In 2024, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 72 nearshore locations in June and 36 offshore locations in July. Nearshore bottom trawls collected 22,190 fish represented by 27 species or morphotypes. The number of species collected at each location ranged from 1 to 12, with a median of 5.5 species. Estimated fish biomass at individual locations ranged from <0.1 to 62.9 kg per ha with a lakewide mean of 3.7 kg per ha. Offshore bottom trawls collected 33,634 fish represented by 12 species or morphotypes. Estimated fish biomass at individual locations ranged from 0.6 to 25.8 kg per ha with a lakewide mean of 8.3 kg per ha, which was the second highest for the period-of-record. Lakewide average densities (fish per ha) of age-1 fish were 1 per ha for Bloater, 5 per ha for Cisco, 1 per ha for Lake Whitefish, 60 per ha for Rainbow Smelt, and 19 per ha for Kiyi. Surface trawling collected 5,177 larval Coregonus individuals which was the third fewest Coregonus larvae collected in a whole lake survey since the larval fish survey began in 2014. Nearshore mean larval Coregonus densities were 176 fish per ha in June 2024 and offshore densities were 7 fish per ha in July 2024. June and July surface water temperatures were near the warmest for the period-of-record.

Lake Superior↗

USGS2024

Origin of the Laurentian Great Lakes fish fauna through upward adaptive radiation cascade prior to the Last Glacial Maximum

The evolutionary histories of adaptive radiations can be marked by dramatic demographic fluctuations. However, the demographic histories of ecologically-linked co-diversifying lineages remain understudied. The Laurentian Great Lakes provide a unique system of two such radiations that are dispersed across depth gradients with a predator-prey relationship. We show that the North American Coregonus species complex (“ciscoes”) radiated rapidly prior to the Last Glacial Maximum (80–90 kya), a globally warm period, followed by rapid expansion in population size. Similar patterns of demographic expansion were observed in the predator species, Lake Charr ( Salvelinus namaycush ), following a brief time lag, which we hypothesize was driven by predator-prey dynamics. Diversification of prey into deep water created ecological opportunities for the predators, facilitating their demographic expansion, which is consistent with an upward adaptive radiation cascade. This study provides a new timeline and environmental context for the origin of the Laurentian Great Lakes fish fauna, and firmly establishes this system as drivers of ecological diversification and rapid speciation through cyclical glaciation.

Communications Biology↗

USGS2024

Factors influencing larval coregonine spatial distribution in Lake Geneva (Europe) and Lake Superior (North America) during a single season near known spawning sites

Survival rate of the larval stage is an important driver of fish recruitment. To understand mechanisms regulating larval survival it is important to understand the relative importance of abiotic and biotic factors that shape larval spatial distributions. We studied larval Coregonus distributions in surface waters (surface to 1 m) by repeatedly sampling study sites in two lakes that varied greatly in trophic state and regional climate. We evaluated the importance of bathymetric depth, Julian Day, edible zooplankton densities (EZ, ind. L −1 ) and wind vectors on larval spatial distributions using generalized additive modeling. In both systems, larval counts declined in a negative exponential fashion with bathymetric depth, indicating shallow depths are critical nursery habitat. The north-south wind vectors and Julian Day (which was positively correlated with surface temperature) influenced larval distributions in Lake Geneva with larval counts related to both variables linearly, whereas the east-west wind vector and EZ were unimportant. Highest larval counts were during an offshore south wind and declined slightly with Julian Day. In Lake Superior, bathymetric depth and the east-west wind vector influenced larval distributions and were unrelated to EZ, Julian Day, and the north-south wind vector. Larval counts were highest when onshore southwest winds preceded sampling. Differences in how wind affected larval distribution (offshore vs. onshore) might be related to larval size with Lake Superior larvae considerably smaller (average length 12.9 mm vs. 15.9 mm); thus, more apt to be subjected to advection. Within coastal waters, Julian Day and wind vectors influence distributions, but their importance seemingly varies lake-to-lake.

Minnesota↗

USGS2024

Trophic ecology of juvenile lean and siscowet lake charr (Salvelinus namaycush) in Lake Superior: Assessing for potential competition

We investigated the spatial overlap, diet, isotopic niche, and growth of juvenile lean and siscowet lake charr ( Salvelinus namaycush ) in Lake Superior to address concerns of potential competition with implications to the study of resource polymorphism. Catch data revealed the greatest levels of sympatry in waters from 40 to 60 m. Juvenile lean and siscowet diet changed ontogenetically with Mysis being the dominant prey item for the smallest lake charr but differentiating with onset of piscivory. As ecotypes increased in size, lean diets became dominated by pelagic prey, whereas siscowets had equal proportions of benthic and pelagic prey. Isotopic niche overlap declined between ecotypes coincident with siscowet lake charr shifting to deeper habitats around 400 mm. Lean and siscowet exhibited different growth trajectories. However, length at age-4 declined in parallel for both ecotypes with no trend in condition suggesting that lake charr growth is sensitive to prey biomass and unlikely related to competition. Our findings indicate minimal evidence of competition and support the concept that multiple sympatric ecotypes of lake charr in Lake Superior are maintained by resource polymorphism.

Lake Superior↗

USGS2024

Paired comparisons with quiet surface drones show evidence of fish behavioral response to motorized vessels during acoustic surveys in Lake Superior

Acoustic surveys are important for fish stock assessments, but fish responses to survey vessels can bias acoustic estimates. We leveraged quiet uncrewed surface vessels (USVs) to characterize potential bias in acoustic surveys. Five conventional motorized ships overtook USVs from astern over 2 km transects at night in Lake Superior in 2022. We examined the difference in acoustic backscatter, average target depth, and average target strength (TS) between USV and motorized vessels. Although sound level measurements from the motorized vessels sometimes exceeded recommendations for scientific vessels, we did not detect differences in acoustic measures among survey vessels. However, the USVs recorded 2 dB higher acoustic backscatter and TS than motorized vessels, leading to ~15% higher fish densities with drones when using in situ TS and echo integration. Differences in fish density would increase to 30-60% if a standard TS value was applied. Target depth did not differ between USVs and motorized ships. These results are consistent with a change in orientation but not depth of insonified fish and limited horizontal avoidance of motorized survey vessels.

Lake Superior↗

USGS2023

Status and trends in the Lake Superior fish community, 2022

In 2022, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 71 nearshore locations in May-June and 35 offshore locations in July, and at 51 Coordinated Science and Monitoring Initiative (CSMI) locations in July-October with bottom trawls, surface trawls, mid-water trawls and acoustics that were previously sampled in 2011 and 2016. Nearshore bottom trawls collected 11,603 fish from 25 species or morphotypes. Nearshore mean biomass was 1.6 kg per ha which was one of the lowest biomass estimates over survey’s 45-year history. Offshore bottom trawls collected 13,876 fish from 11 species or morphotypes. Offshore mean biomass was 5.1 kg per ha, which was less than the annual average since 2011 of 6.5 kg per ha. Recruitment, as measured by age-1 densities, was near zero for Bloater ( Coregonus hoyi ), Cisco ( C. artedi ), and Kiyi ( C. kiyi ), 2 age-1 fish per ha for Lake Whitefish ( C. clupeaformis ) and 77 age-1 fish for Rainbow Smelt ( Osmerus mordax ). All were less than the long-term averages. Sampling at the CSMI locations collected 26 species and morphotypes. The most abundant species’ lakewide were Deepwater Sculpin (all years), young-of-year ciscoe (Bloater, Cisco, and Kiyi, 2022), and Rainbow Smelt (2011 and 2016). Cisco had the highest estimated lakewide biomass in 2011 and 2022 and siscowet Lake Trout had the highest estimated lakewide biomass in 2016. Native species were more abundant than invasive species by numbers (80, 65, and 92%) and biomass (94, 93, 96%) in 2011, 2016, and 2022, respectively. Total lakewide benthic fish biomass declined from 47 thousand metric tons in 2011 to 29 thousand metric tons in 2016 and increased to 33 thousand metric tons in 2022. Total lakewide pelagic fish biomass declined from 61 thousand metric tons in 2011 to 25 thousand metric tons in 2016 and increased to 54 thousand metric tons in 2022. The most unexpected result from our sampling in 2022 was the 2 billion age-0 ciscoe estimate from the mid-water trawl and acoustic sampling in August-October. These fish were broadly distributed across the lake, being collected at 53 of the 54 locations, and their population estimates were highest in the depths >100 m. The factors underlying the survival of these ciscoes into late summer in 2022 as compared to previous years have not been identified, but our annual population surveys of larval ciscoes suggests that lake conditions in June and July may have differed from previous years and enhanced survival. In 2022, ciscoe larval densities in May were lower than average (likely due to a cold winter and spring that delayed hatching), June densities were similar to previous years, and July density estimates were more than double that of any previous year’s estimate.

Lake Superior↗

USGS2023

Status and trends in the Lake Superior fish community, 2023

The U.S. Geological Survey annually conducts fishery surveys across Lake Superior that describe trends in fish species occurrence and relative abundance to inform fisheries management and large lake ecology. In 2023, the Lake Superior fish community was sampled with daytime bottom and surface trawls at 51 nearshore locations in June and 31 offshore locations in July. Nearshore bottom trawls collected 157,804 fish from 25 species or morphotypes. Nearshore mean biomass was 18.3 kg per ha which was the second highest biomass estimate over the survey’s 46-year history. Offshore bottom trawls collected 15,458 fish from 10 species or morphotypes. Offshore mean biomass was 5.2 kg per ha, which was less than the annual average of 6.3 kg per ha. Recruitment, as measured by age-1 densities, was the highest recorded for Bloater, Cisco, and Rainbow Smelt in the nearshore and for Kiyi in the offshore survey’s period-of-records. Lakewide average densities (fish per ha) of age-1 fish were 140 for Bloater, 1,019 for Cisco, 616 for Rainbow Smelt, and 54 for Kiyi, which were the highest estimates for the survey’s period-of-record. Period-of-record averages for these species were 10, 67, and 9 age-1 fish per ha, respectively. Age-1 Lake Whitefish averaged 9 fish per ha which was similar to the long-term average of 8 age-1 Lake Whitefish per ha. If the future can be predicted by past large Bloater, Cisco, and Kiyi (collectively, ciscoe) year-class events, the unprecedented survival of the 2022 ciscoe yearclass will influence the Lake Superior ecosystem for the next 10 to 20-years.

Lake Superior↗

USGS2022

Lake Superior Kiyi reproductive biology

Objective The Lake Superior Kiyi Coregonus kiyi is an understudied species being considered for reintroduction into Laurentian Great Lakes where it no longer occurs. Herein, we provide descriptions of Kiyi reproductive biology with the intention of guiding potential gamete collections for propagation. Methods Data were collected on Kiyi spawning timing, spawning locations, spawning season catch rates, length at sexual maturity, sex ratios, fecundity, egg size, and larval occurrences in Lake Superior from 1996–2021. These data were compared to observations made a century prior in Lakes Michigan, Ontario, and Superior. Result Contemporary Kiyi spawning occurred between late December and late January when surface water temperatures cooled to <4°C. Spawning Kiyi were caught almost exclusively in 38.1-mm stretch mesh, as compared to larger meshes (50.8–76.2 mm). Capture depths for developing, ripe, running, and spent female Kiyi were similar and ranged from 82 to 221 m. Fifty percent of female and male Kiyi were classified as sexually mature at ~150 mm total length. Fecundity estimates ranged from 1,578 to 6,720 eggs/female. Mean diameter of unfertilized eggs was 1.7 mm. Recently hatched larval Kiyi were collected at the surface during May–July at 62 of the 113 locations sampled throughout the lake in 2019. Conclusion Our work suggests that Kiyi gamete collection efforts from mid-December through January using 38.1-mm gill-net panels set at bathymetric depths of at least 100 m would maximize the collection of spawning Kiyi and reduce the bycatch of other Coregonus species. Future research questions include the following: (1) “Do Kiyi form spawning aggregations at specific spawning areas, or do they spawn indiscriminately across the lake?”; (2) “Do Kiyi spawn near the bottom or up in the water column?”; (3) “What is the relationship between fall lake overturn and Kiyi spawn timing?”; and (4) “Could summer larval and age-0 Kiyi collections provide an opportunity for establishing a captive broodstock?”

Lake Superior↗

USGS2022

Concordant patterns of morphological, stable isotope, and genetic variation in a recent ecological radiation (Salmonidae: Coregonus spp.)

Groups of sympatric taxa with low interspecific genetic differentiation, but considerable ecological differences, offer great opportunities to study the dynamics of divergence and speciation. This is the case of ciscoes ( Coregonus spp.) in the Laurentian Great Lakes, which are characterized by a complex evolutionary history and are commonly described as having undergone an adaptive radiation. In this study, morphometrics, stable isotopes and transcriptome sequencing were used to study the relationships within the Coregonus artedi complex in western Lake Superior. We observed general concordance for morphological, ecological and genomic variation, but the last was more taxonomically informative as it showed less overlap among species in multivariate space. Low levels of genetic differentiation were observed between individuals morphologically identified as Coregonus hoyi and C. zenithicus , which could be evidence of incomplete lineage sorting or recent hybridization between the two groups. Transcriptome-based single nucleotide polymorphisms exhibited significant divergence for genes associated with vision, development, metabolism and immunity among species that occupy different habitats. This study highlights the importance of using an integrative approach when studying groups of taxa with a complex evolutionary history, as individual-level analyses of multiple independent data sets can provide a clearer picture of the patterns and processes associated with the origins of biodiversity.

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2022

Variation of cisco egg size among Laurentian Great Lakes populations

Many fish species display inter-population and inter-individual egg size variation. Intra-specific differences in egg size seemingly reflect both energetic experiences of individual spawning fish and long-term population responses to differing ecosystems. Optimal egg size theory implies that selection influences a population’s mean egg size in response to its early-life environment, given the well-established trade-off between egg size and fecundity. Currently, there is strong interest in rehabilitation of Laurentian Great Lakes cisco, Coregonus artedi, which is characterized by inter-population variation of morphological and behavioral traits. However, the extent of cisco egg size variation is under-described. In fall 2018 and 2019, we collected egg samples by stripping ripe females at seven total locations in four Great Lakes. We measured unfertilized egg diameters using imaging software and compared mean egg diameters among locations with and without including maternal total length as a covariate. Lake Michigan females produced the largest eggs overall but were excluded from analyses using the total length covariate because of their significantly larger body sizes. Maternal length had a positive effect on egg size, and when accounting for this effect, females in Lake Huron produced the largest eggs followed by Lake Ontario and Lake Superior. We also found that egg size varied among locations within Lake Superior. These findings aligned with observations of morphological and behavioral differences among populations and suggest that cisco phenotypic variation at a fine spatial scale extends to reproductive biology. Consideration of cisco reproductive traits, such as egg size, may inform restoration strategies, including supplemental stocking.

Journal of Great Lakes Research↗

USGS2022

Genetic population structure of cisco, Coregonus artedi, in the Laurentian Great Lakes

Management of a widely distributed species can be a challenge when management priorities, resource status, and assessment methods vary across jurisdictions. For example, restoration and preservation of coregonine species diversity is a goal of management agencies across the Laurentian Great Lakes. However, management goals and the amount of information available varies across management units, making the focus for management efforts challenging to determine. Genetic data provide a spatially consistent means to assess diversity. Therefore, we examined the genetic stock structure of cisco ( Coregonus artedi ) in the Great Lakes where the species is still extant. Using genotype data from 17 microsatellite DNA loci, we observed low levels of population structure among collections with most contributions to overall diversity occurring among lakes. Cisco from lakes Superior, Michigan, Ontario, and the St. Marys River could be considered single genetic populations while distinct genetic populations were observed among samples from northern Lake Huron. Significant within-lake diversity in Lake Huron is supported by populations found in embayments in northern Lake Huron. The Grand Traverse Bay population in Lake Michigan represents a distinct population with reduced levels of genetic variation when compared to other lakes. The different levels of within lake population structure we observed will be important to consider as future lake-specific management plans are developed.

Great Lakes↗