Complete mitochondrial genome of critically endangered catfish Hemibagrus punctatus (Jerdon, 1849) and comparative analysis for insights into the phylogeny of hemibagrids through mitogenomic approach

BACKGROUND

Hemibagrus punctatus (Jerdon, 1849) is a critically endangered bagrid catfish endemic to the Western Ghats of India, whose population is declining due to anthropogenic activities. The current study aims to compare the mitogenome of H. punctatus with that of other Bagrid catfishes and provide insights into their evolutionary relationships.

METHODS AND RESULTS

Samples were collected from Hemmige Karnataka, India. In the present study, the mitogenome of H. punctatus was successfully assembled, and its phylogenetic relationships with other Bagridae species were studied. The total genomic DNA of samples was extracted following the phenol-chloroform isoamyl alcohol method. Samples were sequenced, and the Illumina paired-end reads were assembled to a contig length of 16,517 bp. The mitochondrial genome was annotated using MitoFish and MitoAnnotator (Iwasaki et al., 2013). A robust phylogenetic analysis employing NJ (Maximum composite likelihood) and ASAP methods supports the classification of H. punctatus within the Bagridae family, which validates the taxonomic status of this species. In conclusion, this research enriches our understanding of H. punctatus mitogenome, shedding light on its evolutionary dynamics within the Bagridae family and contributing to the broader knowledge of mitochondrial genes in the context of evolutionary biology.

CONCLUSIONS

The study's findings contribute to a better understanding of the mitogenome of H. punctatus and provide insights into the evolutionary relationships within other Hemibagrids.

Complete mitochondrial genome of golden variant of freshwater fish Labeo rajasthanicus (Cypriniformes: Cyprinidae): endemic to India

The complete mitochondrial genome of the freshwater fish species Labeo rajasthanicus was obtained, using Illumina NovaSeq 6000 with 2 × 150 bp paired-end sequencing. The mitogenome of L. rajasthanicus is 16,738 bp in length (GenBank accession no.: OQ834146), comprised of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a control region, i.e. D-loop. The arrangement of genes was found to be identical to other Cypriniformes fish mitogenome, available in the NCBI database. The taxonomic status of L. rajasthanicus as a valid species was debated by some researchers and it was considered a synonym of L. boggut. However, phylogenetic analysis in the present study supports the species validity of L. rajasthanicus, as it showed a distinct node well separated from L. boggut and supported by a high bootstrap value. Furtherly, the pairwise genetic divergence among studied species showed the divergence between L. rajasthanicus and L. boggut as 1.6% whereas the minimum divergence was found to be 0.13% with L. dussumieri followed by L. fimbriatus (0.58%) and L. gonius (0.63%). The complete mitogenome of L. rajasthanicus will also be useful as a baseline reference genome for the reconstruction and annotation of the mitogenome of other Labeo species.

Anthropogenic risk assessment of riverine habitat using geospatial modelling tools for conservation and restoration planning: a case study from a tropical river Pranhita, India

The riverine ecosystem provides multiple benefits to human community and contributes to the sustainable development of the ecoregion. The growing dependency on these ecosystems has largely contributed to aggravating the ecological risks, habitat degradation, and loss of ecosystem services. The present study evaluates the ecological risk emanating from nine anthropogenic stressors including river use, hydro-morphology, catchment pollution, and biological stressor on river Pranhita in Godavari Basin of Peninsular India using InVEST (Integrated Valuation of Ecosystem Services and Tradeoffs) Habitat Risk Assessment model. The primary field survey, remote sensing, and secondary data-assisted spatial modelling results revealed low ecological risk (R = 0.65 of 3) in river Pranhita due to anthropogenic activities. Sediment loading, the inflow of nitrogen, and habitat fragmentation were the major stressors with relatively higher risk score (> 1); influence on a sizeable portion of riverine habitat (29-75% of the total area under high-risk zone) indicates the mounting threat from catchment activities. The low-risk value observed in protected river reaches as compared to unprotected areas is likely to be influenced by the abundant presence of intact riparian vegetation which mitigate the catchment stressors and minimal anthropogenic activity within protected areas. This study demonstrates the application of InVEST HRA model for ecological risk assessment of riverine ecosystems and fish assemblages along with their input data generation framework. This has the potential for prioritization of sensitive habitats based on computed ecological risk and stressor identification based on their exposure and consequences for developing appropriate mitigation measures. This model is spatially explicit and accommodates user-defined criteria for ecosystem-level assessment at a regional and national scale to facilitate the resource managers and policymakers for conservation and restoration planning and implementation of targeted management measures for sustainable development.

Quantitative assessment of sediment delivery and retention in four watersheds in the Godavari River Basin, India, using InVEST model - an aquatic ecosystem services perspective

Sediment export and retention are important ecosystem processes in any landscape causing soil erosion and sediment loading in waterways consequently affecting the health of aquatic habitats downstream. The present study quantifies sediment export and retention in four watersheds, viz., Hivra, Satrapur, Konta, and Jagdalpur in the Godavari River Basin, India, using Integrated Valuation of Ecosystem Services and Trade-offs (InVEST) ecosystem service modelling tool. The  results revealed that the sediment export yields ranged from 0.75 (Hivra) to 2.77 t/ha/year (Jagdalpur). The mean deviation between modelled values and observed sediment export yield was - 11.11%, which indicated good prediction by the model. The sediment retention ranged from 16.04 (Hivra) to 101.52 t/ha/year (Konta). Most sediment export and retention occurred on cropland or shrubland land use land cover types in all four watersheds. For decision making on soil conservation, soil loss tolerance limits have been established for these watersheds. For aquatic habitats, sediment concentration is considered more important than the total annual sediment export, since water turbidity is an important determinant of water quality, and the aquatic lives therein. Therefore, the temporal distribution of rainfall and corresponding sediment export becomes important, since these two factors determine the sediment concentration as well as turbidity in the waterbody. In current study, "Precipitation Concentration Index adjusted Sediment Export Yield Index" was developed to account for the effects of the temporal rainfall distribution and its impact on sediment export. The index for four watersheds was quantified (Hivra > Satrapur > Konta > Jagdalpur), which is concordant to the turbidity values reported by respective gauge stations. Thus, the proposed index can efficiently capture the impact of temporal rainfall distribution on sediment export, and consequently its effect on water turbidity. The study revealed the potential of InVEST model to quantify the sediment export and retention in the watersheds studied. Together with the proposed index, it would help the policy makers in making informed decisions for planning conservation strategies for aquatic biodiversity.

