Short Report

Characterization of the complete mitochondrial genome of Mauritian sardinella, Sardinella jussieu (Lacepède, 1803), collected in the Banten Bay, Indonesia

Sinar Pagi Sektiana1,3, Sapto Andriyono1,4,5, Hyun-Woo Kim1,2,*http://orcid.org/0000-0003-1357-5893
Author Information & Copyright
1Interdisciplinary Program of Biomedical, Mechanical and Electrical EngineeringPukyong National University48513BusanRepublic of Korea
2Department of Marine BiologyPukyong National University48513BusanRepublic of Korea
3Aquaculture Technology Study ProgramJakarta Fisheries UniversityJl. AUP Pasar Minggu Jakarta Selatan12520JakartaIndonesia
4Fisheries and Marine FacultyC Campus Jl. Mulyorejo60115SurabayaIndonesia
5Universitas AirlanggaSurabayaEast JavaIndonesia
*82-51-629-5926kimhw@pknu.ac.kr

© The Author(s) 2017. Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Received: Jul 12, 2017; Accepted: Sep 28, 2017

Published Online: Oct 23, 2017

Abstract

Fishes in genus Sardinella are small pelagic species, which plays an important role in marine ecosystem as the first consumer. Those species are also commercially important, whose total catch reaches 278,600 tons in 2011 in Indonesia, but their identification has been difficult for their morphological similarity. In this study, we reported Sardinella jussieu for the first time in Indonesian coastal area (Banten Bay, Indonesia, 6° 0′ 50.00″ S–106° 10′ 21.00″ E). We were able to confirm the species by both its morphological characteristics including the black spot at dorsal fin origin, the dusky pigmentation at caudal fin, 31 total scute numbers, and DNA sequence identity in the GenBank database by the molecular analysis. Its total mitochondrial genome was determined by the combination of next-generation sequencing and typical PCR strategy. The total mitochondrial genome of Sardinella jussieu (16,695 bp) encoded 13 proteins, 2 ribosomal RNAs, 22 transfer RNAs, and the putative control region. All protein-coding genes started with ATG and typical stop codon and ended with TAA or TAG except for ND4 in which AGA is used. Phylogenetic analyses of both COI region and full mitochondrial genome showed that S. jussieu is most closely related to Sardinella albella and Sardinella gibbosa.

Keywords: Barcode; Mitochondrial genome; Sardinella; Indonesia; Next-generation sequencing

Background

