Expansion of Dreissena into offshore waters of Lake Michigan and potential impacts on fish populations

Lake Michigan was invaded by zebra mussels (Dreissena polymorpha) in the late 1980s and then followed by quagga mussels (D. bugensis) around 1997. Through 2000, both species (herein Dreissena) were largely restricted to depths less than 50 m. Herein, we provide results of an annual lake-wide bottom trawl survey in Lake Michigan that reveal the relative biomass and depth distribution of Dreissena between 1999 and 2007 (although biomass estimates from a bottom trawl are biased low). Lake-wide mean biomass density (g/m²) and mean depth of collection revealed no trend between 1999 and 2003 (mean = 0.7 g/m² and 37 m, respectively). Between 2004 and 2007, however, mean lake-wide biomass density increased from 0.8 g/m² to 7.0 g/m², because of increased density at depths between 30 and 110 m, and mean depth of collection increased from 42 to 77 m. This pattern was confirmed by a generalized additive model. Coincident with the Dreissena expansion that occurred beginning in 2004, fish biomass density (generally planktivores) declined 71% between 2003 and 2007. Current understanding of fish population dynamics, however, indicates that Dreissena expansion is not the primary explanation for the decline of fish, and we provide a species-specific account for more likely underlying factors. Nonetheless, future sampling and research may reveal a better understanding of the potential negative interactions between Dreissena and fish in Lake Michigan and elsewhere.     

Additional deepwater sampling of prey fish in Lake Michigan annual bottom trawl survey reveals new insights for depth distribution dynamics

Beginning in 2013, sites at the 128-m bottom depth contour were added to the sampling design of the annual Lake Michigan bottom trawl survey for prey fish, which has been conducted by the U.S. Geological Survey Great Lakes Science Center (GLSC) each fall since 1973, to better assess fish depth distributions in a changing ecosystem. The standard sampling design included bottom depths from 9 to 110 m, although the GLSC also sporadically conducted bottom trawl tows at the 128-m bottom depth contour during 1973–1988. Enactment of this new sampling design in 2013 revealed that mean biomass density of deepwater sculpins (Myoxocephalus thompsonii) at the 128-m depth exceeded the sum of mean biomass densities at shallower depths, indicating that the bulk of the deepwater sculpin population is residing in waters deeper than 110 m. Thus, our findings supported the hypothesis that the depth distribution of the deepwater sculpin population had shifted to deeper waters beginning in 2007, thereby explaining, at least in part, the marked decline in deepwater sculpin abundance since 2006 based on the standard sampling design. In contrast, our results did not support the hypothesis that the slimy sculpin (Cottus cognatus) population had shifted to deeper waters sometime after 2000. A portion of the burbot (Lota lota) population may have also shifted in depth distribution to waters deeper than 110 m after 2007, based on our results. Our findings have served as an impetus to further expand the range of depths sampled in our bottom trawl survey.     

Shell-free Biomass and Population Dynamics of Dreissenids in Offshore Lake Michigan, 2001–2003

The USGS-Great Lakes Science Center has collected dreissenid mussels annually from Lake Michigan since zebra mussels (Dreissena polymorpha) became a significant portion of the bottomtrawl catch in 1999. For this study, we investigated dreissenid distribution, body mass, and recruitment at different depths in Lake Michigan during 2001–2003. The highest densities of dreissenid biomass were observed from depths of 27 to 46 m. The biomass of quagga mussels (Dreissena bugensis) increased exponentially during 2001–2003, while that of zebra mussels did not change significantly. Body mass (standardized for a given shell length) of both species was lowest from depths of 27 to 37m, highest from 55 to 64 m, and declined linearly at deeper depths during 2001–2003. Recruitment in 2003, as characterized by the proportion of mussels < 11 mm in the catch, varied with depth and lake region. For quagga mussels, recruitment declined linearly with depth, and was highest in northern Lake Michigan. For zebra mussels, recruitment generally declined non-linearly with depth, although the pattern was different for north, mid, and southern Lake Michigan. Our analyses suggest that quagga mussels could overtake zebra mussels and become the most abundant mollusk in terms of biomass in Lake Michigan.     

A short-term look at potential changes in Lake Michigan slimy sculpin diets

