Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex.

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Intrinsic variability in shell and soft tissue growth of the freshwater mussel Lampsilis siliquoidea.

Freshwater mussels are ecologically and economically important members of many aquatic ecosystems, but are globally among the most imperiled taxa. Propagation techniques for mussels have been developed and used to boost declining and restore extirpated populations. Here we use a cohort of propagated mussels to estimate the intrinsic variability in size and growth rate of Lampsilis siliquoidea (a commonly propagated species). Understanding the magnitude and pattern of variation in data is critical to determining whether effects observed in nature or experimental treatments are likely to be important. The coefficient of variation (CV) of L. siliquoidea soft tissues (6.0%) was less than the CV of linear shell dimensions (25.1-66.9%). Size-weight relationships were best when mussel width (the maximum left-right dimension with both valves appressed) was used as a predictor, but 95% credible intervals on these predictions for soft tissues were ∼145 mg wide (about 50% of the mean soft tissue mass). Mussels in this study were treated identically, raised from a single cohort and yet variation in soft tissue mass at a particular size class (as determined by shell dimensions) was still high. High variability in mussel size is often acknowledged, but seldom discussed in the context of mussel conservation. High variability will influence the survival of stocked juvenile cohorts, may affect the ability to experimentally detect sublethal stressors and may lead to incongruities between the effects that mussels have on structure (via hard shells) and biogeochemical cycles (via soft tissue metabolism). Given their imperiled status and longevity, there is often reluctance to destructively sample unionid mussel soft tissues even in metabolic studies (e.g., studies of nutrient cycling). High intrinsic variability suggests that using shell dimensions (particularly shell length) as a response variable in studies of sublethal stressors or metabolic processes will make confident identifications of smaller effect sizes difficult.