Outline
Comptes Rendus

Molecular biology and genetics/Biologie et génétique moléculaires
Genetic diversity and phylogenetic analysis of two Tunisian bivalves (Mactridae) Mactra corallina (Linnaeus, 1758) and Eastonia rugosa (Helbling, 1799) based on COI gene sequences
Comptes Rendus. Biologies, Volume 339 (2016) no. 3-4, pp. 115-122.

Abstracts

A partial sequence of mitochondrial cytochrome c oxidase subunit I (COI) was used as a genetic marker for a genetic diversity and phylogenetic analysis (DNA barcoding) of two Mactridae species, Mactra corallina and Eastonia rugosa, collected from the Tunisian coast. These Mactridae species could be distinguished by DNA barcoding techniques and they will be considered as monophyletic clades with the Neighbor-Joining (NJ) tree. The genetic structure detected that E. rugosa presents three haplotypes with a high frequency of HER1 (0.89). However, M. corralina shared 14 haplotypes. The haplotypic diversity (H) was equal to 0.205 and 0.954, respectively, for E. rugosa and M. corallina. While the nucleotide diversity (π) was higher for M. corallina (π = 0.0818), the mismatch distribution showed a unimodal curve for E. rugosa (a recent sudden demographic expansion) and a multimodal distribution for M. corallina (size stability).

Une séquence partielle du cytochrome c oxydase mitochondriale (COI) a été utilisée comme un marqueur génétique pour une étude de la diversité génétique et une analyse phylogénétique (DNA barcoding) de deux espèces de Mactridae, Mactra corallina et Eastonia rugosa, récoltées le long des côtes tunisiennes. Ces deux espèces de Mactridae pouvaient être identifiées par les techniques de DNA barcoding, et elles seront représentées dans des clades monophylétiques au niveau de l’arbre Neighbor-Joining (NJ). La structure génétique a révélé qu’E. rugosa présente trois haplotypes, avec une fréquence élevée d’HER1 (0,89). Cependant, M. corralina montre 14 haplotypes. La diversité haplotypique (H) est égale à 0,205 et 0,954, respectivement, pour E. rugosa et M. corallina. Alors que la diversité nucléotidique (π) est plus élevée chez M. corallina (π = 0,0818), la répartition des non-concordances montre une courbe unimodale pour E. rugosa (une expansion démographique soudaine et récente) et une distribution multimodale pour M. corallina (stabilité de taille).

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crvi.2016.02.001
Keywords: Mactra corallina, Eastonia rugosa, COI gene, Genetic diversity, Phylogenetic analysis
Mots-clés : Mactra corallina, Eastonia rugosa, Gène COI, Diversité génétique, Analyse phylogénétique

Imene Chetoui 1; Françoise Denis 2, 3; Mohamed Boussaid 4; Khoula Telahigue 1; M’Hamed El Cafsi 1

1 Unité de physiologie et d’écophysiologie des organismes aquatiques, département de sciences biologiques, faculté de sciences de Tunis, campus universitaire, 2092, Manar II, Tunis, Tunisie
2 Station de biologie marine, Muséum national d’histoire naturelle, place de la Croix, 29900 Concarneau, France
3 Laboratoire Mer, Molécules, Santé, Université du Maine, EA 2160, avenue O.-Messiaen, 72085 Le Mans cedex 9, France
4 Laboratoire de biotechnologie végétale, Institut national des sciences appliquées de technologie, centre urbain Nord, BP 676, 1080 Tunis cedex, Tunisie
@article{CRBIOL_2016__339_3-4_115_0,
     author = {Imene Chetoui and Fran\c{c}oise Denis and Mohamed Boussaid and Khoula Telahigue and M{\textquoteright}Hamed El Cafsi},
     title = {Genetic diversity and phylogenetic analysis of two {Tunisian} bivalves {(Mactridae)} {\protect\emph{Mactra} corallina {(}Linnaeus,} 1758) and {\protect\emph{Eastonia} rugosa} {(Helbling,} 1799) based on {COI} gene sequences},
     journal = {Comptes Rendus. Biologies},
     pages = {115--122},
     publisher = {Elsevier},
     volume = {339},
     number = {3-4},
     year = {2016},
     doi = {10.1016/j.crvi.2016.02.001},
     language = {en},
}
TY  - JOUR
AU  - Imene Chetoui
AU  - Françoise Denis
AU  - Mohamed Boussaid
AU  - Khoula Telahigue
AU  - M’Hamed El Cafsi
TI  - Genetic diversity and phylogenetic analysis of two Tunisian bivalves (Mactridae) Mactra corallina (Linnaeus, 1758) and Eastonia rugosa (Helbling, 1799) based on COI gene sequences
JO  - Comptes Rendus. Biologies
PY  - 2016
SP  - 115
EP  - 122
VL  - 339
IS  - 3-4
PB  - Elsevier
DO  - 10.1016/j.crvi.2016.02.001
LA  - en
ID  - CRBIOL_2016__339_3-4_115_0
ER  - 
%0 Journal Article
%A Imene Chetoui
%A Françoise Denis
%A Mohamed Boussaid
%A Khoula Telahigue
%A M’Hamed El Cafsi
%T Genetic diversity and phylogenetic analysis of two Tunisian bivalves (Mactridae) Mactra corallina (Linnaeus, 1758) and Eastonia rugosa (Helbling, 1799) based on COI gene sequences
%J Comptes Rendus. Biologies
%D 2016
%P 115-122
%V 339
%N 3-4
%I Elsevier
%R 10.1016/j.crvi.2016.02.001
%G en
%F CRBIOL_2016__339_3-4_115_0
Imene Chetoui; Françoise Denis; Mohamed Boussaid; Khoula Telahigue; M’Hamed El Cafsi. Genetic diversity and phylogenetic analysis of two Tunisian bivalves (Mactridae) Mactra corallina (Linnaeus, 1758) and Eastonia rugosa (Helbling, 1799) based on COI gene sequences. Comptes Rendus. Biologies, Volume 339 (2016) no. 3-4, pp. 115-122. doi : 10.1016/j.crvi.2016.02.001. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/j.crvi.2016.02.001/

