Multiple factors have shaped the phylogeography of Chinese spiny loach Cobitis sinensis in Taiwan as inferred from mitochondrial DNA variation

Population genetic diversity and structure of Rhinogobius candidianus (Gobiidae) in Taiwan: Translocation and release

Rhinogobius candidianus is a freshwater goby distributed in north, northwest, west, and south Taiwan, but this species has been introduced to east Taiwan and became dominant. To investigate its native population genetic diversity and structure and evaluate the sources and diversity of translocated populations, the mitochondrial DNA control region and cytochrome b gene (1981 bp) from 220 specimens were analyzed. These results indicated that (1) the east populations originated from two sources in west Taiwan; (2) translocated populations exist in east Taiwan and south Taiwan; (3) many populations have likely been moved secondarily by human intervention; (4) the effective size of the populations had declined greatly; (5) within the native populations, the ancestral populations colonized Taiwan during the land bridge phase in the Pleistocene through north Taiwan; (6) the landform changes in Taiwan shaped the population structure; and (7) the landforms of the coastline during glaciation also shaped the native range. The low-level genetic diversity, high population differentiation, and population decline greatly suggest the need for resource management and conservation interventions. Four clades (α-δ) should be managed as four distinct evolutionarily significant units, while the translocated populations should be managed as separate management units. Moreover, the translocated populations in east Taiwan should be evaluated and monitored carefully.

Historical connectivity and environmental filtering jointly determine the freshwater fish assemblages on Taiwan and Hainan Islands of China

The biotas of Taiwan and Hainan Islands are of continental origin, but the manner with which historical and ecological factors shaped these insular species is still unclear. Here, we used freshwater fish as a model to fill this gap by quantifying the phylogenetic structure of the insular faunas and disentangling the relative contribution of potential drivers. Firstly, we used clustering and ordination analyses to identify regional species pools. To test whether the insular freshwater fish faunas were phylogenetically clustered or overdispersed, we calculated the net relatedness index (NRI) and the nearest taxon index (NTI). Finally, we implemented logistic regressions to disentangle the relative importance of species attributes (i.e. maximum body length, climatic niche dissimilarity, and diversification) and historical connectivity in explaining the insular faunas. Our results showed that the most possible species pools of Taiwan are Zhejiang and Fujian provinces, and those of Hainan are Guangdong and Guangxi provinces. These insular faunas showed random phylogenetic structures in terms of NRI values. According to the NTI values, however, the Taiwanese fauna displayed more phylogenetic clustering, while the Hainanese one was more overdispersed. Both the standard and phylogenetic logistic regressions identified historical connectivity and climatic niche dissimilarity as the two top explanatory variables for species assemblages on these islands. Our reconstruction of the paleo-connected drainage basins provides insight into how historical processes and ecological factors interact to shape the freshwater fish fauna of the East Asian islands.

Genetic diversity of Rhinogobius delicatus (Perciformes: Gobiidae): origins of the freshwater fish in East Taiwan

Mitochondrial DNA cytochrome b and d-loop sequences (1,984 bp) from 92 specimens of the freshwater goby Rhinogobius delicatus from seven drainages in East Taiwan were identified as two major lineages exhibiting a southern or northern distribution. The existence of low genetic diversity, a pattern of population decline and high population differentiation (F _ST=0.711) support the need for the development of management strategies for the conservation of localized populations. The results of a statistical dispersal-vicariance analysis suggested that the ancestral populations of R. delicatus were widely distributed in East Taiwan. Compared with the phylogeographic patterns of the other endemic eastern Taiwan freshwater fishes, Onychostoma alticorpus, Aphyocypris kikuckii and Hemimyzon taitungensis, our study suggests that the freshwater fishes colonized East Taiwan through northeastern and southwestern Taiwan, although the ancestral populations colonized the island before it reached its present shape.

Are population isolations and declines a threat to island endemic water striders? A lesson from demographic and niche modelling of Metrocoris esakii (Hemiptera: Gerridae)