Development of novel microsatellite markers for population differentiation and detection of natural selection in wild populations of butter catfish, Ompok bimaculatus (Bloch, 1794)

BACKGROUND

Butter catfish (Ompok bimaculatus) is a preferred species in South East Asia, with huge aquaculture potential. However, there is limited information about genetic stock composition due to insufficient markers. The goal of this study was to develop de novo microsatellite markers.

METHODS AND RESULTS

For sequencing, genomic SMRT bell libraries (1.5 Kbp size) were prepared for O. bimaculatus. A total of 114 SSR containing sequences were used for primer designing. Polymorphic loci were validated by genotyping 83 individuals from four distant riverine populations, viz., Brahmaputra, Bichiya, Gomti and Kaveri. A total of 30 microsatellite loci were polymorphic, of which five were found to be associated with functional genes and eight (four positive and four negative) loci were found to be under selection pressure. A total of 115 alleles were detected in all loci and PIC ranged from 0.539 to 0.927 and pair-wise F_ST values from 0.1267 to 0.26002 (p < 0.001), with an overall F_ST value of 0.17047, indicating the presence of population sub-structure. Cross-species transferability of 29 loci (96.67%) was successful in congener species, Ompok pabda.

CONCLUSION

The novel SSR markers developed in this study would facilitate stock characterization of natural populations, to be used in future selection breeding programs and planning conservation strategies in these species. Identified non-neutral markers will give insights into the effect of local adaptation on genetic differentiation in the natural population of this species.

Evaluating the influence of environmental variables on the length-weight relationship and prediction modelling in flathead grey mullet, Mugil cephalus Linnaeus, 1758

Fish stocks that are grown under diverse environmental conditions have different biometric relationships and growth patterns. The biometric length-weight relationship (LWR) is an essential fishery assessment tool, as fish growth is continuous and depends on genetic and environmental factors. The present study attempts to understand the LWR of the flathead grey mullet, Mugil cephalus Linnaeus, 1758, from different locations. The study area encompassed its distribution in the wild across freshwater location (one), coastal habitats (eight locations), and estuaries (six locations) in India to determine the relationship between various environmental parameters. Specimens (n = 476) of M. cephalus were collected from commercial catches and the length and weight of individual specimens were recorded. Monthly data from the study locations were extracted for nine environmental variables from the datasets downloaded from the Physical Oceanography Distributed Active Archive Center (PO.DAAC) and the Copernicus Marine Environment Monitoring Service (CMEMS) over 16 years (2002 to 2017) on the Geographical Information System platform. The parameters of the LWR, intercept 'a' and slope or regression coefficient 'b', varied from 0.005321 to 0.22182 and 2.235 to 3.173, respectively. The condition factor ranged from 0.92 to 1.41. The partial least squares (PLS) score scatter plot matrix indicated differences in the environmental variables between the locations. PLS analysis of the regression coefficient and environment parameters revealed that certain environment variables viz., sea surface temperature, salinity, dissolved oxygen, nitrate, and phosphate, played a positive role. However, chlorophyll, pH, silicate, and iron played a negative role in influencing weight growth across various locations. The results revealed that the M. cephalus specimens from three locations, Mandapam, Karwar, and Ratnagiri, possessed significantly higher fitness to their environment than those from the other six locations. The PLS model can be used to predict weight growth under the various environmental conditions of different ecosystems. The three identified locations are useful sites for the mariculture of this species considering their growth performance, the environmental variables, and their interactions. The results of this study will improve the management and conservation of exploited stocks in regions affected by climate change. Our results will also aid in making environment clearance decisions for coastal development projects and will improve the efficiency of mariculture systems.

Transcriptome Analysis Revealed Osmoregulation Related Regulatory Networks and Hub Genes in the Gills of Hilsa shad, Tenualosa ilisha, during the Migratory Osmotic Stress

Tenualosa ilisha (Hilsa shad), an anadromous fish, usually inhabits coastal and estuarine waters, and migrates to freshwater for spawning. In this study, large-scale gill transcriptome analyses from three salinity regions, i.e., fresh, brackish and marine water, revealed 3277 differentially expressed genes (DEGs), out of which 232 were found to be common between marine vs freshwater and brackish vs freshwater. These genes were mapped into 54 KEGG Pathways, and the most significant of these were focal adhesion, adherens junction, tight junction, and PI3K-Akt signaling pathways. A total of 24 osmoregulatory genes were found to be differentially expressed in different habitats. The gene members of slc16 and slc2 families showed a dissimilar pattern of expressions, while two claudin genes (cldn11 & cldn10), transmembrane tm56b, and voltage-gated potassium channel gene kcna10 were downregulated in freshwater samples, as compared to that of brackish and marine environment. Protein-protein interaction (PPI) network analysis of 232 DEGs showed 101 genes to be involved in PPI, while fn1 gene was found to be interacting with the highest number of genes (36). Twenty-five hub genes belonged to 12 functional groups, with muscle structure development with seven genes, forming the major group. These results provided valuable information about the genes, potentially involved in the molecular mechanisms regulating water homeostasis in gills, during migration for spawning and low-salinity adaptation in Hilsa shad. These genes may form the basis for the bio-marker development for adaptation to the stress levied by major environmental changes, due to hatchery/culture conditions.

