Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

A genomic perspective on an old question: Salmo trouts or Salmo trutta (Teleostei: Salmonidae)?

There are particular challenges in defining the taxonomic status of recently radiated groups due to the low level of phylogenetic signal. Members of the Salmo trutta species-complex, which mostly evolved during and following the Pleistocene, show high morphological and ecological diversity that, along with their very wide geographic distribution, have led to morphological description of 47 extant nominal species. However, many of these species have not been supported by previous phylogenetic studies, which could be partly due to lack of significant genetic differences among them, the limited resolution offered by molecular methods previously used, as well as the often local scale of these studies. The development of next-generation sequencing (NGS) and related analytical tools have enhanced our ability to address such challenging questions. In this study, Genotyping-by-Sequencing (GBS) of 15,169 filtered SNPs and mitochondrial DNA (mtDNA) D-loop sequences were combined to assess the phylogenetic relationships among 166 brown trouts representing 21 described species and three undescribed groups collected from 84 localities throughout their natural distribution in Europe, west Asia, and North Africa. The data were analysed using different clustering algorithms (admixture analysis and discriminant analysis of principal components-DAPC), a Bayes Factor Delimitation (BFD) test, species tree reconstruction, gene flow tests (three- and four-population tests), and Rogue taxa identification tests. Genomic contributions of the Atlantic lineage brown trout were found in all major sea basins excluding the North African and Aral Sea basins, suggesting introgressive hybridization of native brown trouts driven by stocking using strains of the Atlantic lineage. After removing the phylogenetic noise caused by the Atlantic brown trout, admixture clusters and DAPC clustering based on GBS data, respectively, resolved 11 and 13 clusters among the previously described brown trout species, which were also supported by BFD test results. Our results suggest that natural hybridization between different brown trout lineages has probably played an important role in the origin of several of the putative species, including S. marmoratus, S. carpio, S. farioides, S. pellegrini, S. caspius (in the Kura River drainage) and Salmo sp. in the Danube River basin. Overall, our results support a multi-species taxonomy for brown trouts. They also resolve some species in the Adriatic-Mediterranean and Black Sea drainages as members of very closely related genomic clusters that may need taxonomic revision. However, any final conclusions pertaining to the taxonomy of the brown trout complex should be based on an integrative approach combining genomic, morphological, and ecological data. To avoid challenges in taxonomy and conservation of species complexes like brown trouts, it is suggested to describe species based on genomic clusters of populations instead of describing species based only on morphologically differentiated single type populations.

Cytogenetic investigation of Arctic char × brook trout F1, F2 and backcross hybrids revealed remnants of the chromosomal rearrangements.

Arctic char (Salvelinus alpinus) and brook trout (Salvelinus fontinalis) hybridize and their offspring is viable and fertile. This may be a real treat for the native European stocks of Arctic char which gene pools might be unintendedly contaminated with the genetic elements of brook trout. On the other hand, hybrids of these two species are appreciated by customers and have some potential for the aquaculture. Moreover, Salvelinus hybrids and backcross individuals are interesting models in the research focused on influence of hybridization on the genomic organization and chromosome rearrangements. Thus, the main goal of the present study was to examine chromosomes of Arctic char × brook trout F1, F2 hybrids and backcross individuals and compare with genomic information concerning parental species to recognize karyotypic changes provoked by the hybridization events. Application of conventional and molecular (FISH) techniques allow to identify characteristic chromosomes for both parental species in the hybrid progeny and show multiplicity of cytotypes among different types of crosses with variability in structure and number of chromosome (81-85) and chromosome arm (99-101). Chromosome fragment was detected in the karyotype of one F1 and one backcross individual and the presence of one triploid (3n) fish was documented. Occurrence of chromosomes containing internally located telomeric sequences (ITS) inherited after brook trout or both parental species was shown in F1 and backcross progeny. Moreover, additional CMA3-positive signal on chromosome from Arctic char pair no. 2 in F1 fish and interstitially located active NOR visible on subtelo-acrocentric (F2 hybrid) and acrocentric (Sf × H individual) chromosomes were detected. Described polymorphic chromosomes together with specific, interstitial location of CMA3-positive found in F2 and Sf × H hybrids and DAPI-positive regions observed in H × Sa fish at different uniarmed chromosomes pair presumably are remnants of chromosomal rearrangements. Provided results strongly indicate that the hybridization process influenced the genome organization in the Salvelinus hybrid progeny.