Version originale du texte intégral

1 Introduction

Among the classes of mollusks, bivalve species play an important role in marine ecosystems. As filter feeders, they form a significant link between primary producers and act as calcium and carbon accumulators, used for shell construction. Mactra corallina (4.5 to 5.5 cm) and Eastonia rugosa (5.1 to 16.5 cm) are two species belonging to the Mactridae family. They are characterized by an equivalve and inequilateral shell; each valve has two cardinal and lateral teeth, internal ligament, and pallial sinus [1]. They are benthic and sympatric clams living in the infralittoral area of the beach, at depths from 1 m to 100 m. Their geographical distribution extends in the Mediterranean Sea and the Atlantic Ocean [2].

Some species of Mactridae (Macro chinensis and Macra veneriformis) are an important commercial seafood [3–5]; however, M. corallina and E. rugosa are not yet included in the Mediterranean market. The reproductive cycle, biochemical composition and ecology of M. veneriformis [6–8], M. chinensis [9], and Coelomactra antiquata [10] have been investigated. Genetic structure and differentiation were examined for different species of the Mactridae family [11], more specifically in C. antiquate from the Pacific Sea [12] and in two sibling species of the genus Mactra (M. corallina and M. lignaria) [13] from the Adriatic Sea as well as for three Tunisian populations of M. corallina using allozymic and morphological markers [14]. Moreover, scientific studies of E. rugosa species (Helbling, 1799) are rare. Therefore, no genetic marker has been tested on these species.

The identification of Mactridae species is difficult by morphological markers due to the existence of phenotypic plasticity and morphological convergence. Nevertheless, Habe (1977) [15] proposed a classification system for the Mactridae family, but no phylogenetic relationship is still available. Despite this classification, some Mactra species showed confusion. For this reason, the genetic markers might help solve effectively the confusion and distinguish some closely allied Mactridae species.

Since 1990, molecular markers have been used to detect polymorphisms among and within species [16]. They are useful because DNA was more thermostable than isoenzymes and because of the higher probability of variation being neutral. So, DNA was considered as a set of genetic markers involved in 80% of phylogenetic studies [17]. Currently, a fragment of the mitochondrial gene cytochrome c oxidase I (COI) was designed as the molecular barcode “DNA Barcoding” and it was suggested to use it as a universal marker for the identification of the animal's species [18]. Some authors used this reliable method to explain the interspecific relationships between marine bivalves and phylogeographic relationships among Atlantic, Pacific, and Mediterranean species [19–22], so it will be used to identify species and to analyze the structure and phylogeny of populations.

A work carried out on M. corallina [13] resolved the confused and contradictory systematics of the putative species of M. corallina collected from the North Adriatic Sea by analyzing molecular markers (12S, 16S, 18S and COI genes, partial sequences and morphological features). Nevertheless, no studies had been reported about the genetic structure along the South Mediterranean coast.

The aim of the present study is to analyze the genetic variation of two Tunisian Mactridae species, M. corallina and E. rugosa, using the partial sequences of mitochondrial cytochrome c oxidase (COI) to evaluate the potential of DNA barcoding to identify the species of Mactridae and to investigate the phylogenetic relationships with some Mactridae species.

2 Materials and methods

2.1 Specimen sampling

