Mitochondrial DNA variation in natural populations of endangered Indian feather-back fish, Chitala chitala

The concept of "state fishes of India": Current status, knowledge gaps, and strategic plans for conservation and sustainable utilization

The concept of "state fishes of India" highlights the importance and significance of the prioritized fish species distributed across various states within the country. This review article systematically documents the current status of state fishes from various perspectives, identifies the prevailing knowledge gaps, and also highlights the issues and strategic plans essential for the conservation and sustainable utilization of these valuable genetic resources. A total of 8357 publications were checked for the consolidated information on state fish species, and the appropriate items were selected under eight categories: biology, physiology and nutrition, aquaculture, habitat and environmental parameters, genetics and biotechnology, harvest and postharvest, fish health management, and others. The synthesized information was used to present the current status of research and development on state fish species. The knowledge gaps that are to be addressed are also depicted under the perspectives of fisheries management and conservation aquaculture. Based on the findings, strategic plans for the targeted conservation programmes are proposed and discussed under various in situ and ex situ conservation measures. Further, the departmental processes involved in the declaration, the importance of stakeholder involvement, namely, local communities and policymakers, in fostering effective conservation measures, and planning for utilization of these valuable fish genetic resources are also indicated.

Genetic diversity, spatial connectivity, and population structure of Asian silurid catfish Wallago attu (Bloch and Schneider, 1801) in the Ganga River System: insights from mitochondrial DNA analysis

BACKGROUND

The Ganga River System (GRS) is a biodiversity hotspot, its ecological richness is shaped by a complex geological history. In this study, we examined the genetic diversity, spatial connectivity, and population structure of the Asian Silurid catfish, Wallago attu, across seven tributaries of the GRS.

METHODS AND RESULTS

We employed three mitochondrial DNA (mtDNA) regions: cytochrome c oxidase subunit I (COXI), cytochrome b (Cyt b), and control region (CR). Our comprehensive dataset encompassed 2420 bp of mtDNA, derived from 176 W. attu individuals across 19 sampling sites within the seven rivers of GRS. Our findings revealed high gene diversity (Hd:0.99) within W. attu populations. Analysis of Molecular Variance (AMOVA) highlighted that maximum genetic variations were attributed within the populations, and the observed genetic differentiation among the seven populations of W. attu ranged from low to moderate. Network analysis uncovered the presence of three distinct genetic clades, showing no specific association with seven studied rivers. Bayesian skyline plots provided insights into the demographic history of W. attu, suggesting a recent population expansion estimated to have occurred approximately 0.04 million years ago (mya) during the Pleistocene epoch.

CONCLUSIONS

These results significantly enhance our understanding of the genetic diversity and spatial connectivity of W. attu, serving as a vital foundation for developing informed conservation strategies and the sustainable management of this economically valuable resource within the Ganga River System.

Mitochondrial DNA reveals shallow population genetic structure in economically important Cyprinid fish Labeo rohita (Hamilton, 1822) from South and Southeast Asia

BACKGROUND

Assessment of genetic diversity and population genetic structure is important for species that are economically important, threatened, and are at global conservation priority. Analysis of mitochondrial DNA is broadly used in species identification and population genetics studies due to the availability of sufficient reference data and better evolutionary dynamics for phylogeographic investigation. Labeo rohita (Rohu) is an economically important species cultured under carp polyculture systems in Asia. The present study explores the genetic diversity, phylogeography, and population structure of L. rohita from different countries using cytochrome oxidase subunit I (COI) gene.

METHODS AND RESULTS

A total of 17 L. rohita specimens were sampled from River Beas, India. For the genetic study, we amplified and sequenced COI mitochondrial DNA region. The obtained genetic data was combined with 268 COI records available in the NCBI and BOLD databases originating from multiple populations/countries across South and Southeast Asia. As a result, 33 haplotypes were identified that displayed low nucleotide (π = 0.0233) and moderate haplotype diversity (Hd = 0.523). Tajima (D) was found to be negative (P > 0.05), whereas Fu's Fs showed a positive value (P > 0.05). The overall F_ST value between studied populations was 0.481 (P < 0.05).

CONCLUSION

AMOVA analysis indicated higher variation within than among the population examined. The neutrality tests suggested the presence of rare haplotypes and stable demography within studied populations of L. rohita. The Bayesian skyline plot indicated steady population growth until 1 Mya followed by population decline, whereas F_ST values indicated significant genetic differentiation. High heterogeneity was observed in the Pakistan population which could be indicative of long-term isolation and excessive culturing to meet market demands. The present results are the first global comparative analysis of L. rohita and pave the way forward for detailed genomic and ecological studies aimed at the development of improved stock and effective conservation plans. The study also makes recommendations to conserve the genetic integrity of wild species from aquaculture-reared fishes.

