Variation in shell shape and mitochondrial DNA of two species of freshwater mussels in the Kiamichi River, Oklahoma, USA
Data files
Jul 06, 2026 version files 564.90 KB
-
Moore_Vaughn_Franzen_Kiamichi_Morphology.zip
560.12 KB
-
README.md
4.78 KB
Abstract
Numerous freshwater mussel species possess highly variable shell morphologies, though it remains unclear what major factors drive shell morphology patterns in many species. Therefore, basic research of the interactions between mussels and their environment is a pressing research priority. We investigated the shell shape of two mussel species found in the Kiamichi River, USA, to test the hypothesis that river hydrology gradients have an influence on shell traits rather than the alternative that shell trait variation is strictly a product of evolutionary processes. We measured shell length, height, and width of Amblema plicata and Pustulosa pustulosa individuals in the river at four sites, then extracted and sequenced mitochondrial DNA (mtDNA) from each individual to assess genetic variation. After controlling for shell size, we tested whether shell traits varied between sites and varied with hydrological attributes. We constructed phylogenies and generated TCS haplotype networks from mtDNA sequences to examine the distribution of genetic diversity of A. plicata and P. pustulosa in the river and used Mantel tests to test for isolation-by-distance patterns. We found little difference in A. plicata shell shape across sites, whereas P. pustulosa individuals were more compressed in the upper river, likely due to higher stream velocity, and P. pustulosa at downriver sites were more inflated, likely due to increased flow rate. Haplotype networks indicated extensive sharing of mtDNA haplotypes between sites for both species. Our phylogenies showed no indication Kiamichi River A. plicata and P. pustulosa populations are genetically unique, nor did we find evidence of isolation by distance, which suggests shell shape is unlikely to be driven by genetic distinctiveness. Our results show that hydrology likely has a strong influence on shell shape in some mussels, and considering increases in extreme droughts and floods, this study highlights the importance of basic ecological research for informing conservation.
Data for Moore, Vaughn, and Franzen (in revision with journal - Freshwater Mollusk Biology and Conservation)
Description of data
File list
File: Moore_Vaughn_Franzen_Kiamichi_Morphology.zip
README.txt
README file
Haplotype_Network
-
Amblema_network.nex
Alignment used to generate haplotype network in PopArt. Plot was subsequently edited in Illustrator for final touches.
-
Pustulosa_network.nex
Alignment used to generate haplotype network in PopArt. Plot was subsequently edited in Illustrator for final touches.
Hydrology
-
EAM_project_kiamichi_dists.csv
Coordinates for each sampling site and river source in CSV format. site=sampling site; latitude=latitude; longitude=longitude
-
EAM_project_LWH_20250623.csv
List of specimens and their collection localities in CSV format. site=sampling site; latitude=latitude; longitude=longitude; species=species; individual=unique specimen identification number; Length_RV (mm)= max length of right valve in mm Height_RV (mm)= max height of right valve in mm; Width_RV (mm)= max width of right valve in mm
-
kiamichi_distmatrix.csv
River distance matrix. K2=K1, KT=K2, KS=K3, K7=K4, source=source. See EAM_project_all_morphology_habitat_data.csv for additional clarification.
-
Kiamichi_HUC8
Folder containing edited Kiamichi River shapefile.
-
kiamichi_hydro.R
R script used to generate hydrological data for each sampling locality and river distance matrix.
-
kiamichi_net_cleanup.Rdata
R data file generated after using the nhdplustools "cleanup" function.
IBD
-
Apli_pdist.csv
Genetic distance matrix for Amblema plicata in CSV format.
-
IBD.R
R script used to conduct isolation by distance analysis (Mantel tests).
-
kiamichi_distmatrix_Mantel.csv
Geographic distance matrix in CSV format.
-
Ppus_pdist.csv
Genetic distance matrix for Pustulosa pustulosa in csv format.
IQTREE
-
Amblema_plicata_seqs.fasta
Alignment for Amblema plicata phylogenetic analysis.
-
partition.txt
Partition file used for both phylogenetic analyses.
-
Pustulosa_pustulosa_seqs.fasta
Alignment for Pustulosa pustulosa phylogenetic analysis.
Morphometrics
-
EAM_project_all_morphology_habitat_data.csv
Metadata for measured specimens, raw linear measurements, and hydrological data in CSV format.
Description of variables: site=sampling site identification; common_site=alternative site identification; latitude=latitude; longitude=longitude; species=species identification; individual=specimen identification number; Length_RV (mm)= max length of right valve in mm; Height_RV (mm)= max height of right valve in mm; Width_RV (mm)= max width of right valve in mm; comid=the unique, permanent integer used to identify specific water features (like flowlines, catchments, and waterbodies) within the USGS National Hydrography Dataset (NHD) and the NHDPlus geospatial framework; gnis_id=a permanent, unique 10-digit identifier assigned by the USGS to official geographic features; gnis_name=USGS name; lengthkm=length of flowline in kilometers; reachcode=Unique reach identifier; streamleve=Stream level; streamorde=Strahler stream order; totdasqkm=Total cumulative area, in square kilometers; divdasqkm=Divergence-routed cumulative area, in square kilometers; slope=Slope of the flowline from smoothed elevation (unitless); slopelenkm=Flow-line length used to calculate slope, in kilometers; qa_ma=Mean annual flow from runoff; qc_ma=Mean annual flow with reference-gage regression; qe_ma=Mean annual flow from gage adjustment; ve_ma=Velocity from gage adjustment; qe_01=Mean annual flow from gage adjustment for January; qe_02=Mean annual flow from gage adjustment for February; qe_03=Mean annual flow from gage adjustment for March; qe_04=Mean annual flow from gage adjustment for April; qe_05=Mean annual flow from gage adjustment for May; qe_06=Mean annual flow from gage adjustment for June; qe_07=Mean annual flow from gage adjustment for July; qe_08=Mean annual flow from gage adjustment for August; qe_09=Mean annual flow from gage adjustment for September; qe_10=Mean annual flow from gage adjustment for October; qe_11=Mean annual flow from gage adjustment for November; qe_12=Mean annual flow from gage adjustment for December; dist_from_source=distance of point from the river source in kilometers (river kilometers)
-
EAM_project_morphometrics.R
R script used to conduct morphometric analyses. All plots were subsequently edited in Illustrator for final touches.
Access information
Other publicly accessible locations of data:
- NCBI GenBank: Accession