Diporeia hoyi and Mysis relicta are the most important prey items of slimy sculpins (Cottus cognatus) in the Great Lakes. Slimy sculpins were collected from dreissenid-infested bottoms off seven Lake Michigan ports at depths of 27–73 m in fall 2003 to study their lake-wide diets. Relatively large dreissenid biomass occurred at depths of 37- and 46-m. Quagga mussels (Dreissena bugnesis) composed at least 50% of dreissenid biomass at Manistique, Saugatuck, and Sturgeon Bay. Mysis accounted for 82% of the sculpin diet by dry weight at eastern Lake Michigan while Diporeia composed 54–69% of the diet at western Lake Michigan and dominated the diets of slimy sculpins at all sites deeper than 46 m. In northern Lake Michigan, this diet study in new sites showed that slimy sculpin consumed more prey with low energy contents, especially chironomids, than Mysis and Diporeia in shallow sites (depth <55 m). We recommend diet studies on sedentary benthic fishes to be conducted along perimeters of the Great Lakes to observe changes in their diets that may be impacted by changing benthic macroinvertebrate communities.     

Hydroacoustic Estimates of Abundance and Spatial Distribution of Pelagic Prey Fishes in Western Lake Superior

Lake herring (Coregonus artedi) and rainbow smelt (Osmerus mordax) are a valuable prey resource for the recovering lake trout (Salvelinus namaycush) in Lake Superior. However, prey biomass may be insufficient to support the current predator demand. In August 1997, we assessed the abundance and spatial distribution of pelagic coregonines and rainbow smelt in western Lake Superior by combining a 120 kHz split beam acoustics system with midwater trawls. Coregonines comprised the majority of the midwater trawl catches and the length distributions for trawl caught fish coincided with estimated sizes of acoustic targets. Overall mean pelagic prey fish biomass was 15.56 kg ha⁻¹ with the greatest fish biomass occurring in the Apostle Islands region (27.98 kg ha⁻¹), followed by the Duluth Minnesota region (20.22 kg ha⁻¹), and with the lowest biomass occurring in the open waters of western Lake Superior (9.46 kg ha⁻¹). Biomass estimates from hydroacoustics were typically 2–134 times greater than estimates derived from spring bottom trawl surveys. Prey fish biomass for Lake Superior is about order of magnitude less than acoustic estimates for Lakes Michigan and Ontario. Discrepancies observed between bioenergetics-based estimates of predator consumption of coregonines and earlier coregonine biomass estimates may be accounted for by our hydroacoustic estimates.     

Results of the collaborative Lake Ontario bloater restoration stocking and assessment, 2012–2020

Bloater, Coregonus hoyi, are deepwater planktivores native to the Laurentian Great Lakes and Lake Nipigon. Interpretations of commercial fishery time series suggest they were common in Lake Ontario through the early 1900s but by the 1950s were no longer captured by commercial fishers. Annual bottom trawl surveys that began in 1978 and sampled extensively across putative bloater habitat only yielded one individual (1983), suggesting that the species had been locally extirpated. In 2012, a multiagency restoration program stocked bloater into Lake Ontario from gametes collected in Lake Michigan. From 2012 to 2020, 1,028,191 bloater were stocked into Lake Ontario. Bottom trawl surveys first detected stocked fish in 2015, and through 2020 ten bloater have been caught (total length mean = 129 mm, s.d. = 44 mm, range: 96–240 mm). Hatchery applied marks and genetic analyses confirmed the species identification and identified stocking location for some individuals. Trawl capture locations and acoustic telemetry suggested that stocked fish dispersed throughout the main lake within months or sooner, and the depth distribution of recaptured bloater was similar to historic distributions in Lake Ontario and other Great Lakes. Predicted bloater trawl catches, based on modeled population abundance and trawl survey efficiency, were similar to observed catches, suggesting that post-stocking survival is less than 20% and contemporary bottom trawl surveys can quantify bloater abundance at low densities and track restoration.     

Abundance and distribution of benthic macroinvertebrates in offshore soft sediments in Western Lake Huron, 2001–2007

