DNA barcoding of commercially important reef fishes in Weh Island, Aceh, Indonesia

Knowledge on the precise identification of fish resources is critical for sustainable fisheries management. This study employs the DNA barcoding approach to generate a molecular taxonomic catalogue of commercially important reef fishes in the waters of Weh Island (Aceh Province), the most northerly inhabited island in the biodiverse Indonesian Archipelago. The waters not only support artisanal fisheries but also a feeder for the industry in the greater island of Aceh. In total, 230 specimens from 72 species belonging to 32 genera and 17 families were DNA barcoded, representing a major segment of the captured reef fish taxa and a quarter of fish species diversity that had previously been recorded. The sequence read lengths were 639 bp revealing 359 conserved sites, 280 variable sites, 269 parsimony informative and 11 singletons. Our molecular findings paralleled the morphological identification with no evidence of cryptic species or new species discovery. This study is a significant contribution to the fisheries statistics of this area, which would facilitate assessment of species catch composition and hence for strategizing management plans. It is an important input to the DNA barcode library of Indonesian marine fishes and to the global DNA barcode entries in general.


INTRODUCTION
Weh Island (5 • 50 N, 95 • 20 E) is located north of the main island of Sumatra, Indonesia and is within the province of Aceh (widely known for the 2004 Tsunami catastrophe). Surrounded by three great water bodies: the Andaman Sea, the Indian Ocean and the Straits of Malacca, Sumatra Island is rich in biodiversity including in its surrounding marine ecosystem (Brown, 2007). This amazing feature has attracted many scientific investigations. Baird et al. (2012) and Rudi et al. (2012a) reported that the Aceh waters around Weh Island might also be a hotspot for coral biodiversity comparable to the Coral Triangle. In particular, the mean species richness of the Acropora has been estimated to be comparable to regions within the Coral Triangle, an area centered in eastern Indonesia that is considered to be the most diverse marine region in the world (Rudi et al., 2012a). The high richness has been attributed to the area being a source of overlap of three distinct groups of faunas: widespread Indo-Pacific species, Indian Ocean species and Pacific Ocean species (Baird et al., 2012).
The coral ecosystem in Weh Island waters is home to no less than 350 reef fish species community (Campbell et al., 2005;Campbell et al., 2007;Rudi et al., 2008;Rudi et al., 2012b). Allen & Adrim (2003) documented six endemic reef fish species in this area. Indonesia is recognized to have the highest diversity of reef fishes species in the world with more than 2500 reef fish species recorded (Allen & Erdmann, 2012). This biodiversity is a very important source of livelihood to the Weh Island local artisanal fishermen. However, the bigger boats plying these waters transport their catch to other major landing sites in Aceh. The majority of fishes caught in this island are reef fishes, mainly groupers (Rudi et al., 2012b).
To ensure the perpetuity of this biodiversity, particularly with its relationship to the community well-being, it is imperative that management efforts are in place. A reliable species list is one of the most crucial baseline data in developing fisheries management plan and essential for conservation and sustainable exploitation (Basheer et al., 2017;Díaz et al., 2016;Pavan-Kumar et al., 2018;Ward, 2000). Earlier reports on reef fish identification in Weh Island had been based on morphological characters, and several species identification remains inconclusive (Rudi et al., 2012b). Precise species identification of landing catches is very important for capture fisheries management. In most major fisheries countries (including Indonesia), the majority of fish landing data are not recorded based on species names but typically the fishes are registered with a local name. For instance, in Aceh, Indonesia, the groupers are collectively recorded as ''kerapu'' instead of the precise scientific species names. Another example is for some snappers in Malaysia. The Department of Fisheries Malaysia uses the term red snapper (''Merah'') to collectively refer to two species, Lutjanus malabaricus and Lutjanus sebae, while ''remong/kunyit-kunyit'' is applied to species in the yellow-lined group including Lutjanus lutjanus and Lutjanus vitta (Bakar et al., 2018). It is impossible to develop effective management strategies and conservation plans without knowing what species are involved. Inaccurate species identification would undoubtedly affect the accuracy of stock inventories. This has severe implications on the ability to sustainably manage highly exploited species when pooled with species not under threat. Thus, a high proportion of misidentifications in the catch statistics could lead to considerable economic consequences.
Since its introduction in 2003, the DNA barcoding technique has become the golden standard for molecular taxonomy (Hebert et al., 2003). DNA barcoding has been recognized for its effectiveness in species identification for its rapid and robust output based on the cytochrome c oxidase subunit I (COI) gene. This technique has been utilized in species identification of marine fishes globally (Basheer et al., 2017;Bingpeng et al., 2018;Mat Jaafar et al., 2012;Ward et al., 2005;Weigt et al., 2012;Zhang & Hanner, 2011). It has a wide range of applications including identification of species that have similar morphological characters (Bingpeng et al., 2018), individuals at different developmental stages (Hubert et al., 2010), damaged and incomplete specimens (Wang et al., 2018), processed products i.e., food forensics (Chin Chin et al., 2016;Marko et al., 2004), cryptic species and new species discovery (Lim et al., 2016;Farhana et al., 2018).
There have been a few documentations of molecular studies of freshwater fishes of Aceh e.g., DNA barcode of freshwater fishes in Lake Laut Tawar Aceh (Muchlisin, Fadli & Siti-Azizah, 2012), the genetic diversity and structure of the striped snakehead (Tan et al., 2015) DNA barcode of freshwater halfbeak genus Hemirhamphodon from Sundaland (Lim et al., 2016), DNA barcode of tropical eels of Aceh (Muchlisin et al., 2017), and DNA barcode of freshwater halfbeak genus Dermogenys in Southeast Asia (Farhana et al., 2018). However, there is no record of reef fish species identification through DNA barcoding in Aceh. Related studies in the Indonesian waters have been coincidental as part of larger projects over a wider coverage, with limited representative specimens from the Aceh region (Anzani et al., 2019;Madduppa et al., 2017;Madduppa et al., 2016;Prehadi et al., 2015;Sembiring et al., 2015;Wardiatno et al., 2015;Wibowo et al., 2018). Thus, this investigation is the first comprehensive study of DNA barcoding of commercially important reef fishes in Indonesia, with focus on the marine waters of Weh Island. It would also be an important contribution to the global DNA barcode library and will facilitate fisheries management and conservation programmes in Weh Island specifically, and the Aceh waters (and Indonesia) in general.