A new species of trout from the northern drainages of Euphrates River, Turkey (Salmoniformes: Salmonidae).

Salmo fahrettini, a new species, is distributed in the northern tributaries of the Euphrates River. It differs from other Salmo species in adjacent waters by a combination of the following characters: a greyish body; one black spot behind the eye and on the cheek; three to six black spots on the opercle; numerous black spots on the back (missing on the predorsal area), flank and middle part of body, surrounded by a roundish white ring; red spots in the median part of the body, surrounded by a roundish white ring; short and narrow maxilla; increase in the number of black and red spots with an increase in size; adipose fin medium size, no or rarely one red spot at its posterior edge; 109-116 lateral line scales; 27-30 scale rows between dorsal-fin origin and lateral line; 20-23 scale rows between the lateral line and anal-fin origin; maxilla length 8.8-10.0% standard length in males, 8.8-9.6 in females.

Unravelling the stock structure of the Persian brown trout by otolith and scale shape.

The simultaneous use of the scale and otolith morphometry was assessed as a potential tool for the identification of Persian brown trout Salmo trutta stocks of the Lar Lake and five rivers from Lar Basin, Iran. Fourier coefficients (FC) and circularity, rectangularity, roundness, ellipticity and form factor shape indices (SI) were calculated for otolith and scale. Several SIs were significantly different among sites for both structures. Permutational multivariate analysis of variance revealed significant differences between several pairwise comparisons for otolith and scale (FCs and indices separately). Discriminant analysis showed otolith FCs (cross-classification rates: 25-86%) and SI (20-45%) appear to be a relatively acceptable tool to discriminate between several locations. Comparatively, the scale morphometry showed lower discriminatory power (FC = 3-65%; SI = 15-34%), with the exception of SI for Elarm River (60%), Kamardasht River (56%) and Lar Lake (75%). Cross-classification rates improved up to 100% when discriminate analysis incorporating all variables for otolith and scale was performed. The results showed a potential segregation between some water bodies, suggesting that the otolith and scale morphometry could be a useful tool to delimit S. trutta populations in relatively close freshwater environments.

Effects of intracoelomic transmitter implantation on metabolic rate, swimming performance, growth and survival in juveniles of two salmonids.

In this study, we investigated the effects of acoustic tag implantation on standard and routine metabolic rate (SMR and RMR, estimated via oxygen consumption), critical swimming speed (Ucrit ), survival and growth in juveniles of rainbow trout Oncorhynchus mykiss and lake trout Salvelinus namaycush. Tag burdens ranged from 1.8% to 7.5% across the two species. Growth rates in acoustic-tagged fish were equal to or higher than those in other treatments. Acoustic-tagged S. namaycush had a marginally lower Ucrit than controls but that effect was not replicated in the O. mykiss experiment. Tagging did not have clear effects on metabolic rate but there was an interaction whereby SMR and RMR tended to increase with time since surgery in tagged O. mykiss but not in other treatments (the same trend did not occur in S. namaycush). Survival was high across treatments (mean 98% survival among O. mykiss, 97.5% among S. namaycush). There were no statistically significant effects of tag burden (percentage of body mass) except for a weak negative relationship with growth rate (across species) and a weak positive relationship with Ucrit but only in the O. mykiss. Collectively, our findings suggest there were minor, context-dependent effects of acoustic tagging in juvenile S. namaycush and O. mykiss during an eight-week laboratory experiment. Further research will be required to assess whether tagging can cause meaningful behavioural effects in these species in captivity or in the wild and whether there is a tag burden threshold above which deleterious effects consistently occur.

Experimental crossbreeding reveals strain-specific variation in mortality, growth and personality in the brown trout (Salmo trutta).

