Next Article in Journal
Bacterial Diversity and Community in Response to Long-Term Nitrogen Fertilization Gradient in Citrus Orchard Soils
Next Article in Special Issue
Revisiting the Diversity of Barbonymus (Cypriniformes, Cyprinidae) in Sundaland Using DNA-Based Species Delimitation Methods
Previous Article in Journal
Distribution and Potential Limiting Factors of the European Pond Turtle (Emys orbicularis) in Eastern Europe
Previous Article in Special Issue
Phenetic and Genetic Variability of Continental and Island Populations of the Freshwater Copepod Mastigodiaptomus ha Cervantes, 2020 (Copepoda): A Case of Dispersal?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 13
Issue 7
10.3390/d13070281
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mitochondrial Genetic Diversity among Farmed Stocks of Oreochromis spp. (Perciformes, Cichlidae) in Madagascar

by
Nicolas Hubert
1,*,
Elodie Pepey
1,2,
Jean-Michel Mortillaro
1,2,3,
Dirk Steinke
4,
Diana Edithe Andria-Mananjara
3 and
Hugues de Verdal
1,2
1
UMR ISEM (IRD, UM, CNRS), Université de Montpellier, Place Eugène Bataillon, CEDEX 05, 34095 Montpellier, France
2
Cirad, UMR ISEM, Université de Montpellier, 389 Avenue Agropolis, CEDEX 05, 34095 Montpellier, France
3
FOFIFA DRZVP, rue Farafaty, Antananarivo 00101, Madagascar
4
Department of Integrative Biology, Centre for Biodiversity Genomics, University of Guelph, 50 Stone Rd E, Guelph, ON N1G 2W1, Canada
*
Author to whom correspondence should be addressed.
Diversity 2021, 13(7), 281; https://doi.org/10.3390/d13070281
Submission received: 12 May 2021 / Revised: 15 June 2021 / Accepted: 16 June 2021 / Published: 22 June 2021
(This article belongs to the Special Issue Aquatic Organisms Research with DNA Barcodes)
Browse Figures
Versions Notes

Abstract

:
The fast development of aquaculture over the past decades has made it the main source of fish protein and led to its integration into the global food system. Mostly originating from inland production systems, aquaculture has emerged as strategy to decrease malnutrition in low-income countries. The Nile tilapia (Oreochromis niloticus) was introduced to Madagascar in the 1950s, and is now produced nationally at various scales. Aquaculture mostly relies on fry harvested from wild populations and grow-out in ponds for decades. It has recently been diversified by the introduction of several fast-growing strains. Little is known how local genetic diversity compares to recently introduced strains, although high and comparable levels of genetic diversity have previously been observed for both wild populations and local stocks. Our study compares DNA barcode genetic diversity among eight farms and several strains belonging to three species sampled. DNA-based lineage delimitation methods were applied and resulted in the detection of six well differentiated and highly divergent lineages. A comparison of DNA barcode records to sequences on the Barcode of Life Data System (BOLD) helped to trace the origin of several of them. Both haplotype and nucleotide diversity indices highlight high levels of mitochondrial genetic diversity, with several local strains displaying higher diversity than recently introduced strains. This allows for multiple options to maintain high levels of genetic diversity in broodstock and provides more options for selective breeding programs.

1. Introduction

In the light of a global biodiversity decline and the wide-spread depletion of wild fish stocks [1,2,3], aquaculture has become an increasingly important source of animal protein in tropical countries [4,5]. With a global production of more than 160 million tons of fish (in 2015) aquaculture has become an important component of the global food system [5,6], with Asia contributing nearly 90% of the freshwater production. Globally it is now producing nearly twice as much as fisheries [4], mostly through inland aquaculture. It has been frequently integrated in national strategies to prevent malnutrition, particularly in low-income, tropical countries [7].
Inland aquaculture mostly relies on three species or species groups: carp (Cyprinus carpio) which is the most farmed fish worldwide with 13 million tons produced in 2017, followed by tilapia (Oreochromis spp.) and catfishes, such as Pangasius spp., with 5 million tons each produced in 2017 [5]. Their success lies mostly in the fast growth rate of these species and the availability of protocols for farming and captive breeding. As a consequence, they were introduced in numerous countries outside their native distribution ranges [8].
In Madagascar, freshwater fish aquaculture has a long tradition. The first introductions of exotic species for aquaculture trace back to 1857 (giant goramy-Osphronemus goramy), 1861 (goldfish-Carassius auratus), 1912 (carp-Cyprinus carpio) [9], and the 1950s (nile tilapia-Oreochromis niloticus) [10,11]. The successful introduction of these species led to the establishment of multiple populations in the wild, further used as source of young fry for farming [12]. Nile tilapia and carp have been used in semi-intensive productions in association with rice cultivation in paddy fields since the earliest development of aquaculture in Madagascar [13]. Adults and fry are traditionally caught in the wild and introduced to paddy fields when rice plants are sown [12]. Aquaculture in ponds follows similar practices leading to similar genetic diversity found between farmed and wild tilapia populations [14]. More recently, multiple strains of O. niloticus and other Oreochromis species were introduced to upscale tilapia aquaculture by producing fry from imported strains known to have desirable properties, such as fast growth rates [15]. Some of these strains, such as the Genetically Improved Farmed Tilapia (GIFT) strain, are the result of multiple crossings between several farmed strains and wild populations, as well as selective breeding [16,17]. As a result, strains of multiple origins are currently maintained in farmed stocks of tilapia in Madagascar, but their genetic variability is not known.
The present study aims to estimate the genetic diversity of currently farmed strains of Oreochromis in Madagascar, which now include local and recently imported strains of O. niloticus, such as GIFT and JICA, as well as O. mossambicus and O. macrochir [14]. We opted for the use of a standardized mitochondrial marker, a 652 bp fragment of the cytochrome oxidase I gene known as DNA barcode [18,19], for a variety of reasons: (1) it can be easily retrieved thanks to the availability of universal primers for fish [20]; (2) a large number of DNA barcode sequences are available for Oreochromis (>1200 public records for most species) in the Barcode of Life Data System, BOLD [21]; (3) the mitochondrial diversity of farmed tilapia Madagascar is largely unknown.
We generated DNA barcodes for currently farmed strains of tilapia, and applied DNA-based delimitation methods to detect mitochondrial lineages among them and to estimate their diversity. By using new DNA barcode records and published data from BOLD, potential evolutionary origins of revealed lineages are discussed.