Description of a new Pangasius (Valenciennes, 1840) species, from the Cauvery River extends distribution range of the genus up to South Western Ghats in peninsular India

A new species of the genus Pangasius, is described based on 17 specimens collected from the Cauvery River, India. It can be distinguished from its sister species from South and Southeast Asia, by its widely placed, small and rounded vomerine and palatine tooth plates, longer maxillary and mandibular barbels, greater vertebrae count 50 (vs. 44-48), and smaller caudal peduncle depth (6.5-8.2% SL vs. 9.89-13.09% SL). The tooth plates of the new species closely resembles that of Pangasius macronema but can be clearly distinguished from the latter by having lesser gill rakers (16-19 vs. 36-45); a smaller eye (2.4-4.4% SL vs. 5.2-9.6% SL); and larger adipose-fin base (1.5-2.9% SL vs. 0.1-1.2% SL). The mitochondrial cytochrome c oxidase (COI) gene sequence of the new species shows the genetic divergence of 3.5% and 5.1% from P. pangasius and P. silasi respectively, the two sister species found in South Asia and India. The species delimitation approaches, Poisson Tree Processes (PTP) and assemble species by automatic partitioning (ASAP) clearly resolved that the P. icaria is distinct from its sister species. Phylogenetic position of the species with its sister species was evaluated using maximum likelihood and Bayesian analysis. The discovery of this previously unknown species of genus Pangasius from the Cauvery River of peninsular India indicates important biogeographical insight that this genus migrated till the southern division of Western Ghats.

Applications of Sentinel-2 satellite data for spatio-temporal mapping of deep pools for monitoring the riverine connectivity and assessment of ecological dynamics: a case from Godavari, a tropical river in India (2016-2021)

Rivers are dynamic ecosystems with diverse habitats that require ample connectivity to ensure the flow of ecosystem services, thus empowering the sustainable development of an entire basin. Geo-spatial tools offer powerful prospects for monitoring of aquatic ecosystems. The usefulness of Sentinel-2 datasets to assess river connectivity has been explored for an un-gauged seasonal river system. The present study was undertaken in five ecologically unique river reaches viz. Wainganga, Wardha, Pranhita, Godavari-mid and Manair in Godavari Basin in the Indian Deccan Plateau to map water spread dynamics at various time scales, i.e., fortnightly, monthly, seasonal, annual and demi-decadal during 2016-2021. The maximum value of perennial water spread per square kilometre of total floodplain area (2016-2021), determined using Sentinel-2 imageries, was observed in river Wardha (0.18) followed by Pranhita (0.12) and Wainganga (0.11). The water spread showed a decreasing trend, while the number of patches in the river corridor increased over time from post-monsoon to pre-monsoon season. The copious perennial habitat with relatively larger patches, incessant flow in river Pranhita and obstructed flow, large-sized patches reported in river Wardha during summer months, hold importance in terms of providing refuge to aquatic biota. This study provides evidence for the impact of water projects on spatio-temporal water spread dynamics in Godavari Basin. The demonstrated utility of Sentinel-2 imagery coupled with gauge station measurements for river continuity assessment and deep pool mapping would aid in enhancing our understanding on environmental flow at a spatial scale, which in turn would aid in effective river management to achieve the Sustainable Development Goals. The implications of this study for sustainable environmental management and limitations are also discussed.

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 and morphological assessment of a vulnerable large catfish, Silonia silondia (Hamilton, 1822), in natural populations from India

Silonia silondia is a commercially important fish distributed in Asian countries, which is under threat due to overexploitation. This study focuses on the morphological analysis and genetic variation of S. silondia individuals, through truss network and sequencing of two mitochondrial regions, respectively, from six wild populations of the Ganga and Mahanadi river systems in India. A total of 38 haplotypes was observed by analysing combined mitochondrial genes (cytochrome b + ATPase 6/8) in 247 individuals of S. silondia collected from six populations. Average haplotype and nucleotide diversities were 0.8508 and 0.00231, respectively. Genetic structure analysis showed the predominant cause of genetic variation to be within populations. The two clades were observed among the haplotypes and time of divergence from their most probable ancestor was estimated to be around 0.3949 mya. Analysis of combined mitochondrial genes in six populations of S. silondia resulted into three management units or genetic stocks. The truss network analysis was carried out by interconnecting 12 landmarks from digital images of specimens to identify phenotypic stocks. Sixty-five truss morphometric variables were analysed for geometric shape variation which revealed morphological divergence in River Son specimens. The present study presents molecular markers and genetic diversity data which can be critical input for conservation and management of differentiated populations and future monitoring of the genetic bottleneck. The morphological shape analysis clearly shows that variation in the insertion of adipose fin is an important parameter influencing the morphological discrimination.

Microsatellite marker development in Spanish mackerel Scomberomorus commerson using third generation sequencing technology