Phenotypic variation of Chitala chitala (Hamilton, 1822) from Indian rivers using truss network and geometric morphometrics

Chitala chitala (Hamilton, 1822) is an economically important food fish species occurring throughout Indian rivers, which also has ornamental value. This study focuses on morphological variations in C. chitala from seven river basins across India namely; Son, Tons, Ken, Brahmaputra, Ganga, Gomti and Gandak. A truss network was constructed by interconnecting nine landmarks to generate 36 morphometric variables extracted from digital images of specimens sampled from the study locations. Transformed truss measurements were subjected to principal component analysis (PCA), canonical discriminant function analysis (CDFA) and discriminant analyses of principal components (DAPC). DAPC function coefficients performed much better in capturing the variation pattern and discrimination between the rivers which was not achieved using CDFA. Eight truss variables were identified with significant and highest loading for truss variables on principal components and coefficients on discriminant function from DAPC contributing to maximum variation between the rivers. Performance graph and functional distribution of identified truss variables clearly indicated distinction between the rivers. Thin plate spline analysis and procrustes shape analysis further showed the variation in morphology between specimens across the rivers. The significant parameters differentiating specimens from different rivers were linked to dorsal fin origin, the base of the pectoral fin and the perpendicular point on the anal fin from the dorsal fin origin. Variation in the hydrodynamics of the rivers studied might be possibly affecting the fin kinematics and consequently leading to adaption seen as phenotypic variation in C. chitala. The results showcased in the present study shall help in better understanding of intra-specific diversity which is significant for management and conservation of a species.

Genetic diversity in natural population of Notopterus notopterus evaluated through mitochondrial DNA marker ATPase 6/8 regions and cytochrome b gene

BACKGROUND

Population structure and genetic diversity of bronze featherback Notopterus notopterus, fish was not studied yet from Pakistan. So, genetic diversity and population structure of N. notopterus was analysed using two mitochondrial DNA genetic markers, ATPase 6/8 and Cytochrome b.

METHODS AND RESULTS

150 specimens were collected from five different rivers of Pakistan, resulting 56 haplotypes were detected for ATPase 6/8, Cytb and concatenated gene. Haplotype and nucleotide diversity for ATPase 6/8, Cytb and concatenated gene was observed below 1% among five natural populations of N. notopterus. ATPase 6/8 and Cytb genetic variance among populations was 6% and among and within individuals was 94%. Concatenated genetic variance among populations was 11%, among individual 5% and within individuals 84%. Fst value among all population was found 0.091 (p-value=0.02, p<0.05). The combined data set mean coefficient of genetic differentiation (F_ST) was 0.5572. The pair-wise F_ST was 0.000 (Chenab) to 0.88911 (Ravi). Maximum likelihood phylogeographic history of concatenated gene haplotypes showed four distinct diversified clusters. AMOVA, PCoA (Principal Coordinate Analysis) and maximum likelihood tree indicated that the natural populations of N. notopterus were comprised of four genetic stocks among five Pakistani rivers (Chenab, Indus, Jehlum, Ravi and Satluj).

CONCLUSIONS

This study provides the higher level of genetic diversity with confirmatory proofs among genetic stocks of five natural populations of N. notopterus. The information of genetic diversity and genetic variation, from this research can be utilised to help conserve and manage the species in the wild.

Shallow population structure of southern bluefin tuna, Thunnus maccoyii (Castelnau, 1872) between Indian and Atlantic Oceans inferred from mitochondrial DNA sequences

Southern bluefin tuna (Thunnus maccoyii Castelnau, 1872) is distributed across most of the southern temperate ocean and migrates extensively between 30°S and 50°S. Since T. maccoyii has been continually and heavily exploited, it is necessary to investigate the genetic diversity, population structure and demographic history of T. maccoyii for effective management and conservation. Thirty-seven gonad tissues of T. maccoyii were sampled from two locations, which were in the eastern Indian Ocean and the eastern Atlantic Ocean, by scientific observers onboard Korean T. maccoyii longline vessels in 2015. We compared 1240-bp sequences of combined mitochondrial DNA (mtDNA) from the cytochrome c oxidase subunit I (COI, 504-bp) and control region (CR, 736-bp) sequences. The pairwise fixation index (F_ST) and maximum-likelihood tree showed that two clades (A and B) were formed regardless of locations. Clade A occurred more commonly than clade B in both localities: the occurrence ratio of clade A was 69% in the Indian Ocean, and 79% in the Atlantic Ocean, respectively. Our findings suggest that a historic differentiation event may have occurred in T. maccoyii, but recently the connectivity between the two oceans may be possible in T. maccoyii populations.