Although hybridization between populations with low genetic diversity may induce heterosis, it can also lead to reduced fitness of hybrid offspring through outbreeding depression and loss of local adaptations. Using a half-sib mating design, we studied on brown trout (Salmo trutta) how hybridization of migratory hatchery-strain females with males from various strains would affect early mortality, growth and personality in F1 offspring. No differences in mortality or alevin body length were found between the crossing groups by the end of the yolk-sac stage. At later developmental stages, higher mortality and slower growth in one of the geographically distant hybrid groups indicated potential outbreeding depression. The personality component indicating boldness and exploration tendency showed fairly low genetic variation and no phenotypic differences among the crossing groups while the personality component related to freezing behavior indicated stronger freezing responses in the purebred and local cross strain when compared to the two other strains. However, the purebred hatchery strain possessed stronger additive genetic tendency for boldness and explorative behavior, and weaker genetic tendency for freezing behavior, when compared to the wild × hatchery hybrid group. Our results add to the cumulating evidence of risks related to the stocking of fish strains from non-native origins.

An improved environmental DNA assay for bull trout (Salvelinus confluentus) based on the ribosomal internal transcribed spacer I.

The majority of environmental DNA (eDNA) assays for vertebrate species are based on commonly analyzed regions of the mitochondrial genome. However, the high degree of mitochondrial similarity between two species of charr (Salvelinus spp.), southern Dolly Varden (S. malma lordii) and bull trout (Salvelinus confluentus), precludes the development of a mitochondrial eDNA assay to distinguish them. Presented here is an eDNA assay to detect bull trout based on the first ribosomal internal transcribed spacer (ITSI), a nuclear marker. This assay successfully detects bull trout and avoids detecting Dolly Varden as well as brook trout (S. fontinalis), Arctic char (S. alpinus), and lake trout (S. namaycush). In addition, this assay was compared with an extensively used mitochondrial bull trout assay and it was found that the ITSI-based assay produced higher detectability. Our results suggest this assay should out-perform the published mtDNA assay across the range of bull trout, while the added specificity allows reliable bull trout detection in areas where bull trout co-occur with other charr such as Dolly Varden. While clearly a superior assay in this instance, basing assays on ITSI is not without problems. For vertebrates, there are fewer ITSI sequences available than commonly sequenced regions of the mitochondrial genome. Thus, the initial in silico screening of candidate assays must be preceded by much more extensive sampling and sequencing of sympatric or closely related taxa. Further, all copies of the internal transcribed spacers within an individual may not be identical, which can lead to complications. Lastly, the copy number for ITSI varies widely across taxa; the greater detectability associated with this assay cannot be assumed for other species.

The Arctic charr (Salvelinus alpinus) genome and transcriptome assembly.

Arctic charr have a circumpolar distribution, persevere under extreme environmental conditions, and reach ages unknown to most other salmonids. The Salvelinus genus is primarily composed of species with genomes that are structured more like the ancestral salmonid genome than most Oncorhynchus and Salmo species of sister genera. It is thought that this aspect of the genome may be important for local adaptation (due to increased recombination) and anadromy (the migration of fish from saltwater to freshwater). In this study, we describe the generation of a new genetic map, the sequencing and assembly of the Arctic charr genome (GenBank accession: GCF_002910315.2) using the newly created genetic map and a previous genetic map, and present several analyses of the Arctic charr genes and genome assembly. The newly generated genetic map consists of 8,574 unique genetic markers and is similar to previous genetic maps with the exception of three major structural differences. The N50, identified BUSCOs, repetitive DNA content, and total size of the Arctic charr assembled genome are all comparable to other assembled salmonid genomes. An analysis to identify orthologous genes revealed that a large number of orthologs could be identified between salmonids and many appear to have highly conserved gene expression profiles between species. Comparing orthologous gene expression profiles may give us a better insight into which genes are more likely to influence species specific phenotypes.

Genetic and molecular evidence that brown trout Salmo trutta belonging to the Danubian lineage are a single biological species.