2. Materials and Methods

2.1. Sampling, Sequencing and International Repositories

In February 2016, a total of 262 specimens was collected at eight sites maintaining breeding stocks (Figure 1). Collected information, including geocoordinates, are included in the dataset DS-TILMADA (dx.doi.org/10.5883/DS-TILMADA) on BOLD. Specimens were identified to the species level based on the information of the farmers. Specimens were further photographed, individually labeled, and voucher specimens were preserved in a 5% formalin solution. A fin clip or a muscle biopsy were taken from each specimen and fixed in a 96% ethanol solution for further genetic analyses.
Genomic DNA was extracted from fin clip samples using a Qiagen DNeasy 96 tissue extraction kit following manufacturer’s specifications. A 652-bp segment from the 5′ region of the cytochrome oxidase I gene (COI) was amplified using the primer cocktail C_FishF1t1/C_FishR1t1 [20]. PCR amplifications were done on a Veriti 96-well Fast thermocycler (ABI-AppliedBiosystems) with a final volume of 10.0 μL containing 5.0 μL Buffer 2X, 3.3 μL ultrapure water, 1.0 μL each primer (10 μM), 0.2 μL enzyme Phire Hot Start II DNA polymerase (5 U), and 0.5 μL of DNA template (~50 ng). The following thermocycler regime was used: initial denaturation at 98 °C for 5 min followed by 30 cycles denaturation at 98 °C for 5 s, annealing at 56 °C for 20 s, and extension at 72 °C for 30 s, followed by a final extension step at 72 °C for 5 min. PCR products were purified with ExoSap-IT (USB Corporation, Cleveland, OH, USA) and sequenced in both directions. Sequencing reactions were performed at the Centre for Biodiversity Genomics, University of Guelph, Canada, using the BigDye Terminator v3.1 Cycle Sequencing Ready Reaction kit following standard protocols. Sequencing was performed on an ABI 3730xl capillary sequencer (Applied Biosystems). Sequences and collateral information were deposited on BOLD [21], and are available as a public dataset (dx.doi.org/10.5883/DS-TILMADA, Table S1).

2.2. Mitochondrial Lineage Delimitation and Genetic Diversity

Several methods for species delineation based on DNA sequences have been proposed [22,23,24,25,26]. Each of these have different properties, particularly when dealing with singletons (i.e., lineages represented by a single sequence) or heterogeneous speciation rates among lineages [27]. A combination of different approaches is increasingly used to overcome potential pitfalls arising from uneven sampling [28,29,30,31,32]. We used six different sequence-based methods of species delimitation to identify Molecular Operational Taxonomic Units (MOTU): (1) Refined Single Linkage (RESL) as implemented in BOLD and used to generate Barcode Index Numbers (BIN) [25]; (2) Automatic Barcode Gap Discovery (ABGD) [24]; (3) Poisson Tree Process (PTP) in its single (sPTP) and multiple rates version (mPTP) as implemented in the stand-alone software mptp_0.2.3 [26,33]; (4) General Mixed Yule-Coalescent (GMYC) in its single (sGMYC) and multiple threshold version (mGMYC) as implemented in the R package Splits 1.0–19 [34].
The mPTP algorithm and the GMYC both use phylogenetic trees as input file. We reconstructed a maximum likelihood (ML) tree for the former using RAxML [35] with a GTR + I + Γ substitution model. For the GMYC algorithm we calculated an ultrametric, fully resolved the tree using the Bayesian approach implemented in BEAST 2.6.2 [36]. Sequences were collapsed into haplotypes prior to reconstructing the ultrametric tree using RAxML, and Bayesian reconstruction was based on a strict-clock prior of 1.2% per million year [37]. Two Markov chains of 20 million each were ran independently using Yule pure birth and GTR + I + Γ substitution models, other tree priors were used as default. Stability (ESS > 200) and convergence of Markov chains was verified using Tracer 1.7.1 [36]. Trees were sampled every 5000 states, after an initial burn in period of 5 million states. Both runs were combined with trees re-sampled every 20,000 states using LogCombiner 2.6.2, and the maximum credibility tree was constructed using TreeAnnotator 2.6.2 [36].
A final COI gene tree was reconstructed using the SpeciesTreeUCLN algorithm of the StarBEAST2 package [38]. This approach implements a mixed-model including a coalescent component within species and a diversification component between species that allows accounting for variations of substitution rates within and between species [39]. SpeciesTreeUCLN jointly reconstructs gene trees and species trees, and, as such, requires the designation of species, which were determined using the consensus of our species delimitation analyses. The SpeciesTreeUCLN analysis was performed with the same parameters as mentioned above.
Several parameters of genetic diversity were estimated using the R package pegas 1.0 [40], including the number of haplotypes (h), haplotype diversity (Hd) [41], nucleotide diversity (π) [42], genetic diversity based on the number of segregating sites (θ) [43], and Tajima’s D test of neutrality [44]. Kimura 2-parameter (K2P) [45] pairwise genetic distances were calculated using the R package Ape 5.4 [46]. Maximum intraspecific and nearest neighbor genetic distances were calculated from the pairwise K2P distance matrix using the R package Spider 1.5 [47].

