Next Article in Journal
Variability in Clinical Phenotype in TARDBP Mutations: Amyotrophic Lateral Sclerosis Case Description and Literature Review
Previous Article in Journal
Association between NTRK2 Polymorphisms, Hippocampal Volumes and Treatment Resistance in Major Depressive Disorder
Previous Article in Special Issue
Structural Characteristics of Mitochondrial Genomes of Eight Treehoppers (Hemiptera: Membracidae: Centrotinae) and Their Phylogenetic Implications
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 11
10.3390/genes14112038
genes-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

Genetic Analysis Based on Mitochondrial nad2 Gene Reveals a Recent Population Expansion of the Invasive Mussel, Mytella strigata, in China

by
Shaojing Yan
1,
Peizhen Ma
1,2,*,
Chenxia Zuo
1,3,
Yi Zhu
1,4,
Xiaojie Ma
1,4 and
Zhen Zhang
1,4,*
1
Laboratory of Marine Organism Taxonomy & Phylogeny, Qingdao Key Laboratory of Marine Biodiversity and Conservation, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
2
State Key Laboratory of Mariculture Biobreeding and Sustainable Goods, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, China
3
College of Life Sciences, Qingdao University, Qingdao 266000, China
4
University of Chinese Academy of Sciences, Beijing 100049, China
*
Authors to whom correspondence should be addressed.
Genes 2023, 14(11), 2038; https://doi.org/10.3390/genes14112038
Submission received: 27 September 2023 / Revised: 31 October 2023 / Accepted: 1 November 2023 / Published: 3 November 2023
(This article belongs to the Special Issue Molecular Evolution, Mitochondrial Genomics and Mitochondrial Genome Expression in Animals)
Browse Figures
Versions Notes

Abstract

:
Mytella strigata is a highly adaptable invasive alien species that has been established in coastal China since 2014. Mitochondrial DNA (mtDNA) is an important tool for studying the evolution and population genetics of invasive species. In this study, the mitochondrial genome of M. strigata from China was sequenced by Illumina high-throughput sequencing and characterized with 13 protein-coding genes (PCGs). By assessing the selective pressure of 13 PCGs, the nad2 gene had the fastest evolutionary rate and was finally selected for population genetic analysis. A total of 285 nad2 sequences from seven M. strigata populations in China were analyzed and showed obviously T-rich and C-rich characteristics. According to population genetic diversity analysis, all the seven populations had haplotype (gene) diversity (Hd) ≥ 0.5 and nucleotide diversity (Pi) < 0.005. Haplotype networks showed a “star” distribution. Population historical dynamic analyses showed that Fu’s Fs and Tajima’s D values of all populations were negative except the Qukou (QK) and Beihai (BH) populations. The Zhangzhou (ZJ) and Xiamen (XM) populations were unimodal while the other populations were multimodal. These results suggested that the population of M. strigata in China may have passed the bottleneck period and is currently in a state of population expansion.

1. Introduction

Along with the global shipping industry’s vigorous development, invasive marine species are more likely to be introduced into new environments where they begin to establish a population [1]. The establishment and spread of invasive species in new areas could have wide-ranging ecological consequences [2]. At fine scales, invasive species may affect changes in the local gene pool [3]. At larger scales, they may affect the structure, composition, diversity, and biological interactions of the recipient communities [4]. Mollusks and crustaceans have the largest number among marine invasive species [5]. Mytella strigata (Hanley, 1843) (=Mytella charruana (d’Orbigny, 1846)), a marine invasive bivalve [6] native to the Atlantic coast of South America, has been found in many east and southeast Asian countries and regions [7,8]. Studies have shown that M. strigata is mainly transmitted over long distances through human activities such as ships’ ballast water, while ocean currents and coastal currents promote its spread in a small range [9]. Its strong adhesion and ability to survive in high densities enable its attachment to the walls and drainage systems of aquaculture ponds, hulls of boats, bottom sediment, and riverbanks of estuaries [10]. Moreover, it poses a serious competitive threat to native species such as Perna viridis (Linnaeus,1758) [11]. The ecological and social impacts of invasive species can be significant, but their management is often difficult [12]. At the beginning of an invasion, the genetic diversity of an alien species will generally decrease as they need to adapt to the local environment. This stage is called the bottleneck effect. Once the bottleneck is crossed, the genetic diversity of invasive species will increase and their populations will continuously expand [13]. This covert process can lead to an underestimation of invasive species and limit the effectiveness of management [14]. Therefore, studying population genetic structure is very important to understand the population status of invasive species and strengthen relevant control.
Mitochondrial DNA (mtDNA) is an important tool in researching evolution and population genetics of invasive species due to its relatively high mutation rate and apparent simplicity of maternal inheritance [15]. In most eukaryotes, mitochondrial inheritance usually follows strict maternal inheritance (SMI) [16]. However, a unique doubly uniparental inheritance pattern (DUI) has been reported in more than 100 species of bivalves [17]. Mitochondrial DUI is described as follows: the eggs and spermatozoa contain female-lineage mtDNA (F-type) and male-lineage mtDNA (M-type), respectively. Therefore, the zygote is heteroblastic. If the zygote develops into a female, the animal will be homogenous with the F-type in germ cells and somatic tissues [18]. Otherwise, if the zygote develops into a male, it will be heterogeneous with the M-type in germ cells and F-type in somatic tissues [17]. However, a new inheritance type, namely, atypical DUI, has emerged. The M-type is not only found in male germ cells but also in somatic tissues, such as Ruditapes philippinarum (Adams & Reeve, 1850) [19] and Mytilus galloprovincialis (Lamarck, 1819) [20]. In previous studies, M. strigata were also shown to possess atypical DUI, as M-COI was found in the adductor muscle [21,22]. Because the two haplotypes in DUI species are related to sex, the evolution rate of the haplotypes may be different, which may affect the accuracy of population genetic analysis [7,22]. At the same time, population genetic analysis also relies on the relatively high mutation of mitochondrial gene makers. In general, genes with less selective pressure are thought to have more mutations [23] and as a result, the ratios of the nonsynonymous substitution rate (Ka) and synonymous substitution rate (Ks) were commonly employed. Together, it is necessary for more suitable mitochondrial gene makers to be screened and employed for population genetic analysis.
Here, Illumina high-throughput sequencing was performed on the mitochondrial genome of M. strigata from China. Selective pressure was analyzed on a mitochondrial protein-coding gene (PCG) to choose the best mitochondrial marker for population genetic analysis. Then, the population genetic diversity and population history dynamics of seven populations of M. strigata in coastal China were carried out. The aim of this study was to (1) select more suitable mitochondrial PCGs for population genetic analysis, (2) evaluate the invasive status of M. strigata populations in coastal areas of China, and (3) provide molecular biological data for assessing the expansion and determine the need for control measures of M. strigata.

