Molecular phylogeny of the genus Margotrema (Digenea: Allocreadiidae), parasitic flatworms of goodeid freshwater fishes across central Mexico: species boundaries, host-specificity, and geographical congruence

We explored the genetic variation at the intraspecific level of two putative species of Margotrema, M.bravoae and M.guillerminae, to establish the species boundaries. Sequences of a mitochondrial (COI) and nuclear (ITS1) gene were obtained for 125 specimens distributed in 15 localities from seven hydrological systems. An alignment of 750 and 831 bp including gaps of COI and ITS1, respectively, was assembled. We analysed the gene fragments separately as well as together by using Bayesian inference for phylogenetic reconstruction. Based on the phylogenetic analyses, an ultrametric tree was built for each gene, and a general mixed Yule-coalescent model for species delimitation was carried out. A multispecies coalescent analysis was performed in *BEAST using both molecular markers. The results show four independent evolutionary lineages that we interpret as two valid species of Margotrema; the first is represented by an independent lineage, and the second is composed of three lineages. Because the species M.bravoae and M.guillerminae nested together within these three lineages, the validity of M.guillerminae as an independent species is questioned, and it is proposed to represent a junior synonym of M.bravoae. Each lineage shows a congruent geographical pattern with respect to the hydrological system where they occur. Additionally, each lineage shows congruence with respect to the hosts they parasitize, either at the species or at higher taxonomical levels (tribe). The parasite evolutionary history is congruent with the evolutionary and biogeographical history of their hosts.(c) 2013 The Linnean Society of London