Not the same: phylogenetic relationships and ecological niche comparisons between two different forms of Aglaoctenus lagotis from Argentina and Uruguay

We extracted genomic DNA from three species of Aglaoctenus: A. lagotis, and A. castaneus from Uruguay and Argentina. From those we obtained three mitochondrial markers including the 5’ half of cytochrome c oxidase subunit I (cox1), the 3’ half of the 16S rRNA ribosomal subunit plus the complete tRNA-Leu plus 5’ half of the NADH dehydrogenase subunit I (16S+L1+nad1), and a partial fragment of the small ribosomal unit (12S). Additionally, we sequenced the nuclear intron of the gene encoding translation initiation factor 5A (tif5A). The sequences of the 12S, 16S+L1+nad1 and tif5A gene fragments were aligned using the online version of MAFFT v7 using the Q-ins-i algorithm. The alignment of the partial fragment of cox1 sequences was trivial since no insertions/deletions (indels) were observed. The alleles in heterozygous individuals for the tif5A intron were separated using the PHASE algorithm, as implemented in DnaSP v6. 12.03. With this data we inferred gene trees with Maximum Likelihood, Bayesian and statistical parsimony analyses. We perform molecular species delimitation analyses conducted with STACEY, and the species tree and divergence times were co-estimated with *BEAST. Additionally, we build a haplotype network of the nuclear intron tif5A, the concatenated mitochondrial genes cox1+12S+16S+L1+nad1, and a partial fragment of cox1 gene were estimated using statistical parsimony in TCS and implemented in PopART v1.7.