We tested for reproductive isolation between Salmo trutta abanticus, S.t. labrax, S.t. caspius and S.t. fario by conducting crosses to produce F1 and F2 offspring. We also estimated the extent of genetic divergence between all three entities by examining sequence variation across the coI, d-loop and cytb mitochondrial genes. All of the F1 cross-types were successfully produced. After 2 years of culturing, F2 generation were produced as well. Fertilization, hatching and survival rates and hatching performance of F1 and F2 generations were evaluated. F2 generation had similar performance to their parent. Fertilization, hatching, larval survival rate and hatchery performance of F1 and F2 generation were similar except pure bred F2 S.t. abanticus. Purebred F1 individuals shared similar coloration patterns and spots with their parents but direction of the hybridization appeared to be decisive on morphology of hybrids. Some of the hybrids exhibit different morphological characters than their parents. Based on partial alignments of the three genes, phylogenetic analysis showed that these S. trutta are gathering within the same clade and appeared as monophyletic group. We found that there were some morphologic and genetic variation among S. trutta subspecies but the degree of variation does not warrant species level recognition. These findings indicate that the four subspecies constitute a single biological entity, corresponding to different morphs of the Danubian lineage. We therefore recommend that S. trutta belonging to Danubian lineage in Turkey be referred to as Salmo trutta and that strains be named according to location, such as Abant, Caspian, Black Sea and Anatolian.

Connections and containers: Using genetic data to understand how watershed evolution and human activities influence cutthroat trout biogeography.

Species with large geographic distributions often exhibit complex patterns of diversity that can be further complicated by human activities. Cutthroat trout (Oncorhynchus clarkii) are one of the most widely distributed freshwater fish species in western North America exhibiting substantial phenotypic and genetic variability; however, fish stocking practices have translocated populations outside of their native range and may have obscured intraspecific boundaries. This study focuses on cutthroat trout populations representing three distinct evolutionary clades that are found intermixed within a contact zone between the Bonneville and upper Snake River watersheds in the western United States. We used mitochondrial and microsatellite genetic data, as well as historical stocking records, to evaluate whether populations of cutthroat trout in the contact zone are native or are introduced. We found significant genetic differentiation and fine-scale genetic population structure that was organized primarily by watershed boundaries. While we detected increased genetic diversity in some areas in close proximity to the greatest number of stocking events, the highly organized population structure both within and between areas of the contact zone indicates that the populations are native to the watersheds. Intermixing of distinct evolutionary lineages of cutthroat trout appears to be the result of historical connections between paleodrainages. Our analyses provide a context for understanding how genetic data can be used to assess the status of populations as native or introduced.

Hatchery Strain Contributions to Emerging Wild Lake Trout Populations in Lake Huron.

Recent assessments indicate the emergence of naturally produced lake trout (Salvelinus namaycush) recruitment throughout Lake Huron in the North American Laurentian Great Lakes (>50% of fish <7 years). Because naturally produced fish derived from different stocked hatchery strains are unmarked, managers cannot distinguish strains contributing to natural recruitment. We used 15 microsatellite loci to identify strains of naturally produced lake trout (N = 1567) collected in assessment fisheries during early (2002-2004) and late (2009-2012) sampling periods. Individuals from 13 American and Canadian hatchery strains (N = 1143) were genotyped to develop standardized baseline information. Strain contributions were estimated using a Bayesian inferential approach. Deviance information criteria were used to compare models evaluating strain contributions at different spatial and temporal scales. The best performing models were the most complex models, suggesting that hatchery strain contributions to naturally produced lake trout varied spatially among management districts and temporally between time periods. Contributions of Seneca strain lake trout were consistently high across most management districts, with contributions increasing from early to late time periods (estimates ranged from 52% to 94% for the late period across 8 of 9 districts). Strain contributions deviated from expectations based on historical stocking levels, indicating strains differed with respect to survival, reproductive success, and/or dispersal. Knowledge of recruitment levels of strains stocked in different management districts, and how strain-specific recruitment varies temporally, spatially, and as a function of local or regional stocking is important to prioritize strains for future stocking and management of the transition process from primarily hatchery to naturally produced stocks.

Flank spot number and its significance for systematics, taxonomy and conservation of the near-threatened Mediterranean trout Salmo cettii: evidence from a genetically pure population.

