Genomic signals found using RNA sequencing support conservation of walleye (Sander vitreus) in a large freshwater ecosystem.

Abstract

RNA sequencing is an effective approach for studying an aquatic species with little prior molecular information available, yielding both physiological and genomic data, but its genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced messenger RNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (FST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations, and suggest genomic differences possibly underlying observed phenotypic differences. Because RNA sequencing data can yield physiological in addition to genetic information discussed here, we argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.