The studied species M. corallina and E. rugosa were represented by 18 and 19 individuals, respectively. The samples were collected in sandy beaches of Kalaât El Andalous (KA) in the month of Medjerda wadi in the Gulf of Tunis, in northern Tunisia, at latitude 33°41′20′′N and longitude 10°22′40′′E (Fig. 1).

Fig. 1

Map of Tunisia: location of the sampling site.

From each sample, a fragment of adductor muscles was taken and preserved at room temperature in absolute ethanol before DNA extraction.

2.2 DNA extraction, PCR amplification and sequencing

DNA was extracted using a 2% CTAB solution and standard proteinase K [23]. A 591-pb fragment of the cytochrome c oxidase subunit I (COI) of mtDNA was amplified from genomic DNA (Fig. 2). The amplification primers: FuCOI-a (5’GGT-CAA-CAA-ATC-ATA-AAG-ATA-TTG-G3’) and RuCOI-a (5’TAA-ACT-TCA-GGG-TGA-CCA-AAA-AAT-CA3’) were designed according to the method of Folmer et al. [24].

Fig. 2

Electrophoresis gel of the COI marker.

PCR amplification reactions were obtained for a total volume of 25 μL containing pure Taq Amersham Ready-To-Go-PCR Kit, 23 μL H2O, 1 μL of each primer (10 pmol) and 1 μL of ADN extract (25 μg/mL). Amplifications were carried out on a Biometra TGRADIENT Thermal Cycler. The program began by a pre-cycling denaturation cycle at 94 °C for 3 min; 40 cycles of denaturing at 94 °C for 1 min, annealing at 40 °C for 30 s, extension at 72 °C for 2 min, and a post-cycling extension at 72 °C for 10 min. The size and quality of PCR products were visualized on a 1% agarose gel with ethidium bromide on a UV transilluminator. Purification of DNA target fragments was carried out using a QIAquick gel extraction Kit (Qiagen) following the manufacturer's instructions. They were sequenced using the Big Dye TM terminator cycle sequencing chemistry following the manufacturer's protocol. The sequences were recorded with an ABI Prism TM 310 automated sequencer.

2.3 Data analysis

All sequences were aligned using the multiple-alignment program CLUSTAL W [25] as implemented in the software BioEdit v. 7.0.5.3. Sequence alignment was adjusted by visual inspection. The nucleotide sequences were deposited in NCBI with accession numbers from KM673272 to KM673288. From the NCBI (Genbank), we make a blast.

The specific genetic diversity was performed using the programs DNAsp 5.10.01 [26] and ARLEQUIN 3.0 [27]. The number of haplotypes, polymorphic sites, transitions, transversions, haplotype diversity (H), nucleotide diversity (π) and mean number of pairwise differences (k) were evaluated. The interpopulation divergence was determined by the FST index. The significance of FST was examined with a bootstrap analysis with 1000 replicates; P-values of pairwise tests were adjusted with the sequential Bonferroni correction. The genetic relationships among haplotypes for the COI gene sequence data were reconstructed using the Neighbor-Joining tree model using MEGA 4 software [28]. Using program DnaSP 5.10 [26], the null hypothesis of neutral evolution of the COI gene was tested using Fu and Li's D* and F* statistics, [29] and Tajima's D test [30]; the significant D values could be due to factors such as selection population expansion and bottlenecks [31].

The demographic history was inferred by analyzing the mismatch distribution, and the validity of the expansion model was evaluated using the sum of square deviations (SSD) between the observed mismatch and the expected one [32], Ramos-Onsis and Rozas R2 statistics, [33] and raggedness statistics rg [34] were also computed as implemented in DnaSP 5.10.

3 Results

3.1 Phylogenetic analysis

The Neighbor-Joining (NJ) tree (Fig. 3) showed two major clusters corresponding to the two Tunisian Mactridae species (Fig. 4). M. corallina population showed a high heterogeneity, which is characterized by the presence of 3 branches to more than 90% bootstraps. So, Fig. 2 shows the phylogenetic tree for the Mactridae family from COI gene analysis. According to this figure, we assume that all species were monophyletic, with support from 50 to 100%, and the genus Mactra was closer to the Coelomactra species, with support from 35% to 43%. However, E. rugosa was in a different clade, far from the genus Mactra (56%) and more closer to Spisula solidissima and Rangia curneata (80%). This bootstrap value (80%) is supported by the monophyly of the Eastonia and Spisula genus, which are two sister groups.

Fig. 3

Neighbor-joining tree (model Kimura-2 parameter) presenting the phylogenetic relationships among Mactra corallina, Eastonia rugosa and other species of Mactridae family based on the COI gene. Bootstrap (1000 replicates) HMA: Haplotype of M. corallina, and HER: haplotype of E. rugosa.

Fig. 4

Photo of the two Mactridae species.

3.2 Genetic diversity