3. Results

A total of 263 COI sequences were generated from six sites for O. niloticus, and one site for O. mossambicus and O. macrochir. Sequences consisted of mostly full-length sequences (652 bp) and no stop codons were detected, suggesting that the sequences collected represent functional coding regions. DNA-based species delimitation methods resulted in congruent delimitation schemes with 6 MOTUs for ABGD, RESL, and sPTP, 5 MOTUs for mPTP and sGMYC, and 8 MOTUs for mGMYC (Figure 2, Table S1). The final consensus consisted of 6 MOTUs (Figure 2). Maximum distances within MOTU ranged from 0 (BOLD:AAA8513, BOLD:ACR7163, BOLD:ADI0792) to 0.007 (BOLD:AAC9904), and minimum distances between MOTUs ranged from 0.016 (BOLD: AAA8513, BOLD: ADI0792) to 0.065 for BOLD:AAA6537 (Table 1). K2P distances were 7-fold higher on average between than within MOTUs. None of the MOTUs were restricted to a single strain, excepting BOLD:AAA8513 restricted to O. macrochir (Figure 2). The five remaining MOTUs were shared between strains of O. niloticus and O. mossambicus (Figure 2, Table 1).
Shared MOTUs among species of Oreochromis were found on BOLD for several BINs (Table 1). Records associated to: (1) BOLD:AAA6537 members mostly belong to O. niloticus; (2) BOLD:AAA8511 members mostly belong to O. niloticus and O. mossambicus; (3) BOLD:AAA8513 members belong to O. macrochir; (4) BOLD:AAC9904 members mostly belong to O. niloticus; (5) BOLD:ACR7163 members mostly belong to O. urolepis; and (6) BOLD:ADI0792 members mostly belong to O. mossambicus. A substantial proportion of Oreochromis, BOLD records have not been identified to species particularly within the BINs BOLD:AAA6537 and BOLD:AAA8513, resulting in a low frequency of sequences named. Two discrepancies in the proportion of species per MOTU were detected between BOLD records and sequences generated for the present study. BOLD:ACR7163 mostly contains records assigned to O. urolepis in BOLD although this species was never reported from Madagascar. BOLD:ADI0792, mostly represented by O. niloticus, is a new occurrence of this lineage for this species.
As a result of the presence of multiple mitochondrial lineages in breeding stocks, which are estimated to have diverged between 0.6 and 2.7 million years ago (Figure 2), the nucleic diversity in most stocks is high, with π and θw above 0.015 and 8, respectively (Table 2). However, the JICA stock of O. niloticus (Site 4, Table 2) and O. macrochir at Milasoa (Site 8, Table 2) display much lower nucleic diversity. The high diversity estimates of nucleic diversity are accompanied by high haplotype diversity at sites 1, 2, 3, 6, and 7, due to the occurrence of a high number of haplotypes, peaking at 8 in site 2 (Table 2). None of the Tajima’s D tests were significant, indicating that diversity was mostly balanced for the eight breeding stocks, despite some negative values at site 2 (O. niloticus GIFT at Tamatave) and site 4 (O. niloticus JICA at Mahajunga) indicative of a deficit of rare haplotypes.

4. Discussion

The present study highlights high levels of mitochondrial genetic diversity, both at nucleotide and haplotype level, across most Oreochromis stocks assessed. Such high levels of genetic diversity originate from the co-occurrence of multiple, highly divergent MOTUs within most stocks, for both O. niloticus and O. mossambicus. However, high haplotype diversity also indicates the presence of substantial genetic diversity within each MOTU. These likely have multiple origins. In their native range distribution, Oreochromis species are known to easily hybridize after secondary contacts and can result in discordant evolutionary histories between genomes [48]. Multiple cases of introgression of both mitochondrial and nuclear genomes have been reported [49,50,51,52,53]. This likely accounts for the evolutionary success of the group. Hybrids with higher fitness than parental lineages, living outside of the ecological range of their parental species, have been detected in cichlids [54,55] and other perciformes [56]. In the context of fish farming, mitochondrial genomes resulting from ancient introgression events are likely introduced during the building of brood stock or through introductions into the wild. The occurrence of BOLD:ACR7163 and BOLD:AAA6537, mostly assigned to O. urolepis and O. aureus records on BOLD, respectively, seems to support the assumption that heterospecific mitochondrial genomes resulting from ancient introgression events were inadvertently introduced into Madagascar stocks of O. niloticus, as O. urolepis and O. aureus have never been reported as introduced species for Madagascar [57]. Alternatively, inadvertent and unrecorded introductions of O. urolepis and O. aureus individuals within imported batches of O. niloticus cannot be discarded, particularly if young, immature, and morphologically indistinguishable individuals were introduced.
The high mitochondrial diversity and multiple occurrences of MOTUs among species and strains might further reflect local breeding strategies and stock maintenance. Broodstock are known to display similar levels of genetic diversity in comparison with wild populations, which is the result of introductions in the 1950s [14]. Our study suggests tilapia production in Madagascar relies mostly on natural populations, with fry being collected in rivers and grow-out occurrence in ponds. Alternatively, fry are produced in ponds and broodstock genetic diversity is maintained by the regular addition of wild adults [14,58]. Relationships between broodstocks and wild populations depend on the scale of the fish farms. Oswald et al. [58] identified at least three type of fish farms in Madagascar: (1) small-scale, artisanal farms consisting of a single pond for fry production and grow-out with production mostly devoted to local consumption; (2) medium-scale farms where fry production and grow-out are managed in different ponds, and fry production is partly disseminated to small-scale farms for grow-out; and (3) large-scale farms devoted to fry production and distribution, constituting the main source of fry for aquaculture. Fry distribution can occur over large distances within Madagascar, and most small farms have no particular strategy to maintain strains. The similar mitochondrial composition of sampled stocks, and multiple occurrences of MOTUs in O. niloticus and O. mossambicus strains confirms that stocks used to be mixed in the past, and further contributed to increase strain mitochondrial genetic diversity. However, management of the most recently introduced strains seems to depart from this trend as the O. macrochir strain hosts a single and private MOTU. Observed haplotype and nucleotide diversity among imported (GIFT, JICA) and local tilapia strains are mostly similar, however, several local strains display higher mitochondrial genetic diversity. This trend suggests that the initial genetic pool of O. niloticus introduced in the 1950s was already diverse. Its successful introduction likely accounts for the persistence of high mitochondrial genetic diversity in farmed strains, due to their historical reliance on wild populations. Several studies have shown a relationship between individual genetic diversity and fitness in fishes [59,60,61]. As such, this high mitochondrial genetic diversity, also previously reported using nuclear markers [14], suggests local strains constitute good candidates for selective breeding, with genetically diverse wild populations available for maintaining high levels of genetic diversity.