Invasive species have had major impacts on the Great Lakes. This is especially true of exotic dreissenid mussels which are associated with decreased abundance of native macroinvertebrates and changes in food availability for fish. Beginning in 2001, we added a benthic macroinvertebrate survey to the USGS-Great Lakes Science Center's annual fall prey fish assessment of Lake Huron to monitor abundance of macrobenthos. Mean abundance of Diporeia, the most abundant benthic taxon in Lake Huron reported by previous investigators, declined greatly between 2001 and 2007. Diporeia was virtually absent at 27-m sites by 2001, decreased and was lost completely from 46-m depths by 2006, but remained present at reduced densities at 73-m sites. Dreissenids in our samples were almost entirely quagga mussels Dreissena bugensis. Zebra mussels Dreissena polymorpha were virtually absent from our samples, suggesting that they were confined to nearshore areas shallower than we sampled. Loss of Diporeia at individual sites was associated with arrival of quagga mussels, even when mussel densities were low. Quagga mussel density peaked during 2002, then decreased thereafter. During the study quagga mussels became established at most 46-m sites, but remained rare at 73-m sites. Length frequency distributions suggest that initial widespread recruitment may have occurred during 2001–2002. Like other Great Lakes, Lake Huron quagga mussels were associated with decreased abundance of native taxa, but negative effects occurred even though dreissenid densities were much lower. Dreissenid effects may extend well into deep oligotrophic habitats of Lake Huron.     

Impacts of the Introduced Round Goby (Apollonia melanostoma) on Dreissenids (Dreissena polymorpha and Dreissena bugensis) and on Macroinvertebrate Community between 2003 and 2006 in the Littoral Zone of Green Bay,Lake Michigan

We show that the invasion of round gobies (Apollonia melanostoma) in Green Bay, Lake Michigan, has changed the benthic food web in fundamental ways related to their impact on invasive dreissenid mussels. Dreissenid mussels are of specific interest because they are one of the primary dietary items for round gobies. In this study, we collected rocks from each of 10 study sites along approximately 60 km of the eastern shoreline of Green Bay, Lake Michigan, to assess a temporal change in macroinvertebrate abundance related to the northward movement of the round goby invasion front from a point about midway along the shoreline in 2003 to the entire coast in 2006. The pattern of macroinvertebrate abundance in 2003 suggested that round gobies had already caused significant decreases in macroinvertebrate abundances south of the invasion front (interpretation of the data could have been compromised by confounding environmental gradients). In subsequent sampling in 2006 macroinvertebrates were picked off of sampled rocks in the field and underwater transects were videotaped to estimate round goby abundance at each site. Round gobies were collected for stomach analysis to assist in determining which invertebrates would likely be impacted by goby predation. Our results indicated that by 2006, round gobies had become abundant at those sites where they were absent in 2003 and zebra mussels (Dreissena polymorpha), quagga mussels (Dreissena bugensis), isopods, amphipods, trichopterans, and gastropods in the newly invaded sites had significantly decreased at the newly invaded sites.     

Dreissena in Lake Ontario 30 years post-invasion

We examined three decades of changes in dreissenid populations in Lake Ontario and predation by round goby (Neogobius melanostomus). Dreissenids (almost exclusively quagga mussels, Dreissena rostriformis bugensis) peaked in 2003, 13 years after arrival, and then declined at depths <90 m but continued to increase deeper through 2018. Lake-wide density also increased from 2008 to 2018 along with average mussel lengths and lake-wide biomass, which reached an all-time high in 2018 (25.2 ± 3.3 g AFTDW/m²). Round goby densities were estimated at 4.2 fish/m² using videography at 10 to 35 m depth range in 2018. This density should impact mussel populations based on feeding rates, as indicated in the literature. While the abundance of 0–5 mm mussels appears to be high in all three years with measured length distributions (2008, 2013, 2018), the abundance of 5 to 12 mm dreissenids, the size range most commonly consumed by round goby, was low except at >90 m depths. Although the size distributions indicate that round goby is affecting mussel recruitment, we did not find a decline in dreissenid density in the nearshore and mid-depth ranges where goby have been abundant since 2005. The lake-wide densities and biomass of quagga mussels have increased over time, due to both the growth of individual mussels in the shallower depths, and a continuing increase in density at >90 m. Thus, the ecological effects of quagga mussels in Lake Ontario are likely to continue into the foreseeable future.     

Replacement of Zebra Mussels by Quagga Mussels in the Canadian Nearshore of Lake Ontario: the Importance of Substrate,Round Goby Abundance,and Upwelling Frequency

The invasion of the Great Lakes by zebra mussels (Dreissena polymorpha) and quagga mussels (Dreissena bugensis) has been accompanied by tremendous ecological change. In this paper we characterize the extent to which dreissenids dominate the nearshore of the Canadian shoreline of Lake Ontario and examine mussel distribution in relation to environmental factors. We surveyed 27 5-m sites and 25 20-m sites in late August 2003. Quagga mussels dominated all sites (mean: 9,404/m²; range 31–24,270), having almost completely replaced zebra mussels. Round gobies (Neogobius melanostomus) were associated with quagga populations dominated by large mussels. Quagga mussel total mass was low at 5-m sites with high upwelling frequency; we believe this is the first documentation of reduced benthic biomass in areas of upwelling in Lake Ontario. Overall, we estimated 6.32×10¹² quagga mussels weighing 8.13×10¹¹ g dry weight and carpeting ∼66% of the nearshore benthic habitat. Quagga mussels are a dominant and defining feature of the Lake Ontario nearshore, and must be accounted for in management planning.     

