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.

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.

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.

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.

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.

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.

Population structure and genetic diversity of Indian Major Carp, Labeo rohita (Hamilton, 1822) from three phylo-geographically isolated riverine ecosystems of India as revealed by mtDNA cytochrome b region sequences.

The population structure and genetic diversity of Rohu (Labeo rohita Hamilton, 1822) was studied by analysis of the partial sequences of mitochondrial DNA cytochrome b region. We examined 133 samples collected from six locations in three geographically isolated rivers of India. Analysis of 11 haplotypes showed low haplotype diversity (0.00150), nucleotide diversity (π) (0.02884) and low heterogeneity value (0.00374). Analysis of molecular variance (AMOVA) revealed the genetic diversity of L. rohita within population is very high than between the populations. The Fst scores (-0.07479 to 0.07022) were the indication of low genetic structure of L. rohita populations of three rivers of India. Conspicuously, Farakka-Bharuch population pair Fst score of 0.0000, although the sampling sites are from different rivers. The phylogenetic reconstruction of unique haplotypes revealed sharing of a single central haplotype (Hap_1) by all the six populations with a point mutations ranging from 1-25 nucleotides.

Mitochondrial signatures for identification of grouper species from Indian waters.

Groupers are important commercial fish in many parts of the world. Accurate identification is critical for effective conservation assessment and fisheries management. Genetic barcodes provide a simple and reproducible method for the identification of species even in the absence of taxonomic expertise. The generation of reference barcodes from properly identified specimens is an important first step in this direction. Here, 36 species belonging to the subfamily Epinephelinae (Family: Serranidae) were collected from landings on the west coast of India and Port Blair, Andaman, and partial nucleotide sequence data of the mitochondrial cytochrome C oxidase subunit I (COI) gene was generated. Barcodes for 13 species were developed from Indian waters for the first time. Analysis using the COI gene produced phylogenetic trees in concurrence with other multi-gene studies. Epinephelus fasciatus and E. areolatus were found to be a species complex, as hypothesized in other studies. The DNA barcodes developed in the study can be used for identifying species within Epinehelinae, where taxonomic ambiguity still exists.

DNA barcoding reveals species composition of sharks and rays in the Indian commercial fishery.

DNA barcoding was successfully used for the accurate identification of chondrichthyans in the Indian commercial marine fishery. About 528 specimens of 111 chondrichthyan species and 34 families, collected from the Indian EEZ, were barcoded for a 655 bp region of the mitochondrial gene cytochrome c oxidase subunit 1 (COI). Generally, five specimens per species were barcoded, but numbers ranged from 2 to 13. The average Kimura 2 parameter (K2P) distance separating individuals within species was 0.32%, and the average distance separating species within genera was 6.73%. Ten species were suggested as putative new species requiring formal descriptions. Based on the morphology and molecular support, 11 elasmobranch species were confirmed first records for Indian waters. The present study confirms the ability of DNA barcoding for the accurate identification of sharks, rays, and their products from Indian waters.

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

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

Complete mitochondrial genome sequence of Indian medium carp, Labeo gonius (Hamilton, 1822) and its comparison with other related carp species.

In the present study, the complete mitochondrial genome sequence of Labeo gonius is reported using PGM sequencer (Ion Torrent). The complete mitogenome of L. gonius is obtained by the de novo sequences assembly of genomic reads using the Torrent Mapping Alignment Program (TMAP) which is 16 614 bp in length. The mitogenome of L. gonius comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNA genes, and D-loop as control region along with gene order and organization, being similar to most of other fish mitogenomes of NCBI databases. The mitogenome in the present study has 99% similarity to the complete mitogenome sequence of Labeo fimbriatus, as reported earlier. The phylogenetic analysis of Cypriniformes depicted that their mitogenomes are closely related to each other. The complete mitogenome sequence of L. gonius would be helpful in understanding the population genetics, phylogenetics, and evolution of Indian Carps.

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.

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.

