4.1. Comparative morphometric evaluation
The present study used an integrated taxonomic approach to resolve the identity and distribution of Mystus tengara. Comparative morphological evaluation of freshly collected specimens of M. tengara, from West Bengal, showed close similarity to the original description [10], in having four longitudinal stripes separated by 3 pale interspaces, presence of large tympanic spot above pectoral fin, length of four barbels longer than head, and occipital process reaching basal bone of dorsal fin. Specimens of M. tengara collected from Maharashtra also match with the description of Hamilton for M. tengara, except with the presence of small interspace between occipital process and dorsal fin base.
M. tengara is differentiated from M. vittatus (10) by the absence of serrations in dorsal spine (vs. presence), a character which suggested to be an error [27], based on Gunther’s description of “Macrones tengara”. M. tengara further differentiated from M. vittatus [26] by median longitudinal groove reaching to base of occipital process and occipital process reaching basal bone of dorsal fin vs. median longitudinal groove reaching midway behind the hind edge of the eye and base of the occipital process and short interspace between occipital process and basal bone of dorsal fin. In the present study, we observed variations in median longitudinal groove with size and geographical locations.
The re-description of M. tengara [11] to establish and confirm its taxonomic identity and differentiated this species from M. vittatus in having a longer maxillary barbel length (254.5–360.5% HL vs. 214.3–244.9% HL),dorsal spine length (12.3–17.2% SL vs. 10.7–12.2) and length of median longitudinal groove (reaching base of occipital process in juvenile or reaching anterior one third of supraoccipital bone in adult vs. terminating at the anterior border of supraoccipital bone not invading the supraoccipital region). In the present study, specimens of M. tengara, from Eastern India (West Bengal), showed variations from specimens collected from Western India (Maharashtra) in maxillary barbel length, dorsal spine length (14.23–17.79% SL vs. 11.44–13.31) and extent of occipital process (reaching dorsal fin base vs. not reaching) but, without much overall variationsfrom the morphometric descriptionof M. tengara and M. vittatus [11, 24], which shows that the effects of environmental variations in this trait cannot be ruled out.
4.2. DNA barcoding and phylogenetic study
DNA barcodes have been used for confirming identity of species and their distribution [35]. Previous studies have shown that genetic divergence value of 2–3% at DNA barcoding gene (COI) could be used as threshold value to discriminate species [36, 37]. Accordingly, conspecific individuals show genetic divergence value of < 3%, while congeneric species > 3%. During the present study, in Clade I, sequences identified as M. vittatus, M. carcio and M. horai from different geographical locations, were clustered with ‘M. tengara’ (collected in the present study) having a genetic distance of < 3%. This observation suggests that the sequences identified and labelled as M. carico/ M. vittatus/ M. horai in GenBank could be misidentifications of M. tengara. The original description of M. carcio(10) which distinguishes it from M. tengara in the length of maxillary barbel (extending beyond pectoral vs. reaching to end of caudal and serrations on dorsal spine- presence vs absence).
The redescription of M. carcio [19] where it is further distinguish it from M. tengara and M. vittatus based on shorter adipose-fin base length (8.5–11.9% SL vs. 24.0–-31.7 and 21.5–26.9% SL respectively) and posterior fontanel length (reaching base of supraoccipital process vs. not reaching middle of supra occipital bone vs. terminating at the anterior tip of supraoccipital respectively). The description of median longitudinal groove in M. vittatus- terminating at the anterior border of supraoccipital bone, not invading the supraoccipital region [19], is not in agreement with [26] whose description of median longitudinal groove in adult specimens of M. tengara not reaching beyond the middle of supra occipital bone, is also not in agreement to the present study. However, M. carcio, re-described by the authors [19], was distinct in other characters from M. tengara identified in the present study. The re-description of M. carcio [19] was based on specimens collected from Assam, Tripura and Bangladesh but without any molecular evidence. In our phylogenetic analysis, specimens identified as M. carcio from Assam was grouped together with M. tengara sensustricto, whereas, specimens identified as M. carcio, from Bangladesh, grouped together with M. tengara [1] clades with significant genetic variation to be considered distinct.
Further, specimens collected from West Bengal and Maharashtra could be distinguished in morphometric characters such as dorsal spine length (14.23–17.79 vs. 11.44–13.31), pectoral spine length (18.76–20.19 vs. 16.55–17.77), eye diameter (25.05–26.8 vs. 20.06–26.00), maxillary barbel length (reaching posterior tip of anal fin base or to caudal fin base in smaller specimens vs reaching anal fin base) and occipital process (reaching basal bone of dorsal fin vs. a considerable gap present), which reveals geographical variations in these diagnostic phenotypic characters. Similar to our observations, there is a geographical variation in maxillary barbel length [26] in M. tengara from Punjab and Assam (reaching to middle of the pectoral fin vs. reaching to the base of pelvic fin). These finding clearly indicates that these characters may be influenced by size of the fish and the environment, in which they inhabit and cannot be considered as good diagnostic characters for these species [3].
4.3. Molecular evidence reveals cryptic species
In clade II, sequences identified and labelled as M. tengara could likely be misidentifications as the genetic distance of these sequences with those in clade I are higher than 3%. M. tengara recorded from Assam (MH156942) also showed higher genetic distance with sequences in clade I and could be a distinct species. Clade III comprised of M. vittatus and this species was confirmed to be distributed in northern, north-eastern, western and central India. Interestingly, specimens of Mystus cf. tengara, the focus of the present study clustered with M. vittatus recorded from north-east India, and formed clade IV. Though a sister group relationship was observed between clade III and clade IV, the average genetic distance between these two clades was 10.1% suggesting the occurrence of cryptic species in this group. Species names and identities in clade V are also likely to be erroneous due to morphological ambiguities. Studies on generating reference DNA barcodes without morphological taxonomy could often lead to species misidentification [38–39], has been demonstrated recently in hillstream loaches of the Western Ghats [40]. Due to overlapping diagnostic characters and morphological similarities, various authors could have misidentified M. tengara, M. vittatus and M. carcio resulting in the deposition of erroneous sequences in NCBI GenBank.
Based on morphological and genetic evidence of freshly collected M. tengara, from its type locality, we consider sequences that form part of clade I to be M. tengara sensustricto, with the distribution extending from North East India to West Bengal, North India, Central India, Northern peninsular India and Bangladesh. Further confirmation on the identity of M. vittatus and M. carcio, by integrated taxonomic approach based on freshly specimens collected from type locality, is required. The observation of paraphyletic subclade and evaluation of genetic distance between subclades reveals that there could be at least four cryptic species in this group, opening up avenues for future research on the group.