For the two species, a COI fragment of length 591 pb was analyzed. E. rugosa presents three haplotypes with a high frequency of HER1 estimated at 0.894. M. corralina shared 14 haplotypes with frequencies ranging from 0.055 to 0.222 and with the highest frequency for HMA3 (H = 0.222). The nucleotide compositions were characterized by a lower rate of C (13.87% for E. rugosa and 16.29% for M. corallina) and a higher rate of T (37.50% for E. rugosa and 39.15% for M. corallina). So M. corallina was genetically variable as it showed a high number of polymorphic sites (126 substitutions), while E. rugosa was weakly variable, with only three polymorphic sites. Statistical analyzes showed that haplotypic diversity H was equal to 0.205 and 0.954, respectively, for E. rugosa and M. corallina. The nucleotide diversity (π) and the mean number of pairwise differences (k) were higher for M. corallina (π = 0.0818 and k = 48.356) (Table 1).

Table 1

Molecular diversity indices for COI sequences of Mactra corallina and Eastonia rugosa (Mactra rugosa).

Species Mactra corallina Eastonia rugosa
N 18 19
Haplotype number 14 3
Polymorphic sites 126 3
H 0.954 ± 0.0019** 0.205 ± 0.014*
π 0.0818 ± 0.0416* 0.00053 ± 0.00063**
k 48.356 ± 21.987 0.316 ± 0.336
Number of sites with transitions 97 2
Number of sites with transversions 41 1
Number of sites with subtitutions 138 3