Genetic stochasticity and bottlenecking in the course of Pleistocene glaciations have been identified as threatening the survival of local endemics. However, the mechanisms by which local endemic species balance the influences of these two events remain poorly understood. Here, we generated a double-digest restriction site-associated DNA sequencing (ddRAD-seq) data set, mined mitochondrial sequences and constructed ecological niche models for the island endemic water strider Metrocoris esakii (Hemiptera: Gerridae). We found that M. esakii comprised three divergent lineages (i.e., north, central and south) isolated by geographical barriers and generally experienced population declines with the constriction of suitable areas during the Last Glacial Maximum (LGM). Further demographic model testing and stairway plots revealed a history of recent gene flow among the neighbouring lineages and rapid recovery at the end of the LGM, indicating that M. esakii at least had the potential for an adaptive response to population fragmentation and bottlenecking. The northern lineage did not show genetic bottlenecking during the LGM, which was probably due to its large effective population size (N_e ) from migration, which improved its adaptive potential. Relative to the ddRAD-seq data set, the demographic results based on mitochondrial sequences were less conclusive, showing weak differentiation and oversimplified demographic trajectories for the three genetic lineages. Overall, this study provides some degree of optimism for the survival of island endemic water striders from a demographic perspective, but further evaluation of their extinction risk under the impacts of human activities is required.

From mainland to islands: colonization history in the tree frog Kurixalus (Anura: Rhacophoridae)

The origin and colonization history of Kurixalus, a genus of small arboreal tree frogs breeding exclusively in shallow swamps, is under disputed. On the basis of comprehensive sampling program, the evolutionary history of Kurixalus is investigated based on 3 mitochondrial genes. Our results indicate that the genus Kurixalus originated in the Asian mainland and subsequently arrived at its current distribution in Borneo, Taiwan, Ryukyu, and Hainan islands by a series of dispersal events. Moreover, the colonization of Taiwan from mainland Asia has occurred 2 times. The initial colonization of Taiwan occurred at 3.46–8.68 Mya (95% highest posterior density), which rejects the hypothesis that Kurixalus probably originated from Taiwan during the early Oligocene and favors the model of Neogene-origin rather than the model of Quaternary-origin for Taiwanese Kurixalus. Kurixalus eiffingeri has dispersed from Taiwan to the Ryukyus once or 2 times pending more data. Both transoceanic dispersal and landbridge dispersal have played a role in the colonization process; the former resulted in the colonization of Taiwan and the Ryukyus and the latter led to the colonization of Borneo and Hainan.

Geographical and temporal origins of Neocaridina species (Decapoda: Caridea: Atyidae) in Taiwan

BACKGROUND

The freshwater species on Taiwan Island have been documented to have originated from mainland China and the Japanese islands from multiple events and by multiple colonization routes. Moreover, the sequences from the mitochondrial DNA cytochrome c oxidase subunit I (COI) have been used for DNA barcoding to identify the species. This study used the COI sequences to identify Neocaridina species in Taiwan and to examine their geographical and temporal origins.

RESULTS

In total, 479 specimens were collected from 35 localities, which covered almost all rivers in Taiwan. In addition, some sequences were downloaded from GenBank. The maximum likelihood (ML) tree displayed that all sequences were sorted into 13 taxa (clades), and all sequences in Taiwan were sorted into four clades. The Bayesian skyline plots revealed that these four Neocaridina species have declined recently in Taiwan.

CONCLUSIONS

All results support that (1) there are four Neocaridina species in Taiwan, which are N. davidi, N. saccam, N. ketagalan and an undescribed Neocaridina species (N. sp.); (2) these four species colonized Taiwan Island in four colonization events; (3) N. sp. colonized Taiwan first; (4) after the island reached its shape, N. ketagalan and N. saccam colonized Taiwan from the Japanese islands and mainland China, respectively; (5) N. davidi colonized northern Taiwan last; and (6) the cyclic glacial and landform changes in East Asia shaped the colonization events and population structures of the Neocaridina species.

Population structure of Aphyocypris normalis: phylogeography and systematics

