Multi-locus phylogeny of the catfish genus Ictalurus Rafinesque, 1820 (Actinopterygii, Siluriformes) and its systematic and evolutionary implications
BMC Ecology and Evolution volume 23, Article number: 27 (2023)
Ictalurus is one of the most representative groups of North American freshwater fishes. Although this group has a well-studied fossil record and has been the subject of several morphological and molecular phylogenetic studies, incomplete taxonomic sampling and insufficient taxonomic studies have produced a rather complex classification, along with intricate patterns of evolutionary history in the genus that are considered unresolved and remain under debate.
Based on four loci and the most comprehensive taxonomic sampling analyzed to date, including currently recognized species, previously synonymized species, undescribed taxa, and poorly studied populations, this study produced a resolved phylogenetic framework that provided plausible species delimitation and an evolutionary time framework for the genus Ictalurus.
Our phylogenetic hypothesis revealed that Ictalurus comprises at least 13 evolutionary units, partially corroborating the current classification and identifying populations that emerge as putative undescribed taxa. The divergence times of the species indicate that the diversification of Ictalurus dates to the early Oligocene, confirming its status as one of the oldest genera within the family Ictaluridae.
Members of the order Siluriformes, commonly known as catfishes, have been the focus of evolutionary studies given their ancient origin during the early Cretaceous , but also as a result of their wide distribution and species diversity. Catfishes have inhabited salt, brackish, and freshwater areas on all continents [2, 3] and are considered a model system for the research of historical continental relationships and biogeographic patterns [3,4,5,6,7]. The evolutionary history of this group has been explored using different information sources [1, 8,9,10]. Their morphological particularities allow an extensive fossil record that provides evidence of their evolutionary and biogeographic history [11,12,13].
Within Siluriformes, the family Ictaluridae is considered to have evolved since the Eocene (ca. 65 million years ago: Mya) and is the only extant catfish group to occur in North America . Ictaluridae currently comprises approximately 51 extant species belonging to seven genera: Ameiurus Rafinesque 1820, Noturus Rafinesque 1818, Pylodictis Rafinesque 1818, Satan Hubbs and Bailey 1947, Trogloglanis Eigenmann 1919, Prietella Carranza 1954, and Ictalurus Rafinesque 1820 . Sixteen additional fossil species (one of the most complete fossil records of North American freshwater fishes) have also been reported . Although molecular phylogenetic analyses have been conducted for species included in some ictalurid genera; i.e., Noturus, Ameiurus, and Prietella (11, 15–16), those of Ictalurus, the widest distributed genus in North America (from Canada to Belize; 11), have been poorly studied to date. The phylogenetic relationships and species composition within Ictalurus remain uncertain, leading to unstable classification. Several species have been synonymized, and individuals of some populations are even considered undescribed taxa [2, 17].
Lundberg  provided the most comprehensive morphological phylogenetic analysis of the genus Ictalurus in a study that resolved two clades: the furcatus group, which includes Ictalurus furcatus Valenciennes, 1840 and Ictalurus balsanus Jordan and Snyder, 1899, and the punctatus group, which includes Ictalurus australis Meek, 1904, Ictalurus dugesii Bean, 1880, Ictalurus lupus Girard, 1858, Ictalurus mexicanus Meek, 1904, Ictalurus pricei Rutter, 1896, and Ictalurus punctatus Rafinesque, 1818. However, the phylogenetic relationships within the punctatus group were unresolved. Although focused on studying higher-rank relationships, the most recent phylogenetic hypothesis, which included species of Ictalurus based on genetic sequences and morphological data, actually revealed Ictalurus balsanus as the earliest diverged species, followed by the divergence between I. furcatus + I. meridionalis and the punctatus group .
Although based solely on morphological evidence, some studies aimed at analyzing the taxonomic cohesiveness of certain widely distributed species suggest the existence of species complexes. Such is the case of Ictalurus lupus (2, 18–19), which occurs along the extensive Bravo River basin and other independent drainages along the Gulf of Mexico slope, such as the Cuatro Cienegas valley and Soto la Marina River. Further examples include I. pricei, with a distribution range comprising several river drainages along the Pacific slope of the Sierra Madre Occidental , and I. dugesii, for which two differentiated western populations have been proposed . On the other hand, there is some controversy regarding species synonymies; Miller et al.,  consider two synonymized species: Ictalurus ochoterenai de Buen, 1946 with I. dugesii and Ictalurus meridionalis with I. furcatus, while Gilbert  considers Ictalurus meeki Meek, 1902, a synonym of I. pricei; and Lundberg  considers I. australis a synonym of I. punctatus, whereas Fricke et al.  consider I. australis, I. ochoterenai, and I. meridionalis as valid species.
Based on the above, and the most complete taxonomic survey of described and putatively undescribed forms of Ictalurus, the present study infers molecular phylogenetic relationships among Ictalurus species, exploring their biogeographic evolutionary history in North American river drainages and resolving certain genus-related taxonomic uncertainties.
We generated 192 novel sequences of cytb, 135 of coxI, 136 of atpase8/6, and 116 of RAG1 (with no indels), for 97 localities across the distribution range of the Ictalurus species in North America (Table S1). In addition, 63 sequences of cytb, 71 of coxI, 2 of atpase8/6, and 3 of RAG1 (including some species of Ictalurus as well as outgroups) were retrieved from GenBank (Table S1). The final amplicon lengths for cytb, coxI, atpase8/6, and RAG1 were 1092, 609, 861, and 1035 bp, respectively. Assessment of substitution saturation based on the entropy-based index indicates that when the index of substitutional saturation (Iss) is smaller than the critical indices of symmetrical (Iss.cSym) and asymmetrical (Iss.cAsym) substitutional saturation, the gene sequences have experienced little substitution saturation, implying that these sequences are useful (23–24). The results of the saturation test in the present study show that the values of Iss were significantly lower than those of Iss.cSym and Iss.cAsym for the three positions at four loci (P < 0.05; Table S2). This indicates that there was little saturation and that the four data sets correspond to useful sequences.
Using maximum likelihood (ML) and Bayesian inference (BI), the phylogenetic relationships with the genes cytb (Fig. 1), coxI (Fig. S1), and RAG1 (Fig. S2) revealed two well-supported clades, corresponding to the furcatus and punctatus groups. For cytb, the furcatus group recovered two clades corresponding to the recognized species I. balsanus and I. furcatus whereas, with coxI, we recovered three clades corresponding to I. balsanus, I. furcatus, and I. meridionalis. As with cytb, for RAG1, I. meridionalis was nested within I. furcatus (Fig. S2).
For coxI and cytb, the punctatus group produced two main clades, one formed by the populations of I. punctatus and the other grouping the rest of the species/populations. In this latter clade, five supported groups were recovered: (i) populations assigned to I. mexicanus and I. australis, (ii) populations of the Nazas basin, (iii) populations of I. lupus, I. aff. lupus from the Soto la Marina River, I. aff. lupus from Cuatro Cienegas and I. aff. lupus from the Conchos River, (iv) populations assigned to I. dugesii, I. dugesii from Armeria, and I. ochoterenai, and (v) populations assigned to I. pricei, I. aff. dugesii from the Santiago River, Ictalurus sp. from Mezquital, and I. aff. pricei from the Culiacan/San Lorenzo Rivers. Within this clade of the punctatus group, the clade I. mexicanus/I. australis was the earliest diverged lineage (Fig. 1; Fig. S1). RAG1 showed no resolution in the relationships among species, indicating I. punctatus as the early-diverged species, which is related to a clade formed by at least three recovered subclades in a polytomy. The mitochondrial atpase8/6 gene did not yield the same interrelationships. Firstly, the furcatus group was not recovered as a monophyletic group, since I. balsanus emerged as the early diverged clade, related to a clade formed by the remaining species of Ictalurus (Fig. S3); and secondly, low resolution was found within the punctatus group (Fig. S3).
Ictalurus aff. dugesii Armeria was resolved as a monophyletic group with cytb and coxI, but only with the first locus showing reciprocal monophyly in relation to I. dugesii samples (Fig. 1; Fig. S1). The remaining populations of I. dugesii are not supported as reciprocal monophyly based on the three mitochondrial loci (Fig. 1, S1, and S3). Ictalurus ochoterenai (a species previously considered as a synonym of I. dugesii) was not recovered as a monophyletic group, since it was also found to be nested within the I. dugesii clade using all genes (Fig. 1; Fig. S1-S3). The differences between cytb and coxI consisted of the topology based on the first locus: the populations of I. aff. lupus from the Conchos River, I. aff. lupus from Cuatro Cienegas valley, and I. aff. lupus from the Soto la Marina River were recovered as reciprocally monophyletic groups in relation to I. lupus, a clade considered as the lupus complex in the present study (Fig. 1). The coxI phylogeny showed the first two populations to be part of I. lupus and related to I. aff. lupus from the Soto la Marina River (Fig. S1). Regarding the relationships based on the cytb, a clade was formed by the differentiated populations of I. aff. dugesii from the Santiago River, Ictalurus sp. from the Mezquital River, I. aff. pricei from the Culiacan/San Lorenzo Rivers, and the recognized species I. pricei. In this study, we refer to this clade as the pricei complex, which showed strong clade support values and was closely related to the lupus complex (Fig. 1). Relationships based on coxI showed the pricei complex as monophyletic, but the relationships among the monophyletic clades included therein were unresolved (Fig. S1). Another difference from the topology obtained with cytb is that the pricei complex was related to I. dugesii in the cox1 tree (Fig. S1).
With the gene RAG1 (Fig. S2), the following phylogenetic mitochondrial/nuclear (mt/nuc) discordance is highlighted: (1) the two specimens from the Guayalejo River within the Panuco basin, morphologically identified as I. australis but assigned as I. punctatus with mitochondrial loci (Fig. 1; S1; and S3), were nested within I. mexicanus samples (Table 1; Fig. S2); (2) two small morphologically undetermined specimens, one from the Santa Maria River in the Panuco basin, and the other from the Pantepec River in the Tuxpan basin, were assigned to I. punctatus using mitochondrial loci (Fig. 1; S1, and S3) but were nested within I. mexicanus/australis, (Table 1; Fig. S2); and (3) another undetermined specimen from the Sabinas River in the Bravo River basin, confirmed as I. punctatus using mitochondrial loci (Fig. 1; S1, and S3), was nested within the lupus complex (Table 1; Fig. S2).
In the concatenated analysis, the specimens I. lupus and I. aff. pricei from the Culiacan/San Lorenzo Rivers were excluded, including the individuals caused by the above phylogenetic discordance. Thus, the concatenated analysis included the recognized species and six previously differentiated populations (I. aff. dugesii Santiago, Ictalurus sp. Mezquital, Ictalurus sp. Nazas, I. aff. lupus Conchos, I. aff. lupus Cuatro Cienegas, and I. aff. lupus Soto la Marina) and their relationships were resolved with high nodal support (Fig. 2). Overall, the concatenated analysis clearly showed the early diversification of the species in the furcatus and punctatus groups. The phylogenetic tree also evidenced the recovery of two species complexes (those of pricei and lupus) and supported the position of Ictalurus dugesii as the sister taxon of the pricei complex, as suggested by analysis of the independent genes coxI, and RAG1 (Fig. S1 and S2).
Uncorrected p genetic divergences ranged from 2.2 to 11.6% for cytb, 2.6 to 10.9% for coxI, 2.2 to 12.7% for atpase8/6, and 0.2 to 1.8% for RAG1 (Tables S2-S5). In particular, the comparisons between different monophyletic groups presented herein exceeded the genetic divergence range within the species level limit observed for ictalurids using cytb (1.8–3.6%) [16, 25,26,27]. Genetic divergences between previously synonymized species and their respective valid nominal species (i.e., Ictalurus ochoterenai – I. dugesii, I. australis – I. mexicanus and I. meridionalis – I. furcatus) were very low, reaching 0.2, 0.5, and 1.5% for cytb; 0.1, 0.3, and 1.3% for coxI; 0.1, 1.8, and 0.5% for atpase8/6; and 0.2, 0.5, and 0.2% for RAG1, respectively. The genetic divergences of known differentiated populations with respect to closely related species included the following: divergences between I. dugesii and I. aff. dugesii Armeria were 0.9, 0.4, 0.6, and 0.2%, for cytb, coxI, atpase8/6, and RAG1, respectively. Divergences between I. dugesii and I. aff. dugesii Santiago were 2.1, 2.6, 1.9, and 0.2%, for cytb, coxI, atpase8/6, and RAG1, respectively. The genetic divergences among the four monophyletic groups recovered within the pricei complex ranged from 1.6 to 1.9% for cytb, 1.6 to 2% for coxI, 1.4 to 2.1% for atpase8/6, and 0.1 to 0.3% for RAG1 (Tables S2-S5), whereas the divergence among the four monophyletic groups recovered within the lupus complex ranged from 1.9 to 3.4% for cytb, 1.1 to 1.8% for coxI, 1.1 to 4% for atpase8/6, and 0.1 to 0.3% for RAG1 (Tables S2-S5).
Nucleotide substitution rates derived from the time-calibrated tree were 3.054 × 10− 3 for cytb, 2.695 × 10− 3 for coxI, 2.877 × 10− 3 for atpase8/6, and 2.897 × 10− 4 for RAG1 substitution/site/per million years. Although presenting low nodal support, the species tree based on the four loci (Fig. 3) showed the same topology as that yielded by the concatenated analysis (Fig. 2). Likewise, for the time-calibrated tree (Fig. S4), the species tree recorded similar divergence times (Fig. 3). The most recent common ancestor (MRCA) of Ictalurus was dated at ca. 27.9 Mya (23.2–33.9 Mya, 95% high posterior density), during the Oligocene. Most of the diversification events within Ictalurus occurred during the Late Miocene-Early Pliocene (Fig. 3).
Bayesian species delimitation
Two Bayesian Phylogenetics and Phylogeography (BPP) analyses were conducted; A10 (comparison of species delimitation models induced on a given “guide” tree ) and A11 (joint comparison of species delimitation/assignment and tree inference of unguided tree ), both of which yielded similar results. In the A10 analyses, almost all the splits involving the recognized species (Ictalurus pricei, I. dugesii, I. mexicanus, and I. punctatus) and differentiated populations (I. aff. pricei Culiacan/San Lorenzo, Ictalurus sp. Nazas, and I. aff. lupus Soto la Marina) were supported by the four scenarios (Table 2). For the splits involving I. balsanus, I. lupus, and the lineages pair I. aff. dugesii Santiago – Ictalurus sp. Mezquital, only the second (IGθ = 3, 0.002, IG τ0 = 3, 0.004) and fourth (IGθ = 3, 0.002, IG τ0 = 3, 0.4) scenarios were supported when ancestral population size was small (Table 2). Finally, the splits of I. furcatus – I. meridionalis, and the differentiated population pair I. aff. lupus Cuatro Cienegas – I. aff. lupus Conchos, were unsupported (Table 2).
In the A11 analyses, almost all recognized species (except Ictalurus furcatus and I. balsanus) and almost all differentiated populations (except I. aff. lupus Cuatro Cienegas and I. aff. lupus Conchos) were supported in at least the same two scenarios as stated above (second (IGθ = 3, 0.002, IG τ0 = 3, 0.004) and fourth (IGθ = 3, 0.002, IG τ0 = 3, 0.4); Table 3).
This is the first study designed to include all described species of Ictalurus (except for I. meeki) within a molecular phylogenetic analysis, considering not only the currently recognized and synonymized species , but also individuals from well-differentiated populations, regarded as undescribed taxa in several published accounts [2, 17], and from several previously unstudied populations. This could allow improved estimation of the species composition of this genus. Our multilocus and species delimitation approach resolved the phylogenetic relationships among Ictalurus species and revealed inconsistencies with respect to the previous taxonomic scheme. Our results also revealed the existence of genetic lineages that could represent undescribed species and corroborated previous synonymies between species proposed on morphological grounds.
Agreement and conflicts between individual genes, concatenated analysis, and species trees: recovering the evolutionary relationships
Differences among genes were mainly found in the relationships of terminal nodes, which can be attributed to the higher content of variable and informative sites for cytb compared to the other mitochondrial genes, coxI and atpase8/6 [28,29,30], as shown by the number of informative characters and substitution rates. Previous studies found that coxI is more conserved than cytb, in which more than half of all the amino acid sites were invariable across 250 fish species . Moreover, for the nuclear gene RAG1, the lack of resolution in the more derived groups, as was the case with the pricei and lupus complexes, could be associated with incomplete lineage sorting from this nuclear gene. As evidenced by previous studies, the mitochondrial locus showed a high mutation rate, which is related to its smaller effective population size  and is also corroborated by the observed substitution rates.
Although the most informative cytb gene placed I. dugesii as the sister group of the pricei and lupus complexes, the concatenated tree, the coalescent species tree, and the coxI, atpase8/6 and RAG1 genes placed I. dugesii as closely related to the pricei complex. Nevertheless, the phylogenetic conflict among trees derived from individual genes was present in the low support shown in the species tree. This latter relationship is consistent biogeographically since the western limit of the distribution range of I. dugesii (the Chapala-Ameca River) is adjacent to the southern distribution range of the pricei complex, which comprises the river basins of the Sierra Madre Occidental, including the Santiago River (Fig. 4). Similar biogeographic patterns have been reported in studies of other co-distributed fish species, including the Goodeids , Cyprinids  and Catostomids . These freshwater fish groups all feature species that inhabit the central Mexican highlands, the sister species of which are distributed in river drainages of the Sierra Madre Occidental. Previous empirical analyses show that discordances between concatenated and coalescent-based analyses tend to occur when the branches of concatenated analyses are short and weakly supported as a consequence of signal conflict between individual gene trees . Based on this, we regard the species tree as a consistent phylogenetic hypothesis with which to depict the relationships and evolution of the genus Ictalurus (Fig. 3).
Evolutionary relationships and their taxonomic implications
According to the BPP analyses, the second scenario for a small ancestral population and low divergences (IGθ = 3, 0.002, IG τ0 = 3, 0.004) and the fourth scenario for a small ancestral population and deep divergences (IGθ = 3, 0.002, IG τ0 = 3, 0.4), in the BPP A10 and A11 analyses, support almost all the splits and assignments in the species delimitation, except for Ictalurus furcatus and I. meridionalis, and the differentiated populations I. aff. lupus Cuatro Cienegas and I. aff. lupus Conchos, which were either supported by only one scenario or not supported in either of the analyses (Tables 3 and 4). The results highlight the particular case of I. balsanus, which was supported by only one scenario in the A11 analysis, as discussed below. Previous studies, in which different prior parameters were tested, and with values that varied considerably, found unstable support patterns; i.e., higher support for one prior combination and lower support for another [37,38,39]. Indeed, BPP analyses can be sensitive and potentially misleading, especially regarding the age of the root (τ0) that, with increasing values (i.e., 0.04 or 0.4), produced spuriously high posterior probabilities possibly because high values of τ push coalescent events closer to the tips . However, this is not the case in the present study. Despite the increase in τ0 (0.004 to 0.4), under the second and fourth scenarios, the supports were unaffected, thus eliminating any misleading results and/or artifacts of the coalescent events. Moreover, there are cases in which θ prior increases considerably (by at least two orders of magnitude), as occurred in the present case (i.e., 0.002–0.2) while speciation probabilities tend to decrease . This is seemingly associated with the use of an extremely high θ, which is not suitable in all cases . Thus, the high values for θ used in scenarios 1 and 3 (0.2) resulted in an unsuitable set of priors for Ictalurus. This indicates that scenarios 2 and 4, which supported 13 of the 15 proposed groups as potentially well-delimited taxa, are more suitable for discarding a misleading result.
As with the earlier phylogenetic hypothesis based on morphology , a dichotomy between the furcatus and punctatus groups was found in our analyses. Interestingly, our study revealed several differences in species relationships when compared with the most recent phylogenetic hypothesis involving catfishes of the family Ictaluridae, in which the molecular data and morphological characters of extant and fossil species were used .
The furcatus group
The results of the species delimitation analyses do not support the split of Ictalurus meridionalis and I. furcatus in the A10 analysis, and the assignments in the A11 analysis for both of these taxa were highly supported only by the second BPP scenario and, in the case of I. meridionalis, only moderately supported by the third scenario. It is therefore not possible to validate these two taxa and, based on the principle of priority, I. meridionalis is considered a synonym of I. furcatus, corroborating its status as a synonymized species as previously stated by Lundberg . This is consistent with the fact that only coxI supported the reciprocal monophyly (Fig. S1) and that values of genetic divergence of 1.5% were found with cytb between I. meridionalis (Hondo River and Peten samples) and I. furcatus from across a wide geographic range (samples from the Mississippi to Coatzacoalcos). This divergence value lies below the minimum inter-species level within Ictaluridae (1.8-3.6%) [15, 25, 26, 41], and these results indicate that the morphological differentiation reported  could be the result of geographical variation within the widespread I. furcatus. Although the geographic coverage for I. furcatus is incomplete, the samples collected from the extremes of the distribution range of the species (from the Mississippi River in the USA to the Hondo River in Belize) presented two structured populations. Therefore, a future phylogeographic study including optimal sampling along the distribution range of the species is necessary to determine the genetic pattern of these two divergent groups.
Unexpectedly, the assignment of the species Ictalurus balsanus, which is one of those that diverged early within the genus, was only supported by the second scenario in the A11 analysis. However, the deep divergence level of this taxon strongly validates it as a recognized species. Based on the results presented herein, we consider the furcatus group to be formed by the species I. balsanus and I. furcatus, as was originally established , pending a phylogeographic study with a wider geographic range for I. furcatus.
Our results do not entirely support the previous phylogenetic hypothesis , which included molecular data and morphological characters as a total evidence analysis. In this study , I. balsanus was found to be the early divergent species and the sister species of the clade containing all other living and fossil species within the genus. This discrepancy might have been the result of one of the following: the possible influence of the different phylogenetic signals from the genes used (16s, 12s, and RAG2, which were not used in the present study) or the methodological implications associated with the considerable amount of missing data, particularly the lack of complete molecular data for most taxa in the analysis by Arce-H et al. . For instance, only two of the five genes were included in that analysis for I. balsanus. Moreover, it is known that a large quantity of missing data, together with high rates of change, can lead to long-branch attraction . This occurs more frequently using Maximum Parsimony than probabilistic methods  such as those used in Arce-H et al. .
The punctatus group
The punctatus group included the five remaining recognized species of Ictalurus (I. pricei, I. lupus, I. dugesii, I. mexicanus, and I. punctatus) and five well-differentiated populations (Ictalurus sp. Santiago, Ictalurus sp. Mezquital, Ictalurus sp. Culiacan/San Lorenzo, Ictalurus sp. Conchos/Cuatro Cienegas, and Ictalurus sp. Nazas), four of which were previously recognized as undescribed species [2, 17, 19], and the other is a new well-differentiated population revealed in the present study (Ictalurus sp. Soto la Marina). All 11 of these groups were supported by BPP analyses (Tables 3 and 4).
Ictalurus punctatus: The channel catfish presents a native distribution stretching from the north in the Great Lakes, southern Canada, Hudson Bay (Red River drainage), and the Missouri-Mississippi River basins, extending southwards to the Panuco River on the Mexican Atlantic slope [17, 19, 45]. Although the present study lacks broad geographic coverage for I. punctatus, we included samples from both nearby (e.g., individuals from the same locality, Lake Michigan or the Grande River) and distant (e.g., individuals from Canada and the Mississippi to the Panuco River in Mexico) geographical points. This allowed us to cover the genetic variation and differentiation that exists within this species, and to confirm its phylogenetic position within the genus. All the samples identified as I. punctatus (excluding the specimens involved in the mt/nuc discordance) were recovered as a monophyletic assemblage in the concatenated and species tree analyses. In addition, low genetic differentiation between nearby or distant geographic samples was detected with all genes, indicating that all these populations should be treated as a single genetic group or species. Interestingly, a similar pattern was uncovered when dealing with a specific parasite of channel catfishes . The trematode Phyllodistomum lacustri Loewen, 1929, which occurs in the urinary bladder of its hosts, was genetically identical in channel catfishes sampled in Canada, the USA, and northern Mexico . However, a cryptic species complex of P. lacustri was uncovered in populations of catfishes occurring in several river basins of Mexico, where Ictalurus has experienced a great diversification process .
All the specimens of channel catfish sampled were morphologically identified or genotyped (with mitochondrial genes) as I. punctatus; however, we found four cases of mt/nuc discordance in the phylogenetic position, that involved specimens from the Santa Maria River in the Panuco River basin, and the Pantepec River in the Tuxpan River basin. These were assigned to the I. mexicanus clade with the nuclear locus, and a specimen from the Sabinas River within the Bravo River basin was assigned as part of the lupus complex with the nuclear locus.
Due to the high genetic divergence values between I. punctatus and the I. mexicanus – lupus complex (7.5% based on cytb in the former two, and 7.8–8.7% in the latter), as well as the phylogenetic position with both mtDNA and nDNA of most samples from the I. mexicanus – lupus complex, the possibility of an incomplete lineage sorting can be discarded. Another plausible explanation could be the occurrence of a hybridization event. However, a morphological revision of adult individuals would be necessary , as well as an extensive study with nuclear loci, to obtain, for example, heterozygous individuals with nuclear loci, indicating the presence of one allele of each taxon, as has been found for other native species; e.g., I. pricei and I. lupus vs. I. punctatus by Gutierrez-Barragán et al. .
Ictalurus australis and I. mexicanus: Both species were described by Meek , the type locality of I. australis is the Forlon River, a tributary of the Guayalejo River sub-basin of the Panuco River Basin, and its distribution extends along the Panuco, Tuxpan, Cazones, Tecolutla and Nautla Rivers , reaching the Blanco River within the Papaloapan River basin . The type locality of I. mexicanus is the Gallinas River, and it is considered endemic to this tributary in the upper Panuco (17, 52–53). The taxonomic status of I. australis has been questioned and the species has even been considered a junior synonym of I. punctatus (2, 54–55), although several studies still consider I. australis as a valid taxon (19, 56–57). However, according to Miller et al. , I. mexicanus is restricted to the Gallinas River and the populations in the Panuco basin are considered a closely related but undetermined taxon.
We included specimens from the type locality for both species and, according to the low genetic divergence shown between them (< 0.5% with the mitochondrial cytb gene) and the lack of reciprocal monophyly within Ictalurus mexicanus/australis with the cytb and coxI genes, we consider that the samples of I. mexicanus from Gallinas River up- and downstream of the Tamul waterfall, those distributed along the Verde River, and those from the Guayalejo River (according to RAG1 gene) do not represent two distinct lineages, suggesting that I. mexicanus, I. australis and/or Ictalurus sp. sensu Miller et al.  are part of the same taxonomic unit. Thus, based on the principle of priority, and pending the inclusion of more samples from across the distribution range of this taxon and integrative taxonomic analyses, we propose recognition of this population as I. mexicanus. The lack of I. australis samples from the southern basins such as the Tuxpan, Cazones, Tecolutla, Papaloapan, and Nautla Rivers (50–51) indicates the need to include samples from the entire potential distribution range in a phylogeographic analysis to explore the relationship with, and patterns of differentiation from, I. mexicanus of the Panuco basin.
The pricei complex: northwestern pacific group: This monophyletic group is represented by one recognized species, Ictalurus pricei, and three well-differentiated populations with a divergence of 1.8–1.9% at the cytb locus. This value exceeds the minimum divergence found between well-recognized species within Ictaluridae (15, 29–30). These four divergent groups were corroborated by the BPP species delimitation test as independent taxonomic units. The northernmost taxon corresponds to I. pricei sensu stricto from the Yaqui and Fuerte rivers. We did not include samples from other basins such as the Sonora, Mayo, and Casas Grandes Rivers [17, 58], and these populations therefore remain to be analyzed to determine whether or not they can be nested with I. pricei from the Yaqui and Fuerte rivers.
Another differentiated group was represented by the specimens of Ictalurus sp. from the Culiacan/San Lorenzo rivers, confirming previous morphological findings that indicated this population as a potential undescribed taxon . Similarly, the specimens of Ictalurus sp. from the Mezquital River were corroborated as representing a potentially undescribed taxon in our analyses, as previously recognized by other authors (2, 18–19). Finally, the differentiated population of Ictalurus sp. from the Santiago River, which corresponds to the southernmost independent taxonomic unit of the pricei complex, was also corroborated as a potentially undescribed taxon, originally referred to by Miller et al.  as an undescribed species related to I. dugesii. The results of a study by Rosas-Valdez et al.  on the recognition of a cryptic species complex of the trematode Phyllodistomum lacustri may lend further support. The populations of the parasite sampled in catfishes of the Mezquital River basin (identified by the authors as Ictalurus sp.), and those from the Lerma River basin at San Juanico Dam (identified as I. dugesii), but also including introduced specimens of I. punctatus, were retrieved as cryptic species in the molecular phylogenetic analyses of two genes (28s and coxI), with each representing a separate species.
Ictalurus dugesii: The Ictalurus samples from the Lerma-Chapala system, including specimens identified as I. ochoterenai from Chapala Lake and the populations from the Ameca and Armeria River basins, were nested within the I. dugesii clade. Regarding I. ochoterenai and I. dugesii, we found a lack of reciprocal monophyly, as well as a genetic divergence (< 0.5%) below the minimum limit recorded for Ictaluridae using cytb (1.8–3.6%) [15, 25, 26, 41]. Based on the above, the present study confirms I. ochoterenai as a synonym of I. dugesii, as proposed by other authors [2, 59].
An undescribed form of I. dugesii was proposed by Miller et al. , corresponding to the populations of I. dugesii from the Armeria River. In this case, even though we found reciprocal monophyly with cytb, this group was nested within I. dugesii when data were analyzed using the other two mitochondrial genes. It also presented low genetic divergence (0.2%), which suggests a relatively recent isolation event within the Armeria River basin. Thus, the distribution of I. dugesii is confirmed for three disjunct river basins in central Mexico, i.e., the Lerma-Chapala system and the Ameca and Armeria Rivers. This is in accordance with other co-distributed fish species in the region, such as the catostomid Moxostoma austrinum Bean, 1880 , and several species of goodeids (61–62).
The lupus complex: The distribution range of the headwater catfish, I. lupus comprises several rivers of the Gulf of Mexico slope, including the Colorado, Guadalupe, and Nueces drainages in Texas, the Bravo River drainage in the United States and Mexico, as well as the San Fernando, Soto la Marina, and Conchos Rivers, and the endorheic Cuatro Cienegas basin in Coahuila, Mexico . For the cytb analyses, I. lupus sensu stricto from the Pecos and Devil Rivers presented a relatively high genetic divergence with respect to the other populations (Ictalurus sp. Conchos/Cuatro Cienegas, and Ictalurus sp. Soto la Marina), ranging from 2 to 2.8%. Although Ictalurus sp. Conchos and Ictalurus sp. Cuatro Cienegas presented a high divergence value (1.9% in the cytb gene), and the species delimitation test resolved only three groups as potentially independent taxonomic units: I. lupus, Ictalurus sp. Conchos/Cuatro Cienegas, and Ictalurus sp. Soto la Marina (Tables 3 and 4).
The genetically differentiated Ictalurus sp. Conchos seems to correspond to the previously reported divergent population of Ictalurus from the Conchos River, which also occurs in sympatry with I. lupus in the lower Bravo River (2, 17–18, 54). In addition, Ictalurus sp. Soto la Marina was supported as a candidate species in the BPP analyses and presented the highest genetic divergence within the lupus complex, merging as a novel recognized independent evolutionary lineage and putative species. According to our findings, an integrative taxonomic study with broader geographic coverage of I. lupus, as well as a detailed study of the morphoanatomy of this species, are required to test the possible taxonomic independence of the three divergent evolutionary units within the lupus complex: I. lupus sensu stricto, Ictalurus sp. from the Conchos River and Cuatro Cienegas valley, and Ictalurus sp. from the Soto la Marina River.
Ictalurus sp. Nazas: Our study resolved the specimens of Ictalurus from the Nazas River as a highly divergent group within the punctatus complex. All the analyses confirmed this population as a reciprocally monophyletic group, which was apparently the first to diverge within the central-western drainages, indicating this as a potential taxon in the BPP analyses. Our results support previous studies that consider the population of Ictalurus from the Nazas as an undescribed species, although these studies provide no details regarding the morphological characters that distinguish this from other congeners [2, 18, 19, 63] and, consequently, a description of the new species is still pending. Other studies indicate that the Nazas River is an area where independent genetic lineages of fishes occur [64,65,66,67,68].
Divergence times and evolution of Ictalurus
The genus Ictalurus is confirmed as an ancient group of fishes; our dating analysis placed the earliest split between the punctatus and furcatus groups at ca. 33 − 23 Mya (Fig. 3), which is partially consistent with previous studies that estimate the age of the crown group of Ictalurus at 38 − 20 Mya  and 37 − 30 Mya . The considerable species richness in the fossil record, together with the inferred divergence times, indicate that the genus has experienced several extinction events during its evolutionary history, particularly in the case of early-diverging species. The most abundant fossil record for Ictalurus has been found in North America (13, 70–71), and includes the oldest fossil, Ictalurus rhaeas Cope, 1891, discovered at the Cypress Hill Oligocenean Formation, located in southern Saskatchewan, Canada. Hydrographically, this formation is included in the Missouri-Mississippi drainage system . The region originated in the Late Eocene and is one of the most important in terms of explaining the evolution of the North American freshwater fish fauna, including ancient fish fauna elements such as gars (Lepisosteidae), bowfin (Amiidae), and mooneyes (Hiodontidae), as well as the most recent representatives of the North American fish fauna, such as catostomids, salmonids, percids, cyprinids, and catfishes . The first lineages, derived during the evolution of the genus Ictalurus in the early Miocene, show high dispersal capacity and an adaptation to the environmental conditions that allowed them to expand their distribution range along the Atlantic slope of North America, as reflected in the widespread distribution of I. furcatus and I. punctatus, from southern Canada through the US and Mexico southwards to Belize, including the dispersal of the ancestor of I. balsanus to the Balsas River on the Pacific slope of Mexico.
The considerable diversification observed in the punctatus group shows a first split with the cladogenetic event of I. mexicanus in certain tributaries of the Panuco River basin. Although no explanations have been proposed for the origin and isolation of this species in the Verde or Gallinas rivers (currently inhabited by I. mexicanus), this region has served as an area of speciation for different groups, with the occurrence of a significant number of endemic species of cyprinids , goodeids , cichlids (73–74) and poecilids . This high level of endemism seems to be related to the intense volcanism of the Miocene-Pliocene, which promoted sudden subsidence of the graben structure in the basin and the formation of shallow lakes . Ictalurus mexicanus followed the same speciation pattern as other fish species in the area.
The next diversification event was the cladogenesis of Ictalurus sp. Nazas, and the lupus complex distributed in the Northern basins of the Bravo and Soto la Marina Rivers, followed by the formation of the I. dugesii and the pricei complex with wide distribution along the Pacific slope drainages, from the Yaqui River in Sonora southwards to the Lerma River in Michoacan and Guanajuato. Several fish groups experienced the same speciation pattern, including the cyprinids and catostomids [35, 65]. These speciation events associated with the Bravo, Nazas, and Northern pacific slope drainages have been explained by hypothetical connections along the extensive paleo-hydrological system in the Chihuahuan desert region, which dates back to the Oligocene and included the Conchos River, the main channel of the Bravo River, the Nazas River basin, and the headwaters of several western Pacific river drainages (i.e., the Yaqui and Mezquital Rivers) , and the isolation of which is estimated at ca. 5 Mya . The biogeographic role of the Chihuahuan desert paleo-hydrological system is closely associated with the tectonic activity of the Bravo River rift, together with the arid conditions that have been prevalent since the Miocene . This paleo-hydrological system could have had two main roles in the evolution of the genus Ictalurus, in a similar manner to that reported for other fish groups (60, 64, 66–67, 76,77,78): (1) most of the diversification events occurred within or were promoted across the region, and (2) the region acted as a corridor for the punctatus group, allowing them to colonize the Pacific slope, where they subsequently diversified.
In general, derived from a temporary and geographically extensive evolutionary history, most of the speciation events within the genus Ictalurus occurred in the river basins of Mexico. This implies an increase in both species richness and levels of endemism. The results of our study pose a huge challenge for research on the diversity and conservation of this representative and important group of North American freshwater fishes. The most recent evaluation of the IUCN red list for fish species in Mexico categorizes only I. pricei and I. mexicanus in the categories “Endangered” and “Vulnerable”, respectively , mainly based on the premise of the high distribution range of most of the species. However, our research contradicts this notion of the wide distribution of the species and raises the need for a more detailed conservation evaluation of the independent evolutionary units or undescribed species found in the present study to adequately protect the true diversity of the genus Ictalurus in Mexico.
Finally, parasitological data is frequently used, along with the phylogenetic history of hosts, to explain concurrent patterns. In this context, the members of the punctatus group were typical hosts of the trematode Phyllodistomum lacustri, a parasite of the urinary bladder of several species of ictalurids . Molecular phylogenetic analyses of specimens of P. lacustri from Ictalurus punctatus sampled from its natural distribution range across Canada, the USA, and northern Mexico show a pattern congruent with that of the hosts . Interestingly, the presence of a cryptic species complex of P. lacustri that includes other species of Ictalurus, such as I. dugesii in the Santiago basin, and I. pricei (Ictalurus sp. Mezquital in this study) from the Mezquital basin  was also reported. These river basins of western Mexico appear to be important areas for the diversification of I. dugesii and I. pricei and their parasites. We believe that future studies aimed at uncovering diversification patterns of freshwater fishes should, whenever possible, utilize other information sources to corroborate hypotheses. We also consider that the history of the host-parasite association could constitute a robust proxy with which to test such patterns of diversification.
The results obtained in the present study of the systematics, biogeography, and evolution of the genus Ictalurus are based on the most extensive and comprehensive taxonomic sampling to date. Despite the conflicting phylogenetic signals from individual genes, the results of this study conducted using several approaches, including gene trees, a concatenated dataset of nuclear and mitochondrial genes through a multi-locus approach, and the coalescent species tree, recovered a complementary signal, and thus yielded a plausible phylogenetic hypothesis for the genus. Moreover, our study supports the recognition of 13 taxonomically independent units, including seven previously recognized species, I. punctatus, I. mexicanus, I. dugesii, I. lupus, I. pricei, I. furcatus, and I. balsanus, and six putative undescribed species, Ictalurus sp. Culiacan/San Lorenzo, Ictalurus sp. Mezquital, Ictalurus sp. Santiago, Ictalurus sp. Conchos/Cuatro Cienegas, Ictalurus sp. Soto la Marina, and Ictalurus sp. Nazas. These independent evolutionary units require a detailed morphoanatomical study to be formally described as new ictalurid species. We also corroborated previously synonymized species such as I. ochoterenai, I. australis, and I. meridionalis. As in other taxonomically underestimated fish groups (35, 62, 65, 78–79), there has been a lack of morphological studies and/or taxonomic revisions addressing these differentiation patterns in Ictaluridae, a conspicuous, economically important, popular, and important group of fishes. An integrative taxonomy approach is required to describe the evolutionary history of the group and to achieve a better understanding of its species composition and thereby improve its conservation.
Materials and methods
Sampling and sequencing
A total of 187 individuals of the genus Ictalurus were sampled, including members of all recognized species, according to Lundberg  and Miller et al. . Sampling also included several populations recognized as undescribed taxa (2, 17–18, 20), other poorly studied populations, and three previously synonymized species  (Fig. 4; Table 5 and S1). Most of the samples were obtained from three tissue banks: Colección de Peces de la Universidad Michoacana de San Nicolás de Hidalgo, Mexico (CPUM_SEMARNAT-Mich-PEC-227-07-09), Colección Ictiológica del Museo Nacional de Ciencias Naturales, España (MNCN_ICTIO), and Colección de Tejidos del Laboratorio de Genética de la Conservación, CIBNOR, Mexico (LGC_ADN) (Table S1). Other tissue samples were donated by fishermen.
Isolation of genomic DNA was performed using the QIAGEN BioSprint Dneasy Tissue and Blood Kit (Qiagen, Valencia, Ca, USA), following the manufacturer’s instructions. The mitochondrial genes Cytochrome b (cytb), Cytochrome oxidase subunit 1 (coxI), and ATP synthase 6 and 8 (atpase8/6) were amplified using the primers GludG  and H16460 , FISHF1 and FISHR1 , COIII.2, and ATP8.2 , respectively. The nuclear recombination-activation 1 gene (RAG1) was amplified using the primers RAG1F and RAG9R . Amplifications were conducted in a final reaction volume of 25 µl, comprising: 50–100 ng genomic DNA, 1x PCR buffer (containing 1.5mM of MgCl2), 0.2 mM of each dNTP, 0.5 µM of each primer and 1 U of Taq DNA polymerase (Invitrogen). The thermocycler parameters for all amplifications consisted of an initial 2 min denaturation step at 94 ºC. Subsequent cycling parameters for each region were as follows: For cytb, 35 cycles of 45 s at 94 ºC, 60 s at 48 ºC and 90 s at 72 ºC, with a final extension step of 5 min at 72 ºC; for coxI, 35 cycles of 30 s at 94 ºC, 30 s at 52 ºC and 60 s at 72 ºC, with a final extension step of 10 min at 72 ºC; for atpase8/6, 35 cycles of 45 s at 94 ºC, 60 s at 48 ºC and 90 s at 72 ºC, with a final extension step of 5 min at 72 ºC; and for RAG1, 35 cycles of 30 min at 95 ºC, 45 s at 56 ºC and 90 s at 72 ºC, with a final extension step of 7 min at 72 ºC. Positive amplicons were purified using the ExoSAP-IT PCR Product Cleanup Reagent. Sequencing was performed by the Macrogen Korea sequencing service. Sequences were manually aligned in Mega v7.0 . RAG1 alleles were separated using the PHASE algorithm , as implemented in DnaSP v. 6.10.01 . The datasets generated during the current study are available in Genbank (accession numbers ON008558-ON023996). Taxon names and voucher information and associated Genbank numbers are presented in an additional file in Table S1. Moreover, additional sequence data for cytb, coxI, atpase8/6, and RAG1 of the ingroup and the functional outgroups Cranoglanis bouderius Richardson, 1846 and Pangasianodon hypophthalmus Sauvage, 1878 were retrieved from Genbank (Table S1).
Phylogenetic analyses and sequence divergence
We obtained the best-fitting models of nucleotide substitution for each locus (Table 4) considering the Akaike information criteria (AIC) in the program jModelTest2 (88–89). Substitution saturation for the three codon positions at all four loci was assessed with Xia´s method (23–24) using DAMBE 7.3.32 , which estimates an entropy-based index of substitution. To infer the phylogenetic relationships, we conducted a Bayesian Inference (BI) using MrBayes 3.2.3 . We used four simultaneous and independent runs consisting of four chains for 10 million generations, with trees sampled every 1000 generations to calculate posterior probabilities (PP). Convergence and stationary distribution of the runs were verified by the average standard deviation of the split (ASDSF < 0.01) and with a potential scale reduction factor (PSRF) close to 1.0 for all parameters, as indicated in MrBayes. The effective sample size (ESS) for all the parameters was > 200, as visualized in Tracer 1.5 . The first 10% of generations were removed as burn-in. We also performed a Maximum Likelihood (ML) analysis on RaXML v.7.2.8, implementing a GTRCAT model with 1000 bootstrap replicates [93,94,95]. The findings for the best-fitting model of nucleotide substitution and the two phylogenetic reconstruction approaches were run in the CIPRES Science Gateway V.3.3 portal . Both probabilistic methods were run for all four individual loci, cytb, coxI, atpase8/6, and RAG1, as well as for a concatenated matrix including all four loci. Two important points should be noted: firstly, due to the lack of at least two loci for I. lupus and I. aff. pricei Culiacan/San Lorenzo (Table S1), these genetically differentiated groups were excluded from the concatenated matrix, and secondly, based on the findings of reticulate relationships (see results), several individuals morphologically determined as I. australis and I. punctatus (Table S1) were also excluded from the concatenated analyses. Sequence divergences among all recognized, synonymized and undescribed species (Table 5) were estimated using the mean among-groups uncorrected (p) genetic divergences (Table S3-S6) in MEGA 7.
Coalescent-based species tree and divergence time estimation
Using the four genes concatenated matrix (Table 4), we conducted a multispecies coalescent analysis for 95 individuals, representing the six recognized species (except I. lupus), one synonymized species (I. meridionalis), and five differentiated populations (I. aff. lupus Conchos, I. aff. lupus Cuatro Cienegas, I. aff. dugesii Santiago, Ictalurus sp. Mezquital, and Ictalurus sp. Nazas), all recovered as monophyletic groups in the concatenated analysis and presenting a higher cytb genetic divergence than that found between well-recognized species of ictalurids (1.8 − 3.6%) [15, 25, 26, 41]. The synonymized I. meridionalis, which presented a genetic divergence of Dp = 1.5% (a value close to the minimum limit), was also included as an independent taxon in the analysis. Moreover, I. ochoterenai and I. australis, with Dp < 0.6, which is lower than that of their closest relatives I. dugesii and I. mexicanus, were excluded. For this analysis, we used Ameiurus natalis Lesueur, 1891 as an outgroup in *BEAST  implemented in the software Beast v.1.8.3 . To calibrate the species tree, we performed a time-calibrated tree with four fossil constraints, which allowed us to obtain the substitution mean rates for all four loci for subsequent use in the species tree, applying an uncorrelated lognormal relaxed clock prior (see Appendix S2). For species tree priors and population size models, a speciation Yule process and a constant species tree population size were chosen. We performed four independent MCMC runs for 300 million generations, sampling every 3000 generations. Chains convergence was assessed by visualizing the sampled parameter values in Tracer 1.5. We discarded 10% of the generations as burn-in and pooled the estimated parameters using the Log-Combiner module in the BEAST package. The maximum clade credibility species tree was obtained from the Tree-Annotator module in the BEAST package.
Species delimitation analysis
A Bayesian multi-locus species delimitation analysis was conducted using the program Bayesian Phylogenetics and Phylogeography (BPP) v3.4 . In the first approach, an A10 analysis (species delimitation test using a user-specified guide tree [100,101,102] was performed. This method uses a multispecies coalescent model to compare the posterior probabilities of different species delimitation models (100–101) and accommodates lineage sorting due to ancestral polymorphisms . Given that it was possible to perform the analysis with unequal taxon sampling for some loci, the species I. lupus and I. aff. pricei Culiacan/San Lorenzo, which lack the atpase8/6 and RAG1, were included. The topology based on cytb was used as a user-specified guide tree with which to conduct the species delimitation, enabling us to assess all seven recognized species (Table 5), one synonymized species (I. meridionalis), and six differentiated populations (I. aff. lupus Conchos, I. aff. lupus Cuatro Cienegas, I. aff. dugesii Santiago, Ictalurus sp. Mezquital, Ictalurus sp. Nazas, and I. aff. pricei Culiacan/San Lorenzo). Moreover, to validate the outcomes of the A10 analysis of species delimitation, we also performed the A11 analysis (a joint comparison of species delimitation/assignment and species tree estimation) [28, 102] using the same set of priors. Following the strategy of Chan and Grismer  to assess the species delimitation/assignment using A11, the present study conducted this analysis independently in four main groups: the pricei complex, the lupus complex, the punctatus group, and the furcatus group.
The following parameter sets were specified for both the A10 and A11 analyses: a fixed guide tree, a Dirichlet distribution (α = 2) to account for variation in mutation rates among loci and to discern how the effective ancestral population size and the divergence influenced the results, and an inverse-gamma prior (IG) used to specify the population size parameter θ and root-age τ0 of the species tree considering four different scenarios. These used α = 3 for a diffuse prior and adjusted the β to cover the different alternative scenarios: (1) large ancestral population sizes and deep divergences (IGθ = 3, 0.2, IG τ0 = 3, 0.4); (2) small ancestral population sizes and shallow divergences (IGθ = 3, 0.002, IG τ0 = 3, 0.004); (3) large ancestral population sizes and shallow divergences (IGθ = 3, 0.2, IG τ0 = 3, 0.004), and (4) small ancestral population sizes and deep divergences (IGθ = 3, 0.002, IG τ0 = 3, 3, 0.4). Reversible jump (rj) MCMC was run for 500,000 generations, with a burn-in of 8000 and a sampling frequency of two. To assess convergence, we performed the analyses three times to confirm consistency between runs. Posterior probabilities (pp) ≥ 0.95 were considered highly supported, pp ≥ 0.90 and < 0.95 were considered moderately supported, and pp < 0.90 were considered weakly supported.
The datasets used and/or analyzed during the current study are available in the GenBank-NCBI database under the accession numbers ON008558-ON009000, ON009002 and ON023863- ON023993, ON023996.
Kappas I, Vittas S, Pantzartzi CN, Drosopoulou E, Scouras ZG. A time-calibrated mitogenome phylogeny of catfish (Teleostei: Siluriformes). PLoS ONE. 2016;11:12. https://doi.org/10.1371/journal.pone.0166988
Lundberg JG. The phylogeny of Ictalurid catfishes: a synthesis of recent work. In: Mayden RL, editor. Systematics, historical Ecology, and North American Freshwater Fishes. Stanford: Stanford University Press; 1992. pp. 392–420.
Nelson JS, Grande TS, Wilson MVH. Fishes of the world. 5th ed. Hoboken: Jonh Wiley & Sons, Inc; 2016.
Fink SV, Fink WL. Interrelationships of the ostariophysan fishes (Pisces, Teleostei). Zool J Linnean Soc. 1981;72:297–353. https://doi.org/10.1111/j.1096-3642.1981.tb01575.x
Fink SV, Fink WL. Interrelationships of ostariophysan fishes (Teleostei). In: Stiassny MLJ, Parenti LR, Johnson GD, editors. Interrelationships of fishes. New York: Academic Press; 1996. pp. 209–47.
Saitoh K, Miya M, Inoue JG, Ishiguro NB, Nishida M. Mitochondrial genomics of Ostariophysan fishes: perspectives on phylogeny and biogeography. J Mol Evol. 2003;56:464–72. https://doi.org/10.1007/s00239-002-2417-y
Gayet M, Meunier FJ. Palaeontology and palaeobiogeography of catfishes. In: Arratia G, Kapoor BG, Chardon M, Diogo R, editors. Catfishes. Enfield: Science Publishers, Inc; 2003. pp. 491–522.
Diogo R. Higher-level phylogeny of Siluriformes: an overview. In: Arratia G, Kapoor BG, Chardon M, Diogo R, editors. Catfishes. Enfield: Science Publishers, Inc; 2003. pp. 353–84.
Diogo R. Phylogeny, origin and biogeography of catfishes: support for a Pangean origin of modern teleosts’ and reexamination of some mesozoic pangean connections between the Gondwanan and Laurasian supercontinents. Anim Biol. 2004;54:331–51. https://doi.org/10.1163/1570756042729546
Sullivan JP, Lundberg JG, Hardman M. A phylogenetic analysis of the major groups of catfishes (Teleostei: Siluriformes) using RAG1 and rag2 nuclear gene sequences. Mol Phylogenet Evol. 2006;41:636–62. https://doi.org/10.1016/j.ympev.2006.05.044
Grande L, de Pinna MCC. Description of a second species of the catfish Hypsidoris and a reevaluation of the genus and the family Hypsidoridae. J Vertebr Paleontol. 1998;18:451–74. https://doi.org/10.1080/02724634.1998.10011074
Ferraris CJ. Checklist of catfishes, recent and fossil (Osteichthyes: Siluriformes), and catalogue of siluriform primary types. Zootaxa. 2007;1418:1–628. https://doi.org/10.11646/zootaxa.1418.1.1
Arce -HM, Lundberg JG, O’Learya MA. Phylogeny of the north american catfish family Ictaluridae Teleostei: Siluriformes) combining morphology, genes and fossils. Cladistics. 2017;33:406–28. https://doi.org/10.1111/cla.12175
Froese R, Pauly D. FishBase 2000: concepts, design and data sources. In: International Center for Living Aquatic Resources Management (ICLARM) 2000. https://digitalarchive.worldfishcenter.org/bitstream/handle/20.500.12348/2428/WF_311.pdf?sequence1=. Accesed 15 march 2022.
Hardman M. The phylogenetic relationships among bullhead catfishes of the genus Ameirus (Siluriformes: Ictaluridae). Copeia. 2003; 2003: 395–408. https://doi.org/10.1643/0045-8511(2003)003[0020:PRABCO]2.0.CO;2.
Wilcox TP, García de León FJ, Hendrickson DA, Hillis DM. Convergence among cave catfishes: long-branch attraction and a bayesian relative rates test. Mol Phylogenet Evol. 2004;31:1101–13. https://doi.org/10.1016/j.ympev.2003.11.006
Miller RR, Minckley WL, Norris SM. Freshwater fishes of México. 1st ed. Chicago: University of Chicago Press; 2005.
Smith ML, Miller RR. The evolution of the río Grande basin as inferred from its fish fauna. In: Hocutt CH, Wiley EO, editors. The zoogeography of north american freshwater fishes. New York: John Wiley and Sons; 1986. pp. 457–85.
Espinosa PH, Gaspar MTD, Fuentes MP. Listados faunísticos de México III: Los peces dulceacuícolas mexicanos. 1st ed. México D.F.: Universidad Nacional Autónoma de México; 1993.
Ruiz-Campos G, Varela-Romero A, Ceseña-Gallegos D, Ballesteros-Córdova CA, Sánchez-Gonzáles S. Morphometry and meristics of two native catfishes from the Sierra Madre Occidental, México: I. pricei and Ictalurus sp. (Siluriformes: Ictaluridae). Rev Biol Trop. 2020;68:479–91. https://doi.org/10.15517/RBT.V68I2.37041
Gilbert CR. Type catalogue of recent and fossil north american freshwater fishes: families Cyprinidae, Catostomidae, Ictaluridae, Centrarchidae and Elassomatidae. Florida Museum of natural history. Gainesville: University of Florida; 1998.
Fricke R, Eschmeyer WN, Van der Laan R, editors. 2022. eschmeyer’s catalog of fishes: genera, species, references. http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp. Accessed 05 dec.
Xia X, Xie Z, Salemi M, Chen L, Wang Y. An index of substitution saturation and its application. Mol Phylogenet Evol. 2003;26:1–7.
Xia X, Lemey P. Assessing substitution saturation with DAMBE. In: Lemey P, editor. The phylogenetic handbook: a practical Approach to phylogenetic analysis and hypothesis testing. New York, New York: Cambridge University Press; 2009. pp. 615–30.
Egge JJ, Simons AM. The challenge of truly cryptic diversity: diagnosis and description of a new madtom catfish (Ictaluridae: Noturus). Zoolog Scr. 2006;35:581–95. https://doi.org/10.1111/j.1463-6409.2006.00247.x
Egge JJD, Nicholson PW, Stark AW. Morphological and molecular variation in the least madtom Noturus hildebrandi (Siluriformes: Ictaluridae), a Mississippi embayment endemic: evidence for a cryptic lineage in the Hatchie River. J Fish B. 2015;86:493–526. https://doi.org/10.1111/jfb.12574
Flouri T, Jiao X, Rannala B, Yang Z. Species Tree inference with BPP using genomic sequences and the Multispecies Coalescent. Mol Biol Evol. 2018;35:2585–93. https://doi.org/10.1093/molbev/msy147
Meyer A. Phylogenetics relationships and evolutionary processes in east african cichlid fishes. 1993; 8:279–84. https://doi.org/10.1016/0169-5347(93)90255-N
Mandal SD, Chhakchhuak L, Gurusubramanian G, Kumar NS. Mitochondrial markers for identification and phylogenetic studies in insects – A review. DNA Barcodes. 2014;2:1–9. https://doi.org/10.2478/dna-2014-0001
Gordeeva NV, Volkov AA. A new primer set for amplification of the cytochrome b gene in Lantern Fishes (Myctophidae). J Ichtyol. 2016;899–903. https://doi.org/10.1134/S0032945216060059
Satoh TP, Miya M, Mabuchi K, Nishida M. Structure and variation of the mitochondrial genome of fishes. BMC Genom. 2016;17:719. https://doi.org/10.1186/s12864-016-3054-y
Barrow LS, Ralicki HF, Emme SA, Lemmon EM. Species tree estimation of north american chorus frogs (Hylidae: Pseudacris) with parallel tagged amplicon sequencing. Mol Phylogenet Evol. 2014;75:78–90. https://doi.org/10.1016/j.ympev.2014.02.007
Doadrio I, Dominguez O. Phylogenetic relationships within the fish family Goodeidae based on cytochrome b sequence data. Mol Phylogenet Evol. 2004;31:416–30. https://doi.org/10.1016/j.ympev.2003.08.022
Pérez-Rodríguez R, Domínguez-Domínguez O, Pérez-Ponce de León G, Doadrio I. Phylogenetic relationships and biogeography of the genus Algansea Girard (Cypriniformes: Cyprinidae) of central Mexico inferred from molecular data. BMC Evol Biol. 2009;9:223. https://doi.org/10.1186/1471-2148-9-223
Pérez-Rodríguez R, Domínguez-Domínguez O, Mar-Silva AF, Doadrio I, Pérez-Ponce de León G. The historical biogeography of the southern group of the sucker genus Moxostoma (Teleostei: Catostomidae) and the colonization of central Mexico. Zool J Linn Soc. 2016;177:633–47. https://doi.org/10.1111/zoj.12383
Lambert SM, Reeder TW, Wiens JJ. When do species tree and concatenated estimates disagree? An empirical analysis with higher-level scincid lizard phylogeny. Mol Phylogenet Evol. 2015;82:146–55. https://doi.org/10.1016/j.ympev.2014.10.004
Yu G, Rao D, Matsui M, Yang J. Coalescent-based delimitation outperfors distance-based methods for delineating less divergent specie: the case of Kurixalus odontotarsus species group. Sci Rep. 2017;7:16124. https://doi.org/10.1038/s41598-017-16309-1
Luo A, Ling C, Ho SYW, Zhu CD. Comparison of methods for molecular species delimitation across a range of speciation scenarios. Syst Biol. 2018;67:830–46. https://doi.org/10.1093/sysbio/syy011
Chan KO, Grismer LL. To split or not to split? Multilocus phylogeny and molecular species delimitation of southeast asian toads (family: Bufonidae). BMC Evol Biol. 2019;19:95. https://doi.org/10.1186/s12862-019-1422-3
Leaché AD, Fujita MK. Bayesian species delimitation in west african forest geckos (Hemidactylus fasciatus). Proc Royal Soc B. 2010;277:3071–7. https://doi.org/10.1098/rspb.2010.0662
Hardman M. The phylogenetic relationships among Noturus catfishes (Siluriformes: Ictaluridae) as inferred from mitochondrial gene cytochrome b and nuclear recombination activating gene 2. Mol Phylogenet Evol. 2004;30:395–408. https://doi.org/10.1016/s1055-7903(03)00203-3
Rodiles-Hernández RH, Lundberg J, Sullivan JP. Taxonomic discrimination and identification of Extant Blue catfishes (Siluriformes: Ictaluridae: Ictalurus furcatus Group). Proc Acad Nat Sci Philadelphia. 2010;159:67–82. https://doi.org/10.1635/053.159.0105
Wiens JJ. Missing data and design of phylogenetic. J Biomed Inform. 2006;39:34–42. https://doi.org/10.1016/j.jbi.2005.04.001
Philippe H, Zhou Y, Brinkmann H, Rodrigue N, Delsuc F. Heterothachy amd long-branch attraction in phylogenetics. BMC Ecol Evol. 2005;5:50. https://doi.org/10.1186/1471-2148-5-50
Fuller P, Neilson M. Ictalurus punctatus (Rafinesque, 1818). In: U.S. Geological Survey. Nonindigenous Aquatic Species Database. Gainesville, Florida. 2017 https://nas.er.usgs.gov/queries/FactSheet.aspx?speciesID=740. Accessed 18march 2022.
Rosas-Valdez R, Choudhury A. Pérez-Ponce de León G. Molecular prospecting for cryptic species in Phyllodistomum lacustri (Platyhelminthes, Gorgoderidae). Zool. Scr. 2011. 40; 296–305. https://doi.org/10.1111/j.1463-6409.2011.00472.x
Bean PT, Jackson JT, McHenry DJ, Bonner TH, Fostner MRJ. Rediscovery of the headwater catfish Ictalurus lupus (Ictaluridae) in a western gulf-slope drainage. Southwest Nat. 2011;56:285–9. http://www.jstor.org/stable/23028182
Gutiérrez-Barragán A, García-De León FJ, Varela-Romero A, Ballesteros-Córdova CA, Grijelva-Chon JM, De la Re-Vega E. Evidence of hybridization between Yaqui catfish Ictalurus pricei (rutter, 1896) and channel catfish Ictalurus punctatus (Rafinesque, 1818) in north-west México revealed by analysis of mitochondrial and nuclear genes. Aquat Conserv : Mar Freshw Ecosyst. 2021;1–8. https://doi.org/10.1002/aqc.3709
Meek SE. The fresh-water of the Mexico north of the isthmus of Tehuantepec. Zoological series. Volume 5. Chicago: Field Columbian Museum Publications; 1904.
Soto-Galera E. Fichas de Especies del Proyecto W040. 2003. http://www.conabio.gob.mx/conocimiento/ise/fichasnom/Ictalurusaustralis00.pdf. Accessed 19 March 2022.
Cordero GR, Espinosa PH. Ictalurus australis (Meek, 1904). In: Ceballos G, Díaz-Pardo E, Martínez-Estevéz L, Espinosa-Pérez H, editors. Los Peces dulceacuícolas de México en peligro de Extinción. Ciudad de México: Fondo de Cultura Económica; 2016. pp. 201–2.
Soto-Galera E. Fichas de Especies del Proyecto W040. 2003. http://www.conabio.gob.mx/conocimiento/ise/fichasnom/Ictalurusmexicanus00.pdf. Accessed 19 March 2022.
Espinosa-Pérez H, Huidobro-Campos L, Ictaurus mexicanus. (Meek, 1904): Bagre del Verde. In: Ceballos G, Díaz-Pardo E, Martínez-Estévez L, Espinosa-Pérez H, editors. Los Peces dulceacuícolas de México en Peligro de Extinción. Ciudad de México: Fondo de Cultura Económica; 2016. P. 207–210.
Miller RR. Composition and derivation of the freshwater fish fauna of México. An Esc Nac Cien Biol. 1986;30:121–53.
Obregón-Barbosa H, Contreras-Balderas S, Lozano-Vilano ML. The fishes of northern and central Veracruz, México. Hydrobiologia. 1994;286:79–95.
Cordero RG, Espinosa-Pérez H. Ictaurus australis (Meek, 1904): Bagre del Panuco. In: Ceballos G, Díaz-Pardo E, Martínez-Estévez L, Espinosa-Pérez H, editors. Los Peces dulceacuícolas de México en peligro de Extinción. Ciudad de México: Fondo de Cultura Económica; 2016. pp. 201–2.
García-De León FJ, Hernández-Sandoval AI, Contreras-Catala F, Sánchez-Velasco L, Ruiz-Campos G. Distribution of fishes in the Río Guayalejo-Río Tamesí system and relationships with environmental factors in northeastern Mexico. Environ Biol Fish. 2017;101:167–80. https://doi.org/10.1007/s10641-017-0689-8
Varela-Romero A, Hendrickson DA, Yepiz-Plascencia G. Status of the Yaqui Catfish (Ictalurus pricei) in the United States and Northwestern Mexico. Southwest Nat. 2011;56:278–86. https://doi.org/10.2307/23028181
Álvarez J. Contribución al conocimiento de los bagres fósiles de Chapala y Zacoalco. Jalisco, México DF: Instituto Nacional de Antropología e Historia: Paleontología I.; 1966.
Pérez-Rodríguez R, Domínguez-Domínguez O, Doadrio I, Cuevas-García E, Pérez-Ponce de León G. Comparative historical biogeography of three groups of Nearctic freshwater fishes across central Mexico. J Fish Biol. 2015;86:993–1015. https://doi.org/10.1111/jfb.12611
Domínguez-Domínguez O, Pedraza-Lara C, Gurrola-Sánchez N, Perea S, Pérez-Rodríguez R, Israde-Alcántara I, Garduño-Monroy VH, Doadrio I. Pérez-Ponce de León G, Brooks DR. Historical Biogeography of the Goodeinae (Cyprinodontiforms). In: Uribe-Aranzabal MC, Grier HJ, editors. Viviparous fishes II. Homestead: New Life Publications; 2010. p. 13–30.
Beltrán-López RG, Pérez-Rodríguez R, Montañez-García OC, Artigas-Azas JM, Köck M, Mar-Silva AF, Domínguez-Domínguez O. Genetic differentiation in the genus Characodon: implications for conservation and taxonomy. PeerJ. 2021;9:e11492. https://doi.org/10.7717/peerj.11492
Miller RR, Smith ML. Origin and geography of the fishes of central Mexico. In: Hocutt CH, Wiley EO, editors. The zoogeography of north american freshwater fishes. New York: John Wiley and Sons; 1986. pp. 413–55.
Domínguez-Domínguez O, Vila M, Pérez-Rodríguez R, Remón N, Doadrio I. Complex evolutionary history of the mexican stoneroller Campostoma ornatum Girard, 1856 (Actinopterygii: Cyprinidae). BMC Evol Biol. 2011;11:153. https://doi.org/10.1186/1471-2148-11-153
Schönhuth S, Doadrio I, Domínguez-Domínguez O, Hillis DM, Mayden RL. Molecular evolution of southern North American Cyprinidae (Actinopterygii), with the description of the new genus tampichthys from central Mexico. Mol Phylogenet Evol. 2008;47:729–56. https://doi.org/10.1016/j.ympev.2007.11.036
Schönhuth S, Blum MJ, Lozano-Vilano L, Neely DA, Romero- Varela A, Espinosa-Pérez H, Perdices A, Mayden RL. Interbasin exchange and repeated headwater capture across the Sierra Madre Occidental inferred from the phylogeography of mexican stonerollers. J Biogeogr. 2011;38:1406–21. https://doi.org/10.1111/j.1365-2699.2011.02481.x
Schönhuth S, Perdices A, Lozano-Vilano L, Garcia- De León FJ, Espinosa-Pérez H, Mayden RL. Phylogenetic relationships of north american western chubs of the genus Gila (Cyprinidae, Teleostei), with emphasis on southern species. Mol Phylogenet Evol. 2014;70:210–30. https://doi.org/10.1016/j.ympev.2013.09.021
Corona-Santiago DK, Domínguez-Domínguez O, Tovar-Mora L, Pardos-Blas JR, Herrerías-Diego Y, Pérez-Rodríguez R, Doadrio I. Historical biogeography reveals new independent evolutionary lineages in the Pantosteus plebeius-nebuliferus species-group (Actinopterygii: Catostomidae). BMC Evol Biol. 2018;18:206. https://doi.org/10.1186/s12862-018-1321-z
Hardman M, Hardman LM. The relative importance of body size and paleoclimatic change as explanatory variables influencing lineage diversification rate: an evolutionary analysis of Bullhead catfishes (Siluriformes: Ictaluridae). Syst Biol. 2008;57:116–30. https://doi.org/10.1080/10635150801902193
Lundberg JG. The fossil catfishes of North America. Univ. Mich. Mus. Paleont. Pap. Paleont. 1975. 11: 1–51.
De Sant´ Anna V, Collette BB, Godfrey SJ. Belone countermani, a new Miocene needlefish (Belonidae) from the St. Marys Formation of Valvert Cliffs, Maryland. Proc. Biol. Soc. Wash. 2013; 126:137–150. https://doi.org/10.2988/0006-324X-126.2.137
Cavin L. 2017. Freshwater fishes: 250 million years of evolutionary history. 1rst ed. London: ISTE Press-Elsevier. 2017.
De la Maza-Benignos M, Ornelas-García CP, Lozano-Vilano MDL, García-Ramírez ME, Doadrio I. Phylogeographic analysis of genus Herichthys (Perciformes: Cichlidae), with descriptions of Nosferatu new genus and H. tepehua n. sp. Hydrobiologia. 2015;748:201–31. https://doi.org/10.1007/s10750-014-1891-8
Pérez-Miranda F, Mejía O, Soto‐Galera E, Espinosa‐Pérez H, Piálek L, Říčan O. Phylogeny and species diversity of the genus Herichthys (Teleostei: Cichlidae). J Zool System Evol Res. 2018;56:223–47. https://doi.org/10.1111/jzs.12197
Galloway WE, Whiteaker TL, Ganey-Curry P. History of Cenozoic North American drainage basin evolution, sediment yield, and accumulation in the Gulf of Mexico basin. Geosphere. 2011;7:938–73. https://doi.org/10.1130/GES00647.1
Echelle AA, Carson EW, Echelle AF, Bussche RAVD, Dowling TE, Meyer A. Historical biogeography of the New-World pupfish genus Cyprinodon (Teleostei:Cyprinodontidae). Copeia; 2005: 220–239. https://doi.org/10.1643/CG-03-093R3
Schönhuth S, Hillis DM, Neely DA, Lozano-Vilano L, Perdices A, Mayden RL. Phylogeny, diversity, and species delimitation of the North American Round-Nosed Minnows (Teleostei: Dionda), as inferred from mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2012;62:427–46. https://doi.org/10.1016/j.ympev.2011.10.011
Schönhuth S, Lozano-Vilano L, Perdices A, Espinosa H, Mayden RL. Phylogeny, genetic diversity, and phylogeography of the genus Codoma (Teleostei, Cyprinidae). Zool Scr. 2014;44:11–28. https://doi.org/10.1111/zsc.12083
Lyons TJ, Máiz-Tomé L, Tognelli M, Daniels A, Meredith C, Bullock R, Harrison I, editors, Contreras-MacBeath, Hendrickson T, Arroyave DA, Mercado Silva J, Köck N, Domínguez M, Domínguez O, Valdés González A, Espinosa Pérez H, Gómez Balandra MA, Matamoros W, Schmitter-Soto JJ, Soto-Galera E, Rivas González JM, Vega-Cendejas ME, Ornelas-García CP, Norris S, Mejía Guerrero HO. The status and distribution of freshwater fishes in Mexico. Cambridge, UK and Albuquerque, New Mexico, USA: IUCN and ABQ BioPark. 2020.
Palumbi S, Martin A, Romano S, Macmillan WO, Stice L, Grabowski G. The simple fool’s guide to PCR, Version 2.0. Honolulu: Univ. Hawaii; 1991.
Perdices A, Bermingham E, Montilla A, Doadrio I. Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America. Mol Phylogenet Evol. 2002;25:172–89. https://doi.org/10.1016/S1055-7903(02)00224-5
Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN. DNA barcoding Australia’s fish species. Philos T R Soc B. 2005;360:1847–57. https://doi.org/10.1098/rstb.2005.1716
Perdices A, Doadrio I. The Molecular Systematics and Biogeography of the european Cobitids based on mitochondrial DNA sequences. Mol Phylogenet Evol. 2001;19:468–78. https://doi.org/10.1006/mpev.2000.0900
Quenouille B, Bermingham E, Planes S. Molecular systematics of the damselfishes (Teleostei: Pomacentridae): bayesian phylogenetic analyses of mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2004;31:66–88. https://doi.org/10.1016/S1055-7903(03)00278-1
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4. https://doi.org/10.1093/molbev/msw054
Stephens M, Donelly PA. Comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003;73:1162–9. https://doi.org/10.1086/379378
Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, Sánchez-Gracia A. Mol Biol Evol. 2017;34:3299–302. https://doi.org/10.1093/molbev/msx248. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets.
Guindon S, Gascuel O. A simple, fast and accurate method to estimate large phylogenies by maximum-likelihood. Syst Biol. 2003;52:696–704. https://doi.org/10.1080/10635150390235520
Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772. https://doi.org/10.1038/nmeth.2109
Xia X. DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Mol Biol Evol. 2013;30:1720–8. https://doi.org/10.1093/molbev/mst064
Ronquist F, Teslenko M, Mark PVD, Ayres DL, Darling AS, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42. https://doi.org/10.1093/sysbio/sys029
Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarisation in bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018;67:901–4. https://doi.org/10.1093/sysbio/syy032
Stamatakis A. Phylogenetic models of rate heterogeneity: a high performance computing perspective. In: Proceedings 20th IEEE International Parallel & Distributed Processing Symposium. 2006. http://www.cecs.uci.edu/~papers/ipdps06/index.html. Accessed 25 Mar 2022.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90. https://doi.org/10.1093/bioinformatics/btl446
Stamatakis A. A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008;57:758–71. https://doi.org/10.1080/10635150802429642
Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Proceedings of the Gateway Computing Environments Workshop (GCE). 2010. https://www.computer.org/csdl/proceedings/gce/2010/12OmNy7h3cn. Accessed 25 Mar 2022.
Heled J, Drummond AJ. Bayesian inference of species trees from multilocus data. Mol Biol Evol. 2010;27:570–80. https://doi.org/10.1093/molbev/msp274
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73. https://doi.org/10.1093/molbev/mss075
Yang Z. A tutorial of BPP for species tree estimation and species delimitation. Curr Zool. 2015;61:854–65. https://doi.org/10.1093/czoolo/61.5.854
Yang Z, Rannala B. Bayesian species delimitation using multilocus sequence data. Proc. Natl. Acad. Sci. U.S.A. 2010; 107: 9264–9269. https://doi.org/10.1073/pnas.0913022107
Rannala B, Yang Z. Improved reversible jump algorithms for bayesian species delimitation. Genetics. 2013;194:245–53. https://doi.org/10.1534/genetics.112.149039
Yang Z, Rannala B. Unguided species delimitation using DNA sequence data from multiple loci. Mol Bio l Evol. 2014;31:3125–35. https://doi.org/10.1093/molbev/msu279
We are grateful to Rocio Rodiles Hernández for providing us with tissue samples of Ictalurus meridionalis. We also wish to thank Dr. Francisco J. García de León for providing us with tissue samples for some ictalurid species used as outgroups.
This study was funded by the projects “Sistemática y taxonomía del género Ictalurus (Siluriformes:Pisces), en México (PRODEP-SEP2016/UMSNH-PTC-384)” and “Delimitación de especies dentro del género de peces dulceacuícolas Ictalurus (Siluriformes:Pisces) (CIC-UMSNH-2016-2017)”. This study was partly supported by the Spanish Ministry of Science and Innovation (PID2019-103936GB-C22). The funding body played no role in the design of the study and collection, analysis, interpretation of data, and in writing the manuscript.
Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.
Ethics approval and consent to participate
We did not conduct any animal experiments in this study. As stated in the Material and methods section, most of the tissue samples used in the present study were obtained from three tissue banks, while the few specimens sampled were donated by fishermen
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
About this article
Cite this article
Pérez-Rodríguez, R., Domínguez-Domínguez, O., Pedraza-Lara, C. et al. Multi-locus phylogeny of the catfish genus Ictalurus Rafinesque, 1820 (Actinopterygii, Siluriformes) and its systematic and evolutionary implications. BMC Ecol Evo 23, 27 (2023). https://doi.org/10.1186/s12862-023-02134-w
- Freshwater fishes
- North America
- Taxonomic classification