The number of spots on the flank of the Mediterranean trout Salmo cettii population from Tellesimo Stream in Sicily, Italy, varied between seven and 37, with a median and mode of 21, numbers much lower than those reported in the literature and different from those of other Italian trout populations. This finding could be ascribed to the phenotypic plasticity of the species or, alternatively, could provide evidence to support the different origin of Sicilian autochthonous trout. Given the ease of use of this character for the identification of native fish, it could be particularly important for conservation and management purposes.

Using Linkage Maps as a Tool To Determine Patterns of Chromosome Synteny in the Genus Salvelinus.

Next generation sequencing techniques have revolutionized the collection of genome and transcriptome data from nonmodel organisms. This manuscript details the application of restriction site-associated DNA sequencing (RADseq) to generate a marker-dense genetic map for Brook Trout (Salvelinus fontinalis). The consensus map was constructed from three full-sib families totaling 176 F1 individuals. The map consisted of 42 linkage groups with a total female map size of 2502.5 cM, and a total male map size of 1863.8 cM. Synteny was confirmed with Atlantic Salmon for 38 linkage groups, with Rainbow Trout for 37 linkage groups, Arctic Char for 36 linkage groups, and with a previously published Brook Trout linkage map for 39 linkage groups. Comparative mapping confirmed the presence of 8 metacentric and 34 acrocentric chromosomes in Brook Trout. Six metacentric chromosomes seem to be conserved with Arctic Char suggesting there have been at least two species-specific fusion and fission events within the genus Salvelinus In addition, the sex marker (sdY; sexually dimorphic on the Y chromosome) was mapped to Brook Trout BC35, which is homologous with Atlantic Salmon Ssa09qa, Rainbow Trout Omy25, and Arctic Char AC04q. Ultimately, this linkage map will be a useful resource for studies on the genome organization of Salvelinus, and facilitates comparisons of the Salvelinus genome with Salmo and Oncorhynchus.

A real-time polymerase chain reaction method for the identification of four commercially important salmon and trout species.

A duplex quantitative real-time PCR (qPCR) assay was developed for rapid and accurate identification of four commercially important salmon and trout species (Oncorhynchus keta, Oncorhynchus nerka, Oncorhynchus mykiss, and Salmo salar) commonly used for production process of fish in China. The assays targeting the mitochondrial control region (CR) and 16S rRNA gene were able to simultaneously discriminate four target species and the family Salmonidae from processed as well as fresh fish. The qPCR efficiency of each reaction was calculated according to the standard curve, and the method was validated by amplification DNA extracted from single or artificial mixtures prepared with the reference salmon and trout species. Testing of 11 commercial salmon and trout products by the established qPCR assay demonstrated that it was really a useful and academic technique to identify four commercially important salmon and trout species.

Temporal stability and rates of post-depositional change in geochemical signatures of brown trout Salmo trutta scales.

This study investigates temporal stability in the scale microchemistry of brown trout Salmo trutta in feeder streams of a large heterogeneous lake catchment and rates of change after migration into the lake. Laser-ablation inductively coupled plasma mass spectrometry was used to quantify the elemental concentrations of Na, Mg, Mn, Cu, Zn, Ba and Sr in archived (1997-2002) scales of juvenile S. trutta collected from six major feeder streams of Lough Mask, County Mayo, Ireland. Water-element Ca ratios within these streams were determined for the fish sampling period and for a later period (2013-2015). Salmo trutta scale Sr and Ba concentrations were significantly (P < 0·05) correlated with stream water sample Sr:Ca and Ba:Ca ratios respectively from both periods, indicating multi-annual stability in scale and water-elemental signatures. Discriminant analysis of scale chemistries correctly classified 91% of sampled juvenile S. trutta to their stream of origin using a cross-validated classification model. This model was used to test whether assumed post-depositional change in scale element concentrations reduced correct natal stream classification of S. trutta in successive years after migration into Lough Mask. Fish residing in the lake for 1-3 years could be reliably classified to their most likely natal stream, but the probability of correct classification diminished strongly with longer lake residence. Use of scale chemistry to identify natal streams of lake S. trutta should focus on recent migrants, but may not require contemporary water chemistry data.

