Next Article in Journal
Combined Use of Job Stress Models and the Incidence of Glycemic Alterations (Prediabetes and Diabetes): Results from ELSA-Brasil Study
Previous Article in Journal
Developing Youth Environmental Health Literacy and Civic Leadership through Community Air Monitoring in Imperial County, California
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 17
Issue 5
10.3390/ijerph17051538
ijerph-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
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Population Genetic Structure Reveals Two Lineages of Amynthas triastriatus (Oligochaeta: Megascolecidae) in China, with Notes on a New Subspecies of Amynthas triastriatus

by
Yan Dong
1,2,
Jibao Jiang
1,2,
Zhu Yuan
1,2,
Qi Zhao
1,2 and
Jiangping Qiu
1,2,*
1
School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai 200240, China
2
Key Laboratory of Urban Agriculture (South), Ministry of Agriculture, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(5), 1538; https://doi.org/10.3390/ijerph17051538
Submission received: 21 January 2020 / Revised: 22 February 2020 / Accepted: 24 February 2020 / Published: 27 February 2020
(This article belongs to the Section Environment and Applied Ecology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Amynthas triastriatus (Oligochaete: Megascolecidae) is a widely distributed endemic species in Southern China. To shed light on the population genetic diversity and to elucidate the population differentiation and dispersal of A. triastriatus, a population genetic structure study was undertaken based on samples from 35 locations collected from 2010 to 2016. Two exclusive lineages within A. triastriatus—lineage A and lineage B—were revealed. Lineage A was mainly distributed at high altitudes while lineage B was mainly distributed at low altitudes in Southeast China. The genetic diversity indices indicated that the populations of A. triastriatus had a strong genetic structure and distinct dispersal histories underlying the haplogroups observed in this study. Combined with morphological differences, these results indicated a new cryptic subspecies of A. triastriatus. Lineage A was almost degenerated to parthenogenesis and lineage B had a trend to parthenogenesis, which suggested that parthenogenesis could be an internal factor that influenced the differentiation and dispersal of A. triastriatus. The divergence time estimates showed that A. triastriatus originated around Guangxi and Guangdong provinces and generated into two main lineages 2.97 Ma (95%: 2.17–3.15 Ma) at the time of Quaternary glaciation (2.58 Ma), which suggested that the Quaternary glaciation may have been one of main factors that promoted the colonization of A. triastriatus.

1. Introduction

Understanding the evolutionary history of a species is a fundamental goal of evolutionary biology [1] and many approaches are based on population genetic structure analysis, both in fauna and flora [2,3]. Population genetic structure studies combined with geographical distribution, time of differentiation and ancient geographical events are commonly used to determine the differentiation and dispersal of populations and to analyze their influencing factors.
Earthworms are very important macro-invertebrate detritivores that tend to remain in the same areas during long periods of time and show low ability to cross mountains; for these reasons they are often an ideal model for phylogeographical study [4,5,6,7,8]. Earthworms play important roles in the agroecosystem as key organisms that influence soil structure formation, soil carbon dynamics and biogeochemical cycles [9,10,11], and hence failure to recognize accurate species boundaries within earthworms compromises many aspects of applied ecological, biodiversity, systematic, and evolutionary studies [12,13,14].
Their structural simplicity makes it difficult to distinguish earthworm species accurately using only morphological characteristics [15]. Molecular tools could be helpful to solve this problem, and analysis of population genetic structure and colonization history could reveal differentiation and dispersal of species and provide a solid base for taxonomy. Recent studies on this aspect for earthworms have expanded rapidly. Phylogenetic analysis using different substitution rates for the CO1 gene was carried out to illustrate the differentiation and diffusion of earthworms [6,7,16]. Parthenogenesis, geological events, climatic changes and human beings are the factors that exert significant influence on the differentiation and diffusion of earthworms [16,17,18,19,20]. Additionally, phylogenetic analysis constructed by single gene or different gene combinations have led to an-increasing discovery of cryptic earthworm species [6,12,15,21,22], which are taxonomically important to earthworms.
In recent research, mitochondrial COI, COII, NDI, 12S rRNA and 16S rRNA have been widely used to study the phylogenetic evaluation, mitochondrial diversity and population genetic diversity of earthworm using multivariate concatenation [18,21,23,24,25,26,27,28,29]. Complete mitochondrial genomes and genomic features of four pheretimoid earthworms (Amynthas longisiphonus, Amynthas corticis, Amynthas gracilis and Metaphire californica) revealed that mitochondrial ATP8 exhibited the highest evolutionary rate and the NADH gene was second to it [30]. Nuclear 18S rRNA has been more helpful within recognized families, rather than to elucidate deep annelid relationships [31,32]. The entire 28S rRNA appears suitable for within-genus and sometimes even within-species complexity if the sequences are analyzed with other genes [15,17,24,33]. The focus of this study was the intraspecies level using mitochondrial genes.
Amynthas triastriatus (Chen 1946) has two pairs of spermathecal pores in segments 7/8–8/9 and belongs to the aeruginosus-group. It is named after the triangle formed by three genital papillae in each male pore zone and is a widely distributed endemic species in Southern China [34]. There are many populations of A. triastriatus that are similar in morphology due to a tendency to parthenogenesis; for example, some populations have a thin and lustreless seminal chamber and prostate glands absent, and others have a plump and glossy seminal chamber and small prostate glands. As a result, their classification status and the genetic relationships are obscure.
A previous study based on only seven individuals indicated that the ancestor of A. triastriatus originated in the vicinity of Guangxi and Guangdong provinces then split into two lineages around 4.31 Ma (3.24–5.37 Ma) [34]—this was a little earlier than the Quaternary glaciation in China. Therefore, we assume that A. triastriatus presents two lineages, one displays parthenogenesis and the other is degenerate towards parthenogenesis. This raised general questions as whether there are cryptic species, whether the population genetic structure are useful to clear the classification status of A. triastriatus and whether genetic isolation occurred between the two lineages due to population expansion and reduction along with ancient geographical events. In the present study, a fragment of mitochondrial gene COI of 65 individuals collected from 35 locations and one Aporrectodea species, and 4691bp segments comprising mitochondrial genes COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI from 35 locations and one Aporrectodea species were used to determine: (1) the population genetic diversity and the classification status of A. triastriatus; (2) the population differentiation and dispersal of A. triastriatus; (3) the time of divergence and colonization histories of A. triastriatus, as well the main factors impacting these.

2. Materials and Methods

2.1. Earthworm Sampling

A total of 232 individuals of A. triastriatus were collected by digging at 35 different locations (Figure 1, Table 1) from 17 provinces in China from 2010 to 2016. One individual of Aporrectodea trapezoides was collected as outgroup in Yunnan Province, China. Collected individuals were washed by clean water, anaesthetized in 10% ethanol solution, then fixed in anhydrous alcohol and preserved at 4 °C. The tail muscle issue (±0.02 g) was removed and cleaned carefully under a stereomicroscope for DNA extraction, and the rest of the individual was used for taxonomic identification.

2.2. DNA Extraction, Amplification and Sequencing

Total genomic DNA was extracted using the E.Z.N.A.® Mollusc DNA kit (Omega Bio-Tek, Norcross, GA, USA) according to the manufacturer’s instructions. The quality of DNA samples was checked by electrophoresis using a 1% agarose gel and stored at −20 °C.
The mitochondrial genes were amplified using 50 μL of reaction mixture with 1 μL template DNA, 0.6 μL Trans Taq® HIFI DNA Polymerase (Trans, Beijing, China), 4 μL dNTP, 5 μL HiFi Buffer, 2 μL primer and 35.4 μL ddH2O. The PCR of genes COI, COII, ATP8, 12S rRNA, 16S rRNA and NDI were carried out for 5 min at 94 °C followed by 32 cycles of 94 °C for 30 s, 50 °C for 30 s and 72 °C for 1 min, with an extension of 10 min at 72 °C. The PCR of gene ND6 was carried out for 5 min at 94 °C followed by 40 cycles of 94 °C for 30 s, 51.9 °C for 30 s and 72 °C for 1 min, with an extension of 10 min at 72 °C. The primers of mitochondrial COI, COII, 12S rRNA, 16S rRNA and NDI are listed in Table 2. The primers of genes ATP8 and ND6 were designed in this paper.
The PCR products were visualized on a 1% agarose gel and subsequently purified then sequenced by the Beijing Genomics Institute (Shanghai, China) using an ABI 3730 DNA analyzer. The sequences were compared with known earthworm sequences in GenBank using the BLAST search algorithm [35,36,37,38].

2.3. Sequence Alignment, Population Genetic Structure and Genetic Diversity, Time of Divergence and Colonization History

Population genetic structure and genetic diversity studies were conducted using single-gene datasets of COI with one Aporrectodea species as an outgroup. Time of divergence and colonization history were conducted using a combination of COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI datasets with one Aporrectodea species as an outgroup. Sequences of COI and a combined dataset were aligned using the Clustal X 1.8 program [39].
Models for sequence evolution and corresponding parameters were estimated using jModeltest v.2.1.7 [40]. The best-fitting models were used in the phylogenetic estimation for the maximum likelihood (ML) and Bayesian phylogeny estimation (BI) analyses. jModeltest revealed TPM3uf+G to be the best-fitting model for the COI gene (679 bp) with a gamma shape parameter of 0.2360, and GTR+G to be the best-fitting model for 4691 bp segments comprising mitochondrial genes COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI with a gamma shape parameter of 0.3950.
The phylogenetic trees were reconstructed using different methods based on single-gene COI datasets of 65 A. triastriatus individuals and one individual of Aporrectodea trapezoides as outgroup. ML analyses were performed using PhyML v.3.0 [41]; clade support was assessed using bootstrapping with 1000 pseudo replicates. BI was performed using MRBAYES v.3.2 [42]; the parameters were set to 10,000,000 with trees sampled every 1000 generations, discarding the first 10% as burn-in. Posterior probabilities and bootstrap support for each branch were calculated from the sampled trees.
Population structure was evaluated in Arlequin v.3.5.2 [43] via AMOVA with 10,000 permutations for statistical confidence. A haplotype network was constructed using the median joining (MJ) method Network v.4.6.1.6 [44] to infer the relationships among haplotypes and their geographical distribution. Nucleotide mismatch distribution analyses of lineages A and B were calculated with DnaSP 5.0, and neutrality tests were used to test population equilibrium.
Basic statistics of mtDNA diversity, including nucleotide and haplotype diversity, Tajima’s D [45] and Fu’s Fs [46] neutrality tests and nucleotide mismatch distribution analysis with and within lineages were calculated with DnaSP 5.0 [47]. Bayesian Skyline Plot (BSP) was constructed with BEAST 1.8.2 [48] and Tracer 1.7 [49] to infer the timing of the population events.
The Bayesian tree and the time of divergence of A. triastriatus were estimated using BEAST 1.8.2 [48] based on 4691 bp segments comprising mitochondrial genes COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI. These were from 18 individuals from 18 locations of lineage A, 17 individuals from 17 locations of lineage B and one individual of Aporrectodea trapezoides as an outgroup. A previous study has estimated 2.1 My−1 as the average substitution rate for COI, COII, ATP8, ND6, NDI, 12S and 16S genes [50], and this value was applied in the present study. Since earthworm lack a fossil record, the age of the calibration clade between A. triastriatus and Aporrectodea trapezoides was set at a mean age of 212 Ma (95%: 200–230 Ma); this represented the minimum age of Lumbricoidea and Megascolecoidea in our tree and is considered to be when the breakup and drift of Gondwana occurred [50,51,52,53]. Five replicate runs were performed with 10,000,000 MCMC steps, and the first 10% served as burn-in. The output was visualized in Tracer v1.6 to determine whether convergence and suitable effective sample sizes were achieved for all parameters. The maximum clade credibility tree was accessed using Tree Annotator v1.8.2 [48]. The tree was visualized and edited using Fig Tree v2.1.4.

3. Results

3.1. Population Genetic Structure and Genetic Diversity

The posterior probabilities from the Bayesian analysis and the bootstrap values from the ML analysis were plotted on the ML tree (Figure 2). A. triastriatus fell into two exclusive lineages, A and B, and this result was strongly supported in both the BI tree (1.00) and ML tree (100%). Lineage A consisted of four clades (clades 1–4) and lineage B consisted of two clades (clades 5, 6). Clade 1 was mainly distributed in Guizhou, Hubei, Chongqing and Sichuan provinces (10 of 15 individuals), clades 2 and 3 were mainly distributed in Anhui and Fujian provinces (six of 11 individuals) and clade 4 was mainly distributed in Guangxi, Guangdong and Hunan provinces (six of eight individuals). Clade 5 was mainly distributed in East China, such as Zhejiang, Anhui province, east of Jiangxi and Fujian province (18 of 22 individuals) and clade 6 was mainly distributed in the west of Southern China, such as Sichuan, Hunan, Guangdong, Guangxi and Guizhou provinces (six of nine individuals).
The haplotype network was constructed using a 679 bp fragment of COI from 65 individuals from the 35 locations. It exhibited a certain level of population genetic structure indicating that A. triastriatus had two exclusive lineages, which was in concordance with the phylogenetic trees (Figure 3). The topological relationship of lineage A fell into four clades. Clade 1 included Hap_1-8, clade 2 included Hap_9-11, clade 3 included Hap_12, two individuals from Anhui province and one individual from Zhejiang province of Hap_13, and clade 4 included the remainder of Hap_13 as well as Hap_14. Lineage B fell into two clades. Clade 5 included 22 individuals of Hap_15 and clade 6 included the remainder of Hap_15 as well as Hap_16 and Hap_17.
All 65 individuals of the 35 locations yielded high-quality DNA and were successfully sequenced for the COI gene. Sequences were deposited in GenBank under accession numbers shown in Table S1 of the Supplementary Materials. The mean base composition of the fragment showed a strong bias of A+T (A: 35.03%; C: 18.96%; G: 16.02%; T: 29.99%), as commonly found in invertebrate mitochondrial genomes [54]. The network showed that there were 31 nucleotide changes between lineages A and B, one within lineage B and 1–2 within lineage A. The sequenced region contained 16 polymorphic sites defining 14 haplotypes of lineage A, two polymorphic sites defining three haplotypes of lineage B and 49 polymorphic sites defining 17 haplotypes of A. triastriatus (Table 3). Haplotype and nucleotide diversities were used as measures of genetic diversity between and within lineages. General haplotype diversity was 0.871 for lineage A but only 0.123 for lineage B, and was 0.758 for A. triastriatus. The nucleotide diversity was low in lineage A and extremely low in lineage B. The genetic distances were 0.5% for lineage A, 0.1% for lineage B and 6.3% for A. triastriatus. These results showed that lineages A and B had exclusive haplotypes and suggested a great genetic differentiation exists between lineages.

3.2. Time of Divergence and Colonization History

As shown in Figure 4, the root of the tree, corresponding to the diversification of A. triastriatus and Aporrectodea trapezoides, was estimated to occur about 215.28 Ma (95%: 200.41–229.89 Ma). Amynthas triastriatus became two main lineages about 2.97 Ma (95%: 2.17–3.15 Ma) (Figure 4, Table 4).
Hierarchical population structure in A. triastriatus was studied using AMOVA to check structure at the lineage level based on a fragment of the COI gene (Table 5). The percentage of variation among lineages was 93.89%, and within lineages was just 6.11%. The FST value was 0.9389 (Nm = 0.02), indicating significant population genetic structure. The value of pairwise comparisons of ΦST among lineages was significantly high, suggesting a great genetic differentiation exists between lineages A and B.
Tajima’s D values were negative within lineages, but positive and significant between lineages, indicating that the COI gene was possibly not a neutral evolution. Fu’s Fs test is more sensitive to detect population expansion [46]; Fu’s Fs test values were negative and nonsignificant within lineages, indicating that lineages A and B had experienced demographic expansion. Fu’s Fs test values were positive and nonsignificant between lineages, indicating that A. triastriatus had experienced demographic expansion. The results of nucleotide mismatch distribution analysis based on fragments of COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI genes showed the following: a transition between unimodal and bimodal distributions for COI and 12S, the same as expected for COII, 16S, ND1, ATP8 and ND6 of lineage A; unimodal distribution for COI and ND1, the same as expected for COII, 16S, ATP8 and ND6; and bimodal distribution for the 12S gene of lineage B. The results revealed that lineage A (Figure 5a) and lineage B (Figure 5b) were not steady-state lineages, suggesting an older expansion. BSP revealed an explicit demographic history for the populations of A. triastriatus (Figure 6): a flat population history during previous climate fluctuations of last glacial period (MIS 2-4), a sudden population decline during last glacial period, a rapid population expansion after last glacial period.

4. Discussion

4.1. Application of Molecular Techniques in Taxonomy and Differentiation of Earthworms

Earthworms are important soil animals, but their simple structure and burrowing nature makes their taxonomy challenging and has led to a gross underestimation of the true level of earthworm biodiversity [55]. DNA taxonomy and associated molecular tools might be the only way to reveal the true level of biodiversity [56]. In recent years there has been increasing research use of sequence data along with morphological characters to distinguish earthworm species and to erect cryptic species [6,57,58].
Through analysising the population differentiation and diffusion process, population genetic structure studies and phylogenetic studies could provide an important basis to illustrate the classification statue of species. In previous study, phylogenetic trees and high haplotypes diversity indicated two distinct clades of Drawida japonica Michaelsen, 1892 collected in Shandong and Liaodong peninsulas [17]. A high genetic diversity suggests the presence of five cryptic allopatric species collected from the central Iberian Peninsula [15]. Lineages of A. triastriatus had a strong genetic structure and distinct demographic histories underlying the haplogroups observed in this study. The phylogenetic trees showed that A. triastriatus had two exclusive lineages, a differentiation was also shown by AMOVA, high ΦST index and the high percentage of variation between lineages. Taken together, these results indicated a high degree of isolation and a great genetic differentiation exists between two lineages.
Taxonomically speaking, COI gene divergence of Aporrectodea caliginosa, Lumbricus castaneus, Lumbricus terrestris and Satchellius mammalis was below 4% and ranged between 0.0% and 2.0% within Aporrectodea longa [12,59]. The interspecific COI gene divergence values of invasive Asian earthworms in the northeast United States ranged from 15.84% to 24.03%, and intraspecific COI gene divergence values ranged from 0.01% to 0.4% [29]. In this study, the COI gene divergence between the two lineages of A. triastriatus was 6.3%, which is neither in the intraspecific or interspecific range. However, there were also certain morphological differences between the two lineages. Thus, it is proposed to erect a new subspecies of A. triastriatus (Chen 1946), named Amynthas triastriatus usualis subsp. nov.; the description for this proposed subspecies is included as the final part of this paper.

4.2. Parthenogenesis and the Demographic and Dispersal of A. triastriatus

Parthenogenesis is a unisexual reproduction pattern corresponding to bisexual reproduction that is present in a few animal and plant species. Generally, the loss of sexual reproduction has been thought a dead end in evolution, leading to early extinction [60,61]. But in contrast to this perception, parthenogenesis has been observed in many species of earthworm (Lumbricidae and Megascolecidae) and parthenogenesis in earthworms is often related to polyploidy or aneuploidy [62,63,64,65,66,67,68,69,70]. Muldal (1952) points out that parthenogenesis is important as it makes the retention of polyploidy possible, and also favors the spread of polyploid forms into new areas, since even a single parthenogenetic individual may establish a population. The association between parthenogenesis and high polyploidy in earthworms produces an unexpected level of heterozygosity, an advantageous condition that provides resistance to environmental stress [65,69].
Parthenogenesis often results in polymorphism in earthworms [67], with morphological variability mainly relating to the reduction of reproductive structures such as seminal vesicles, spermathecae, prostates, and an empty seminal chamber. [69]. In the case of A. triastriatus, lineage A with a thin and lustreless seminal chamber and no prostate gland observed was almost degenerated to parthenogenesis, while lineage B with a plump and glossy seminal chamber and small prostate glands had a tendency to parthenogenetic reproduction. Compared to their sexual relatives, parthenogenetics often occur at high latitudes, high altitudes, on islands or island-like habitats, in xeric environments or in disturbed habitats [71,72,73,74,75]. Lineage A was mainly distributed at high altitudes in Southwest China, whereas lineage B was mainly distributed at low altitudes in Southeast China. Unisexual lineages can be highly ecologically successful, with broad geographical ranges that overlap and/or are more extensive than those of the sexual ancestor [76]. For example, the locations of A4 and B5, A8 and B7, A18 and B14 overlapped, and the AMOVA FST values indicated that lineages A and B had a high degree of isolation and a significant genetic differentiation exists with lineage A having greater genetic diversity than lineage B. These findings supported those of Simon et al. (2003), that unisexual lineages with high genetic diversity can live in wider geographical ranges than their sexual ancestor.
Based on phylogenetic studies, parthenogenetic lineages usually occupy terminal nodes of phylogenetic trees; therefore, unisexual forms are usually assumed to arise from sexual congeners [76]. The phylogenetic trees produced from 321 species of Megascolecidae earthworms from China showed that all of these species fell into 14 groups, with A. triastriatus located at the terminal node of the terminal clade of the fourteenth group [52]. In the present study, lineage A was at the terminal node too. Additionally, the time of divergence for lineage A (0.48 Ma; 95%: 0.22–0.53 Ma) was later than for lineage B (0.72 Ma; 95%: 0.29–1.47 Ma), indicating that these unisexual forms did indeed arise from sexual congeners, which was in accordance with previous findings [28,53,77].

4.3. Differentiation and Colonization of A. triastriatus

The divergence time estimate of 215.28 Ma for outgroup Aporrectodea trapezoides (95%: 200.41–229.89 Ma) was not significantly different from previous estimates of the diversification between Lumbricidae and Megascolecidae at about 200–220 Ma [50,51,52,53]. Combining the results of population genetic structure and phylogenetic analysis, as well as paleogeographic events and paleoclimate, it was attempted to determine the differentiation and colonization routes of lineages A and B. In a previous study, the Megascolecidae family of China was separated into 14 groups after pheretima spread into China from the Indo-China Peninsula, then spread to the east or north [52]. Additionally, the split between A. triastriatus and GX201201-11A (an unknown species in the same cluster as A. triastriatus collected in Guangxi province, China) was estimated to have taken place 6.49 Ma (95%: 5.23–7.88 Ma) in Guangxi province [34].
It is proposed that A. triastriatus then formed into two main lineages (A and B) at around 2.97 Ma (95%: 2.17-3.15 Ma), which happened together with the Quaternary glaciation (2.58 Ma). The diversification of lineage A was divided into three stages.
The first stage (node 18 in Figure 4) occurred at the end of Zhonglianggan glaciation (MIS 12) and resulted from the Kunlun-Yellow River tectonic movement (1.2–0.6 Ma BP) that occurred at the transition from the Early to Middle Pleistocene [78,79]. Lineage A (Hap_14) dispersed to Guangdong province (Figure 7). The Heishiding Mountains in Guangdong province was a possible glacial refugium for Liquidambar species during the Quaternary glaciation [80], and the network analysis showed that hap_14 (GD56, collected in Heishiding mountains) was the hap that connected to lineage B. And BSP analysis shown that the Quaternary glaciation have influence on the diffusion of A. triastriatus, significantly several refuge should exist during this period. Therefore, the Heishiding Mountains is proposed as a potential refugium for lineage A.
The second stage (nodes 19, 20, 21 and 25) occurred during MIS 7, 8, 9 and 10 during the Penultimate Glacial Period and before MIS 6, during which glaciation remained widespread with the maximum temperature decrease experienced in Western China and alpine areas [78,81,82,83,84]. Lineage A dispersed north to the mountainous terrain in the Sichuan, Chongqing and Guizhou provinces. Guizhou province, with its comparably stable ecological conditions during environmental fluctuations, is a refugium for Ginkgo biloba [85,86], and lineage A populations from Guizhou province contained 4 haplotypes with high levels of genetic diversity. Therefore, the Sichuan, Chongqing and Guizhou provinces are proposed as potential refugia for lineage A.
The third stage occurred during the last glacial period (MIS 2-4), in which the Gonghe Movement (~0.15 Ma BP) occurred between the Middle and Late Pleistocene to raise the Tibetan Plateau to its present height, High-Asia experienced several temperature plunges, and glaciation existed in most of the mountains in western China [87,88,89,90]. Lineage A dispersed to the east to Fujian province, followed by dispersal to Guangxi province and north to Hunan and Jiangxi provinces. Humid and suitable climate conditions in South China promoted its colonization.
Unlike lineage A, lineage B mainly dispersed to the east. The diversification of lineage B was also divided into three stages. The first stage (nodes 2, 3, 4 and 6 in Figure 4) occurred during the Wangkun Glaciation (MIS 16-20) in which Kunlun-Yellow River tectonic movement not only strongly uplifted the Tibetan Plateau but also the surrounding mountains. Coverage of the Tibetan Plateau by ice sheet peaked during MIS 16 at 18 times the present glacial coverage [91,92]; the mid-Pleistocene Revolution (~0.9 Ma), the most significant transition boundary, led to a significant drop in global temperature and extension of the ice sheet [92]. In this period, lineage B dispersed to Guangdong, Fujian and Zhejiang provinces, far from the Tibetan Plateau. Guangdong (Heishiding), Fujian (Wuyi Mountain) and Zhejiang (Tianmu Mountain) were proposed to be refugia of flora and fauna [86,93]. A lineage B was mainly distributed in Zhejiang province and the surrounding area, it is suggested that this glaciation period blocked the colonization of additional locations by lineage B, resulting in refugia in Zhejiang province.
The second stage (nodes 11 and 13) occurred during MIS 6 and 8, two glacial periods of the Penultimate Glaciation, with lineage B still dispersed in Zhejiang province. The third stage, as for lineage A, occurred during the last glacial period (MIS 2-4). One route of lineage B dispersal was to Anhui and Jiangxi provinces and another route was north to Guizhou and Guangxi provinces. Geological and climatic changes during this period had a profound effect on the colonization of lineage B, and it is also possible that the distant dispersal of B12 was affected by other factors such as birds or human beings.
The population structure of Drawida japonica Michaelsen, 1892 were analysed using 79 samples obtained from the Shandong and Liaodong peninsulas. Results show that the intraspecific genetic diversity of D. japonica has been influenced by geographic isolation [17]. The results of population genetic structure analysis of hormogastrid earthworms (Hormogaster elisae Álvarez 1977) collected from the central Iberian Peninsula and some invasive species (Amynthas corticis (Kinberg, 1867), Amynthas gracilis (Kinberg, 1867)) collected from a volcanic island show that environmental factors may have some influence on earthworms’ genetic evolution [15,94]. The quaternary glaciation resulted in repeated drastic environmental changes that profoundly shaped the current distributions and genetic structure of many plants and animal species in temperate zones of the Northern Hemisphere [95,96,97]. As the youngest major glacial period, Quaternary glaciation ended in a recessional event extending from 0.015 to 0.01 Ma at mid-latitudes [98,99]. The present-day geographical distribution pattern of A. triastriatus formed about 0.01 Ma, with its dispersal process totally covered by the Quaternary glaciation period. This period had a major impact on climate and vegetation in Southeast Asia through the interaction of temperature, rainfall and topography [100,101]. We suggest that the Quaternary glaciation period promoted the differentiation and colonization of A. triastriatus. Drawida japonica, H. elisae, A. corticis, A. gracilis and A. triastriatus are wildspread species, ancient geography event, climate changes, geographic isolation (island, sea, river) and environment factors should play important role in the differentiation and diffusion of wildspread earthworm species.
The phylogenetic evaluation of the genus Amynthas was carried out using 77 Amynthas species from South China. It has been determined that at least one branch of Amynthas spread to the southeast of China, and another spread to the southwest of China [28]. In the present study, the spread of A. triastriatus was consistent with the research of Sun et al., with lineage A widely distributed in the southwest of China, and lineage B widely distributed in the southeast of China. Overall, the colonization routes of A. triastriatus were from south to north, which is in accord with the known dispersal direction of earthworms in China [52,102].

5. Conclusions

Amynthas triastriatus (Oligochaete: Megascolecidae) of the aeruginosus-group is a widely distributed endemic earthworm species in Southern China. Population genetic structure exhibited a very high degree of isolation and two exclusive lineages based on existing specimen and gene data. Lineage A was almost degenerated to parthenogenesis, lineage B had a tendency to parthenogenesis. There were distinguish genetic differentiations exist between lineages and lineage A had high haplotype diversity and strong genetic structure, additionally, lineage A mainly distributied at high alitude areas. These results indicated a new cryptic species of A. triastriatus; combined with morphological differences, a new subspecies of A. triastriatus is proposed. And the Quaternary glaciation may be one of main factors that promoted the colonization of A. triastriatus. In future, there may be more studies of differentiation and dispersal of A. triastriatus including more populations based on this study, and the population genetic structure analysis should play an important role in the earthworm taxonomy.

6. Description of a New Subspecies of Amynthas Triastriatus

Based on the above results, to unravel the classification puzzle of Amynthas triastriatus it is proposed that a new subspecies be established that encompasses the populations of lineage B.
Material. Holotype: 1 clitellate (C-FJ201111-04A): China, Fujian Province, Meihua Mountain Nature Reserve (25.31278° N, 116.89306° E), 1207 m asl, black soil beside road, 13th July, 2011, J. B. Jiang, J. Sun, X. D. Lei and H. W. Feng coll. Paratypes: 85 clitellates: 2 clitellates (C-FJ201111-04B) same date as for holotype. 4 clitellates (C-FJ201004-01): China, Guangxi Province, Youxi Country (26.11244° N, 118.42387° E), 952 m asl, black soil beside road, 24th Aug, 2010, J. B. Jiang coll. 1 clitellate (C-GD201115-04): China, Guangdong Province, Lianping county (23.46417° N, 114.91444° E), 410 m asl, brown soil under eucalyptus trees, 04th Aug, 2011, J. B. Jiang coll. 12 clitellates (C-GX201319-06): China, Guangxi Province, Jiuwan Mountain Nature Reserve (25.20595° N, 108.67555° E), 1226 m asl, yellow soil, 18th May, 2013, J. P. Qiu, Y, Hong, J. B. Jiang, L. L. Zhang and Y. Dong coll. 4 clitellates (C-GZ201609-02): China, Guizhou Province, Fanjing Mountain National Reserve (27.52480° N, 108.59300° E), 1316 m asl, yellow soil under bryophyte beside road, 21th May, 2013, J. P. Qiu, J. B. Jiang, L. L. Zhang and Y. Dong coll. 2 clitellates (C-HU201309-02): China, Hunan Province, Bamian Mountain Nature Reserve (27.25680° N, 112.72433° E), 1034 m asl, black soil under willow trees, 21th Jun, 2013, J. P. Qiu, J. B. Jiang, L. L. Zhang and Y. Dong coll. 2 clitellates (C-JX201302-04): China, Jiangxi Province, Longhu mountain (28.03037° N, 117.09680° E), 185 m asl, black soil beside road, 15th Jun, 2013, J. B. Jiang, L. L. Zhang, Y. Dong and S. Lei coll. 19 clitellates (C-JX201311-01): China, Jiangxi Province, Jinggang Mountain National Reserve (26.55213° N, 114.10773° E), 1300 m asl, black soil under bamboo forest beside road, 17th Jun, 2013, J. B. Jiang, L. L. Zhang, Y. Dong and S. Lei coll. 4 clitellates (C-SC201004-06): China, Sichuan Province, Gongga Mountain National Reserve (29.65556° N, 102.11667° E), 1558 m asl, brown soil beside river, 06th Aug, 2010, J. P. Qiu, J. B. Jiang coll. 1 clitellate (C-ZJ201101-02): China, Zhejiang Province, Tianmu Mountain National Reserve (30.32886° N, 102.11667° E), 1052 m asl, black soil beside road, 02nd Nov, 2011, J. B. Jiang, J. Sun, L. L. Zhang coll. 1 clitellate (C-ZJ201511-01): China, Zhejiang Province, Jinhua City (28.83361° N, 119.24639° E), 160 m asl, black soil beside road, 08th July, 2015, Y. Dong, M. S. Chen, X. Gao and C. Peng coll. 4 clitellates (C-ZJ201515-03): China, Zhejiang Province, Yunhe City (28.10444° N, 119.63472° E), 290 m asl, yellow soil beside road, 09th July, 2015, Y. Dong, M. S. Chen, X. Gao and C. Peng coll. 1 clitellate (C-ZJ201517-02): China, Zhejiang Province, Wenzhou City (27.84556° N, 119.77306° E), 710 m asl, brown soil beside road, 09th July, 2015, Y. Dong, M. S. Chen, X. Gao and C. Peng coll. 6 clitellates (C-ZJ201518-01): China, Zhejiang Province, Yandang mountain (28.37639° N, 121.10972° E), 120 m asl, brown soil beside road, 10th July, 2015, Y. Dong, M. S. Chen, X. Gao and C. Peng coll. 10 clitellates (C-AH201601-01): China, Anhui Province, Huangshan City (29.95601° N, 118.09112° E), 325 m asl, brown soil, 08th May, 2016, J. B. Jiang, J. Sun, Y. Dong and Y. Zheng coll. 2 clitellates (C-AH201609-02): China, Anhui Province, Huangshan City (30.31276° N, 118.02561° E), 120 m asl, yellow soil, 09th May, 2016, J. B. Jiang, J. Sun, Y. Dong and Y. Zheng coll. 10 clitellates (C-AH201610-01): China, Anhui Province, Huangshan City (30.35924° N, 117.93609° E), 146 m asl, brown soil, 09th May, 2016, J. B. Jiang, J. Sun, Y. Dong and Y. Zheng coll.
Diagnosis. Dimensions 120–150 mm by 4.9–6.8 mm at clitellum, segments 108–111. First dorsal pore in 11/12. Setae numbering 20–26 at III, 24–30 at V, 30–34 at VIII, 40–48 at XX and 55–60 at XXV; 12–16 between male pores; seta numbering between spermathecal pores 10–14 at VIII. Two pairs of spermathecal pores in 7/8–8/9, 1/3 circumference ventrally apart from each other. One pair of male pores in XVIII, 1/3 circumference apart ventrally, each on the top of a raised, round porophore surrounded by 3–5 circular ridges. Two presetal and two postsetal collapse-topped genital papillae present on XVIII. Ampulla oval-shaped, stout duct as long as 1/5 ampulla. Diverticulum as long as 1/3 main pouch, terminal 1/2 dilated into an oval-shaped glossy seminal chamber. Prostate glands small or only one side observed.
Description. External characters: Brown pigment on dorsum, no pigment on ventrum. Dimensions 125 mm by 5.1 mm at clitellum, segments 111. Prostomium 1/2 epilobous. First dorsal pore in 11/12. Setae numbering 22 at III, 26 at V, 32 at VIII, 44 at XX and 58 at XXV; 14 between male pores; 14 at VIII between spermathecal pores. Setal formula: AA = 1.2 − 1.4AB, ZZ = 1.6 − 2.0ZY. Clitellum annular, taupe on dorsal, in XIV–XVI, setae not visible external. Two pairs of spermathecal pores in 7/8–8/9, 1/3 circumference apart ventrally. Two pairs of postsetal genital papillae present on VII, two pairs presetal and two pairs postsetal on VIII, three pairs presetal on IX, 1/3 circumference ventrally apart from each other (paratypes: numbers of genital papillae variable from 14–18). One pair of male pores in XVIII, 1/3 circumference apart ventrally, each on the top of a raised, round porophore surrounded by 3–5 circular ridges, two presetal and two postsetal genital papillae on median side. Two big collapse-topped genital papillae present on XVIII, 1/4 circumference ventrally apart from each other. (Figure 8A). Single female pore in XIV, milky.
Internal characters. Septa 5/6–7/8, thick and muscular, 10/11–12/13 slightly thickened, 8/9–9/10 absent. Gizzard bucket-shaped, in IX–X. Intestine enlarged distinctly from XV. One pair of intestinal caeca in XXVII, brown, simple, dorsal margin smooth, 17 short pointed saccules in ventral margin, extending anteriorly to XXIII. Four pairs of esophageal hearts in X–XIII. Ovaries in XIII. Two pairs of spermathecae in VIII–IX, about 1.7–2.0 mm long. Ampulla oval-shaped, stout duct as long as 1/5 ampulla. Diverticulum as long as 1/3 main pouch, terminal 1/2 dilated into an oval-shaped glossy seminal chamber. (Figure 8B). Holandric: two pairs of testis sacs, in X–XI undeveloped. Two pairs of seminal vesicles in XI and XII, the second pair more developed than the first pair. Holotype with one pair of thick, small, nubby lobate prostate glands and S-curved duct. Paratypes with one pair of prostate glands or only one side observed and other side S-curved duct (Figure 8C).
Ethymology. The species is named after the existence of prostate glands.
Remarks. There are certain differences between the new subspecies and A. triastriatus. For instance, the new subspecies has small prostate glands that are absent in A. triastriatus, and the seminal chamber of the new subspecies is plump and glossy while thin and lustreless in A. triastriatus. Additionally, the first dorsal pore of the new subspecies is located at 11/12, but at 10/11 in A. triastriatus. Furthermore, the new subspecies has more papillae in both the spermathecal and male pore regions than A. triastriatus, and more finger sacs in the ventral intestinal caeca than A. triastriatus.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/17/5/1538/s1, Sequences of A. triastriatus were deposited in GenBank under accession numbers shown in Table S1.

Author Contributions

J.Q. initiated the project. Y.D. and J.Q. designed the study. Y.D. drafted the manuscript. Y.D., J.J., Z.Y. and Q.Z. collected the data. Y.D. and J.Q. performed molecular analyses and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was founded by National Nature Science Foundation of China grant No. 41701272, No. 41771279, No. 41471204, No. 31272265, No. 31401967, No. 31172360, and a National Science & Technology Fundamental Resources Investigation Program of China 2018FY100300.

Acknowledgments

We would like to thank Prof. Hong Yong from the College of Agriculture, Chonbuk National University; Jing Sun from Sichuan Institute of Industrial Technology; Mingsheng Chen from Shanghai Jiao Tong University; Yali Wang from Henan University of Technology; Liangliang Zhang and Zhimin Tao for their valuable help with the collections, Weibin Jiang from Shanghai Normal University for software support, and Graeme Brown for English language editing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hoffman, E.A.; Blouin, M.S. Evolutionary history of the northern leopard frog: Reconstruction of phylogeny, phylogeography, and historical changes in population demography from mitochondrial DNA. Evolution 2004, 58, 145–159. [Google Scholar] [CrossRef]
  2. Nyakaana, S.; Arctander, P. Population genetic structure of the African elephant in Uganda based on variation at mitochondrial and nuclear loci: Evidence for male-biased gene flow. Mol. Ecol. 1999, 8, 1105–1115. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, H.L.; Chen, H.; Hu, L.L.; Wang, L.X.; Wang, S.H.; Wang, M.L.; Cheng, X.Z. Genetic diversity a population structure analysis of a accessions in the Chinese cowpea [Vigna unguiculata (L.) Walp.] germplasm collection. Crop J. 2017, 5, 363–372. [Google Scholar] [CrossRef]
  4. Winding, A.; Ronn, R.; Hendriksen, N.B. Bacteria and protozoa in soil microhabitats as affected by earthworms. Biol. Fert. Soils 1997, 24, 133–140. [Google Scholar] [CrossRef]
  5. Fernández, R.; Almodóvar, A.; Novo, M.; Gutiérrez, M.; Díaz Cosín, D.J. A vagrant clone in a peregrine species: Phylogeography, high clonal diversity and geographical distribution in the earthworm Aporrectodea trapezoides (Dugès, 1828). Soil Biol. Biochem. 2011, 43, 2085–2093. [Google Scholar] [CrossRef]
  6. Novo, M.; Almodovar, A.; Diaz-Cosin, D.J. High genetic divergence of hormogastrid earthworms (Annelida, Oligochaeta) in the central Iberian Peninsula: Evolutionary and demographic implications. Zool. Scr. 2009, 38, 537–552. [Google Scholar] [CrossRef]
  7. Shekhovtsov, S.V.; Berman, D.I.; Peltek, S.E. Phylogeography of the earthworm Eisenia nordenskioldi nordenskioldi (Lumbricidae, Oligochaeta) in northeastern Eurasia. Dokl. Biol. Sci. 2015, 461, 85–88. [Google Scholar] [CrossRef]
  8. Fernández, R.; Almodóvar, A.; Novo, M.; Gutiérrez, M.; Díaz Cosín, D.J. Earthworms, good indicators for palaeogeographical studies? Testing the genetic structure and demographic history in the peregrine earthworm Aporrectodea trapezoides (Dugès, 1828) in southern Europe. Soil Biol. Biochem. 2013, 58, 127–135. [Google Scholar] [CrossRef]
  9. Thakuria, D.; Schmidt, O.; Finan, D.; Egan, D.; Doohan, F.M. Gut wall bacteria of earthworms: A natural selection process. ISME J. 2010, 4, 357–366. [Google Scholar] [CrossRef] [Green Version]
  10. Groffman, P.M.; Fahey, T.J.; Fisk, M.C.; Yavitt, J.B.; Sherman, R.E.; Bohlen, P.J.; Maerz, J.C. Earthworms increase soil microbial biomass carrying capacity and nitrogen retention in northern hardwood forests. Soil Biol. Biochem. 2015, 87, 51–58. [Google Scholar] [CrossRef]
  11. Braga, L.P.P.; Yoshiura, C.A.; Borges, C.D.; Horn, M.A.; Brown, G.G.; Drake, H.L.; Tsai, S.M. Disentangling the influence of earthworms in sugarcane rhizosphere. Sci. Rep. 2016, 6, 38923. [Google Scholar] [CrossRef] [PubMed]
  12. King, R.A.; Tibble, A.L.; Symondson, W.O.C. Opening a can of worms: Unprecedented sympatric cryptic diversity within British lumbricid earthworms. Mol. Ecol. 2008, 17, 4684–4698. [Google Scholar] [CrossRef] [PubMed]
  13. Perez-Losada, M.; Eiroa, J.; Mato, S.; Dominguez, J. Phylogenetic species delimitation of the earthworms Eisenia fetida (Savigny, 1826) and Eisenia andrei Bouche, 1972 (Oligochaeta, Lumbricidae) based on mitochondrial and nuclear DNA sequences. Pedobiologia 2005, 49, 317–324. [Google Scholar] [CrossRef]
  14. Domínguez, J.; Bohlen, P.J.; Parmelee, R.W. Earthworms Increase Nitrogen Leaching to Greater Soil Depths in Row Crop Agroecosystems. Ecosystems 2004, 7, 672–685. [Google Scholar] [CrossRef]
  15. Novo, M.; Almodovar, A.; Fernandez, R.; Trigo, D.; Cosin, D.J.D. Cryptic speciation of hormogastrid earthworms revealed by mitochondrial and nuclear data. Mol. Phylogenet. Evol. 2010, 56, 507–512. [Google Scholar] [CrossRef] [PubMed]
  16. Minamiya, Y.; Yokoyama, J.; Fukuda, T. A phylogeographic study of the Japanese earthworm, Metaphire sieboldi (Horst, 1883) (Oligochaeta: Megascolecidae): Inferences from mitochondrial DNA sequences. Eur. J. Soil Biol. 2009, 45, 423–430. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, Y.F.; Zhang, D.H.; Xu, Y.L.; Zhang, G.S.; Sun, Z.J. Effects of fragmentation on genetic variation in populations of the terrestrial earthworm Drawida japonica Michaelsen, 1892 (Oligochaeta, Moniligastridae) in Shandong and Liaodong peninsulas, China. J. Nat. Hist. 2012, 46, 1387–1405. [Google Scholar] [CrossRef]
  18. Aspe, N.M.; James, S.W. Molecular phylogeny and biogeographic distribution of pheretimoid earthworms (clitellata: Megascolecidae) of the Philippine archipelago. Eur. J. Soil Biol. 2018, 85, 89–97. [Google Scholar] [CrossRef]
  19. Dupont, L.; Gresille, Y.; Richard, B.; Decaens, T.; Mathieu, J. Dispersal constraints and fine-scale spatial genetic structure in two earthworm species. Biol. J. Linn. Soc. 2015, 114, 335–347. [Google Scholar] [CrossRef]
  20. Novo, M.; Almodovar, A.; Fernandez, R.; Giribet, G.; Cosin, D.J.D. Understanding the biogeography of a group of earthworms in the Mediterranean basin-The phylogenetic puzzle of Hormogastridae (Clitellata: Oligochaeta). Mol. Phylogenet. Evol. 2011, 61, 125–135. [Google Scholar] [CrossRef]
  21. James, S.W.; Porco, D.; Decaens, T.; Richard, B.; Rougerie, R.; Erseus, C. DNA Barcoding Reveals Cryptic Diversity in Lumbricus terrestris L., 1758 (Clitellata): Resurrection of L. herculeus (Savigny, 1826). PLoS ONE 2010, 5, e15629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Buckley, T.R.; James, S.; Allwood, J.; Bartlam, S.; Howitt, R.; Prada, D. Phylogenetic analysis of New Zealand earthworms (Oligochaeta: Megascolecidae) reveals ancient clades and cryptic taxonomic diversity. Mol. Phylogenet. Evol. 2011, 58, 85–96. [Google Scholar] [CrossRef] [PubMed]
  23. Richard, B.; Decaens, T.; Rougerie, R.; James, S.W.; Porco, D.; Hebert, P.D.N. Re-integrating earthworm juveniles into soil biodiversity studies: Species identification through DNA barcoding. Mol. Ecol. Resour. 2010, 10, 606–614. [Google Scholar] [CrossRef] [PubMed]
  24. Chang, C.H.; James, S. A critique of earthworm molecular phylogenetics. Pedobiologia 2011, 54, S3–S9. [Google Scholar] [CrossRef]
  25. Csuzdi, C.; Chang, C.H.; Pavlicek, T.; Szederjesi, T.; Esopi, D.; Szlavecz, K. Molecular phylogeny and systematics of native North American lumbricid earthworms (Clitellata: Megadrili). PloS ONE 2017, 12, 8–e0181504. [Google Scholar] [CrossRef] [Green Version]
  26. James, S.W.; Davidson, S.K. Molecular phylogeny of earthworms (Annelida:Crassiclitellata) based on 28S, 18S and 16S gene sequences. Invertebr. Syst. 2012, 26, 213–229. [Google Scholar] [CrossRef]
  27. Schubert, M.; Csuzdi, C.; Chang, C.-H.; Pavlícek, T.; Szederjesi, T.; Esopi, D.; Szlávecz, K. Molecular phylogeny and systematics of native North American lumbricid earthworms (Clitellata: Megadrili). PLoS ONE 2017, 12, e0181504. [Google Scholar]
  28. Sun, J.; James, S.W.; Jiang, J.B.; Yao, B.; Zhang, L.L.; Liu, M.Q.; Qiu, J.P.; Hu, F. Phylogenetic evaluation of Amynthas earthworms from South China reveals the initial ancestral state of spermathecae. Mol. Phylogenet. Evol. 2017, 115, 106–114. [Google Scholar] [CrossRef]
  29. Schult, N.; Pittenger, K.; Davalos, S.; McHugh, D. Phylogeographic analysis of invasive Asian earthworms (Amynthas) in the northeast United States. Invertebr. Biol. 2016, 135, 314–327. [Google Scholar] [CrossRef]
  30. Zhang, L.L.; Jiang, J.B.; Dong, Y.; Qiu, J.P. Complete mitochondrial genome of four pheretimoid earthworms (Clitellata: Oligochaeta) and their phylogenetic reconstruction. Gene 2015, 574, 308–316. [Google Scholar] [CrossRef]
  31. Halanych, K.M.; Janosik, A.M. A review of molecular markers used for annelid phylogenetics. Integr. Comp. Biol. 2006, 46, 533–543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Nygren, A.; Sundberg, P. Phylogeny and evolution of reproductive modes in Autolytinae (Syllidae, Annelida). Mol. Phylogenet. Evol. 2003, 29, 235–249. [Google Scholar] [CrossRef]
  33. Perez-Losada, M.; Ricoy, M.; Marshall, J.C.; Dominguez, J. Phylogenetic assessment of the earthworm Aporrectodea caliginosa species complex (Oligochaeta: Lumbricidae) based on mitochondrial and nuclear DNA sequences. Mol. Phylogenet. Evol. 2009, 52, 293–302. [Google Scholar] [CrossRef] [PubMed]
  34. Jiang, J.B. Taxonomy and Molecular Phylogeny of the Family Megascolecidae Earthworms from China. Ph.D. Thesis, Shanghai Jiao Tong University Shanghai, Shanghai, China, 2016. [Google Scholar]
  35. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  36. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar] [PubMed]
  37. Bely, A.E.; Wray, G.A. Molecular phylogeny of naidid worms (Annelida: Clitellata) based on cytochrome oxidase I. Mol. Phylogenet. Evol. 2004, 30, 50–63. [Google Scholar] [CrossRef]
  38. Hillis, D.M.; Mable, B.K.; Larson, A.; Davis, S.K.; Zimmer, E.A. Nucleic acids, IV: Sequencing and cloning. In Molecular Systematics; Hillis, D.M., Moritz, C., Mable, B.K., Eds.; Sinauer: Sutherland, MA, USA, 1996; pp. 321–381. [Google Scholar]
  39. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef] [Green Version]
  40. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [Green Version]
  41. Guindon, S.; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [Green Version]
  42. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Hohna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  43. Excoffier, L.; Laval, G.; Schneider, S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol. Bioinform. Online 2005, 1, 47–50. [Google Scholar] [CrossRef] [Green Version]
  44. Bandelt, H.J.; Forster, P.; Rohl, A. Median-Joining Networks for Inferring Intraspecific Phylogenies. Mol. Biol. Evol. 1999, 16, 37–48. [Google Scholar] [CrossRef] [PubMed]
  45. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphisms. Genetics 1989, 123, 585–595. [Google Scholar] [PubMed]
  46. Fu, Y.X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 1997, 147, 915–925. [Google Scholar] [PubMed]
  47. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  48. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [Green Version]
  49. Rambaut, A.; Drummond, A.J.; Dong, X.; Baele, G.; Suchard, M.A. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef] [Green Version]
  50. Yuan, Z.; Jiang, J.B.; Dong, Y.; Zhao, Q.; Gao, X.; Qiu, J.P. The dispersal and diversification of earthworms (Annelida: Oligochaeta) related to paleogeographical events in the Hengduan Mountains. Eur. J. Soil Biol. 2019, 94, 103–118. [Google Scholar] [CrossRef]
  51. Bouche, M. L’eolution spatiotemporelle des lombrixiens. Doc. Pedozool. Integrol. 1998, 3, 1–28. [Google Scholar]
  52. Jiang, J.B.; Qiu, J.P. Origin and evolution of earthworms belonging to the family Megascolecidae in China. Biodivers. Sci. 2018, 26, 1074–1082. [Google Scholar] [CrossRef] [Green Version]
  53. Dominguez, J.; Aira, M.; Breinholt, J.W.; Stojanovic, M.; James, S.W.; Perez-Losada, M. Underground evolution: New roots for the old tree of lumbricid earthworms. Mol. Phylogenet. Evol. 2015, 83, 7–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Boore, J.L.; Brown, W.M. Complete DNA sequence of the mitochondrial genome of the annelid worm Lumbricus terrestris. Genet. Mol. Res. 1995, 141, 305–319. [Google Scholar]
  55. Sket, B. The nature of biodiversity in hypogean waters and how it is endangered. Biodivers. Conserv. 1999, 8, 1319–1338. [Google Scholar] [CrossRef]
  56. Proudlove, G.; Wood, P.J. The blind leading the blind: Cryptic subterranean species and DNA taxonomy. Trends Ecol. Evol. 2003, 18, 272–273. [Google Scholar] [CrossRef]
  57. Klarica, J.; Kloss-Brandstatter, A.; Traugott, M.; Juen, A. Comparing four mitochondrial genes in earthworms—Implications for identification, phylogenetics, and discovery of cryptic species. Soil Biol. Biochem. 2012, 45, 23–30. [Google Scholar] [CrossRef]
  58. Decaens, T.; Porco, D.; Rougerie, R.; Brown, G.G.; James, S.W. Potential of DNA barcoding for earthworm research in taxonomy and ecology. Appl. Soil Ecol. 2013, 65, 35–42. [Google Scholar] [CrossRef]
  59. Martinsson, S.; Rhodén, C.; Erséus, C. Barcoding gap, but no support for cryptic speciation in the earthworm Aporrectodea longa (Clitellata: Lumbricidae). Mitochondrial DNA A DNA Mapp. Seq. Anal. 2017, 28, 147–155. [Google Scholar]
  60. Welch, D.M.; Meselson, M. Evidence for the evolution of bdelloid rotifers without sexual reproduction or genetic exchange. Science 2000, 288, 1211–1215. [Google Scholar] [CrossRef] [Green Version]
  61. Graham, B. The Masterpiece of Nature: The Evolution and Genetics of Sexuality; University of Califormia Press: Berkeley, CA, USA, 1982; p. 635. [Google Scholar]
  62. Muldal, S. The chromosomes of the earthworms: I. The evolution of polyploidy. Helv. Chim. Acta 1952, 6, 56–76. [Google Scholar] [CrossRef] [Green Version]
  63. Casellato, S. On polyploidy in oligochaetes with particular reference to Lumbricids. In On Earthworms; Omodeo, P., Ed.; Bonvicini Pagliai AM: Modena, Italy, 1987; pp. 77–87. [Google Scholar]
  64. Viktorov, A.G. Diversity of polyploid races in the family Lumbricidae. Soil Biol. Biochem. 1997, 29, 217–221. [Google Scholar] [CrossRef]
  65. Edwards, C.; Bohlen, P. Biology and Ecology of Earthworms, 3rd ed.; Chapman and Hall: London, UK, 1996; p. 426. [Google Scholar]
  66. Omodeo, P. Specie e razze poliploidi nei lombrichi. Conv. Genet. 1953, 23, 43–49. [Google Scholar]
  67. Gates, G. Contributions to North American earthworms. No. 3. IV. The trapezoides species group. Bull. Tall Timbers Res. Stn. 1972, 12, 146. [Google Scholar]
  68. Qiu, J.; Bouché, M. Contribution to the taxonomy of Hormogastridae (Annelida: Oligochaeta) with description of new species from Spain. Doc. Pedozool. Integrol. 1998, 4, 164–177. [Google Scholar]
  69. Cosín, D.J.D.; Novo, M.; Fernández, R. Reproduction of Earthworms: Sexual Selection and Parthenogenesis. In Biology of Earthworms; Springer: Berlin/Heidelberg, Germany, 2011; pp. 69–86. [Google Scholar]
  70. Terhivuo, J.; Saura, A. Dispersal and clonal diversity of North-European parthenogenetic earthworms. Biol. Invasions 2006, 8, 1205–1218. [Google Scholar] [CrossRef]
  71. Ghiselin, M.T. The Economy of Nature and the Evolution of Sex; University of California Press: Berkeley, CA, USA, 1974; p. 346. [Google Scholar]
  72. Glesener, R.R.; Tilman, D. Sexuality and the Components of Environmental Uncertainty: Clues from Geographic Parthenogenesis in Terrestrial Animals. Am. Nat. 1978, 112, 659–673. [Google Scholar] [CrossRef]
  73. Hamilton, W.D.; Henderson, P.A.; Moran, N.; Alexander, R.; Tinkle, D.; Press, C. Fluctuation of environment and coevolved antagonist polymorphism as factors in the maintenance of sex. In Natural Selection and Social Behaviour; Alexander, R.D., Tinkle, D.W., Eds.; Chiron: New York, NY, USA, 1981; pp. 363–381. [Google Scholar]
  74. Lynch, M. Destabilizing Hybridization, General-Purpose Genotypes and Geographic Parthenogenesis. Q. Rev. Biol. 1984, 59, 527–590. [Google Scholar] [CrossRef]
  75. Kearney, M.; Blacket, M.J.; Strasburg, J.L.; Moritz, C. Waves of parthenogenesis in the desert: Evidence for the parallel loss of sex in a grasshopper and a gecko from Australia. Mol. Ecol. 2006, 15, 1743–1748. [Google Scholar] [CrossRef] [PubMed]
  76. Simon, J.C.; Delmotte, F.; Rispe, C.; Crease, T. Phylogenetic relationships between parthenogens and their sexual relatives: The possible routes to parthenogenesis in animals. Biol. J. Linn. Soc. 2003, 79, 151–163. [Google Scholar] [CrossRef]
  77. Minamiya, Y.; Hayakawa, H.; Ohga, K.; Shimano, S.; Ito, M.T.; Fukuda, T. Variability of sexual organ possession rates and phylogenetic analyses of a parthenogenetic Japanese earthworm, Amynthas vittatus (Oligochaeta: Megascolecidae). Genes Genet. Syst. 2011, 86, 27–35. [Google Scholar] [CrossRef] [Green Version]
  78. Zhou, S.Z.; Li, J.J.; Zhang, S.Q. Quaternary glaciation of the Bailang River Valley, Qilian Shan. Quat. Int. 2002, 97–98, 103–110. [Google Scholar] [CrossRef]
  79. Jiang, Z.Y.; Wu, Y.Q.; Cui, Z.J. Kunlun-Yellow River Tectonic Motion and formation of modern physical geography pattern of China. J. Beijing Normal Univ. 2005, 41, 85–88. [Google Scholar]
  80. Wu, W. Natural Hybridization, PhyIogeography and Speciation Patterns of Altingiaceae. Ph.D. Thesis, Sun Yat-sen University Guangdong, Guangdong, China, 2009. [Google Scholar]
  81. Jörg, M.S.; Silvio, T.; Zhizhong, Z.; Xihao, W.; Susan, I.-O.; Rainer, W.; Heinrich, B.; Peter, W.K.; Christian, S. The limited influence of glaciations in Tibet on global climate over the past 170,000 yr. Earth Planet. Sci. Lett. 2002, 194, 287–297. [Google Scholar]
  82. Zhao, J.; Zhou, S.Z.; He, Y.Q.; Ye, Y.G.; Liu, S.Y. ESR dating of glacial tills and glaciations in the Urumqi River headwaters, Tianshan Mountains, China. Quat. Int. 2006, 144, 61–67. [Google Scholar] [CrossRef]
  83. Owen, L.A.; Finkel, R.C.; Barnard, P.L.; Haizhou, M.; Asahi, K.; Caffee, M.W.; Derbyshire, E. Climatic and topographic controls on the style and timing of Late Quaternary glaciation throughout Tibet and the Himalaya defined by 10Be cosmogenic radionuclide surface exposure dating. Quat. Sci. Rev. 2005, 24, 1391–1411. [Google Scholar]
  84. Zhao, J.; Liu, S.; He, Y.; Song, Y. Quaternary glacial chronology of the Ateaoyinake River Valley, Tianshan Mountains, China. Geomorphology 2009, 103, 276–284. [Google Scholar] [CrossRef]
  85. Shen, L.; Chen, X.Y.; Zhang, X.; Li, Y.Y.; Fu, C.X.; Qiu, Y.X. Genetic variation of ginkgo biloba l. (ginkgoaceae) based on cpdna pcr-rflps: Inference of glacial refugia. Heredity 2004, 94, 396–401. [Google Scholar] [CrossRef] [Green Version]
  86. Gong, W.; Chen, C.; Dobeš, C.; Fu, C.-X.; Koch, M. Phylogeography of a living fossil: Pleistocene glaciations forced Ginkgo biloba L. (Ginkgoaceae) into two refuge areas in China with limited subsequent postglacial expansion. Mol. Phylogenet. Evol. 2008, 48, 1094–1105. [Google Scholar] [CrossRef]
  87. Li, J.J.; Fang, X.M.; Pan, B.T.; Zhao, Z.J.; Song, Y.G. Late cenozoic intensive uplift of qinghai-xizang plateau and its impacts on environments in surrounding area. Quat. Sci. 2001, 21, 381–391. [Google Scholar]
  88. Thompson, L.G.; Yao, T.; Davis, M.E.; Henderson, K.A.; MosleyThompson, E.; Lin, P.N.; Beer, J.; Synal, H.A.; ColeDai, J.; Bolzan, J.F. Tropical climate instability: The last glacial cycle from a Qinghai-Tibetan ice core. Science 1997, 276, 1821–1825. [Google Scholar] [CrossRef] [Green Version]
  89. Yi, C.L.; Zhu, Z.Y.; Wei, L.; Cui, Z.J.; Zheng, B.X.; Shi, Y.F. Advances in numerical dating of Quaternary glaciations in China. Z. Geomorphol. 2007, 51, 153–175. [Google Scholar] [CrossRef]
  90. Yang, J.Q.; Zhang, W.; Cui, Z.J.; Yi, C.L.; Chen, Y.X.; Xu, X.G. Climate change since 11.5 ka on the Diancang Massif on the southeastern margin of the Tibetan Plateau. Quat. Res. 2010, 73, 304–312. [Google Scholar] [CrossRef]
  91. Cui, Z.J.; Wu, Y.Q.; Liu, G.N.; Ge, A.K.; Pang, Q.Q.; Xu, Q.H. On Kunlun-Yellow River tectonic movement. Sci. China 1998, 41, 592–600. [Google Scholar] [CrossRef]
  92. Shi, Y.F. Evolution of the cryosphere in the Tibetan Plateau, China, and its relationship with the global change in the mid-Quaternary. J. Glaciol. Geocryol. 1998, 20, 197–208. [Google Scholar]
  93. Chen, C.D. Biodiversity in the Wuyi Mountains and its importance in China. Chin. Biodivers. 1999, 7, 320–326. [Google Scholar]
  94. Novo, M.; Cunha, L.; Maceda-Veiga, A.; Talavera, J.A.; Hodson, M.E.; Spurgeon, D.; Bruford, M.W.; Morgan, A.J.; Kille, P. Multiple introductions and environmental factors affecting the establishment of invasive species on a volcanic island. Soil Biol. Biochem. 2015, 85, 89–100. [Google Scholar] [CrossRef] [Green Version]
  95. Hewitt, G.M. Genetic consequences of climatic oscillations in the Quaternary. Philos. Tans. R. Soc. Lond. B Biol. Sci. 2004, 359, 183–195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Qiu, Y.X.; Fu, C.X.; Comes, H.P. Plant molecular phylogeography in China and adjacent regions: Tracing the genetic imprints of Quaternary climate and environmental change in the world’s most diverse temperate flora. Mol. Phylogenet. Evol. 2011, 59, 225–244. [Google Scholar] [CrossRef]
  97. Ren, G.P.; Mateo, R.G.; Liu, J.Q.; Suchan, T.; Alvarez, N.; Guisan, A.; Conti, E.; Salamin, N. Genetic consequences of Quaternary climatic oscillations in the Himalayas: Primula tibetica as a case study based on restriction site-associated DNA sequencing. New Phytol. 2017, 213, 1500–1512. [Google Scholar] [CrossRef] [Green Version]
  98. Wright, H.E. Late Pleistocene climate of Europe: A review. GSA Bull. 1961, 72, 933–984. [Google Scholar] [CrossRef]
  99. Dawson, A.G. Ice Age Earth; Routledge: London, UK, 1992; p. 293. [Google Scholar]
  100. Flenley, J.R. Quaternary Vegetation. (Book Reviews: The Equatorial Rain Forest. A Geological History). Science 1979, 206, 550. [Google Scholar]
  101. Verstappen, H. Quaternary Climatic Changes and Natural Environment in SE Asia. GeoJournal 1980, 4, 45–54. [Google Scholar] [CrossRef]
  102. Sun, J. Taxonomy and Molecular Phylogeny of Amynthas Earthworms from China. Ph.D. Thesis, Shanghai Jiao Tong University Shanghai, Shanghai, China, 2013. [Google Scholar]