Aphyocyprisnormalis (Teleostei: Cyprinidae) is an endemic species in South China, but little is known about its genetic structure. This study examined the population structure of A.normalis using sequences of the mitochondrial DNA control region and cytochrome b gene (2,086 bp). In total, 107 specimens were collected from nine populations. All 105 mtDNA haplotypes were identified as belonging to two allopatric phylogroups. The results of a statistical dispersal-vicariance analysis (S-DIVA) suggested that the ancestral populations of A.normalis were distributed widely on Hainan Island and east of the Leizhou Peninsula. A comparison of the fixation indices N_ST (0.532) and G_ST (0.004) revealed that the phylogeny and geography had a significant relationship. Our study found that (1) the Wuzhishan and Yinggeling Mountain Range was an important barrier limiting gene exchange between populations on both sides; (2) cyclic climate changes may have shaped migrations and population differentiations; and (3) different colonization times caused different population diversities between codistributed species. In addition, the inter- and intraspecific diversities of the genus Aphyocypris were estimated.

Genetic variation of the land-locked freshwater shrimp Caridina pseudodenticulata (Decapoda: Atyidae: Caridina) in Taiwan

Caridina pseudodenticulata is a land-locked freshwater shrimp, which is endemic to North and West Taiwan. It is listed as Vulnerable under criteria A1ce in the IUCN Red List. Our study used mitochondrial DNA fragment sequences of the cytochrome oxidase subunit I (COI) gene to examine its genetic structure. In total, 204 sequences were analyzed from 22 populations, and all 39 haplotypes were identified. Our study (1) found that the genetic diversity within population of this species was similar to that of other land-locked freshwater prawn in Taiwan; (2) displayed that during glaciation, the Taiwan Strait was largely above water, an event which created opportunities for colonization Taiwan and dispersal widely; (3) showed that although the phylogenetic analysis showed the lack of a population genetic structure, the AMOVA and haplotypes distribution pattern revealed geographically divided; (4) indicated that C. pseudodenticulata could across to neighbor rivers through the mouths of rivers, which confluent with each other during floods.

Phylogeography and genetic structure of the oriental river prawn Macrobrachium nipponense (Crustacea: Decapoda: Palaemonidae) in East Asia

The oriental river prawn (Macrobrachium nipponense) is mainly distributed in East Asia. The phylogeography, population genetic structure and historical demography of this species in the East Asia were examined by using partial sequences of the cytochrome oxidase subunit I (COI) and 16S rRNA in mitochondrial DNA. Ten populations that included 239 individuals were collected from Taiwan (Shihmen Reservoir, SMR, Mingte Reservoir, MTR and Chengching Lake Reservoir, CLR), mainland China (Taihu Lake, TLC, Min River, MRC, Jiulong River, JRC and Shenzhen Reservoir, SRC), Japan (Biwa Lake, BLJ and Kasumigaura Lake, KLJ) and Korea (Han River, HRK). The nucleotide diversity (π) of all individuals was 0.01134, with values ranging from 0.0089 (BLJ, Japan) to 0.01425 (MTR, Taiwan). A total of 83 haplotypes were obtained, and the haplotypes were divided into 2 main lineages: lineage A included the specimens from BLJ, KLJ, CLR, MTR, TLC, MRC and JRC, and lineage B comprised the ones from HRK, SRC, SMR, MTR, TLC, MRC and JRC. Lineage A could be further divided two sub-lineages (A₁ and A₂). Individuals of lineage A₂ were only from TLC. Demographic expansion was observed in each lineage, starting within the second-to-latest interglacial period for lineage A and within the last glacial period for lineage B. All F_ST values among the ten populations were significantly different, except for the values between MRC and JRC, and SMR and SRC. The phylogeography and genetic structure of M. nipponense in East Asia might be influenced by Pleistocene glacial cycles, lake isolation and human introduction. The possible dispersal routes of M. nipponense in the East Asia were also discussed.

The population history of Garra orientalis (Teleostei: Cyprinidae) using mitochondrial DNA and microsatellite data with approximate Bayesian computation

Background

The South China landmass has been characterized by a complex geological history, including mountain lifting, climate changes, and river capture/reversal events. To determine how this complexity has influenced the landmass’s phylogeography, our study examined the phylogeography of Garra orientalis, a cyprinid widely distributed in South China, using sequences from the mitochondrial DNA control region and cytochrome b gene (1887 bp) and polymorphisms of thirteen microsatellite loci.

Results