Sardinella is a genus of fish in the family Clupeidae found in the Atlantic, Indian, and the Pacific Ocean. The paddle-shaped supramaxilla bones are major characteristics, which help distinguish Sardinella from other genera. Morphological characters distinguish Sardinella from all other clupeoid genera with the presence of two fleshy outgrowths on the hind margin of the gill opening (Whitehead 1985). According to FishBase (http://www.fishbase.org/), there are currently 22 recognized species in the genus Sardinella. Sardinella is important not only in marine food webs as a base consumer supporting tuna, seabirds, and marine mammals (Willette et al. 2011) but also in industry as the protein source with a low cost using as a bait for large fish or a feed in aquaculture.

Seven species in the genus Sardinella are currently known in Indonesian waters including Sardinella fimbriata, Sardinella gibbosa, Sardinella lemuru, Sardinella albella, Sardinella atricauda, Sardinella branchysoma, and Sardinella melanura, whose total catch in Indonesia reaches 278,600 tons in 2011 (MMAF 2012). Morphological identification in Sardinella is mainly characterized by their gill raker, pelvic scute, scales, and otolith (Homayuni et al. 2013; Bräger and Moritz 2016; Begg and Waldman 1999). However, species identification in the genus Sardinella is often hard for its broad geographical ranges, overlapping distributions (Willette et al. 2011) and morphological similarities (Sivakumaran et al. 1987) especially in larval stages (Ditty et al. 1994), which makes it difficult to manage the Sardinella resources in Indonesia.

In addition to the traditional morphological identification, the genetic information is now alternatively used for the species identification for its fast and exact results. The most widely used genetic markers are partial mitochondrial DNA sequences such as cytochrome C oxidase I (COI) or cytochrome B (CytB) (Palumbi et al. 1991; Ward et al. 2005; Vrijenhoek 1994). However, full mitochondrial genome sequences provide more information about its biogeographical or evolutional information than those fragmental sequences. Therefore, more than 5000 mitochondrial genomes have been deposited in GenBank database (www.ncbi.nlm.nih.gov) from 33,500 species identified based on morphological characteristics (www.fishbase.org).

In this study, we report the Mauritian sardinella, Sardinella jussieu, for the first time in Indonesian coastal waters, which was collected from the Banten Bay. S. jussieu was previously reported only in the Western Indian Ocean, Taiwan, Hong Kong, and Vietnam (www.fishbase.org). Morphological characteristics of Sardinella jussieu are distinguished within other Sardinella species with the presence of black spot at dorsal fin origin and dusky pigmentation at caudal fin, total scute measurement which is 31, and vertical striae on a scale not meeting at center and no perforation on hind part (Whitehead 1985). After confirmation of the species by the molecular COI markers, its total mitochondrial genome sequence was determined by the combination of the traditional PCR methods and next-generation sequencing (NGS) techniques.

Methods

Sample collection and morphological measurement

Five individuals of S. jussieu were collected in the Banten Bay, Indonesia (6° 0′ 50.00″S 106° 10′ 21.00″ E), in January 2016 as the part of the regular fish survey (Fig. 1). Collected fish were directly stored in 96% ethanol and kept at − 20 °C until the further analysis (Knebelsberger and Stöger 2012). Morphological identification was made by their body shape, type of scale, fin feature, morphometric (i.e., standard length, body width, and head length), and meristic characteristic (total number of scutes) (Whitehead 1985; Strauss and Bond 1990).

fas-20-0-26-g1
Fig. 1. Sampling location in Banten Bay, Indonesia (red square)
Download Original Figure

Genomic DNA extraction and next-generation sequencing

Genomic DNA was extracted using an AccuPrep® Genomic DNA Extraction Kit (Bioneer) according to the manufacturer’s instruction. A small portion of tail fin was dissected, which was further homogenized by the TissueLyser II (Qiagen). Purified genomic DNA was quantified with nanoDrop (Thermofisher Scientific D1000), aliquoted, and stored at the − 70 °C for further analysis.

Two universal primer sets targeting cytochrome c oxidase I (COI) region, Fish F1 and Fish R1 (Ward et al. 2005), and targeting cytochrome b (cyt-B) region, GLUDG-L and CB2-H (Palumbi et al. 1991), were used to obtain the partial sequences of each gene, respectively (Table 1). The quality of all the primers used in this experiment was analyzed by the OligoAnalyzer 3.1 (http://sg.idtdna.com/calc/analyzer) and commercially synthesized by Bioneer Co. (Korea). Each PCR mixture (20 μL) contained 12.8 μL ultrapure water, 1 μL primer (0.5 μM, forward and reverse), 0.2 μL Ex Taq DNA polymerase (TaKaRa, Japan), 2 μL 10× Buffer, 2 μL dNTPs (1 μM, TaKaRa, Japan), and 100 ng genomic DNA as template. PCR was carried out under the following condition: initial denaturation step at 95 °C for 3 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 50 °C for 30 s, and extension at 72 °C for 45 s (COI target sequence) or 30 s (Cyt-B target sequence). The process was completed with a final extension at 72 °C for 10 min. Two PCR products targeting partial sequences of COI and Cyt B were then purified with AccuPrep Gel purification kit (Bioneer, Korea) and ligated into a cloning vector (Promega, USA), sequenced in both directions.

Table 1. Primers used for the mitochondrial genome of Sardinella jussieu

Name

Sequence (5' to 3')

Product size (bp)

Fish F1

TCAACCAACCACAAAGACATTGGCAC

652

Fish R1

TAGACTTCTGGGTGGCCAAAGAATCA

GLUDG-L

TGACTTGAARAAYCGTTG

451

GB2-H

CCCTCAGAATGATATTTGTCCTCA

CYTB-F

GCCTACGAAAAACCCACCCGCTCC

8.2 k

CO1-R

GTAAGGTCTACGGATGCCCCTGCG

CYTB-R

AACGGAGGAGAAAGCGGTTGCGATG

8.7 k

CO1-F

CTTCCTGCTTCTCCTGGCCTCCTC

SARD F

TTAAAGTCCTCCCTGAGGCCC

683

SARD R

TTAGGAGGGAGTCGTCAAATGC

Download Excel Table

In order to obtain two large PCR products (~ 8 kb), two pairs of sequence-specific primer sets (CYT-F and CO1-R and CO1-F and CYT B-R) were designed based on the obtained partial sequences of each region (Table 1). Each PCR reaction (30 μL) contained 19.7 μL ultrapure water, 1 μL of each primer (0.5 μM), 0.3 μL Ex Taq Hot Start Version DNA polymerase (TAKARA, Japan), 3 μL 10× Buffer, 3 μL dNTPs (1 mM, Takara, Japan), and 100 ng genomic DNA as template. PCR was carried out with two-step PCR protocol for long PCR under the following condition: initial denaturation step at 94 °C for 3 min, followed by 30 cycles of denaturation at 98 °C for 10s, and annealing and extension at 68 °C for 10 min. The process was completed with a final extension at 72 °C for 10 min. Two large PCR products were pooled together in equal concentration and fragmented to ~ 350 bp in length by Covaris M220 (Covaris Inc.). TruSeq® sample preparation kit V2 (Illumina, USA) was used for the construction of a library from fragmented sequence and quality and quantity of the constructed library was measured using 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Sequencing was performed by Illumina Miseq platform (2 × 300 bp pair ends) (Illumina, USA).

Assembly of mitochondrial genome by the bioinformatic analysis

Raw reads from MiSeq sequencer, with under Qv 20 and more than ambiguous nucleotides, were removed from raw read using CLC Genomic Workbench v 7.5 (CLC BIO Aarhus, Denmark). Mothür software was used to pairing forward and reverse sequence with more than 7 bp overlapped and without any mismatch. Paired sequence then assembled using Geneious R8 with minimum 20 bp of overlapping sequence and 100% overlap identity. Ambiguous sequences of the D-loop region were reconfirmed by the typical end-point PCR and with sequence-specific primers (Sard_F and Sard_R) and DNA sequencing of its PCR products by Sanger sequencing method (Table 1).

Results and discussion

Morphological and molecular identification of Sardinella jussieu

As the result of morphometric measurements, we determined that the collected five fish were S. jussieu. Among the morphologically similar fish Sardinella species including S. albella, S. atricauda, S. fimbriata, S. marquesensis, S. sindensis, and S. gibbosa, the scale and pigmentation patterns are useful characteristics to identify species (Bräger and Moritz 2016; Strauss and Bond 1990). The average ratio of body depth (BD) to standard length (SL) of the collected samples was 27.5%, and total scute numbers were 31 (Table 2). Vertical striae on scales did not meet at center with no perforations on hind part of the scale, and the pigmented dorsal and caudal fins were also identified (Fig. 2). Those morphological characteristics suggested that the collected samples were S. jussieu. The most closely related Sardinella species, S. albella and S. gibbosa, are distinguished from S. jussieu in the presence of scale perforations (Table 3). Molecular identification of five Sardinella samples confirmed the morphological identification. The COI region of five individuals (652 bp) exhibited 100% sequence identity to Sardinella sp. (GenBank accession number: KJ566769) collected from the coastal water in Thailand and 99% to S. jussieu (GenBank accession no.: HQ231358) collected from the Philippines (Quilang et al. 2011). Based on the morphological characteristics and DNA sequence identity, we concluded that five Sardinella samples collected in the Banten Bay, Indonesia, were Mauritian sardinella, Sardinella jussieu.

Table 2. General morphometric and meristic (total scute) of S. jussieu

Sample

Measurement

Standard length/SL (mm)

Body depth/BD (mm)

Head length/HL (mm)

Eye diameter/ED (mm)

BD/SL (%)

HL/SL (%)

ED/SL (%)

HL/ED (%)

Total scute

1

50

13.5

12

3.5

27.0

24.0

7.0

29.2

30

2

47

12.5

11

3

26.6

23.4

6.4

27.3

31

3

45

12.5

11

3

27.8

24.4

6.7

27.3

31

4

41

11.5

11

3

28.0

26.8

7.3

27.3

31

5

52

14.5

13

3.5

27.9

25.0

6.7

26.9

32

Average

47

12.9

11.6

3.2

27.5

24.7

6.8

27.6

31

Download Excel Table
fas-20-0-26-g2
Fig. 2. Mauritian sardinella (S. jussieu) collected from Banten Bay, Indonesia (a). The fish scale of S. jussieu presents no perforations and vertical striated with not meeting at center (b) according to Whitehead (1985) (c). Black scale bar = 1 cm
Download Original Figure
Table 3. Comparison of morphological characteristic of seven Sardinella species

Name

Scale

Fin

Striae connected/overlapped

Perforations

Dark spot at dorsal fin origin

Dark spot at caudal fin

S. fimbriata

S. gibbosa

S. albella

S. atricauda

S. brachysoma

S. melanura

S. jussieu

Download Excel Table

Complete mitochondrial genome of the Sardinella jussieu

In order to have additional information of S. jussieu, the complete mitochondrial genome sequence was determined by the NGS and bioinformatic sequence assembly. Its mitochondrial genome was 16,695 bp in length consisting of 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, and the putative control region (Fig. 3). The base composition was 4415 A (26%), 4132 T (25%), 4900 C (29%), and 3248 G (19%). The purines and pyrimidines are A+T content (51%) slightly higher than G+C content (49%). The highest A+T content was observed in the putative control region (66%), which is similar to the other previous studies. The H strands encode 28 genes while the L strands encode 9 genes (Table 4). Among the protein-coding genes, three overlap nucleotides up to 10 bp, ATP8–ATP6, ND4L–ND4, and ND5–ND6, were detected. The transfer RNA gene pair tRNA−Ile–tRNA−Gln and t RNA−Thr–tRNA−Pro overlaps 1 bp as well. A total of 1292 bp of noncoding nucleotides are apparent in the S. jussieu with 1029 bp at putative control region, and 263 remains spread over 11 intergenic nucleotides; 68.3% (11.397 bp) of total mitochondrial genome sequence encoded 13 proteins and the size of each gene ranged from 168 bp (ATP8) to 1836 bp (ND5). Except for ND6, all protein-coding genes were encoded by H strand (Fig. 3). Although all 13 genes begin with typical start codon, ATG, there were several stop codons including typical ones such as TAA (CO1, COII, ATP8, ATP6, COIII, ND4L, ND5, CYTB), TAG (ND2, ND3, ND6, ND1), and exceptional AGA in ND4 gene (Table 4). Overlapping nucleotides were identified in three pairs of protein-coding genes (10 nucleotides for ATP8 and ATP6, seven for ND4L and ND4, and four for ND5 and ND6).

fas-20-0-26-g3
Fig. 3. Mitochondrial genomic organization of Sardinella jussieu
Download Original Figure
Table 4. Organization of the full-length mitochondrial genome of Sardinella jussieu

Feature

Position numbers

Size (bp)

Strand

Intergenic nucleotides

Codon

Anticodon/position

Start

Stop

tRNA–Phe

1–68

68

H

GAA/31–33

12S rRNA

69–1019

951

H

0

tRNA–Val

1020–1091

72

H

0

TAC/1053–1055

16S rRNA

1092–2777

1686

H

0

tRNA–Leu

2779–2853

75

H

1

TAA/2814–2816

ND1

2854–3828

975

H

0

ATG

TAG

tRNA–Ile

3837–3908

72

H

8

GAT/3869–3871

tRNA–Gln

3908–3978

71

L

− 1

TTG/3944–3946

tRNA–Met

3978–4046

69

H

− 1

CAT/4008–4010

ND2

4020–5093

1074

H

− 27

ATG

TAG

tRNA–Trp

5092–5163

72

H

− 2

TCA/5125–5127

tRNA–Ala

5165–5233

69

L

1

TGC/5173–5175

tRNA–Asn

5236–5308

73

L

2

GTT/5273–5275

tRNA–Cys

5345–5410

66

L

36

GCA/5366–5368

tRNA–Tyr

5414–5484

71

L

3

GTA/5450–5452

COX1

5678–7036

1359

H

193

ATG

TAA

tRNA–Ser

7037–7107

71

L

0

TGA/7073–7075

tRNA–Asp

7112–7180

69

H

4

GTC/7142–7144

COII

7193–7897

705

H

12

ATG

TAA

tRNA–Lys

7884–7957

74

H

− 14

TTT/7918–7920

ATP8

7959–8126

168

H

1

ATG

TAA

ATP6

8117–8800

684

H

− 10

ATG

TAA

COIII

8800–9585

786

H

− 1

ATG

TAA

tRNA–Gly

9585–9656

72

H

0

TCC/9618–9620

ND3

9600–10007

408

H

− 57

ATG

TAG

tRNA–Arg

10006–10075

70

H

− 2

TCG/1037–10039

ND4L

10076–10372

297

H

0

ATG

TAA

ND4

10366–11751

1386

H

− 7

ATG

AGA

tRNA–His

11747–11815

69

H

− 5

GTG/11747–11815

tRNA–Ser

11816–11883

68

H

0

GCT/11842–11844

tRNA–LEu

11884–11955

72

H

0

TAG/11916–11918

ND5

11956–13791

1836

H

0

ATG

TAA

ND6

13788–14309

522

L

− 4

ATG

TAG

tRNA–Glu

14310–14378

69

L

0

TTC/14346–14348

CYTB

14385–15581

1197

H

6

ATG

TAA

tRNA–Thr

15526–15597

72

H

− 56

TGT/15558–15560

tRNA–Pro

15597–15666

70

L

− 1

TGG/15620–15622

Control region

15667–16695

1029

H

0

Download Excel Table

The mitochondrial genome of S. jussieu contained 22 tRNA genes (Fig. 4), which showed the difference in their sizes from 68 bp (tRNA–Phe) to 71 (tRNA–Gln). Fourteen tRNA genes encode in H strand and 8 genes encoded in L strand (Fig. 3). The 12S rRNA gene (951 bp) of S. jussieu was located between the tRNA–Phe and tRNA–Val, whereas 1686 bp of 16S rRNA was between tRNA–Val and tRNA–Leu. Twenty-one tRNA structures were predicted to have typical three arms except for tRNAser, which showed two arms. That result was also identified in the other Sardinella species (Lavoué et al. 2007). The putative control region of S. jussieu (1029 bp) was longest among three other Sardinella species including S. longiceps (958 bp) (GenBank accession number: NC033407), S. albella (986 bp) (GenBank accession number: NC016726), and S. maderensis (986 bp) (GenBank accession number: NC009587).

fas-20-0-26-g4
Fig. 4. Putative secondary structure tRNA genes in mitochondrial genomic. Proposed structure of 22 tRNA genes encoded in the mitochondrial of Sardinella jussieu
Download Original Figure

Total mitochondrial DNA sequence of S. jussieu showed 84–93% identity with those of currently known three other Sardinella species among which S. albella is the most closely related to S. jussieu (Fig. 5a). In order to know the better evolutional relationship of S. jussieu, its COI sequence was compared with those of the other 12 Sardinella species (Fig. 5b). As shown in the analysis by the full mitochondrial genomes, S. jussieu showed the most closely related to S. albella with 96% sequence identity. In fact, DNA sequence identity of two species S. albella and Sardinella gibbosa was too high to be distinct with each other in the COI region (Figs. 5b). Although morphological keys to discriminate two species were proposed, the post-pelvic ventral scutes and gill rakers number on a lower limb (Stern et al. 2016), S. albella and S. gibbosa frequently misidentified as shown in the COI barcodes. From the reason, it is required to compare full-length mitochondrial sequences of two species for the better classification. As the lowest sequence identity to other Sardinella species, control region of S. jussieu mitochondrial genome would be the good candidate to discriminate them.

fas-20-0-26-g5
Fig. 5. a Phylogenetic tree of mitochondrial genome of four species belonging to Sardinella. The phylogenetic tree was constructed using molecular evolutionary genetic analysis ver.6.0 (MEGA 6, MEGA Inc. Englewood, NJ), program with the minimum evolution algorithm, the evolutionary distances were computed using Kimura 2-Parameter method and b Phylogenetic tree of CO1 sequences of 18 species belonging to genus Sardinella. The phylogenetic tree was constructed using molecular evolutionary genetic analysis ver.6.0 (MEGA 6, MEGA Inc. Englewood, NJ), program with the minimum evolution algorithm, the evolutionary distances were computed using Kimura 2-Parameter method
Download Original Figure

In this study, we identified that S. jussieu inhabits in Java island, Indonesia, as well as the two previously known Sardinella species, S. albella and S. gibbosa. Although S. jussieu is originally distributed in the western Indian Ocean from the western coast of southern India from Bombay South to Sri Lanka also Madagascar and Mauritius, the recent information is also caught in Taiwan (Hu et al. 2015), Hong Kong (Leung 1997), and the Philippines (Quilang et al. 2011). The result strongly supported that S. jussieu is more widely distributed than we have thought and the large-scale survey should be made to know the spatiotemporal distribution of four Sardinella species in Indonesia. We here reported the full-length mitochondrial genome sequence of S. jussieu collected from Java island, which would provide the important information for the scientific management of Sardinella species in Indonesia. We expect that more Sardinella species may exist in Java island and more information about the mitochondrial genome of the other unreported Sardinella species such as S. gibbosa would be a useful information for the molecular biological tools to discriminate different Sardinella species in Indonesia.

Conclusion

This study determined the complete mitochondrial DNA (mtDNA) sequence of S. jussieu in Java Island, Indonesia, for the first time. The mtDNA sequence is 16.695 bp in length and comprises the typical set of 2 rRNAs, 22 tRNA genes, 13 protein-coding genes, and putative control region. Mitochondrial genome structure of S. jussieu was identical to those in other Sardinella genus. Phylogenetic analysis using full mitochondrial genome exhibits that S. jussieu were most closely related to S. albella. However, comparison in the COI region showed that relationship between S. albella and S. gibbosa was ambiguous and determination of the complete mitochondrial DNA sequence of S. gibbosa is required for the better understanding of evolutional relationship between S. jussieu and those species. Those information would provide the basic information for the scientific management of Sardinella species in Indonesia.

Acknowledgements

All authors thank the funding of grant from the Pukyong National University, and authors also thank the Jakarta fisheries University students who helped for collecting the fish samples.

Funding

This work was supported by a grant from the Pukyong National University in 2016.

Availability of data and materials

The dataset(s) supporting the conclusions of this article is(are) included within the article.

Authors’ contributions

SPS carried out the whole process, participated in the whole experiment, and drafted the manuscript. SA added and drafted the final phylogenetic analysis. HWK participated in the design of the study and edited the manuscript. All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1.

Begg GA, Waldman JR. An holistic approach to fish stock identification. Fish Res. 1999; 43(1):35-44

2.

Bräger Z, Moritz T. A scale atlas for common Mediterranean teleost fishes. Vertebr Zool. 2016;66

3.

Ditty JG, Houde ED, Shaw RF. Egg and larval development of Spanish sardine, Sardinella aurita (Family Clupeidae), with a synopsis of characters to identify clupeid larvae from the northern Gulf of Mexico. Bull Mar Sci. 1994; 54(2):367-380.

4.

Homayuni H, Marjani M, Mousavi-Sabet H. Descriptive key to the otoliths of three Sardinella species (Pisces, Clupeidae) from the northern Oman Sea. AACL Bioflux. 2013;6(3):211-21.

5.

Hu W, et al. Study on fish life history traits and variation in the Taiwan Strait and its adjacent waters. Acta Oceanol Sin. 2015; 34(2):45

6.

Knebelsberger T, Stöger I. DNA extraction, preservation, and amplification. DNA barcodes: methods and protocols. 2012; p. 311-338

7.

Lavoué S, et al. Phylogenetic relationships among anchovies, sardines, herrings and their relatives (Clupeiformes), inferred from whole mitogenome sequences. Mol Phylogenet Evol. 2007; 43(3):1096-1105

8.

Leung A. The epibenthic ichthyofauna of Tolo Harbour and Hong Kong’s northeastern waters: a long term record of change. Proceedings of the Eighth International Marine Biological Workshop: The Marine Flora and Fauna of Hong Kong and Southern China, Hong Kong University Press, Hong Kong. 1997.

9.

MMAF. Capture fisheries statistics of Indonesia, 2011. 2012; Jakarta: Directorate General of Capture Fisheries. p. 190.

10.

Palumbi S, et al. The simple fool’s guide to PCR, version 2.0, privately published document compiled by S. Palumbi. Dept. 1991; Honolulu: Zoology, Univ. Hawaii. p. 96822.

11.

Quilang JP, et al. DNA barcoding of the Philippine endemic freshwater sardine Sardinella tawilis (Clupeiformes: Clupeidae) and its marine relatives. Philipp Agric Sci. 2011;94(3)

12.

Sivakumaran K, Manickasundaram M, Ramaiyan V. Problems of identification among species of Sardinella. CMFRI Bulletin National Symposium on Research and Development in Marine Fisheries Sessions I & II 1987. 1987 Kochi: CMFRI. .

13.

Stern N, Rinkevich B, Goren M. Integrative approach revises the frequently misidentified species of Sardinella (Clupeidae) of the Indo-West Pacific Ocean. J Fish Biol. 2016; 89(5):2282-2305

14.

Strauss RE, Bond CE. Taxonomic methods: morphology. Methods for fish biology. 1990; Bethesda: American Fisheries Society. p. 109-140.

15.

Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994; 3(5):294-299

16.

Ward RD, et al. DNA barcoding Australia's fish species. Philos Trans R Soc Lond B Biol Sci. 2005; 360(1462):1847-1857

17.

Whitehead P. FAO species catalogue vol. 7. Clupeoid fishes of the world (suborder Clupeoidei): an annotated and illustrated catalogue of the herrings, sardines, pilchards, sprats, shads, anchovies and wolf-herrings. Rome: Food and Agriculture Organization of the United Nations; 1985.

18.

Willette D, Santos M, Aragon M. First report of the Taiwan sardinella Sardinella hualiensis (Clupeiformes: Clupeidae) in the Philippines. J Fish Biol. 2011; 79(7):2087-2094

Abbreviations

COI region

Cytochrome c oxidase subunit 1 region

Cyt-B

Cytochrome B subunit

mtDNA

Mitochondrial DNA

ND4

NADH dehydrogenase subunit 4

ND5

NADH dehydrogenase subunit 5

ND6

NADH dehydrogenase subunit 6

NGS

Next-generation sequencing