5. Conclusions

The present study shows that tilapia stocks in Madagascar have genetically diverse mitochondrial genomes, likely resulting from an intricate history of ancient introgressions, introduction of genetically diverse populations in the 1950s, and local practices. The detection of higher mitochondrial genetic diversity in several local strains in comparison to introduced strains of O. niloticus, suggests local practices are beneficial in maintaining a high level of genetic diversity. In the context of fast development of modern aquaculture, this high genetic diversity, also previously observed for nuclear markers [14], certainly constitutes an asset. With genetically diverse wild populations, multiple options are available for maintaining high levels of genetic diversity in broodstock and several strategies are further available for selective breeding programs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/d13070281/s1. Table S1: Results of the genetic species delimitation, including the majority-rule consensus.

Author Contributions

Conceptualization, E.P., H.d.V. and N.H.; methodology, E.P., H.d.V. and N.H.; software, E.P., H.d.V. and N.H.; validation, J.-M.M., E.P., D.S., D.E.A.-M., H.d.V. and N.H.; formal analysis, N.H.; investigation, J.-M.M., E.P., D.S., D.E.A.-M., H.d.V. and N.H.; resources, E.P., D.E.A.-M., H.d.V., D.S. and N.H.; data curation, E.P. and N.H.; writing—original draft preparation, N.H., E.P., J.-M.M. and H.d.V.; writing—review and editing, D.E.A.-M. and D.S.; visualization, J.-M.M., E.P., D.S., D.E.A.-M., H.d.V. and N.H.; supervision, E.P. and N.H.; project administration, E.P.; funding acquisition, E.P., D.S. and N.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the CFREF to the University of Guelph’s Food from Thought program, by the “Institut de Recherche pour le Développement” through annual funding and “Fonds d’Amorçage” (226F2ABIOS) and through the AMPIANA project funded by the European Union (EuropeAid/135-182/DD).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in BOLD at [dx.doi.org/10.5883/DS-TILMADA], reference number DS-TILMADA.

Acknowledgments