Twenty five years of changes in Dreissena spp. populations in Lake Erie

Lake Erie has the longest history of colonization by both Dreissena polymorpha and Dreissena rostriformis bugensis in North America and is therefore optimal for the study of long-term dynamics of dreissenid species. In addition, the morphometry of Lake Erie basins varies dramatically from the shallow western to the deep eastern basin, making this waterbody a convenient model to investigate patterns of Dreissena distribution, as well as interspecies interactions among dreissenids. We compare our data on the distribution, density and wet biomass of both dreissenid species in Lake Erie collected in 2009 and 2011–2012 with previous data. We found that Dreissena spp. distribution in Lake Erie varied depending on the time since the initial invasion, collection depth, and lake basin. In 2009–2012, zebra mussels were smaller than in 1992 and were consistently smaller than quagga mussels. During 2009–2012, quagga mussels were found at all depths and in all basins, while zebra mussels were common in the western basin only, and in the central and eastern basins were limited to shallow depths, resulting in an almost complete replacement of D. polymorpha with D. rostriformis bugensis. In the shallowest western basin of Lake Erie, zebra mussels represented > 30% of the combined dreissenid density even after more than 20 years of coexistence, providing strong evidence that, even in lakes as large as Lake Erie (or at least its western basin), D. polymorpha may sustain a significant presence for decades without being displaced by quagga mussels.     

Veliger density and environmental conditions control quagga mussel colonization rates in two perialpine lakes

Quagga and zebra mussels (Dreissena bugensis and D. polymorpha) are spreading across lakes in Europe and North America. In particular, quagga mussels colonize lakes to great depths (>200 m). To better understand the colonization pattern of quagga mussels in deep lakes, we studied the settlement of quagga mussels along a depth gradient on colonization plates at multiple depths (1–140 m) in the pelagic zone of two recently invaded perialpine lakes, Lake Constance and Lake Geneva. We measured colonization rates every three months over one year on colonization plates deployed in both lakes at defined depths. We also assessed long-term population dynamics from abundance and size distribution using repeated photogrammetry of colonization plates. Highest colonization rates and largest mussel sizes occurred above 8 m depth, and almost no zebra mussels were found. Colonization rates decreased to almost zero below 30 m. Colonization rates on plates were associated with variation in environmental conditions as well as veliger densities in the plankton across season and depth. Temperature was the most important environmental parameter that influenced colonization. Our results will help to better understand the seasonal colonization patterns of invasive quagga mussels in deep lakes.     

Pelagic Occurrence of Benthic Animals Near Shore in Lake Michigan

Collections with towed nets showed that presumably sedentary benthic animals occurred continually in the water above bottom in small numbers over depths of 6 m and 9 m in southeastern Lake Michigan. Densities were usually higher at night and in the months from June through September. Pelagic individuals at both depths consisted mainly of two chironomids, Saetheria tylus and Chironomus sp., and two naidids, Stylaria lacustris and Nais pardalis, but a total of 71 species were collected over a seven-month period. About 0.1% of the benthic animals were present above bottom at any time during a night in July 1974. The dominant pelagic species were the same as the benthic dominants at 6 m, but were much less important in the benthos at 9 m. We propose that this behavior, analogous to drifting of stream invertebrates, serves to maintain benthic populations in the frequently disturbed sediments of shallow (6 m deep or less) areas in Lake Michigan, and probably many other large lakes.     

Slimy sculpin depth shifts and habitat squeeze following the round goby invasion in the Laurentian Great Lakes

The collapse of Diporeia spp. and invasions of dreissenid mussels (zebra, Dreissena polymorpha; quagga, D. bugensis) and round goby (Neogobius melanostomus) have been associated with declines in abundance of native benthic fishes in the Great Lakes, including historically abundant slimy sculpin (Cottus cognatus). We hypothesized that as round goby colonized deeper habitat, slimy sculpin avoided habitat competition, predation, and aggression from round goby by shifting to deeper habitat. Accordingly, we predicted increased depth overlap of slimy sculpin with both round goby and deepwater sculpin (Myoxocephalus thompsonii) that resulted in habitat squeeze by both species. We used long-term bottom trawl data from Lakes Michigan, Huron, and Ontario to evaluate shifts in slimy sculpin depth and their depth overlap with round goby and deepwater sculpin. Lake Huron most supported our hypotheses as slimy sculpin shifted to deeper habitat coincident with the round goby invasion, and depth overlap between slimy sculpin and both species recently increased. Slimy sculpin depth trends in Lakes Michigan and Ontario suggest other ecological and environmental factors better predicted sculpin depth in these lakes.     