Mitochondrial DNA markers reveal genetic connectivity among populations of Osteoglossiform fish Chitala chitala

Genetic diversity and population structure in Indian featherback fish, Chitala chitala (Hamilton, 1822) was investigated by combined analyses of two full mitochondrial genes, ATPase 6/8 and Cytochrome b. A total of 403 individuals, collected from 14 rivers yielded 61 haplotypes. Hierarchical partitioning analysis identified 19.01% variance 'among' and 80.99% variance 'within groups and populations'. The mean coefficient of genetic differentiation (F_ST) was observed to be significant 0.26 (p < 0.05). Mantel tests rejected the hypothesis that genetic and geographic distances were correlated. The patterns of genetic differentiation, AMOVA and principal coordinate analyses indicated that natural populations were sub-structured and comprised of four genetic stocks of C. chitala in Indian rivers. The results also supported the higher resolution potential of concatenated gene sequences. The knowledge of genetic variation and divergence, from this study, can be utilized for its scientific conservation and management in the wild.

Microsatellite marker set for genetic diversity assessment of primitive Chitala chitala (Hamilton, 1822) derived through SMRT sequencing technology

In present study, single molecule-real time sequencing technology was used to obtain a validated set of microsatellite markers for application in population genetics of the primitive fish, Chitala chitala. Assembly of circular consensus sequencing reads resulted into 1164 sequences which contained 2005 repetitive motifs. A total of 100 sequences were used for primer designing and amplification yielded a set of 28 validated polymorphic markers. These loci were used to genotype n = 72 samples from three distant riverine populations of India, namely Son, Satluj and Brahmaputra, for determining intraspecific genetic variation. The microsatellite loci exhibited high level of polymorphism with PIC values ranging from 0.281 to 0.901. The genetic parameters revealed that mean heterozygosity ranged from 0.6802 to 0.6826 and the populations were found to be genetically diverse (F_st 0.03-0.06). This indicated the potential application of these microsatellite marker set that can used for stock characterization of C. chitala, in the wild. These newly developed loci were assayed for cross transferability in another notopterid fish, Notopterus notopterus.

Molecular insights into the genetic and haplotype diversity among four populations of Catla catla from Madhya Pradesh revealed through mtDNA cyto b gene sequences

In the present investigation, the genetic structure of four populations of Catla catla, sequences of mitochondrial gene, cytochrome b (cyto b) from four populations were sequenced and analyzed. The sequences of mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The lowest 249 polymorphic sites and 0.00 parsimony informative sites were detected in populations of Fish Federation Pond (CCFFB) whereas highest 330 polymorphic sites and 56 parsimony informative sites were detected in populations of Narmada River (CCNRH) in the cyto b gene sequences in Catla catla populations. The twelve different haplotypes were detected among the four populations studied, lowest population specific haplotype as 2.00 was observed in Fish Federation Pond (CCFFB) and highest was in Population of Narmada River and Tighra reservoir. Sequencing of cyto b gene revealed 12 number of haplotypes (h) with haplotype (gene) diversity (Hd) 0.8736 and nucleotide diversity (π) 0.6474. These data clearly indicated that, feral/wild population showing highest values of polymorphisms, parsimony, haplotype diversity showing good, healthy habitat is lotic water (Narmada River) and lentic water body (Tighra reservoir). The results also concluded that the partial cyto b is polymorphic and can be a potential marker to determine ecological habitat based genetic differentiation among the populations.

Genetic variations and population structure in three populations of beardless barb, Cyclocheilichthys apogon (Valenciennes, 1842) inferred from mitochondrial cytochrome b sequences

This study aimed to investigate the relation of dispersal barrier and genetic diversity, population structure, and demographic history of 46 samples of beardless barb, Cyclocheilichthys apogon (Valenciennes, 1842) collected from three different locations in North-eastern Thailand. The analysis of molecular variance (AMOVA) was employed in order to determine the genetic variability within and among populations of this fish. The neutrality tests and mismatch distribution analysis were additionally applied to assess the neutrality and demographic expansion of the populations, respectively. Contiguous sequences within range from 1100 to1140 bp were obtained with varying 16 different haplotypes with high-haplotype diversity (0.8773 ± 0.0327) and low-nucleotide diversity (0.00215 ± 0.00020). The variations within and among populations accounted for 98.98 and 1.02% of the total variation, respectively. The low level of pairwise Fst estimations indicated a possible gene flow among populations and a suggestion of genetic homogeneity at this geographical range. A supportive idea of having a single-maternal lineage and past demographic expansion or selection experiencing has distinctly appeared among these populations. The current data suggests that all three populations distinctly exist as a single stock unit and that is an important factor in identifying genetic variation of C. apogon in this geographical area to be used in establishing effective plans and strategies for a conservation and management.