Ijerph 17 01538 g001 550
Figure 1. The 35 sampling locations in China. The circles represent 18 locations of lineage A. The crosses represent 17 locations of lineage B. Location details are shown in Table 1.
Figure 1. The 35 sampling locations in China. The circles represent 18 locations of lineage A. The crosses represent 17 locations of lineage B. Location details are shown in Table 1.
Ijerph 17 01538 g001
Ijerph 17 01538 g002 550
Figure 2. Phylogenetic trees reconstructed using different methods (ML and BI) based on the COI gene. Branch lengths are proportional to the amount of change along the branches. Bootstrap proportions (if ≥70%) and Bayesian posterior probabilities (if ≥0.95) are shown above and below the branches, respectively. Aporrectodea trapezoides was used as the outgroup taxa.
Figure 2. Phylogenetic trees reconstructed using different methods (ML and BI) based on the COI gene. Branch lengths are proportional to the amount of change along the branches. Bootstrap proportions (if ≥70%) and Bayesian posterior probabilities (if ≥0.95) are shown above and below the branches, respectively. Aporrectodea trapezoides was used as the outgroup taxa.
Ijerph 17 01538 g002
Ijerph 17 01538 g003 550
Figure 3. Unrooted network of 17 haplotypes in 65 A. triastriatus individuals from 35 locations (lineage A: blue; lineage B: green) based on 679 bp of a fragment of the mitochondrial COI gene. Each circle in the haplotype network corresponds to one haplotype, and its size is proportional to its frequency among samples. The numbers along branches represent the number of mutations between haplotypes. Dark dots are median vectors that represent intermediate haplotypes that were not sampled or are extinct.
Figure 3. Unrooted network of 17 haplotypes in 65 A. triastriatus individuals from 35 locations (lineage A: blue; lineage B: green) based on 679 bp of a fragment of the mitochondrial COI gene. Each circle in the haplotype network corresponds to one haplotype, and its size is proportional to its frequency among samples. The numbers along branches represent the number of mutations between haplotypes. Dark dots are median vectors that represent intermediate haplotypes that were not sampled or are extinct.
Ijerph 17 01538 g003
Ijerph 17 01538 g004 550
Figure 4. Divergence time estimate for A. triastriatus using a combined COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI database based on a relaxed molecular clock (uncorrected lognormal). Details of numbers in each node are shown in Table 4.
Figure 4. Divergence time estimate for A. triastriatus using a combined COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI database based on a relaxed molecular clock (uncorrected lognormal). Details of numbers in each node are shown in Table 4.
Ijerph 17 01538 g004
Ijerph 17 01538 g005a 550Ijerph 17 01538 g005b 550
Figure 5. Graphs of the mismatch distribution of lineage A (a) and lineage B (b) based on mitochondrial genes COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI. The abscissa represents the number of pairwise distances and the ordinate represents the number of observations. The lines represent the expected (smoothed) and observed (dashed) frequency of pairwise differences under the sudden lineage expansion model.
Figure 5. Graphs of the mismatch distribution of lineage A (a) and lineage B (b) based on mitochondrial genes COI, COII, ATP8, ND6, 12S rRNA, 16S rRNA and NDI. The abscissa represents the number of pairwise distances and the ordinate represents the number of observations. The lines represent the expected (smoothed) and observed (dashed) frequency of pairwise differences under the sudden lineage expansion model.
Ijerph 17 01538 g005aIjerph 17 01538 g005b
Ijerph 17 01538 g006 550
Figure 6. Bayesian skyline plots of population dynamic for populations of A. triastriatus. The shade areas in light blue are within the 95% highest posterior density interval. The dotted vertical lines indicate important population-size changes of A. triastriatus.
Figure 6. Bayesian skyline plots of population dynamic for populations of A. triastriatus. The shade areas in light blue are within the 95% highest posterior density interval. The dotted vertical lines indicate important population-size changes of A. triastriatus.
Ijerph 17 01538 g006
Ijerph 17 01538 g007 550
Figure 7. Dispersal routes of lineage A (left) and lineage B (right). Reconstructions are based on the network and genetic analysis. The arrows represent possible spread directions but not actual dispersal routes. Different colored arrows represent different branches in Figure 4.
Figure 7. Dispersal routes of lineage A (left) and lineage B (right). Reconstructions are based on the network and genetic analysis. The arrows represent possible spread directions but not actual dispersal routes. Different colored arrows represent different branches in Figure 4.
Ijerph 17 01538 g007
Ijerph 17 01538 g008 550
Figure 8. (A) Ventral view showing spermathecal pores, female pores and male pores of A. triastriatus usualis; (B) spermathecae of A. triastriatus usualis; (C) prostate glands of A. triastriatus usualis.
Figure 8. (A) Ventral view showing spermathecal pores, female pores and male pores of A. triastriatus usualis; (B) spermathecae of A. triastriatus usualis; (C) prostate glands of A. triastriatus usualis.
Ijerph 17 01538 g008
Table
Table 1. Collection data and information on A. triastriatus sequenced in this study using the COI gene. The code in bold formatting was sequenced using 7 genes.
Table 1. Collection data and information on A. triastriatus sequenced in this study using the COI gene. The code in bold formatting was sequenced using 7 genes.
TaxonLabelLocation
Code		Collection LocationGPS CoordinatesCollection DateElevation (m)Gene CodeNumber of Individuals
AC-AH201608-011Anhui30.12324° N, 118.09639° E05/201693AH86 AH60 AH72 AH81 AH91 AH9419
C-CQ201304-022Chongqing29.01467° N, 107.08238° E06/20131162CQ162
C-FJ201106-013Fujian26.10627° N, 118.42997° E07/2011342FJ101
C-FJ201111-014Fujian25.31278° N, 116.89306° E07/20111207FJ6518
C-GD201109-035Guangdong23.46417° N, 111.91444° E10/2011184GD564
C-GX201210-056Guangxi25.59723° N, 109.93519° E06/20121536GX3513
C-GX201318-057Guangxi24.15622° N, 110.21279° E05/20131198GX194 GX172 GX18920
C-GX201319-048Guangxi25.20595° N, 108.67555° E05/20131226GX198 GX2102
C-GZ201101-019Guizhou27.73722° N, 108.81611° E10/2011370GZ01 GZ41 GZ4427
C-GZ201103-0210Guizhou28.23472° N, 107.17444° E10/20111534GZ224
C-GZ201610-0111Guizhou27.53520° N, 108.59300° E10/2016604GZ142 GZ135 GZ13714
C-GZ201611-0512Guizhou27.39100° N, 108.11610° E10/2016890GZ153 GZ151 GZ1583
C-HB201101-0613Hubei31.47611° N, 110.38055° E09/20111399HB061
C-HU201304-0114Hunan25.94700° N, 113.70489° E06/2013982HN154
C-HU201311-0415Hunan27.27608° N, 112.70904° E06/2013500HN431
C-JX201313-0216Jiangxi26.54490° N, 114.18954° E06/2013571JX36 JX33 JX4010
C-SC201006-0217Sichuan29.64400° N, 102.12433° E08/20101625SC182
C-ZJ201511-0218Zhejiang28.83361° N, 119.24639° E07/2015160ZJ521
BC-AH201601-011Anhui29.95601° N, 118.09112° E05/2016325AH50 AH67 AH6910
C-AH201609-022Anhui30.31267° N, 118.02561° E05/2016120AH922
C-AH201610-013Anhui30.35924° N, 117.93609° E05/2016146AH100 AH101 AH10410
C-FJ201004-014Fujian26.11244° N, 118.42387° E08/2010952FJ094
C-FJ201111-045Fujian25.31278° N, 116.89306° E07/20111207FJ683
C-GD201115-046Guangdong23.46417° N, 111.91444° E08/2011410GD831
C-GX201319-067Guangxi25.20595° N, 108.67555° E05/20131226GX200 GX6312
C-GZ201609-028Guizhou27.52480° N, 108.59300° E10/2016731GZ136 GZ143 GZ1444
C-HU201309-029Hunan27.25680° N, 112.72433° E06/2013172HN34 HN042
C-JX201302-0410Jiangxi28.03037° N, 117.09680° E06/2013185JX062
C-JX201311-0111Jiangxi26.55213° N, 114.10773° E06/20131300JX26 JX50 JX60 JX62 JX8919
C-SC201004-0612Sichuan29.65556° N, 102.11667° E08/20101558SC154
C-ZJ201101-0213Zhejiang30.32886° N, 119.42599° E11/20111052ZJ021
C-ZJ201511-0114Zhejiang28.83361° N, 119.24639° E07/2015160ZJ511
C-ZJ201515-0315Zhejiang28.10444° N, 119.63472° E07/2015290ZJ83 ZJ874
C-ZJ201517-0216Zhejiang27.84556° N, 119.77306° E07/2015710ZJ931
C-ZJ201518-0117Zhejiang28.37639° N, 121.10972° E07/2015120ZJ98 ZJ1036
Table
Table 2. Primers used for PCR and sequencing.
Table 2. Primers used for PCR and sequencing.
GenePrimerSequenceSource
CO1LCO14905′-GGTCAACAAATCATAAAGATATTGG-3′[36]
HCO21985′-TAAACTTCAGGGTGACCAAAAAATCA-3[36]
CO1-E5′-TATACTTCTGGGTGTCCGAAGAATCA-3′[37]
CO2rRNA-Asn-CO2-tRNA-Asp: LumbF15′-GGCACCTATTTGTTAATTAGG-3′[33]
rRNA-Asn-CO2-tRNA-Asp: LumbR25′-GTGAGGCATAGAAATACACC-3′[33]
ATP8ATP8-F5′-GYTTAGTTRCCAMCYGGTGTATTTC-3′This study
ATP8-R5′-CTTYYTACTTGGAAGGTARRTGTAC-3′This study
ND6ND6-F5′-TGTATGGYGCACACRGGYHTTTGAA-3′This study
ND6-R5′-RATDGCTGGRGTNGGYTTAAACATA-3′This study
12S rRNA12S-tRNA-Val-16S-LumbF15′-CTTAAAGATTTTGGCGGTGTC-3′[33]
12S-tRNA-Val-16S-LumbR15′-CCTTTGCACGGTTAGGATAC-3′[33]
16S rRNA16Sar5′-CGCCTGTTTATCAAAAACAT-3′[38]
16Sbr5′-CCGGTCTGAACTCAGATCACGT-3′[38]
ND1rRNA-Leu-ND1-LumbF25′-GAATAGTGCCACAGGTTTAAAC-3′[33]
rRNA-Leu-ND1-LumbR1b5′-TTAACGTCATCAGAGTTATC-3′[33]
Table
Table 3. Summary of genetic diversity measures and neutrality tests between and within lineages based on the COI gene, and the mean number of pairwise distances based on a fragment of the COI gene and on a combination of seven genes.
Table 3. Summary of genetic diversity measures and neutrality tests between and within lineages based on the COI gene, and the mean number of pairwise distances based on a fragment of the COI gene and on a combination of seven genes.
TaxonNumber of IndividualsNumber HaplotypesPolymorphic SitesHaplotype DiversityNucleotide DiversityTajima’s DFu’s FSMean Number of Pairwise Distance
COI 7Gene
A (within lineage)3414160.8660.0052−1.028−4.6370.5%0.6%
B (within lineage)31320.1270.0002−1.506−2.3970.1%0.3%
Total (with lineage)6517490.7690.03092.303 *7.5256.3%6.7%
* p < 0.05.
Table
Table 4. Divergence times of A. triastriatus, estimated by the Bayesian method.
Table 4. Divergence times of A. triastriatus, estimated by the Bayesian method.
NodeEstimated Age in Millions of Years Ago (Ma)
AgeLower 95% Confidence LimitUpper 95% Confidence Limit
12.972.173.15
20.720.291.47
30.690.271.42
40.600.261.35
50.080.010.13
60.420.260.41
70.090.040.14
80.070.030.11
90.030.010.06
100.020.000.04
110.260.030.31
120.020.000.04
130.150.020.19
140.050.010.05
150.030.010.04
160.020.000.03
170.020.010.04
180.480.220.53
190.330.190.39
200.240.170.28
210.270.150.35
220.040.010.06
230.020.000.04
240.010.000.03
250.230.140.25
260.090.040.11
270.040.010.06
280.020.000.04
290.060.030.09
300.040.010.07
310.020.000.03
320.050.020.07
330.030.010.05
340.010.000.03
Table
Table 5. ANOVA results between and within lineages based on COI a gene fragment. d.f. = degrees of freedom, ΦST = ΦStatistic, FST = genetic differences among population.
Table 5. ANOVA results between and within lineages based on COI a gene fragment. d.f. = degrees of freedom, ΦST = ΦStatistic, FST = genetic differences among population.
Source of Variationd.f.Sum of SquaresVariance ComponentsPercentage of VariationΦST
Between lineages1629.61719.37529 Va93.890.9396
Within lineages6379.4451.26103 Vb6.11(p < 0.0001)
Total64709.06220.63632
FST0.9389