2. Materials and Methods

2.1. Mussels and Sampling

All experiments were performed in accordance with the relevant guidelines and regulations. In April 2023, a total of 252 individuals were collected from seven populations in China: Zhanjiang (ZJ), Beihai (BH), Shanwei (SW), Zhangzhou (ZZ), Qukou (QK), Gangbei (GB), and Xiamen (XM) (Figure 1, Table 1). All the samples were directly collected from the attachment substrates, including rocks, concrete, oyster shells, and net cages, in intertidal zones. All samples were identified as M. strigata according to the morphological descriptions (Figure S1) [10,11,24]. All samples were preserved in 95% alcohol and deposited in the Laboratory of Marine Organism Taxonomy and Phylogeny, Qingdao Key Laboratory of Marine Biodiversity and Conservation, Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China.

2.2. DNA Extraction, Sequencing, Assembly, and Annotation for Mitochondrial Genomes

One sample from ZJ was chosen for mitochondrial genome sequencing. Total DNA was extracted from the adductor muscle using a TIANamp Marine Animals DNA kit (DP324-03, Tiangen Biotech Co., Ltd., Beijing, China) according to the reagent instructions. The genomic library was constructed with the whole-genome shotgun approach and sequenced on the Illumina NovaSeq platform (Illumina, San Diego, CA, USA) at Shanghai Personal Biotechnology Co., Ltd. using the 2 × 150 bp paired-end sequencing mode. A5-miseq v20150522 software [25] and SPAdes v3.9.0 software [26] were used for the de novo assembly of high-quality data. The complete mitogenome sequences were uploaded to the MITOS2 web server for functional annotation [27]. The genetic code selection was set to “Invertebrate” and the other settings were adjusted according to the default parameters. The boundaries of PCGs were determined by an ORF finder (https://www.ncbi.nlm.nih.gov/orffinder, accessed on 13 October 2023), and then they were corrected manually and annotated by comparison with the previous mitochondrial genomes of M. strigata (Genbank accession nos.: MT800514 and MT991018).

2.3. Selective Pressure (Ka/Ks) Analysis

PhyloSuite v1.2.2 software was employed to perform the preparation of the PCG sequences from mitochondrial genomes sequenced in this study [28]. The remaining two mitochondrial genomes used for selective pressure analysis came from GenBank (GenBank accession nos.: MT800514 and MT991018). Multiple Alignment Using Fast Fourier Transformation (MAFFT v7.505) software was used to align PCGs independently [29]. The codon mode for PCGs was selected using the invertebrate mitochondrial code table. Gblocks was employed here to conduct the initial alignments and retain only the conserved regions with clear positional homology [30]. DnaSP v6.12.03 software [31] was then used to calculate the Ka and Ks of each PCG from mitochondrial genomes of M. strigata. Based on the selective pressure results, the PCG with highest variation, nad2, was selected as the marker gene for genetic analysis.

2.4. Amplification and Sequencing of Marker Gene

The total DNA of individuals from all seven populations was extracted from the adductor muscle using a TIANamp Marine Animals DNA kit (DP324-03, Tiangen Biotech Co., Ltd., Beijing, China). Then, PCR primers were designed using primer premier v5.5.0 software for the mitochondrial nad2 gene: forward primer 5′-GGTGATAAGAATAATAATAAGACCTTTG-3′, reverse primer 5′-CCACTCTTATGACCACTATTATC-3′. The PCR reaction mixture contained 0.5 µL of each primer, 1.0 µL of DNA, 12.5 µL of 2× Taq PCR MasterMix (PC1120; Beijing Solarbio Science & Technology Co., Ltd., Beijing, China), and 10.5 µL of ddH2O in a total volume of 25 µL. The PCR procedure involved initial denaturation at 95 °C for 3 min; 32 cycles of 95 °C for 30 s, 48 °C for 1 min, and 72 °C for 1 min; and a final extension at 72 °C for 5 min. Finally, agarose gel electrophoresis was used to detect the purity of PCR products and the Applied Biosystems 3730xl Genetic Analyzer (Beijing Tsingke Biotech Co., Ltd., Beijing, China) was used for sequencing.

2.5. Data Analysis

MEGA v7.0.26 software was employed to analyze the base composition and sequence variation of the nad2 gene sequence. DnaSP v6.12.03 [31] software was used to analyze the number of haplotypes (h), number of polymorphic (segregating) sites (S), haplotype (gene) diversity (Hd), nucleotide diversity (Pi), and the average number of nucleotide differences (K) of nad2 sequences from the seven populations. MEGA v7.0.26 software was used to align and evaluate the inter- and intra-population genetic distances of nad2 sequences from the seven populations with the Maximum Composite Likelihood model. The software of Arlequin v3.5 was used to perform the analysis of molecular variance (AMOVA) of the populations. PopART v1.7 [32] and Arlequin v3.5 [33] were used to construct haplotype networks. The genetic differentiation coefficient (F-statistics, FST) of inter-population and its significance, as well as Tajima’s D and Fu’s Fs, were also calculated using Arlequin v3.5 [33]. Finally, DnaSP v6.12.03 [31] software was used to perform mismatch distribution.

3. Results

3.1. Mitochondrial Genome of M. strigata

The mitochondrial genome measured in this study had a double-stranded circular structure, with GenBank accession No. OR666116. The mitochondrial genome consisted of 2 rRNA genes (12s and 16s), 22 tRNA genes, and 13 PCGs, including atp6, cox1, cox2, cox3, cytb, nad1, nad2, nad3, nad4, nad4l, nad5, nad6, and atp8.
The results of selective pressure analysis showed that all the PCGs had Ka/Ks < 1 (Table 2). This means that these genes exhibit a purifying selection. Among them, the selective pressure of nad2 was Ka/Ks = 0.15321, the maximum of these genes.

3.2. Base Composition and Sequence Variation of nad2

A total of 285 nad2 sequences were obtained from seven populations in this study and the sequence length was 691 bp. The percent identity of these sequences was found to be as high as 99.71% by aligning with M. strigata sequences in GenBank. The average contents of A, T, C, and G of the nad2 sequence were 13.67%, 34.88%, 41.44%, and 10.02%, respectively, exhibiting obviously T-rich and C-rich characteristics (Table 3). Totally, 31 variable sites were detected, accounting for 4.49% of the total number of sites. Among them, there were 18 parsimony informative sites and 13 singleton variable sites. In addition, there were two transition sites and one transversion site in mitochondrial nad2 sequences, and the average transition to transversion ratio was 3.9.

3.3. Population Genetic Diversity Analysis

All 285 nad2 gene sequences were used for genetic diversity analysis. Results showed that a total of 31 haplotypes (GenBank accession Nos. OR604589-OR604619) and 31 polymorphic (segregating) sites were charactered (Table 4). Among them, the ZJ population (19 h and 19 S) had the most haplotypes and polymorphic (segregating) sites, while the BH population (7 h) and QK population (9 S) had the fewest haplotypes and polymorphic (segregating) sites, respectively. The Hd, Pi, and K of all populations were 0.869, 0.00384, and 2.652, respectively. Similarly, the ZJ population had the highest Hd (0.898), Pi (0.00403), and K (2.783). The difference was that the BH population had the lowest Pi (0.00301) and K (2.077), while the GB population occupied the lowest Hd (0.712).
The population genetic distance and genetic differentiation analysis of M. strigata is shown in Table 5. The genetic distances between populations ranged from 0.00328 (BH and XM) to 0.00467 (ZZ and GB). And within populations, the BH population had the smallest genetic distance of 0.00303 and the ZJ population had the largest one of 0.00407. In addition, larger genetic differentiation (0.15~0.25) was observed at GB-ZJ (0.16360), GB-BH (0.15309), GB-SW (0.17637), and GB-ZZ (0.22865) (p < 0.05). In QK-ZJ (0.07061), QK-BH (0.05826), QK-SW (0.07391), QK-GB (0.13784), and GB-XM (0.14683), there was a moderate genetic differentiation (0.05~0.15) (p < 0.05). The results of AMOVA (Table 6) showed that 7.23% of the genetic variation occurred among populations and 92.77% occurred within populations.
The haplotype network based on the nad2 gene showed a “star” distribution with Hap_20 as the center (Figure 2). Among the 31 haplotypes, Hap_1, Hap_3, Hap_5, and Hap_6 were the dominant haplotypes with 196 individuals. Hap_2, Hap_3, and Hap_6 were shared by all the seven populations. Of all 20 unique haplotypes, the ZJ population possessed 10 haplotypes (Hap_16~Hap_19), while the remaining populations had another 10 haplotypes (Hap_11, Hap_12, and Hap_14) except the BH population. Significantly, the ZJ population had the most diverse and numerous haplotypes (19 and 43, respectively).

3.4. Neutrality Test and Mismatch Distribution Analysis

The results of the neutrality test (Table 7) showed that the Fu’s Fs values of all populations except the BH and QK populations were negative. The Fu’s Fs tests of the ZJ and SW populations were statistically significant (p < 0.05) and the mismatch distribution analysis (Figure 3) showed that the ZJ population displayed a unimodal distribution. The XM population also presented as a unimodal distribution but there was no statistical significance (p > 0.05). The mismatch distribution of the BH, ZZ, QK, and GB populations showed bimodal or multimodal distribution. In addition, the Tajima’s D values of the QK population were positive, while those of other populations were negative. There were no statistically significant values of Tajima’s D test in all populations (p > 0.05).

4. Discussion

4.1. Screening and Characterization of the Molecular Marker

Compared to the more conserved nuclear genome, the mitochondrial genome has a faster rate of genetic mutation. In primates and mollusks, the mutation rate of mtDNA genomes is about 5–10 times and 4 times than nuclear DNA, respectively [34,35]. In bivalves, although the ratio of the mtDNA to nuclear DNA mutation rate is lower (about 1.8), it is still higher than nuclear DNA [35]. At the same time, studies have shown that the mitochondrial mutation rate of Pteriomorphia, which M. strigata belongs to, is faster than that of other infraclasses in bivalves [36]. On the single gene level, the commonly used mitochondrial genes such as cox1 tend to have more genetic diversity than the nuclear genes such as its1 [22,37,38]. Due to the high rate of evolution, mtDNA could be a useful molecule for the analysis of evolutionary processes [39]. Therefore, mitochondrial genome analysis was performed on M. strigata taken from coastal areas of China and 13 PCGs were identified. Selective pressure analysis showed that the selective pressure of all 13 PCGs was Ka/Ks < 1, indicating that all the PCGs are subject to purifying selection. Purifying selection acts to weed out harmful mutations from the population and helps maintain the fitness of the species [40]. The maximum selective pressure (Ka/Ks = 0.15321) was found in nad2, suggesting that nad2 is subject to the least selective pressure and has the highest variation. Therefore, the nad2 gene was selected as the most appropriate genetic marker to analyze the population genetics of M. strigata. The only molecular marker used in the previous population genetic analysis of M. strigata was cox1 [7,21,22]. However, in contrast to cox1, the nad2 gene has no sex-related haplogroup, which is the second reason it was chosen.
In this study, the A + T base content was slightly lower than the G + C base content in the mitochondrial nad2 gene of M. strigata. The result showed a CG bias. Base bias is a characteristic of sequence evolution, which may result from natural mutation [41]. Although high AT base content is more common in invertebrate mtDNA [42], the nucleotide mutation rate of different fragments is inconsistent, which leads to different GC content in different parts [43]. Meanwhile, the difference in base content is species-specific [44]. In general, transition is more likely to occur than transversion because transition does not change the type of base [45]. In the results, the average transition to transversion ratio was 4.6. This means that the sequence variation of nad2 occurs mostly between a purine and a purine or between a pyrimidine and a pyrimidine.

4.2. Population Genetic Diversity and Structure of M. strigata

In the results, all the seven populations had Hd ≥ 0.5 and Pi < 0.005, suggesting that M. strigata has crossed the bottleneck effect in coastal China [46]. It is currently undergoing a population expansion with a rapid population growth and accumulation of mutations [47]. The invasion of M. strigata in China began around 2014 [7]. Zuo et al. hypothesized that M. strigata in China had experienced or was currently experiencing a population bottleneck during 2020~2021 [22], and this study further proves that the Chinese population of M. strigata in 2023 is in the stage of population expansion after the bottleneck effect. The phenomenon of the population expansion of M. strigata in China could also be seen in practical production. In 2019, a survey revealed that M. strigata has been present in farming ponds of Meretrix petechialis (Lamarck, 1818) in Taiwan and caused serious impact [10]. In 2020, M. strigata was found attached to an oyster farming facility in Hainan [7]. In 2022, the density of M. strigata surpassed that of the native species P. viridis in Donghai island [9]. In addition, the ZJ population, where M. strigata was first found in mainland China, has the highest genetic diversity. This corresponds to the greatest genetic distance within the ZJ population
Based on nad2 sequences, the genetic distance between populations of M. strigata ranged from 0.00328 to 0.00467. The genetic distance between the BH population and the XM population was the nearest, indicating a closer relationship [48], while that of the ZZ population and the GB population was the furthest. AMOVA analysis results showed that the genetic variation of M. strigata was mainly from within the population rather than between the populations. This was also evidenced by genetic differentiation. It is generally believed that genetic differentiation between populations is small, moderate, and larger when FST is 0~0.05, 0.15~0.25, and 0.15~0.25, respectively [49]. Except for the QK and GB populations, the genetic differentiation among other populations was not obvious. This may be related to the geographical location of the QK and GB populations on Hainan Island. In addition, it is found that Hap_20 acts as the center and diverges to other haplotypes according to the haplotype network. Therefore, Hap_20 is presumed to be an ancient haplotype [50]. The “star” distribution of haplotype networks also supported that the different populations of M. strigata began to expand after experiencing bottleneck [51]. However, the ZJ population and the ZZ population do not exist in the ancient haplotype Hap_20. From this, we hypothesize that the ZJ population and the ZZ population may have been introduced by other populations. This might also explain why there was a large genetic distance within the ZJ population and a small genetic distance between the ZJ population and other populations. The ZJ and SW populations had the most unique haplotypes, indicating that these two populations could have undergone adaptive evolution [7].

4.3. Population Historical Dynamics of M. strigata

Neutrality tests and mismatch distribution are two effective tests to study the historical dynamics of populations [52]. A population is identified as stable when the value of the neutrality test approaches 0. In another case, a population is considered to have expanded in history when the value is negative and statistically significant [53]. In our results, Fu’s Fs and Tajima’s D values of all populations except the QK and BH populations were negative. However, only the ZJ population and the SW population had statistically significant Fu’s Fs values. Tajima’s D test is more likely to reveal the history of expansion of an ancient population, while Fu’s Fs test is more sensitive to the recent expansion of a population [54,55]. Therefore, it is speculated that there have been recent population expansion events in the ZJ population and the SW population. For a stable population, mismatch distribution analyses will display bimodal or multimodal distribution. However, if a population has experienced expansion or sustained growth, the results will show a unimodal distribution. According to the results, the ZJ and XM populations showed a unimodal distribution, but the neutrality test of the XM population was not significant. This once again supported that recent expansion and sustained growth had occurred in the ZJ population [56].

5. Conclusions

In this study, gene nad2 was ultimately selected as the best mitochondrial marker for population genetic analysis by assessing the selective pressure of 13 mitochondrial PCGs. Population genetic analysis showed that the seven M. strigata populations in China could be expanding after the bottleneck. Among them, recent expansion and sustained growth events had occurred in the ZJ population. Therefore, relevant measures, including removal (physical removal and biological control) and utilization (human food and feed), need to be taken to prevent the formation of large and stable invasive populations of M. strigata in the future.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes14112038/s1, Figure S1: Shell morphology of Mytella strigata.

Author Contributions

Investigation, C.Z.; original draft, S.Y. and P.M.; collected samples, Y.Z. and X.M.; data analysis, S.Y.; reviewed and edited the manuscript, P.M. and Z.Z.; in charge of the methodology, funding provision, Z.Z. and P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China: 2022YFD2401204 and National Natural Science Foundation of China: 42006080.

Institutional Review Board Statement

Ethical review and approval were waived for this study because the invasive mussel in this study is an invertebrate with no sense or subjective experience.

Informed Consent Statement

Not applicable.

Data Availability Statement

Relevant information has been added in the article.

Acknowledgments

We would like to thank the reviewers for their constructive comments and helpful suggestions that improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Teem, J.L.; Alphey, L.; Descamps, S.; Edgington, M.P.; Edwards, O.; Gemmell, N.; Harvey Samuel, T.; Melnick, R.L.; Oh, K.P.; Piaggio, A.J.; et al. Genetic Biocontrol for Invasive Species. Front. Bioeng. Biotechnol. 2020, 8, 452. [Google Scholar] [CrossRef] [PubMed]
  2. Sherpa, S.; Després, L. The evolutionary dynamics of biological invasions: A multi-approach perspective. Evol. Appl. 2021, 14, 1463–1484. [Google Scholar] [CrossRef]
  3. Sax, D.F.; Stachowicz, J.J.; Brown, J.H.; Bruno, J.F.; Dawson, M.N.; Gaines, S.D.; Grosberg, R.K.; Hastings, A.; Holt, R.D.; Mayfield, M.M.; et al. Ecological and evolutionary insights from species invasions. Trends Ecol. Evol. 2007, 22, 465–471. [Google Scholar] [CrossRef]
  4. Renault, D.; Hess, M.C.M.; Braschi, J.; Cuthbert, R.N.; Sperandii, M.G.; Bazzichetto, M.; Chabrerie, O.; Thiébaut, G.; Buisson, E.; Grandjean, F.; et al. Advancing biological invasion hypothesis testing using functional diversity indices. Sci. Total Environ. 2022, 834, 155102. [Google Scholar] [CrossRef] [PubMed]
  5. Molnar, J.L.; Gamboa, R.L.; Revenga, C.; Spalding, M.D. Assessing the global threat of invasive species to marine biodiversity. Front. Ecol. Environ. 2008, 6, 485–492. [Google Scholar] [CrossRef]
  6. Lim, C.S.; Tay, T.S.; Tan, K.S.; Teo, S.L. Removal of larvae of two marine invasive bivalves, Mytilopsis sallei (Récluz, 1849) and Mytella strigata (Hanley, 1843), by water treatment processes. Mar. Pollut. Bull. 2020, 155, 111154. [Google Scholar] [CrossRef]
  7. Ma, P.Z.; Li, H.M.; Liu, Y.M.; Li, C.; Zhang, Z.; Wang, H.Y. First confirmed occurrence of the invasive mussel Mytella strigata (Hanley, 1843) in Guangdong and Hainan, China, and its rapid spread in Indo-West Pacific regions. BioInvasions Rec. 2022, 11, 947–963. [Google Scholar] [CrossRef]
  8. Xu, X.; Zhang, X.; Peng, J.; Deng, Y.; Liu, Y.; Jiang, L.; Zhao, L. Survival and physiological energetics of highly invasive mussels exposed to heatwaves. Mar. Environ. Res. 2023, 187, 105948. [Google Scholar] [CrossRef]
  9. Ma, P.Z.; Zuo, C.X.; Li, H.M.; Wang, H.Y.; Wang, Q.H.; Zhang, Z. The basic biology and biological invasion of Mytella strigata. Acta Ecol. Sin. 2023, 43, 5251–5259. [Google Scholar] [CrossRef]
  10. Huang, Y.C.; Li, Z.K.; Chen, W.L.; Chan, C.C.; Hsu, H.Y.; Lin, Y.T.; Huang, Y.S.; Han, Y.S. First record of the invasive biofouling mussel Mytella strigata (Hanley, 1843) (Bivalvia: Mytilidae) from clam ponds in Taiwan. BioInvasions Rec. 2021, 10, 304–312. [Google Scholar] [CrossRef]
  11. Sanpanichl, K.; Wells, F.E. Mytella strigata (Hanley, 1843) emerging as an invasive marine threat in Southeast Asia. BioInvasions Rec. 2019, 8, 343–356. [Google Scholar] [CrossRef]
  12. Dahlberg, A.D.; Waller, D.L.; Hammond, D.; Lund, K.; Phelps, N.B.D. Open water dreissenid mussel control projects: Lessons learned from a retrospective analysis. Sci. Rep. 2023, 13, 10410. [Google Scholar] [CrossRef]
  13. Li, H.Y. Genetic Diversity and Differentiation among Populations of Flaveria bidentis in Hebei province. Master’s Thesis, Hebei Agricultural University, Baoding, China, June 2009. [Google Scholar]
  14. Jarić, I.; Heger, T.; Castro Monzon, F.; Jeschke, J.M.; Kowarik, I.; McConkey, K.R.; Pyšek, P.; Sagouis, A.; Essl, F. Crypticity in Biological Invasions. Trends Ecol. Evol. 2019, 34, 291–302. [Google Scholar] [CrossRef] [PubMed]
  15. White, D.J.; Wolff, J.N.; Pierson, M.; Gemmell, N.J. Revealing the hidden complexities of mtDNA inheritance. Mol. Ecol. 2008, 17, 4925–4942. [Google Scholar] [CrossRef]
  16. Birky, C.W. The inheritance of genes in mitochondria and chloroplasts: Laws, mechanisms, and models. Annu. Rev. Genet. 2001, 35, 125–148. [Google Scholar] [CrossRef]
  17. Ghiselli, F.; Gomes Dos Santos, A.; Adema, C.M.; Lopes Lima, M.; Sharbrough, J.; Boore, J.L. Molluscan mitochondrial genomes break the rules. Philos. Trans. R. Soc. London Ser. B. 2021, 376, 20200159. [Google Scholar] [CrossRef] [PubMed]
  18. Sato, K.; Sato, M. Multiple ways to prevent transmission of paternal mitochondrial DNA for maternal inheritance in animals. J. Biochem. 2017, 162, 247–253. [Google Scholar] [CrossRef]
  19. Ghiselli, F.; Milani, L.; Passamonti, M. Strict sex-specific mtDNA segregation in the germ line of the DUI species Venerupis philippinarum (Bivalvia: Veneridae). Mol. Biol. Evol. 2011, 28, 949–961. [Google Scholar] [CrossRef]
  20. Obata, M.; Kamiya, C.; Kawamura, K.; Komaru, A. Sperm mitochondrial DNA transmission to both male and female offspring in the blue mussel Mytilus galloprovincialis. Dev. Growth Differ. 2006, 48, 253–261. [Google Scholar] [CrossRef]
  21. Alves, F.A.d.S.; Beasley, C.R.; Hoeh, W.R.; Rocha, R.M.D.; Tagliaro, C.H. Detection of mitochondrial DNA heteroplasmy suggests a doubly uniparental inheritance pattern in the mussel Mytella charruana. Braz. J. Biosci. 2012, 10, 176–185. [Google Scholar]
  22. Zuo, C.X.; Zhang, T.T.; Zhang, C.C.; Zhao, D.T.; Zhu, Y.; Ma, X.J.; Wang, H.Y.; Ma, P.Z.; Zhang, Z. Evaluating the utility of five gene fragments for genetic diversity analyses of Mytella strigata populations. Fishes 2023, 8, 34. [Google Scholar] [CrossRef]
  23. Ngandu, N.K.; Scheffler, K.; Moore, P.; Woodman, Z.; Martin, D.; Seoighe, C. Extensive purifying selection acting on synonymous sites in HIV-1 Group M sequences. Virol. J. 2008, 5, 160. [Google Scholar] [CrossRef] [PubMed]
  24. Lim, J.Y.; Tay, T.S.; Lim, C.S.; Lee, S.S.C.; Teo, S.L.M.; Tan, K.S. Mytella strigata (Bivalvia: Mytilidae): An alien mussel recently introduced to Singapore and spreading rapidly. Molluscan Res. 2018, 38, 170–186. [Google Scholar] [CrossRef]
  25. Coil, D.; Jospin, G.; Darling, A.E. A5-miseq: An updated pipeline to assemble microbial genomes from Illumina MiSeq data. Bioinformatics 2015, 31, 587–589. [Google Scholar] [CrossRef] [PubMed]
  26. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  27. Donath, A.; Jühling, F.; Al-Arab, M.; Bernhart, S.H.; Reinhardt, F.; Stadler, P.F.; Middendorf, M.; Bernt, M. Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes. Nucleic Acids Res. 2019, 47, 10543–10552. [Google Scholar] [CrossRef]
  28. Zhang, D.; Gao, F.; Jakovlić, I.; Zou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 2020, 20, 348–355. [Google Scholar] [CrossRef]
  29. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef]
  30. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef]
  31. Rozas, J.; Ferrer Mata, A.; Sánchez DelBarrio, J.C.; Guirao Rico, S.; Librado, P.; Ramos Onsins, S.E.; Sánchez Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  32. Leigh, J.W.; Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  33. Excoffier, L.; Lischer, H.E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef] [PubMed]
  34. Brown, W.M.; Prager, E.M.; Wang, A.; Wilson, A.C. Mitochondrial DNA sequences of primates: Tempo and mode of evolution. J. Mol. Evol. 1982, 18, 225–239. [Google Scholar] [CrossRef] [PubMed]
  35. Allio, R.; Donega, S.; Galtier, N.; Nabholz, B. Large Variation in the Ratio of Mitochondrial to Nuclear Mutation Rate across Animals: Implications for Genetic Diversity and the Use of Mitochondrial DNA as a Molecular Marker. Mol. Biol. Evol. 2017, 34, 2762–2772. [Google Scholar] [CrossRef]
  36. Mortz, M.; Levivier, A.; Lartillot, N.; Dufresne, F.; Blier, P.U. Long-Lived Species of Bivalves Exhibit Low MT-DNA Substitution Rates. Front. Mol. Biosci. 2021, 8, 626042. [Google Scholar] [CrossRef]
  37. Lee, Y.; Ni, G.; Shin, J.; Kim, T.; Kern, E.M.A.; Kim, Y.; Kim, S.C.; Chan, B.; Goto, R.; Nakano, T.; et al. Phylogeography of Mytilisepta virgata (Mytilidae: Bivalvia) in the northwestern Pacific: Cryptic mitochondrial lineages and mito-nuclear discordance. Mol. Phylogenet. Evol. 2021, 157, 107037. [Google Scholar] [CrossRef]
  38. Peretolchina, T.E.; Bukin Iu, S.; Sitnikova, T.; Shcherbakov, D. Genetic differentiation of the endemic Baikalian mollusk Baicalia carinata (Mollusca: Caenogastropoda). Genetika 2007, 43, 1667–1675. [Google Scholar] [CrossRef]
  39. Brown, W.M.; George, M.; Wilson, A.C. Rapid evolution of animal mitochondrial DNA. Proc. Natl. Acad. Sci. USA 1979, 76, 1967–1971. [Google Scholar] [CrossRef]
  40. De Silva, D.R.; Nichols, R.; Elgar, G. Purifying selection in deeply conserved human enhancers is more consistent than in coding sequences. PLoS ONE 2014, 9, e103357. [Google Scholar] [CrossRef]
  41. Li, Z.X.; Li, J.Q.; Xu, W.J.; Fan, L.X. Mitochondrial Genome Characterization and Phylogeny of Diplodiscus nigromaculati. Sichuan J. Zool. 2022, 41, 248–255. [Google Scholar] [CrossRef]
  42. Mao, Y.L. Studies on molecular phylogenetics of genus mytilus and population morphology and genetics of three mussels. Master’s Thesis, Ocean University of China, Qingdao, China, June 2011. [Google Scholar]
  43. Zhong, D.; Zhao, G.J.; Zhang, Z.S.; Xun, A.L. Advance in the Entire Balance and Local Unbalance of Base Distribution in Genome. Hereditas 2002, 24, 351–355. [Google Scholar] [CrossRef] [PubMed]
  44. Sueoka, N. On the genetic basis of variation and heterogeneity of DNA base composition. Proc. Natl. Acad. Sci. USA 1962, 48, 582–592. [Google Scholar] [CrossRef]
  45. Stoltzfus, A.; Norris, R.W. On the Causes of Evolutionary Transition: Transversion Bias. Mol. Biol. Evol. 2016, 33, 595–602. [Google Scholar] [CrossRef] [PubMed]
  46. Ren, Y.L.; Dai, J.C.; Zhao, N.H.; Zhao, H.; Wei, J.; Nie, Z.L. Genetic diversity analysis of Schizothorax biddulphi based on mitochondrial DNA COII and ND4 genes. J. Fish. China 2022, 46, 2274–2285. [Google Scholar] [CrossRef]
  47. Grant, W.S.; Bowen, B.W. Shallow population histories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation. J. Hered. 1998, 89, 415–426. [Google Scholar] [CrossRef]
  48. Yin, X.Y.; Chen, Y.T.; Zhuang, E.J.; Zhao, J.J.; Dong, L.; Li, H.L. Genetic diversity analysis of Taenia asiatica in Dali based on mitochondrial nad1 gene. J. Trop. Med. 2023, 23, 301–304. [Google Scholar] [CrossRef]
  49. Zheng, Z.H.; Wu, S.S.; Wei, Y.; Zhong, D.B.; Zheng, X.L. Analysis of the genetic diversity of 15 Aedes albopictus populations in Guangzhou based on the mitochondrial COI gene. Chin. J. Zoonoses. 2021, 37, 985–994. [Google Scholar] [CrossRef]
  50. Guo, Q.; Zhu, Q.D.; Zhou, Z.J.; Shi, F.M. Phylogeography and Genetic Structure of the Bush Cricket Decma fissa (Orthoptera, Tettigoniidae) in Southern China. Zool. Stud. 2023, 62, e32. [Google Scholar] [CrossRef]
  51. Slatkin, M.; Hudson, R.R. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 1991, 129, 555–562. [Google Scholar] [CrossRef]
  52. Mousset, S.; Derome, N.; Veuille, M. A test of neutrality and constant population size based on the mismatch distribution. Mol. Biol. Evol. 2004, 21, 724–731. [Google Scholar] [CrossRef]
  53. Guo, J.; Li, J.; Massart, S.; He, K.; Francis, F.; Wang, Z. Analysis of the Genetic Diversity of Two Rhopalosiphum Species from China and Europe Based on Nuclear and Mitochondrial Genes. Insects 2023, 14, 57. [Google Scholar] [CrossRef] [PubMed]
  54. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar] [CrossRef]
  55. Fu, Y.X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 1997, 147, 915–925. [Google Scholar] [CrossRef] [PubMed]
  56. Rogers, A.R.; Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 1992, 9, 552–569. [Google Scholar] [CrossRef] [PubMed]
