Pangasiidae, Bleeker, 1858

Duong, Thuy Yen, Pham, Linh Thi Khanh, Le, Xuyen Thi Kim, Nguyen, Ngoc Tran Thi, Nor, Siti Azizah Mohd & Le, Thanh Hoa, 2023, Mitophylogeny of Pangasiid Catfishes and its Taxonomic Implications for Pangasiidae and the Suborder Siluroidei., Zoological Studies 62 (48), pp. 1-17 : 7-11

Phylogenetic relationships within Pangasiidae

The complete mtDNA sequence datasets (concatenated 13 PCGs) recovered Pangasius and Pangasianodon as sister subclades and the Pangasiidae as a monophyletic clade with a 100% bootstrap value ( Fig. 2 ). The ML tree clarified the monophyly of Pangasianodon and Pangasius with 88% bootstrap support. The ML tree clarified the monophyly of Pangasianodon , with 88% bootstrap support. Pangasianodon gigas is identified as a sister taxon to a group of P. sanitwongsei and Pn. hypophthalmus (5 mitogenomes), and Pangasius bocourti (as named in GenBank under no. MN842723) is positioned in the Pn. hypophthalmus group with 100% bootstrap support. Within the Pangasius , P. mekongensis was recovered as a sister taxon to the Indian Pangasius pangasius Hamilton-Buchanan, 1822 ( Hossain et al. 2009) with a very high bootstrap (100%), while P. krempfi was placed in between P. larnaudii and the two abovementioned Pangasius species ( P. mekongensis and P. pangasius ) with a moderate bootstrap value (68%). The most concerning feature is that P. sanitwongsei was positioned as a sister taxon to a subgroup of all the Pn. hypophthalmus sequences (bootstrap 98%), and was in between this subgroup and Pn. gigas ( Fig. 2 ). Several samples may have been misidentified and, therefore, phylogenetically misplaced. In fact, the sequence named “ Pangasianodon_hypophthalmus _(BaijinCo-Foshan)-China-MZ286355” is Pangasius larnaudii , and the “ Pangasius bocourti (QingyuanGD) -China-MN842723” sample is Pangasianodon hypophthalmus . These pangasiid sequences were grouped into their corresponding phylogenetic clades ( Fig. 2 ).

Furthermore, using single-gene datasets, we investigated the close phylogenetic relationships between Pangasius and Pangasianodon . Multiple cox 1 and cyt B barcode sequences are available in GenBank and in previous publications; therefore, we downloaded all cox 1 (551 bp) and cyt B (634 bp) sequences and extracted cox 1 and cyt B, respectively, from the complete mitogenomes (listed in Table S3). The cox 1 and cyt B topologies also revealed that the Pangasiidae were a sister group to the Austroglanididae . All three families ( Pangasiidae , Cranoglanididae , and Austroglanididae ) together with the Ictaluridae formed a large group that was always a sister group to the Ariidae , as discovered in the mitophylogeny constructed based on the complete mitogenome data in our current study and as in the previously reported analysis ( Schedel et al. 2022).

The cox 1 phylogenetic tree ( Fig. 3 ) indicated that P. mekongensis (4 sequences) is a sister taxon to P. pangasius (8 sequences) with relatively high bootstrap support (88%). In the cyt B tree ( Fig. 4 ), this species is sistered with P. pangasius with 100% nodal support. Pangasius krempfi , in the cox 1 tree, was shown to be close to the Helicophagus species ( Helicophagus leptorhynchus and Helicophagus waandersii ) with a low bootstrap of 44%, and was placed as a sister taxon in the cyt B tree in a non-stable phylogenetic relationship with a 63% bootstrap, to P. macronema and P. polyuranodon . The partial cox 1 datasets (four sequences, consisting of two from South Africa, one from China, and one from Cambodia) and cyt B datasets (two sequences from Thailand and one from China) of the correctly identified P. sanitwongsei species placed this taxon into the Pangasius clade, although with low support (bootstrap 36% and 29%, respectively) ( Figs. 3 and 4 ). It should be noted that several P. sanitwongsei sequences, including Psan-(PB2)-( China)-JN020073 and Psan-(PB1)-CN-JN020086, as well as a P. bocourti sequence, Pboc-(QingyuanGD)-China-MN842723, were placed in the Pn. hypophthalmus cluster, which could be attributed to missampling or misidentification ( Figs. 3 and 4 ).

0

.

05

In both the cox 1 and cyt B phylogenies, the Pangasianodon was clearly resolved as a sister to the Pangasius with medium nodal support (bootstrap: 86% for cox 1 and 52% for cyt B). While the monophyletic Pangasianodon was resolved with its two members ( Pn. gigas was the sister taxon to Pn. hypophthalmus in both the cox 1 and the cyt B phylogenetic topologies), the Pangasius cluster included not only Pangasius but Helicophagus and Pseudolais species as well. Though they had low bootstrap values (28% for cox 1 and 29% for cyt B phylogeny), Helicophagus and Pseudolais species were grouped as sister taxa to some of the Pangasius species.

