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 <em>Aglaoctenus:</em> <em>A. lagotis, </em>and<em> A. castaneus </em>from Uruguay and Argentina. From those we obtained three mitochondrial markers including the 5’ half of cytochrome c oxidase subunit I (<em>cox1</em>), 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 (<em>12S</em>). Additionally, we sequenced the nuclear intron of the gene encoding translation initiation factor 5A (<em>tif5A</em>). The sequences of the <em>12S</em>, <em>16S</em>+<em>L1+nad1 </em>and <em>tif5A </em>gene fragments were aligned using the online version of MAFFT v7 using the Q-ins-i algorithm. The alignment of the partial fragment of <em>cox1 </em>sequences was trivial since no insertions/deletions (indels) were observed. The alleles in heterozygous individuals for the <em>tif5A</em> 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 <em>tif5A</em>, the concatenated mitochondrial genes <em>cox1</em>+<em>12S</em>+<em>16S</em>+<em>L1</em>+<em>nad1</em>, and a partial fragment of <em>cox1 </em>gene were estimated using statistical parsimony in TCS and implemented in PopART v1.7.