In total, 157 specimens were collected from eight populations. All 88 mtDNA haplotypes were identified as belonging to three major lineages, and these lineages were almost allopatric in their distributions. The results of a statistical dispersal-vicariance analysis suggested that the ancestral populations of G. orientalis were distributed south of the Yunkai Mountains, including on Hainan Island. The mtDNA data revealed a strong relationship between phylogeny and geography. In the microsatellite analysis, a total of 339 alleles with an average of 26 alleles per locus were observed across thirteen microsatellite loci. A clustering algorithm for microsatellite data revealed an admixture-like genetic structure. Although the mtDNA and microsatellite data sets displayed a discordant population structure, the results of an approximate Bayesian computation approach showed that these two markers revealed congruent historical signals. The population history of G. orientalis reflects vicariance events and dispersal related to the complex geological history of South China.

Conclusion

Our results (i) found that the discordances between mtDNA and microsatellite markers were accounted for by admixtures; (ii) showed that the Wuzhishan and Yinggeling mountain ranges and Qiongzhou Strait were important barriers limiting gene exchange between populations on both sides; (iii) indicated that during glaciation and inter-glacial periods, the strait and continental shelves were exposed and sank, which contributed with the dispersion and differentiation of populations; and (iv) displayed that the admixtures between lineages took place in coastal populations and then colonized the tributaries of the Pearl River.

Origins of Semisulcospira libertina (gastropoda: semisulcospiridae) in Taiwan

The most accepted hypothesis has suggested that the fauna in Taiwan Island originated from South China, but some studies supported the Japan, Ryukyu Archipelago, and Taiwan Islands as a unique biogeographical district. This study examines whether the populations of freshwater snail Semisulcospira libertina in Taiwan are closer to those in Japan based on the mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. Our study shows the populations in North Taiwan originated from Japan and the cyclic glacial caused the migrations among islands and continent repeatedly; the populations in South Taiwan might originate from South China or South Asia. Our results will not only affect the conclusions in phylogeography of freshwater species in Taiwan but also change the sampling plans in the future studies about evolutionary of freshwater species in East Asia.

A comparative phylogeographic study reveals discordant evolutionary histories of alpine ground beetles (Coleoptera, Carabidae)

Taiwan, an island with three major mountain ranges, provides an ideal topography to study mountain-island effect on organisms that would be diversified in the isolation areas. Glaciations, however, might drive these organisms to lower elevations, causing gene flow among previously isolated populations. Two hypotheses have been proposed to depict the possible refugia for alpine organisms during glaciations. Nunatak hypothesis suggests that alpine species might have stayed in situ in high mountain areas during glaciations. Massif de refuge, on the other hand, proposes that alpine species might have migrated to lower ice-free areas. By sampling five sympatric carabid species of Nebria and Leistus, and using two mitochondrial genes and two nuclear genes, we evaluated the mountain-island effect on alpine carabids and tested the two proposed hypotheses with comparative phylogeographic method. Results from the phylogenetic relationships, network analysis, lineage calibration, and genetic structure indicate that the deep divergence among populations in all L. smetanai, N. formosana, and N. niitakana was subjected to long-term isolation, a phenomenon in agreement with the nunatak hypothesis. However, genetic admixture among populations of N. uenoiana and some populations of L. nokoensis complex suggests that gene flow occurred during glaciations, as a massif de refuge depicts. The speciation event in N. niitakana is estimated to have occurred before 1.89 million years ago (Mya), while differentiation among isolated populations in N. niitakana, N. formosana, L. smetanai, and L. nokoensis complex might have taken place during 0.65-1.65 Mya. While each of the alpine carabids arriving in Taiwan during different glaciation events acquired its evolutionary history, all of them had confronted the existing mountain ranges.

Molecular Evidence for Multiple Origins of the European Spined Loaches (Teleostei, Cobitidae)