Genes 14 02038 g001
Figure 1. Sampling station maps for the seven populations of Mytella strigata. Black circles represent sampling points. The size of the black circle indicates the number of samples.
Figure 1. Sampling station maps for the seven populations of Mytella strigata. Black circles represent sampling points. The size of the black circle indicates the number of samples.
Genes 14 02038 g001
Genes 14 02038 g002
Figure 2. Haplotype TCS networks in seven populations of Mytella strigata based on mitochondrial nad2 gene. The size of the circle represents the frequency of the haplotype. Different colors represent different populations. The number of segments of the connecting lines between haplotypes represents the number of nucleotide replacements.
Figure 2. Haplotype TCS networks in seven populations of Mytella strigata based on mitochondrial nad2 gene. The size of the circle represents the frequency of the haplotype. Different colors represent different populations. The number of segments of the connecting lines between haplotypes represents the number of nucleotide replacements.
Genes 14 02038 g002
Genes 14 02038 g003
Figure 3. Mismatch distribution of Mytella strigata populations. The bar chart represents the distribution observed, and the line chart represents the distribution predicted by the model. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen.
Figure 3. Mismatch distribution of Mytella strigata populations. The bar chart represents the distribution observed, and the line chart represents the distribution predicted by the model. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen.
Genes 14 02038 g003
Table 1. Information about the sampling stations. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen.
Table 1. Information about the sampling stations. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen.
LocationCollection DateGeographical CoordinatesQuantity
ZJ2 April 2023110°31′43″ E, 20°55′53″ N43
BH3 April 2023109°11′36″ E, 21°24′42″ N38
SW5 April 2023115°36′28″ E, 22°51′25″ N45
ZZ6 April 2023117°58′0″ E, 24°9′14″ N38
QK1 April 2023110°35′7″ E, 19°59′7″ N40
GB1 April 2023110°30′50″ E, 18°53′12″ N43
XM6 April 2023118°4′41″ E, 24°28′1″ N38
Total -285
Table 2. The evolutionary constraint (Ka/Ks) analyses of 13 mitochondrial protein-coding genes (PCGs). Ka: nonsynonymous substitution rate, Ks: synonymous substitution rate calculations.
Table 2. The evolutionary constraint (Ka/Ks) analyses of 13 mitochondrial protein-coding genes (PCGs). Ka: nonsynonymous substitution rate, Ks: synonymous substitution rate calculations.
Mitochondrial PCGsKaKsKa/Ks
atp60.001260.012330.10219
atp80.000000.037500.00000
cox10.001130.016140.07001
cox20.002790.00000-
cox30.000000.013510.00000
cytb0.000000.003750.00000
nad10.000000.014530.00000
nad20.000930.006070.15321
nad30.000000.016330.00000
nad40.001340.019380.06914
nad4l0.000000.00000-
nad50.002020.016500.12242
nad60.001870.017440.10722
Table 3. Base composition of nad2 gene sequence by population (%). ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen. A: Adenine, T: Thymine, C: Cytosine, G: Guanine.
Table 3. Base composition of nad2 gene sequence by population (%). ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen. A: Adenine, T: Thymine, C: Cytosine, G: Guanine.
Nucleotide CompositionATCG
ZJ13.5034.6641.2610.58
BH13.3335.3141.519.85
SW13.1935.5641.0810.18
ZZ13.5035.4041.269.85
QK13.9534.0342.269.76
GB14.7834.6640.749.83
XM13.4134.5541.9410.10
Average13.6734.8841.4410.02
Table 4. Genetic diversity analysis based on nad2 molecular markers. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen. h, number of haplotypes; S, number of polymorphic (segregating) sites; Hd, haplotype (gene) diversity; Pi, nucleotide diversity; K, average number of nucleotide differences.
Table 4. Genetic diversity analysis based on nad2 molecular markers. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen. h, number of haplotypes; S, number of polymorphic (segregating) sites; Hd, haplotype (gene) diversity; Pi, nucleotide diversity; K, average number of nucleotide differences.
hSHdPiK
ZJ19190.8980.004032.783
BH7100.8110.003012.077
SW14180.8670.003772.606
ZZ8110.7720.003452.385
QK890.8040.003642.517
GB9130.7120.003682.545
XM11150.8820.003492.410
Total31310.8690.003842.652
Table 5. Genetic distance (diagonal and below diagonal) and genetic differentiation coefficient (F-statistics, FST) (above diagonal) between populations of Mytella strigata based on nad2 gene sequences. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen. A: Adenine, T: Thymine, C: Cytosine, G: Guanine.
Table 5. Genetic distance (diagonal and below diagonal) and genetic differentiation coefficient (F-statistics, FST) (above diagonal) between populations of Mytella strigata based on nad2 gene sequences. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen. A: Adenine, T: Thymine, C: Cytosine, G: Guanine.
ZJBHSWZZQKGBXM
ZJ0.00407−0.00912−0.003600.015630.07061 *0.16360 *0.01903
BH0.003510.00303−0.001920.024290.05826 *0.15309 *0.00331
SW0.003920.003410.003800.003820.07391 *0.17637 *0.01665
ZZ0.003830.003330.003650.003480.04562 *0.22865 *0.01477
QK0.004170.003560.004040.003750.003670.13784 *0.00103
GB0.004650.003990.004570.004670.004290.003720.14683 *
XM0.003860.003280.003720.003550.003600.004240.00351
* indicates significant level p < 0.05.
Table 6. AMOVA results of Mytella strigata among seven populations. d.f.: degrees of freedom.
Table 6. AMOVA results of Mytella strigata among seven populations. d.f.: degrees of freedom.
Source of Variationd.f.Sum of SquaresVariance ComponentsPercentage of Variation
Among populations631.1100.09690 Va7.23
Within populations278345.4241.24253 Vb92.77
Table 7. Fu’s Fs and Tajima’s D test from different populations of Mytella strigata. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen.
Table 7. Fu’s Fs and Tajima’s D test from different populations of Mytella strigata. ZJ: ZhangJiang, BH: Beihai, SW: Shanwei, ZZ: Zhangzhou, QK: Qukou, GB: Gangbei, XM: Xiamen.
Fu’s Fsp-ValueTajima’s Dp-Value
ZJ−10.533990.000 *−1.189430.114
BH0.050770.531−0.383890.410
SW−4.452610.028 *−1.174570.128
ZZ−0.255690.495−0.272260.419
QK0.019030.5410.553550.737
GB−0.537390.439−0.471090.376
XM−2.589370.101−1.045820.167
* indicates significant level p < 0.05. Tajima’s D and Fu’s Fs are both neutrality test values.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, S.; Ma, P.; Zuo, C.; Zhu, Y.; Ma, X.; Zhang, Z. Genetic Analysis Based on Mitochondrial nad2 Gene Reveals a Recent Population Expansion of the Invasive Mussel, Mytella strigata, in China. Genes 2023, 14, 2038. https://doi.org/10.3390/genes14112038

AMA Style

Yan S, Ma P, Zuo C, Zhu Y, Ma X, Zhang Z. Genetic Analysis Based on Mitochondrial nad2 Gene Reveals a Recent Population Expansion of the Invasive Mussel, Mytella strigata, in China. Genes. 2023; 14(11):2038. https://doi.org/10.3390/genes14112038

Chicago/Turabian Style

Yan, Shaojing, Peizhen Ma, Chenxia Zuo, Yi Zhu, Xiaojie Ma, and Zhen Zhang. 2023. "Genetic Analysis Based on Mitochondrial nad2 Gene Reveals a Recent Population Expansion of the Invasive Mussel, Mytella strigata, in China" Genes 14, no. 11: 2038. https://doi.org/10.3390/genes14112038

APA Style

Yan, S., Ma, P., Zuo, C., Zhu, Y., Ma, X., & Zhang, Z. (2023). Genetic Analysis Based on Mitochondrial nad2 Gene Reveals a Recent Population Expansion of the Invasive Mussel, Mytella strigata, in China. Genes, 14(11), 2038. https://doi.org/10.3390/genes14112038

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