Population genetic structure of Macrobrachium rosenbergii (Palaemonidae) from Indian waters using mitochondrial ATPase 6/8 gene

Macrobrachium rosenbergii, giant freshwater prawn, is one of the most commercially important crustaceans. In the present study, primers for ATPase 6/8 region of mt-DNA were designed and successfully amplified (827 bp) in the species. The nucleotide variation in ATPase 6/8 gene revealed the population structuring in natural populations of M. rosenbergii in Indian waters. A total of 35 haplotypes were observed in 93 individuals collected from different locations. Low nucleotide diversity and high haplotype diversity were noticed for the ATPase 6/8 gene. Significant pairwise F_ST and, haplotype network indicated occurrence of distinct populations. Observed mismatch distribution and Tajima's D test suggested demographical stability of giant freshwater prawn. The genetic stock structure revealed in this study will be helpful for conservation and management of stocks of M. rosenbergii in Indian waters.

Genetic characterization of Bagarius species using cytochrome c oxidase I and cytochrome b genes

In this study, we first inferred the genetic variability of two Bagarius bagarius populations collected from Ganges and Brahmaputra rivers of India using two mtDNA markers. Sequence analysis of COI gene did not show significant differences between two populations whereas cytochrome b gene showed significant differences between two populations. Followed by, genetic relationship of B. bagarius and B. yarrielli was analyzed using COI and cytochrome b gene and the results showed a higher level genetic variation between two species. The present study provides support for the suitability of COI and cytochrome b genes for the identification of B. bagarius and B. yarrielli.

Characterization of mitochondrial ATPase 6/8 genes in wild Labeo calbasu (Hamilton, 1822) and mapping of natural genetic diversity

We characterized mitochondrial ATP synthase (ATPase) 6 and 8 genes in Labeo calbasu (Hamilton, 1822) and determined genetic variation in wild populations across the natural distribution in Indian rivers. A total of 206 individuals were sampled from 11 riverine sites belonging to distinct geographical locations covering five major river basins. Sequencing of 842 base pairs of ATPase 6/8 revealed 21 haplotypes with haplotype diversity ranging from 0.1250 (River Satluj) to 0.8846 (River Bhagirathi). Analysis of molecular variance (AMOVA) of mitochondrial DNA (mtDNA) data revealed significant genetic differentiation among sites (FST = 0.192, p < 0.0001) which was indicative of moderate level of genetic structuring in the wild L. calbasu populations. The patterns of genetic divergence and haplotype network of mtDNA revealed distinct clades present in Indian rivers. The analysis of data demonstrated the potential of ATPase 6/8 genes in determining the genetic diversity and indicated considerable sub-structuring in wild calbasu populations present in different rivers.

Genetic variation in Labeo fimbriatus (Cypriniformes: Cyprinidae) populations as revealed by partial cytochrome b sequences of mitochondrial DNA

Labeo fimbriatus, a medium sized carp is assessed as a commercially important aquaculture species in Indian subcontinent. In the present study, the genetic diversity and population structure of four Indian riverine populations of L. fimbriatus have been evaluated using partial cytochrome b sequences of mitochondrial DNA. Sequencing and analysis of this gene from 108 individuals defined 7 distinct haplotypes. Haplotype diversity (Hd) and nucleotide diversity (π) ranged from 0.067 to 0.405 and 0.00023 to 0.03231, respectively. The Mahanadi population had the highest π level. Analysis of molecular variance (AMOVA) indicated that 47.36% of genetic variation contained within population and 53.76% of genetic variation among groups. Pairwise FST analysis indicated that there was little or no genetic differentiation among populations (-0.0018 to 04572) from different geographical regions except Mahanadi population. The Mahanadi population can be considered as a separate stock from rest three riverine populations. Accordingly, the genetic information generated from this study can be implemented while taking decision in formulating base population for the sustainable selective breeding programs of this species.

Genetic divergence in wild population of endangered yellowtail catfish Pangasius pangasius (Hamilton-Buchanan, 1822) revealed by mtDNA