Spanish mackerel S. commerson belonging to family Scombridae, represent a group of highly commercial marine fisheries with an ever-growing demand world over. Analysing the genetic diversity of this species is of utmost importance and necessary for conservation purposes. Microsatellites are molecular tools with advantages that are ideal for population analyses. This study provides the first multiplex panel set of species-specific microsatellite loci for S. commerson that can be applied when assessing both intra- and inter population genetic variation. Microsatellite marker panels were developed in S. commerson, using Third Generation Sequencing technology in PacBio RSII, based on Single-Molecule Real-Time (SMRT). Thirty- two microsatellite loci were isolated and characterized for S. commerson, by genotyping 20 individuals each obtained from the Kochi and Veraval in the Arabian sea and Chennai along Bay of Bengal coast (n = 3). The number of alleles per locus in S. commerson varied from 4 to 17, while the mean observed and expected heterozygosities ranged from 0.656 to 0.753. The Polymorphic Information Content (PIC) were highly informative, 85% loci with PIC value 0 > 0.75. This suite of markers provides the first species specific nuclear multiplex microsatellite marker panels (32 loci) for S. commerson and thus allows assessment of different populations structures of the species across its distribution range, with more specificity. These newly developed loci have also been validated for cross transferability in another scomberid fish Scomberomorus guttatus.

Transcriptome analysis reveals immune pathways underlying resistance in the common carp Cyprinus carpio against the oomycete Aphanomyces invadans

Infection with Aphanomyces invadans is a serious fish disease with major global impacts. Despite affecting over 160 fish species, some of the species like the common carp Cyprinus carpio are resistant to A. invadans infection. In the present study, we investigated the transcriptomes of head kidney of common carp experimentally infected with A. invadans. In time course analysis, 5288 genes were found to be differentially expressed (DEGs), of which 731 were involved in 21 immune pathways. The analysis of immune-related DEGs suggested that efficient processing and presentation of A. invadans antigens, enhanced phagocytosis, recognition of pathogen-associated molecular patterns, and increased recruitment of leukocytes to the sites of infection contribute to resistance of common carp against A. invadans. Herein, we provide a systematic understanding of the disease resistance mechanisms in common carp at molecular level as a valuable resource for developing disease management strategies for this devastating fish-pathogenic oomycete.

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.

Resumption of donor-origin spermatogenesis in senescent goldfish Carassius auratus (Linnaeus, 1758) following spermatogonial cell therapy

Stem cell research has come into prominence because of its applications in assisted reproductive technology and the treatment of deadly diseases. In teleost fishes, spermatogonial stem cells have been effectively used to produce surrogate gametes and progeny through germ cell transplantation technique. The present study is the first report of an innovative application of stem cell therapy in fish species for revitalising the reproductive competence of senescent individuals. Senescent male goldfish, Carassius auratus, approximately 10 years of age were procured from a fish-breeding farm and were reared locally in the lab for an additional two years. The senescence of the individuals was then evaluated and confirmed using histological analysis, gonadal index assessment, and germ-cell specific vasa gene expression. Analyses revealed absence of spermatogonial cells and other germ cells in the testes of the senescent fish (n = 5). Spermatogonial cells from sexually immature C. auratus male donor were isolated using discontinuous percoll gradients, labelled with the fluorescent dye PKH-26, and transplanted into the gonads of senescent C. auratus males through urogenital papilla. Six months after the transplant, spermatozoa were collected through applying gentle manual pressure on the abdomen and were observed under a microscope. All C. auratus males with the transplant had produced spermatozoa from the transplanted cells. This was confirmed by the retention of PKH-26 in the spermatozoa and diagnostic SSR locus. Gravid C. auratus females were artificially inseminated with the spermatozoa of those senescent males and natural spawning was allowed. As a result viable progeny were produced. These observation suggests that the reproductive competence of senescent male fishes can be revitalised through spermatogonial stem cell therapy to produce functional gametes.

Non-invasive Method for Collection of Clarias magur Spermatozoa for Breeding

Induced spawning is more inefficient in Clarias magur than in other fish species such as cyprinids and salmonids. Ovulation can be induced in the female C. magur by using pituitary extract or synthetic hormones. However, milt from the male C. magur cannot be obtained by hand stripping because the volume of seminal fluid in the testes is extremely low. Notably, similar problems are observed in other male catfishes such as C. gariepinus and C. lazera. Because milt from the males cannot be obtained for use in artificial fertilization of eggs, males are invariably sacrificed, and the testicular tissue is excised and macerated to obtain spermatozoa.•We developed an alternative approach that allowed harvesting of C. magur spermatozoa through surgery for artificial fertilization without sacrificing male fish.•The surgically obtained spermatozoa were used to inseminate C. magur eggs; the cross resulted in healthy progeny with a fertilization rate of 80%-98% and hatching of up to 85% of fertilized embryos; similar to those obtained using the conventional sacrificial approach (hatching percentage range of 45%-85%) [1].•This indicated the viability of partial surgical harvest of testicular tissue in seed production in C. magur for aquaculture without sacrificing male fish.

Draft genome assembly of Tenualosa ilisha, Hilsa shad, provides resource for osmoregulation studies

This study provides the first high-quality draft genome assembly (762.5 Mb) of Tenualosa ilisha that is highly contiguous and nearly complete. We observed a total of 2,864 contigs, with 96.4% completeness with N₅₀ of 2.65 Mbp and the largest contig length of 17.4 Mbp, along with a complete mitochondrial genome of 16,745 bases. A total number of 33,042 protein coding genes were predicted, among these, 512 genes were classified under 61 Gene Ontology (GO) terms, associated with various homeostasis processes. Highest number of genes belongs to cellular calcium ion homeostasis, followed by tissue homeostasis. A total of 97 genes were identified, with 16 GO terms related to water homeostasis. Claudins, Aquaporins, Connexins/Gap junctions, Adenylate cyclase, Solute carriers and Voltage gated potassium channel genes were observed to be higher in number in T. ilisha, as compared to that in other teleost species. Seven novel gene variants, in addition to claudin gene (CLDZ), were found in T. ilisha. The present study also identified two putative novel genes, NKAIN3 and L4AM1, for the first time in fish, for which further studies are required for pinpointing their functions in fish. In addition, 1.6 million simple sequence repeats were mined from draft genome assembly. The study provides a valuable genomic resource for the anadromous Hilsa. It will form a basis for future studies, pertaining to its adaptation mechanisms to different salinity levels during migration, which in turn would facilitate in its domestication.

