Comparative study of DNA extraction to initiate harmonized protocol for a simple method of species identification: Fresh and canned Tuna case study

ABSTRACT Seven different methods of DNA extraction were performed and compared in order to have the better DNA recovery in terms of quality and quantity. A PCR assay was realised based on the Cytb mitochondrial gene to identify 53 commercial food products of wide range diversity (Scombridae Thunnus species) commercialised in the Tunisian market, and we checked whether these products were correctly labelled under European and national legislations. Phylogenetic analyses on Cytb mitochondrial gene were used to study the relationships among the considered species. Tuna species was easily differentiated and confirmed by direct sequencing. The results showed a high deficiency in labelling: 96% of canned tuna did not specify the species, although the name of the species had to be mentioned. Partial least squares regression (PLS) showed relatively high values (45%) in all analysed data, suggesting substantial difference between the extraction procedures and matrices with SDS/PCI method being the most convenient. Graphical abstract


Introduction
Seafood and aquaculture products are among the most vulnerable produces to food fraud. Therefore, traceability in these relative sectors is essential to prevent fraud, guarantee compliance with food safety regulations and enforce sustainable fisheries management. Traceability is also critical for ensuring the quality of fish products and minimising health risks for consumers (FAO, 2018;Van Ruth & Granato, 2017;Velasco et al., 2020). Tuna fisheries and processing represent economic activities of primordial importance around the world. They are greatly demanded commodities, and most of these products are internationally traded for human consumption. However, not all tuna products have the equal market price as for certain species, such as Thunnus thynnus though not always available, there are consumers who are even prepared to pay huge sums, while other species are quite affordable (for instance, Skip jack). Therefore, specious labelling has been commonly detected (Sotelo et al., 2018). The best instrument used to combat seafood fraud is through legislation (Ferrito et al., 2019). For instance, the European labelling regulation aims to provide consumers information on seafood product including the commercial and scientific names, to ensure their traceability and identification along the value chain (EU 1379(EU /2013. However, such regulation excludes prepared and processed aquatic product from the application of the mandatory labelling and therefore remains a gap to be considered by the legislator (D'Amico et al., 2016). For example, in the reglementation (EEC 1536/92) where it is specified that preserved tuna and bonito must be prepared exclusively from certain fish species (i.e. tuna cans should contain only any Thunnus species or Katsuwonus pelamis), mixing of species is not allowed in each tin (unless muscular structure has disappeared), only commercial names are required.
To overcome these problems, there are many recent genetic methodologies available for this purpose (Gordoa et al., 2017;Mata et al., 2020;Sotelo et al., 2018). Most of the genetic approaches used for the determination of species identity are based on the amplification of a region of mitochondrial (mtDNA) DNA by polymerase chain reaction (PCR; Piskata et al., 2017;Xiong et al., 2016;Yang et al., 2016). Cytochrome b is one of the familiar mitochondrial genes used in DNA barcoding. The mitochondrial cytochrome b gene is widely used in systematic studies to resolve divergences at many taxonomic levels. In addition, cytochrome b has some advantages such as its resistance against high temperature, abundancy in cell, following maternal heredity, as well as its high mutation. The authentication of tuna fresh, steak, as well as canned tuna using cytochrome b is reported in previous study (Besbes et al., 2011;Pedrosa-Gerasmio et al., 2012;Abdullah & Rehbein, 2014;Wulansari;Wulansari & Nurilmala, 2015).
Several protocols for DNA extraction from fish samples have been described including chemical reagents and commercial kit (Chapela et al., 2007;Besbes et al., 2011;Piskata et al., 2017). In general, fresh tissue of fish that has not been thermally treated has the best probability of yielding adequate DNA. However, the detection of DNA in processed products can be adversely affected by DNA fragmentation during the processing steps and the use of ingredients that may inhibit the PCR reaction. The critical step remains the extraction of high-quality DNA in sufficient large quantities from heterogeneous food matrices. In fact, DNA amplification on which rely the effectiveness of the method for species identification requires expensive reagents, which are not always available in laboratories especially in developing countries. Therefore, the choice of the extraction method should consider the reagent efficiency in term of DNA quantity/purity as well as their relative costs . In this study, we aimed in the first step to extract DNA using different protocols to initiate harmonised method for species identification among Mediterranean area including Tunisian and Spanish laboratories. In a second step, PCR of a 250-bp fragment in the mitochondrial Cytb region followed by a phylogenetic analysis were applied for the identification of species in 53 canned tuna brand commercialised in Tunisian market.

Commercial food products
The fresh tuna samples used in this study were caught along the Tunisian coast and identified morphologically by researchers; Tuna-type products were purchased at local supermarkets or directly from Tunisian companies specialised in the commercialisation of aquatic food products. The scientific names are given in Table 1. All samples were frozen at −80°C upon arrival in the laboratory and canned tuna was kept at room temperature. Upon arrival in the laboratory, all samples were stored at −80°C until analysis.

DNA extraction
The DNA of tuna samples was extracted with the sterile scalpel from tissue of each sample according to the different protocols (EX1-EX7).

Kit fast DNA (EX1)
The DNA was extracted following the Kit -manufacturer's instructions. To the weighed subsample (0.1 g); 800 µL Genomic Lyse buffer was added (to lyse the cells and solubilise the proteins, DNA and other cellular constituents) in the presence 0.1 mg (3.2 U) of Proteinase K (to digest the protein) and incubated at 65°C for 10 to 30 minutes. The DNA-containing solution was mixed with the genomic binding buffer (100 µL), then loaded on a column which binds the DNA.

Method sing the Phenol-chloroform-isoamyl Alcohol: SDS/PCI (EX2)
The extraction of DNA was performed according to (Besbes et al., 2011) with the following modifications. As tuna tissue contains fewer lipids, a one-step DNA extraction was performed. Sample (0.1 g) was homogenised with 500 µL of the buffer A (50 mM Tris-HCl, 100 mM EDTA and 1% SDS with pH 8.0) and 0.1 mg (3.2 U) of proteinase K (Sigma Aldrich, France). The mixture was incubated at 55°C with continuous shaking during 4 h, cooled on ice and then centrifuged at 12,000g for 10 min. The supernatant was extracted with phenol-chloroform-isoamylic alcohol in a 25:24:1 ratio, then precipitated with 600 µL of isopropanol at −20°C and mixed by inversion. Using ethanol (70%), the pellet was washed and dried, then resuspended in 50 µL TE buffer.

Method Using the SDS/PCI/NaCl (EX3)
In this case, the same procedure was applied to 0.1 g of sample using the SDS/PCI method (Hsieh et al., 2005), but with the addition of NaCl to the buffer A (final concentration of 0.2 M NaCl).

Method Using chelex (EX4)
The sample (0.1 g) was homogenised in 300 µL of a Chelex water solution (5% -Chelex 100 resin Qiagen) mixed with 20 µL of proteinase K (0.1 mg) and 30 µL TE buffer; incubated at 56°C during 4 h and followed by a heating (100°C) for 15 min (Jérôme et al., 2003). The solution was centrifuged (12,000 g for 10 min), the DNA pellet for PCR amplification was stored at 4°C until use.

Method Using TRITON/CI Chloroform-isoamyl Alcohol (EX6)
The muscle tissue (0.1 g) was homogenised in 500 µL of pH 8.0-buffer B (2 mM EDTA,10 mM Tris-HCl), to which was added 50 µL ß-mercaptoethanol 0.2% at final concentration and 50 µL Triton X-100 at 1% final concentration and 0.1 mg (3.2 U) of proteinase K (Besbes et al., 2011). The mixture was incubated at 65°C for 12 h and then centrifuged at 12,000 g (4°C-10 min). The suspended DNA was mixed vigorously with an equal volume of chloroform-isoamylic alcohol and then centrifuged at 12,000 g (10 min). For DNA precipitation, ethanol and sodium acetate were added to the supernatant, the pellets were then resuspended in 50 µL TE buffer.

Method Using dneasy tissue kit (EX7)
DNA was extracted by means of a commercial kit (DNeasy Tissue kit, Qiagen, Darmstadt, Germany), as previously described (Prado et al., 2005), 25 mg of sample were homogenised with 180 µL of lysis buffer ATL and 20 µL proteinase K. The homogenate was incubated during 3 h at 56°C and then Vortex for 15 s. 200 µL Buffer AL were then added to the sample, and thoroughly mixed by vortexing. The supernatant, obtained after centrifugation at 6000 g for 5 min at room temperature, was added to 200 µl ethanol. The mixture (400 µL) was then loaded on the DNeasy Mini spin column where DNA was bonded to the silica-gel membrane during the centrifugation steps. The bounded DNA-DNeasy Mini spin membrane was washed using two different washing buffers and then eluted using the elution buffer.

DNA quantification and purity
Extracted DNA was estimated by absorbance at 260 nm. One optical density 'DO' at 260 nm equalling 50 µg/mL DNA. Purity was checked by determining the ratio A260/A280 (Sambrock et al., 1989).

Primer design
About 53 tuna sample Cytb gene sequences were obtained from GenBank. In addition, Cytb gene sequences of tuna species purchased for this study were determined by a direct sequencing analysis of the amplified products using the sequencing primer pairs, namely, TunaF1 (forward primer 5ʹ-TACTGAGGCGAAGGCTGATT-3ʹ) and TunaF2 (reverse primer 5ʹ-CAAAGAAGATTACGCTGTTATCCC-3ʹ). This primer pairs based on the highly conserved sequence among the target species, for detecting tuna, were designed to amplify the Cytb gene. For all species, one set of primer pairs TunaF1/TunaF2 were used in the PCR amplifications to obtain 250bp DNA fragments of Cytb.

PCR amplification
PCR amplifications were carried out using Bio-Rad MyCycler (Bio-Rad Amplification conditions were as follows: first an initial denaturation step of 94°C for 3 min, followed by 35 cycles of denaturation (30 s at 94°C), annealing (30 s at 55°C) and extension (40 s at 72°C), and a final extension for 5 min at 72°C was also included. The set of newly designed primers used for PCR amplification were designed as follows TUNAF1: (forward primer 5ʹ-GCGGTATTTTGACCGTGCGA AGG-3ʹ) TUNAR2: (reverse primer 5ʹ-CCTTAATTCAACATCGAGG TCGC-3ʹ).

Cleanup and sequencing of PCR products
Before sequencing, double-stranded PCR products were purified by filtration through a DNA Clean & Concentrator-25 of ZYMO RESEARCH CORP according to the manufacturer's protocol.
PCR fragments were used for direct cycle sequencing with the dye terminator cycle sequencing kit (Beckman). Sequencing analysis was performed with a Beckman Coulter CEQ 8000 DNA sequencer in both directions with the primers used for PCR amplification

Sequence alignment, genetic distances, and phylogenetic analyses
The DNA sequences were edited with BioEdit software (Hall, 1999). In addition, the sequences were also compared to the Genbank data base by BLAST (Altschul et al., 1990) using the high similarity algorithm. The sequences were aligned by Clustal W, then the generated cytochrome b sequences were identified in GenBank (www.ncbi.nlm.nih.gov) using the BLASTn search tool.
Phylogenetic and molecular evolutionary analyses were conducted using MEGA X (Kumar et al., 2018(Kumar et al., , 2004. Nucleotide divergences were computed using the Tamura-Nei model, (Tamura & Nei, 1993), which takes substitutional rate biases and the inequality of base frequencies into account. Phylogenetic trees were constructed using the neighbor-joining (NJ) method, (Saitou & Nei, 1987) and the robustness of topology nodes was tested by the bootstrap method with 1500 iterations were shown in the branches (Altschul et al., 1990). Genetic distance was measured using the Tamura-Nei model with number of nucleotide substitutions per site (Bartlett & Davidson, 1992).

Multivariate analysis
For each extraction, the results were presented as mean ± standard deviation (n = 6). All the data were statistically treated using ANOVA with an SPSS Statistical Software System version 11.01 (SPSS Chicago, IL, USA). Tukey's test was used to determine the possible significant differences among mean values at the 5% level. The results were also analysed by partial least squares (PLS) regression using multivariate statistical software (The Unscrambler version 9.8, CAMO Software AS, Oslo, Norway). Leverage correction and cross-validation of all the data were applied. The variables were weighted with the inverse of the standard deviation of all the objects in order to compensate for the different scales of the variables. Martens Uncertainty test was used to examine the influence of considered variables on the PCR amplification.
Finally, (DO A260/A280, DNA YIELD and tuna sample type) were also analysed by principal component analysis (PCA) using Past Software (Past 3.20, University of Oslo). Variables were scaled by columns so they have a mean of zero and a variance of one in order to compensate for their different scales.

Evaluation of DNA yield, purity and amplificability
The different samples used to evaluate the DNA yield and purity, included the tissue of fresh tuna and processed products (Table 1). Gel electrophoresis and spectrophotometry analysis revealed significant differences (p < .05) in DNA yield and purity according to the procedure of extraction and sample matrix ( Table 2). In the process of DNA extraction, three steps are necessary including breaking cell walls, nucleic material release from the nuclei into the media and DNA precipitation. In this study, the highest precipitated amount of DNA was obtained following the PCI/SDS protocol and that, independently of fresh or processed tuna. However, adding NaCl to the PCI/SDS medium induced a significant decrease in DNA levels in all analysed samples. Such result could be due to the inadequate concentration of NaCl used in this experiment as it has either dispersive or binding effects on DNA fragment causing its precipitation along with proteins (Bauer et al., 2003;Besbes et al., 2011). However, no optimisation of NaCl concentration was performed in this study as sufficient quantity of DNA level with good purity (A260/ A280 ratio >1.81) was obtained in all cases.
Significant lower DNA contents were obtained for canned tuna when using specific method such as the Chelex method. Moreover, for such processed samples, both methods yielded poor DNA quality (A260/A280 < 1.21) which may be attributed to interferences with substances such as protein and flavouring aromatic compounds added to tuna preparation.
Similarly, CTAB methods revealed considerable low DNA quantities with poor DNA quality (A260/A280 = 1.51-1.67, respectively). In addition, poor reproducibility of DNA extraction was noticed as generally high standard errors were found for DNA quantification for fresh and processed tuna.
Among all trials, Fast DNA kit method yielded the lowest DNA recovery in all samples with values ranging from 20.8 to 42.8 ng DNA/µL of extract. Such kit was found inappropriate to DNA precipitation and recovery in tuna matrix.
The quantity and quality of the extracted DNA from tuna clearly depended on the type of matrix and on the level of processing.
Several parameters including the degree of DNA damage occurring during processing and the extraction procedure, as well as the average length of DNA fragments, affect considerably the success of DNA amplification. In this study, the DNA extracted from the different tuna matrices were differently amplified with a set of primers amplifying the 250 bpcytochrome b fragments as shown in Table 1 and Figure 1. The amplification was successful particularly in fresh tuna, which was used as reference, but it was difficult in processed tuna; proving that processing may damage the DNA. Thus DNA-fragments from fresh tuna tissues were amplified by PCR independently of the extracting methods (Table 2). Concerning the DNA from the various processed tuna matrices, only DNA extracted with the PCI/SDS and DNeasy Tissue kit gave a PCR product of the 250 bp-expected size for all the analysed samples (Figure 1).

PCR amplification of thon mtDNA from commercial food products
Several parameters including the degree of DNA damage occurring during processing and the extraction procedure, as well as the average length of DNA fragments, affect considerably the success of DNA amplification. In this study, the DNA-fragments from fresh tuna tissues were amplified by PCR independently of the extracting methods. Concerning the DNA from the various processed tuna matrices, only DNA extracted with the PCI/SDS and DNeasy Tissue kit gave a PCR product of the 250 bp-expected size for all the analysed samples (Figure 1). The PCR products were individually amplified by primers TUNAF1 and TUNAR2. Including 2 primers, the length of each consensus mitochondrial cytochrome b gene fragment was 250 bp, as shown in Figure 1.
In this study, a total of 53 different commercial food products of TUNA were subjected to DNA extraction and purification. DNA amplification was successful in all primers and all samples analysed (Figure 1). Accordingly, the size of the molecular target 250 bp proved to be accurate to achieve PCR amplification even in processed foods that included tuna as ingredients and that have been subjected to a variety of technological treatments such as boiling or freezing.

Species identification of commercial TUNA food product
The best labelling that can be found for the different fresh tuna found in the local market was the species common name such as the Thunnus thynnus. The other canned tuna products were labelled tuna or did not have any species or genera specification on their labels. In 43 of such processed products, the tuna species included were generically referred to as "tuna" ( Table 1). The DNA amplification products were sequenced and aligned against the Cytb sequences of tuna obtained from the NCBI GenBank (Table 1).
The sequences acquired in this study were deposited in GenBank with accession numbers: (MT240963-MT241007/ MT274604-MT274604) and all the sequences of the tunafood products compared and identified by BLAST (Table 1).
Authentication of Thunnus species with Cytb as a "barcode gene" by BLAST showed the high similarity, this is exemplified by reference sample of K.pelamis (TKP), which gave the best hit (98% identity) to KP669125.1 (T. albacares), and a percent ID of 99% for the next following tuna species (T. thynnus, KJ623830.1).
The new sequences, which were equally deposited in the NCBI GenBank, were aligned against the other scombridae species sequences allowing the construction of a neighbourjoining tree (NJ; Figure 2).
The European labelling regulations establish compulsory information such as the scientific name and the corresponding commercial denomination according to the official list of names of each species tuna.
Thus, the constructed phylogram showed two distinct family clades (tuna and bonito), for the scombridae families, each associated with each reference species.
We found that the 24 canned tuna products were in the same cluster with K. pelamis and these tuna products exhibiting the highest values of intraspecific homology (99%) and the bootstrap value were high. Additionally, when we compared the sequences of the tuna-food products by BLAST (Figure 2), 24 products were also identified as K.pelamis.
The rest of tuna-processed products were grouped together with three cluster Thunnus albacores (seven canned tuna products were grouped in the same cluster of T.albacares reference with a bootstrap values at the nodes of 84%), Euthynnus alleteratus (five canned tuna products were grouped in the same cluster of E.alleratus reference with a bootstrap values at the nodes of 99%) and Thunnus thynnus (nine canned tuna products were grouped in the same cluster of T.thynnus reference with a bootstrap values at the nodes of 99%) in the tree (Figure 2).
From the NJ tree (Figure 2), we found two tuna type products (T43 MT240996 and T46 MT240997) were mislabelled. The alignment results against available online database in GenBank showed a despite the high similarity of the Thunnus maccoyii and Thunnus tonggol, respectively, thus confirming mislabelling of samples.
Species divergence is another tool to evaluate similarity of an unknown sequence with a pool of reference samples. In this study, the genetic distances of Auxis rochei was easily differentiated from the studied Katsowonus pelamis, Thunnus thynnus and Euthunnus alleteratus species was 0.96, 0.94 and 1.04, respectively (Supplementary material Table 3).
The degree of unlabelled or unspecific cases, similar to the findings of other studies on tuna, was observed from the beginning of the commercial chain. These results suggested that labelling is overlooked by control inspections (Gordoa et al., 2017).
The results showed that the use of the unspecific generic name hamper the detection of species substitution. This is also covered by the attitude of consumers, who rarely ask about the tuna species at the time of buying, because in Tunisia the public are usually unaware of the existence of other tuna species and are not well informed about fish products and policy. In this scenario of incomplete information and consumer ignorance, deception conditions are favourable (Gordoa et al., 2017), as was proved in this study.
Consumers are misinformed about tuna species sold, which makes them unable to differentiate between species and has a negative financial impact on them. We conclude that the Tunisian regulations on labelling and on bluefin catch certificates have proven to be ineffective, highlighting the need for policy enforcements and a review of the accepted commercial names. Several actions should be taken, efforts should be made to have more effective control on labelling and labelling authenticity, including genetic testing, and on the required documentation, and it needs to be accompanied by dissuasive sanctions, proportionate to deception turnover and its associated risks and economic losses. There is also an urgent need for information campaigns addressed to Tunisian consumers: species substitution must be fought from above with efficient control and from below through well-informed and conscientious consumers.
This study showed that DNA barcoding is a fast and high accuracy tool to identify the species in order to avoid economic fraudulent such as illegal trading and mislabelling as reported in previous study (Nurilmala et al., 2013;Sotelo et al., 2018).

Partial least square regression analysis
Due to the complexity of data originating from the different extraction methods and the different matrices, multivariate statistical procedure such as principal component analysis (PCA) and Partial least squares regression models (Martens & Naes, 1989) were run to identify the variable most significant to the PCR amplification of tuna product In this study, partial least squares regression models (PLS) were performed to identify the most significant variable relative to the different extraction method to PCR amplification of fresh and tuna products (Martens & Naes, 1989). The PCR amplification was used as Y-matrix; all the others parameters (DO, A260/A280, DNA YIELD, extraction methods and tuna sample type) as X-variables. The loading plot of tuna (Figure 3), corresponded to a variance in PCR result of 32%, of which 30% was along PC1 and 2% along PC2. Of the total variation of all other variables of 30%, 15% was along PC1 and 15% along PC2. Such result highlighted that fresh tissue (right upper region) is a discriminating parameter. The comparison of the different Tuna products revealed that DO A260/A280, DNA YIELD, fresh tuna (right hand side) and processed oil tuna (left hand side) contribute the most to important variation on PC1. The PC2 was mainly describing processed tuna product.
When comparing the extraction methods, results showed that DO, A260/A280, DNA YIELD, EX2, EX7 and EX6 (right hand side) and EX1, EX3, EX4 and EX5 (left hand), contribute the most to the significant variation on PC1. Results showed that EX2, EX7 and EX6 were positively correlated with PCR amplification, while EX1, EX3, EX4 and EX5 were negatively correlated with results of PCR. Such data confirmed results previously recorded for the DNA yield/purity and amplicon observed on agarose gel.
The quality of DNA, with the highest amplicon length obtained from the fresh and the different tuna food product, was best preserved using the SDS/PCI precipitation method, with the advantage of being the fastest and simplest protocol and is less expensive than the commercial kit.
The relatively high values of the PLS of all the data (45%) suggest substantial difference between the extraction procedures and matrices when all the studied parameters were compared. Such statistical data confirm results previously recorded. Thus, the SDS/PCI method can be a good choice for processed tuna products.
Based on the result obtained from this study, detection of cyt b gene of species included in a particular processed food may largely dependent on the processing technique applied. The extraction of DNA that is of a high quantity and quality is generally recognised as the most important factor for the success of molecular techniques, such as PCR and DNA sequencing. Overall, this study identified the PCI/SDS method as the most suitable method for the extraction of DNA from tuna, due to the high DNA yields delivered the ease of use and is less expensive than the commercial kit and the applicability for high throughput extraction from multiple specimens.
The partial least squares regression models (PLS) suggest substantial difference between the extraction procedures and matrices when all the studied parameters were compared. Thus, the SDS/PCI method can be a good choice for highly processed foodstuffs.
Finally, tuna samples were also analysed by principal component analysis (PCA) to examine the influence of sample matrix (type of tuna) and DNA extraction method on the PCR result, DNA yield and DO. Variables were scaled by column so they have a mean of zero and a variance of one to compensate for their different scales (Figure 4). Three grouping variables were used to analyse the results: type  3. Gráficos de carga de la regresión de mínimos cuadrados parciales realizada en todas las variables X DO: A260/A280, rendimiento de ADN, EX1: Método de ADN rápido del kit, EX2: Método PCI/SDS, EX3: Método SDS/PCI/NaCl, EX4: Método CHELEX, EX5: Método CTAB/CI, EX6: Método TRITON/CI, EX7: Método del kit DNeasy Tissue y todos los alimentos del atún y la variable Y (resultado de la RCP: amplificación de la RCP). Todos los datos de las proyecciones del PCA representan el valor medio de seis réplicas para cada parámetro.  of DNA extraction method (1-7), product type (fresh, frozen and canned) and processing (none, olive oil and vegetable oil). The tuna sample type was used as Y-matrix; all the other parameters (DO A260/A280, DNA YIELD and PCR RESULTS) as X-variables. In the scores and loadings plot (Figure 4), the explained variance of the two first principal components was of 82.1%, of which 51.6% was along PC1 and 30.5% along PC2. Such result highlighted that extraction method 2 yielded the higher DNA yield in all cases (upper right region). According to the product type, fresh tuna samples produce better DO and PCR results than frozen or canned tuna. Comparing all products, the fresh tuna yielded the highest amounts of DNA independently of the extraction methods; while matrices with high degree of processing such as canned tuna, yielded the lowest amounts of DNA.

Conclusions
A label indicating authentication of imported products aimed at benefiting producers from developing countries will surely be well-accepted and demanded by an increasing sector of consumers. Additionally, it would encourage traders to avoid fraudulent commercialisation of products from developing countries, even if they have to pay for authentication certificates. Genetic tools for application in food fish traceability are abundant and affordable even for modest budgets, and now it is time to enforce their use in routine controls of international tuna trade. In this study, the processed tuna fish found in 53 samples was confirmed by DNA barcoding method to identify its authenticity. The amplification of extracted DNA was successfully performed. Considering the present results, initiation of method harmonisation and blockchain application should be reinforced to control international tuna trade.