Circular Rep encoding single stranded (CRESS) DNA virus-like sequences detected in quagga mussels (Dreissena rostriformis bugensis) and sediments from the central Lake Michigan benthos

Circular Rep Encoding Single Stranded (CRESS) DNA viruses are a diverse group of viruses that have been identified in both terrestrial and aquatic ecosystems. Recent work in the Laurentian Great Lakes characterized the ecology and diversity of CRESS-DNA viruses associated with amphipods, Diporeia spp. In the last 20 years the Lake Michigan benthos has changed considerably with drastic population declines of Diporeia spp. concurrent with an increase in invasive quagga mussel (Dreissena rostriformis bugensis) abundance. The purpose of this study was to characterize CRESS-DNA virus-like elements (VLES; which could represent complete, partial/defective, endogenized or satellite viruses) associated with both invasive quagga mussels and in sediment collected in central Lake Michigan. Viral metagenomic libraries were prepared from two size classes (>25 mm and <15 mm shell length) of Lake Michigan quagga mussels and for two different sediment layers (136 mm and 290 mm below the lake floor) in a sediment core extracted from the lake. Viral metagenomes were different between quagga mussels and sediment cores. Nine VLE sequences were present in both the quagga mussel tissues and the sediment core layers analyzed. Cs¹³⁷ radiometric dating results indicate that these VLEs were present in the sediment prior to arrival of quagga mussels in Lake Michigan. These data suggest quagga mussels may interact with CRESS-DNA VLEs and algal DNA VLEs historically present in the Lake Michigan benthos. Overall, these data suggest that quagga mussels interacted with CRESS-DNA VLES present in Lake Michigan benthos since at least 1952.     

Impact of Zebra Mussels (Dreissena polymorpha) on Fingernail Clams (Sphaeriidae) in Extreme Southern Lake Michigan

The demography of the zebra mussel (Dreissena polymorpha) and its impacts on native fingernail clams (Sphaeriidae) in Lake Michigan near Michigan City, Indiana were studied from 1992 to 1997 at 5, 10, and 15 m depths. Zebra mussel densities ranged from 0 to 6,209/m², and were greatest at deeper stations. Size-frequency distributions suggest that initial colonization of the Michigan City area occurred in 1991, with adults typically living for 2 to 3 years. Abundance of adult populations may be limited by the sandy habitat typical of the areas studied and the occasional high-energy wave action in the shallow (5 m) near-shore zone of this area. Densities of fingernail clam ranged from 0 to 1,312/m² and were also greatest at deeper stations. Approximately half of the fingernail clams shells between 1 and 2 mm in size were used as substrate for attachment by juvenile and adult zebra mussels, while over 90% of the clams > 2 mm showed attachment by zebra mussels. Overall median densities of Sphaeriidae decreased from 832/m² to 13/m² at the 15 m depth and from 234/m² to 0/m² at the 10 m depth during the study. It appears that zebra mussel colonization caused a dramatic reduction of Sphaeriidae density by 1997 that may eventually result in their loss from the area.     

Measurement Error Associated With Surveys of Fish Abundance in Lake Michigan

In fisheries, imprecise measurements in catch data from surveys add uncertainty to the results of fishery stock assessments. The USGS Great Lakes Science Center (GLSC) began to survey the fall fish community of Lake Michigan in 1962 with bottom trawls. The measurement error was evaluated at the level of individual tows for nine fish species collected in this survey by applying a measurementerror regression model to replicated trawl data. It was found that the estimates of measurement-error variance ranged from 0.37 (deepwater sculpin, Myoxocephalus thompsoni) to 1.23 (alewife, Alosa pseudoharengus) on a logarithmic scale corresponding to a coefficient of variation = 66 to 156%. The estimates appeared to increase with the range of temperature occupied by the fish species. This association may be a result of the variability in the fall thermal structure of the lake. The estimates may also be influenced by other factors, such as pelagic behavior and schooling. Measurement error might be reduced by surveying the fish community during other seasons and/or by using additional technologies, such as acoustics. Measurement-error estimates should be considered when interpreting results of assessments that use abundance information from USGS-GLSC surveys of Lake Michigan and could be used if the survey design was altered. This study is the first to report estimates of measurement-error variance associated with this survey.     