Depletion of endogenous germ cells in striped catfish Pangasianodon hypophthalmus (Sauvage, 1878) by heat-chemical treatments

Germ cell (GC) transplantation (GCT) is a proven powerful reproductive technique to enhance the production efficiency of domesticated animals and aid to the recovery of endangered germ lines. In mammals, several methods have been adopted for the eradication of GCs such as treatment with cytotoxic drugs, irradiation, cold ischaemia and hyperthermic treatment. Some of these methods have also been tried in fishes, and conditions for sterilization of gonads have been established. Here, we report the production of GC-depleted male striped catfish Pangasianodon hypophthalmus in 12 weeks by the combination of heat and chemical treatments. The cytotoxic drug busulphan (40 mg/kg) was intraperitoneally injected into the animals at 2-week intervals (six doses in total) and maintained in water at 38°C between weeks 1 and 12. The effectiveness of the treatment was assessed using gonadal index and histology. At the end of 12 weeks, very severe gonadal degeneration was observed in fish treated with the heat-chemical combination, and 100% of sampled fish (n = 5) were found devoid of endogenous GCs. On contrary, high temperature alone caused minor gonadal degeneration. Results obtained in this study suggest that endogenous GCs of large-bodied fish such as P. hypophthalmus can also be sterilized by heat and chemical treatments within a considerably short period.

Tissue specific alpha-2-Macroglobulin (A2M) splice isoform diversity in Hilsa shad, Tenualosa ilisha (Hamilton, 1822)

The present study, for the first time, reported twelve A2M isoforms in Tenualosa ilisha, through SMRT sequencing. Hilsa shad, T. ilisha, an anadromous fish, faces environmental stresses and is thus prone to diseases. Here, expression profiles of different A2M isoforms in four tissues were studied in T. ilisha, for the tissue specific diversity of A2M. Large scale high quality full length transcripts (>0.99% accuracy) were obtained from liver, ovary, testes and gill transcriptomes, through Iso-sequencing on PacBio RSII. A total of 12 isoforms, with complete putatative proteins, were detected in three tissues (7 isoforms in liver, 4 in ovary and 1 in testes). Complete structure of A2M mRNA was predicted from these isoforms, containing 4680 bp sequence, 35 exons and 1508 amino acids. With Homo sapiens A2M as reference, six functional domains (A2M_N,A2M_N2, A2M, Thiol-ester_cl, Complement and Receptor domain), along with a bait region, were predicted in A2M consensus protein. A total of 35 splice sites were identified in T. ilisha A2M consensus transcript, with highest frequency (55.7%) of GT-AG splice sites, as compared to that of Homo sapiens. Liver showed longest isoform (X1) consisting of all domains, while smallest (X10) was found in ovary with one Receptor domain. Present study predicted five putative markers (I-212, I-269, A-472, S-567 and Y-906) for EUS disease resistance in A2M protein, which were present in MG2 domains (A2M_N and A2M_N2), by comparing with that of resistant and susceptible/unknown response species. These markers classified fishes into two groups, resistant and susceptible response. Potential markers, predicted in T. ilisha, placed it to be EUS susceptible category. Putative markers reported in A2M protein may serve as molecular markers in diagnosis of EUS disease resistance/susceptibility in fishes and may have a potential for inclusion in the marker panel for pilot studies. Further, challenging studies are required to confirm the role of particular A2M isoforms and markers identified in immune protection against EUS disease.

Ontogeny of the digestive tract in stinging catfish, Heteropneustes fossilis (Bloch) larvae

Heteropneustes fossilis (Bloch) is an important candidate species for diversification of freshwater aquaculture in India. However, high mortality rate during larval rearing is the most serious bottleneck in commercial production of this species. A proper understanding of the ontogenic development of digestive system provides the basis to understand the nutritional physiology of larvae and develop appropriate feeding strategies. In the present study, the ontogenical development of the digestive tract in H. fossilis larvae was studied from hatching until 30 day post-hatching (dph) at 29 °C. At hatching (2.8 ± 0.2 mm standard length, SL), the digestive tract was undifferentiated and attached dorsally to the yolk sac. At 1 dph (2.9 ± 0.2 mm SL), the mouth opened and oral valves were visible. At 2 dph (3.0 ± 0.3 mm SL), goblet cells were observed in the buccoparyngaeal cavity. At this age, exogenous feeding started and the intestine was differentiated into the anterior and posterior regions, and the rudimentary liver and pancreas were also seen. Small supranuclear vacuoles were observed in the enterocytes of the posterior intestine at 2 dph. Zymogen granules were observed in acinar cells of pancreas by 3 dph, and islets of Langerhans were visible at 4 dph (3.5 ± 0.1 mm SL). At the same age, most of the yolk sac reserves were consumed, whereas they were completely exhausted by 5 dph (3.9 ± 0.5 mm SL). Between 4 and 6 dph, the liver elongated in size and started to accumulate lipids in the hepatocytes. Gastric glands were detected at 4 dph, and the pyloric sphincter was completely differentiated at 9 dph (6.1 ± 0.4 mm SL) as an epithelial fold that separated stomach from the anterior intestine. By 13 dph (8.6 ± 0.2 mm SL), profuse gastric glands were visible inside longitudinal mucosal folds of the stomach. The formation of gastric glands and their development were noticed as the last events in the development of the digestive tract in H. fossilis. This indicated the end of the larval period and the commencement of the juvenile stage. Considering these observations, it is suggested that H. fossilis larvae have a morphologically complete digestive tract by 13 dph. The findings of the study on the development of the digestive system in H. fossilis may help in synchronising the larval stage of development and feeding strategies and would be helpful in improving larval rearing techniques for catfish species.

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 and morphological evidences resolve taxonomic ambiguity between Systomus sarana sarana (Hamilton, 1822) and S. sarana subnasutus (Valenciennes) and suggest elevating them into distinct species

