This study was established to test the hypothesis of whether the codon optimization of fish growth hormone gene (FGH) based on P. pastoris preferred codon will improve the quantity of secreted rFGH in culture supernatant that can directly be used as fish feed supplements. The optimized FGH coding sequence (oFGH) and native sequence (nFGH) of giant grouper fish (Epinephelus lanceolatus) were cloned into P. pastoris expression vector (pPICZαA) downstream of alcohol oxidase gene (AOX1) for efficient induction of extracellular rFGH by adding 1% of absolute methanol. The results showed that recombinant P. pastoris was able to produce 2.80 ± 0.27 mg of oFGH compared to 1.75 ± 0.25 of nFGH in one litre of culture supernatant. The total body weight of tiger grouper fingerlings fed with oFGH increased significantly at third (P < 0.05) and fourth weeks (P < 0.01) of four-week experiment period compared to those fed with nFGH. Both oFGH and nFGH significantly enhanced the final biomass and fish survival percentage. In conclusion, codon optimization of FGH fragment was useful to increase rFGH quantity in the culture supernatant of P. pastoris that can be directly used as fish feed supplements. Further studies are still required for large scale production of rFGH and practical application in aquaculture production.
Morphological identification of fish taxa can sometimes prove difficult because phenotypic variation is either being affected by environmental factors, phenotypic characters are highly conserved or marker selection has been inappropriate. DNA based markers especially neutral mitochondrial DNA (mtDNA) have been used widely in recent times to provide better resolution of systematic relationships among vertebrate taxa. The Asian Arowana (Scleropages formosus) is a high value ornamental fish belonging to the family Osteoglossidae with a number of different colour variants distributed geographically across different locations around Southeast Asia. Systematic relationships among colour variants still remain unresolved. Partial sequences of the Cytochrome B (Cyt B) and DNA barcoding gene, Cytochrome C Oxidase I (COI) were used here to assess genetic relationships among colour variants and as a tool for molecular identification for differentiating among colour variants in this species. Results of the study show that in general, colour pattern shows no relationship with extent of COI or Cyt B mtDNA differentiation and so cannot be used to identify taxa. Partial sequences of the mtDNA genes were sufficient however, to identify S. formosus from a closely related species within the order Osteoglossidae.
A phylogeographic study of an economically important freshwater fish, the striped snakehead, Channa striata in Sundaland was carried out using data from mtDNA ND5 gene target to elucidate genetic patterning. Templates obtained from a total of 280 individuals representing 24 sampling sites revealed 27 putative haplotypes. Three distinct genetic lineages were apparent; 1)northwest Peninsular Malaysia, 2)southern Peninsular, east Peninsular, Sumatra and SW (western Sarawak) and 3) central west Peninsular and Malaysian Borneo (except SW). Genetic structuring between lineages showed a significant signature of natural geographical barriers that have been acting as effective dividers between these populations. However, genetic propinquity between the SW and southern Peninsular and east Peninsular Malaysia populations was taken as evidence of ancient river connectivity between these regions during the Pleistocene epoch. Alternatively, close genetic relationship between central west Peninsular Malaysia and Malaysian Borneo populations implied anthropogenic activities. Further, haplotype sharing between the east Peninsular Malaysia and Sumatra populations revealed extraordinary migration ability of C. striata (>500 km) through ancient connectivity. These results provide interesting insights into the historical and contemporary landscape arrangement in shaping genetic patterns of freshwater species in Sundaland.
Climatic differences across a taxon's range may be associated with specific bioenergetic demands and may result in genetics-based metabolic adaptation, particularly in aquatic ectothermic organisms that rely on heat exchange with the environment to regulate key physiological processes. Extending down the east coast of Australia, the Great Dividing Range (GDR) has a strong influence on climate and the evolutionary history of freshwater fish species. Despite the GDR acting as a strong contemporary barrier to fish movement, many species, and species with shared ancestries, are found on both sides of the GDR, indicative of historical dispersal events. We sequenced complete mitogenomes from the four extant species of the freshwater cod genus Maccullochella, two of which occur on the semi-arid, inland side of the GDR, and two on the mesic coastal side. We constructed a dated phylogeny and explored the relative influences of purifying and positive selection in the evolution of mitogenome divergence among species. Results supported mid- to late-Pleistocene divergence of Maccullochella across the GDR (220-710 thousand years ago), bringing forward previously reported dates. Against a background of pervasive purifying selection, we detected potentially functionally relevant fixed amino acid differences across the GDR. Although many amino acid differences between inland and coastal species may have become fixed under relaxed purifying selection in coastal environments rather than positive selection, there was evidence of episodic positive selection acting on specific codons in the Mary River coastal lineage, which has consistently experienced the warmest and least extreme climate in the genus.
The Japanese scad Decapterus maruadsi (Carangidae) is an economically important marine species in Asia but its exploitation shows signs of overfishing. To document its stock structure, a population genetic and phylogeographic study of several populations of this species from the central part of the Indo-West Pacific region was conducted using the mitochondrial cytochrome b gene. Genetic homogeneity within the Sundaland region's population, including Rosario (the Philippines) and Ranong (Andaman Sea) populations was revealed with low nucleotide diversity (π = 0.001-0.003) but high haplotype diversity (h = 0.503-0.822). In contrast, a clear genetic structure was observed between this group and the northern Vietnam populations as revealed by FST, AMOVA and SAMOVA, while the central Vietnam population of Khanh Hoa is an admixed group between the two differentiated regional populations. The neutrality and mismatch distribution analyses supported a demographic expansion of D. maruadsi in between last Pleistocene to early Holocene period which influenced present day distribution pattern. Contemporary factors such as oceanic currents and different life history traits are also believed to play significant roles in the observed population structure and biogeographical pattern. Based on these results, recommendations on how stocks of the Japanese scad should be managed are offered.
In this study, we first conducted a genome survey assay for Sillago sihama by Illumina sequencing platform, and then developed 15 polymorphic microsatellite loci in a wild population. A total of 129.46 Gb raw data were obtained, of which 115.07 Gb were clean data, with a sequencing depth of 179.3-folds. This genome was estimated to be 522.6 Mb in size, with the heterozygosity, repeat content and GC content being 0.63%, 21% and 44%. A total of 630,028 microsatellites were identified from the genome, of which, dinucleotide repeat was the most abundant (56.80%), followed by mononucleotide repeat (30.23%). Furthermore, 60 pairs of primers were designed and synthesized based on microsatellite sequences, of which 15 were polymorphic in a wild population. A total of 91 alleles were found, with an average of 6.07 per locus. Number of alleles, observed and expected heterozygosity per locus ranged from two to 13, from 0.250 to 0.862, and from 0.396 to 0.901, respectively. Twelve loci were highly informative (PIC > 0.5), and the others were medium informative (0.25
Invasive alien fish species have become a silent treat towards the ecosystem especially the native fish population in Malaysia. There has been a need to develop rapid identification methods that can aid management teams in identifying fish species that are not native to our ecosystem. Current visual identification methods are highly tedious and require time, delaying action towards curbing the invasion. The LAMP assay successfully identified six popular invasive fish species in Malaysia. None of the LAMP assays showed false positives and the Limit of Detection of the LAMP primers were highly sensitive and could detect DNA samples up to 1 × 10-15 ng/μl. The LAMP primers designed were highly specific to the target species and did not amplify non target species. DNA sequencing was done to ensure the accuracy of LAMP assay results. This study demonstrates that LAMP is a suitable tool in species identification efforts of invasive fish species in Malaysia.
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.
Mislabelling in fish products is a highly significant emerging issue in world fish trade in terms of health and economic concerns. DNA barcoding is an efficient sequencing-based tool for detecting fish species substitution but due to DNA degradation, it is in many cases difficult to amplify PCR products of the full-length barcode marker (~650 bp), especially in severely processed products. In the present study, a pair of universal primers targeting a 198 bp sequence of the mitochondrial 16s rRNA gene was designed for identification of fish species in the processed fish products commonly consumed in Malaysia. The specificity of the universal primers was tested by both in-silico studies using bioinformatics software and through cross-reaction assessment by practical PCR experiments against the DNA from 38 fish species and 22 other non-target species (animals and plants) and found to be specific for all the tested fish species. To eliminate the possibility of any false-negative detection, eukaryotic endogenous control was used during specificity evaluation. The developed primer set was validated with various heat-treated (boiled, autoclaved and microwaved) fish samples and was found to show high stability under all processing conditions. The newly developed marker successfully identified 92% of the tested commercial fish products with 96-100% sequence similarities. This study reveals a considerable degree of species mislabelling (20.8%); 5 out of 24 fish products were found to be mislabelled. The new marker developed in this work is a reliable tool to identify fish species even in highly processed products and might be useful in detecting fish species substitution thus protecting consumers' health and economic interests.
The migration of anadromous fish in heterogenic environments unceasingly imposes a selective pressure that results in genetic variation for local adaptation. However, discrimination of anadromous fish populations by fine-scale local adaptation is challenging because of their high rate of gene flow, highly connected divergent population, and large population size. Recent advances in next-generation sequencing (NGS) have expanded the prospects of defining the weakly structured population of anadromous fish. Therefore, we used NGS-based restriction site-associated DNA (NextRAD) techniques on 300 individuals of an anadromous Hilsa shad (Tenualosa ilisha) species, collected from nine strategic habitats, across their diverse migratory habitats, which include sea, estuary, and different freshwater rivers. The NextRAD technique successfully identified 15,453 single nucleotide polymorphism (SNP) loci. Outlier tests using the FST OutFLANK and pcadapt approaches identified 74 and 449 SNPs (49 SNPs being common), respectively, as putative adaptive loci under a divergent selection process. Our results, based on the different cluster analyses of these putatively adaptive loci, suggested that local adaptation has divided the Hilsa shad population into two genetically structured clusters, in which marine and estuarine collection sites were dominated by individuals of one genetic cluster and different riverine collection sites were dominated by individuals of another genetic cluster. The phylogenetic analysis revealed that all the riverine populations of Hilsa shad were further subdivided into the north-western riverine (turbid freshwater) and the north-eastern riverine (clear freshwater) ecotypes. Among all of the putatively adaptive loci, only 36 loci were observed to be in the coding region, and the encoded genes might be associated with important biological functions related to the local adaptation of Hilsa shad. In summary, our study provides both neutral and adaptive contexts for the observed genetic divergence of Hilsa shad and, consequently, resolves the previous inconclusive findings on their population genetic structure across their diverse migratory habitats. Moreover, the study has clearly demonstrated that NextRAD sequencing is an innovative approach to explore how dispersal and local adaptation can shape genetic divergence of non-model anadromous fish that intersect diverse migratory habitats during their life-history stages.
Snakehead fishes of the family Channidae are predatory freshwater teleosts from Africa and Asia comprising 38 valid species. Snakeheads are important food fishes (aquaculture, live food trade) and have been introduced widely with several species becoming highly invasive. A channid barcode library was recently assembled by Serrao and co-workers to better detect and identify potential and established invasive snakehead species outside their native range. Comparing our own recent phylogenetic results of this taxonomically confusing group with those previously reported revealed several inconsistencies that prompted us to expand and improve on previous studies. By generating 343 novel snakehead coxI sequences and combining them with an additional 434 coxI sequences from GenBank we highlight several problems with previous efforts towards the assembly of a snakehead reference barcode library. We found that 16.3% of the channid coxI sequences deposited in GenBank are based on misidentifications. With the inclusion of our own data we were, however, able to solve these cases of perpetuated taxonomic confusion. Different species delimitation approaches we employed (BIN, GMYC, and PTP) were congruent in suggesting a potentially much higher species diversity within snakeheads than currently recognized. In total, 90 BINs were recovered and within a total of 15 currently recognized species multiple BINs were identified. This higher species diversity is mostly due to either the incorporation of undescribed, narrow range, endemics from the Eastern Himalaya biodiversity hotspot or the incorporation of several widespread species characterized by deep genetic splits between geographically well-defined lineages. In the latter case, over-lumping in the past has deflated the actual species numbers. Further integrative approaches are clearly needed for providing a better taxonomic understanding of snakehead diversity, new species descriptions and taxonomic revisions of the group.
Adaptive differences across species' ranges can have important implications for population persistence and conservation management decisions. Despite advances in genomic technologies, detecting adaptive variation in natural populations remains challenging. Key challenges in gene-environment association studies involve distinguishing the effects of drift from those of selection and identifying subtle signatures of polygenic adaptation. We used paired-end restriction site-associated DNA sequencing data (6,605 biallelic single nucleotide polymorphisms; SNPs) to examine population structure and test for signatures of adaptation across the geographic range of an iconic Australian endemic freshwater fish species, the Murray cod Maccullochella peelii. Two univariate gene-association methods identified 61 genomic regions associated with climate variation. We also tested for subtle signatures of polygenic adaptation using a multivariate method (redundancy analysis; RDA). The RDA analysis suggested that climate (temperature- and precipitation-related variables) and geography had similar magnitudes of effect in shaping the distribution of SNP genotypes across the sampled range of Murray cod. Although there was poor agreement among the candidate SNPs identified by the univariate methods, the top 5% of SNPs contributing to significant RDA axes included 67% of the SNPs identified by univariate methods. We discuss the potential implications of our findings for the management of Murray cod and other species generally, particularly in relation to informing conservation actions such as translocations to improve evolutionary resilience of natural populations. Our results highlight the value of using a combination of different approaches, including polygenic methods, when testing for signatures of adaptation in landscape genomic studies.
The Merbok Estuary comprises one of the largest remaining mangrove forests in Peninsular Malaysia. Its value is significant as it provides important services to local and global communities. It also offers a unique opportunity to study the structure and functioning of mangrove ecosystems. However, its biodiversity is still partially inventoried, limiting its research value. A recent checklist based on morphological examination, reported 138 fish species residing, frequenting or subject to entering the Merbok Estuary. In this work, we reassessed the fish diversity of the Merbok Estuary by DNA barcoding 350 specimens assignable to 134 species initially identified based on morphology. Our results consistently revealed the presence of 139 Molecular Operational Taxonomic Units (MOTUs). 123 of them are congruent with morphology-based species delimitation (one species = one MOTU). In two cases, two morphological species share the same MOTU (two species = one MOTU), while we unveiled cryptic diversity (i.e. COI-based genetic variability > 2%) within seven other species (one species = two MOTUs), calling for further taxonomic investigations. This study provides a comprehensive core-list of fish taxa in Merbok Estuary, demonstrating the advantages of combining morphological and molecular evidence to describe diverse but still poorly studied tropical fish communities. It also delivers a large DNA reference collection for brackish fishes occurring in this region which will facilitate further biodiversity-oriented research studies and management activities.
The complete mitogenome of the ray Taeniura lymma was recovered from genome skimming using the HiSeq sequencing system. The T. lymma mitogenome has 17,652 base pairs (59.13% A + T content) made up of 13 protein-coding genes, 2 ribosomal subunit genes, 22 transfer RNAs and a 1906 bp non-coding AT-rich region. This mitogenome sequence is the second for a ray from Australian waters, the first for the genus Taeniura and the ninth for the family Dasyatidae.
For centuries, morphology-based fish identification has been applied without molecular evaluation. Many studies showed that specimens with a similar morphology are frequently found to be quite genetically distinct. One of the fish species that still remains taxonomically problematic is a commercial snapper species, Lutjanus johnii. Because of morphological ambiguities among local fish taxonomists in Malaysia, we examined the ability of the cytochrome oxidase I (COI) gene to genetically examine the taxonomic status of L. johnii. A 626-base pair COI region was successfully amplified and aligned with conspecific sequences that were retrieved from GenBank. The phylogenetic tree obtained showed two major clusters; the first cluster consists of L. johnii from Straits of Malacca, Thailand, Australia, and China while the second cluster comprises L. johnii from China and India. The latter group showed sequence divergence greater than 3.5%. After observing this, we suspected that there might be a cryptic species between the South China Sea and Indian Ocean. This is the first molecular report concerning the commercial species of snapper, L. johnii, in Malaysia, which had only gained provisional recognition from morphological examination.
Conservation is imperative for the Asian snakeheads Channa striata, as the species has been overfished due to its high market demand. Using maternal markers (mitochondrial cytochrome c oxidase subunit 1 gene (COI)), we discovered that evolutionary forces that drove population divergence did not show any match between the genetic and morphological divergence pattern. However, there is evidence of incomplete divergence patterns between the Borneo population and the populations from Peninsular Malaysia. This supports the claim of historical coalescence of C. striata during Pleistocene glaciations. Ecological heterogeneity caused high phenotypic variance and was not correlated with genetic variance among the populations. Spatial conservation assessments are required to manage different stock units. Results on DNA barcoding show no evidence of cryptic species in C. striata in Malaysia. The newly obtained sequences add to the database of freshwater fish DNA barcodes and in future will provide information relevant to identification of species.
To satisfy increasing demands for fish as food, progress must occur towards greater aquaculture productivity whilst retaining the wild and farmed genetic resources that underpin global fish production. We review the main selection methods that have been developed for genetic improvement in aquaculture, and discuss their virtues and shortcomings. Examples of the application of mass, cohort, within family, and combined between-family and within-family selection are given. In addition, we review the manner in which fish genetic resources can be lost at the intra-specific, species and ecosystem levels and discuss options to best prevent this. We illustrate that fundamental principles of genetic management are common in the implementation of both selective breeding and conservation programmes, and should be emphasized in capacity development efforts. We highlight the value of applied genetics approaches for increasing aquaculture productivity and the conservation of fish genetic resources.
Malaysian arowana (dragonfish; Scleropages formosus) is an ancient osteoglossid fish from southeast Asia. Due to the high demand of the ornamental fish trade and because of habitat loss, the species is close to extinction. We isolated and characterized 10 polymorphic microsatellites of this species, using 5'-anchored PCR. The number of alleles at the 10 microsatellite loci ranged from 2 to 28, with a mean of 7.8/locus. The observed heterozygosity ranged from 0.03 to 0.93 (mean: 0.39), whereas the expected heterozygosity ranged from 0.03 to 0.94 (mean: 0.46). Seven microsatellites deviated from Hardy-Weinberg equilibrium, and three conformed to Hardy-Weinberg equilibrium and were in linkage equilibrium. These 10 novel microsatellites should facilitate studies of genetic diversity and population structure of arowana to help plan actions for the conservation of the indigenous Malaysian arowana.
The most economically important form of aquaculture is fish farming, which is an industry that accounts for an ever increasing share of world fishery production. Molecular markers can be used to enhance the productivity of the aquaculture and fish industries to meet the increasing demand. Molecular markers can be identified via a DNA test regardless of the developmental stage, age or environmental challenges experienced by the organism. The application of 16s and cytochrome b markers has enabled rapid progress in investigations of genetic variability and inbreeding, parentage assignments, species and strain identification and the construction of high resolution genetic linkage maps for aquaculture fisheries. In this review, the advantages of principles and potential power tools of 16s and cytochrome b markers are discussed. Main findings in term of trend, aspects and debates on the reviewed issue made from the model of aquatic species for the benefit of aquaculture genomics and aquaculture genetics research are discussed. The concepts in this review are illustrated with various research examples and results that relate theory to reality and provide a strong review of the current status of these biotechnology topics.
A total of 143 microsatellites were isolated from Mystus nemurus using a 5' anchored polymerase chain reaction technique or the random amplified hybridization microsatellite method, the first set of microsatellite markers developed for the Southeast Asian river catfish. Twenty polymorphic microsatellite loci were used as markers for population characterization of M. nemurus from six different geographical locations in Malaysia (Perak, Kedah, Johor, UPM, Sarawak and Terengganu). The number of alleles per locus ranged from 2 to 11 with 6.3 as the average number of alleles per locus. Characterization of the populations showed relatively high levels of genetic variation compared with previous studies using allozyme markers. The highest genetic similarity was found between Perak and Kedah, while the highest genetic distance was found between Terengganu and Kedah. The majority of clustering was in accordance with geographical locations and the histories of the populations. Microsatellite analysis indicated that the Sarawak population might be genetically closer to the Peninsular Malaysian populations than has been previously shown by other molecular marker studies.