Authors are most grateful to Lionel Dabbadie (Cirad) and Tojoharivelo Rakotomalala (FOFIFA) for the identification of the different sites studied and for organizing field sampling. The authors are also grateful to Marc Canonne (Cirad), the NGO APDRA Pisciculture Paysanne and FOFIFA for technical and logistic assistance before or during field work. Finally, the authors are very grateful to André Ramilison and Zarre, fish farmers, who agreed to make their fish available for this study and the staff of the Ivoloina fish farming station. This is publication ISEM 2021-134 SUD.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Worm, B.; Barbier, E.B.; Beaumont, N.; Duffy, E.; Folke, C.; Halpern, B.S.; Jackson, J.B.C.; Lotze, H.K.; Micheli, F.; Palumbi, S.R.; et al. Impacts of biodiversity loss on ocean ecosystem services. Science 2006, 314, 787–790. [Google Scholar] [CrossRef] [Green Version]
  2. Newbold, T.; Hudson, L.N.; Arnell, A.P.; Contu, S.; De Palma, A.; Ferrier, S.; Hill, S.L.L.; Hoskins, A.J.; Lysenko, I.; Phillips, H.R.P. Has land use pushed terrestrial biodiversity beyond the planetary boundary? A global assessment. Science 2016, 353, 288–291. [Google Scholar] [CrossRef]
  3. Thiault, L.; Mora, C.; Cinner, J.E.; Cheung, W.W.L.; Graham, N.A.J.; Januchowski-Hartley, F.A.; Mouillot, D.; Sumaila, U.R.; Claudet, J. Escaping the perfect storm of simultaneous climate change impacts on agriculture and marine fisheries. Sci. Adv. 2019, 5, eaaw9976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. FAO. Des Pêches Et De L’aquaculture; FAO: Rome, Italy, 2018; ISBN 9789251306925. [Google Scholar]
  5. Naylor, R.L.; Hardy, R.W.; Buschmann, A.H.; Bush, S.R.; Cao, L.; Klinger, D.H.; Little, D.C.; Lubchenco, J.; Shumway, S.E.; Troell, M. A 20-year retrospective review of global aquaculture. Nature 2021, 591, 551–563. [Google Scholar] [CrossRef] [PubMed]
  6. Tacon, A.G.J. Trends in global aquaculture and aquafeed production: 2000–2017. Rev. Fish. Sci. Aquac. 2020, 28, 43–56. [Google Scholar] [CrossRef]
  7. Lobell, D.B.; Burke, M.B.; Tebaldi, C.; Mastrandrea, M.D.; Falcon, W.P.; Naylor, R.L. Prioritizing Climate Change Adaptation Needs for Food Security in 2030. Science 2008, 319, 607–610. [Google Scholar] [CrossRef] [PubMed]
  8. Froese, R.; Pauly, D. Fishbase. Available online: http://www.fishbase.org (accessed on 14 January 2021).
  9. Kiener, A. Intérêt et perspectives de la pisciculture de la carpe à Madagascar. Bull. Madag. 1958, 8, 693–702. [Google Scholar]
  10. Kiener, A. Poissons, pêche et pisciculture à Madagascar. Centrotechique For. Trop. Fr. 1963, 405, 599–606. [Google Scholar]
  11. Therezien, Y. L’introduction de poissons d’eau douce à Madagascar, leur influence sur la modification du biotope. Bull. Français Piscic. 1960, 199, 45–61. [Google Scholar] [CrossRef]
  12. Duchaufour, H.; Razafimbelo-Andriamifidy, T.; Rakotoarisoa, J.; Ramamonjisoa, B. ; Rakotondravao. Recherche Interdisciplinaire pour le Développement Durable Application à Différentes Thématiques de Territoire et la Biodiversité des Espaces Ruraux Malgaches; Cirad: Antananarivo, Madagascar, 2016. [Google Scholar]
  13. Bentz, B.; Oswald, M. Respective roles of national institutions and farmers groups in the implementation of an innovation enabling smallholders to reproduce Carp inside their rice fields in betafo (Madagascar). In Proceedings of the Symposium Innovation and Sustainable Development in Agriculture and Food, ISDA 2010, Montpellier, France, 28 June 2010. [Google Scholar]
  14. Ravakarivelo, M.; Pepey, E.; Benzie, J.; Raminosoa, N.; Rasamoelina, H.; Mikolasek, O.; De Verdal, H. Genetic variation in wild populations and farmed stocks of Nile tilapia (Oreochromis niloticus) in Madagascar. Revue d’Elevage et de Médecine Vétérinaire des Pays Tropicaux 2019, 72, 101. [Google Scholar] [CrossRef]
  15. Ponzoni, R.W.; Hamzah, A.; Tan, S.; Kamaruzzaman, N. Genetic parameters and response to selection for live weight in the GIFT strain of Nile tilapia (Oreochromis niloticus). Aquaculture 2005, 247, 203–210. [Google Scholar] [CrossRef]
  16. Ponzoni, R.W.; Nguyen, N.H.; Khaw, H.L.; Hamzah, A.; Bakar, K.R.A.; Yee, H.Y. Genetic improvement of Nile tilapia (Oreochromis niloticus) with special reference to the work conducted by the WorldFish Center with the GIFT strain. Rev. Aquac. 2011, 3, 27–41. [Google Scholar] [CrossRef]
  17. Gupta, M.V.; Acosta, B.O. From drawing board to dining table: The success story of the GIFT project. NAGA WorldFish Cent. Q. 2004, 27, 4–14. [Google Scholar]
  18. Ward, R.D.; Hanner, R.H.; Hebert, P.D.N. The campaign to DNA barcode all fishes, FISH-BOL. J. Fish Biol. 2009, 74, 329–356. [Google Scholar] [CrossRef]
  19. Hubert, N.; Hanner, R.; Holm, E.; Mandrak, N.E.; Taylor, E.; Burridge, M.; Watkinson, D.; Dumont, P.; Curry, A.; Bentzen, P.; et al. Identifying Canadian freshwater fishes through DNA barcodes. PLoS ONE 2008, 3. [Google Scholar] [CrossRef] [Green Version]
  20. Ivanova, N.V.; Zemlak, T.S.; Hanner, R.H.; Hébert, P.D.N. Universal primers cocktails for fish DNA barcoding. Mol. Ecol. Notes 2007, 7, 544–548. [Google Scholar] [CrossRef]
  21. Ratnasingham, S.; Hebert, P.D.N. BOLD: The Barcode of Life Data System (www.barcodinglife.org). Mol. Ecol. Notes 2007, 7, 355–364. [Google Scholar] [CrossRef] [Green Version]
  22. Hebert, P.D.N.; deWaard, J.R.; Zakharov, E.; Prosser, S.W.J.; Sones, J.E.; McKeown, J.T.A.; Mantle, B.; La Salle, J. A DNA “barcode blitz”: Rapid digitization and sequencing of a natural history collection. PLoS ONE 2013, 8, e68535. [Google Scholar] [CrossRef]
  23. Pons, J.; Barraclough, T.G.; Gomez-Zurita, J.; Cardoso, A.; Duran, D.P.; Hazell, S.; Kamoun, S.; Sumlin, W.D.; Vogler, A.P. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst. Biol. 2006, 55, 595–606. [Google Scholar] [CrossRef] [Green Version]
  24. Puillandre, N.; Lambert, A.; Brouillet, S.; Achaz, G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 2012, 21, 1864–1877. [Google Scholar] [CrossRef]
  25. Ratnasingham, S.; Hebert, P.D.N. A DNA-based registry for all animal species: The barcode index number (BIN) system. PLoS ONE 2013, 8, e66213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kapli, P.; Lutteropp, S.; Zhang, J.; Kobert, K.; Pavlidis, P.; Stamatakis, A.; Flouri, T. Multi-rate Poisson Tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo. Bioinformatics 2017, 33, 1630–1638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Luo, A.; Ling, C.; Ho, S.Y.W.; Zhu, C.-D. Comparison of methods for molecular species delimitation across a range of speciation scenarios. Syst. Biol. 2018, 67, 830–846. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kekkonen, M.; Mutanen, M.; Kaila, L.; Nieminen, M.; Hebert, P.D.N. Delineating species with DNA barcodes: A case of taxon dependent method performance in moths. PLoS ONE 2015, 10, e0122481. [Google Scholar] [CrossRef] [Green Version]
  29. Kekkonen, M.; Hebert, P.D.N. DNA barcode-based delineation of putative species: Efficient start for taxonomic workflows. Mol. Ecol. Resour. 2014, 14, 706–715. [Google Scholar] [CrossRef]
  30. Shen, Y.; Hubert, N.; Huang, Y.; Wang, X.; Gan, X.; Peng, Z.; He, S. DNA barcoding the ichthyofauna of the Yangtze River: Insights from the molecular inventory of a mega-diverse temperate fauna. Mol. Ecol. Resour. 2019, 19, 1278–1291. [Google Scholar] [CrossRef]
  31. Sholihah, A.; Delrieu-Trottin, E.; Sukmono, T.; Dahruddin, H.; Risdawati, R.; Elvyra, R.; Wibowo, A.; Kustiati, K.; Busson, F.; Sauri, S.; et al. Disentangling the taxonomy of the subfamily Rasborinae (Cypriniformes, Danionidae) in Sundaland using DNA barcodes. Sci. Rep. 2020, 10, 2818. [Google Scholar] [CrossRef] [Green Version]
  32. Limmon, G.; Delrieu-Trottin, E.; Patikawa, J.; Rijoly, F.; Dahruddin, H.; Busson, F.; Steinke, D.; Hubert, N. Assessing species diversity of Coral Triangle artisanal fisheries: A DNA barcode reference library for the shore fishes retailed at Ambon harbor (Indonesia). Ecol. Evol. 2020, 10, 3356–3366. [Google Scholar] [CrossRef]
  33. Zhang, J.; Kapli, P.; Pavlidis, P.; Stamatakis, A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics 2013, 29, 2869–2876. [Google Scholar] [CrossRef] [Green Version]
  34. Fujisawa, T.; Barraclough, T.G. Delimiting species using single-locus data and the generalized mixed Yule coalescent approach: A revised method and evaluation on simulated data sets. Syst. Biol. 2013, 62, 707–724. [Google Scholar] [CrossRef] [Green Version]
  35. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  36. Bouckaert, R.; Heled, J.; Kühnert, D.; Vaughan, T.; Wu, C.H.; Xie, D.; Suchard, M.A.; Rambaut, A.; Drummond, A.J. BEAST 2: A software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 2014, 10, 1–6. [Google Scholar] [CrossRef] [Green Version]
  37. Bermingham, E.; McCafferty, S.; Martin, A.P. Fish biogeography and molecular clocks: Perspectives from the Panamanian isthmus. In Molecular Systematics of Fishes; Kocher, T.D., Stepien, C.A., Eds.; CA Academic Press: San Diego, CA, USA, 1997; pp. 113–128. [Google Scholar]
  38. Ogilvie, H.A.; Bouckaert, R.R.; Drummond, A.J. StarBEAST2 brings faster species tree inference and accurate estimates of substitution rates. Mol. Biol. Evol. 2017, 34, 2101–2114. [Google Scholar] [CrossRef]
  39. Ho, S.Y.W.; Larson, G. Molecular clocks: When times are a-changin’. TRENDS Genet. 2006, 22, 79–83. [Google Scholar] [CrossRef]
  40. Paradis, E. pegas: An {R} package for population genetics with an integrated-modular approach. Bioinformatics 2010, 26, 419–420. [Google Scholar] [CrossRef] [Green Version]
  41. Nei, M.; Tajima, F. DNA polymorphism detectable by restriction endonucleases. Genetics 1981, 97, 145–163. [Google Scholar] [CrossRef]
  42. Nei, M. Molecular Evolutionary Genetics; Columbia University Press: Chichester, NY, USA, 1987; ISBN 0231063210. [Google Scholar]
  43. Watterson, G.A. On the number of segregating sites in genetical models without recombination. Theor. Popul. Biol. 1975, 7, 256–276. [Google Scholar] [CrossRef]
  44. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar] [CrossRef]
  45. Kimura, M. A Simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide-sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef]
  46. Paradis, E.; Schliep, K. ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 2019, 35, 526–528. [Google Scholar] [CrossRef]
  47. Brown, S.D.J.; Collins, R.A.; Boyer, S.; Lefort, C.; Malumbres-Olarte, J.; Vink, C.J.; Cruickshank, R.H. Spider: An R package for the analysis of species identity and evolution, with particular reference to DNA barcoding. Mol. Ecol. Resour. 2012, 12, 562–565. [Google Scholar] [CrossRef]
  48. Ford, A.G.P.; Bullen, T.R.; Pang, L.; Genner, M.J.; Bills, R.; Flouri, T.; Ngatunga, B.P.; Rüber, L.; Schliewen, U.K.; Seehausen, O. Molecular phylogeny of Oreochromis (Cichlidae: Oreochromini) reveals mito-nuclear discordance and multiple colonisation of adverse aquatic environments. Mol. Phylogenet. Evol. 2019, 136, 215–226. [Google Scholar] [CrossRef] [Green Version]
  49. Nyingi, D.W.; Agnèse, J. Recent introgressive hybridization revealed by exclusive mtDNA transfer from Oreochromis leucostictus (Trewavas, 1933) to Oreochromis niloticus (Linnaeus, 1758) in Lake Baringo, Kenya. J. Fish Biol. 2007, 70, 148–154. [Google Scholar] [CrossRef]
  50. D’Amato, M.E.; Esterhuyse, M.M.; Van Der Waal, B.C.W.; Brink, D.; Volckaert, F.A.M. Hybridization and phylogeography of the Mozambique tilapia Oreochromis mossambicus in southern Africa evidenced by mitochondrial and microsatellite DNA genotyping. Conserv. Genet. 2007, 8, 475–488. [Google Scholar] [CrossRef]
  51. Gregg, R.E.; Howard, J.H.; Snhonhiwa, F. Introgressive hybridization of tilapias in Zimbabwe. J. Fish Biol. 1998, 52, 1–10. [Google Scholar] [CrossRef]
  52. Mojekwu, T.O.; Cunningham, M.J.; Bills, R.I.; Pretorius, P.C.; Hoareau, T.B. Utility of DNA barcoding in native Oreochromis species. J. Fish Biol. 2021, 98, 498–506. [Google Scholar] [CrossRef]
  53. Firmat, C.; Alibert, P.; Losseau, M.; Baroiller, J.-F.; Schliewen, U.K. Successive invasion-mediated interspecific hybridizations and population structure in the endangered cichlid Oreochromis mossambicus. PLoS ONE 2013, 8, e63880. [Google Scholar] [CrossRef]
  54. Seehausen, O. Hybridization and adaptive radiation. Trends Ecol. Evol. 2004, 19, 198–207. [Google Scholar] [CrossRef]
  55. Selz, O.M.; Seehausen, O. Interspecific hybridization can generate functional novelty in cichlid fish. Proc. R. Soc. B 2019, 286, 20191621. [Google Scholar] [CrossRef] [Green Version]
  56. Englebrecht, C.; Freyhof, J.; Nolte, A.; Rassman, K.; Schliewen, U.; Tautz, D. Phylogeography of the bullhead Cottus gobio (Pisces: Teleostei: Cottidae) suggests a pre-Pleistocene origin of the major central European populations. Mol. Ecol. 2000, 9, 709–722. [Google Scholar] [CrossRef]
  57. Fricke, R.; Mahafina, J.; Behivoke, F.; Jaonalison, H.; Léopold, M.; Ponton, D. Annotated checklist of the fishes of Madagascar, southwestern Indian Ocean, with 158 new records. FishTaxa 2018, 3, 1–432. [Google Scholar]
  58. Oswald, M.R.; Ravakarivelo, M.; Mikolasek, O.; Rasamoelina, H.; de Verdal, H.; Bentz, B.; Pepey, E.; Al, E. Combining a comprehensive approach to fish-farming systems with assessment of their genetics—From planning to realization. In Actes Projet FSP PARRUR, Recherche Interdisciplinaire pour le Développement Durable et la Biodiversité des Espaces Ruraux Malgaches. Application à Différentes Thématiques de Territoire; Parrur, Ed.; Cirad: Antananarivo, MG, USA, 2016; pp. 219–267. [Google Scholar]
  59. Borrell, Y.J.; Pineda, H.; McCarthy, I.; Vazquez, E.; Sanchez, J.A.; Lizana, G.B. Correlations between fitness and heterozygosity at allozyme and microsatellite loci in the Atlantic salmon, Salmo salar L. Heredity 2004, 92, 585–593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Lieutenant-Gosselin, M.; Bernatchez, L. Local heterozygosity-fitness correlations with global positive effects on fitness in threespine stickleback. Evolution 2006, 60, 1658–1668. [Google Scholar] [CrossRef] [PubMed]
  61. Borrell, Y.J.; Carleos, C.E.; Sánchez, J.A.; Vázquez, E.; Gallego, V.; Asturiano, J.F.; Blanco, G. Heterozygosity–fitness correlations in the gilthead sea bream Sparus aurata using microsatellite loci from unknown and gene-rich genomic locations. J. Fish Biol. 2011, 79, 1111–1129. [Google Scholar] [CrossRef]