Taxonomic ambiguity exists in genus Systomus and recently many new species were described under this genus. Systomus sarana subnasutus is considered a valid subspecies of S. sarana sarana although revisions have been done by some researchers. We employed a combination of morpho-meristics and molecular tools (Cytochrome c oxidase I, 16S and Cytochrome b genes of mitochondrial genome) to resolve the two species. Three morpho-meristic characters, head length/maxillary barbel length (HL/MxBL), Lateral Line Scales (LLSs) as well as two truss-based characters, had discernible variation between the two taxa. The sequence analysis (2353 nucleotides) depicted a separate clad of S. sarana subnasutus with high bootstrap support. The findings from combined use of morphology, meristics and mitogenes were concordant. The corroborative results suggest the possibility of two different species. The results suggest to adopt suitable management measures, accordingly.

Genetic divergence in natural and farm populations of Pengba fish, Osteobrama belangeri (Valenciennes, 1844), an endemic fish of North-East India derived from mtDNA ATPase 6/8 gene

Pengba fish, Osteobrama belangeri is a freshwater inhabitant, highly endemic, threatened and economically important minor carp for its food value. In the present investigation, population structure of O. belangeri was examined using mitochondrial ATPase 6/8 gene from geographically distinct locations along the Indo-Burma biodiversity hotspot. A total of 80 individuals were collected belonging to natural and farm populations. The hierarchical analysis of molecular variance and conventional Fst values (0.825 in ATPase 6/8, p < .05) indicated significant genetic structure among populations. The result showed that ATPase 6/8 genes are potential marker in determining the genetic divergence between natural and farm populations of O. belangeri from North-East India.

The complete mitochondrial genome of olive barb, Systomus sarana sarana (Hamilton, 1822) and its phylogenetic status

Systomus sarana sarana, commonly known as olive barb, is an important food and ornamental fish species with wide distribution in South Asia. Here, the complete mitogenome was sequenced on HiSeq 2500. With 16,590 nucleotides, the base composition was 32.9% (A), 26% (C), 15.4% (G) and 25.7% (T), comparable with other carps. The clustering pattern depicted the monophyly of S. sarana sarana with sister cyprinids.

Characterization of threatened endemic fish Osteobrama belangeri (Valenciennes, 1844) and related species from North-East India based on morphological and molecular analysis

Four species of the genus Osteobrama collected from five different sampling sites of Chindwin and Barak river basin of North-East India, namely O. belangeri, O. cotio, O. cunma and O. feae, were characterized. Meristic study showed differences in anal fin rays count and lateral line scales. Morphometric analysis revealed significant differences among Osteobrama species. In Truss analysis, a clear pattern of differentiation was observed among the four species with discriminant function analysis assigning 100% correctly to the particular species. In molecular analysis, four Osteobrama species collected were barcoded with COI and 16S rRNA sequences and phylogenetically these four species formed two distinct clusters, O. belangeri form one separate cluster, from the other three species i.e. O. cotio, O. cunma and O. feae. High-resolution clusters generated (NJ trees) aided the groupings of species corresponding to their genera and families which are in confirmation to the values generated by COI Automatic Barcode Gap Discovery bioinformatics platform.

Hypoxic stress -responsive genes in air breathing catfish, Clarias magur (Hamilton 1822) and their possible physiological adaptive function

The Indian catfish, Clarias magur (previous name C. batrachus) is an air breathing fish, inhabitant of aquatic bodies characterized by low dissolved oxygen levels. It is exposed to hypoxic conditions in its natural habitat. Thus, it can be useful model to study the mechanism of hypoxia stress tolerance. In C. magur, molecular processes facilitating its adaptation to hypoxia stress remain largely unexplored, in part due to unavailability of genomic resources. The suppression subtractive hybridization technique (SSH) was employed to compare the differential expression of transcripts under experimental hypoxic conditions, to that of normoxic conditions. Twelve subtracted cDNA libraries (six each forward and reverse) were constructed from brain, heart, liver, muscle, spleen and head kidney tissues. A total of 2020 clones were screened and sequenced, resulting into 1805 high quality expressed sequence tags (ESTs). Annotation of these differentially expressed ESTs resulted into the identification of genes involved in vast majority of pathways/processes affecting metabolism, cellular processes, signal transduction and/or immune functions. Additionally, 18 potential novel genes expressed in hypoxia stress exposed fish were also identified. The study had catalogued the differentially expressed genes from hypoxia stress induced C. magur, where most of them are reported for the first time in a hypoxia-tolerant fish species. The results not only provided insights for the hypoxia stress altered cellular functions in C. magur, but also generated a valuable functional genomics resource to assist targeted studies on functional genomics and future genome projects.

Expression profile and in silico characterization of novel RTF2h gene under oxidative stress in Indian catfish, Clarias magur (Hamilton 1822)