[Genetic Connectivity Between Sympatric Populations of Closely Related Char Species, Dolly Varden Salvelinus malma and White Char Salvelinus albus].

The closely related chars Salvelinus malma and Salvelinus albus, which sympatrically inhabit the Kamchatka River basin and Kronotsky Lake (Kamchatka), attract the attention of the researchers because of their debated origin and taxonomic status. Previous studies of sympatric populations of these chars revealed small but statistically significant genetic differences between these species at a number of molecular markers, suggesting the presence of the genetic exchange and hybridization. In this study, based on genotypic characterization of nine microsatellite loci, a considerable level of historical and contemporary genetic migration between sympatric populations of these chars was demonstrated. At the individual level a high degree of hybridization was observed, mainly among the Dolly Varden individuals from the studied populations. The obtained evidence on the genetic connectivity between sympatric S. malma and S. albus do not support the separate species status of S. albus.

Hiding in Plain Sight: A Case for Cryptic Metapopulations in Brook Trout (Salvelinus fontinalis).

A fundamental issue in the management and conservation of biodiversity is how to define a population. Spatially contiguous fish occupying a stream network have often been considered to represent a single, homogenous population. However, they may also represent multiple discrete populations, a single population with genetic isolation-by-distance, or a metapopulation. We used microsatellite DNA and a large-scale mark-recapture study to assess population structure in a spatially contiguous sample of Brook Trout (Salvelinus fontinalis), a species of conservation concern. We found evidence for limited genetic exchange across small spatial scales and in the absence of barriers to physical movement. Mark-recapture and stationary passive integrated transponder antenna records demonstrated that fish from two tributaries very seldom moved into the opposite tributary, but movements between the tributaries and mainstem were more common. Using Bayesian genetic clustering, we identified two genetic groups that exhibited significantly different growth rates over three years of study, yet survival rates were very similar. Our study highlights the importance of considering the possibility of multiple genetically distinct populations occurring within spatially contiguous habitats, and suggests the existence of a cryptic metapopulation: a spatially continuous distribution of organisms exhibiting metapopulation-like behaviors.

The complete mitogenome of brown trout (Salmo trutta fario) and its phylogeny.

The complete mitochondrial genome of Salmo trutta fario, commonly known as brown trout, was sequenced using NGS technology. The mitochondrial genome size was determined to be 16 677 bp and composed of 13 protein-coding gene (PCG), 22 tRNAs, 2 rRNA genes, and 1 putative control region. The overall mitogenome composition of S. trutta fario is A: 28.13%, G: 16.44%, C: 29.47%, and T: 25.96% with A + T content of 54.09% and G + C content of 45.91%. The gene arrangement and the order are similar to other vertebrates. The phylogenetic tree constructed using 42 complete mitogenomes of Salmonidae fishes confirmed the position of the present species under the genus Salmo of subfamily Salmoninae. NGS platform was proved to be a rapid and time-saving technology to reveal complete mitogenomes.

Morphometric variation and molecular characterization of snow trout species from Kashmir valley, India.

There is a significant taxonomic ambiguity among snow trout species due to their morphometric similarities. In view of this, a morphometric and molecular study was conducted on five different species of genus Schizothorax that have been reported from Kashmir valley. Morphometric data analyzed using multivariate statistics (Principal component analysis and cluster analysis) indicated the significant grouping of species to individual clusters. Mitochondrial DNA cytochrome oxidase I (COI) gene analysis revealed 0.2%-4.5% genetic divergence among the five species. This study confirms that utility of cytochrome oxidase I in species delineation along with morphometric data. Phylogenetic tree obtained using Neighbor-Joining method revealed that all the five species represented distinct species group. The Schizothorax genus formed two distinct clades; one containing S. niger, S. curvifrons and S. plagiostomus, while other clade containing S. esocinus and S. labiatus. This phlogeny trend was also supported by cluster analysis of morphometric characters. The phylogenetic analysis with other published COI sequences revealed distinct nature of these five species. The study may aid in the taxonomic identification of snow trout species in India. This may further increase the knowledge of the ichthyologists in planning conservation and management strategies for these important fish species along with their natural habitat.