Assessment of Genetic Diversity, Gene Flow and Demographic History of Frigate Tuna (Auxis thazard) Populations in Tanzanian Marine Waters Using Mitochondrial DNA Control Region

Aim: To investigate genetic diversity, gene flow and demographic history of frigate tuna                   (Auxis thazard) populations found in Tanzanian marine waters. Study Design: The study used a descriptive research design whereby fish samples were randomly collected from four locations and the genetic variation within and among the four populations analyzed using mitochondrial DNA (mtDNA) control region. Place and Duration of the Study: Fish samples were collected from landing sites in Tanga, Dar es Salaam, Mtwara and Zanzibar, Tanzania. The study was conducted between July 2020 and June 2021. Methodology: A total of 100 frigate tuna were randomly sampled from small-scale fishermen at the landing sites of Dar-es-Salaam (20), Tanga (30), Mtwara (30), and Zanzibar (20). For each fish, 50 g muscle tissue was obtained and put in a vial containing 95% ethanol. DNA was extracted from the muscle using a commercial DNA Kit (Quick-DNATM Miniprep Plus Kit, Zymo Research Corp.) according to the instructions of the manufacturer and a fragment of 432 bp of the mtDNA control region was amplified and sequenced using ABI PRISIM™ 3100 Genetic Analyzer (Applied Biosystems). Haplotype and nucleotide diversity, gene flow and historic demographic were estimated from 92 fish samples. Results: A total of 88 haplotypes were identified in all fish samples. The highest haplotype diversity was found in Zanzibar (1.000 ± 0.017) and Mtwara (1.000 ± 0.010) populations while the lowest was observed in Tanga population (0.992 ± 0.012). Tanga population had the highest nucleotide diversity (0.078 ± 0.018) while Dar es Salaam had the lowest (0.016 ± 0.009). The highest genetic differentiation (FST) was found between Tanga and Dar-es-Salaam (0.178) and the lowest was observed between Mtwara and Zanzibar (0.016). The genetic distances between pairs of populations were small. The phylogenetic tree revealed two main clusters. Cluster 1 which consisted of nine haplotypes was dominated by Tanga population with seven haplotypes while the remaining two haplotypes were from Mtwara and the reference sequences of Euthynnus affinis (Kawakawa tuna). Cluster II had 84 haplotypes of individuals from all four populations, with no population specific subcluster. The number of immigrants per generation was highest between Mtwara and Zanzibar (Nm=18.310) and lowest between Tanga and Dar-es-Salaam (Nm=1.180). The neutrality test indicated negative values, suggesting a recent population expansion. Conclusion: There is high genetic diversity within the populations, but there is no significant genetic differentiation among the four frigate tuna populations, suggesting that the four populations comprise a single panmictic population.

[1]  Zhe Liu,et al.  Genetic diversity of Gymnocypris chilianensis (Cypriniformes, Cyprinidae) unveiled by the mitochondrial DNA D-loop region , 2021, Mitochondrial DNA. Part B, Resources.

[2]  Y. Sung,et al.  Recent population expansion of longtail tuna Thunnus tonggol (Bleeker, 1851) inferred from the mitochondrial DNA markers , 2020, PeerJ.

[3]  Jianguo Du,et al.  Genetic diversity and population structure of Terapon jarbua (Forskål, 1775) (Teleostei, Terapontidae) in Malaysian waters , 2020, ZooKeys.

[4]  R. Jehle,et al.  High genetic diversity and lack of pronounced population structure in five species of sympatric Pacific eels. , 2018, Fisheries management and ecology.

[5]  Min-Yun Liu,et al.  Genetic Stock Structure of Terapon jarbua in Taiwanese Waters , 2015 .

[6]  W. Hong,et al.  Population genetic structure and demographic history of the mudskipper Boleophthalmus pectinirostris on the northwestern pacific coast , 2015, Environmental Biology of Fishes.

[7]  P. M. Galetti,et al.  Genetic Multipartitions Based on D-Loop Sequences and Chromosomal Patterns in Brown Chromis, Chromis multilineata (Pomacentridae), in the Western Atlantic , 2014, BioMed research international.

[8]  Koichiro Tamura,et al.  MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. , 2013, Molecular biology and evolution.

[9]  R. Meena,et al.  Mitochondrial DNA analysis reveals three stocks of yellowfin tuna Thunnus albacares (Bonnaterre, 1788) in Indian waters , 2013, Conservation Genetics.

[10]  R. Meena,et al.  Three genetic stocks of frigate tuna Auxis thazard thazard (Lacepede, 1800) along the Indian coast revealed from sequence analyses of mitochondrial DNA D-loop region , 2012 .

[11]  E. Lavergne Estuarine fish biodiversity of Socotra Island (N.W. Indian Ocean): from the fish community to the functioning of Terapon jarbua populations , 2012 .

[12]  L. Excoffier,et al.  Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows , 2010, Molecular ecology resources.

[13]  N. Galtier,et al.  Mitochondrial DNA as a marker of molecular diversity: a reappraisal , 2009, Molecular ecology.

[14]  Pablo Librado,et al.  DnaSP v5: a software for comprehensive analysis of DNA polymorphism data , 2009, Bioinform..

[15]  E. González,et al.  Genetic structuring and migration patterns of Atlantic bigeye tuna, Thunnus obesus (Lowe, 1839) , 2008, BMC Evolutionary Biology.

[16]  D. Hurwood,et al.  Evidence for fine geographical scale heterogeneity in gene frequencies in yellowfin tuna (Thunnus albacares) from the north Indian Ocean around Sri Lanka , 2008 .

[17]  Rodrigo Lopez,et al.  Clustal W and Clustal X version 2.0 , 2007, Bioinform..

[18]  N. Taniguchi,et al.  Genetic variation in skipjack tuna Katsuwonus pelamis(L.) using PCR-RFLP analysis of the mitochondrial DNA D-loop region , 2006 .

[19]  D. R. Robertson,et al.  PHYLOGEOGRAPHY OF THE TRUMPETFISHES (AULOSTOMUS): RING SPECIES COMPLEX ON A GLOBAL SCALE , 2001, Evolution; international journal of organic evolution.

[20]  Y. Fu,et al.  Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. , 1997, Genetics.

[21]  M. Slatkin,et al.  Estimation of levels of gene flow from DNA sequence data. , 1992, Genetics.

[22]  F. Tajima Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. , 1989, Genetics.

[23]  M Slatkin,et al.  Gene flow and the geographic structure of natural populations. , 1987, Science.

[24]  S. Wright THE INTERPRETATION OF POPULATION STRUCTURE BY F‐STATISTICS WITH SPECIAL REGARD TO SYSTEMS OF MATING , 1965 .

[25]  J. Maguire The state of world highly migratory, straddling and other high seas fishery resources and associated species , 2006 .

[26]  L. Lim Introduction to Conservation Genetics , 2001 .