The mismatch frequency spectra for the species studied are shown in Fig. 5. The observed distribution was unimodal for the E. rugosa population, indicating a demographic expansion in the past. This result was supported by the significant negative values in Tajima's D test (D = –1.71, P < 0.05), Fu's Fu test (Fu's Fu = –1.084, P < 0.05) Fu and Li's D* test (Fu D* = –2.34, P < 0.001) and Fu & Li's F* test (Fu's F* = –2.24, P < 0.001). So, the values of the sum of squared differences (SSD = 0.008, P > 0.1), the raggedness index (rg = 0.49, P > 0.1) and the Ramos-Onsins and Rozas test (R2 = 0.162, P > 0.1) were non-significant. The hypothesis of selective neutrality is rejected. So it confirms a recent sudden expansion population scenario. For M. corallina, the observed multimodal mismatch distribution underlined a constant population size. This population showed non-significant values of neutrality tests: Tajima's D test (D = 1.35, P > 0.1), Fu's Fu test (Fu's Fu = 2.25, P > 0.1), and Fu F* test (Fu F* = 0.26, P > 0.1); they suggested a contraction of population size due to the balanced selection. Moreover, the value of the sum of squared differences (SSD) was high (SSD = 0.075, P < 0.1), which rejects the contraction hypothesis and confirms the stability size of the M. corallina population (Table 2).

Fig. 5

Pairwise mismatch distributions; observed (dashed line) and expected (solid line) curves for two species of Mactridae. X-axis: number of pairwise differences; Y-axis: frequency of mismatches; Exp: expected and Obs: observed.

Table 2

Neutrality test, mismatch distribution fitting to a sudden expansion model and estimates of the expansion statistic R2 for Mactra corallina and Eastonia rugosa (Mactra rugosa). P-values were obtained by coalescent simulations with 100 replicates.

Species N Test D (Tajima) Test D* (Fu and Li) Test F* (Fu and Li) Fu's Fs R 2 rg Mismatch SSD
Mactra corallina 18 1.353NS –0.236 0.269 0.250NS 0.184 0.04 48.366 0.075*
Eastonia rugosa 19 –1.719* –2.345** –2.498** –1.084* 0.162 0.49 0.316 0.008NS

4 Discussion

The development of molecular techniques and tools could to reduce the problems of confusion of similar morphological sibling species and to identify them.. A new approach based on COI sequences could overcome these difficulties. It has been efficient in many animal species (e.g. birds [35], fishes [36], crustaceans, [37] and mollusks [38]. Also, the COI marker is considered a good genetic marker because of its power for identifying specimens [39], for evaluating biodiversity [40], for discriminating the species genetically close [41], and for defining the cryptic species [42]. But this marker presents two criteria for species delimitations: (1) the unequal ratio of intraspecific variations and interspecific divergence, and (2) the reciprocal monophyly of species [35,43]. So, in our study, we were interested first to analyze the phylogenetic relationship between haplotypes for two clam species (M. corallina and E. rugosa). Some different sequences of other species of the Mactridae family from GenBank were used for comparative purposes (Table 3). According to the N-J tree, we observed as well that the sequences of M. Corallina from our study and sequences of M. corallina coralllina and M. Corallina lignaria from gene bank formed a single clade with a high bootstrap support (63% to 99%). We assumed that all species of the Mactridae family were monophyletic (having a common ancestor) and M. corallina species was considered a paraphyletic group. On the other hand, we considered that the genus Coelomactra was closer to Mactra (43%) than to other genus of Rangia (33%), Spisula, and Eastonia (21%). So our results showed D = 0.24 between M. corallina and Mactra species and D = 0.27 between M. corallina and Coelomactra species). This value of D fell between the limits established for these species by [11]. Only the genus (Mactra) was reciprocally monophyletic with 99% bootstrap support. Several authors have confirmed these hypotheses after studying the phylogeny of Mactridae species based on genetic markers (16S and COI) [11] and on two other molecular markers, 12S and18S [14].

Table 3

Species and GenBank accession numbers.

Species GenBank accession numbers
Mactra veneriformis JN674621
Mactra corallina corallina GQ166585 and FJ830444
Mactra corallina lignaria GQ166586 and FJ830437
Mactra alta JN674615
Mactra chinensis JN674630
Coelomactra cumingui JN674610
Lutraria australis JN674600
Lutraria arcuata JN674601
Spisula solodissima AY707811
Rangia cuneata U47652
Hexaplex trunculus ( outgroups) EU391577

E. rugosa was a paraphyletic species (with a separate group); it was closer to the genus Spisula and Rangia and showed low genetic distances (D) equal to 0.33 and 0.32, respectively.

The subdivision of the M. corallina population (MA) into two distinct clades (NJ) shows the high genetic divergence within this population. This large genetic divergence can be caused by the freshwater of Medjerda wadi. So, this genetic differentiation has been demonstrated by allozymic markers in M. Corallina in the same locality, whose difference in salinity could probably explain the presence of genetic differences (allelic frequencies difference and absence of some alleles) between populations. As well Kong et al. (2007) [44] suggest that the freshwater outflow of the Yangtze River and upwellings of the Zhejiang province are probably the causes of the high genetic divergence in the C. antiquata populations.

In a second part, we analyzed the genetic variation of M. corallina and E. rugosa using the cytochrome c oxidase subunit I (COI) region of mtDNA. The results demonstrate that DNA sequence polymorphism in the COI was effective in the characterization of interspecific structures.

The value of haplotype diversity (H) of mtDNA of M. corallina (0.954) appears to be relatively higher than those reported in other bivalve species [21,22,45,46] and show a great similarity with those recorded in Spissula solidissima [47], while this value for E. rugosa was lower than for the species previously cited. There are several possible explanations for such reduced genetic diversity of E. rugosa: it could result from population life history traits (after demographic crashes or selective events) or be related to the current population's size. But the high haplotype and nucleotide diversity of M. corallina are related to the stability of the population. So, haplotype and nucleotide diversity depend essentially on the pelagic larval phase and on their dispersal ability. Both biological phenomena have an impact on population size [48].

Moreover, as COI is considered more polymorphic than other mitochondrial genes [49], it was frequently used to study the population genetic structure of marine species [17,50] and could be very useful for phylogenetic reconstructions [22] for close species. Therefore, this marker showed a variable number of polymorphic sites for the two species of Mactridae. It demonstrated a range of nucleotide substitution rates that can be used for different evolutionary analyses. So we can say that the M. corallina population was more heterogeneous. However, the E. rugosa population seems to be homogenous, with a low number of polymorphic sites equal to 3.

For the E. rugosa population, Tajima's D neutrality tests and Fu's Fu test revealed significant negative values. A unimodal model mismatch distribution and high haplotype diversity (0.205) were observed. We can conclude to a sudden demographic expansion model of Rogers and Harpending [34]. These negatives values of neutrality tests will be more effective to reject the null hypothesis of constant size and sustain a recent population expansion scenario due to recent selective sweeps. The low and non-significant values of the sum of squared differences (SSD) confirm this hypothesis and reject the bottleneck event. However, we noted that the expansion model of the E. rugosa population was represented by Tajima's D and Fu's Fu tests, but not by the R2 test, due to small sample size [51–53]. According to previous studies, the same signal in favor of a sudden and recent population expansion was found in bivalves, such as Pinna nobilis from the Mediterranean Sea [22], Cerastoderma glaucum in shallow waters [45], S. solidissima from the Atlantic coast [47], Ruditapes philippinarum in the Pacific coast [20], which was related to the recent colonization of the present area and the historical episodes corresponding to the major warm interstadial period during the last glacial cycle [54].

For the M. corallina population, we tested two hypotheses in order to determine the demographic history of this population. Firstly, non-significant and positives values of neutrality tests (Tajima's D and Fu's Fu tests) have suggested a contraction of the population size due to the balanced selection (multiple haplotypes with low and high frequencies). Secondly, the high and significant value of the sum of squared differences (SSD) and nucleotide diversity (π) confirmed the stability size of the M. corallina population. So, we can deduce that this stability of M. corallina population was significantly related to the power of the sum of squared differences (SSD) test. The same example was observed in cephalopods (Octopus vulgaris) from the Mediterranean Sea [55]. They demonstrated the demographic stability of these populations, determined by non-significant values of Fu's [33] FS and Ramos-Onsins and Rozas’ R2 [34] tests, knowing that they observed multimodal mismatch distributions for all populations.

5 Conclusion

The mitochondrial COI has demonstrated its power firstly for identifying Mactridae species (M. corallina and E. rugosa) in this study, secondly for explaining the genetic diversity and demographic history for both studied species.

Nevertheless, it is necessary to improve the genetic study of these two species, first by referring to other populations (samples) collected all along the Mediterranean coast, and then by using a vast other genetic markers such as (microsatellites).

Disclosure of interest

The authors declare that they have no competing interest.

Acknowledgements

Special thanks are extended to different members of the “Station de biologie marine” of Concarneau.


References

[1] R.C. Moore; C. Teichert Treatise on invertebrate paleontology, Part N, Mollusca 6, Bivalvia, Geological Society of America and University of Kansas, Lawrence, KS, USA, 1969

[2] FAO. Fiches FAO d’identification des espèces pour les besoins de la pêche en Méditerranée et mer Noire, zone de pêche, 1987 (760 p)

[3] D.K. Ryou Studies on the population dynamics of surf clam, Mactra veneriformis Reeve (Bivalvia) in the coast of Kunsan, Korea I. Production, Korean. J. Malacol., Volume 13 (1997), pp. 193-201

[4] Q.Z. Qi; R.C. Gao; W.R. Qiu; X.Q. Huang The life history of Coelomactra antiquata, J. Fujian Norm. Univ., Volume 11 (1995), pp. 82-88

[5] J.F. Wu; H.H. Zhang; C.Y. Liang; S.W. Chen Resource and conservation strategy of Coelomactra antiquata in coast of Guangdong province, J. Ocean Univ. China, Volume 22 (2002), pp. 68-69

[6] K. Qiaozhen; Li. Qi Annual dynamics of glycogen, lipids, and proteins during the reproductive cycle of the surf clam Mactra veneriformis from the north coast of Shandong Peninsular, Invert. Reprod. Dev., Volume 57 (2013) no. 1, pp. 49-60

[7] H. Kazuyoshi; T. Norihisa; K. Kunio; N. Yasuo; M. Masatoshi Life history characteristics of the surf clam Mactra veneriformis (Bivalvia: Veneroida: Mactridae) on a sandy tidal flat in Tokyo Bay, Pac. Sci., Volume 66 (2012) no. 3, pp. 335-346

[8] T. Nakano; J. Nasuda; Y. Agata; T. Yurimoto; Y. Maeno; Y. Nakamura; F. Yamada; A. Tamaki Life history and population dynamics of the surf clam, Mactra veneriformis (Bivalvia: Mactridae), on an estuarine intertidal sandflat in western Kyushu, Molluscan Res., Volume 32 (2012) no. 3, pp. 159-176

[9] Li. Qi; Y. Lin; K. Qiaozhen; K. Lingfeng Gametogenic cycle and biochemical composition of the clam Mactra chinensis (Mollusca: Bivalvia): implications for aquaculture and wild stock management, Mar. Biol. Res., Volume 7 (2011), pp. 407-415

[10] D.J. Liu; K.E. Xie Reproductive biology of Coelomactra antiquata, J. Chin. Zool., Volume 38 (2003), pp. 10-15

[11] Ni. Lehai; Li. Qi; K. Lingfeng; H. Shuqing; Li. Lujing DNA barcoding and phylogeny in the family Mactridae (Bivalvia: Heterodonta): evidence for cryptic species, Biochem. Syst. Ecol., Volume 44 (2012), pp. 164-172

[12] K. Lingfeng; Li. Qi Genetic evidence for the existence of cryptic species in an endangered clam Coelomactra antiquata, Mar. Biol., Volume 156 (2009), pp. 1507-1515

[13] I. Guarniero; F. Plazzi; A. Bonfitto; A. Rinaldi; M. Trentini; M. Passamonti The bivalve mollusc Mactra corallina: genetic evidence of existing sibling species, J. Mar. Biol. Assoc. UK, Volume 90 (2010) no. 3, pp. 633-644

[14] I. Chetoui; M. El Cafsi; M. Boussaid Allozymic and morphological variation in three populations of surf clam Mactra corallina (Bivalvia: Mactridae) from Tunisian sandy beaches, Cah. Biol. Mar., Volume 53 (2012), pp. 409-417

[15] T. Habe; N. Lehai; L. Qi; K. Lingfeng; H. Shuqing; L. Lujing DNA barcoding and phylogeny in the family Mactridae (Bivalvia: Heterodonta): Evidence for cryptic species, Biochem Syst Ecol, Volume 44 (1977), pp. 164-172 (in Japanese)

[16] T. Sugita; A. Nishikawa; T. Shinoda Identification of Trichosporon asahii by PCR based on sequences of the internal transcribed spacer regions, J. Clin. Microbiol., Volume 36 (1998), pp. 2742-2744

[17] J.C. Avise Phylogeography. The history and formation of species, Harvard University Press, Cambridge, MA, USA, 2000 (464 p)

[18] P.D.N. Hebert; A. Cywinska; S.L. Ball; J.R. deWaard Biological identifications through DNA barcodes, Proc. R. Soc. Lond. B, Volume 270 (2003), pp. 313-321

[19] A. Wood; S. Apte; E.S. MacAvoy; J.P.A. Gardner A molecular phylogeny of the marine mussel genus Perna (Bivalvia: Mytilidae) based on nuclear (ITS1&2) and mitochondrial (COI) DNA sequences, Mol. Phylogenet. Evol., Volume 44 (2007), pp. 685-698

[20] M. Yangli; G. Tianxiang; Y. Takashi; X. Yongshuang Fishery molecular phylogeography of Ruditapes philippinarum in the Northwestern Pacific Ocean based on COI gene, J. Exp. Mar. Biol. Ecol., Volume 407 (2011), pp. 171-181

[21] A. Gharbi; N. Chatti; K. Said; A.V. Wormhoudt Genetic variation and population structure of the carpet shell clam Ruditapes decussatus along the Tunisian coast inferred from mtDNA and ITS1 sequence analysis, Biologia, Volume 65 (2010) no. 4, pp. 688-696

[22] L. Rabaoui; R. Mejri; S. Tlig-Zouari; L. Bahri; O.M. Ben Hassine; C.T. Tsigenopoulos Genetic variation among populations of the endangered fan mussel Pinna nobilis (Mollusca: Bivalvia) along the Tunisian coastline, Hydrobiologia, Volume 678 (2011), pp. 99-111

[23] B. Winnepennickx; T. Backeljan; R. Dewachter Extraction of high molecular weight DNA from mollusks, Trends Genet., Volume 9 (1993) no. 12, p. 407

[24] O. Folmer; M. Black; R. Lutz; R. Vrijenhoek DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from metazoan invertebrates, Mol. Mar. Biol. Biotechnol., Volume 3 (1994), pp. 294-299

[25] J.D. Thompson; D.G. Higgins; T.J. Gibson; W. Cllustal Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice, Nucleic Acids Res., Volume 12 (1994), pp. 4673-4680

[26] J. Rozas; P. Librado; J.C. Sanchez-DelBarrio; X. Messeguer; R. Rozas DNA sequence polymorphism, version 5.10.01, 2010 (A software package for the analysis of nucleotide polymorphism from aligned DNA sequence data)

[27] L. Excoffier; A. Estoup; J.-M. Cornuet Bayesian analysis of an admixture model with mutations and arbitrarily linked markers, Genetics, Volume 169 (2005), pp. 1727-1738

[28] K. Tamura; M. Nei Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees, Mol. Biol. Evol., Volume 10 (1993), pp. 512-526

[29] Y.X. Fu; W.H. Li Statistical tests of neutrality of mutations, Genetics, Volume 133 (1993), pp. 693-709

[30] F. Tajima Statistical method for testing the neutral mutation hypothesis by DNA polymorphism, Genetics, Volume 123 (1989), pp. 585-595

[31] F. Tajima The effect of change in population size on DNA polymorphism, Genetics, Volume 123 (1989), pp. 597-601

[32] Y.X. Fu Statistical tests of neutrality of mutations against population growth, hitchhiking, and background selection, Genetics, Volume 147 (1997), pp. 915-925

[33] S.E. Ramos-Onsins; J. Rozas Statistical properties of new neutrality tests against population growth, Mol. Biol. Evol., Volume 19 (2002), pp. 2092-2100

[34] A.R. Rogers; H. Harpending Population growth makes waves in the distribution of pairwise genetic differences, Mol. Biol. Evol., Volume 9 (1992), pp. 552-569

[35] P.D.N. Hebert; M.Y. Stockle; T.S. Zemlak; C.M. Francis Identification of birds through DNA barcodes, PLoS Biol., Volume 2 (2004), p. e312

[36] R.D. Ward; T.S. Zemlak; B.H. Innes; P.R. Last; P.D.N. Hebert DNA barcoding Australia's fish species, Philos. Trans. R. Soc. Lond. B, Volume 360 (2005), pp. 1847-1857

[37] T. Lefebure; C.J. Douady; M. Gouy; J. Gibert Relationship between morphological taxonomy and molecular divergence within Crustacea: proposal of a molecular threshold to help species delimitation, Mol. Phylogenet. Evol., Volume 40 (2006), pp. 435-447

[38] Y. Feng; Q. Li; L. Kong; X. Zheng COI-based DNA barcoding of Arcoida species (Bivalvia: Pteriomorphia) along the coast of China, Mol. Ecol. Resour., Volume 11 (2011), pp. 435-441

[39] W.S. Lakra; M. Goswami; A. Gopalakrishnan; D.P. Singh; A. Singh; N.S. Nagpure Genetic relatedness among fish species of genus Channa using mitochondrial DNA genes, Biochem. Syst. Ecol., Volume 38 (2010), pp. 1212-1219

[40] M.A. Smith; B.L. Fisher; P.D.N. Hebert DNA barcoding for effective biodiversity assessment of a hyperdiverse arthropod group: the ants of Madagascar, Philos. Trans. R. Soc. B, Volume 360 (2005), pp. 1825-1834

[41] P.D.N. Hebert; S. Ratnasingham; J.R. deWaard Barcoding animal life: cytochrome c oxidase subunit I divergences among closely related species, Proc. R. Soc. Lond. B, Volume 270 (2003), p. S96-S99

[42] A. Bucklin; P.H. Wiebe; S.B. Smolenack; N.J. Copley; J.G. Beaudet; K.G. Bonner; J. Farber-Lorda; J.J. Pierson DNA barcodes for species identification of euphausiids (Euphausiacea, Crustacea), J. Plankton Res., Volume 29 (2007), pp. 483-493

[43] J.J. Wiens; T.A. Penkrot Delimiting species using DNA and morphological variation and discordant species limits in spiny lizards (Sceloporus), Syst. Biol., Volume 51 (2002), pp. 69-91

[44] L. Kong; Q. Li; Z. Qi Genetic and morphological differentiation in the clam Coelomactra antiquata (Bivalvia: Veneroida) along the coast of China, J. Exp. Mar. Biol. Ecol., Volume 343 (2007), pp. 110-117

[45] R. Nikula; R. Vainola Phylogeography of Cerastoderma glaucum (Bivalvia: Cardiidae) across Europe: a major break in the Eastern Mediterranean, Mar. Biol., Volume 143 (2003), pp. 339-350

[46] E.S. Cho; C.G. Joung; S.G. Sohn; C.D. Kim; S.J. Han Population genetic structure of the ark shell Scapharca broughtonii Schrenck from Korea, China, and Russia based on COI gene sequences, Mar. Biotechnol., Volume 9 (2007), pp. 203-216

[47] M.P. Hare; J.R. Weinberg Phylogeography of surfclams, Spisula solidissima, in the western North Atlantic based on mitochondrial and nuclear DNA sequences, Mar. Biol., Volume 146 (2005), pp. 707-716

[48] A.J. Bohonak Dispersal, gene flow and population structure, Q. Rev. Biol., Volume 74 (1999), pp. 21-45

[49] V. Katsares; A. Tsiora; S. Galinou-Mitsoudi; A. Imsiridou Genetic structure of the endangered species Pinna nobilis (Mollusca: Bivalvia) inferred from mtDNA sequences, Biologia, Volume 63 (2008), pp. 412-417

[50] C. Saccone; C. De Giorgi; C. Gissi; G. Pesole; A. Reyes Evolutionary genomics in Metazoa: the mitochondrial DNA as a model system, Genetics, Volume 238 (1999), pp. 195-209

[51] K.L. Simonsen; G.A. Churchill; C.F. Aquardo Properties of statistical tests of neutrality for DNA polymorphism data, Genetics, Volume 141 (1995), pp. 413-429

[52] F. Depaulis; S. Mousset; M. Veuille Haplotype tests using coalescent simulations conditional on the number of segregating sites, Mol. Biol. Evol., Volume 18 (2001), pp. 1136-1138

[53] S.R. Palumbi Population genetics, demographic connectivity, and the design of marine reserved, Ecol. Appl., Volume 13 (2003) no. 1, pp. 146-158

[54] K. Lambeck; T.M. Esat; E.K. Potter Links between climate and sea levels for the past three million years, Nature, Volume 419 (2002), pp. 199-206

[55] K. Fadhlaoui-Zid; L. Knittweis; D. Aurelle; C. Nafkha; O. Jarboui; F. Maltagliati Genetic structure of Octopus vulgaris (Cephalopoda, Octopodidae) in the central Mediterranean Sea inferred from the mitochondrial COIII gene, C. R. Biologies, Volume 335 (2012), pp. 625-636


Comments - Policy