## 5 Conclusions
(b) Tree inferred exactly as for (a) with a constraint tree representing the known morphospecies imposed on the search space; log-likelihood = \(-\)1184.068. (c) Maximum clade credibility tree inferred using *BEAST, with models of site evolution estimated simultaneously to tree estimation using bModefTest and a birth-death process for the tree prior; 95% HPDI of tree log-likelihood = \([-\)1123.086, \(-\)1102.800], mean log-likelihood = \(-\)1106.792. Table S1: Location of field sites from which adult Dermaptera were collected. Numbers refer to labels on Figure S1. Table S2: Summary statistics for morphometric variables for the Anisolabididae morphospecies. \"Forceps asymmetry\" is an index calculated as the sum of the absolute differences between the heights and widths of a specimen's two forceps, \"Parameter ratio\" is the ratio of a specimen's right parameter's length to its width, and \"Forceps ratio\" is the ratio of the maximum forceps length to the maximum forceps width. These three measurements were dimensionless; all others were in mm. Individual cells show the mean, the standard deviation in brackets, and the minimum and maximum values below them. N refers to the number of individuals of each morphospecies measured. Note that only the three dimensionless variables plus the basal width of the forceps were used for the morphometric and classification analysis. Table S3: Details of linear discriminant function analysis. \"Classification\" refers to the grouping variable used (i.e., \"Genus\" refers to the discrimination of the two gene, \"_Anisolabis_\" to that between _Anisolabis_ spp., etc.). Function refers to the individual discriminant functions generated, of which there can be a maximum of k \(-\) 1, where k is the number of levels of the grouping variable (morphospecies or genera in this case). Functions are numbered in order of the proportion of variance explained by each, which is shown in brackets. \"Forceps asymmetry\" is an index calculated as the sum of the absolute differences between the heights and widths of a specimen's two forceps, \"Paramere ratio\" is the ratio of a specimen's right parameter's length to its width, and \"Forceps ratio\" is the ratio of the maximum forceps length to the maximum forceps width. These three measurements were dimensionless; the basal width of the forceps was in mm. Table S4: Uncorrected p-distances within and between earwigs species on the Barcode of Life Data Systems public repository's cytochrome oxidase c subunit 1 barcodes (507-680 bp). Only records with species names attached and with sequences > 500 bp were retrieved. Interspecific means and variances are shown on the left and right of the diagonal. Intraspecific mean shown on diagonal with variance in brackets. Variances were calculated using 1000 bootstraps. N refers to total individuals retrieved. Forficula auricularia A and B refers to cryptic sibling species [58]. **Author Contributions:** Conceptualization, O.P.S., M.B., P.A.U. and A.A.H. ; methodology, O.P.S., M.B., P.A.U., J.H., D.S., M.N., T.H., M.v.H. and A.A.H. ; formal analysis, O.P.S. ; investigation, O.P.S., M.B., P.A.U., J.H., D.S., M.N., T.H., M.v.H. and A.A.H. ; resources, P.A.U., J.H., D.S., M.N., T.H., M.v.H. and A.A.H. ; data curation, O.P.S. ; writing--original draft preparation, O.P.S. ; writing--review and editing, O.P.S., M.B., P.A.U., J.H., D.S., M.N. and A.A.H. ; visualization, O.P.S. ; supervision, M.B., P.A.U. and A.A.H. ; project administration, P.A.U., J.H., D.S., M.N., T.H., M.v.H. and A.A.H. ; funding acquisition, P.A.U., J.H., D.S., M.N., T.H., M.V.H. and A.A.H. **Funding:** This research was funded by the Grains Research and Development Corporation, grant number CSE00059. The APC was funded by the Pest and Environmental Adaptation Research Group. **Acknowledgments:** Gratitude is extended to the property owners who granted permission to collect invertebrate specimens. We acknowledge [PERSON] for his assistance with identifying morphogenesis and navigating the literature.