Using a novel library of 5637 expressed sequence tags (ESTs) from the brain tissue of the Asian seabass (Lates calcarifer), we first characterized the brain transcriptome for this economically important species. The ESTs generated from the brain of L. calcarifer yielded 2410 unique transcripts (UTs) which comprise of 982 consensi and 1428 singletons. Based on database similarity, 1005 UTs (41.7%) can be assigned putative functions and were grouped into 12 functional categories related to the brain function. Amongst others, we have identified genes that are putatively involved in energy metabolism, ion pumps and channels, synapse related genes, neurotransmitter and its receptors, stress induced genes and hormone related genes. Subsequently we selected a putative preprocGnRH-II precursor for further characterization. The complete cDNA sequence of the gene obtained was found to code for an 85-amino acid polypeptide that significantly matched preprocGnRH-II precursor sequences from other vertebrates, and possesses structural characteristics that are similar to that of other species, consisting of a signal peptide (23 residues), a GnRH decapeptide (10 residues), an amidation/proteolytic-processing signal (glycine-lysine-argine) and a GnRH associated peptide (GAP) (49 residues). Phylogenetic analysis showed that this putative L. calcarifer preprocGnRH-II sequence is a member of the subcohort Euteleostei and divergent from the sequences of the subcohort Otocephalan. These findings provide compelling evidence that the putative L. calcarifer preprocGnRH-II precursor obtained in this study is orthologous to that of other vertebrates. The functional prediction of this preprocGnRH-II precursor sequence through in silico analyses emphasizes the effectiveness of the EST approach in gene identification in L. calcarifer.
There is very little information on the capacity of freshwater carnivorous fish to biosynthesize highly unsaturated fatty acids (HUFA). The striped snakehead fish (Channa striata) is a carnivorous species cultured inland of several Southeast Asian countries due to its pharmaceutical properties in wound healing enhancement. We described here the full-length cDNA cloning of a striped snakehead fatty acid desaturase (fads), which is responsible for desaturation of unsaturated fatty acids in the HUFA biosynthesis. Bioinformatics analysis reveals a protein coding region with length of 445 amino acids containing all characteristic features of desaturase enzyme, including a cytochrome b5-domain with the heme-binding motif, two transmembrane domains and three histidine-rich regions. The striped snakehead fads amino acid sequence shares high similarity with known fads of other teleosts. The mRNA expression of striped snakehead fads also showed an ontogenic-related increase in expression in 0-20 days after hatch larva. Using ISH, we localized the presence of fads in larva brain, liver and intestinal tissues.
Mackerel (Scombridae; Rastrelliger) are small commercially important pelagic fish found in tropical regions. They serve as a cheap source of animal protein and are commonly used as live bait. By using a truss morphometrics protocol and RAPD analysis, we examined morphological and genetic variation among 77 individual mackerel that were caught using long lines and gillnets at 11 locations along the west coast of Peninsular Malaysia. Nineteen morphometric traits were evaluated and genetic information was estimated using five 10-base RAPD random primers. Total DNA was extracted from muscle tissue. Morphometric discriminant function analysis revealed that two morphologically distinct groups of Rastrelliger kanagurta and a single group of R. brachysoma can be found along the west coast of Peninsular Malaysia. We also found that the head-related characters and those from the anterior part of the body of Rastrelliger spp significantly contribute to stock assessment of this population. RAPD analysis showed a trend similar to that of the morphometric analysis, suggesting a genetic component to the observed phenotypic differentiation. These data will be useful for developing conservation strategies for these species.
Haruan (Channa striatus) is in great demand in the Malaysian domestic fish market. In the present study, mtDNA cyt b was used to investigate genetic variation of C. striatus among populations in Peninsular Malaysia. The overall population of C. striatus demonstrated a high level of haplotype diversity (h) and a low-to-moderate level of nucleotide diversity (π). Analysis of molecular variance (AMOVA) results showed a significantly different genetic differentiation among 6 populations (F(ST)=0.37566, P=0.01). Gene flow (Nm) was high and ranged from 0.32469 to infinity (∞). No significant relationship between genetic distance and geographic distance was detected. A UPGMA tree based on the distance matrix of net interpopulation nucleotide divergence (d(A)) and haplotype network of mtDNA cyt b revealed that C. striatus is divided into 2 major clades. The neutrality and mismatch distribution tests for all populations suggested that C. striatus in the study areas had undergone population expansion. The estimated time of population expansion in the mtDNA cyt b of C. striatus populations occurred 0.72-6.19 million years ago. Genetic diversity of mtDNA cyt b and population structure among Haruan populations in Peninsular Malaysia will be useful in fisheries management for standardization for Good Agriculture Practices (GAP) in fish-farming technology, as well as providing the basis for Good Manufacturing Practices (GMP).
Cryptocaryoniasis (also known as marine white spot disease) is mediated by Cryptocaryon irritans. This obligate ectoparasitic protozoan infects virtually all marine teleosts, which includes Lates calcarifer, a highly valuable aquaculture species. Little is known about L. calcarifer-C. irritans interactions. This study was undertaken to gain an informative snapshot of the L. calcarifer transcriptomic response over the course of C. irritans infection. An in-house fabricated cDNA microarray slides containing 3872 probes from L. calcarifer liver and spleen cDNA libraries were used as a tool to investigate the response of L. calcarifer to C. irritans infection. Juvenile fish were infected with parasites for four days, and total RNA was extracted from liver tissue, which was harvested daily. We compared the transcriptomes of C. irritans-infected liver to uninfected liver over an infection period of four days; the comparison was used to identify the genes with altered expression levels in response to C. irritans infection. The greatest number of infection-modulated genes was recorded at 2 and 3 days post-infection. These genes were mainly associated with the immune response and were associated in particular with the acute phase response. Acute phase proteins such as hepcidin, C-type lectin and serum amyloid A are among the highly modulated genes. Our results indicate that an induced acute phase response in L. calcarifer toward C. irritans infection is similar to the responses observed in bacterial infections of teleosts. This response demonstrates the importance of first line defenses in teleost innate immune responses against ectoparasite infection.
DNA barcodes, typically focusing on the cytochrome oxidase I gene (COI) in many animals, have been used widely as a species-identification tool. The ability of DNA barcoding to distinguish species from a range of taxa and to reveal cryptic species has been well documented. Despite the wealth of DNA barcode data for fish from many temperate regions, there are relatively few available from the Southeast Asian region. Here, we target the marine fish Family Carangidae, one of the most commercially-important families from the Indo-Malay Archipelago (IMA), to produce an initial reference DNA barcode library.
Cryptocaryon irritans causes Cyptocaryonosis or white spot disease in a wide range of marine fish including Lates calcarifer (Asian seabass). However, the immune response of this fish to the parasite is still poorly understood. In this study, quantitative polymerase chain reaction (qPCR) was performed to assess the expression profile of immune-related genes in L. calcarifer infected by C. irritans. A total of 21 immune-related genes encoding various functions in the fish immune system were utilized for the qPCR analysis. The experiment was initiated with the infection of juvenile fish by exposure to theronts from 200 C. irritans cysts, and non-infected juvenile fish were used as controls. Spleen, liver, gills and kidney tissues were harvested at three days post-infection from control and infected fish. In addition, organs were also harvested on day-10 post-infection from fish that had been allowed to recover from day-4 up to day-10 post-infection. L. calcarifer exhibited pathological changes on day-3 post-infection with the characteristic presence of white spots on the entire fish body, excessive mucus production and formation of a flap over the fish eye. High quality total RNA was extracted from all tissues and qPCR was performed. The qPCR analysis on the cohort of 21 immune-related genes of the various organs harvested on day-3 post-infection demonstrated that most genes were induced significantly (p
This study provides a first description of the phylogeographic patterns and evolutionary history of two species of the mudskipper genus Periophthalmus. These amphibious gobies are distributed throughout the whole Indo-Pacific region and Atlantic coast of Africa, in peritidal habitats of soft-bottom coastal ecosystems. Three sequence datasets of two widely distributed species, Periophthalmus argentilineatus and P. kalolo, were obtained by amplifying and sequencing two mtDNA markers (D-loop and 16S rDNA) and the nDNA rag1 region. The three datasets were then used to perform phylogeographic, demographic and population genetic analyses. Our results indicate that tectonic events and past climatic oscillations strongly contributed to shape present genetic differentiation, phylogeographic and demographic patterns. We found support for the monophyly of P. kalolo, and only shallow genetic differentiation between East-African and Indo-Malayan populations of this species. However, our collections of the morphospecies P. argentilineatus include three molecularly distinct lineages, one of them more closely related to P. kalolo. The presence of Miocenic timings for the most recent common ancestors of some of these morphologically similar clades, suggests the presence of strong stabilising selection in mudskippers' habitats. At population level, demographic analyses and palaeoecological records of mangrove ecosystems suggest that Pleistocene bottlenecks and expansion plus secondary contact events of the studied species were associated with recurrent sea transgressions during interglacials, and sea regressions or stable regimes during glacials, respectively.
CXCR3 is a CXC chemokine receptor 3 which binds to CXC ligand 4 (CXCL4), 9, 10 and 11. CXC chemokine receptor 3a (CXCR3a) is one of the splice variants of CXCR3. It plays crucial role in defense and other physiological processes. In this study, we report the molecular cloning, characterization and gene expression of CXCR3a from striped murrel Channa striatus (Cs). The full length CsCXCR3a cDNA sequence was obtained from the constructed cDNA library of striped murrel by cloning and sequencing using an internal sequencing primer. The full length sequence is 1425 nucleotides in length including an open reading frame of 1086 nucleotides which is encoded with a polypeptide of 361 amino acids (mol. wt. 40 kDa). CsCXCR3a domain analysis showed that the protein contains a G protein coupled receptor between 55 and 305 along with its family signature at 129-145. The transmembrane prediction analysis showed that CsCXCR3a protein contains 7 transmembrane helical regions at 34-65, 80-106, 113-146, 154-181, 208-242, 249-278 and 284-308. The 'DRY' motif from CsCXCR3a protein sequence at (140)Asp-(141)Arg-(142)Tyr which is responsible for G-protein binding is also highly conserved with CXCR3 from other species. Phylogenetic tree showed that the CXC chemokine receptors 3, 4, 5 and 6, each formed a separate clade, but 1 and 2 were clustered together, which may be due to the high similarity between these receptors. The predicted 3D structure revealed cysteine residues, which are responsible for 'CXC' motif at 116 and 198. The CsCXR3a transcript was found to be high in kidney, further its expression was up-regulated by sodium nitrite acute toxicity exposure, fungal, bacterial and poly I:C challenges. Overall, these results supported the active involvement of CsCXCR3a in inflammatory process of striped murrel during infection. However, further study is necessary to explore the striped murrel chemokine signaling pathways and their roles in defense system.
In this study, we reported a molecular characterization of a novel proto-type galectin-1 from the striped murrel Channa striatus (named as CsGal-1). The full length CsGal-1 was identified from an established striped murrel cDNA library and further we confirmed the sequence by cloning. The complete cDNA sequence of CsGal-1 is 590 base pairs (bp) in length and its coding region encoded a poly peptide of 135 amino acids. The polypeptide contains a galactoside binding lectin domain at 4-135. The domain carries a sugar binding site at 45-74 along with its signatures (H(45)-X-Asn(47)-X-Arg(49) and Trp(69)-X-X-Glu(72)-X-Arg(74)). CsGal-1 shares a highly conserved carbohydrate recognition domain (CRD) with galectin-1 from other proto-type galectin of teleosts. The mRNA expressions of CsGal-1 in healthy and various immune stimulants including Aphanomyces invadans, Aeromonas hydrophila, Escherchia coli lipopolysaccharide and poly I:C injected tissues of C. striatus were examined using qRT-PCR. CsGal-1 mRNA is highly expressed in kidney and is up-regulated with different immune stimulants at various time points. To understand its biological activity, the coding region of CsGal-1 gene was expressed in an E. coli BL21 (DE3) cloning system and its recombinant protein was purified. The recombinant CsGal-1 protein was agglutinated with mouse erythrocytes at a concentration of 4μg/mL in a calcium independent manner. CsGal-1 activity was inhibited by d-galactose at 25mM(-1) and d-glucose and d-fructose at 100mM(-1). The results of microbial binding assay showed that the recombinant CsGal-1 protein agglutinated only with the Gram-negative bacteria. Interestingly, we observed no agglutination against Gram-positive bacteria. Overall, the study showed that CsGal-1 is an important immune gene involved in the recognition and elimination of pathogens in C. striatus.
The endogenous production of long-chain polyunsaturated fatty acids (LC-PUFA) in carnivorous teleost species inhabiting freshwater environments is poorly understood. Although a predatory lifestyle could potentially supply sufficient LC-PUFA to satisfy the requirements of these species, the nutrient-poor characteristics of the freshwater food web could impede this advantage. In this study, we report the cloning and functional characterisation of an elongase enzyme in the LC-PUFA biosynthesis pathway from striped snakehead (Channa striata), which is a strict freshwater piscivore that shows high deposition of LC-PUFA in its flesh. We also functionally characterised a previously isolated fatty acyl desaturase cDNA from this species. Results showed that the striped snakehead desaturase is capable of Δ4 and Δ5 desaturation activities, while the elongase showed the characteristics of Elovl5 elongases. Collectively, these findings reveal that striped snakehead exhibits the genetic resources to synthesise docosahexaenoic acid (DHA; 22:6n-3) from eicosapentaenoic acid (EPA; 20:5n-3). Both genes are expressed at considerable levels in the brain and the liver. In liver, both genes were up-regulated by dietary C18 PUFA, although this increase did not correspond to a significant rise in the deposition of muscle LC-PUFA. Brain tissue of fish fed with plant oil diets showed higher expression of fads2 gene compared to fish fed with fish oil-based diet, which could ensure DHA levels remain constant under limited dietary DHA intake. This suggests the importance of DHA production from EPA via the ∆4 desaturation step in order to maintain an optimal reserve of DHA in the neuronal tissues of carnivores.
Taxonomic confusion exists within the genus Epinephelus due to the lack of morphological specializations and the overwhelming number of species reported in several studies. The homogenous nature of the morphology has created confusion in the Malaysian Marine fish species Epinephelus fuscoguttatus and Epinephelus hexagonatus. In this study, the partial DNA sequence of the 16S gene and mitochondrial nucleotide sequences of two gene regions, Cytochrome Oxidase Subunit I and III were used to investigate the phylogenetic relationship between them. In the phylogenetic trees, E. fuscoguttatus was monophyletic with E. hexagonatus species and morphology examination shows that no significant differences were found in the morphometric features between these two taxa. This suggests that E. fuscoguttatus is not distinguishable from E. hexagonatus species, and that E. fuscoguttatus have been identified to be E. hexagonatus species is likely attributed to differences in environment and ability to camouflage themselves under certain conditions. Interestingly, this finding was also supported by Principal Component Analysis on Attenuated Total Reflectance-Fourier-transform Infrared (ATR-FTIR) data analysis. Molecular, morphological and meristic characteristics were combined with ATR-FTIR analysis used in this study offer new perspectives in fish species identification. To our knowledge, this is the first report of an extensive genetic population study of E. fuscoguttatus in Malaysia and this understanding will play an important role in informing genetic stock-specific strategies for the management and conservation of this highly valued fish.
Phylogeographic patterns and population structure of the pelagic Indian mackerel, Rastrelliger kanagurta were examined in 23 populations collected from the Indonesian-Malaysian Archipelago (IMA) and the West Indian Ocean (WIO). Despite the vast expanse of the IMA and neighbouring seas, no evidence for geographical structure was evident. An indication that R. kanagurta populations across this region are essentially panmictic. This study also revealed that historical isolation was insufficient for R. kanagurta to attain migration drift equilibrium. Two distinct subpopulations were detected between the WIO and the IMA (and adjacent populations); interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the IMA and WIO regions suggest historical isolation as a result of fluctuations in sea levels during the late Pleistocene. This occurrence resulted in the evolution of a phylogeographic break for this species to the north of the Andaman Sea.
The mitogenome of the Australian freshwater blackfish, Gadopsis marmoratus was recovered coverage by genome skimming using the MiSeq sequencer (GenBank Accession Number: NC_024436). The blackfish mitogenome has 16,407 base pairs made up of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs, and a 819 bp non-coding AT-rich region. This is the 5th mitogenome sequence to be reported for the family Percichthyidae.
Genetic variability and differences in wild striped snakehead Channa striata from Malaysia were analysed by genotyping nine novel nuclear microsatellite loci. Analysis revealed moderate-to-high genetic diversity in most of the populations, indicative of large effective population sizes. The highly diversified populations are admixed populations and, therefore, can be recommended as potential candidates for selective breeding and conservation since they each contain most of the alleles found in their particular region. Three homogenous groups of the wild populations were identified, apparently separated by effective barriers, in accordance with contemporary drainage patterns. The highest population pairwise FST found between members of the same group reflects the ancient population connectivity; yet prolonged geographical isolation resulted in adaptation of alleles to local contemporary environmental change. A significant relationship between genetic distance and geographical isolation was observed (r = 0·644, P < 0·01). Anthropogenic perturbations indicated apparent genetic proximity between distant populations.
Giant grouper (Epinephelus lanceolatus) is a commercially important species, but its wild population has recently been classified as vulnerable. This species has significant potential for use in aquaculture, though a greater understanding of population genetics is necessary for selective breeding programs to minimize kinship for genetically healthy individuals. High-throughput pyrosequencing of genomic DNA was used to identify and characterize novel tetra- and trinucleotide microsatellite markers in giant grouper from Sabah, Malaysia. In total, of 62,763 sequences containing simple sequence repeats (SSRs) were obtained, and 78 SSR loci were selected to possibly contain tetra- and trinucleotide repeats. Of these loci, 16 had tetra- and 8 had trinucleotide repeats, all of which exhibited polymorphisms within easily genotyped regions. A total of 143 alleles were identified with an average of 5.94 alleles per locus, with mean observed and expected heterozygosities of 0.648 and 0.620, respectively. Among of them, 15 microsatellite markers were identified without null alleles and with Hardy-Weinberg equilibrium. These alleles showed a combined non-exclusion probability of 0.01138. The probability of individual identification (PID) value combined with in descending order 12 microsatellite markers was 0.00008, which strongly suggests that the use of the microsatellite markers developed in this study in various combinations would result in a high resolution method for parentage analysis and individual identification. These markers could be used to establish a broodstock management program for giant grouper and to provide a foundation for genetic studies such as population structure, parentage analysis, and kinship selection.
Despite the potential of vegetable oils as aquafeed ingredients, a major drawback associated with their utilization is the inferior level of beneficial n-3 long-chain polyunsaturated fatty acids (LC-PUFA). Echium oil (EO), which is rich in stearidonic acid (SDA, 18:4n-3), could potentially improve the deposition of n-3 LC-PUFA as the biosynthesis of LC-PUFA is enhanced through bypassing the rate-limiting ∆6 desaturation step. We report for the first time an attempt to investigate whether the presence of a desaturase (Fads2) capable of ∆4 desaturation activities and an elongase (Elovl5) will leverage the provision of dietary SDA to produce a higher rate of LC-PUFA bioconversion. Experimental diets were designed containing fish oil (FO), EO or linseed oil (LO) (100FO, 100EO, 100LO), and diets which comprised equal mixtures of the designated oils (50EOFO and 50EOLO) were evaluated in a 12-week feeding trial involving striped snakeheads (Channa striata). There was no significant difference in growth and feed conversion efficiency. The hepatic fatty acid composition and higher expression of fads2 and elovl5 genes in fish fed EO-based diets indicate the utilization of dietary SDA for LC-PUFA biosynthesis. Collectively, this resulted in a higher deposition of muscle eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) compared to LO-based diets. Dietary EO improved the ratio of n-3 LC-PUFA to n-6 LC-PUFA in fish muscle, which is desirable for human populations with excessive consumption of n-6 PUFA. This study validates the contribution of SDA in improving the content of n-3 LC-PUFA and the ratio of EPA to arachidonic acid (ARA, 20:4n-6) in a freshwater carnivorous species.
The mitogenome of an Australian sample of the mudskipper, Periophthalmus minutus, was recovered from partial sequencing using the MiSeq sequencer. This mudskipper has a mitogenome of 16,506 base pairs (55% A + T content) made up of two ribosomal subunit genes, 13 protein-coding genes, 22 transfer RNAs, and a 838 bp non-coding AT-rich region. This is the first sequenced mitogenome for the genus Periophthalmus and the fifth for the subfamily Oxudercinae.
In this study, the complete mitogenome sequence of the Clarion angelfish, Holacanthus clarionensis (Perciformes: Pomacanthidae) has been sequenced by next-generation sequencing method. The length of the assembled mitogenome is 16,615 bp, including 13 protein coding genes, 22 transfer RNAs, and two ribosomal RNAs genes. The overall base composition of Clarion angelfish is 28.3% for A, 29.3% for C, 16.5% for G, 25.9% for T and show 85% identities to flame angelfish Centropyge loriculus. The complete mitogenome of the Clarion angelfish provides essential and important DNA molecular data for further phylogeography and evolutionary analysis for marine angelfish phylogeny.
In this study, the complete mitogenome sequence of the Blue-face angelfish, Pomacanthus xanthometapon (Perciformes: Pomacanthidae) has been sequenced by the next-generation sequencing method. The assembled mitogenome consisting of 16,533 bp includes 13 protein coding genes, 22 transfer RNAs, and two ribosomal RNAs genes. The overall base composition of Blue-face angelfish is 28.7% for A, 28.9% for C, 15.9% for G, 26.6% for T and show 84% identities to flame angelfish Centropyge loriculus. The complete mitogenome of the Blue-face angelfish provides essential and important DNA molecular data for further phylogeography and evolutionary analysis for marine angelfish phylogeny.