During the investigation of genes involved in the hypoxia tolerance, novel transcript, Replication Termination Factor 2 homologue (RTF2h), was found to be differentially expressed in brain of Clarias magur (previous name, C. batrachus) whose function was still undefined. Thus, present study was aimed to examine the transcriptional response of novel RTF2h gene, for its possible involvement in hypoxia tolerance in C. magur. Novel transcripts expressed under hypoxic stress were identified from ESTs obtained from SSH libraries of C. magur. Homology analysis of novel transcript (JK487668) revealed it to have RING superfamily signature and was closely related with RTF2 homologue. To characterize the expression pattern of CmRTF2h gene in different tissues of C. magur, qRT-PCR analysis was performed which showed significant increased expression, in spleen following short-term hypoxia exposure (12 H; 2.33-fold), while after long-term hypoxia exposure, significant up-regulation was found three tissues: in spleen (6.57-fold), liver (2.31-fold) and head kidney (2.99-fold) and was down-regulated in brain (3.52-fold). Further, the consensus sequence, obtained from SSH EST sequence and transcripts from hypoxia induced transcriptome data, through multiple sequence alignments, homology modeling and phylogentic analysis together, confirmed it as a novel CmRTF2h gene. In-silico protein-protein interaction and docking studies suggested its closely related function to that of RTF2. The present study reports the expression pattern of CmRTF2h under hypoxia. The up-regulated expression of CmRTF2h under hypoxic conditions may contribute to defense mechanism against induced DNA damage, which in turn may of adaptive significance for hypoxia tolerance in C. magur.

Complete mitogenome of Indian mottled eel, Anguilla bengalensis bengalensis (Gray, 1831) through PacBio RSII sequencing

Complete mitogenome sequence for Anguilla bengalensis bengalensis (family Anguillidae) was generated through third-generation sequencing platform. The 16 714 bp mitgenome sequence contained 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and a non-coding (control) region. The gene order was identical to that observed in most of the other vertebrates. The comparison of complete mitogenome sequence of Indian mottled eel generated during this study with two other subspecies did not agree with the taxonomic status of the three subspecies and considered as one species.

DNA barcoding Indian freshwater fishes

DNA barcoding is a promising technique for species identification using a short mitochondrial DNA sequence of cytochrome c oxidase I (COI) gene. In the present study, DNA barcodes were generated from 72 species of freshwater fish covering the Orders Cypriniformes, Siluriformes, Perciformes, Synbranchiformes, and Osteoglossiformes representing 50 genera and 19 families. All the samples were collected from diverse sites except the species endemic to a particular location. Species were represented by multiple specimens in the great majority of the barcoded species. A total of 284 COI sequences were generated. After amplification and sequencing of 700 base pair fragment of COI, primers were trimmed which invariably generated a 655 base pair barcode sequence. The average Kimura two-parameter (K2P) distances within-species, genera, families, and orders were 0.40%, 9.60%, 13.10%, and 17.16%, respectively. DNA barcode discriminated congeneric species without any confusion. The study strongly validated the efficiency of COI as an ideal marker for DNA barcoding of Indian freshwater fishes.

Hypoxia induced altered expression of heat shock protein genes (Hsc71, Hsp90α and Hsp10) in Indian Catfish, Clarias batrachus (Linnaeus, 1758) under oxidative stress

Heat shock proteins (Hsps) are typically associated with stress response and tolerance. The Indian catfish, Clarias batrachus, is a freshwater air-breathing hypoxia tolerant teleost and is potentially important catfish species for aquaculture and for its economic value as food. The present study aimed at determining the transcriptional response of three Hsps, CbHsc71, CbHsp90α and CbHsp10 in hypoxia tolerant Indian catfish, C. batrachus, under experimental and natural hypoxia. The expression profile of above three genes were studied under different periods of hypoxia, through qRT-PCR. Primers were designed from ESTs obtained through SSH libraries constructed from hypoxia treated fishes. The Hsp10 ESTs and deduced protein was in silico characterized for its ORF and for its physical and chemical properties, respectively, using GeneScan, blastp, scanprosite, superfamily and other softwares. A phylogenetic tree was constructed based on deduced amino acid sequences of Hsc71, Hsp90α, Hsp90β of Homo sapiens and other fishes along with CbHsp10 protein in MEGA4. The deduced protein sequences of CbHsp10 was found to have characteristic Hsp10 family signatures, and it is proposed for inclusion of methionine in the consensus sequences of Hsp10 family signature, after the "proline" residue. At transcription level, these genes were found to be differentially regulated under hypoxia stress, in different tissues of C. batrachus. The CbHsc71 and CbHsp90α were up-regulated after short and long-term hypoxia, whereas CbHsp10 was significantly down-regulated after short-term hypoxia. The differential expression of these Hsps may play a role in protection and survival under hypoxia induced oxidative stress in C. batrachus.

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.

Genealogy and phylogeography of Cyprinid fish Labeo rohita (Hamilton, 1822) inferred from ATPase 6 and 8 mitochondrial DNA gene analysis

ATPase 6/8 gene (842 bp) of mitochondrial DNA was sequenced in Labeo rohita samples (n = 253) collected from nine rivers belonging to four river basins; Indus, Ganges, Brahmaputra and Mahanadi. Analysis revealed 44 haplotypes with high haplotype diversity (Hd) 0.694 and low nucleotide diversity (π) 0.001. The within population variation was larger (83.44%) than among population differences (16.56%). The mean FST value (0.166; P &lt; 0.05) for overall populations revealed moderate level of genetic structuring in the wild L. rohita populations. The haplotype network presented a single clade for wild L. rohita population, from different rivers. Negative values for Fu’s index (FS), mismatch distribution analysis indicated period of expansion in L. rohita population. The time after recent expansion was estimated for each population, between 0.042 to 0.167 mya. The pattern of Isolation by Distance (IBD) was not significant (r = -0.113, P &lt; 0.287), when all the sampling locations were compared (Mantel test), however, when an outlier (Indus, Brahmaputra and Mahanadi) was removed from the whole population set, a clear positive correlation between pairwise FST and geographic distance (Km) was seen. The analysis of data demonstrated that ATPase6/8 gene polymorphism is a potential marker to understand genetic population structure of wild L. rohita existing in different rivers. The study identified population substructure in wild L. rohita with common ancestral origin [Current Zoology 60 (4): 460–471, 2014].

Cross-priming of microsatellite loci in subfamily cyprininae (family Cyprinidae): their utility in finding markers for population genetic analysis in three Indian major carps

This study is aimed to identify polymorphic microsatellite markers and establish their potential for population genetics studies in three carp (family cyprinidae; subfamily cyprininae) species, Labeo rohita, Catla catla and Cirrhinus mrigala through use of cyprinid primers. These species have high commercial value and knowledge of genetic variation is important for management of farmed and wild populations. We tested 108 microsatellite primers from 11 species belonging to three different cyprinid subfamilies, Cyprininae, Barbinae and Leuciscinae out of which 63 primers (58.33%) successfully amplified orthologous loci in three focal species. Forty-two loci generated from 29 primers were polymorphic in these three carp species. Sequencing of amplified product confirmed the presence of SSRs in these 42 loci and orthologous nature of the loci. To validate potential of these 42 polymorphic loci in determining the genetic variation, we analyzed 486 samples of three focal species collected from Indus, Ganges and Brahmaputra river systems. Results indicated significant genetic variation, with mean number of alleles per locus ranging from 6.80 to 14.40 and observed heterozygosity ranging from 0.50 to 0.74 in the three focal species. Highly significant (P < 0.00001) allelic homogeneity values revealed that the identified loci can be efficiently used in population genetics analysis of these carp species. Further, thirty-two loci from 19 primers were useful for genotyping in more than one species. The data from the present study was compiled with cross-species amplification data from previous results on eight species of subfamily cyprininae to compare cross-transferability of microsatellite loci. It was revealed that out of 226 heterologous loci amplified, 152 loci that originated from 77 loci exhibited polymorphism and 45 primers were of multispecies utility, common for 2-7 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.

Physiological responses to acute experimental hypoxia in the air-breathing Indian catfish, Clarias batrachus (Linnaeus, 1758)

With an aim to study the mechanism of adaptation to acute hypoxic periods by hypoxia-tolerant catfish, Clarias batrachus, the mass-specific metabolic rate (VO2) along with its hematological parameters, metabolic response and antioxidant enzyme activities were studied. During progressive hypoxia, C. batrachus was found to be an oxyconformer and showed a steady decline in its aquatic oxygen consumption rate. When C. batrachus was exposed for different periods at experimental hypoxia level (0.98 +/- 0.1 mg/L, DO), hemoglobin and hematocrit concentrations were increased, along with decrease in mean cellular hemoglobin concentration, which reflected a physiological adaptation to enhance oxygen transport capacity. Significant increase in serum glucose and lactate concentration as well as lactate dehydrogenase activity was observed. Antioxidant enzymes were found to operate independently of one another, while total glutathione concentration was unaffected in any of the tissues across treatments. These observations suggested that hypoxia resulted in the development of oxidative stress and C. batrachus was able to respond through increase in the oxygen carrying capacity, metabolic depression and efficient antioxidant defense system to survive periods of acute hypoxia.

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.

Complete mitochondrial genome sequences of two endangered Indian catfish species, Clarias batrachus and Pangasius pangasius

Complete mitogenome sequences for Clarias batrachus (magur) and Pangasius pangasius (family Claridae and Pangasiidae, respectively) were generated, which were 16,511 and 16,476 bp, respectively. Both the mitogenomes contained 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs and 1 non-coding (control) region. The order of the genes was identical to that observed in most other vertebrates. The genome resource of complete mitogenome sequence of Indian catfish species generated during this study would be useful in phylogenetic studies.

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.

Molecular characterization and expression analysis of three hypoxia-inducible factor alpha subunits, HIF-1α, -2α and -3α in hypoxia-tolerant Indian catfish, Clarias batrachus [Linnaeus, 1758]

The present study aimed at characterization of three HIF-α subunits, HIF-1α -2α and -3α from hypoxia-tolerant Clarias batrachus, as well as to elucidate their expression pattern under short and long-term hypoxic conditions and identification of biomarker candidate. The complete cDNAs of HIF-1α, -2α and -3α were 2,833, 4,270 and 3,256 bp in length, encoding 774, 818 and 628 amino acid residues, respectively. In C. batrachus, HIF-α subunits were structurally similar in DNA binding, dimerization, degradation and transcriptional activation domains, but differed in their oxygen-dependent degradation domains. Presence of c-Jun N-terminal kinase binding domain in HIF-α subunits was reported here for the first time in fish. In adult C. batrachus, three HIF-α mRNAs were detected in different tissues under normoxic conditions, however HIF-1α was highly expressed in all the tissues studied, in comparison to HIF-2α and -3α. Short-term hypoxia exposure caused significant increase in three HIF-α transcripts in brain, liver and head kidney, while after long-term hypoxia exposure, significant up-regulation of HIF-1α in spleen and -2α in muscle was observed and HIF-3α significantly down-regulated in head kidney. These observations suggest that the differential expression of HIF-α subunits in C. batrachus was hypoxic time period dependent and may play specialized roles in adaptive response to hypoxia. HIF-2α, with its highly elevated expression in muscle tissues, can be a robust biomarker candidate for exposure to hypoxic environment.

Complete sequence and characterization of mitochondrial genome in great snakehead, Channa marulius (Hamilton, 1822)

The complete nucleotide sequence of mitogenome of the great snakehead, Channa marulius (Channidae), was determined and found to be 16,569 base pairs in length. The content and arrangement of different genes on the mitogenome was found similar to other typical teleosts. The overall base composition of the L-strand was found to be T (19.1%), C (31.5%), A (34.8%) and G (14.6%). The control region was 915 nt long and without any repetitive region. The mitogenome sequence data would be useful for studying phylogenetic relationship of C. marulius with other perciform species.