Pangasius pangasius, an endangered freshwater fish species, is an important component of capture fishery from Indian rivers. Samples collected through commercial catches from three riverine populations were analyzed with cytb (307 bp) and ATPase6&8 (842 bp) regions for population variation and differentiation. The sequences of the both the mitochondrial regions revealed high haplotype and low nucleotide diversity. Shallow genetic diversity based on ATPase6&8 was observed, however its haplotypes network clearly indicated two distinct mitochondrial lineages. Mismatch distribution suggested population bottlenecks followed by expansion in Mahanadi population. The present study indicated the ATPase6&8 to be a potential mitochondrial marker for studying the population sub-structuring in the wild population of P. pangasius.

Genetic divergence in natural populations of bronze featherback, Notopterus notopterus (Osteoglossiformes: Notopteridae) from five Indian rivers, analyzed through mtDNA ATPase6/8 regions

The present study characterized 842 bp fragment of mitochondrial ATP synthase 6 and 8 (ATPase6/8) genes in Notopterus notopterus. In all, 97 samples of N. notopterus were collected from five distant rivers; viz Satluj, Gomti, Yamuna, Brahmaputra and Mahanadi representing 4 river basins in India. The analysis of variation revealed presence of 23 haplotypes in ATPase6/8 gene with haplotype diversity (Hd) of 0.899 and nucleotide diversity (π) of 0.00336. The within population variation which was 41.78% of the total variation of 58.22% was found among population. The Fst value of 0.582 (P < 0.05) of the total population was found significant. The results concluded that the polymorphism in ATPase6/8 gene is a potential marker that is important for determining genetic divergence of wild N. notopterus populations. The findings reveal common ancestry of mahanadi population with the populations in rivers of Indo-Gangetic region. However, long evolutionary isolation must be responsible for the high genetic divergence between N. notopterus in Mahanadi and other regions.

Characterization of polymorphic microsatellite markers and genetic diversity in wild bronze featherback, Notopterus notopterus (Pallas, 1769)

Six polymorphic microsatellite DNA loci were identified in the primitive fish, bronze featherback, Notopterus notopterus for the first time and demonstrated significant population genetic structure. Out of the six primers, one primer (NN90) was specific to N. notopterus (microsatellite sequence within the RAG1 gene) and five primers were product of successful cross-species amplification. Sixty-four primers available from 3 fish species of order Osteoglossiformes and families Notopteridae and Osteoglossidae were tested to amplify homologous microsatellite loci in N. notopterus. Fifteen primer pairs exhibited successful cross-priming PCR product. However, polymorphism was detected only at five loci. To assess the significance of these six loci (including NN90) in population genetic study, 215 samples of N. notopterus from five rivers, viz Satluj, Gomti, Yamuna, Brahmaputra and Mahanadi were analyzed. The five sample sets displayed different diversity levels and observed heterozygosity ranged from 0.6036 to 0.7373. Significant genotype heterogeneity (P < 0.0001) and high FST (0.2205) over all loci indicated that the samples are not drawn from the same genepool. The identified microsatellite loci are promising for use in fine-scale population structure analysis of N. notopterus.

Cytochrome b gene sequence divergence of seven sisorid species of catfish genus Glyptothorax (siluriformes, sisoridae) from India

Samples of seven Sisorid catfish species of the genus Glyptothorax (G. garhwali, G. dakpathari, G. brevipinnis, G. ngapang, G. granulus, G. ventrolineatus, and G. davissinghi) were collected from the Himalayan region and the Western Ghats of India. They were analyzed for the mitochondrial cytochrome b gene (Cyt b). Out of 1152 nucleotide positions analyzed, 269 (23.3%) were found to be variable and 235 (20.3%) were parsimoniously informative. The sequences showed 111 (9.6%) fourfold degenerate sites. The overall transition/transversion bias was R = 3.457. The average proportion of base substitutions measured as P-distance for all sequences of seven Glyptothorax species and other five species comprising four comparison groups was intraspecies, P = 0.17 ± 0.05%, intragenus, P = 10.75 ± 0.48%, intrafamily, P = 20.07 ± 1.43%, and intraorder, P = 21.10 ± 0.45%. Within the Glyptothorax genus maximum divergence was obtained among G. brevipinnis sequences, whereas the least divergence was obtained within G. davissinghi. The phylogenetic trees for 193 and 47 sequences of Sisorid catfishes together were developed using the Cyt b gene and four different analytical approaches: Bayesian (BA), neighbour-joining (NJ), maximum parsimony (MP), and maximum likelihood (ML). The analysis revealed a monophyletic origin for the all investigated representatives of the genus Glyptothorax (99, 100, 99, and 73% support level in our BA, NJ, MP, and ML analyses respectively) and with some reservations for Sisoridae, which is the principal family investigated. The monophyletic origin of the two subfamilies of Sisorid catfish defined in the literature was partly also supported by molecular phylogenetic data.