Normalized HMW DNA yield (nY) was calculated as follows using data obtained from the Agilent Technologies TapeStation 2200 DNA Analyzer and extract tissue weights modified by the correction ratio explained above:
Data for each taxon, time interval and response variable (%R and nY) were analyzed separately. A Shapiro-Wilk normality test was used to determine whether data were normally distributed. For each taxon, the effect of storage solution on DNA quality was assessed using a repeated measures design: normally distributed data were analyzed using a parametric repeated measures ANOVA and non-normal data were analyzed using a non-parametric Friedman χ2 test (RStudio v. 1.2.1335). Post-hoc tests for repeated measures ANOVAs were performed with the ‘nlme’ package in RStudio (Version 3.1–137). We used the Friedman Conover test as a post-hoc test for Friedman χ2 analyses and performed them with a Bonferroni correction in the ‘PMCMR’ package (Version 4.3).