Non-native Dreissena associated with increased native benthic community abundance with greater lake depth

Although the typical interaction between non-native invasive species and native species is considered to be negative, in some cases, non-native species may facilitate native species. Zebra and quagga mussels (Dreissena spp.) are aggressive invaders in freshwater systems, and they can alter energy flow by diverting nutrients from pelagic to benthic food-webs. In the last two decades, quagga mussels have largely replaced zebra mussels in shallow regions of the Laurentian Great Lakes and colonized deeper waters previously devoid of all dreissenids. Here, we aim to characterize potential positive effects of dreissenids in relation to depth on the benthic community in lakes Michigan and Huron. For this study, we used benthic survey data collected from Lake Michigan in 2015 and Lake Huron in 2017 and annual U.S. EPA Great Lakes National Program Office Long-term Biology Monitoring Program data for both lakes from 1998 to 2019. Benthic species richness and abundance (excluding dreissenids) in both lakes were almost three-fold higher in the nearshore (<70 m) compared to offshore (>70 m) communities. We found that, even though abundance of benthic invertebrates decreased with increased depth, total benthos density and biomass were higher in the presence than in the absence of quagga mussels in both lakes. Moreover, increased quagga mussel density and biomass with depth offset the lower benthos density and biomass at deeper depths, and samples with dreissenids had high densities of oligochaetes in both nearshore and offshore communities. These patterns are consistent with facilitative effects of quagga mussels on both shallow and deep-water benthic communities.     

Lake Champlain offshore benthic invertebrate community before and after zebra mussel invasion

Benthic invertebrates are important bio-indicators of water quality and play a significant role in aquatic systems. Lake Champlain has limited benthic invertebrate data which hinders development of food web models, assessment of invasive species impacts, and evaluation of management actions. In June 2016, we assessed benthic invertebrates along three transects in the main basin of Lake Champlain ranging from 5 to 100?m, and then compared results to densities from a limited survey in 1991 prior to the zebra mussel (Dreissena polymorpha) invasion. In 2016, total biomass and density were 1–2 orders of magnitude greater at 5?m than at 20–100?m. Zebra mussels, chironomids, oligochaetes, and gastropods were dominant at 5?m, and oligochaetes and sphaeriids were dominant at 20–100?m. Total density at the 5-m site was 94% lower in 2016 compared to 1991, but similar at the 100-m site. Diporeia, while abundant in many freshwater bodies, is historically rare in Lake Champlain and was not detected in our sampling. Because Lake Champlain benthic invertebrate densities are low and display dissimilar distributions to the Great Lakes, we hypothesize the offshore fish community is likely much more reliant on pelagic rather than benthic production. Although the current composition and biomass suggest the benthic community in Lake Champlain may not be greatly impacted by an invasion of quagga mussel (D. rostriformis bugensis), the potential for quaggas to re-route energy from pelagic to benthic habitats, as it has in the Great Lakes, could limit the Lake Champlain offshore fish community.     

White and longnose sucker diet composition in Lake Michigan and Saginaw Bay,Lake Huron

Sucker species (Catostomidae) are a common benthic invertivore in North American waterbodies yet are understudied in the Laurentian Great Lakes. Despite large biomass and potential competition with more economically valuable fish species, the diets of Great Lakes suckers are poorly described. We explored the gut contents of adult white suckers (Catostomus commersonii) in Lake Michigan and Saginaw Bay, Lake Huron, and longnose suckers (Catostomus catostomus) in Lake Michigan. Chironomidae was a primary prey for white suckers in both lakes, along with Amphipoda in Saginaw Bay. For longnose suckers in Lake Michigan, gut contents were dominated by Amphipoda and Isopoda. Relative to other benthic invertebrate taxa, white suckers positively selected Chironomidae as preferred prey (indexed via Chesson’s alpha). Moreover, the average length of consumed Chironomidae and Amphipoda increased with white sucker length, suggesting preferential size selection of prey as suckers grow. Although dreissenid mussels are overwhelmingly abundant in the benthos of Lake Michigan and Lake Huron, relatively few dreissenid mussels were consumed by either sucker species. Thus, suckers are unlikely to contribute to the control of these invasive invertebrates.