We present a phylogenetic investigation of the Northern Clade, the major monophyletic clade within the freshwater fish family Cobitidae, one of the most prominent families of freshwater fishes found in Asian and European waters. Phylogenetic reconstructions based on the cytochrome b and RAG-1 genes show the genera Microcobitis, Sabanejewia, Koreocobitis and Kichulchoia as monophyletic groups. These reconstructions also show a Cobitis sensu lato and a Misgurnus sensu lato group. The Cobitis sensu lato group includes all species of Cobitis, Iksookimia, Niwaella and Kichulchoia, while the Misgurnus sensu lato group includes Misgurnus, Paramisgurnus and Koreocobitis. Although the monophyly of both the Cobitis sensu lato and Misgurnus sensu lato groups is supported, relationships within the groups are incongruent with current generic definitions. The absence of monophyly of most genera included in the Cobitis sensu lato group (Cobitis, Iksookimia and Niwaella) or their low genetic differentiation (Kichuchoia) supports their consideration as synonyms of Cobitis. Molecular phylogenies indicate that the Asian species of Misgurnus experienced a mitochondrial introgression from a lineage of Cobitis. We also find two nuclear haplotypes in the same Cobitis species from the Adriatic area that, in the absence of morphological differentiation, may indicate molecular introgression. Most lineages within the Northern Clade consist of species found in East Asia. However, some lineages also contain species from Europe and Asia Minor. The phylogenetic relationships presented here are consistent with previous studies suggesting an East Asian origin of the Northern Clade. According to the current distributions and phylogenetic relationships of the Misgurnus sensu lato and Cobitis clade lineages, particularly of M. fossilis and C. melanoleuca, the range expansion of East Asian species into Europe was most likely via Siberia into Northern and Central Europe. Phylogenetic analyses also show that the Cobitis sensu lato group consists of two clear subgroups (I and II), each presenting geographical differences. Subgroup I is distributed exclusively in East Asian drainages with an Eastern European offshoot (C. melanoleuca), whereas Subgroup II includes species widespread throughout Europe (including the Mediterranean), Asia Minor, the Black Sea and the Caucasus, with some lineages related to species restricted to East Asia.

Population Structure and Historical Demography of the Oriental River Prawn (Macrobrachium nipponense) in Taiwan

The oriental river prawn (Macrobrachium nipponense) is a non-obligatory amphidromous prawn, and it has a wide distribution covering almost the entire Taiwan. Mitochondrial DNA fragment sequences of the cytochrome oxidase subunit I (COI) and 16S rRNA were combined and used to elucidate the population structure and historical demography of oriental river prawn in Taiwan. A total of 202 individuals from six reservoirs and three estuaries were separately collected. Nucleotide diversity (π) of all populations was 0.01217, with values ranging from 0.00188 (Shihmen Reservoir, SMR, northern Taiwan) to 0.01425 (Mingte Reservoir, MTR, west-central Taiwan). All 76 haplotypes were divided into 2 lineages: lineage A included individuals from all sampling areas except SMR, and lineage B included specimens from all sampling locations except Chengching Lake Reservoir (CLR) and Liyu Lake Reservoir (LLR). All F_ST values among nine populations were significantly different except the one between Jhonggang River Estuary (JGE, west-central Taiwan) and Kaoping River Estuary (KPE, southern Taiwan). UPGMA tree of nine populations showed two main groups: the first group included the SMR and Tamsui River Estuary (TSE) (both located northern Taiwan), and the second one included the other seven populations (west-central, southern and eastern Taiwan). Demographic analyses implied a population expansion occurred during the recent history of the species. The dispersal route of this species might be from China to west-central and west-southern Taiwan, and then the part individuals belonging to lineage A and B dispersed southerly and northerly, respectively. And then part individuals in west-central Taiwan fell back to and stay at estuaries as the sea level rose about 18,000 years ago.

Differentiation in stag beetles, Neolucanus swinhoei complex (Coleoptera: Lucanidae): four major lineages caused by periodical Pleistocene glaciations and separation by a mountain range

