<title>Background</title> <p>Population genetic structure is dependent on the interaction of the biology of a species and the environment in which it resides. Marine organisms generally show low levels of genetic differentiation over large geographic distances [##UREF##0##1##,##UREF##1##2##]. In open marine systems it is evident that higher dispersal potential during planktonic egg, larval, or adult stages coupled with the absence of physical barriers to movement seems to greatly facilitate extensive gene flow among populations of organisms [##REF##10879524##3##,##UREF##2##4##]. However, populations of freshwater fish species from different basins often show significant genetic differentiation resulting from isolation, while populations within a basin often show no or low levels of genetic differentiation [##UREF##3##5##, ####UREF##4##6##, ##UREF##5##7##, ##REF##15964769##8##, ##UREF##6##9####6##9##]. There is evidence that for some species instream barriers can restrict gene flow and induce high levels of genetic differentiation between populations [##REF##11091310##10##,##UREF##7##11##]. Even in the absence of restricted gene flow, multiple unrelated founding events could also shape the same pattern of population genetic structure [##UREF##2##4##].</p> <p>Due to periodic climatic oscillations during the Pleistocene, range contractions, range expansions and changes in interconnectivity within and between drainages during the periodic climatic oscillations of the Pleistocene have greatly influenced the distribution of many fish species, have significantly altered the amount and distribution of intraspecific genetic variation in many species and have remodelled the population structures of temperate fish faunas [##UREF##0##1##,##UREF##3##5##]. The major climatic oscillations occurred during the past 800 kyr (with a 100 kyr dominant cycle), and the last glacial maximum (LGM) occurred about 18 kyr ago with a decline of the sea levels of about 120–140 m [##REF##12226674##12##]. Each of the severe climatic shifts and the change in the sea levels could have produced great changes in the freshwater fish species' geographical distribution and abundance [##UREF##0##1##,##UREF##3##5##,##REF##10434419##13##,##REF##10679153##14##]. The legacy of these changes in the phylogeographic distribution of genetic diversity has been traced in several species [##UREF##3##5##,##UREF##8##15##,##REF##15522807##16##]. In Cyprinids, it has been shown that glacial, interglacial and postglacial oscillations have resulted in major changes of distributions and have shaped current distributions by changes of river courses and river confluences in response to falling sea levels [##REF##10434419##13##,##REF##10679153##14##].</p> <p>The Family Salangidae, which is classified as Protacanthopterygii, Order Osmeriformes, Suborder Osmeroidei [##REF##15737583##17##,##REF##12079643##18##], comprises six genera and approximately 17 species [##UREF##9##19##,##UREF##10##20##]. The Salangid, <italic>Neosalanx taihuensis</italic>, is one of the most economically important species in this Family. This species is endemic to China and is restricted to large bodied freshwater systems such as inland lakes, out-flowing rivers and estuaries, covering a continuous distribution from the temperate (Yellow River basin) to the subtropical zone (Pearl River basin) [##UREF##11##21##]. Also, this species exhibits annualism poor swimming ability and has some neotenic features such as transparent miniaturized bodies, cartilaginous endoskeletons, and notochords [##UREF##10##20##,##UREF##12##22##]. The combination of relative large spatial distribution encompassing discrete drainages in tropical and subtropical regions and the low vagility constitutes a good model system to study the palaeoclimatic effects on population genetic structure and phylogeography of freshwater [##UREF##13##23##]. In recent years, because of water pollution, overfishing and habitat destruction, the population of <italic>N. taihuensis </italic>has undergone the lost of appropriate habitat and the decrease in population size [##UREF##14##24##,##UREF##15##25##]. In order to protect the resource of <italic>N. taihuensis</italic>, the habitat detection has to be implemented immediately, and the study on the <italic>N. taihuensis </italic>population structure and its phylogeoraphic pattern should be enhanced since the failure to detect population units will lead to local overfishing and ultimately to severe declines [##REF##12535089##26##]. The Salangidae family has been studied for more than 200 years, but most of studies have focused on taxonomy [##UREF##16##27##,##UREF##17##28##], biology [##UREF##18##29##], and molecular phylogeny [##UREF##19##30##]. Previous studies in China have shown that in two unrelated species,<italic>Opsariichthys bidens </italic>(Teleostei, Cyprinidae) and <italic>Brachymystax lenok </italic>Pallas (Salmoninae, Salmonidae), both had genetic patterns coinciding with contemporary drainage structure [##REF##15964769##8##,##UREF##6##9##], but the related study on <italic>N. taihuensis </italic>is not conducted till now.</p> <p>The basins of Yangtze River, Huai River and Pearl River, in which <italic>N. taihuensis </italic>mainly distributes, are important river systems in China. Although a strong uplifting of the Tibetan plateau occurred during the Pleistocene in southwest China [##UREF##20##31##], there is no evidence to indicate that large-scale river rearrangement occurred after the Pliocene in the mid-lower reaches of these basins [##UREF##21##32##]. The Yangze River and the Huai River basin are connected by the Peking Hangzhou Canal, while there are no connection between Pearl River basin and the others[##UREF##22##33##]. Given the intolerance of the <italic>N. taihuensis </italic>to waters with some salt content, members of freshwater fish division cannot disperse by river connections [##UREF##23##34##]. Therefore, we expect vicariance among basins might play an important role in shaping population genetic structure of <italic>N. taihuensis</italic>. In addition, the glaciations during the Pleistocene, especially the last glacial maximum (LGM) with sea levels lowered (120–140 m below present sea level), should have severely affected the genetic structure of this species. In order to detect the population structure and phylogeographic pattern of <italic>N. taihuensis </italic>at large geographical scales, samples of <italic>N. taihuensis </italic>were collected from 13 populations among the Pearl River basin, the Yangtze River basin, and the Huai River basin (Figure. ##FIG##0##1##). The complete mtDNA cytochrome <italic>b </italic>gene (<italic>cyt b</italic>) sequence was used to analyze population genetic differentiation of this species, population dynamics and demographic history, and to show how the glaciations influenced the phylogeographic patterns of freshwater fish, and how this species responded to the severe climatic oscillations and decline of the sea levels during the Pleistocene ice ages.</p>
<title>Methods</title> <title>Sample Collection</title> <p>We collected a total of 354 individuals of <italic>N. taihuensis </italic>from 13 populations in the Pearl River, Yangtze River, and Huai River basins from March in 2004 to September in 2005, 44 came from two populations of the Pearl River basin in 2005 (P1–P2), 196 from seven populations of the Yangtze River basin in 2004 (Y1–Y7), and 114 from four populations of the Huai River basin in 2005 (H1–H4) (Figure. ##FIG##0##1## and Table ##TAB##0##1##). To eliminate the possibility that these samples came from the related individuals (schools of fishes), we selected many patches from one location and resampled 2–4 times in same population. All samples (whole fish) were stored in 95% ethanol.</p> <title>DNA extraction, PCR amplification and Sequencing</title> <p>Total DNA was extracted from the muscle tissue using standard phenol/chloroform procedures [##UREF##41##70##] or following the method of Asahida <italic>et al</italic>. [##UREF##42##71##]. PCR was performed according to the method of Zhang <italic>et al</italic>. [##UREF##19##30##]. 1141 base pairs (bp) of <italic>cytb </italic>DNA were amplified using the primers L14321 (5'-CAGTGACTTCAAAAACCACCG-3') and L15634 (5'-CTTAGCTTTGGGAGTTAAGGGT-3') in a 50 μl reaction volume mixture containing 25 μl <italic>Premix Taq </italic>(1.25 U EX Taq polymerase, 0.4 mM of each dNTP Mixture, 4 mM Mg2+; TaKaRa, Tokyo, Japan), 1.0 μM of each primer, approximately 50–100 ng total genomic DNA template and 18–22 μl of sterile distilled water. A PE 9700 Thermal Cycler (Applied Biosystems) was used for the following cycles: pre-denaturation at 94°C for 5 min; 30 cycles of denaturation at 94°C for 30 s, annealing at 54°C for 45 s, and extension at 72°C for 1 min 10 s, plus a final extension at 72°C for 10 min. Purified PCR products were directly sequenced using primers L14321 and L15634 from both ends using the Big Dye Terminator Sequencing Kit (Perkin-Elmer, Norwalk, CT) in a semi-automated DNA analyzer (3700; Applied Biosystems). To avoid the errors in amplification and sequencing, all singletons among polymorphism sites were verified by an additional amplification and sequencing. Sequences were aligned using the Clustal X program [##REF##9396791##72##] and rechecked against the inferred reading frame for the corresponding protein. Haplotypes were identified using Clustal X, and the <italic>cytb </italic>haplotype sequences of <italic>N. taihuensis </italic>were deposited in GenBank (<ext-link ext-link-type="gen" xlink:href="EU376454">EU376454</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU376489">EU376489</ext-link>).</p> <title>Data Analysis</title> <title>Genetic Diversity and Population Structure</title> <p>In order to test whether there is recombination within <italic>cyt b </italic>gene, linkage disequilibrium (LD) test was conducted using DnaSP 4.10 [##REF##14668244##73##] and the average degree of LD, or non-random association between nucleotide variants at different polymorphic sites was tested using ZZ statistics [##REF##11454763##74##], which could be used for detecting intragenic recombination with the null hypothesis no recombination. Also, neutrality test for <italic>cyt b </italic>was conducted in MEGA4.0 [##REF##17488738##75##] to check selection by Z-test using the null hypothesis H0: d<sub>N </sub>= d<sub>S </sub>[##REF##3444411##76##,##REF##9520431##77##], where d<sub>S </sub>and d<sub>N </sub>is the average number of synonymous substitutions per synonymous site and the average number of nonsynonymous substitutions per nonsynonymous site, respectively.</p> <p>Genetic diversity was measured for all samples and for each basin grouping using haplotype diversity (<italic>h</italic>) and nucleotide diversity (<italic>π</italic>) [##UREF##43##78##]. Values for the numbers of haplotypes (<italic>H</italic>), polymorphic sites (<italic>S</italic>) and the mean numbers of pair-wise differences among sequences (<italic>K</italic>) were also estimated. These diversity indices were computed using the software DnaSP 4.10 [##REF##14668244##73##]. As unequal numbers of individuals varying from 13 to 34 per population in this study, Pearson's correlation tests were conducted to test values independence of the number of haplotypes (<italic>H</italic>), haplotype diversity (<italic>h</italic>) and nucleotide diversity (<italic>π</italic>) against the number of individuals (<italic>N</italic>). Population structure and genetic variation of <italic>N. taihuensis </italic>were characterized and compared with Arlequin version 3.11 [##REF##19325852##79##]. Analysis of molecular variance (AMOVA) was used to asses the population configuration and the geographical pattern of population subdivision. In this study, populations were grouped according to different geographical hierarchies using <italic>Φ</italic>-statistics [##UREF##44##80##]. Three hierarchical levels of subdivision were obtained: <italic>Φ</italic><sub><italic>CT</italic></sub>, the degree of differentiation among all basins, <italic>Φ</italic><sub><italic>SC</italic></sub>, the degree of differentiation among populations within basins, and <italic>Φ</italic><sub><italic>ST</italic></sub>, the degree of differentiation among all populations. We tested whether the derived indices were significantly different from zero using a nonparametric permutation method (10,000 Permutations). Moreover, we estimated pair-wise genetics differentiation between populations with <italic>Φ</italic><sub><italic>ST </italic></sub>[##UREF##44##80##,##REF##1644282##81##] that includes information on haplotypes frequency and information on haplotype sequences. Null hypothesis of genetic homogeneity was assessed by 10,000 replications and sequential Bonferroni corrections [##UREF##45##82##] for multiple comparisons were applied to all pairs comparisons. To verify the hypothesis of isolation-by-distance (IBD), correlation between pair-wise linearized Φst/(1-Φst) and river distances (Km) between populations were analyzed using the Mantel test [##REF##6018555##83##] with 10,000 permutations. Mantel test were conducted for all 13 populations among basins, for 7 populations within Yangze River basin and for 4 populations within Huai River basin, respectively. The river distances between populations were determined based on river courses by AcrView 3.2, which was also used in NCA analysis.</p> <title>Phylogenetic and Phylogeographic Analysis</title> <p>For phylogenetic analysis, we constructed unrooted maximum-likelihood (ML) tree using the program PAUP* 4.0 [##UREF##46##84##]. The TrN+I (I = 0.8148) model was selected as the best-fit model for analysis using MODELTEST 3.06 [##REF##9918953##85##].</p> <p>Intraspecific data typically consists of many similar sequences, some of which may be ancestral and its phylogenetic relationship are often more clearly and accurately represented by a network [##REF##11146143##86##]. In the present study, we constructed a haplotype network based on statistical parsimony [##REF##1385266##87##] using the program TCS version 1.18 [##REF##11050560##88##]. A distance matrix for all pair-wise haplotype comparisons was constructed and the maximum number of mutational differences justified by the parsimony limit of 0.95 was estimated. The network first constructed haplotypes that differed by a single change and added increasingly more distant haplotypes until either all were included or the maximum number of mutational steps was reached [##REF##11050560##88##]. To check the evolutionary mechanisms responsible for the spatial distribution of genetic variation, we conducted the Nested Clade Analysis (NCA). The resulting cladogram calculated by the TCS program was converted into a nested design using the nesting rules [##REF##8100789##89##]. Reticulations, or equally parsimonious connections within the network, were resolved with two steps. First, alternative connections between haplotypes were broken following a series of rules based on coalescence theory [##REF##8349118##90##]. Second, synapomorphies in the form of nonsynonymous substitutions in the <italic>cytb </italic>gene, a conservative class of substitutions, was used to resolve all cases in which assignment to a nested series was ambiguous [##REF##8100789##89##]. The current and historic patterns of phylogenetic and geographic associations were statistically tested using NCA implemented in GeoDis 2.0 [##REF##10736051##91##], under the null hypothesis of no geographic association among haplotypes using 10,000 permutations. The most suitable model explaining the historical geographical distribution was identified using the inference key [##REF##15012756##51##]. It is appropriate to address the caveats of NCA analysis, the power of this method has been criticized for years [##REF##12453245##92##, ####UREF##47##93##, ##REF##17542853##94####17542853##94##], but Templeton defended NCA on both theoretical and empirical grounds, and insisted this method an extensively validated method for strong phylogeographic inference [##REF##15012756##51##,##REF##18346121##95##]. However, as no other better all-encompassing method can offer the ability to explore patterns relating to complex historical scenarios at present, the NCA method can still provide useful complementary information for phylogeographical analysis [##REF##18284568##96##].</p> <title>Demographic Analysis</title> <p>To detect mechanisms responsible for observed mtDNA patterns on various spatial and time scales, we developed different data analysis methods. Due to the limits and pitfalls of each method, the combined results might provide the best approximation of the population structure and history.</p> <p>First, changes over time in nonparametric estimates of the effective population size of <italic>N. taihuensis </italic>were evaluated with the skyline plots [##REF##11719579##38##,##REF##10880500##39##] using R 2.6.1 [##UREF##48##97##]. This method can provide a rough profile of population demographic fluctuation, and can be used as a model selection tool. The classic skyline plot typically produces 'noisy' plots that display the stochastic variability inherent in the coalescent process. While using the Akaike Information Criterion [##UREF##49##98##], the generalized skyline plot could reduce this noise, and thus produce smoother estimated population size plots [##REF##11719579##38##]. As the power of skyline plots depends on the numbers of variable nucleotide sites [##REF##11719579##38##,##REF##15703244##55##], we conducted skyline analysis for the entire region, the Yangtze River basin and the Huai River basin with exception of the Pearl River basin. Second, for the entire region and for each basin, we examined the observed distribution of pair-wise differences between sequences (mismatch distribution) [##REF##1316531##53##]. Theoretical studies show that mismatch distributions are usually ragged or multimodal for populations at stationary demographic equilibrium, but are typically smoother or unimodal for populations that have recently undergone a demographic expansion [##REF##1316531##53##]. Occasionally, a pronounced rate of heterogeneity is reported even among closely related lineages [##REF##16503171##59##], so the pair-wise relative rate test (pRRT) implemented in HyPhy [##UREF##50##99##] was used to estimate and verify molecular clock constancy at the intraspecific level based on 1141 bp of 36 haplotypes and <italic>P. chinensis </italic>(DQ191115) [##UREF##19##30##] before performing the mismatch distribution analysis. One haplotype (H31) was identified to evolve at a significantly different pace from the others. After removing this haplotype (H31), no other violation of the molecular clock constancy was found, so the other 35 haplotypes were used for further analysis. Mismatch analysis was conducted using Arlequin 3.11 [##REF##19325852##79##] under a model of population expansion. The overall validity of the estimated demographic model was evaluated by the tests of raggedness index (<italic>Hri</italic>) [##REF##8088750##100##] and the sum of squared differences (<italic>SSD</italic>) [##REF##15836643##101##]. Significance of <italic>Hri </italic>and <italic>SSD </italic>were assessed by parametric bootstraps (10,000 replicates), and the significant value was taken as evidence for departure from the estimated demographic model of sudden population expansion. Third, Fu and Li's <italic>D</italic>* [##REF##8454210##42##], Fu's <italic>Fs </italic>[##REF##9335623##40##] and Ramos-Onsins and Rozas's <italic>R</italic><sub>2 </sub>[##REF##12446801##41##] tests for mutation/drift equilibrium were performed in DnaSP 4.10 [##REF##14668244##73##] and Arlequin 3.11 [##REF##19325852##79##] with 10,000 simulations. Although these methods are commonly used to test the selective neutrality of genetic markers, these estimators are also sensitive to demographic processes such as recent population expansion or bottleneck. Because we were interested in discriminating between demographic expansion and contraction, we chose two class test statistics, each with particular sensitivity to one demographic scenario. Fu and Li's <italic>D</italic>* is designed to detect an excess of old mutations, characteristic of a population that has experienced a historical reduction in effective population size [##REF##8454210##42##,##REF##8722804##102##]. In contrast, Fu's <italic>F</italic><sub><italic>S </italic></sub>and <italic>R</italic><sub>2 </sub>are sensitive to an excess of recent mutations [##REF##9335623##40##,##REF##12446801##41##]. Finally, we used the coalescent-based method implemented in FLUCTUATE 1.4 [##REF##9584114##43##] to estimate exponential growth rate (<italic>g</italic>) for this species. All runs employed the following strategy: 10 short chains of 4,000 steps and five long chains of 400,000 steps, sampling every 20th step; random starting trees; empirical nucleotide frequencies; initial <italic>g </italic>value of 0.0, starting <italic>θ</italic>-value from Watterson's estimate [##REF##1145509##103##]; and a 2.424 transition/transversion (ti/tv) rate determined from MEGA4 [##REF##17488738##75##] under the Tamura-Nei model [##REF##8336541##104##] using only ingroup sequences. Runs were repeated five times to ensure consistency of estimates. As the estimates of <italic>g </italic>are biased upwards [##REF##9584114##43##], a conservative approach in testing for significance was adopted in our study, with values larger than three standard deviations (<italic>SD</italic>) of <italic>g </italic>regarded as significant.</p> <p>For stepwise expansion population, we used <italic>τ </italic>to calculate the expansion time. The mismatch distribution analysis provided a rough estimate of <italic>τ</italic>, with the starting time of the expansion in units of 1/(2 ut) generations, where u is the mutation rate per locus per generation. The relationship with the absolute time in years (t), is t = <italic>τ</italic>/2 uT. The value of T is generation time (1 year for <italic>N. taihuensis</italic>), and the value of u is derived from u = μk, where μ is the mutation rate per nucleotide per generation, and k is the number of nucleotides in the sequence. For other populations, however, we used GENETREE [##REF##10792974##44##] to estimate the time to the most recent common ancestor (TMRCA). GENETREE is based on a coalescent method to simulate gene trees conditional on their topology, which can be used to determine the distribution of TMRCA. In present study, the most probable ancestral haplotype was analyzed by TCS program, and the TMRCA for a constant population was calculated in GENETREE by first finding the maximum-likelihood estimate [##REF##16689893##105##] of <italic>θ </italic>following the method of Joy <italic>et al</italic>. [##REF##12690197##106##]. Each run was repeated three times with a different starting seed number and 1,000,000 coalescent simulations. To get a more detailed picture of the mutational pattern from population expansion, GENETREE was also used to construct a gene tree to determine the distribution of each mutation for <italic>N. taihuensis</italic>.</p> <p>Since there is no fossil record for <italic>N. taihuensis </italic>to analyze their cartilage character, the genetic time clock is difficult to calibrate. However, it has been convincingly shown that most groups of fish have slower evolutionary rates compared with those of other species, for example primates [##REF##7807548##35##]. A calibration rate of 0.8–1% per Myr per generation has been accepted for mitochondrial protein-coding genes in Salmoniformes or Osmeriformes in many studies [##UREF##41##70##,##UREF##51##107##,##REF##11123621##108##]. For convenience, we used the mutation rate μ = 1% in this study.</p>
<title>Results</title> <title>Sequence Variation and Genetic Diversity</title> <p>For 1141 bp of the complete <italic>cyt b </italic>gene analyzed in 354 individuals, there were 38 variable nucleotide sites, and no insertions or deletions in any of the sequences. The average nucleotide composition for all individuals was highest for C (0.334), followed by T (0.277), A (0.214), and G (0.175). Nucleotide composition showed an anti-G bias (G = 17.5%), which is characteristic for the mitochondrial genome [##REF##7807548##35##]. Since there were no stop codons when the <italic>cytb </italic>sequence was translated into amino acid sequences, and based on the criteria of Zhang &amp; Hewitt [##UREF##24##36##], we found no evidence for nuclear mitochondrial pseudogenes (Numts) in our PCR sequences. Linkage disequilibrium (LD) test and selective neutrality test were implemented for <italic>cyt b </italic>sequences in <italic>N. taihuensis</italic>, which indicated that there was no signs of recombination (<italic>P </italic>= 0.925) and the neutral hypothesis (no selection) could not be rejected (<italic>P </italic>= 0.172) for these <italic>cyt b </italic>sequences. Thus, the <italic>cyt b </italic>was suitable marker for analyzing the evolutionary history of this species.</p> <p>Thirty-six different haplotypes were identified in the 354 samples analyzed (Table ##TAB##0##1## and Table ##TAB##1##2##). The number of haplotypes ranged from 5 to 25 for each basin. The populations with the greatest number of haplotypes were the Taihu lake population (8) and Chaohu lake population (8) in Yangze River basin, and the Hongzehu lake population (8) in Huai River basin. The haplotype frequency distribution was strongly skewed, with the vast majority of haplotypes found only once (24 out of 36) and restricted to a single population. H36, the most common haplotype found in 12 sampled populations with the exception of the Xujiahe reservoir population, occupied 49.15% of all individuals among all populations (56.82% for the Pearl River basin, 37.43% for the Yangtze River basin, 72.57% for the Huai River basin). Other common haplotypes found in all three basins were H27, H33 and H25, with 16.38%, 10.73% and 3.11% of all individuals, respectively (Table ##TAB##1##2##).</p> <p>Diversity indices (average ± standard deviation), <italic>h </italic>and <italic>π </italic>, are summarized in Table ##TAB##0##1##. The total of <italic>h </italic>was 0.713 ± 0.022 ranging from 0.473 ± 0.057 to 0.781 ± 0.017, and <italic>π </italic>was 0.0022 ± 0.0001 ranging from 0.0010 ± 0.0002 to 0.0027 ± 0.0001 among there basins. There are no evidences that heterogeneous levels of genetic variability were produced by unequal sample size for each population through Pearson's correlation test between diversity indices and sample size (r = 0.082 and <italic>P </italic>= 0.790 for numbers of haplotypes, r = -0.02 and <italic>P </italic>= 0.996 for haplotype diversity, and r = 0.141 and <italic>P </italic>= 0.645 for nucleotide diversity, respectively). These results showed a medium/high haplotype diversity and a low nucleotide diversity, and that the genetic diversity for three near-coastal populations (the Taihu lake population and the Chaohu lake population within the Yangze River basin, and the Hongzehu lake population within the Huai River basin) were higher than others (Table ##TAB##0##1##).</p> <title>Population Structure and Phylogeography</title> <p>Unrooted ML tree revealed that there were no distinct haplotype groups with high bootstrap support (&lt;60 %) and there was no geographical structure among haplotypes (Figure. ##FIG##1##2##), which was supported by network analysis (Figure. ##FIG##2##3##). In this study, only a few missing haplotypes were detected, which indicated that we gathered enough samples (both numbers of individuals and locations within the distribution) for phylogeographical analysis. Each haplotype was connected with others by 1 to 3 mutational steps, suggesting that there was no deep branching among halotypes for this species. In network, H36 yielded the highest outgroup weight (0.2941) followed by H27 (0.0588), H33 (0.0588), and H35 (0.0441), indicating H36 was the most probable ancestral haplotype [##REF##8075830##37##]. In fact, the high frequency of shared haplotypes (i.e., H36, H27, H33, H29) were found among all basins, which produced a pattern of population relationships associated with the ancestry of the haplotypes found within each population rather than with the hydrological network (Figure. ##FIG##1##2##, Figure. ##FIG##2##3##). Thus, we suggested that there was no obvious phylogeographic pattern within <italic>N. taihuensis</italic>.</p> <p>Pair-wise comparisons <italic>Φ</italic><sub><italic>ST </italic></sub>test among 13 populations showed that the differences among the most of populations were irrelevant to their distributions (Table ##TAB##2##3##). AMOVA analysis revealed that there were significant subdivisions within and among populations of <italic>N. taihuensis</italic>, but no significant genetic variance was found among basins, with only 8.43% of the genetic variance was found among basins (Table ##TAB##3##4##), which provided evidences that the high level of population structure was not related to the hydrological pattern. The Mantel test indicated no significant relationship between Φst/(1-Φst) and geographic distance among populations, with <italic>P </italic>= 0.59 for all 13 populations among basins, <italic>P </italic>= 0.82 for 7 populations within Yangze River basin and <italic>P </italic>= 0.19 for 4 populations within Huai River basin, respectively.</p> <p>The NCA analysis (see additional file ##SUPPL##0##1##) provided some historical information for the pattern of genetic differentiation for <italic>N. taihuensis</italic>. The widespread clade 1–6 with significantly small <italic>Dc </italic>and large <italic>Dn</italic>, were indicative long-distance colonization (and/or past fragmentation). Clade 1–2, clade 1–4, and two second cladistic levels of clade 2-1, clade 2-2 indicated restricted gene flow with isolation by distance, and the overall picture of NCA indicated historical contiguous range expansion.</p> <title>Demographic Analysis</title> <p>An examination of demographic histories revealed the marked differences among the basins under study (Table ##TAB##4##5##, Figure ##FIG##3##4##, Figure ##FIG##4##5##). The skyline plots [##REF##11719579##38##,##REF##10880500##39##]of <italic>N. taihuensis </italic>in the entire region and the Yangze River basin showed a sudden stepwise expansion. Fu's <italic>Fs </italic>[##REF##9335623##40##], Ramos-Onsins and Rozas's <italic>R</italic><sub>2 </sub>[##REF##12446801##41##] tests for the entire region and Fu's <italic>Fs </italic>test for Yangze River basin were statistically significant negative, supported this point. Although both <italic>Fs </italic>and <italic>R</italic><sub>2 </sub>tests could not effectively tell the causes of bottleneck or expansion, Fu and Li's <italic>D</italic>* [##REF##8454210##42##] tests were not significant (<italic>P </italic>&gt; 0.05) for these populations, suggesting that there were no historical reduction in effective population size in these regions. The mtDNA mismatch analyses for the entire region and the Yangze River basin showed bimodal profile which might result from constant population size among an old population or an admixture population (Figure ##FIG##4##5##). However, the mismatch distribution goodness of fit test (Table ##TAB##4##5##) was not significant, which indicated that there were no severe departure from the estimated demographic model. Moreover, the results of the coalescent based analysis using FLUCTUATE1.4 [##REF##9584114##43##] showed high growth rate, g = 2648.97 ± 556.26 for the entire region and <italic>g </italic>= 3268.21 ± 873.30 for the Yangze River basin, providing powerful evidences of population expansion in these regions. In contrast to that of the entire region and the Yangze River basin, the skyline plot for Huai River basin provided evidence that the exponential growth model fitted the population demographic fluctuation well (Figure ##FIG##3##4##). Significant Fu's <italic>Fs </italic>(<italic>P </italic>&lt;0.05) and Ramos-Onsins and Rozas's <italic>R</italic><sub>2 </sub>(<italic>P </italic>&lt; 0.05) test, as well as high growth rate (<italic>g </italic>= 2083.14 ± 529.35) confirmed population growth in Huai River basin, which was also supported by the test of <italic>Hri </italic>(<italic>P </italic>&gt; 0.95), although the bimodal mismatch distribution (Figure. ##FIG##4##5##) and significance of <italic>SSD </italic>(<italic>P </italic>&lt; 0.05) indicated a poor fit for the stepwise growth model (Table ##TAB##4##5##). For the Pearl River basin, the bimodal mismatch distribution (Figure. ##FIG##4##5##) and significance of <italic>SSD </italic>(<italic>P </italic>&lt; 0.05) values indicated a relative constant population size. Fu and Li's <italic>D</italic>* test, Fu's <italic>Fs </italic>and <italic>R</italic><sub>2 </sub>test for Pearl River basin were not significant, which also rejected population expansion/bottleneck model (Table ##TAB##4##5##). In addition, a relative constant population size was confirmed by a relative low growth rate (g = -309.64 ± 955.56) with the approximate 95% confidence interval of <italic>g </italic>included zero in this region.</p> <p>According to the tau value (τ) and the mutation rates for the <italic>cyt b </italic>gene (1% per nucleotide per generation per million years), the estimated expansion time was found to be 245.5 kyr ago for the entire region and 221.9 kyr ago for the Yangze River basin. Using the coalescent method, an estimate of TMRCA was 190.0 kyr for this species in the Huai River basin and 73.53 kyr in the Pearl River basin, respectively (Table ##TAB##4##5##).</p> <p>To get a more detailed picture of the mutational pattern for population expansion of <italic>N. taihuensis </italic>within the Yangze River basin and the Huai River basin, GENETREE [##REF##10792974##44##] was used to construct a gene tree for determining the distribution of mutation time. To make the sequence data suitable for the infinite-sites model, four haplotypes (H06, H07, H10, and H13) were excluded in this analysis. It is worth to note that the possible removal of data can have important effects on inference, and can lead to unusual conclusions [##REF##10792974##44##], but the results from this procedure can provide a rough profile for the distribution of mutations. As shown in Figure. ##FIG##5##6##, it was evident that the vast majority of mutations occurred near the tips of the gene tree, with 61.54% mutations occurring within 15.6 kyr in the Yangze River basin and 66.67% mutations occurring within 16.2 kyr in the Huai River basin.</p>
<title>Discussion</title> <title>Population Genetic Structure and Phylogeography</title> <p>Our results revealed low nucleotide diversity (<italic>π </italic>= 0.0022 ± 0.0001) and medium/high haplotype diversity (<italic>h </italic>= 0.713 ± 0.022) in <italic>N. taihuensis </italic>(Table ##TAB##0##1##), which could reflect a recent expansion or a short evolutionary history of the population [##UREF##0##1##,##UREF##25##45##, ####UREF##26##46##, ##UREF##27##47####27##47##]. This scenario was also supported by the absence of deep branching among haplotypes in the genealogical tree and the "star-like" shaped network (Figure. ##FIG##1##2## and Figure. ##FIG##2##3##) as well as NCA and demographic history analysis (Table ##TAB##4##5##).</p> <p>Phylogenetic inference based on unrooted ML tree revealed that there were no distinct clades with high bootstrap support (&lt;60%) and individuals from three basins did not partition into distinct clades (Figure. ##FIG##1##2##), which indicated no geographical structure among haplotypes. The statistical parsimony network displayed a "star-like" shape with H36 being the most probable ancestral haplotype (outgroup weight: 0.2941) (Figure. ##FIG##2##3##), which also revealed a lack of geographical structure.</p> <p>Beyond the prediction that vicariance among basins should play an important role in shaping population genetic structure, the results of AMOVA showed that the contemporary drainage structure served as a poor model in explaining the high levels of genetic differentiation, with only 8.43% variance among basins, it was suggested that there must be other factors shaping the genetic structure of <italic>N. taihuensis</italic>. Our results were in contrast to the 'Stream Hierarchy Model' [##UREF##28##48##], which showed that the distribution of genetic variation would follow different basins. Similar results were found in many studies of freshwater fish species [##UREF##3##5##,##UREF##6##9##,##REF##16702764##49##,##UREF##29##50##], including <italic>Opsariichthys bidens </italic>(Teleostei, Cyprinidae) conducted in the Pearl River and the Yangze River in China [##REF##15964769##8##]. However, some researchers [##REF##11091310##10##,##UREF##7##11##,##REF##15522807##16##] have reported population genetic structure in some freshwater fish species being analogous to our study, and these provided evidence that the contemporary drainage structure did not coincide with the genetic relationships among populations, or that the instream barriers and complex population histories played an important role in the determination of the population genetic structure.</p> <p>Pair-wise comparisons between <italic>Φ</italic><sub><italic>ST </italic></sub>and AMOVA analysis confirmed the existence of some degree of restricted gene flow, which could be one of the important factors determining the genetic structure among populations of <italic>N. taihuensis</italic>. However, the Mantel test indicated no significant relationship between genetic differentiation and geographic distance among populations. The discrepancy between these results indicated that the pattern of genetic differentiation for <italic>N. taihuensis </italic>was complex and could not be interpreted simply using the IBD model.</p> <p>In contrast, NCA analysis revealed that multiple unrelated long-distance founding events could be important in determining the patterns of genetic diversity. In NCA, the large <italic>Dn </italic>and small <italic>Dc </italic>suggests long distance dispersal [##REF##15522807##16##]. The NCA inferences for widespread clade 1–6, including the most probable ancestral haplotype H36 with a significantly small <italic>Dc </italic>and large <italic>Dn</italic>, were indicative of range expansion and long-distance colonization (and/or past fragmentation) [##REF##15012756##51##]. The observed population relationships agreed with the population structure predicted by normal dispersal models of recently expanding populations, with a portion of individuals migrating beyond neighbouring demes because of rare long-distance migrants founding pocket populations under this dispersal scenario [##REF##15522807##16##]. Also, theoretical study has proven that even though no effectively restricting gene flow exists under these conditions, the effects on these founder-events can persist for thousands of generations, which results in higher levels of divergence within than between rivers [##UREF##2##4##]. In addition, it was also suggested that historical basin rearrangements might form the analogous population genetic structure mentioned above [##REF##11091310##10##,##REF##15522807##16##]. Although a strong uplifting of the Tibetan plateau occurred during the Pleistocene in southwest China [##UREF##20##31##], there is no evidence to indicate that large-scale river rearrangement occurred after the Pliocene in our study area [##UREF##21##32##]. It may imply that historical basin rearrangements are not critical in determining the current population genetic structure of this species. Accordingly, we suggested that the current complex population genetic structure of this species was mainly constituted by multiple unrelated founding events (long-distance colonization) of dispersal, and also affected by contiguous population expansion and some degree of restricted gene flow.</p> <title>Demographic history</title> <p>Statistical analysis of sequences was originally developed to test selective neutrality of mutations and has recently been used to test for population expansion. Demographic events are also commonly analyzed by using the distribution of pair-wise differences or mismatch distributions of non-recombining DNA sequence such as mtDNA [##REF##1743491##52##,##REF##1316531##53##]. Comparative studies on the test power of various statistics suggested the superiority of <italic>Fs </italic>and <italic>R</italic><sub>2 </sub>over other statistics [##REF##9335623##40##,##REF##12446801##41##]. The power of these tests also depends on sample size, and the <italic>R</italic><sub>2 </sub>test is more powerful for a small sample while <italic>Fs </italic>is more powerful for a large sample [##REF##12446801##41##]. In this study, <italic>Fs </italic>tests were significant for the entire region, Yangze River basin and Huai River basin. However, <italic>R</italic><sub>2 </sub>tests were significant for the entire region and Huai River basin, but not significant for Yangze River basin (Table ##TAB##4##5##). The different results for population expansion in this study may reflect the differences in sample size requirements of the tests as well as source information.</p> <p>Mismatch distribution analysis was a common method to infer population demography; however, some of inferences of the mismatch analysis were not consistent with the results of the other indices in this study (Table ##TAB##4##5##). Besides a relative lower statistical power of demography analysis [##REF##12446801##41##], other factors might also induce such scenario. First, the single stepwise expansion model in mismatch analysis might be inadequate for some populations [##UREF##30##54##], and the use of an incorrect demographic model can also lead to biased and invalid estimates of demographic history and other evolutionary parameters [##REF##15703244##55##]. As the case for Huai River basin in this study, the skyline plot provided evidence that the exponential growth model could be more realistic than that of stepwise expansion (Figure ##FIG##3##4##). Second, mismatch analysis needed a pair of genes being chosen at random from the population. In populations having gone through a recent and still large expansion, the internal branches should be very short due to the star-like structure of the tree, and a very few mutations would accumulate on those branches [##REF##9335623##40##,##REF##1743491##52##]. In this case, pairs drawn from the sample would be not independent due to the shared portions of their genealogy [##REF##10388826##56##]. In this study, both unrooted ML tree and the statistical parsimony network displayed "star-like" shape with no deep branching among halotypes (Figure ##FIG##1##2##, Figure ##FIG##2##3##). So the inadequate information inferred from mismatch analysis might attribute to the sign in the data set being swamped by non-independence in this study. Moreover, although the unrooted ML tree indicated that there were no distinct haplotype groups with high bootstrap support, the results of our phylogeographic analysis revealed that there were high levels of genetic differentiation among populations, and that the genetic pattern of this species was mainly constituted by multiple unrelated founding events. So cases of founder events not containing all genetic types in some regions might be unavoidable, and could severely affect the shape of the mismatch distribution. Finally, other factors such as high frequency of the ancestral haplotypes, population substructure, or inbreeding could also affect the shape of the mismatch distribution, but to an extent these effects had not been quantified [##REF##10388826##56##].</p> <p>However, the methods of mismatch distribution analysis and neutrality test did not make full use of the data [##UREF##31##57##]. The coalescent-based method, by incorporating information from the genealogical tree structure of DNA sequences, used more information in the data and seemed to be more adequate to estimate demographic growth [##REF##9584114##43##,##UREF##31##57##,##REF##14526072##58##]. In Yangze River and Huai River basin, three different coalescent-based methods, the skyline plot [##REF##11719579##38##,##REF##10880500##39##], the estimate of population growth rate (<italic>g</italic>) using FLUCTUATE 1.4 [##REF##9584114##43##] and the distribution of mutations conducted by GENETREE [##REF##10792974##44##], which provided the similar information of population expansion in these regions, suggested that the superiority of the coalescent-based method. In conclusion, despite of some inconsistency in these results, it was evident that the combined demographic analysis could confirm population expansion in Yangze River and Huai River basin, while a relative constant population size in Pearl River basin.</p> <p>During the last glacial maximum (LGM, about 18 kyr ago), with sea levels lowered (120–140 m below present sea level) [##REF##12226674##12##], the present distribution range of <italic>N. taihuensis </italic>was almost completely exposed and eradicated, and this species should have been severely impacted by the past glaciations. Our estimate of TMRCA in the Pearl River and the Huai River basin, and the population expansion time for <italic>N. taihuensis </italic>in the Yangze River basin were much older than that of LGM. It was difficult to link the population demographic history to any particular Pleistocene paleo-climatic event. However, the gene tree (Figure. ##FIG##5##6##) provided evidence that the Pleistocene ice ages had a great effect on the demographic history of this species. That the vast majority of mutations occurred near the tips of the gene tree suggests that population expansion events mainly occurred after LGM. Therefore, we proposed that in the Yangze River basin and the Huai River basin, the <italic>N. taihuensis </italic>should have been the most severely impacted by the past glaciations, surviving in some refugium, and then post-glacial colonizations and population expansions occurred with the climate warming and sea levels uplifting after the late Quaternary (especially the LGM) [##REF##16503171##59##]. Further, three near-coastal populations (the Taihu lake population and the Chaohu population in the Yangze River basin, the Hongzehu Lake population in the Huai River basin) contained higher numbers of haplotypes, haplotype diversity (<italic>h</italic>) and nucleotide diversity (<italic>π</italic>) (Table ##TAB##0##1##), suggesting that a glacial refugium could exist near the estuaries. However, <italic>N. taihuensis </italic>has been exposed to pollution and habitat destruction for decades; its population size has decreased drastically in out-flowing rivers and estuaries of these basins under study [##UREF##14##24##,##UREF##15##25##]. No sample from these locations could be collected to detect the probable glacial refugium. Further studies should be conducted in future. In addition, it is worthy to note that many attached lakes formed in these areas with sea levels uplifting after the LGM. For example, Chaohu Lake and Taihu Lake in the Yangze River basin shaped about 12 and 2.6 kyr ago, respectively, and Hongzehu lake and Weishanhu lake in the Huai River basin was shaped about 1.4 and 0.6 kyr ago, respectively [##UREF##32##60##]. This indicated that the formation of appropriate habitations facilitated historical population expansions. This inference was compatible with the "r-strategy" characteristic of <italic>N. taihuensis </italic>with high relative fertility and a shorter generation time [##UREF##18##29##], which also suggests its population size could be severely affected by climate and habitation. Thus, it is not surprising that <italic>N. taihuensis </italic>has maintained a relatively constant population size in the Pearl River basin, because this area is a subtropical zone, with average 2 – 3°C higher than in the Yangze River basin or the Huai River basin during the Pleistocene [##UREF##33##61##]. There was no large-scale creation of appropriate habitations (attached lakes) after the LGM in this area.</p> <title>Implications for Conservation</title> <p>The Salangid, <italic>Neosalanx taihuensis</italic>, is one of the most important commercial species in this family. The dry or frozen salangids are good materials in Chinese food, and they are also exported to Japan and other southern-east Asia countries. But, it has to be noticed that the population size of <italic>N. taihuensis </italic>and its fishing yield have decreased rapidly in recent years [##UREF##14##24##,##UREF##15##25##]. For example, the yield of Salangids per year in Poyanghu lake decreased from 6 × 10<sup>5 </sup>kg in 1960s to 10<sup>4 </sup>kg in 1987 [##UREF##34##62##], and the fishery production of <italic>N. taihuensis </italic>was completely collapsed in Poyanghu lake thereafter [##UREF##35##63##]. The yield of Salangids per year in the Yellow River has also decreased since 1987, and the lowest yield per year less than 10<sup>3 </sup>kg occurred in 1990 [##UREF##18##29##]. According to our investigation being from 2004 to 2005 and the information from the local fishery management offices, both adults and juveniles are now rarely found in the Yellow River. Some authors thought that over-fishing was a serious problem [##UREF##36##64##], but You <italic>et al</italic>. suggested that over-fishing was not the most important threat to this species [##UREF##37##65##]. On the other hand, the water pollution should be also the one of the most serious threat to <italic>N. taihuensis </italic>[##UREF##15##25##,##UREF##35##63##,##UREF##38##66##]. The environment of <italic>N. taihuensis </italic>has been polluted because of the waste water from urban, industry and agriculture [##UREF##15##25##]. Finally, the threat to <italic>N. taihuensis </italic>is also from the habitat fragmentation. Since mid 1950s, the aquatic networks have been fragmented dramatically with the construction of dykes and floodgates across the river-lake connections [##UREF##32##60##,##UREF##38##66##]. Most of the lakes in the mid-lower reaches lost their connection with the main river watercourse by the mid 1980s [##UREF##32##60##]. Consequently, most of the small lakes nearly disappeared and the mid-large lakes shrank considerably resulting in unsuitable habitat for this species [##UREF##38##66##,##UREF##39##67##]. Up to the present, no specific protect action has been taken for <italic>N. taihuensis</italic>, so it is urgent to make an effective plan to protect this economically important species.</p> <p>The study on population structure and its phylogeoraphic pattern could provide particularly information for us to more broadly consider conservation options and more accurately assign appropriate fishery management [##REF##16911199##68##]. According to the model proposed by Moritz [##UREF##40##69##], Evolutionary Significant Units (ESUs) are designated on the basis of reciprocal monophyly at mitochondrial markers. In the present study, phylogenetic inference based on ML and NCA methods revealed that it was no geographical structure among haplotypes in <italic>N. taihuensis</italic>, suggesting that there was not a genetically distinct and independent population that could be considered as an ESU. However, AMOVA analysis indicated that significant genetic variance resulted from genetic differences among populations and within basins. In addition, based on low vagility, geographic isolation among most populations and rapidly population decline over the recent decades for this species [##UREF##14##24##,##UREF##15##25##,##UREF##35##63##], we proposed that all populations should be managed and conserved separately. Moreover, special attention must be paid on habitat protection and pollution control because this species is susceptive to appropriate habitat change.</p>
<title>Conclusion</title> <p>Our results showed that this species presented a very random genetic pattern and was not compatible with the present hydrological structure. This genetic pattern might be formed mainly by multiple unrelated long distance founding events, combining with population expansion and the instream barriers restricting gene flow. Demographic analysis revealed that the populations in the Pearl River basin (TMRCA approximately 73.53 kyr) has kept a relatively constant population size of <italic>N. taihuensis</italic>, but populations in the Yangze River basin and the Huai River basin expanded about 221.9 and 190.0 kyr ago, respectively, with the majority of mutations occurring after the LGM. We concluded that the more favourable climate and the formation of appropriate large-scale habitations facilitated this expansion. We suggested that it was necessary to develop a suitable management plan for <italic>N. taihuensis </italic>based on our results.</p>
<title>Background</title> <p>The Salangid icefish <italic>Neosalanx taihuensis </italic>(Salangidae) is an economically important fish, which is endemic to China, restricted to large freshwater systems (e.g. lakes, large rivers and estuaries) and typically exhibit low vagility. The continuous distribution ranges from the temperate region of the Huai and Yellow River basins to the subtropical region of the Pearl River basin. This wide ranging distribution makes the species an ideal model for the study of palaeoclimatic effects on population genetic structure and phylogeography. Here, we aim to analyze population genetic differentiation within and between river basins and demographic history in order to understand how this species responded to severe climatic oscillations, decline of the sea levels during the Pleistocene ice ages and tectonic activity.</p> <title>Results</title> <p>We obtained the complete mtDNA cytochrome <italic>b </italic>sequences (1141 bp) of 354 individuals from 13 populations in the Pearl River, the Yangze River and the Huai River basin. Thirty-six haplotypes were detected. Haplotype frequency distributions were strongly skewed, with most haplotypes (n = 24) represented only in single samples each and thus restricted to a single population. The most common haplotype (H36) was found in 49.15% of all individuals. Analysis of molecular variance (AMOVA) revealed a random pattern in the distribution of genetic diversity, which is inconsistent with contemporary hydrological structure. Significant levels of genetic subdivision were detected among populations within basins rather than between the three basins. Demographic analysis revealed that the population size in the Pearl River basin has remained relatively constant whereas the populations in the Yangze River and the Huai River basins expanded about 221 and 190 kyr ago, respectively, with the majority of mutations occurring after the last glacial maximum (LGM).</p> <title>Conclusion</title> <p>The observed complex genetic pattern of <italic>N. taihuensis </italic>is coherent with a scenario of multiple unrelated founding events by long-distance colonization and dispersal combined with contiguous population expansion and locally restricted gene flow. We also found that this species was likely severely impacted by past glaciations. More favourable climate and the formation of large suitable habitations together facilitated population expansion after the late Quaternary (especially the LGM). We proposed that all populations should be managed and conserved separately, especially for habitat protection.</p>
<title>Abbreviations</title> <p>AMOVA: Analysis of molecular variance; ML: maximum likelihood; NCA: nested clade analysis; IBD: isolation by distance; LGM: last glacial maximum; TMRCA: the most recent common ancestor.</p> <title>Authors' contributions</title> <p>LZ collected samples, isolated genomic DNA, and conducted the amplification and sequencing of the complete Cyt <italic>b </italic>gene. He collected, analyzed and summarized the data, and drafted the manuscript; JZ collected some samples, contributed some sequencing and contributed some to its design; ZL participated in analysis and interpretation of data, and helped draft the manuscript; SMF revised the manuscript and conducted part of analysis; FW and MX revised the draft manuscript and gave their some comments; ML conceived the study and participated in its design and data interpretation, and preparing the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This project was supported by the <bold>N</bold>ational Basic Research Program of China (973 Program: 2007CB411600) and the National Natural Science Foundation of China (No. 30570256). The authors would like to thank Dr. L.F. Zhu, and Dr. B. W. Zhang for their technical guidance and data analysis, L. Yan and Y. L. Hao for their laboratory assistance and useful ideas, and some staffs for their help in collecting samples. Special thanks are given to Prof. Y. Tao for polishing the English, Drs. T. Meng and D. W. Qi for providing geographic information and maps, and some anonymous reviewers for their constructive comments.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p>Map of the distribution of <italic>N. taihuensis </italic>and the sampling locations for this study.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p><bold>Unrooted Maximum-likelihood tree for all 36 haplotypes of <italic>N. taihuensis</italic></bold>. Labels are haplotype identification numbers shown in Table 2. Values indicate bootstrap support for each node based on maximum likelihood inference.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p><bold>Minimum spanning network based on Templeton <italic>et al</italic>. (1987, 1993) statistical parsimony</bold>. Nodes contain the haplotype number and are proportional to the haplotype frequency. White nodes indicate undetected intermediate haplotype states separated by one mutational step. Boxes indicate one-step to two-step nesting levels for the nested clade analysis in NCA.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p><bold>Skyline plots of changes in effective population size on time</bold>. The <italic>x </italic>axis is time of mutations per site before present; the <italic>y </italic>axis is the effective population size equal to <italic>N</italic><sub><italic>e</italic></sub><italic>μ</italic>. Classic skyline was shown as thin line, generalized skyline as thick line, and dot line as the expected demographic history.</p></caption></fig> <fig id="F5" position="float"><label>Figure 5</label><caption><p><bold>Observed and expected mismatch distributions showing the frequencies of pairwise differences</bold>. The observed distributions (bars) are compared for their goodness-of-fit to a Poisson distribution under a model of sudden expansion illustrated by the overlaid curve (black dots and solid lines). <italic>X</italic>-axis: number of pair-wise differences, <italic>Y</italic>-axis: frequency.</p></caption></fig> <fig id="F6" position="float"><label>Figure 6</label><caption><p><bold>Gene tree of the <italic>N. taihuensis </italic>mitochondrial <italic>cytb </italic>gene</bold>. The tree is based on 1,000,000 coalescent simulations. The absolute time scale on the right shows the TMRCA (kyr) in generations using a mutation rate of 1% × 10<sup>-6 </sup>per year. The tree shows the ancestral distribution of mutations and events (TMRCA, recent rapid expansion) in the population history.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Descriptive statistics of <italic>N. taihuensis </italic>phylogroups based on mitochondrial <italic>cytb </italic>sequence data</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Population</th><th align="center">location</th><th align="center"><italic>N</italic></th><th align="center"><italic>H</italic></th><th align="center"><italic>h</italic></th><th align="center"><italic>S</italic></th><th align="center"><italic>k</italic></th><th align="center"><italic>π</italic></th></tr></thead><tbody><tr><td align="center"><bold>Pearl River basin</bold></td><td/><td align="center">44</td><td align="center">5</td><td align="center">0.622 ± 0.064</td><td align="center">8</td><td align="center">1.261</td><td align="center">0.0011 ± 0.0003</td></tr><tr><td align="center">P1</td><td align="center">Luofu River</td><td align="center">31</td><td align="center">5</td><td align="center">0.443 ± 0.105</td><td align="center">8</td><td align="center">1.213</td><td align="center">0.0011 ± 0.0004</td></tr><tr><td align="center">P2</td><td align="center">Pearl River</td><td align="center">13</td><td align="center">3</td><td align="center">0.513 ± 0.144</td><td align="center">2</td><td align="center">0.744</td><td align="center">0.0007 ± 0.0002</td></tr><tr><td align="center"><bold>Yangze River basin</bold></td><td/><td align="center">197</td><td align="center">25</td><td align="center">0.781 ± 0.017</td><td align="center">28</td><td align="center">3.080</td><td align="center">0.0027 ± 0.0001</td></tr><tr><td align="center">Y1</td><td align="center">Taihu lake</td><td align="center">31</td><td align="center">8</td><td align="center">0.798 ± 0.039</td><td align="center">11</td><td align="center">3.174</td><td align="center">0.0028 ± 0.0002</td></tr><tr><td align="center">Y2</td><td align="center">Chaohu lake</td><td align="center">25</td><td align="center">8</td><td align="center">0.763 ± 0.068</td><td align="center">10</td><td align="center">1.380</td><td align="center">0.0012 ± 0.0003</td></tr><tr><td align="center">Y3</td><td align="center">Bohu lake</td><td align="center">29</td><td align="center">5</td><td align="center">0.650 ± 0.078</td><td align="center">10</td><td align="center">2.606</td><td align="center">0.0023 ± 0.0004</td></tr><tr><td align="center">Y4</td><td align="center">Poyanghu lake</td><td align="center">22</td><td align="center">5</td><td align="center">0.407 ± 0.128</td><td align="center">8</td><td align="center">1.589</td><td align="center">0.0014 ± 0.0005</td></tr><tr><td align="center">Y5</td><td align="center">Dongtinghu lake</td><td align="center">27</td><td align="center">5</td><td align="center">0.393 ± 0.114</td><td align="center">6</td><td align="center">0.576</td><td align="center">0.0005 ± 0.0002</td></tr><tr><td align="center">Y6</td><td align="center">Tianrezhou</td><td align="center">34</td><td align="center">5</td><td align="center">0.570 ± 0.078</td><td align="center">8</td><td align="center">1.210</td><td align="center">0.0011 ± 0.0003</td></tr><tr><td align="center">Y7</td><td align="center">Xujiahe Reservoir</td><td align="center">29</td><td align="center">2</td><td align="center">0.296 ± 0.093</td><td align="center">3</td><td align="center">0.887</td><td align="center">0.0008 ± 0.0002</td></tr><tr><td align="center"><bold>Huai River basin</bold></td><td/><td align="center">113</td><td align="center">15</td><td align="center">0.464 ± 0.057</td><td align="center">18</td><td align="center">1.110</td><td align="center">0.0010 ± 0.0002</td></tr><tr><td align="center">H1</td><td align="center">Hongzehu lake</td><td align="center">26</td><td align="center">8</td><td align="center">0.686 ± 0.088</td><td align="center">11</td><td align="center">2.348</td><td align="center">0.0021 ± 0.0003</td></tr><tr><td align="center">H2</td><td align="center">Weishanhu lake</td><td align="center">31</td><td align="center">4</td><td align="center">0.553 ± 0.060</td><td align="center">6</td><td align="center">0.796</td><td align="center">0.0007 ± 0.0002</td></tr><tr><td align="center">H3</td><td align="center">Wabuhu lake</td><td align="center">30</td><td align="center">2</td><td align="center">0.067 ± 0.061</td><td align="center">1</td><td align="center">0.067</td><td align="center">0.0001 ± 0.0001</td></tr><tr><td align="center">H4</td><td align="center">Chengdonghu lake</td><td align="center">26</td><td align="center">6</td><td align="center">0.354 ± 0.119</td><td align="center">9</td><td align="center">0.834</td><td align="center">0.0007 ± 0.0003</td></tr><tr><td align="center"><bold>Entire region (all samples)</bold></td><td/><td align="center">354</td><td align="center">36</td><td align="center">0.713 ± 0.022</td><td align="center">38</td><td align="center">2.528</td><td align="center">0.0022 ± 0.0001</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Summary of mtDNA <italic>cytb </italic>region haplotype distributions</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Haplotypes\populations</th><th align="center" colspan="2">Pearl River basin</th><th align="center" colspan="7">Yangze River basin</th><th align="center" colspan="4">Huai River basin</th><th align="center">Total</th></tr><tr><th/><th colspan="13"><hr/></th><th/></tr><tr><th/><th align="center">P1</th><th align="center">P2</th><th align="center">Y1</th><th align="center">Y2</th><th align="center">Y3</th><th align="center">Y4</th><th align="center">Y5</th><th align="center">Y6</th><th align="center">Y7</th><th align="center">H1</th><th align="center">H2</th><th align="center">H3</th><th align="center">H4</th><th/></tr></thead><tbody><tr><td align="center"><bold>H01</bold></td><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H02</bold></td><td align="center">3</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">3</td></tr><tr><td align="center"><bold>H03</bold></td><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H04</bold></td><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H05</bold></td><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H06</bold></td><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H07</bold></td><td/><td/><td align="center">2</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">2</td></tr><tr><td align="center"><bold>H08</bold></td><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H09</bold></td><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H10</bold></td><td/><td/><td/><td align="center">4</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">4</td></tr><tr><td align="center"><bold>H11</bold></td><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H12</bold></td><td/><td/><td/><td/><td align="center">2</td><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">2</td></tr><tr><td align="center"><bold>H13</bold></td><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H14</bold></td><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H15</bold></td><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H16</bold></td><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H17</bold></td><td/><td/><td/><td/><td/><td/><td align="center">3</td><td/><td/><td/><td/><td/><td/><td align="center">3</td></tr><tr><td align="center"><bold>H18</bold></td><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H19</bold></td><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H20</bold></td><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H21</bold></td><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H22</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H23</bold></td><td/><td/><td/><td align="center">1</td><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td align="center">2</td></tr><tr><td align="center"><bold>H24</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H25</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">2</td><td/><td/><td/><td align="center">2</td></tr><tr><td align="center"><bold>H26</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H27</bold></td><td/><td/><td align="center">9</td><td/><td align="center">16</td><td/><td/><td align="center">3</td><td align="center">24</td><td align="center">5</td><td align="center">1</td><td/><td/><td align="center">58</td></tr><tr><td align="center"><bold>H28</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td/><td align="center">1</td></tr><tr><td align="center"><bold>H29</bold></td><td align="center">1</td><td align="center">9</td><td/><td/><td/><td align="center">2</td><td/><td align="center">8</td><td/><td/><td align="center">11</td><td/><td/><td align="center">31</td></tr><tr><td align="center"><bold>H30</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td/><td align="center">1</td></tr><tr><td align="center"><bold>H31</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>H32</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>H33</bold></td><td/><td/><td align="center">7</td><td align="center">1</td><td align="center">6</td><td align="center">17</td><td/><td/><td align="center">5</td><td align="center">1</td><td/><td/><td align="center">1</td><td align="center">38</td></tr><tr><td align="center"><bold>H34</bold></td><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>H35</bold></td><td align="center">3</td><td align="center">2</td><td/><td align="center">5</td><td/><td/><td/><td/><td/><td/><td/><td/><td align="center">1</td><td align="center">11</td></tr><tr><td align="center"><bold>H36</bold></td><td align="center">23</td><td align="center">2</td><td align="center">9</td><td align="center">11</td><td align="center">4</td><td align="center">1</td><td align="center">21</td><td align="center">21</td><td/><td align="center">14</td><td align="center">18</td><td align="center">29</td><td align="center">21</td><td align="center">174</td></tr><tr><td align="center"><bold>Sample Number</bold></td><td align="center">31</td><td align="center">13</td><td align="center">31</td><td align="center">25</td><td align="center">29</td><td align="center">22</td><td align="center">27</td><td align="center">34</td><td align="center">29</td><td align="center">26</td><td align="center">31</td><td align="center">30</td><td align="center">26</td><td align="center">354</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Pairwise <italic>Φ</italic><sub><italic>ST </italic></sub>among <italic>N. taihuensis </italic>populations</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">populations</th><th align="center" colspan="2">Pearl River basin</th><th align="center" colspan="7">Yangze River basin</th><th align="center" colspan="4">Huai River basin</th></tr><tr><th/><th colspan="13"><hr/></th></tr><tr><th/><th align="center">P<sub>1</sub></th><th align="center">P<sub>2</sub></th><th align="center">Y<sub>1</sub></th><th align="center">Y<sub>2</sub></th><th align="center">Y<sub>3</sub></th><th align="center">Y<sub>4</sub></th><th align="center">Y<sub>5</sub></th><th align="center">Y<sub>6</sub></th><th align="center">Y<sub>7</sub></th><th align="center">H<sub>1</sub></th><th align="center">H<sub>2</sub></th><th align="center">H<sub>3</sub></th><th align="center">H<sub>4</sub></th></tr></thead><tbody><tr><td align="center">P<sub>1</sub></td><td/><td align="center">61.9</td><td align="center">1335.0</td><td align="center">1155.4</td><td align="center">953.7</td><td align="center">808.1</td><td align="center">808.5</td><td align="center">869.7</td><td align="center">1081.9</td><td align="center">1429.4</td><td align="center">1532.0</td><td align="center">1221.0</td><td align="center">1193.8</td></tr><tr><td align="center">P<sub>2</sub></td><td align="center">0.3146*</td><td/><td align="center">1360.6</td><td align="center">1176.6</td><td align="center">974.9</td><td align="center">829.3</td><td align="center">777.5</td><td align="center">838.7</td><td align="center">1074.7</td><td align="center">1450.6</td><td align="center">1553.2</td><td align="center">1215.0</td><td align="center">1242.2</td></tr><tr><td align="center">Y<sub>1</sub></td><td align="center">0.3081*</td><td align="center">0.4161*</td><td/><td align="center">331.0</td><td align="center">471.7</td><td align="center">619.1</td><td align="center">983.8</td><td align="center">822.0</td><td align="center">658.2</td><td align="center">278.9</td><td align="center">604.6</td><td align="center">380.5</td><td align="center">423.6</td></tr><tr><td align="center">Y<sub>2</sub></td><td align="center">0.0309</td><td align="center">0.2873*</td><td align="center">0.3292*</td><td/><td align="center">194.9</td><td align="center">344.0</td><td align="center">645.9</td><td align="center">484.1</td><td align="center">328.8</td><td align="center">261.2</td><td align="center">462.1</td><td align="center">108.8</td><td align="center">125.1</td></tr><tr><td align="center">Y<sub>3</sub></td><td align="center">0.5349*</td><td align="center">0.6046*</td><td align="center">0.0511</td><td align="center">0.5442*</td><td/><td align="center">142.3</td><td align="center">513.0</td><td align="center">351.1</td><td align="center">333.4</td><td align="center">470.5</td><td align="center">600.8</td><td align="center">268.1</td><td align="center">259.2</td></tr><tr><td align="center">Y<sub>4</sub></td><td align="center">0.7094*</td><td align="center">0.7592*</td><td align="center">0.2791*</td><td align="center">0.6947*</td><td align="center">0.2818*</td><td/><td align="center">413.4</td><td align="center">291.8</td><td align="center">389.6</td><td align="center">618.0</td><td align="center">720.6</td><td align="center">409.6</td><td align="center">382.4</td></tr><tr><td align="center">Y<sub>5</sub></td><td align="center">0.0420</td><td align="center">0.4304*</td><td align="center">0.4070*</td><td align="center">0.0633*</td><td align="center">0.6241*</td><td align="center">0.7834*</td><td/><td align="center">156.4</td><td align="center">419.6</td><td align="center">885.8</td><td align="center">951.7</td><td align="center">658.0</td><td align="center">626.7</td></tr><tr><td align="center">Y<sub>6</sub></td><td align="center">0.0340</td><td align="center">0.1611</td><td align="center">0.3255*</td><td align="center">0.0817*</td><td align="center">0.5451*</td><td align="center">0.7073*</td><td align="center">0.0768*</td><td/><td align="center">271.5</td><td align="center">723.9</td><td align="center">789.8</td><td align="center">496.1</td><td align="center">464.8</td></tr><tr><td align="center">Y<sub>7</sub></td><td align="center">0.7476*</td><td align="center">0.8327*</td><td align="center">0.2340*</td><td align="center">0.7560*</td><td align="center">0.0675</td><td align="center">0.5831*</td><td align="center">0.8353*</td><td align="center">0.7463*</td><td/><td align="center">560.0</td><td align="center">500.8</td><td align="center">333.8</td><td align="center">276.9</td></tr><tr><td align="center">H<sub>1</sub></td><td align="center">0.0753</td><td align="center">0.2961</td><td align="center">0.0996</td><td align="center">0.1331</td><td align="center">0.3043*</td><td align="center">0.5349*</td><td align="center">0.1874*</td><td align="center">0.1016</td><td align="center">0.5247*</td><td/><td align="center">316.1</td><td align="center">224.9</td><td align="center">291.3</td></tr><tr><td align="center">H<sub>2</sub></td><td align="center">0.1196*</td><td align="center">0.1137</td><td align="center">0.4003*</td><td align="center">0.1459*</td><td align="center">0.6103*</td><td align="center">0.7622*</td><td align="center">0.1541*</td><td align="center">-0.0008</td><td align="center">0.8109*</td><td align="center">0.1988*</td><td/><td align="center">351.8</td><td align="center">338.1</td></tr><tr><td align="center">H<sub>3</sub></td><td align="center">0.0582</td><td align="center">0.6489*</td><td align="center">0.4499*</td><td align="center">0.1012*</td><td align="center">0.6725*</td><td align="center">0.8403*</td><td align="center">0.0300</td><td align="center">0.0927*</td><td align="center">0.8882*</td><td align="center">0.2283*</td><td align="center">0.2123*</td><td/><td align="center">53.6</td></tr><tr><td align="center">H<sub>4</sub></td><td align="center">-0.0008</td><td align="center">0.3672*</td><td align="center">0.3368*</td><td align="center">0.0295</td><td align="center">0.5673*</td><td align="center">0.7339*</td><td align="center">0.0169</td><td align="center">0.0383</td><td align="center">0.7941*</td><td align="center">0.1157</td><td align="center">0.1212*</td><td align="center">0.0190</td><td/></tr></tbody></table></table-wrap> <table-wrap id="T4" position="float"><label>Table 4</label><caption><p>The results of AMOVA for <italic>N. taihuensis </italic>mtDNA <italic>cytb </italic>estimated using <italic>Φ</italic>-statistics</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="center">Source of variation</th><th align="center">d.f.</th><th align="center">Sum of squares</th><th align="center">Variance components</th><th align="center">% variation</th><th align="center">Φ statistics (P-value)</th></tr></thead><tbody><tr><td align="center"><bold>Among basins</bold></td><td align="center">2</td><td align="center">55.055</td><td align="center">0.1143</td><td align="center">8.43</td><td align="center">Φ<sub>CT </sub>= 0.0843</td></tr><tr><td align="center"><bold>Among populations within basins</bold></td><td align="center">10</td><td align="center">158.543</td><td align="center">0.5593</td><td align="center">41.26</td><td align="center">Φ<sub>SC </sub>= 0.4506**</td></tr><tr><td align="center"><bold>Within populations</bold></td><td align="center">341</td><td align="center">232.530</td><td align="center">0.6819</td><td align="center">50.31</td><td align="center">Φ<sub>ST </sub>= 0.4970**</td></tr><tr><td align="center"><bold>Total</bold></td><td align="center">353</td><td align="center">446.127</td><td align="center">1.3555</td><td/><td/></tr></tbody></table></table-wrap> <table-wrap id="T5" position="float"><label>Table 5</label><caption><p>Statistical test for neutrality, mismatch analysis and the estimate of demographic parameters for <italic>N. taihuensis </italic>based on mitochondrial <italic>cytb </italic>sequence data</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th align="center">Entire region</th><th align="center">Pearl River basin</th><th align="center">Yangze River basin</th><th align="center">Huai River basin</th></tr></thead><tbody><tr><td align="center"><italic>Hri </italic>(<bold><italic>P </italic>value</bold>)</td><td align="center">0.0675(0.4687)</td><td align="center">0.1418(0.0608)</td><td align="center">0.0518(0.2609)</td><td align="center">0.1446(0.999)</td></tr><tr><td align="center"><italic>SSD </italic>(<bold><italic>P </italic>value</bold>)</td><td align="center">0.0388 (0.3002)</td><td align="center">0.0239(0.0411) *</td><td align="center">0.0336(0.1401)</td><td align="center">0.0131(0.001)**</td></tr><tr><td align="center"><italic>D</italic>* (<bold><italic>P </italic>value</bold>)</td><td align="center">-0.1047(0.5730)</td><td align="center">-0.1036(0.5910)</td><td align="center">-0.0991(0.5590)</td><td align="center">-0.1276(0.5060)</td></tr><tr><td align="center"><italic>R</italic><sub>2 </sub>(<bold><italic>P </italic>value</bold>)</td><td align="center">0.0726(0.0160)*</td><td align="center">0.1177(0.1794)</td><td align="center">0.0825(0.1430)</td><td align="center">0.0859(0.0020)**</td></tr><tr><td align="center"><italic>Fs </italic>(<bold><italic>P </italic>value</bold>)</td><td align="center">-18.9472 (0.0002)**</td><td align="center">0.4106(0.6238)</td><td align="center">-7.4204 (0.0279) *</td><td align="center">-9.2412 (0.0019) **</td></tr><tr><td align="center"><italic>g </italic>± 3*SD</td><td align="center">2648.97 ± 556.26</td><td align="center">-309.64 ± 955.56</td><td align="center">3268.21 ± 873.30</td><td align="center">2083.14 ± 529.35</td></tr><tr><td align="center"><bold><italic>τ </italic></bold>(<bold>90%CI</bold>)</td><td align="center">5.602 (0.619–10.090)</td><td/><td align="center">5.064 (1.334–8.439)</td><td/></tr><tr><td align="center">TMRCA (<bold>kyr</bold>)</td><td/><td align="center">73.53</td><td/><td align="center">190.0</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Nested clade analysis. Statistical analysis of the current and historic patterns of phylogenetic and geographic associations.</p></caption></supplementary-material>
<table-wrap-foot><p>Note:<italic>N</italic>, number of individuals; <italic>H</italic>, number of haplotypes; <italic>S</italic>, number of segregating sites; <italic>h</italic>, gene diversity (± Standard Deviation); <italic>k</italic>, mean pair-wise nucleotide; <italic>π</italic>, nucleotide diversity (± Standard Deviation)</p></table-wrap-foot> <table-wrap-foot><p><italic>Φ</italic><sub><italic>ST </italic></sub>values are reported below, and the river distances (km) between populations above, the diagonal. Asterisks indicate significant values after Bonferroni correction.</p></table-wrap-foot> <table-wrap-foot><p>(*<italic>P </italic>&lt; 0.05, ** <italic>P </italic>&lt; 0.01)</p></table-wrap-foot> <table-wrap-foot><p>(*<italic>P </italic>&lt; 0.05, ** <italic>P </italic>&lt; 0.01)</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-201-1"/> <graphic xlink:href="1471-2148-8-201-2"/> <graphic xlink:href="1471-2148-8-201-3"/> <graphic xlink:href="1471-2148-8-201-4"/> <graphic xlink:href="1471-2148-8-201-5"/> <graphic xlink:href="1471-2148-8-201-6"/>
<media xlink:href="1471-2148-8-201-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Avise"], "given-names": ["JC"], "source": ["Phylogeography: The History and Formation of Species"], "year": ["2000"], "publisher-name": ["USA: Harvard University Press"]}, {"surname": ["Palumbi"], "given-names": ["SR"], "article-title": ["Genetic divergence, reproductive isolation, and marine speciation"], "source": ["Annu Rev Ecol Syst"], "year": ["1994"], "volume": ["25"], "fpage": ["547"], "lpage": ["572"], "pub-id": ["10.1146/annurev.es.25.110194.002555"]}, {"surname": ["Kamal", "Richard", "Hewitt"], "given-names": ["MI", "AN", "GM"], "article-title": ["Spatial patterns of genetic variation generated by different forms of dispersal during range expansion"], "source": ["Heredity"], "year": ["1996"], "volume": ["77"], "fpage": ["282"], "lpage": ["291"], "pub-id": ["10.1038/hdy.1996.142"]}, {"surname": ["Ward", "Woodwark", "Skibinski"], "given-names": ["R", "M", "O"], "article-title": ["A comparison of genetic diversity levels in marine, freshwater and anadromous fishes"], "source": ["Journal of Fish Biology"], "year": ["1994"], "volume": ["44"], "fpage": ["213"], "lpage": ["232"], "pub-id": ["10.1111/j.1095-8649.1994.tb01200.x"]}, {"surname": ["Gyllensten"], "given-names": ["U"], "article-title": ["The genetic structure of fish: differences in the intraspecific distribution of biochemical genetic variation between marine, anadromous, and freshwater species"], "source": ["Journal of Fish Biology"], "year": ["1985"], "volume": ["26"], "fpage": ["691"], "lpage": ["699"], "pub-id": ["10.1111/j.1095-8649.1985.tb04309.x"]}, {"surname": ["Hansen", "Mensberg"], "given-names": ["MM", "K-LD"], "article-title": ["Genetic differentiation and relationship between genetic and geographical distance in Danish sea trout (Salmo trutta L.) populations"], "source": ["Heredity"], "year": ["1998"], "volume": ["81"], "fpage": ["493"], "lpage": ["508"], "pub-id": ["10.1046/j.1365-2540.1998.00408.x"]}, {"surname": ["Xia", "Chen", "Sheng"], "given-names": ["YZ", "YY", "Y"], "article-title": ["Phylogeographic Structure of Lenok ("], "italic": ["Brachymystax lenok "], "source": ["Zoological Studies"], "year": ["2006"], "volume": ["45"], "fpage": ["190"], "lpage": ["200"]}, {"surname": ["Currens", "Schreck", "Li"], "given-names": ["K", "C", "WH"], "article-title": ["Allozyme and morphological divergence of rainbow trout (Oncorhynchus mykiss) above and below waterfalls in the Deschutes River, Oregon"], "source": ["Copeia"], "year": ["1990"], "volume": ["1990"], "fpage": ["730"], "lpage": ["746"], "pub-id": ["10.2307/1446439"]}, {"surname": ["Avise", "Arnold", "Ball", "Bermingham", "Lamb"], "given-names": ["JC", "J", "RM", "E", "T"], "article-title": ["Intraspecific Phylogeography: The Mitochondrial DNA Bridge Between Population Genetics and Systematics"], "source": ["Annual Review of Ecology and Systematics"], "year": ["1987"], "volume": ["18"], "fpage": ["489"], "lpage": ["522"]}, {"surname": ["Zhang", "Qiao"], "given-names": ["YL", "XG"], "article-title": ["Study on phylogeny and zoogeography of fishes of the family Salangidae"], "source": ["Acta Zoologica Taiwanica"], "year": ["1994"], "volume": ["5"], "fpage": ["95"], "lpage": ["115"]}, {"surname": ["Nelson"], "given-names": ["J"], "source": ["Fishes of the World"], "year": ["1994"], "edition": ["3"], "publisher-name": ["New York: John Wiley & Sons, Inc."]}, {"surname": ["Cheng", "Zheng"], "given-names": ["QT", "BS"], "source": ["Systematic Synopsis of Chinese Fishes"], "year": ["1987"], "publisher-name": ["Beijing: Science Press"]}, {"surname": ["Roberts"], "given-names": ["T"], "article-title": ["Skeletal anatomy and classification of the neotenic Asian salmoniform superfamily Salangoidea (icefishes or noodlefishes)"], "source": ["Proceedings of the California Academy of Science"], "year": ["1984"], "volume": ["43"], "fpage": ["179"], "lpage": ["220"]}, {"surname": ["Bernatchez", "Wilson"], "given-names": ["L", "CC"], "article-title": ["Comparative phylogeography of Nearctic and Palearctic fishes"], "source": ["Molecular Ecology"], "year": ["1998"], "volume": ["7"]}, {"surname": ["Wang", "Lu", "Hu", "Zhou", "Xu", "Lei"], "given-names": ["ZS", "C", "HJ", "Y", "CR", "GC"], "article-title": ["Freshwater icefishes (Salangidae) in the Yangtze River basin of China: Spatial distribution patterns and environmental determinants"], "source": ["Environmental Biology of Fishes"], "year": ["2005"], "volume": ["73"], "fpage": ["253"], "pub-id": ["10.1007/s10641-005-2146-3"]}, {"surname": ["Li"], "given-names": ["SF"], "source": ["A study on biodiversity and its conservation of major fishes in the Yangtze River"], "year": ["2001"], "publisher-name": ["Shanghai: Shanghai Scientific and Technical Publishers"]}, {"surname": ["Regan"], "given-names": ["C"], "article-title": ["Description of new fishes from Lake Candidius, Formosa, collected by Dr. A. Moltrecht"], "source": ["Annals and Magazine of Natural History Series"], "year": ["1908"], "volume": ["8"], "fpage": ["358"], "lpage": ["360"]}, {"surname": ["Wakiya", "Takahasi"], "given-names": ["Y", "N"], "article-title": ["Study on fishes of the family Salangidae"], "source": ["Journal of the College of Agriculture, Imperial University of Tokyo"], "year": ["1937"], "volume": ["14"], "fpage": ["265"], "lpage": ["296"]}, {"surname": ["Dou", "Chen"], "given-names": ["SZ", "DG"], "article-title": ["Taxonomy, biology and abundance of Icefishes Salangidae in the Yellow River estuary of the Bohai Sea, China"], "source": ["Journal of Fish Biology"], "year": ["1994"], "volume": ["45"], "fpage": ["737"], "lpage": ["748"], "pub-id": ["10.1111/j.1095-8649.1994.tb00940.x"]}, {"surname": ["Zhang", "Li", "Xu", "Takita", "Wei"], "given-names": ["J", "M", "MQ", "T", "FW"], "article-title": ["Molecular phylogeny of icefish Salangidae based on complete mtDNA cytochrome b sequences, with comments on estuarine fish evolution"], "source": ["Biological Journal of the Linnean Society"], "year": ["2007"], "volume": ["91"], "fpage": ["325"], "lpage": ["340"], "pub-id": ["10.1111/j.1095-8312.2007.00785.x"]}, {"surname": ["Kutzbach", "Ruddiman", "Prell"], "given-names": ["J", "W", "W"], "person-group": ["WF R"], "article-title": ["Tectonic Uplift and Climate Change"], "source": ["New York"], "year": ["1997"], "fpage": ["149"], "lpage": ["170"]}, {"surname": ["Luo", "Xing", "Wang", "Shen", "Ren", "Zeng", "Zhu"], "given-names": ["LX", "JM", "NQ", "YC", "ME", "ZX", "SD"], "source": ["The Physical Geography of China"], "year": ["1981"], "publisher-name": ["Beijing: Science Press"]}, {"surname": ["Huang", "Su", "Mei"], "given-names": ["xQ", "FC", "HA"], "source": ["Rivers in China"], "year": ["1995"], "publisher-name": ["Beijing: Commercial Press"]}, {"surname": ["Myers"], "given-names": ["GS"], "article-title": ["Salt-tolerance of fresh-water fish groups in relation to zoogeographical problems"], "source": ["Bijdr Dierk"], "year": ["1949"], "volume": ["28"], "fpage": ["315"], "lpage": ["322"]}, {"surname": ["Zhang", "Hewitt"], "given-names": ["D", "G"], "article-title": ["Nuclear integrations challenges for mitochondrial DNA markers"], "source": ["Trends of Ecology and Evolution"], "year": ["1996"], "volume": ["11"], "fpage": ["247"], "lpage": ["251"], "pub-id": ["10.1016/0169-5347(96)10031-8"]}, {"surname": ["Cook", "Sinclair", "Stefansson"], "given-names": ["RM", "A", "G"], "article-title": ["Potential collapse of North Sea cod stocks"], "source": ["Nature"], "year": ["1997"], "volume": ["385"], "fpage": ["521"], "lpage": ["522"], "pub-id": ["10.1038/385521a0"]}, {"surname": ["Grant", "Bowen"], "given-names": ["WAS", "BW"], "article-title": ["Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation"], "source": ["J Hered"], "year": ["1998"], "volume": ["89"], "fpage": ["415"], "lpage": ["426"], "pub-id": ["10.1093/jhered/89.5.415"]}, {"surname": ["Frankham"], "given-names": ["R"], "article-title": ["Relationship of genetic variation to population size in wildlife"], "source": ["Conservation Biology"], "year": ["1996"], "volume": ["10"], "fpage": ["1500"], "lpage": ["1508"], "pub-id": ["10.1046/j.1523-1739.1996.10061500.x"]}, {"surname": ["Meffe", "Vrijenhoek"], "given-names": ["G", "R"], "article-title": ["Conservation genetics in the management of desert fishes"], "source": ["Conservation Biology"], "year": ["1988"], "volume": ["2"], "fpage": ["157"], "lpage": ["169"], "pub-id": ["10.1111/j.1523-1739.1988.tb00167.x"]}, {"surname": ["Hynes", "Ferguson", "McCann"], "given-names": ["R", "A", "M"], "article-title": ["Variation in mitochondrial DNA and post-glacial colonisation of north-west Europe by. brown trout (Salmo trutta L.)"], "source": ["Journal of Fish Biology"], "year": ["1996"], "volume": ["48"], "fpage": ["54"], "lpage": ["67"]}, {"surname": ["Polanski", "Kimmel", "Chakraborty"], "given-names": ["A", "M", "R"], "article-title": ["Application of a time-dependent coalescence process for inferring the history of population size changes from DNA sequence data"], "source": ["Proceedings of the National Academy of Sciences"], "year": ["1998"], "volume": ["95"], "fpage": ["5456"], "lpage": ["5461"], "pub-id": ["10.1073/pnas.95.10.5456"]}, {"surname": ["Emerson", "Paradis", "Thebaud"], "given-names": ["BC", "E", "C"], "article-title": ["Revealing the demographic histories of species using DNA sequences"], "source": ["Trends in Ecology & Evolution"], "year": ["2001"], "volume": ["16"], "fpage": ["707"], "pub-id": ["10.1016/S0169-5347(01)02305-9"]}, {"surname": ["Wang", "Dou"], "given-names": ["SM", "HS"], "source": ["Lakes in China"], "year": ["1998"], "publisher-name": ["Beijing: Science Press"]}, {"surname": ["Xia"], "given-names": ["ZK"], "source": ["The Pleistocene Envioment"], "year": ["1997"], "publisher-name": ["Beijing: Beijing University Press"]}, {"surname": ["Chen", "Zhang"], "given-names": ["GH", "B"], "article-title": ["study on the growth of salangid fishes in Poyang lake"], "source": ["JiangXi science"], "year": ["1991"], "volume": ["9"], "fpage": ["225"], "lpage": ["232"]}, {"surname": ["Wang", "Fu", "Lui"], "given-names": ["ZS", "CZ", "GC"], "article-title": ["Biodiversity of Chinese icefishes (Salangidae) and their consering strategies"], "source": ["Biodiversity Science"], "year": ["2002"], "volume": ["10"], "fpage": ["416"], "lpage": ["424"]}, {"surname": ["Wang", "Lu", "Hu", "Xu", "Lei"], "given-names": ["ZS", "C", "HJ", "CR", "GC"], "article-title": ["Dynamics of icefish (Salangidae) stocks in Nanyi Lake, eastern China: Degradation and overfishing"], "source": ["J Freshwater Ecol"], "year": ["2004"], "volume": ["19"], "fpage": ["271"], "lpage": ["278"]}, {"surname": ["You", "You"], "given-names": ["Y", "Q"], "article-title": ["A study on bioecological and economic effect of "], "italic": ["Neosalanx Taihuensis "], "source": ["Journal of Southwest Nationalities College"], "year": ["1999"], "volume": ["25"], "fpage": ["269"], "lpage": ["273"]}, {"surname": ["Zeng"], "given-names": ["XC"], "source": ["Fishery resources in the Yangtze River basin"], "year": ["1990"], "publisher-name": ["Beijing: Marine Press"]}, {"surname": ["Xie", "Chen"], "given-names": ["P", "Y"], "article-title": ["Threats to biodiversity in Chinese inland waters"], "source": ["Ambio"], "year": ["1999"], "volume": ["28"], "fpage": ["674"], "lpage": ["681"]}, {"surname": ["Moritz"], "given-names": ["C"], "article-title": ["Applications of mitochondrial DNA analysis in conservation: a critical review"], "source": ["Molecular Ecology"], "year": ["1994"], "volume": ["3"], "fpage": ["401"], "lpage": ["411"], "pub-id": ["10.1111/j.1365-294X.1994.tb00080.x"]}, {"surname": ["Sheldon", "Robert", "Michael"], "given-names": ["JM", "HD", "JS"], "article-title": ["Phylogeny of Pacific salmon and trout based on growth hormone type-2 and mitochondrial NADH dehydrogenase subunit 3 DNA sequences"], "source": ["Can J Fish Aquat Sci"], "year": ["1996"], "volume": ["53"], "fpage": ["1165"], "lpage": ["1176"], "pub-id": ["10.1139/cjfas-53-5-1165"]}, {"surname": ["Asahida", "Kobayashi", "Saito", "Nakayama"], "given-names": ["T", "T", "K", "I"], "article-title": ["Tissue. preservation and total DNA extraction from fish stored at ambient temperature using buffers containing high concentration of urea"], "source": ["Fisheries Science"], "year": ["1996"], "volume": ["62"], "fpage": ["727"], "lpage": ["730"]}, {"surname": ["Nei"], "given-names": ["M"], "source": ["Molecular Evolutionary Genetics"], "year": ["1987"], "publisher-name": ["New York: Columbia University Press"]}, {"surname": ["Weir", "Cockerham"], "given-names": ["B", "C"], "article-title": ["Estimating F-statistic for the analysis of population structure"], "source": ["Evolution"], "year": ["1984"], "volume": ["38"], "fpage": ["1358"], "lpage": ["1370"], "pub-id": ["10.2307/2408641"]}, {"surname": ["Rice"], "given-names": ["WR"], "article-title": ["Analyzing tables of statistical tests"], "source": ["Evolution"], "year": ["1989"], "volume": ["43"], "fpage": ["223"], "lpage": ["225"], "pub-id": ["10.2307/2409177"]}, {"surname": ["Swofford"], "given-names": ["D"], "source": ["paup*: Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4.0b10"], "year": ["2002"], "publisher-name": ["Massachusetts: Sinauer As sociates, Sunderland"]}, {"surname": ["Petit"], "given-names": ["RJ"], "article-title": ["On the falsifiability of the nested clade phylogeographic analysis method"], "source": ["Molecular Ecology"], "year": ["2008"], "volume": ["17"], "fpage": ["1404"], "lpage": ["1404"], "pub-id": ["10.1111/j.1365-294X.2008.03692.x"]}, {"surname": ["Ihaka", "Gentleman"], "given-names": ["R", "R"], "article-title": ["R: A Language for Data Analysis and Graphics"], "source": ["Journal of Computational and Graphical Statistics"], "year": ["1996"], "volume": ["5"], "fpage": ["299"], "lpage": ["314"], "pub-id": ["10.2307/1390807"]}, {"surname": ["Akaike"], "given-names": ["H"], "article-title": ["A new look at the statistical model identification"], "source": ["IEEE Transactions on Automatic Control"], "year": ["1974"], "volume": ["AC-19"], "fpage": ["716"], "lpage": ["723"], "pub-id": ["10.1109/TAC.1974.1100705"]}, {"surname": ["Pond", "Frost", "Muse"], "given-names": ["SLK", "SDW", "SV"], "article-title": ["HyPhy: hypothesis testing using phylogenies"], "year": ["2005"], "volume": ["21"], "fpage": ["676"], "lpage": ["679"]}, {"surname": ["Smith"], "given-names": ["GR"], "article-title": ["Introgression in Fishes: Significance for paleontology, cladistics, and evolutionary rates"], "source": ["Systematic Biology"], "year": ["1992"], "volume": ["41"], "fpage": ["41"], "lpage": ["57"], "pub-id": ["10.2307/2992505"]}]
BMC Evol Biol. 2008 Jul 14; 8:201