Figures

Fig. 2. PhyML-phylogeny of the order Siluriformes, including 32 catfish families (117 sequences) of three suborders, Siluroidei, Loricarioidei, and Diplomystoidei based on the complete concatenated nucleotide sequences of all 13 mitochondrial protein coding genes (about 11,408 bp in length) (Table S2). Two sequences of Gonorynchiformes were used as an outgroup. The alignment was performed by MAFFT (Katoh and Standley 2013), curated by BMGE v1.12 (Criscuolo and Gribaldo 2010), the tree was reconstructed in PhyML 3.3 (Guindon et al. 2010) using a maximum likelihood method and 1000 bootstrap resamplings, and the output Newick tree was extracted and visualized using FigTree v1.4.4 (Rambaut 2018). The nodal bootstrap support values (shown at each node) were interpreted from the concurrently constructed tree using the above MAFFT-BMGE alignment by MEGA X (Kumar et al. 2018). The basal nodes of the three suborders (Diplomystoidei, Loricarioidei, and Siluroidei) as well as the two major “Big Asia” and “Big Africa” groups (background highlighted) are shown by arrows. The Pangasius mekongensis, Pangasianodon hypophthalmus, and Pangasius krempfi sequences in this study are indicated by stars and with the associated families’ background highlighted. The taxa were presented with their full names. The abbreviations of the isolates are given in brackets, including the geographical origin or voucher records of each sequenced specimen (where available), which were retrieved from the previous studies (Saitoh et al. 2003; Nakatani et al. 2011; Kappas et al. 2016; Zhang et al. 2021; Schedel et al. 2022). The country of origin or geographical regions where the sample was reported are given in full name, if available. Accession numbers are given at the end of each sequence label. The scale bar represents the number of substitutions per site.

Fig. 3. Detailed PhyML-phylogeny based on the analysis of the partial cox1 sequences (551 bp) showing the detailed relationships of the family Pangasiidae and related families (Austroglanididae, Ictaluridae, and Cranoglanididae). In total, 83 sequences, including 81 from Pangasius and Pangasianodon and 2 outgroup sequences from the order Gymnotiformes, were included (Table S3). The alignment was performed by MAFFT (Katoh and Standley 2013), curated by BMGE v1.12 (Criscuolo and Gribaldo 2010), the tree was reconstructed in PhyML 3.3 (Guindon et al. 2010) using a maximum likelihood method and 1000 bootstrap resamplings, and the output Newick tree was extracted and visualized using FigTree v1.4.4 (Rambaut 2018). The basal nodes of the Pangasiidae and two sister groups (Pangasianodon and (Pangasius + Helicophagus + Pseudolais)) are shown by arrows. The Pangasius mekongensis, Pangasianodon hypophthalmus, and Pangasius krempfi sequences in this study are bolded. The taxonmisidentified sequences were added with a question mark at the end. The taxa from Pangasiidae were shortened and those from other related families were presented with their full names. The abbreviations of the isolates are given in brackets, including the geographical origin or voucher records of each sequenced specimen (where available), which were retrieved from the previous studies (Karinthanyakit and Jondeung 2012; Tran and Duong 2019; Schedel et al. 2022). The country of origin or where the sample was reported is given in full or in brackets, if available. Accession numbers are given at the end of each sequence label. The scale bar represents the number of substitutions per site.

Fig. 4. Detailed PhyML-phylogeny based on the analysis of the partial cytB sequences (634 bp) showing the detailed relationships of the family Pangasiidae and related families (Austroglanididae, Ictaluridae, and Cranoglanididae). In total, 80 sequences, including 78 from Pangasius and Pangasianodon and 2 outgroup sequences from the order Clupeiformes, were included (Table S3). The alignment was performed by MAFFT (Katoh and Standley 2013), curated by BMGE v1.12 (Criscuolo and Gribaldo 2010), the tree was reconstructed in PhyML 3.3 (Guindon et al. 2010) using a maximum likelihood method and 1000 bootstrap resamplings, and the output Newick tree was extracted and visualized using FigTree v1.4.4 (Rambaut 2018). The basal nodes of the Pangasiidae and two sister groups (Pangasianodon and (Pangasius + Helicophagus + Pseudolais)) are shown by arrows. The Pangasius mekongensis, Pangasianodon hypophthalmus, and Pangasius krempfi sequences in this study are bolded. The taxonmisidentified sequences were added with a question mark at the end. The taxa from Pangasiidae were shortened and those from other related families were presented with their full names. The abbreviations of the isolates are given in brackets, including the geographical origin or voucher records of each sequenced specimen (where available), which were retrieved from the previous studies. The country of origin or where the sample was reported is given in full or in brackets, if available. Accession numbers are given at the end of each sequence label. The scale bar represents the number of substitutions per site.