In Silico Mining of Conserved miRNAs of Indian Catfish Clarias batrachus (Linnaeus, 1758) from the Contigs, ESTs, and BAC End Sequences.

MicroRNAs are small non-coding RNAs which play significant role in RNA interference. The present work deals with mining of the conserved miRNA and their target genes from the contigs, ESTs, and BAC end sequences of commercially important catfish, Clarias batrachus, from India. A total of 138, 1 and 1 conserved pre-miRNA sequences, were mined from the contigs, ESTs, and BAC end sequences, respectively. The analysis of families of the conserved pre-miRNA revealed conservation of the fish-specific family mir-430 and other important families, such as mir-455, let-7, mir-133, and mir-137. The mir-455 is involved in hypoxia signaling, let-7 family represents potential anti-tumor molecules involved in human cancer therapy, whereas mir-133 and mir-137 have high therapeutic potentials. Using an alternate computational in silico approach, mining of mature miRNAs resulted in identification of 210 mature miRNAs from contigs, 1 from EST, and 2 each from forward as well as reverse BAC end sequences. Target prediction of these putative miRNAs resulted in the identification of 66,758 and 18,747 target genes in C. batrachus and Danio rerio, respectively. Functional annotation of these miRNAs indicated their involvement in diverse biological functions. The findings of the present study can serve as a valuable resource for further functional genomics studies in C. batrachus.

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.

Development and characterization of genic SSR markers from low depth genome sequence of Clarias batrachus (magur).

Indian magur (Clarias batrachus) is an important freshwater catfish, which is listed as endangered under A3cde+4acde ver. 3.1 categories by the IUCN (2015) due to decreasing population trend. Microsatellites or short sequence repeats (SSRs) tagged to genes have been utilized as gene marker. In the present study, 31,814 SSRs of C. batrachus (magur) were identified using microsatellite identification tool programme from the next generation sequencing data generated on Roche 454 and Ion Torrent platforms. A bioinformatics pipeline, with stringent criteria resulted in selection of 1672 microsatellite loci falling in the genic region. Initially, a total of 30 loci were selected for primer development; and of these 14 were successfully amplified and five were found to be polymorphic in 30 individuals of C. batrachus (magur). The observed as well as expected heterozygosity ranged from 0.038 to 0.526 and 0.434 to 0.784, respectively, and the number of observed alleles ranged from three to five. The study reported the application of next generation sequencing technologies for rapid development of microsatellite loci in Indian catfish species, C. batrachus (magur).

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.

Detection of goldfish haematopoietic necrosis herpes virus (Cyprinid herpesvirus-2) with multi-drug resistant Aeromonas hydrophila infection in goldfish: First evidence of any viral disease outbreak in ornamental freshwater aquaculture farms in India.

This outbreak report details of a mortality event where Cyprinid herpes virus-2 (CyHV-2) was detected in association with multidrug-resistant Aeromonas hydrophila infection in goldfish, Carassius auratus, from commercial farms. The goldfish exhibited large scale haemorrhages on the body, fins and gills, lepidorthosis, necrosed gills, protruded anus and shrunken eyes. White nodular necrotic foci in spleen and kidneys were noticed, along with necrosis and fusion of gill lamellae. Transmission electron microscopy of affected tissues revealed the presence of mature virus particles. Involvement of CyHV-2 was confirmed by PCR, sequencing and observed cytopathic effect in koi carp fin cell line along with experimental infection study. A bacterium isolated from the internal organs of affected fish was found to be pathogenic Aeromonas hydrophila having resistance to more than 10 classes of antibiotics. We postulate that CyHV-2 was the primary etiological agent responsible for this outbreak with secondary infection by A. hydrophila. The experimental infection trials in Labeo rohita and koi carp by intraperitoneal challenge with CyHV-2 tissue homogenates failed to reproduce the disease in those co-cultured fish species. This is the first report of a viral disease outbreak in organised earthen ornamental fish farms in India and bears further investigation.