Share and Cite

MDPI and ACS Style

Dong, Y.; Jiang, J.; Yuan, Z.; Zhao, Q.; Qiu, J. Population Genetic Structure Reveals Two Lineages of Amynthas triastriatus (Oligochaeta: Megascolecidae) in China, with Notes on a New Subspecies of Amynthas triastriatus. Int. J. Environ. Res. Public Health 2020, 17, 1538. https://doi.org/10.3390/ijerph17051538

AMA Style

Dong Y, Jiang J, Yuan Z, Zhao Q, Qiu J. Population Genetic Structure Reveals Two Lineages of Amynthas triastriatus (Oligochaeta: Megascolecidae) in China, with Notes on a New Subspecies of Amynthas triastriatus. International Journal of Environmental Research and Public Health. 2020; 17(5):1538. https://doi.org/10.3390/ijerph17051538

Chicago/Turabian Style

Dong, Yan, Jibao Jiang, Zhu Yuan, Qi Zhao, and Jiangping Qiu. 2020. "Population Genetic Structure Reveals Two Lineages of Amynthas triastriatus (Oligochaeta: Megascolecidae) in China, with Notes on a New Subspecies of Amynthas triastriatus" International Journal of Environmental Research and Public Health 17, no. 5: 1538. https://doi.org/10.3390/ijerph17051538

APA Style

Dong, Y., Jiang, J., Yuan, Z., Zhao, Q., & Qiu, J. (2020). Population Genetic Structure Reveals Two Lineages of Amynthas triastriatus (Oligochaeta: Megascolecidae) in China, with Notes on a New Subspecies of Amynthas triastriatus. International Journal of Environmental Research and Public Health, 17(5), 1538. https://doi.org/10.3390/ijerph17051538

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