## Materials and Methods
### Phylogeny  
The generic classifications after [PERSON] [6] were tested by gathering more sequence data. We obtained a 358 bp fragment of cytochrome \(b\) (_cob_), a 442 bp fragment of 28S ribosomal DNA (_28S_), and a 466 bp fragment of the tubulin alpha 1 gene (_tuba1_) for selected individuals from the Anisolabididae morphospecies. For each gene fragment, we aimed to obtain at least two individuals for each morphospecies. Cytochrome \(b\) was amplified and sequenced under the same reaction conditions and cycling profile as _cox1_, using the CB3/CB4 primer pair [73]. The _28S_ fragment was amplified with primers designed using Primer3 [74; 75] implemented in Geneious using the default settings. Suitable primer binding sites were identified for a fragment > 400 bp with a multiple sequence alignment of four Anisolabididae _28S_ sequences from [PERSON] [76] in GenBank [77] (accession numbers: AB119545, AB119544, AB119542, AB119549). Two forward (An16F: GAGAATCCGAATATCTGAA and An17F: AGAAAATCCGAATCTGAAG) and two reverse primers (An481R: AATAAATTGCCAACAAATGC and An485R: TGAGAATAAATTGCCAACA) were selected and tested, amplifying a fragment 452-462 bp. PCR was performed using each possible pair, and the combination that produced the strongest bands after staining and electrophoresis was selected for sequencing. _28S_ was amplified with a cycle of 93 \({}^{\circ}\)C denaturation for three minutes, followed by 30 cycles of 93 \({}^{\circ}\)C denaturation for 60 s, 55 \({}^{\circ}\)C annealing for 90 s, and 72 \({}^{\circ}\)C extension for 75 s, and then 72 \({}^{\circ}\)C for four minutes for the final extension. For _tuba1_, the primer pair DDVtabAF/DDVtabAR [78] was used with a cycle of 94 \({}^{\circ}\)C denaturation for three minutes, followed by 35 cycles of 94 \({}^{\circ}\)C denaturation for 30 s, 50 \({}^{\circ}\)C annealing for 30 s, and 72 \({}^{\circ}\)C extension for 50 s, followed by 72 \({}^{\circ}\)C for six minutes for the final extension. Nuclear genes were amplified, visualised, and sequenced using the same reagent mixtures as mitochondrial ones. Primary sequence data were assembled from forward and reverse strands in Geneious as above. For both _28S_ and _tuba1_, ambiguous reads at the ends of the assembly could not be resolved and ends were dropped starting from the first ambiguous read. Sequences were aligned in MEGA X as above, and _cob_ and _tuba1_ sequences were translated to assure that sequences corresponded to an amino acid sequence without stop codons. In addition to their phylogenetic use, _cob_ sequences were also used analogously to _cox1_ sequences, estimating the within- and between-species genetic distance for Anisolabididae morphospecies based on the fragment. Gene trees for each fragment were estimated using the maximum likelihood method in RAxML version 8.2.10 [79] with the standard hill-climbing algorithm, a general time-reversible model of site evolution, and a gamma-distributed rate parameter which was estimated empirically. Mitochondrial fragments (_cox1_ and _cob_) were combined into a single mitochondrial data matrix of 1016 bp. Log-likelihoods of estimated trees using rate heterogeneity parameters under gamma and CAT models were compared before selecting the gamma parameter [80]. Node support was estimated with 1000 bootstraps using RAxML's rapid bootstrapping algorithm. Following this, the sequence alignments were concatenated with SequenceMatrix version 1.8 [81]. A partition homogeneity test (incongruence length difference test) was performed in PAUP* version 4 [82] on the concatenated data to determine whether concatenating the fragments was likely to improve phylogenetic signal. The test failed to reject the null hypothesis of congruence at \(\alpha=0.05\). The concatenated dataset contained gaps (e.g., _28S_ could be obtained for one individual and _tub1_ for another, but both fragments could not be obtained for both individuals). Individuals for which only one fragment was obtained were removed prior to the final analysis.