Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf http://zoobank.org/urn:lsid:zoobank.org:pub:A758A612-C3B3-4C32-8EE8-986F8530706C COLONIZATION OF LAKE ERIE TRIBUTARIES BY ALLOCAPNIA RECTA (CAPNIIDAE) Alison L. Yasick, Robert A. Krebs1, and Julie A. Wolin Department of Biological, Geological and Environmental Sciences, Cleveland State University, Cleveland, OH 44115 U.S.A. E-mail: a.yasick@csuohio.edu Department of Biological, Geological and Environmental Sciences, Cleveland State University, Cleveland, OH 44115 U.S.A. E-mail: r.krebs@csuohio.edu Department of Biological, Geological and Environmental Sciences, Cleveland State University, Cleveland, OH 44115 U.S.A. E-mail: j.wolin@csuohio.edu 1corresponding author information: Robert A. Krebs, Department of Biological, Geological, and Environmental Sciences, Cleveland State University, 2121 Euclid Ave., Cleveland, Ohio 44115-2406 USA, r.krebs@csuohio.edu, 216-523-7553 ABSTRACT Plecoptera display a variety of dispersal mechanisms. The most reduced of these, walking, is common for winter stoneflies in the family Capniidae. To examine dispersal in a winter stonefly, Allocapnia recta, we sequenced a fragment of the 16s rRNA mitochondrial gene from four adjacent Lake Erie drainages. Nineteen unique haplotypes were obtained from 107 specimens. The four watersheds explained 37% of the haplotype diversity and all pairwise contrasts among streams were statistically significant; not one haplotype was common to all four watersheds. This extreme level of divergence among populations separated briefly in both space and time suggest minimal if any measurable gene flow and therefore the possibility that geological processes may better explain northern colonization to the Great Lakes than may adult migration. One possibility is for post-glacial change in topography, in particular isostatic adjustment to raise land, which led to flow reversal at the headwaters of rivers along what now is Ohio's Northern Divide, separating the Laurentian Great Lakes from the Ohio/Mississippi River system. Keywords: Allocapnia recta, distribution, haplotypes, migration, mtDNA, population structure INTRODUCTION Historical dispersal may be inferred by the presence of a species (Hynes 1988; Zwick 2000), but continued movement of individuals is critical to sustain genetic connectivity within populations, and that may be better measured from genetic markers (Avise 2000). Directly observing dispersal in aquatic insects is rarely practical (Miller et al. Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 41 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf 2009), thus genetic assessment and gene flow provides a useful proxy to assess dispersal in groups of smaller insect species (Finn & Adler 2006; Krosch et al. 2012) including Plecoptera (Kauwe et al. 2004; Schultheis et al. 2012; Pessino et al. 2014). Movement of aquatic invertebrates is expected to be both passive, by downstream drift, or by active flight, but is often limited by habitat (Finn et al. 2006). Flight in stoneflies can generally be limited by a continuum of extremes in wing size, ranging from full wing and full flight capabilities to apterous, flightless species. In the most rudimentary forms of adult flight, surface skimming may be the primary use of wings (Malmqvist 2000; Winterbourn 2005; Marden 2008). While the adult terrestrial stage in plecopterans is also short lived, some species, like Pteronarcys californica Newport (Pteronarcyidae), may possess a wide distribution, suggesting that duration of the adult stage is not necessarily a limit to dispersal (Kauwe et al. 2004). But, what of the many winter emerging species? The Capniidae presently compose 166 species, contributing 23% of the total North American fauna (Stark et al. 2014), and many have small or even rudimentary wings (Briers et al. 2003). A common requirement for cool streams, little silt and generally good water quality (Palma & Figueroa 2008) may further restrict their movement. Within the winter stoneflies, the distribution of A. recta (Claassen) is geographically broad (Stark et al. 2014), and the species is often locally abundant (Grubbs et al. 2012). Yasick et al. (2007) assessed potential physical limits on dispersal by comparing genetic variability between Allocapnia recta and Leuctra tenuis (Pictet), a long-winged summer emerging species (Stark et al. 1998). Two populations of L. tenuis were more genetically similar than those of A. recta, a pattern also found by Watanabe et al. (2010) in Trichoptera. Here we extend that work to test whether genetic variation is correlated most strongly with linear distance along waterways or direct distance overland between watersheds within the glaciated zone. An alternative hypothesis is that geographic distance is not correlated with genetic diversity, a possibility if each stream was colonized separately from the south after glacial retreat. MATERIALS AND METHODS Adult specimens of Allocapnia recta (Family Capniidae) were collected from small streams of four Lake Erie drainages (Fig. 1) in northeast Ohio between 2004 to 2007: 1. Rocky River. Summit County, East Branch Rocky River (41.2115, -81.6831); 2. Cuyahoga River. Summit County, tributary Boston Run 1 (41.2314, -81.5086) and Boston Run 2 (41.2335, -81.5021); 3. Chagrin River, Geauga County, tributary East Branch Chagrin River 1 (41.5960, -81.2512) and East Branch Chagrin River 2 (41.6071, -81.2875); 4. Grand River. Lake County, Talcott Creek (41.7217, -81.0830). Samples were obtained during the peak of annual winter adult emergence between November and February. The best collecting days occurred when temperatures exceeded 0°C with few to no clouds (Yasick, personal obs.). Specimens were collected from tree trunks, on snow and ice, or using a beating sheet on low-lying vegetation within 15 m of the stream. Individual A. recta were placed in a 1.5 ml microcentrifuge tube containing 95% ethanol. If a male and female were captured in copula, the mating pair was placed in the same tube. Specimens were identified to species based on the structure of the male genitalia according to Ross and Ricker (1971). The distal portions of the abdomen of specimens were retained as vouchers, while the remainder was used for DNA extraction. Only females captured in copula were used for DNA extraction. Insect preparation, mtDNA extraction, and amplification methods were conducted as in Yasick et al. (2007) with forward primer (SF_arF) 5'- TCG AAC AGA CCT AAA CTT TG -3' and the reverse primer (SF arR) 5'- AAT AAT TTA AAG TCT GAC CTG C -3'. Samples were sequenced at the Cleveland Clinic Genomic Core Facility initially from the forward primer, which gave a 467 base-pair sequence, reduced slightly from 492 bp in Yasick et al. (2007). Improved sequencing quality has generally eliminated the need for redundant sequencing of most samples, although chromatograms were assessed for any ambiguities by viewing them in Sequencher® (v. 4.10.1, Gene Codes Corp.). Any samples where the haplotype could not be assigned with certainty were sequenced in the reverse direction. Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 42 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf Fig. 1. Collecting sites of Allocapnia recta within the Rocky, Cuyahoga, Chagrin, and Grand rivers of northeast Ohio, USA. Inset figure is the Laurentian Great Lakes region of eastern North America, with a line indicating the Wisconsinan glacial maximum, approximately 15,000 years ago. Sampling adequacy was tested by rarefaction in PAST 3.0 (Hammer et al. 2001). Diversity indices (H and n) were calculated using Arlequin v. 3.5 (Excoffier & Lischer 2010). Haplotypes were defined and formatted in DnaSP 5.10 (Librado & Rozas 2009) to produce minimum spanning trees in Network v. 4.6 (Röhl 2004). To verify monophyly and species identification of all samples, a maximum likelihood phylogeny was constructed using MEGA V5 (Tamura et al. 2011) comparing A. recta to sequences of other Capniidae from Terry (2004) used with permission (GenBank accession numbers EF623127-31) and its sister family, Leuctridae (EF623173-81). Analysis of Molecular Variance (AMOVA) and pairwise Fst analyses were also run in ARLEQUIN, using 30,000 permutations for tests of significance. Hedrick (2005) and Merimans and Hedrick (2011) propose G'st as a standardized method of measuring genetic variation between highly variable populations with multiple alleles, as compared to Fst, because a G'st score of 1 indicates haplotypes frequencies are completely different, while a score of 0 is indicative of identical haplotypes. RESULTS Haplotypes were obtained from 107 specimens of A. recta from four watersheds in Northeast Ohio; the Grand (n = 25), Cuyahoga (n = 25; pooled from two sites), Chagrin (n = 27; pooled from two sites), and the Rocky River (n = 30). Nineteen haplotypes were identified across the four watersheds using a Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 43 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf Table 1. Haplotype frequencies for Allocapnia recta in the four watersheds, the Rocky, Cuyahoga, Chagrin and Grand rivers of northeast Ohio. Sites H1 H2 H3 H4 H5 H6 H7 H8 H9 H10 H12 H16 H17 H19 Minor Rocky 0 2 3 11 0 0 2 2 2 3 2 0 0 0 3 Cuyahoga 7 3 13 0 0 2 0 0 0 0 0 0 0 0 0 Chagrin 0 2 10 10 3 1 0 0 0 0 0 0 0 0 1 Grand 0 0 0 6 0 0 0 0 0 0 0 3 5 10 1 Sum 7 7 26 27 3 3 2 2 2 3 2 3 5 10 5 Frequency 0,07 0,07 0,24 0,25 0,03 0,03 0,02 0,02 0,02 0,03 0,02 0,03 0,05 0,09 0,05 fragment of the 16S rRNA gene (GenBank accession numbers, KC881036-KC881054). Given a small set of rare haplotypes, rarefaction curves (not shown) indicated sufficient sampling overall, with only one more haplotype expected to be found for each additional 20 individuals sequenced. Genetic variation, whether measured by the frequency of haplotypes or by n, the amount of nucleotide variation, was not significantly different among rivers (Table 2). Table 2. Sample size (N), haplotypes, unique haplotypes, gene diversity (H) and nucleotide diversity (n) estimated for the four populations of A. recta and overall in the region. SE is one standard error. Sites N H H (se) 71X100 7i (se) Rocky 30 0,85 0,06 0,37 0,18 Cuyahoga 25 0,66 0,07 0,19 0,12 Chagrin 27 0,73 0,07 0,25 0,14 Grand 25 0,76 0,08 0,32 0,17 Totals 107 0,86 0,02 0,47 0,18 While 14 of the 19 haplotypes occurred in multiple copies within the study, only four haplotypes were shared among watersheds, H2, H3, H4 and H6, and none occurred in all streams sampled (Table 1). The two most common haplotypes were represented by 26 and 27 individuals, respectively, but each was missing from one of the rivers, and each watershed possessed at least one of the 10 private haplotypes, those present only in one watershed and occurring in multiple copies: five in Rocky River, three in the Grand River and one each in the Cuyahoga and Chagrin Rivers. Thus, partitioning genetic differences among the four watersheds explained 37% (P <0.05) of the haplotype diversity (Table 3), and all pairwise contrasts were statistically significant (Table 4). A contrast of the two largest watersheds by area, the Cuyahoga and Grand Rivers, gave a G'st of 1.0, indicating complete isolation; the two most similar, yet significantly different watersheds, the Rocky and Chagrin Rivers, were not adjacent. Despite that variation could be significantly partitioned among these four rivers, all 19 haplotypes formed one large monophyletic clade. Geographically, the network supported separation of subgroups within watersheds (Fig. 2), and not just random isolation of different haplotypes. Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 44 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf Table 3. Molecular Analysis of Variance of 16s mtDNA for four populations of A. recta. Source of variation Df Sum of squares Variance component % Variation Among populations 3 58,9 0,69*** 36,7 Within population 103 122,9 1,19 63,3 Total 106 181,7 1,88 *** P< 0.001 Table 4. Pairwise population differences based on location from west to east in the region presenting both Fst (above the diagonal) and G'st (below the diagonal). Rocky Cuyahoga Chagrin Grand Rocky 0,47 *** 0,07 * 0,34 *** Cuyahoga 0,86 0,35 ** 0,62 *** Chagrin 0,31 0,45 0,29 *** Grand 0,68 1 0,72 * P< 0.05, ** P<0.01, *** P< 0.001 DISCUSSION The initial hypothesis stated that dispersal distance drives patterns of intraspecific variation in A. recta across neighboring watersheds. That expectation was only partly supported, because whether watersheds were adjacent or distant (up to 75 km away), levels of differences were similarly large. No single haplotype was found in all four watersheds, and haplotypes from the two largest rivers were completely different (G'ST of 1.0). Thus, little connectivity among populations, or conversely, little gene flow, even between neighboring watersheds, was observed. Concurrently, similar levels of genetic diversity were present in each of the four rivers, but this diversity arises from different sets of haplotypes. One plausible cause is that many haplotypes arrived during initial colonization events within the region, and genetic drift has led to some haplotype losses. The Capniidae lineage appears very diverse for haplotypes, yet members of the family present some of the shortest 16S rRNA gene branch lengths (Terry, 2004). That result is consistent with high levels of within-species variation and incomplete lineage sorting, leading to poor differentiation among species for any one gene (Kutschera et al. 2014). Such an outcome is theoretically more likely in lineages where population sizes have been historically large (van Velzen et al. 2012). At present, migration in A. recta appears to be slow, as males of A. recta are brachypterous, the species emerges in winter, and females tend to remain in the natal stream to oviposit (Malmquist 2000; Alp et al. 2012). Furthermore, the current distribution of A. recta regionally suggests that this is not a vagile species. This stonefly inhabits cool streams from the highlands of Alabama and Georgia, northeast to Maine, largely within the range of the Appalachian Mountains (Earle 2009). From this historic base, A. recta has only sporadically dispersed into regions once covered by the Wisconsin glacial lobes (Ross & Ricker 1971); the species does not occur in Michigan (Stark et al. 1986; Grubbs et al. 2012) and inhabits mainly the unglaciated regions of Illinois and Indiana Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 45 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf Fig. 2. A network of 16S mt rRNA haplotypes for Allocapnia recta collected in four watersheds of NE Ohio, USA. The diameter of each circle is proportional to haplotype frequency and levels of shading denote location increasing in shading level from west to east: Rocky River (white), Cuyahoga River (light grey), Chagrin River (dark grey) and Grand River (black). Numbers denote haplotypes and hash marks indicate base changes where more than one nucleotide difference separates the connected haplotypes. (DeWalt & Grubbs 2011) where a dispersal route was formed by the Ohio River watershed. This distributional limit is important for considering colonization into the Laurentian Great Lakes region. One important path for many aquatic species after the last glacial retreat 15,000 years ago was the Wabash River (Teller 2004), which drains extensive plains before entering the Ohio River near its confluence with the Mississippi River. Prior to the opening of the St. Lawrence river, a developing Lake Erie would have covered the lake plains of NW Ohio (White & Totten 1982). After a northern outlet opened, flow reversed to create a slow warm-water river, the modern Maumee (Lewis et al. 2012), the largest tributary of Lake Erie. This environment favored migration of amphibians, fish, and molluscs into Lake Erie (Mandrak & Crossman 1992; Austin et al. 2002; Graf 2002, Borden & Krebs 2009), but failed to provide habitat for silt sensitive cold-water species like A. recta (Meador & Goldstein 2003; Allan 2004). A failure to collect this species along the Maumee River or in the lake itself (Grubbs et al. 2012; DeWalt et al. 2012) is suggestive of a physiological barrier that limits or slows dispersal. That cool streams flowing from the escarpment of northeastern Ohio provide habitat (Szabo et al. 1988; DeWalt et al. 2012) may implicate Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 46 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf temperature effects, but several related species, like A. granulata, A. pygmaea and A. minima also colonized glacial areas (Frison 1942; DeWalt & Grubbs 2011; Stark et al. 2014). The extensive southern distribution of the species in highlands [note the removal of A. recta from taxa lists for Mississippi (Nations et al. 2007)] may suggest that sediment composition is a more important factor (Pessino et al. 2014). Ross & Ricker (1971) recorded A. recta along the boundary of the once-glaciated region of Ohio, and DeWalt et al. (2012) extended records to include the counties along Lake Erie for which we report genetic data, as well as the neighboring counties to the south. These provide a link for A. recta to the Lake Erie watershed by stream capture following isostatic rebounding (Coffey 1958; Bishop 1995), which was also proposed to explain the disjunct distribution of genetic variants above waterfalls in NE Ohio streams for two unionid mussels, Pyganodon grandis (Krebs 2004) and Lampsilis siliquoidea (Krebs et al. 2013), as well as smallmouth bass, Micropterus dolomieu (Borden & Krebs 2009). The high levels of mtDNA haplotype diversity strengthens support for a geological cause of observed patterns rather than a biological one, because waves of rare adult dispersal are likely to reduce variation to few or even single haplotypes in founding populations (Theissinger et al. 2013; Krebs et al. 2015). As the glacial ice melted, re-colonization by A. recta may have occurred first into the Rocky River, and expanded eastward. This post-glacial population dynamic may explain why specimens of A. recta collected from the Rocky River are the most diverse, and why only some haplotypes are shared among watersheds. Similar limits are predicted for caddisflies (Lehrian et al. 2010). By contrast, a long-winged warm-weather species like Acroneuria frisoni shows low levels of divergence in nearby watersheds in comparisons made across its range (Pessino et al. 2014), as similarly reported for Leuctra tenuis between the Chagrin and Grand rivers (Yasick et al. 2007). However, all three of these species, A. frisoni, L. tenuis and A. recta, are characterized by high haplotype diversity in the region north of the Ohio River. Variation within streams may historically have been even higher, with genetic drift creating the differences observed today. Current gene flow in A. recta remains minimal, while levels of movement even along streams remain unknown. A prime region in which to extend this research are streams draining south from Ohio's northern divide. ACKNOWLEDGEMENTS We thank Holden Arboretum, Lake County Metroparks, and the Cuyahoga Valley National Park for permission to collect on public land. This work was made possible by funding from the Shoemaker Foundation. REFERENCES Allan, J.D. 2004. Landscapes and riverscapes: the influence of land use on stream ecosystems. Annual Review of Ecology, Evolution, and Systematics, 35:257-284. Alp, M., I. Keller, A.M. Westram, & C. Robinson. 2012. How river structure and biological traits influence gene flow: a population genetic study of two stream invertebrates with differing dispersal abilities. Freshwater Biology, 57:969981. Austin, J.D., S.C. Lougheed, L. Neidrauer, A.A. Chek, & P.T. Boag. 2002. Cryptic lineages in a small frog: the post-glacial history of the spring peeper, Pseudacris crucifer (Anura: Hylidae). Molecular Phylogenetics and Evolution, 25:316329. Avise, J.C. 2000. Phylogeography: The History and Formation of Species. Harvard University Press, Cambridge, Massachusetts. Bishop, P. 1995. Drainage rearrangement by river capture, beheading and diversion. Progress in Physical Geography, 19:449-473. Borden, W.C. & R.A. Krebs. 2009. Phylogeography and postglacial dispersal of smallmouth bass (Micropterus dolomieu, Centrarchidae) into the Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences, 66:2142-2156. Briers, R.A., H.M. Cariss, & J.H.R. Gee. 2003. Flight activity of adult stoneflies in relation to weather. Ecological Entomology, 28:21-40. Coffey, G.N. 1958. Major glacial drainage changes in Ohio. The Ohio Journal of Science, 58:43-49. DeWalt, R.E. & S.A. Grubbs. 2011. Updates to the Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 47 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf stonefly fauna of Illinois and Indiana. Illiesia, 7 (3):31-50. Available from: http://www2.pms-lj.si/illiesia/papers/Illiesia07-03.pdf. DeWalt, R.E., Y. Cao, T. Tweddale, S.A. Grubbs, L. Hinz, M. Pessino, & J.L. Robinson. 2012. Ohio USA stoneflies (Insecta, Plecoptera): species richness estimation, distribution of functional niche traits, drainage affiliations, and relationships to other states. ZooKeys, 178:1-26. Earle, J.I. 2009. New state stonefly (Plecoptera) records for Pennsylvania, with additional records and information on rare species. Illiesia, 5:169-181. Available online: http://www2.pms-lj.si/illiesia/papers/Illiesia05-16.pdf. Excoffier, L. & H.E.L. Lischer. 2010. Arlequin suite V. 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 10:564567. Finn, D.S. & P.H. Adler. 2006. Population genetic structure of a rare high-elevation black fly, Metacnephia coloradensis, occupying Colorado lake outlet streams. Freshwater Biology, 51:2240-2251. Finn, D.S., D.M. Theobald, W.C. Black IV, & N.L. Poff. 2006. Spatial population genetic structure and limited dispersal in a Rocky Mountain alpine stream insect. Molecular Ecology, 15:3553-3566. Frison, T.H. 1942. Studies of North American Plecoptera with special reference to the fauna of Illinois, 22:235-350. Graf, D.L. 2002. Historical biogeography and late glacial origin of the freshwater pearly mussel (Bivalvia: Unionidae) faunas of Lake Erie. North American Museum of Comparative Zoology, Occasional Papers on Molluscs, 6:175-210. Grubbs, S.A., M. Pessino, & R.E. DeWalt. 2012. Michigan Plecoptera (Stoneflies): distribution patterns and an updated species list. Illiesia, 8:162-173. Available online: http://www2.pms-lj.si/illiesia/papers/Illiesia08-18.pdf Hammer, 0., D.A.T. Harper, & P.D. Ryan. 2001. PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica 4, 9 Pp. (Available from: http://folk.uio.no/ohammer/past/). Hedrick, P.W. 2005. A standardized genetic differentiation measure. Evolution. 59:16331638. Hynes, H.B.N. 1988. Biogeography and origins of the North American stoneflies (Plecoptera). Memoirs of the Entomological Society of Canada, 144:31-37. Kauwe, J.S., D.K. Shiozawa, & R.P. Evans. 2004. Phylogeographic and nested clade analysis of the stonefly Pteronarcys californica (Plecoptera: Pteronarcyidae) in the western USA. Journal of the North American Benthological Society, 23:824-838. Krebs, R.A. 2004. Combining paternally and maternally inherited mitochondrial DNA for analysis of population structure in mussels. Molecular Ecology, 13:1701-1705. Krebs, R.A., W.C. Borden, N.M. Evans, & F.P. Doerder. 2013. Differences in population structure between maternally- and paternally-inherited forms of mitochondria in Lampsilis siliquoidea (Bivalvia: Unionidae). Biological Journal of the Linnean Society, 109:229-240. Krebs, R.A., B.D. Allen, N.M. Evans, and D.T. Zanatta. 2015. Mitochondrial DNA structure of Pyganodon grandis (Bivalvia: Unionidae) from the Lake Erie watershed and selected locations in its northern distribution. American Malacological Bulletin, (in press). Krosch, M.N., A.M. Baker, P.B. Mather, & P.S. Cranston. 2012. Comparison of intraspecific genetic structure among related chironomids (Diptera) from New Zealand and Patagonia: disparity between potential and realized dispersal. Freshwater Science, 31:1105-1120. Kutschera, V.E., T. Bidon, F. Hailer, J.L. Rodi, S.R. Fain, & A. Janke. 2014. Bears in a forest of gene trees: phylogenetic inference is complicated by incomplete lineage sorting and gene flow. Molecular Biology and Evolution, 31:2004-2017. Lehrian, S., M. Balint, P. Haase, & S.U. Pauls. 2010. Genetic population structure of an autumn-emerging caddisfly with inherently low dispersal capacity and insights into its phylogeography. Journal of the North American Benthological Society, 29:1100-1118. Lewis, C.F.M., G.D.M. Cameron, T.W. Anderson, C.W. Heil Jr., & P.L. Gareau. 2012. Lake levels in the Erie Basin of the Laurentian Great Lakes. Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 48 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf Journal of Paleolimnology, 47:493-511. Librado, P. & J. Rozas. 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25:14511452. Malmqvist, B. 2000. How does wing length relate to distribution patterns of stoneflies (Plecoptera) and mayflies (Ephemeroptera)? Biological Conservation, 93:271-276. Mandrak, N.E. & E.J. Crossman. 1992. Postglacial dispersal of freshwater fishes into Ontario. Canadian Journal of Zoology, 70:2247-2259. Marden, J.H. 2008. Quantitative and evolutionary biology of alternative splicing: how changing the mix of alternative transcripts affects phenotypic plasticity and reaction norms. Heredity, 100:111-120. Meador, M.R. & R.M. Goldstein. 2003. Assessing water quality at large geographic scales: relations among land use, water physicochemistry, riparian condition, and fish community structure. Environmental Management, 31:504-517. Merimans, P.G. & P.W. Hedrick. 2011. Assessing population structure: Fst and related measures. Molecular Ecology Resources, 11:5-18. Miller, M.P., D.W. Blinn, & P. Keim. 2009. Correlations between observed dispersal capabilities and patterns of genetic differentiation in populations of four aquatic insect species from the Arizona White Mountains, USA. Freshwater Biology, 47:16601673. Nations, T.M., B.P. Stark, & M.B. Hicks 2007. The Winter Stoneflies (Plecoptera: Capniidae) of Mississippi. Illiesia, 3:70-94. Available online: http://www2.pms-lj.si/illiesia/Illiesia03-09.pdf Palma, A. & R. Figueroa. 2008. Latitudinal diversity of Plecoptera (Insecta) on local and global scales. Illiesia, 4:81-90. Available online: http://www2.pms-lj.si/illiesia/Illiesia04-08.pdf Pessino, M., E.T. Chabot, R. Giordano, & R.E. DeWalt. 2014. Refugia and postglacial expansion of Acroneuria frisoni Stark & Brown (Plecoptera: Perlidae) in North America. Freshwater Science, 33:232-249. DOI: 10.1086/675306 Ross H.H. & W.E. Ricker. 1971. The classification, evolution, and dispersal of the winter stonefly genus Allocapnia. Illinois Biological Monographs, University of Illinois Press, Urbana, Illinois. Available at https://www.ideals.illinois.edu/handle/2142/27332. Röhl, A. 2004. NETWORK version 4.5.1.0. Shareware phylogenetic network software presented by Fluxus Technology Ltd. Available at http://www.fluxus-engineering.com. April 2013. Schultheis, A.S., J.Y. Booth, L.R. Perlmutter, J.E. Bond, & A.L. Sheldon. 2012. Phylogeography and species biogeography of montane Great Basin stoneflies. Molecular Ecology, 21:33253340. Stark, B.P., S.W. Szczytko, & R.W. Baumann.1986. North American Stoneflies (Plecoptera): systematics, distribution, and taxonomic references. The Great Basin Naturalist, 46:383397. Stark, B.P., S.W. Szczytko, and C.R. Nelson. 1998. American Stoneflies: A Photographic Guide to Plecoptera. The Caddis Press, Columbus, Ohio. Stark, B.P., R.W. Baumann, M. Maehr, & R.E. DeWalt. 2014. Valid Stonefly Names for North America. Plecoptera Species File Online, Version 1.0/4.1. 21 October 2014. Available at http://plecoptera.speciesfile.org Szabo, J.P., C.H. Carter, P.W. Bruno, & E.J. Jones. 1988. Glacial and post-glacial deposits of Northeast Ohio. The Ohio Journal of Science, 88:66-74. Tamura, K., D. Peterson, N. Peterson, G. Stecher, M. Nei, & S. Kumar. 2011. MEGA 5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Molecular Biology and Evolution, 28:2731-2739. Teller, J.T. 2004. Controls, history, outbursts, and impact of large late-Quaternary proglacial lakes in North America. Developments in Quaternary Science. Pages 45-61. In A.R. Gillespie, S.C. Porter & B.F. Atwater [editors]. The Quaternary Period in the United States. INQUA Anniversary Volume, Elsevier, New York. Terry, M.D. 2004. Phylogeny of the polyneopterous insects with emphasis on Plecoptera: molecular and morphological evidence. unpublished Ph.D. dissertation, Brigham Young University, Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 49 Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2015. Colonization of Lake Erie tributaries by Allocapnia recta (Capniidae). Illiesia, 11(05):41-50. Available online: http://illiesia.speciesfile.org/papers/Illiesia11-05.pdf Provo, Utah. Available at http://contentdm.lib.byu.edu/ ETD/image/etd376.pdf Theissinger, K., M. Balint, K.A. Feldheim, P. Haase, J. Johannesen, I. Laube, & S.U. Pauls. 2013. Glacial survival and post-glacial recolonization of an arctic-alpine freshwater insect (Arcynopteryx dichroa, Plecoptera, Perlodidae) in Europe. Journal of Biogeography, 40:236-248. van Velzen, R.,E. Weitschek, G. Felici, & F.T. Bakker. 2012. DNA barcoding of recently diverged species: relative performance of matching methods. PLOS ONE 7:e30490. Available at http: 10.1371/journal.pone.0030490 Watanabe, K., M.T. Monaghan, Y. Takemon, & T. Omura, 2010 Dispersal ability determines the genetic effects of habitat fragmentation in three species of aquatic insect. Aquatic Conservation: Marine and Freshwater Ecosystems, 20:574-579. White, G.W. & S.M. Totten. 1982. Glacial geology of northeastern Ohio. Ohio Geological Survey Bulletin, 68:1-75. Winterbourn, M.J. 2005. Dispersal, feeding, and parasitism of adult stoneflies (Plecoptera) at a New Zealand forest stream. Aquatic Insects, 27:155-166. Yasick, A.L., R.A. Krebs, and J.A. Wolin. 2007. The effects of dispersal ability in winter and summer stoneflies on their genetic differentiation. Ecological Entomology, 32:399404. Zwick, P. 2000. Phylogenetic system and zoogeography of the Plecoptera. Annual Review of Entomology, 45:709-746. Received 10 September 2014, Accepted 15 January 2015, Published 10 February 2015 Illiesia - http://illiesia.speciesfile.org Volume 11 - Number 5 - Page 50