Diversity 13 00281 g001 550
Figure 1. Sampling sites in Madagascar for the 262 DNA barcode records of Oreochromis spp. analyzed in this study. Site 1, O. niloticus GIFT strain at Matera (−18.0433, 49.3682); Site 2, O. niloticus GIFT strain at Tamatave (−18.449, 49.3987); Site 3, O. niloticus local strain at Ampefy (−19.0408, 46.7363); Site 4, O. niloticus JICA strain at Mahajunga (−15.6932, 46.3193); Site 5, O. niloticus local strain at Fenerive (−17.3666, 49.3961); Site 6, O. niloticus local strain at Brickaville (−18.8172, 49.0674); Site 7, O. mossambicus local strain at Fenerive (−17.3632, 49.3992); Site 8, O. macrochir local strain at Milasoa (−18.2224, 47.1461).
Figure 1. Sampling sites in Madagascar for the 262 DNA barcode records of Oreochromis spp. analyzed in this study. Site 1, O. niloticus GIFT strain at Matera (−18.0433, 49.3682); Site 2, O. niloticus GIFT strain at Tamatave (−18.449, 49.3987); Site 3, O. niloticus local strain at Ampefy (−19.0408, 46.7363); Site 4, O. niloticus JICA strain at Mahajunga (−15.6932, 46.3193); Site 5, O. niloticus local strain at Fenerive (−17.3666, 49.3961); Site 6, O. niloticus local strain at Brickaville (−18.8172, 49.0674); Site 7, O. mossambicus local strain at Fenerive (−17.3632, 49.3992); Site 8, O. macrochir local strain at Milasoa (−18.2224, 47.1461).
Diversity 13 00281 g001
Diversity 13 00281 g002 550
Figure 2. Mitochondrial gene tree for the 262 DNA barcodes of Oreochromis spp. inferred with SpeciesTreeUCLN, including 95% HPD interval for node age estimates, genetic lineage delimitation results for five methods (mGMYC discarded) and their 50% consensus, BOLD Barcode Index Numbers (BIN) for each MOTU, and distribution of farmed strain individuals.
Figure 2. Mitochondrial gene tree for the 262 DNA barcodes of Oreochromis spp. inferred with SpeciesTreeUCLN, including 95% HPD interval for node age estimates, genetic lineage delimitation results for five methods (mGMYC discarded) and their 50% consensus, BOLD Barcode Index Numbers (BIN) for each MOTU, and distribution of farmed strain individuals.
Diversity 13 00281 g002
Table
Table 1. Summary of genetic distances and species composition of each MOTUs including their BOLD Barcode Index Number (BIN), number of individuals in present study, number of publicly available records in BOLD (not including newly generated records here), maximum intraspecific and minimum interspecific K2P genetic distances, and relative proportions of species within MOTUs according to the present study and in publicly available BOLD records.
Table 1. Summary of genetic distances and species composition of each MOTUs including their BOLD Barcode Index Number (BIN), number of individuals in present study, number of publicly available records in BOLD (not including newly generated records here), maximum intraspecific and minimum interspecific K2P genetic distances, and relative proportions of species within MOTUs according to the present study and in publicly available BOLD records.
NK2P DistanceO. niloticus (Freq.)O. mossambicus (Freq.)O. macrochir (Freq.)O. aureus (Freq.)O. urolepis (Freq.)
BINPresent StudyBOLDWithin (Max)Between (Min)Present StudyBOLDPresent StudyBOLDPresent StudyBOLDPresent StudyBOLDPresent StudyBOLD
BOLD:AAA653743040.0040.0650.750.330.250.0001---0.3--
BOLD:AAA8511222340.0040.0350.99950.140.00050.52------
BOLD:AAA8513371800.016----10.17----
BOLD:AAC99041724170.0070.0500.99980.730.00020.0001------
BOLD:ACR7163162400.0360.56-0.440.01-----0.55
BOLD:ADI0792112800.0160.9-0.10.29------
Table
Table 2. Summary statistics of genetic diversity per MOTU. N, number of individuals; h, number of haplotypes; Hd, haplotypic diversity; π, nucleotide diversity; θw, theta; Tajima’s D test including D value and significance of the test (* significant at 0.025 threshold, none of the tests were significant).
Table 2. Summary statistics of genetic diversity per MOTU. N, number of individuals; h, number of haplotypes; Hd, haplotypic diversity; π, nucleotide diversity; θw, theta; Tajima’s D test including D value and significance of the test (* significant at 0.025 threshold, none of the tests were significant).
SpeciesStrainLocalitySiteNhHdπθwDp-Value
O. niloticusGIFTMatera11630.5750.03118.0820.4340.664
O. niloticusGIFTTamatave23180.5390.01815.519−1.2620.207
O. niloticusLocalAmpefy35870.6220.0168.2090.5540.579
O. niloticusJICAMahajunga44640.243<0.0010.455−1.4810.127
O. niloticusLocalFenerive53030.3930.0088.582−1.7190.086
O. niloticusLocalBrickaville63160.7850.03714.2682.0860.037
O. mossambicusLocalFenerive71360.7180.03723.520.0670.947
O. macrochirLocalMilasoa83720.054<0.0010.240−1.1310.258
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hubert, N.; Pepey, E.; Mortillaro, J.-M.; Steinke, D.; Andria-Mananjara, D.E.; de Verdal, H. Mitochondrial Genetic Diversity among Farmed Stocks of Oreochromis spp. (Perciformes, Cichlidae) in Madagascar. Diversity 2021, 13, 281. https://doi.org/10.3390/d13070281

AMA Style

Hubert N, Pepey E, Mortillaro J-M, Steinke D, Andria-Mananjara DE, de Verdal H. Mitochondrial Genetic Diversity among Farmed Stocks of Oreochromis spp. (Perciformes, Cichlidae) in Madagascar. Diversity. 2021; 13(7):281. https://doi.org/10.3390/d13070281

Chicago/Turabian Style

Hubert, Nicolas, Elodie Pepey, Jean-Michel Mortillaro, Dirk Steinke, Diana Edithe Andria-Mananjara, and Hugues de Verdal. 2021. "Mitochondrial Genetic Diversity among Farmed Stocks of Oreochromis spp. (Perciformes, Cichlidae) in Madagascar" Diversity 13, no. 7: 281. https://doi.org/10.3390/d13070281

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop