Transcriptome analysis revealed hub genes for muscle growth in Indian major carp, Catla catla (Hamilton, 1822).

Catla (Catla catla) is the fastest growing Indian major carp species and forms an important component of the freshwater aquaculture systems in the Indian sub-continent. The molecular mechanisms of growth of the species in response to seasonal water temperature variations hitherto are still unknown. In the current study, high-throughput transcriptome sequencing was used to study the differential gene expression pattern in catla muscle tissues between pre-winter and post-winter fingerling groups and fast-growing table size fish. Transcriptome analysis identified 1677 differentially expressed genes (DEGs) in three different growth stages and 236 common DEGs between fingerling at low temperature and table fish post-winter, including four genes under GH/IGF1 axis and 163 genes under signature for compensatory muscle growth. Molecular pathways for the mapped genes identified 42 KEGG pathways and the critical pathways under Environmental Information Processing identified were PI3K-Akt signaling, AMPK signaling pathway, Calcium signaling pathway and MAPK signaling pathway. In this study, 14 differentially expressed potential regulatory hub genes for growth were identified, for the first time and categorized into three major GO groups: unfolded protein binding, rNA processing and biogenesis and muscle development and differentiation. These regulatory hub genes, except acta1, were found to be upregulated in fast-growing table size and post-winter fingerling groups. The results provided valuable information about the key genes, with potential to be used as biomarkers of growth in breeding programs and contributed to our understanding of the molecular mechanisms and pathways regulating muscle growth, in response to temperature fluctuations and different growth stages in C. catla.

Discovery of alternatively spliced isoforms and long non-coding RNA in full length brain transcriptomes of anadromous Hilsa shad, Tenualosa ilisha (Hamilton, 1822).

BACKGROUND: Full length transcriptomes, achieved through long-read sequencing, along with the isoform analysis can reveal complexities in the gene expression profiles, as well as annotate the transcriptomes of non-model organisms. METHODS AND RESULT: Full length transcripts of brain transcriptome of Tenualosa ilisha, Hilsa shad, were generated through PacBio single molecule real-time sequencing and were characterized. A total of 8.30 Gb clean reads were generated, with PacBio RSII, which resulted in 57,651 high quality consensus transcripts. After removing redundant reads, a total of 19,220 high-quality non-redundant transcripts and 17,341 full length ORF transcripts were classified to 7522 putative ortholog groups. Genes involved in various neural pathways were identified. In addition, isoform clusters and lncRNAs were discovered, along with Hilsa specific transcripts with coding frames and 29,147 SSRs in 944 transcripts (1141 annotated). CONCLUSION: The present study provided, for the first time, a comprehensive view of the alternative isoforms of genes and transcriptome complexity in Hilsa shad brain and forms a rich resource for functional studies in brain of this anadromous fish.

Emergence of carp edema virus in cultured ornamental koi carp, Cyprinus carpio koi, in India.

A disease outbreak was reported in adult koi, Cyprinus carpio koi, from a fish farm in Kerala, India, during June 2015. The clinical signs were observed only in recently introduced adult koi, and an existing population of fish did not show any clinical signs or mortality. Microscopic examination of wet mounts from the gills of affected koi revealed minor infestation of Dactylogyrus sp. in a few koi. In bacteriological studies, only opportunistic bacteria were isolated from the gills of affected fish. The histopathological examination of the affected fish revealed necrotic changes in gills and, importantly, virus particles were demonstrated in cytoplasm of gill epithelial cells in transmission electron microscopy. The tissue samples from affected koi were negative for common viruses reported from koi viz. cyprinid herpesvirus 3, spring viraemia of carp virus, koi ranavirus and red sea bream iridovirus in PCR screening. However, gill tissue from affected koi carp was positive for carp edema virus (CEV) in the first step of nested PCR, and sequencing of PCR amplicons confirmed infection with CEV. No cytopathic effect was observed in six fish cell lines following inoculation of filtered tissue homogenate prepared from gills of affected fish. In bioassay, the symptoms could be reproduced by inoculation of naive koi with filtrate from gill tissue homogenate of CEV-positive fish. Subsequently, screening of koi showing clinical signs similar to koi sleepy disease from different locations revealed that CEV infection was widespread. To our knowledge, this is the first report of infection with CEV in koi from India.

De novo development and characterization of polymorphic microsatellite markers in a schilbid catfish, Silonia silondia (Hamilton, 1822) and their validation for population genetic studies.

The stock characterization of wild populations of Silonia silondia is important for its scientific management. At present, the information on genetic parameters of S. silondia is very limited. The species-specific microsatellite markers were developed in current study. The validated markers were used to genotype individuals from four distant rivers. To develop de novo microsatellite loci, an enriched genomic library was constructed for S. silondia using affinity-capture approach. The markers were validated for utility in population genetics. A total number of 76 individuals from four natural riverine populations were used to generate data for population analysis. The screening of isolated repeat sequences yielded eleven novel polymorphic microsatellite loci. The microsatellite loci exhibited high level of polymorphism, with 6-24 alleles per locus and the PIC value ranged from 0.604 to 0.927. The observed (Ho) and expected (He) heterozygosities ranged from 0.081 to 0.84 and 0.66 to 0.938, respectively. The AMOVA analysis indicated significant genetic differentiation among riverine populations (overall FST = 0.075; P < 0.0001) with maximum variation (92.5%) within populations. Cross-priming assessment revealed successful amplification (35-38 %) of heterologous loci in four related species viz. Clupisoma garua, C. taakree, Ailia coila and Eutropiichthys vacha. The results demonstrated that these de novo polymorphic microsatellite loci are promising for population genetic variation and diversity studies in S. silondia. Cross-priming results indicated that these primers can help to get polymorphic microsatellite loci in the related catfish species of family Schilbidae.

A new fish cell line derived from the caudal fin of freshwater angelfish Pterophyllum scalare: development and characterization.

In this study, a new cell line derived from the caudal fin of the freshwater angelfish Pterophyllum scalare was developed and characterized. The cell line was designated angelfish fin (AFF) and subcultured 44 times since its development. These cells grew well in Leibovitz's -15 medium supplemented with 10% foetal bovine saline (FBS) at 28° C and the modal chromosome number (2n) was 48. The AFF cell-line is mainly comprised of epithelial cells as confirmed by immunocytological technique using anti-cytokeratin antibodies, an epithelial cell marker. This cell line was tested for growth in a temperatures range from 20 to 37° C and at various FBS concentrations of 5-20% at 28° C. The cell line was cryopreserved at different passage levels and revived successfully with 80% survival rate. Polymerase chain reaction amplification and sequencing of partial mitochondrial 16s rRNA and coI genes confirmed that the AFF cell-line originated from angelfish. Mycoplasma sp. contamination was not detected in AFF cells and checked by Hoechst 33258 fluorescence staining. At the 42nd passage the cells were transfected with 2 µg of pAcGFP1-N1 expression vector. The AFF cells exhibited cytotoxic effects when exposed to the bacterial extra cellular products from Serratia marcescens and Proteus hauseri. The AFF cells and cells from kidney and brain did not show cytopathic effect when exposed to cyprinid herpes virus2 and viral nervous necrosis virus. The newly developed AFF cell line will be useful for the isolation of viruses affecting angelfishes, such as iridoviruses, in the future.

Molecular phylogenetics of three species of the genus Rastrelliger using mitochondrial DNA markers.

In the present study three species of mackerel, Rastrelliger present in Indian waters were taken for genetic identification using cytochrome c oxidase subunit I and 16s rRNA sequences. Accurate identification of these species is important for fishery management as its morphological characters are very similar. In this study, the sequences of COI and 16S rRNA were determined from 19 individuals of three Rastrelliger species, Rastrelliger kanagurta, Rastrelliger brachysoma and Rastrelliger faughni from Andamans and Indian mainland to study the phylogenetic relationship. The intraspecies and interspecies genetic distance ranged from 0.000 to 0.002 and 0.007 to 0.015 respectively based on 16S rRNA sequences. Using COI data analysis, the intraspecies genetic distance ranged from 0.000 to 0.012, while it varied from 0.039 to 0.086 for interspecies. The present study clearly demarcates three species of mackerel based on the mitochondrial genetic sequences and also showed a non-descriptive genetic distance of R. kanagurta from mainland and Andaman Islands.

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.

Mitochondrial ATPase 6/8 genes reveal genetic divergence in the Coilia dussumieri (Valenciennes, 1848) populations of north east and northwest coasts of India.

The golden anchovy, Coilia dussumieri, though possessing discontinuous distribution along northeast and northwest coasts of India, it is being managed as unit stock for fishery assessment purposes. By considering the need for stock specific management of the species, mitochondrial ATP synthase 6 and 8 (ATPase 6/8) genes were analyzed for delineating genetic stock structure of the species. Sequence analysis revealed a total of 34 haplotypes across four populations from both the east and west coasts of India. Haplotype diversity (h) was found in the range of 0.7421-0.9368. Similarly, nucleotide diversity (π) varied from 0.0012 to 0.0025. AMOVA results indicated a high total variance of 72.66% between east and west coast populations and less (1.34%) among populations within the respective coast. Phylogenetic tree constructed using pair wise FST also indicated the genetic divergence of populations of east and west coasts of India. The findings of the present study will be helpful in developing stock specific management measures for conservation and sustainable utilization of the species.

Microsatellite markers to determine population genetic structure in the golden anchovy, Coilia dussumieri.

Coilia dussumieri (Valenciennes, 1848) commonly called as golden anchovy, constitutes a considerable fishery in the northern part of both the west and east coasts of India. Despite its clear-cut geographic isolation, the species is treated as a unit stock for fishery management purposes. We evaluated 32 microsatellite primer pairs from three closely related species (resource species) belonging to the family Engraulidae through cross-species amplification in C. dussumieri. Successful cross-priming was obtained with 10 loci, which were sequenced for confirmation of repeats. Loci were tested for delineating the genetic stock structure of four populations of C. dussumieri from both the coasts of India. The number of alleles per locus ranged from 8 to 18, with a mean of 12.3. Results of pairwise F ST indicated genetic stock structuring between the east and west coast populations of India and also validated the utilization of identified microsatellite markers in population genetic structure analysis.

Plectranthias alcocki, a new anthiine fish species (Perciformes: Serranidae) from the Arabian Sea, off southwest India.

A new species of anthiine fish, Plectranthias alcocki n. sp. is described and illustrated based on two specimens, (63.7-72.5 mm SL), recently collected from deep-waters of the Arabian Sea, off Kollam, Kerala, India. The following combination of characters distinguishes it from all other congeners: Dorsal-fin rays X, 15; anal-fin rays III, 7; pectoral-fin rays 14, all unbranched; pelvic-fin rays I, 5; lateral-line complete, the pored lateral-line scales 28; scales above lateral line to origin of dorsal fin 1; scales dorsally on head extending to posterior nostrils; no scales on maxilla or chin; gill rakers 5 + 11 (2 + 7 developed); circumpeduncular scales 10; fourth dorsal spine longest, 2.8 (2.6) in head length (HL), longest dorsal-fin soft ray (second)  2.4 (2.7) in head length; body depth 34.4 (35)% SL; head length 46 (49.8)% SL; orbital length 8.6 in SL; margin of preopercle finely serrate, the serrae 33 (28), ventral edge without antrorse spines; dorsal fin continuous and notched; first anal-fin spine 4.9 (5.6) in HL, second anal-fin spine 2.2 (2.6) in HL; pelvic fins relatively short, 4.0-4.3 in SL; the dorsal fin with a black blotch at base of fourth to eighth spines, one at base of the last three spines, and two at base of soft portion of fin, the dark pigment extending onto adjacent body.

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.

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.

The sequence and de novo assembly of the genome of the Indian oil sardine, Sardinella longiceps.

The Indian oil sardine, Sardinella longiceps, is a widely distributed and commercially important small pelagic fish of the Northern Indian Ocean. The genome of the Indian oil sardine has been characterized using Illumina and Nanopore platforms. The assembly is 1.077 Gb (31.86 Mb Scaffold N50) in size with a repeat content of 23.24%. The BUSCO (Benchmarking Universal Single Copy Orthologues) completeness of the assembly is 93.5% when compared with Actinopterygii (ray finned fishes) data set. A total of 46316 protein coding genes were predicted. Sardinella longiceps is nutritionally rich with high levels of omega-3 polyunsaturated fatty acids (PUFA). The core genes for omega-3 PUFA biosynthesis, such as Elovl 1a and 1b,Elovl 2, Elovl 4a and 4b,Elovl 8a and 8b,and Fads 2, were observed in Sardinella longiceps. The presence of these genes may indicate the PUFA biosynthetic capability of Indian oil sardine, which needs to be confirmed functionally.

Complete mitochondrial genome sequence of Catla catla and its phylogenetic consideration.

Complete nucleotide sequence of mitochondrial genome (mitogenome) of the Catla catla (Ostariophysi: Cypriniformes: Cyprinidae) was determined in the present study. Its length is 16,594 bp and contains 13 protein coding genes, 22 transfer RNAs, two ribosomal RNAs and one non-coding control region. Most of the genes were encoded on the H-strand, while the ND6 and eight tRNA (Gln, Ala, Asn, Cys, Tyr, Ser (UCN), Glu and Pro) genes were encoded on the L-strand. The reading frames of two pair of genes overlapped: ATPase 8 with 6 and ND4L with ND4 by seven nucleotides each. The main non-coding region was 929 bp, with three conserved sequence blocks (CSB-I, CSB-II, and CSB-III) and an unusual simple sequence repeat, (TA)(7). Phylogenetic analyses based on complete mitochondrial genome sequences were in favor of the traditional taxonomy of family Cyprinidae. In conclusion present mitogenome of Catla catla adds more information to our understanding of diversity and evolution of mitogenome in fishes.

Development of an ES-like cell culture system (RESC) from rohu, Labeo rohita (Ham.).

An embryonic stem (ES)-like cell culture system RESC from a commercially important freshwater carp, Labeo rohita, was developed using blastula stage embryos. The cells were cultured in Leibovitz-15 (L-15) medium in gelatin-coated cell culture flask supplemented with 15 % fetal bovine serum along with 10 ng ml(-1) basic fibroblast growth factor at 28 °C under feeder-free conditions. The ES-like cells were characterized by their unique morphology, alkaline phosphatase activity, embryoid body formation tendency, expression of transcription factor Oct4, and consistent chromosome count. The RESC cells when treated with retinoic acid differentiated into cells of different lineages. The RESC developed from mid-blastula embryos of L. rohita would be a useful tool for cellular differentiation and gene expression studies.

Ontogeny of the digestive tract in butter catfish Ompok bimaculatus (Bloch) larvae.

The ontogeny of the digestive tract was studied histologically in butter catfish Ompok bimaculatus from hatching to 30 days post-hatching (dph). At hatching, the digestive tract of butter catfish consisted of a straight tube with a smooth lumen dorsally attached to the yolk sac. Between 1 and 2 dph, the mouth opened, oral valves were visible and canine-like teeth and taste buds were detected. During this period, intestine was differentiated into the anterior and posterior intestine, and the digestive accessory glands were also developed. Exogenous feeding started at 2 dph, and there was a 2-day mixed endogenous-exogenous feeding period. Most of the yolk sac reserves were consumed between 2 and 3 dph, and by 5 dph, the yolk sac was completely depleted and no longer visible in histological sections. Between 3 and 4 dph, several vacuoles (neutral lipids) were observed in the intestine and also in hepatocytes, indicating a functional absorption of nutrients from food. At 8 dph, differentiation of gastric glands was noticed, and by 9-11 dph, there were abundant gastric tubular glands arranged along numerous longitudinal folds. During the same period, pyloric sphincter appeared as an epithelial fold that separated the stomach from the anterior intestine. From 12 dph to the end of the study at 30 dph, no noticeable histological modifications were observed. The development of gastric glands is considered as the last major events in digestive tract development and their presence designates the end of larval period and the onset of the juvenile period. Hence, it is suggested that, butter catfish larvae have a morphologically complete digestive tract by 12 dph. These findings on the development of the digestive system in butter catfish may lead to a better understanding of the ontogeny and would be useful to improve the larval rearing techniques of this promising catfish species for freshwater aquaculture diversification.

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.

Muscle transcriptome resource for growth, lipid metabolism and immune system in Hilsa shad, Tenualosa ilisha.

The information on the genes involved in muscle growth, lipid metabolism and immune systems would help to understand the mechanisms during the spawning migration in Hilsa shad, which in turn would be useful in its future domestication process. The primary objective of this study was to generate the transcriptome profile of its muscle through RNA seq. The total RNA was isolated and library was prepared from muscle tissue of Tenualosa ilisha, which was collected from Padma River at Farakka, India. The prepared library was then sequenced by Illumina HiSeq platform, HiSeq 2000, using paired-end strategy. A total of 8.68 GB of pair-end reads of muscle transcriptome was generated, and 43,384,267 pair-end reads were assembled into 3,04,233 contigs, of which 23.99% of assembled contigs has length ≥ 150 bp. The total GO terms were categorised into cellular component, molecular function and biological process through PANTHER database. Fifty-three genes related to muscle growth were identified and genes in different pathways were: 75 in PI3/AKT, 46 in mTOR, 76 in MAPK signalling, 24 in Janus kinase-signal transducer and activator of transcription, 45 in AMPK and 27 in cGMP pathways. This study also mined the genes involved in lipid metabolism, in which glycerophospholipid metabolism contained highest number of genes (32) and four were found to be involved in fatty acid biosynthesis. There were 58 immune related genes found, in which 31 were under innate and 27 under adaptive immunity. The present study included a large genomic resource of T. ilisha muscle generated through RNAseq, which revealed the essential dataset for our understanding of regulatory processes, specifically during the seasonal spawning migration. As Hilsa is a slow growing fish, the genes identified for muscle growth provided the basic information to study myogenesis. In addition, genes identified for lipid metabolism and immune system would provide resources for lipid synthesis and understanding of Hilsa defense mechanisms, respectively.

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 N50 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.