Hierarchical Clustering Analysis (HCA) for the imputed (N = 105) dataset using ‘hclust’ command in R v3.5.1. Euclidean distances among the specimens were calculated using the ‘dist’ command in R. In both cases the normalized data (H/D%, T/D%, U/D% and H/T) were used for Euclidean distance calculation. Note the numbers are the last numbers of specimen number in the text (e.g., 1 for KSKV2019Jara/1), followed by Type 1, 2, 3, or 4 after a forward slash, which is followed by ‘m’ for microconch or ‘M’ for macroconch.