Sample collection
A total of 230 coral reef fish specimens were collected from fish landing sites and fish markets located throughout Weh Island (except for the west coast of Weh Island since this is a protected forest area, with no fish landing sites permitted) from May 2015 to November 2016 ( Fig. 1). Reef fishes were selected from the total landings based on the following taxonomic keys: Heemstra & Randall (1993), Craig, Mitcheson & Heemstra (2011), and Froese & Pauly (2019. Each sample was documented according to the requirements of the BOLD Systems, including digital specimen photograph. Fin biopsy was performed at the collection sites. Tissue samples were fixed with 96% ethanol placed in individual two mL microcentrifuge tube. Sample sizes of each morphospecies ranged from 1 to 9 individuals with an average of 3 specimens classified into 72 presumed species (53 species were represented by more than one specimen while 19 species by only a single individual) depending on the availability of the species. Field experiments were approved by the Faculty of Marine and Fisheries, Syiah Kuala University, Banda Aceh, Indonesia (16030/A2.1/KP/2016). The committee of Research Ethics and Animal Care of Syiah Kuala University provided full approval for this research (Ethic Code No: 958/2015). Voucher specimens have been deposited at the Faculty of Marine and Fisheries, Syiah Kuala University at Banda Aceh, Indonesia. No specific authorization were required for sampling since samples were obtained from commercial fishers and fish markets. In many cases, tissues were contributed by the fishers. For those samples which lack morphological vouchers (big fishes), ''photographic voucher'' were obtained following the Fish-BOL collaborator's protocol (Steinke & Hanner, 2011). All DNA sequences have been deposited  Table 1.

DNA extraction, primer, and PCR assay
The modified Cetyltrimethyl Ammonium Bromide (CTAB) protocol was used to isolate genomic DNA (Bakar et al., 2018). The extracted DNA was quantified using the Nano Drop 2000c Spectrophotometer (Thermo Scientific, Waltham, MA). The COI gene was amplified using the universal primer sets of Ward et al. (2005). PCR was set up in a 25 µL mixture reaction and was amplified in a BIO-RAD T100 Thermal Cycler (BioRad Laboratories Inc., USA). The total reaction volume of 25 µL contained 2.0 µL DNA template, 0.5 µL of each primer, 2.5 µL of 10x i-Taq TM plus PCR Buffer, 2.0 of 25 mM MgCl 2 , 1.0 µL of dNTP, 0.25µL of i-Taq TM plus DNA Polymerase and 16.25 µL of Milli-Q water. The thermal conditions consisted of an initial denaturation at 95 • C for 2 min followed by 30 cycles of denaturation at 94 • C for 45 s; annealing at 47.9-60 • C for 45 s; elongation at 72 • C for 1 min and a final extension of 72 • C for 10 min before termination of the reaction at 4 • C.

Gel electrophoresis, staining, and DNA sequencing
Amplicons were separated and visualized on a 1.7% agarose gel stained with 2 to 2.5 µL of RedSafe TM Nucleic Acid Staining Solution (IntRON Biotechnology, Gyeonggi-do, Korea). A 2.0 µL volume of PCR product was loaded onto the agarose gel and electrophoresed at 100V for 30 min to assess its quality. Satisfactory PCR products were sent to the First BASE

Data analysis
Forward and reverse COI sequences were aligned in Clustal W and edited in MEGA 6.06 (Tamura et al., 2013). To ensure accurate alignment and detection of stop codons if present, the aligned sequences were translated into proteins. To confirm their genetic identity, the aligned sequences were compared with previously known sequences and closely related species in BLAST (Basic Local Alignment Search Tool; http://www.ncbi.nlm.nih.gov/BLAST) and the BOLD Identification System (IDS) (http://www.boldsystems.org/). Pairwise genetic distances at different taxonomic levels (within species, within genus, within family) were estimated based on the Kimura-2-parameter (K2P) model using MEGA 6.06 (Kimura, 1980;Tamura et al., 2013) and Sequence Analysis Engine of BOLD (http://www.boldsystems.org/). The MEGA 6.06 software was also used to cluster COI haplotypes into a Neighbour-Joining (NJ) phylogeny, employing 1000 bootstrap replicates. Furthermore, DnaSP 5.10 was used to summarize haplotype distributions (Librado & Rozas, 2009;Rozas et al., 2003).
Comparisons at the species level of the maximum intraspecific genetic distance with the minimum distance to the nearest neighbour were performed by applying the BOLD's 'Barcoding Gap Analysis' tool (http://www.boldsystems.org). The Automatic Barcode Gap Discovery (ABGD) (with default settings and the K2P model) was used for species delimitation to determine the number of operational taxonomic units (OTUs) based on pairwise sequence distances between individuals within the dataset (https://bioinfo.mnhn.fr/abi/public/abgd/abgdweb.html) (Puillandre et al., 2012).
Further to the molecular work, the conservation status of the identified reef fishes was determined based on the International Union for Conservation of Nature (IUCN) website (http://www.iucnredlist.org, accessed on 13 April 2020). The IUCN criteria are widely accepted for assessing species conservation status (Hoffmann et al., 2008) including in the evaluation of the conservation status of plants (Callmander, Schatz & Lowry, 2005), corals (Carpenter et al., 2008), sea cucumbers (Madduppa et al., 2017) and several fish species (Morris, Roberts & Hawkins, 2000;Sadovy de Mitcheson et al., 2013), among others. The IUCN categorizes the conservation status into four groups namely: Threatened (T), Near Threatened (NT), Data Deficient (DD) and Least Concern (LC) (further description of each category can be found in http://www.iucnredlist.org). Besides their conservation status, their trade status was also determined using CITES (Convention on International Trade in Endangered Species) online at http://www.cites.org/ (accessed on 13 April 2020).

Cytochrome oxidase subunit I diversity assessment
All 230 specimens were successfully sequenced generating 146 haplotypes. The sequence read lengths were 639 bp with average nucleotide composition of A = 23.78%, T = 29.58%,  C = 28.39%, G = 18.25%. There were 359 conserved sites and 280 variable sites, of which 269 were parsimony informative and 11 singletons. The base composition showed that the AT content (53.36%) was higher than the GC content (46.64%). The mean T content was the highest, and the mean G content was the lowest. The GC content decreased in the order of first, second, to the third position codon with mean values of 57.09%, 42.70% and 40.14% respectively. Of the 230 COI sequences, no insertions, deletions or stop codons were detected in the samples. In addition, intra-taxon K2P genetic distances increased from lower to higher taxonomic levels: 0% to 1.55% within species, 3.79% to 20.73% between species within genus and 9.79% to 23.94% between genera within family. The mean K2P genetic distance between specimens was 0.24% within species, 11.16% between species within genus and 19.01% between genera within family (Fig. 2, Table 2). The average K2P genetic distance between species within genus was 46.5 fold higher than the average K2P genetic distance within species.

Species delimitation
The assessment of species identities with previously known sequences and closely related species in BLAST and BOLD databases yielded 98-100% identities indicating the efficacy of COI sequences to provide species-level resolution. In addition, Barcoding Gap Analysis showed that all putative species had a maximum intraspecies distance of less than 2%. The mean distance to the nearest neighbour (NN) was 11.13%. The latter analysis also showed that all specimens exhibited high distance values to their nearest neighbour indicating the presence of ''barcode gap'' among the 72 observed species (Fig. 3, Table 3). However, it should also be noted that there were several single-specimen species (not applicable for this analysis) and numerous NN that were not congenerics. In concordance with Barcoding Gap Analysis, ABGD also generated 72 operational taxonomic units (OTUs) with the initial partition at a prior intraspecific divergence (P) (P = 0.0077-0.03594) (Fig. 4). Furthermore, the NJ tree showed that all the recognized species formed monophyletic clusters without any overlap between species. All assemblages of conspecific individuals had a bootstrap support of 99% and 1.00 posterior probability (Fig. 5). Congeneric samples were also well resolved.

Conservation and trade status
Based on the IUCN categories, of the 72 species sampled in Weh Island, 61 species are categorised as Least Concern (85%), one species (1%) is Data Deficient, ten species (14%) categorized as Not Evaluated (NE) while none (0%) are deemed to be under Near Threatened category. However, two species showed decreasing population trend (L. decussatus and V. albimarginata) (http://www.iucnredlist.org, accessed 13 April 2020). None of the 72 identified species has been evaluated by CITES (http://www.cites.org/, accessed 13 April 2020).

Establishment of DNA barcode for Weh Island Reef Fishes
This study has demonstrated the utility of DNA Barcoding to complement the morphological identification of 72 reef fish species from Weh Island. The establishment of a reliable DNA barcode reference records for commercial reef fishes in Weh island has been initiated, contributing to the DNA barcode library of Indonesian marine fishes and the worldwide DNA barcode entries in general. This baseline database is critical for the future fisheries management on this island. Database sequence similarity and genetic distance comparisons with voucher reference supported the accurate identification of the 72 putative species. The exact or near matches (98% to 100%) identity with reference DNA libraries both in BLAST and the BOLD Identification System (IDS) is regarded as firm evidence in the success of the DNA barcoding approach (Alcantara & Yambot, 2016;Bhattacharjee et al., 2012;Cerutti-Pereyra et al., 2012;Filonzi et al., 2010;Joly et al., 2014;Kress & Erickson, 2012;Ratnasingham & Hebert, 2007;Ward, 2009;Zhang & Hanner, 2012). Species delineation was also supported by ABGD and the monophyletic clustering on the NJ tree of all morphologically recognized species.
As expected, and in concordance with previous studies a marked elevation of conspecifics versus congeneric genetics distances were observed although the values varied considerably.  (31X). However, the current study showed a higher value of 46.5-fold. This may be attributable to several factors. Firstly, the conspecifics were also 'con-population' ie. from the same population. Although sampled from different landing sites, the populations were most likely a single population due to the proximity of all the landing sites. Previous studies may have represented conspecifics sampling from different geographical populations, with higher potential for genetic variation, consequently increasing the conspecific distance, reflected by the narrower gap between the conspecific versus congeneric distances. However, this factor does not explain the discrepancy as the conspecific distances in our study were of the same magnitude in the current and previous studies. Thus, the elevated comparisons are likely attributable to higher between-species distances i.e., our study sampled a lower number of closely related species.
The ''barcoding gap'' analysis also supported the delineation of the samples into 72 taxa. Even though there is no universal threshold to delineate species, Meier, Zhang & Ali (2008) defined a barcoding gap as the difference between the smallest congeneric and maximum conspecific divergence. The barcode gap results offer yet another explanation for the elevated congeneric to conspecific differences in our study. As observed (Table 3), when the nearest neighbor is not a closely related species (intra-genera) but belonged to another genus or family, the genetic distance was high (above the mean) and therefore caused biasness in the expected NN values. Furthermore, since the number of species sampled were much lower (72 from the 350 morphologically documented fish species) than the record from previous studies (Campbell et al., 2005;Campbell et al., 2007;Rudi et al., 2008;Rudi et al., 2012b), it is likely that many ''true'' nearest neighbours were not sampled, thus inflating the congeneric mean.
Reef fishes are among one the most taxonomically challenging groups due to their astonishing diversity of morphological characteristics and colouration, often environmentally influenced (Duarte et al., 2017;Hubert et al., 2010). Their morphometric traits may also change during the stages of ontogenetic development (Bingpeng et al., 2018). Furthermore, traditional morphological classification methods (especially colour) depend on well-preserved samples which are frequently lacking from samples collected from landing sites and fish markets. The majority of fishermen in Weh Island rely on traditional chilling methods for fish preservation (e.g., icing the fish with ice flakes or ice tubes) which often result in loss of colouration and often incomplete/damaged morphological characters, thus making taxonomic identification more complicated. However, there appears to be no issue in accurate morphological identification of the samples. One obvious reason for the perfect match in morphological and molecular analysis is the familiarity with these common commercial species of the island.

Conservation and trade status
At a glance, the conservation status for each of the 72 barcoded species does not raise any cause for ''alarm'' with 85%, 14% and 1% categorized into LC, NE and DD, respectively. The majority (LC) do not appear to require any additional protection as required for Critically Endangered, Endangered, Vulnerable or Near Threatened categories (IUCN, 2012). However, sustainable considerations are needed in exploiting the species listed in NE and DD categories as these have no or limited biological, ecological or distributional information for a precise assessment to be made using the IUCN criteria (Sadovy de Mitcheson et al., 2013). Until such data is available and conservation status verified, it would be prudent to accord these group careful attention, at least until their status are evaluated (IUCN, 2012). Disregarding the management of LC category is also ''unsafe'' as they make up 80% of the bony fish and therefore the fisheries of this area (IUCN, 2019). In coral reef fisheries, ''fishing through the food web'' is a common fishing practice and Weh Island fisheries is no exception. The fishermen catch both high and low trophic species since all fish species have a market value (Darling & D'Agata, 2017). The dominant families landed are: Epinephelidae, Scaridae, Holocentridae, Labridae, and Lutjanidae. Epinephelidae (groupers) and Lutjanidae (snappers) are at the highest trophic level and the most valuable. In Aceh, the price of one kg grouper may reach 7-10 USD and therefore highly prized by the fishermen. They are being threatened on a global scale due to overfishing (Sadovy de Mitcheson et al., 2013). Being predators, their declination from the ocean may lead to imbalance in the trophic chain. Scaridae (parrotfish) are sold at a more affordable price (one kg may cost 2-3 USD) which makes it more accessible to the general population and thus also at a risk of overharvesting. They are low trophic herbivore species feeding on turf algae that grow on reef substrate and hence function to clear new sites for coral larval settlement (Bonaldo & Bellwood, 2009). A decline in herbivores such as the parrotfish would hinder coral replenishment. Bozec et al. (2016) suggested a 10% harvest limitation and a minimum size of 30 cm of parrotfish catches might sustain the health of reefs in the Carribean. Without a sound management plan and immediate implementation, there is a possibility that reef fish catches in this region will go into further decline due to the overexploitation, causing economic losses and long term ecological damage.
In future, DNA barcoding in the different fields of fisheries in Indonesia could be routinely utilized as fishery managers familiarize with the technique and it becomes affordable. Ardura, Planes & Garcia-Vazquez (2013) suggested that DNA barcoding be included in routine surveys and periodical monitoring of catch diversity. Through precise species identification, DNA barcoding could be applied in ecosystem management and conservation, ichthyoplankton identification, seafood forensics, biosecurity, and invasive species detection, predator-prey relationship, seed identification, fingerling survival and health management (Pavan-Kumar et al., 2018). However, with all its advantages, DNA barcoding also has several limitations, among others: this method requires a reliable DNA library so that newly generated sequences could be compared with the existing databases (Taylor & Harris, 2012). Nuclear mitochondrial pseudogenes (numts) that may be coamplified with mtDNA when using conserved universal PCR primers could lead to erroneous estimates of species diversity (Song et al., 2008). Furthermore, the method is inefficient for hybrid detection (Dudu et al., 2016). Occurrence of hybridization has been documented among a broad range of teleost fishes, both freshwater (Scribner, Page & Bartron, 2000) and marine (Montanari et al., 2012;Qu et al., 2018;Yaakub et al., 2006). The COI gene is of mitochondrial origin, typically derived only from the maternal donor. For this reason, any hybrid fish will be inevitably and erroneously identified as its maternal species (Bhattacharya et al., 2016;Linacre et al., 2011). Thus, if the specimen is a hybrid, the COI barcoding result will enable only the identification of the female parent rather than male (Kochzius et al., 2010). In order to allow the precise identification of the parental pair of a hybrid, additional nuclear markers are suggested to be employed (Bhattacharya et al., 2016;Kochzius et al., 2010;Linacre et al., 2011). Montanari et al. (2012) confirmed the hybrid status of the butterflyfish, genus Chaetodon through the utilisation of 20 genus-specific microsatellite loci to complement the mtDNA results. Yaakub et al. (2006) utilised a multi-loci approach based on mitochondrial COI gene, nuclear S7 (intron 1) and nuclear copy of mitochondrial (numt) D-loop region to corroborate the identity of Thalassoma quinquevittatum as the maternal and T. jansenii as the paternal contributor in a study on hybrid specimens (Labridae: Thalassoma) collected from Holmes Reef, an isolated oceanic atoll in the Coral Sea. The nuclear copy of the mitochondrial genome (numt) was detected in all putative parental species and hybrids, further supporting the evidence from the nuclear S7 intro 1. In another study (Qu et al., 2018) the nuclear RYR3 marker was found to be an efficient marker to identify first-generation grouper hybrids to assist barcoding in family Epinephelidae.
Finally, results of the present study have contributed significant records for the molecular taxonomy of commercially important reef fishes from Weh Island. Sequences from the study have been deposited in the BOLD Systems websites, establishing a molecular biodiversity catalogue of reef fishes from Weh Island and the Indonesian Archipelago in general. The data could facilitate improved monitoring, conservation, and management of fisheries in this area. When funds permit, more intensive studies should be conducted with a wider coverage and through direct sampling.

CONCLUSIONS
The establishment of a DNA barcode reference library of important reef fishes from Weh Island in the Indonesian Archipelago has been achieved through this study. In total, 230 specimens from 72 species belonging to 32 genera in 17 families were barcoded. There are no overlaps of pairwise genetic variations between conspecific and interspecific comparisons. Epinephelidae (25%), Scaridae (13%), Holocentridae (11%), Labridae (7%) and Lutjanidae (7%) were the dominant family found in Weh Island landing sites, respectively. The study has a significant contribution to the DNA barcode library of Indonesian marine fishes and the worldwide barcode entries in general. It will also contribute towards achieving more efficient monitoring, conservation, and management of fisheries and has a potential role in the aquaculture industry of this region.

ADDITIONAL INFORMATION AND DECLARATIONS Funding
The research was funded by LPSDM Aceh (Contract number: 189/S3.LN/LPSDMA/X/2015) and Universiti Malaysia Terengganu (Addressing Fisheries and Aquaculture Challenges in the Grouper (Family Epinephelidae) through Molecular Genetics with Special Focus on E. coioides -UMT/RMIC/STRA/53280). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Grant Disclosures
The following grant information was disclosed by the authors:

Competing Interests
The authors declare there are no competing interests.

Author Contributions
• Nur Fadli and Siti Azizah Mohd Nor conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft.
• Ahmad Sofiman Othman analyzed the data, authored or reviewed drafts of the paper, and approved the final draft.
• Hizir Sofyan analyzed the data, prepared figures and/or tables, and approved the final draft.
• Zainal A. Muchlisin performed the experiments, prepared figures and/or tables, and approved the final draft.

Animal Ethics
The following information was supplied relating to ethical approvals (i.e., approving body and any reference numbers): The Committee of Research Ethics and Animal Care of Syiah Kuala University provided full approval for this research (Ethic Code No: 958/2015).

Field Study Permissions
The following information was supplied relating to field study approvals (i.e., approving body and any reference numbers): Field experiments were approved by the Faculty of Marine and Fisheries, Syiah Kuala University, Banda Aceh, Indonesia (16030/A2.1/KP/2016).