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.
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.
One of the most iconic Australian fish is the Murray cod, Maccullochella peelii (Mitchell 1838), a freshwater species that can grow to ∼1.8 metres in length and live to age ≥48 years. The Murray cod is of a conservation concern as a result of strong population contractions, but it is also popular for recreational fishing and is of growing aquaculture interest. In this study, we report the whole genome sequence of the Murray cod to support ongoing population genetics, conservation, and management research, as well as to better understand the evolutionary ecology and history of the species. A draft Murray cod genome of 633 Mbp (N50 = 109 974bp; BUSCO and CEGMA completeness of 94.2% and 91.9%, respectively) with an estimated 148 Mbp of putative repetitive sequences was assembled from the combined sequencing data of 2 fish individuals with an identical maternal lineage; 47.2 Gb of Illumina HiSeq data and 804 Mb of Nanopore data were generated from the first individual while 23.2 Gb of Illumina MiSeq data were generated from the second individual. The inclusion of Nanopore reads for scaffolding followed by subsequent gap-closing using Illumina data led to a 29% reduction in the number of scaffolds and a 55% and 54% increase in the scaffold and contig N50, respectively. We also report the first transcriptome of Murray cod that was subsequently used to annotate the Murray cod genome, leading to the identification of 26 539 protein-coding genes. We present the whole genome of the Murray cod and anticipate this will be a catalyst for a range of genetic, genomic, and phylogenetic studies of the Murray cod and more generally other fish species of the Percichthydae family.
Management of wild fisheries resources requires accurate knowledge on which species are being routinely exploited, but it can be hard to identify fishes to species level, especially in speciose fish groups where colour patterns vary with age. Snappers of the genus Lutjanus represent one such group, where fishes can be hard to identify and as a result fisheries statistics fail to capture species-level taxonomic information. This study employs traditional morphological and DNA barcoding approaches to identify adult and juvenile Lutjanus species harvested in Malaysian waters. Our results reveal a suite of species that differs markedly from those that have previously been considered important in the Malaysian wild-capture fishery and show that official fisheries statistics do not relate to exploitation at the species level. Furthermore, DNA barcoding uncovered two divergent groups of bigeye snapper ('Lutjanus lutjanus') distributed on either side of the Malay Peninsula, displaying a biogeographical pattern similar to distributions observed for many co-occurring reef-distributed fish groups. One of these bigeye snapper groups almost certainly represents an unrecognized species in need of taxonomic description. The study demonstrates the utility of DNA barcoding in uncovering overlooked diversity and for assessing species catch composition in a complicated but economically important taxonomic group.
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.
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.
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.
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.
The complete mitochondrial genome of the conservationally significant Macquarie perch (Macquaria australasica) was obtained from low-coverage shotgun sequencing using the MiSeq sequencer. The M. australasica mitogenome has 16,496 base pairs (55% A + T content) 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 first mitogenome sequence for the genus Macquaria, and the third to be reported for the family Percichthyidae.
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.
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.
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.
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.
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.
There is a lack of understanding on how the environment and trophic niche affect the capability of long-chain polyunsaturated fatty acids (LC-PUFA) in freshwater carnivorous teleost. In this present study, we isolated and functionally characterised a fatty acyl desaturase (Fads) from the striped snakehead Channa striata. Sequence comparison and phylogenetic analysis suggested a Fads2 protein that is closely related to previously characterised Fads2 proteins from freshwater carnivorous and marine herbivorous fish species. We further demonstrated the capacity of Δ6 and Δ5 desaturation activities for this particular desaturase, with highest activities towards the conversion of omega-3 (n-3) polyunsaturated fatty acids (PUFA). Low Δ4 desaturation activity was also detected, although the significance of this at a physiological level remains to be studied. The expression of this striped snakehead Δ6/Δ5 fads2 gene was highest in brain, followed by liver and intestine. In liver, diet fortified with high LC-PUFA concentration impeded the expression of Δ6/Δ5 fads2 gene compared to vegetable oil (VO) based diets. The discovery of Δ6/Δ5 Fads2 desaturase here complements the previous discovery of a Δ4 Fads2 desaturase and an Elovl5 elongase, lending proof to the existence of all the required enzymatic machinery to biosynthesise LC-PUFA from C18 PUFA in a freshwater carnivorous species.
A broad-scale comparative phylogeographic and phylogenetic study of pennah croakers, mainly Pennahia anea, P. macrocephalus, and P. ovata was conducted to elucidate the mechanisms that may have driven the diversification of marine organisms in Southeast Asian waters. A total of 316 individuals from the three species, and an additional eight and six individuals of P. argentata and P. pawak were employed in this study. Two genetically divergent lineages each of P. argentata and P. anea (lineages L1 and L2) were respectively detected from the analyses based on mitochondrial cytochrome b gene data. Historical biogeography analysis with a multi-gene dataset revealed that Pennahia species most likely originated in the South China Sea and expanded into the eastern Indian Ocean, East China Sea, and northwestern Pacific Ocean through three separate range expansions. The main diversifications of Pennahia species occurred during Miocene and Pliocene periods, and the occurrences of lineage divergences within P. anea and P. argentata were during the Pleistocene, likely as a consequence of cyclical glaciations. The population expansions that occurred after the sea level rise might be the reason for the population homogeneity observed in P. macrocephalus and most P. anea L2 South China Sea populations. The structure observed between the two populations of P. ovata, and the restricted distributions of P. anea lineage L1 and P. ovata in the eastern Indian Ocean, might have been hampered by the northward flowing ocean current at the Malacca Strait and by the distribution of coral reefs or rocky bottoms. While our results support S. Ekman's center-of-origin hypothesis taking place in the South China Sea, the Malacca Strait serving as the center of overlap is a supplementary postulation for explaining the present-day high diversity of pennah croakers centered in these waters.
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.
We examine genetic structuring in three commercially important species of the teleost family Carangidae from Malaysian waters: yellowtail scad Atule mate, bigeye scad Selar crumenophthalmus and yellowstripe scad Selaroides leptolepis, from the Indo-Malay Archipelago. In view of their distribution across contrasting habitats, we tested the hypothesis that pelagic species display less genetic divergence compared with demersal species, due to their potential to undertake long-distance migrations in oceanic waters. To evaluate population genetic structure, we sequenced two mitochondrial (mt)DNA [650 bp of cytochrome oxidase I (coI), 450 bp of control region (CR)] and one nuclear gene (910 bp of rag1) in each species. One hundred and eighty samples from four geographical regions within the Indo-Malay Archipelago including a population of yellowtail from Kuwait were examined. Findings revealed that the extent of genetic structuring among populations in the semi-pelagic and pelagic, yellowtail and bigeye were lower than demersal yellowstripe, consistent with the hypothesis that pelagic species display less genetic divergence compared with demersal species. The yellowtail phylogeny identified three distinct clades with bootstrap values of 86%-99% in mtDNA and 63%-67% in rag1. However, in bigeye, three clades were also observed from mtDNA data while only one clade was identified in rag1 dataset. In yellowstripe, the mtDNA tree was split into three closely related clades and two clades in rag1 tree with bootstraps value of 73%-99% and 56% respectively. However, no geographic structure appears in both mtDNA and rag1 datasets. Hierarchical molecular variance analysis (AMOVA), pair wise FST comparisons and the nearest-neighbour statistic (Snn ) showed significant genetic differences among Kuwait and Indo-Malay yellowtail. Within the Indo-Malay Archipelago itself, two distinct mitochondrial lineages were detected in yellowtail suggesting potential cryptic species. Findings suggests varying degrees of genetic structuring, key information relevant to management of exploited stocks, though more rapidly evolving genetic markers should be used in future to better delimit the nature and dynamics of putative stock boundaries.
Benthic species, though ecologically important, are vulnerable to genetic loss and population size reduction due to impacts from fishing trawls. An assessment of genetic diversity and population structure is therefore needed to assist in a resource management program. To address this issue, the two-spined yellowtail stargazer (Uranoscopus cognatus) was collected within selected locations in the Indo-West Pacific (IWP). The partial mitochondrial DNA cytochrome c oxidase subunit 1 and the nuclear DNA recombination activating gene 1 were sequenced. Genetic diversity analyses revealed that the populations were moderately to highly diversified (haplotype diversity, H = 0.490-0.900, nucleotide diversity, π = 0.0010-0.0034) except sampling station (ST) 1 and 14. The low diversity level, however was apparent only in the matrilineal marker (H = 0.118-0.216; π = 0.0004-0.0008), possibly due to stochastic factors or anthropogenic stressors. Population structure analyses revealed a retention of ancestral polymorphism that was likely due to incomplete lineage sorting in U. cognatus, and prolonged vicariance by the Indo-Pacific Barrier has partitioned them into separate stock units. Population segregation was also shown by the phenotypic divergence in allopatric populations, regarding the premaxillary protrusion, which is possibly associated with the mechanism for upper jaw movement in biomechanical feeding approaches. The moderate genetic diversity estimated for each region, in addition to past population expansion events, indicated that U. cognatus within the IWP was still healthy and abundant (except in ST1 and 14), and two stock units were identified to be subjected to a specific resource management program.
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