Taxonomic debates on Neolucanus swinhoei complex consisting of N. swinhoei, N. doro doro, N. doro horaguchii, and N. euganiae, distributed exclusively in Taiwan, have been ongoing for several decades because of their overlapping morphological characters. To clarify their taxonomic status and phylogeographical history, we analyzed nine morphological characteristics and four molecular amplicons. Phylogenetic inferences based on COI+16S rDNA+wingless showed one eastern and three western lineages, with the latter consisting of one low-hill and two montane lineages. Intermingled DNA sequences from different populations within each lineage, many low FST values, and a high variance component between lineages indicate the possibility of gene flow among populations. However, positive relationships were observed between the genetic divergences of 16S rDNA and its FST values with geographic distance. A divergence estimation based on COI+16S revealed that these beetles might have originated from Asian mainland and differentiated into western and eastern lineages ca. 1Mya, with the differentiation of the western lineages occurring approximately 0.50-0.75Mya. Isolation by mountain ranges and limited flying capability of these beetles as well as populations retreat to and expansion from refugia in response to glaciation cycles have resulted in the current distribution of N. swinhoei complex. Although most morphological characters are variable and undistinguishable, multi-dimensional scaling analysis based on measurable characteristics could recognize hill N. swinhoei as a cluster distinct from the others. However, based on the realities of genetic admixture, shared phylogeographical history and overlapping characteristics, all of these stag beetles should be regarded as Neolucanus swinhoei Bates, 1866.

Phylogeographical analysis on Squalidus argentatus recapitulates historical landscapes and drainage evolution on the island of Taiwan and mainland China

Phylogeographical analyses on Squalidus argentatus samples from thirteen localities within mainland China and Taiwan were conducted for biogeographic studies, as their dispersal strictly depends on geological evolution of the landmasses. A total of 95 haplotypes were genotyped for mtDNA cyt b gene in 160 specimens from nine river systems. Relatively high levels of haplotype diversity (h = 0.984) and low levels of nucleotide diversity (π = 0.020) were detected in S. argentatus. Two major phylogenetic haplotype groups, A and B, were revealed via phylogenetic analysis. The degree of intergroup divergence (3.96%) indicates that these groups diverged about 4.55 myr (million years) ago. Haplotype network and population analyses indicated significant genetic structure (F_ST = 0.775), largely concordant with the geographical location of the populations. According to SAMOVA analysis, we divided these populations into four units: Yangtze-Pearl, Qiantang-Minjiang, Jiulong-Beijiang and Taiwan groups. Mismatch distribution analysis, neutrality tests and Bayesian skyline plots indicated a significant population expansion for lineage A and B, approximately dated 0.35 and 0.04 myr ago, respectively. We found strong geographical organization of the haplotype clades across different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow.

Lineage-specific late pleistocene expansion of an endemic subtropical gossamer-wing damselfly, Euphaea formosa, in Taiwan

BACKGROUND

Pleistocene glacial oscillations have significantly affected the historical population dynamics of temperate taxa. However, the general effects of recent climatic changes on the evolutionary history and genetic structure of extant subtropical species remain poorly understood. In the present study, phylogeographic and historical demographic analyses based on mitochondrial and nuclear DNA sequences were used. The aim was to investigate whether Pleistocene climatic cycles, paleo-drainages or mountain vicariance of Taiwan shaped the evolutionary diversification of a subtropical gossamer-wing damselfly, Euphaea formosa.

RESULTS

E. formosa populations originated in the middle Pleistocene period (0.3 Mya) and consisted of two evolutionarily independent lineages. It is likely that they derived from the Pleistocene paleo-drainages of northern and southern Minjiang, or alternatively by divergence within Taiwan. The ancestral North-central lineage colonized northwestern Taiwan first and maintained a slowly growing population throughout much of the early to middle Pleistocene period. The ancestral widespread lineage reached central-southern Taiwan and experienced a spatial and demographic expansion into eastern Taiwan. This expansion began approximately 30,000 years ago in the Holocene interglacial period. The ancestral southern expansion into eastern Taiwan indicates that the central mountain range (CMR) formed a barrier to east-west expansion. However, E. formosa populations in the three major biogeographic regions (East, South, and North-Central) exhibit no significant genetic partitions, suggesting that river drainages and mountains did not form strong geographical barriers against gene flow among extant populations.

CONCLUSIONS

The present study implies that the antiquity of E. formosa's colonization is associated with its high dispersal ability and larval tolerance to the late Pleistocene dry grasslands. The effect of late Pleistocene climatic changes on the subtropical damselfly's historical demography is lineage-specific, depending predominantly on its colonization history and geography. It is proposed that the Riss and Würm glaciations in the late Pleistocene period had a greater impact on the evolutionary diversification of subtropical insular species than the last glacial maximum (LGM).