Molecular phylogeography of endangered sharp-snouted pitviper (Deinagkistrodon acutus; Reptilia, Viperidae) in Mainland China

Genomic analysis reveals deep population divergence in the water snake Trimerodytes percarinatus (Serpentes, Natricidae)

Although several phylogeographic studies of Asian snakes have been conducted, most have focused on pitvipers, with non-venomous snakes, such as colubrids or natricids, remaining poorly studied. The Chinese keelback water snake (Trimerodytes percarinatus Boulenger) is a widespread, semiaquatic, non-venomous species occurring in China and southeastern Asia. Based on mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, we explored the population genetic structure, genetic diversity, and evolutionary history of this species. MtDNA-based phylogenetic analysis showed that T. percarinatus was composed of five highly supported and geographically structured lineages. SNP-based phylogenetic analysis, principal component analysis, and population structure analysis consistently revealed four distinct, geographically non-overlapping lineages, which was different from the mtDNA-based analysis in topology. Estimation of divergence dates and ancestral area of origin suggest that T. percarinatus originated ~12.68 million years ago (95% highest posterior density: 10.36-15.96 Mya) in a region covering southwestern China and Vietnam. Intraspecific divergence may have been triggered by the Qinghai-Xizang Plateau uplift. Population demographics and ecological niche modeling indicated that the effective population size fluctuated during 0.5 Mya and 0.002 Mya. Based on the data collected here, we also comment on the intraspecific taxonomy of T. percarinatus and question the validity of the subspecies T. p. suriki.

Comparative analysis of Deinagkistrodon acutus venom from Taiwan and China utilizing chromatographic, electrophoretic, and bioinformatic approaches, along with ELISA employing a monospecific antivenom

Deinagkistrodon acutus is a medically important pitviper inhabiting mainly South China and Taiwan. The hemorrhagic effects of its envenoming are compatible to its venom, which is abundant in metalloproteases (svMPs) and C-type lectin-like proteins. In this study, we investigated geographic variations in the venom of D. acutus collected from Taiwan and four Mainland Chinese provinces: Fujian, Jiangxi, Anhui, and Hunan. The variations were assessed through high-performance liquid chromatography, non-metric multidimensional scaling analysis, gel electrophoresis, and enzyme-linked immunosorbent assay (ELISA) with a monospecific antivenom (DaMAV) generated against the Taiwanese D. acutus venom, and discussed based on venom-protein sequences in databases and literature related to D. acutus venom. Additionally, the cross-reactivity of DaMAV against Crotalus horridus and Calloselasma rhodostoma venoms was investigated. We noted differential abundances of D. acutus venom metalloproteases, C-type lectin-like proteins, and phospholipase A2, along with point mutations and selective expression of serine protease isoforms. The ELISA results revealed that the venom from Taiwan was more reactive toward Taiwanese DaMAV than the four Mainland Chinese venoms, consistent with chromatographic profile differences, whereas C. horridus venom presented moderate cross-reactivity with DaMAV. The observed immunoreactivities of these venom with DaMAV can be attributed to the high prevalence of their PIII-svMPs, which are the dominant antigens, and the conservation of PIII-svMP epitopes.

Sharp-nosed Pit Viper (Deinagkistrodon acutus) from Taiwan and China: A comparative study on venom toxicity and neutralization by two specific antivenoms across the Strait

In East Asia, the Sharp-nosed Pit Viper (Deinagkistrodon acutus) is a medically important venomous snake in Taiwan and China, two geographical areas long separated by the Taiwan Strait. Yet, snake venom variation is little known between specimens found across the Strait. This study thus investigated the intra-species variation of D. acutus venoms from Taiwan (Da-Taiwan) and China (Da-China) in their profiles of gel electrophoresis, toxicity, immunoreactivity and neutralization effect by antivenom. Da-China venom exhibited higher procoagulant, hemorrhagic and lethal activities than Da-Taiwan venom, presumably attributed to the higher abundance of moderate-to-high molecular weight toxins (procoagulants and hemorrhagins) in the venom. The mono-specific antivenoms produced in Taiwan (DaMAV-Taiwan) and China (DaMAV-China) were immunoreactive toward both venoms, and were able to neutralize the venom toxicity to different extents. DaMAV-Taiwan was more efficacious in neutralizing the venom procoagulant and lethal effects, while DaMAV-China was more potent against hemorrhagic effect. The discrepancy in efficacy between the two antivenoms could be due to varying proportions of neutralizing antibodies in the respective products, influenced by techniques of antibody raising and purification. Further study is warranted to elucidate variation in the proteome and antigenicity of D. acutus venom between snakes from Taiwan and China.

Population Genetic Structure of Chlorops oryzae (Diptera, Chloropidae) in China

Frequent outbreaks have made Chlorops oryzae one of the major pests of rice in some regions. In order to understand the ecological adaptation of C. oryzae at the molecular level, and provide a scientific basis for formulating management strategies, we used two molecular markers, COI and ITS1 sequences, to systematically analyze the genetic structure of 31 populations. The higher haplotype diversity and lower nucleotide diversity indicated that the C. oryzae populations experienced rapid expansion after a “Bottleneck effect”. The results of the mismatch distribution, neutrality test (Fu’s Fs < 0, p < 0.001), and haplotype network analysis suggested that the population has recently undergone an expansion. Although genetic differentiation among C. oryzae populations was found to have existed at low/medium levels (Fst: 0.183 for COI, 0.065 for ITS1), the frequent gene flow presented as well (Nm: 2.23 for COI, 3.60 for ITS1) was supposed to be responsible for frequent local outbreaks.

Recombinase polymerase amplification combined with fast DNA extraction for on–spot identification of Deinagkistrodon acutus, a threatened species

This study addresses the use of recombinase polymerase amplification combined with fast DNA extraction for on–spot identification of Deinagkistrodon acutus, a snake species threatened due to over–exploitation and habitat destruction. For its conservation, an efficient species identification method is urgently neededto fight against illegal capture and trade. Fourteen individuals representing 12 snake species (including D. acutus and other snake species) were collected from mountainous regions in Southern China. Genomic DNA was extracted within five minutes by a modified alkaline lysis method. Species–specific primers for recombinase polymerase amplification (RPA) were designed based on the sequences of cytochrome C oxidase subunit I (COI) barcode region, and an optimized RPA assay system was set up. Specificity and sensitivity of the assay were checked, and the assay was validated by identifying 10 commercial Qi She crude drug samples derived from D. acutus. Under optimized RPA conditions, a distinct single band of 354 bp was amplified only for D. acutus but not for the related snake species. The entire procedure can be completed in 30 min at room temperature. Commercial Qi She crude drug identification validated effectiveness of the established assay system. Using a recombinase polymerase amplification (RPA) assay with rapid DNA extraction, we established an on–spot D. acutus identification method with good specificity and sensitivity. This method could become an efficient tool for rigorous supervision of illegal D. acutus capture and trade.

Cryptic diversity and phylogeography of the Rhabdophis nuchalis group (Squamata: Colubridae)

Previous studies, have found that the rapid uplift of the Tibetan plateau accelerated the diversification of species. However, there are few relevant biogeographic data for the Colubridae in this region. We conducted a comprehensive study of the Rhabdophis nuchalis Group, which presently contains four nominal species, R. nuchalis, R. pentasupralabialis, R. leonardi, and R. chiwen. Building upon previous studies with specimens we have recently examined, greater interspecific and intraspecific diversity has been revealed. Here we address three questions: (1) Do the intraspecific differences represent only geographic variation within lineages, or are there cryptic species? (2) What are the interspecific relationships among members of the R. nuchalis Group? (3) What has been the biogeographic history of this species group? To resolve these questions we used four mitochondrial gene sequences and one nuclear sequence to investigate the molecular phylogenetic and geographic relationships among populations. Our molecular analysis reveals cryptic species diversity within the R. nuchalis Group, and seven clades were identified in the analysis. Ancestral area estimation suggests that the R. nuchalis Group originated in the Hengduan Mountains approximately 6.24 Mya and expanded its range northward to the Qinling-Daba Mountains. The Sichuan Basin appears to have been a barrier to migration. Species divergence seems to have been related to the rapid uplift of the Qinghai-Tibet Plateau.

Climate warming will increase chances of hybridization and introgression between two Takydromus lizards (Lacertidae)

Coexisting species may experience population and range changes alone or jointly in response to environmental change. Here, we used six climate variables and ten modeling algorithms to predict the distribution of two Takydromus species (T. septentrionalis and T. sexlineatus) in China. We identified the sympatric and allopatric areas by comparing projections between the two species based on habitat suitability under present and future climate scenarios. We constructed the hypervolumes of six climate variables for the two species and then evaluated overlaps between hypervolumes. From this study, we know the following. First, minimum temperature of coldest month contributes the most to the prediction of habitat suitability. Second, habitats suitable for the two species will shift northward in response to climate warming. Third, the range of T. sexlineatus will expand across the four future time intervals before 2,100, namely the 2021-2040, 2041-2060, 2061-2080, and 2081-2100 intervals, under both Shared socioeconomic pathway (SSP) 245 and SSP585 scenarios, and the range of T. septentrionalis will also expand in the future except at the 2081-2100 interval under the SSP585 scenario. Fourth, the sympatric areas will contract or expand under the SSP245 scenario and expand across the four future time intervals before 2,100 under the SSP585 scenario. Fifth, the niche hypervolumes of the two species partially overlapped, and the differences in niche centroid show some degree of niche differentiation between the two species. These results allow to conclude that climate warming will not only drive the northward drift of sympatric areas but also increase the size of these areas if nothing is done to limit the emission of greenhouse gases. Given the existence of hybridization and introgression between T. septentrionalis and T. sexlineatus in the field where they coexist, we also conclude that climate warming will increase chances of hybridization and introgression between the two species.

Protective Effects of Carbon Dots Derived from Phellodendri Chinensis Cortex Carbonisata against Deinagkistrodon acutus Venom-Induced Acute Kidney Injury

BACKGROUND

As an emerging nanomaterial, carbon dots (CDs) have been the focus of tremendous attention for biomedical applications. However, little information is available on their bioactivity of inhibiting acute kidney injury (AKI) induced by snake venom.

METHODS

This study reports the development of a green, one-step pyrolysis process to synthesize CDs using Phellodendri Chinensis Cortex (PCC) as the sole precursor, and their potential application as a protectant against Deinagkistrodon acutus (D. acutus) venom-induced AKI was investigated for the first time. The AKI model was established by injecting D. acutus venom into the abdominal cavity of mice and the potential protective effects of PCC Carbonisata-CDs (PCCC-CDs) on renal abnormalities including dysfunction, inflammatory reactions, tissue damage, and thrombocytopenia at six time points (1, 3, and 12 h, and 1, 2, and 5 days) were investigated.

RESULTS

These results demonstrated that PCCC-CDs significantly inhibited the kidney dysfunction (reduced serum creatinine (SCR), blood urea nitrogen (BUN), urinary total protein (UTP), and microalbuminuria (MALB) concentrations) and the production of chemoattractant (monocyte chemotactic protein 1 (MCP-1)), proinflammatory cytokines (interleukin (IL)-1β), and anti-inflammatory cytokine (IL-10) in response to intraperitoneal injection of D. acutus venom. The beneficial effect of PCCC-CDs on the envenomed mice was similar to that on the change in renal histology and thrombocytopenia.

CONCLUSIONS

These results demonstrated the remarkable protective effects of PCCC-CDs against AKI induced by D. acutus venom, which would not only broaden the biomedical applications of CDs but also provide a potential target for the development of new therapeutic drugs for AKI induced by D. acutus snakebite envenomation.

Complete mitochondrial genome of the sharp-snouted pitviper Deinagkistrodon acutus (Reptilia, Viperidae)

The sharp-snouted pitviper, Deinagkistrodon acutus, belongs to the monotypic genus Deinagkistrodon of the family Viperidae. It is one of the important species of snakes for its commercial and medicinal value. Complete mitochondrial genome sequence of the D. acutus was sequenced by next-generation sequencing technology. The total mitochondrial genome is 17,538 bp in length, containing 13 protein-coding genes, 22 tRNA genes, 2 ribosome RNA genes, and 2 control regions. Most of the genes of D. acutus are encoded on the H-strand, except for the ND6 subunit gene and eight tRNA genes which distribute on the L-strand. Phylogenetic reconstruction result indicated that our newly determined mitochondrial genome sequence could meet the demands. The complete mitochondrial genome sequence presented here will be useful to study the evolutionary relationships and genetic diversity of D. acutus.

Genetic structure of two sympatric gudgeon fishes (Xenophysogobio boulengeri and X. nudicorpa) in the upper reaches of Yangtze River Basin

Background

Xenophysogobio boulengeri and X. nudicorpa are the only two species within the genus Xenophysogobio (Cyprinidae, Cypriniformes), and both are endemic to the upper reaches of the Yangtze River. In recent years, due to human activities, the natural resources available to both species have declined sharply. Sympatric species with overlapping niches inevitably compete for their habitats, and genetic structure and diversity can reflect population history and their potential for adaptation to changing environments, which is useful for management decisions.

Methods

In the present study, microsatellite DNA and mitochondrial DNA (mtDNA) markers were used to investigate the patterns of population genetic structure for X. boulengeri and X. nudicorpa. Microsatellite DNA data, jointly with traditional summary statistics including F _ST and F _is, were used to assess the population genetic structure by structure analysis. The mtDNA sequences were then used to examine these patterns through time to detect demographic history.

Results

Xenophysogobio boulengeri and X. nudicorpa exhibited high levels of genetic diversity in Yangtze River populations, except for two populations of X. nudicorpa in the Jinsha River, which were low in mtDNA diversity. X. boulengeri showed genetic homogeneity among populations, whereas X. nudicorpa appeared to have significant geographic genetic divergence. Both species experienced a late-Pleistocene sudden population expansion in Yangtze River populations, but not in the Jinsha River populations of X. nudicorpa.

Discussion

The genetic homogeneity of X. boulengeri populations might result from similar population expansion events and environment features. The geographic genetic subdivision for X. nudicorpa between the Jinsha and Yangtze Rivers might be caused by the geographic isolation in the middle Pliocene, as well as climate and environmental heterogeneity.

Multilocus phylogeography of the brown-spotted pitviper Protobothrops mucrosquamatus (Reptilia: Serpentes: Viperidae) sheds a new light on the diversification pattern in Asia

Understanding the influence of geographical events and climate changes on genetic diversity is essential in explaining current patterns of genetic structure and geographic distribution of organisms. We inferred phylogenetic relationships, investigated historical demography, explored the evolutionary history, and clarified intraspecific taxonomy of Protobothrops mucrosquamatus, which is one of the commonest and most wide-ranging Asian pitvipers. A total of 184 samples from 54 localities were sequenced and analyzed for two mitochondrial gene fragments and two nuclear genes. Phylogenetic reconstruction based on mtDNA sequences revealed the existence of a minimum of five geographically structured and well-supported lineages within P. mucrosquamatus. Based on the mtDNA gene tree, and the geographic relationship between populations allied by matrilineal lineages, a complex longitudinal and latitudinal diversification pattern was uncovered in P. mucrosquamatus. The estimated date of the origin of the species (about 5.3 Ma) and divergence of the intraspecific lineages match the rapid uplifting of Qinghai-Xizang Plateau, and is also consistent with those of some other co-distributed Asian pitvipers. Formation of the two island lineages (Taiwan and Hainan) was generally congruent with the first isolation of the islands, but the two lineages showed different relationships with the continental Asian populations in comparison with some other pitvipers. Population historical demographic analyses, based on three methods, showed that all lineages have experienced slight population expansion in and around the Dali Glacial. Tests of intraspecific taxonomy indicated that no cryptic taxon is present within this widely distributed snake.

Three-dimensional computed tomography reconstructive diagnosis of snakebite-induced cerebral infarction

Despite common injury caused by snakebite, snakebite-induced ischemic stroke is rare. We reported on a patient who incurred a large cerebral infarction after being bitten by a Deinagkistrodon acutus, one of the most poisonous snakes in the southwestern of China. Applying 3D computed tomography (CT) of head combined with cerebral angiography examinations showed a large cerebral infarction, hernia in the right brain, developmental abnormalities of the right middle cerebral artery and cerebral artery of right brain. In conclusion, head CT imaging combined with cerebral angiography provides an efficient approach in diagnosis of stroke caused by snakebites.

Phylogeographic structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China

BACKGROUND

Species with broad distributions frequently divide into multiple genetic forms and may therefore be viewed as "cryptic species". Here, we used the mitochondrial cytochrome b (Cytb) and 12 nuclear DNA loci to investigate phylogeographic structures of the sharpbelly (Hemiculter leucisculus) in rivers in southern China and explored how the geological and climatic factors have shaped the genetic diversity and evolutionary history of this species.

RESULTS

Our mitochondrial phylogenetic analysis identified three major lineages (lineages A, B, and C). Lineages B and C showed a relatively narrower geographic distribution, whereas lineage A was widely distributed in numerous drainages. Divergence dates suggested that H. leucisculus populations diverged between 1.61-2.38 Ma. Bayesian species delimitation analysis using 12 nuclear DNA loci indicated the three lineages probably represented three valid taxa. Isolation-with-migration (IM) analysis found substantial gene flow has occurred among the three lineages. Demographic analyses showed that lineages B and C have experienced rapid demographic expansion at 0.03 Ma and 0.08 Ma, respectively.

CONCLUSIONS

Hemiculter leucisculus populations in drainages in southern China comprise three mtDNA lineages, and each of which may represent a separate species. Intense uplift of the Qinghai-Tibetan Plateau, evolution of Asian monsoons, changes in paleo-drainages, and poor dispersal ability may have driven the divergence of the three putative species. However, gene flow occurs among the three lineages. Climatic fluctuations have a prominent impact on the populations from the lineages B and C, but exerted little influence on the lineage A.

Phylogenetic analysis of the king cobra, Ophiophagus hannah in Thailand based on mitochondrial DNA sequences

Background: Snakes possess adaptive characteristics of morphology that may result in incorrect reconstruction of phylogeny. Molecular approaches have become the major source of new information for advancing our understanding of evolutionary, genetic relationships, and species identification.

Objective: We studied DNA sequences of Ophiophagus hannah in different parts of Thailand and compared them with those of O. hannah from other countries.

Materials and Methods: We studied the mitochondrial NADH dehydrogenase subunit 2 (ND2) and control region (CR) sequences of 12 individuals O. hannah from different localities across Thailand. Phylogenetic analysis of DNA sequences were compared to the published sequences of O. hannah deposited in NCBI GenBank database from other countries.

Results: O. hannah could be categorized into 2 Clades, 5 haplotypes, and 4 localities based on 43 different nucleotide positions from the 736 bp of ND2 and 673 bp of CR. Clade A was haplotype A from southern Thailand. Clade B consisted of haplotypes B, C, D, and E. Haplotype B and C came from northern Thailand. Haplotype D came from western Thailand, while haplotype E was O. hannah from central Thailand. The DNA sequences of Clade B were similar to the sequences of O. hannah in Myanmar and southern China that are already deposited in NCBI GenBank database.

Conclusion: We found a different genotype of O. hannah from southern Thailand and suggest that this may be a new species of O. hannah.

Huangshan population of Chinese Zacco platypus (Teleostei, Cyprinidae) harbors diverse matrilines and high genetic diversity

Six main mitochondrial DNA (mtDNA) lineages have been described in minnow (Zacco platypus) samples obtained from northern, western and southern China. Perdices et al. (2004) predicted that further sampling of other tributaries might discover more lineages of this species. In this study, we collected 26 Zacco platypus individuals in the Huangshan area of eastern China and determined the cytochrome b (cytb) sequence variations. Combined with reported data in GenBank, we identified ten matrilines (Zacco A-J) in a total of 169 samples, with relatively high molecular divergence found among them. The Huangshan population had the greatest genetic variation among all sampled regions and hosted six of the ten matrilines. Our results highlight the significance of the Huangshan area for the conservation of Zacco platypus.

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri, Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.

The phylogeographical pattern and conservation of the Chinese cobra (Naja atra) across its range based on mitochondrial control region sequences

The vulnerable Chinese cobra (Naja atra) ranges from southeastern China south of the Yangtze River to northern Vietnam and Laos. Large mountain ranges and water bodies may influence the pattern of genetic diversity of this species. We sequenced the mitochondrial DNA control region (1029 bp) using 285 individuals collected from 23 localities across the species' range and obtained 18 sequences unique to Taiwan from GenBank for phylogenetic and population analysis. Two distinct clades were identified, one including haplotypes from the two westernmost localities (Hekou and Miyi) and the other including haplotypes from all sampling sites except Miyi. A strong population structure was found (Φst = 0.76, P<0.0001) with high haplotype diversity (h = 1.00) and low nucleotide diversity (π = 0.0049). The Luoxiao and Nanling Mountains act as historical geographical barriers limiting gene exchange. In the haplotype network there were two “star” clusters. Haplotypes from populations east of the Luoxiao Mountains were represented within one cluster and haplotypes from populations west of the mountain range within the other, with haplotypes from populations south of the Nanling Mountains in between. Lineage sorting between mainland and island populations is incomplete. It remains unknown as to how much adaptive differentiation there is between population groups or within each group. We caution against long-distance transfers within any group, especially when environmental differences are apparent.

Phylogeography of the Rickett's big-footed bat, Myotis pilosus (Chiroptera: Vespertilionidae): a novel pattern of genetic structure of bats in China

Background

China is characterized by complex topographic structure and dramatic palaeoclimatic changes, making species biogeography studies particularly interesting. Previous researchers have also demonstrated multiple species experienced complex population histories, meanwhile multiple shelters existed in Chinese mainland. Despite this, species phylogeography is still largely unexplored. In the present study, we used a combination of microsatellites and mitochondrial DNA (mtDNA) to investigate the phylogeography of the east Asian fish-eating bat (Myotis pilosus).

Results

Phylogenetic analyses showed that M. pilosus comprised three main lineages: A, B and C, which corresponded to distinct geographic populations of the Yangtze Plain (YTP), Sichuan Basin (SCB) and North and South of China (NSC), respectively. The most recent common ancestor of M. pilosus was dated as 0.25 million years before present (BP). Population expansion events were inferred for populations of Clade C, North China Plain region, Clade B and YunGui Plateau region at 38,700, 15,900, 4,520 and 4,520 years BP, respectively. Conflicting results were obtained from mtDNA and microsatellite analyses; strong population genetic structure was obtained from mtDNA data but not microsatellite data. The microsatellite data indicated that genetic subdivision fits an isolation-by-distance (IBD) model, but the mtDNA data failed to support this model.

Conclusions

Our results suggested that Pleistocene climatic oscillations might have had a profound influence on the demographic history of M. pilosus. Spatial genetic structures of maternal lineages that are different from those observed in other sympatric bats species may be as a result of interactions among special population history and local environmental factors. There are at least three possible refugia for M. pilosus during glacial episodes. Apparently contradictory genetic structure patterns of mtDNA and microsatellite could be explained by male-mediated gene flow among populations. This study also provides insights on the necessity of conservation of M. pilosus populations to conserve this genetic biodiversity, especially in the areas of YTP, SCB and NSC regions.

Ongoing speciation in the Tibetan plateau Gymnocypris species complex

Local adaptation towards divergent ecological conditions often results in genetic differentiation and adaptive phenotypic divergence. To illuminate the ecological distinctiveness of the schizothoracine fish, we studied a Gymnocypris species complex consisting of three morphs distributed across four bodies of water (the Yellow River, Lake Qinghai, the Ganzi River and Lake Keluke) in the Northeast Tibetan Plateau. We used a combination of mitochondrial (16S rRNA and Cyt b) and nuclear (RAG-2) genetic sequences to investigate the phylogeography of these morphs based on a sample of 277 specimens. Analysis of gill rakers allowed for mapping of phenotypic trajectories along the phylogeny. The phylogenetic and morphological analyses showed that the three sparsely rakered morphs were present at two extremes of the phylogenetic tree: the Yellow River morphs were located at the basal phylogenetic split, and the Lake Keluke and Ganzi River morphs at the peak, with the densely rakered Lake Qinghai morphs located between these two extremes. Age estimation further indicated that the sparsely rakered morphs constituted the oldest and youngest lineages, whereas the densely rakered morph was assigned to an intermediate-age lineage. These results are most compatible with the process of evolutionary convergence or reversal. Disruptive natural selection due to divergent habitats and dietary preferences is likely the driving force behind the formation of new morphs, and the similarities between their phenotypes may be attributable to the similarities between their forms of niche tracking associated with food acquisition. This study provides the first genetic evidence for the occurrence of convergence or reversal in the schizothoracine fish of the Tibetan Plateau at small temporal scales.

Genealogy and Demographic History of a Widespread Amphibian throughout Indochina

Relatively little is known about spatial patterns of cryptic diversity in tropical species and the processes that generate them. Few studies examine the geographic distribution of genetic lineages in Southeast Asia, an area hypothesized to harbor substantial cryptic diversity. We investigated the evolutionary history of Asian tree frogs of the Polypedates leucomystax complex (n = 172) based on 1800 bp of the mtDNA genes ND1 and cytochrome b and tested hypotheses pertaining to climate, geology, and dispersal patterns. Analyses revealed substantial genetic diversity and lineage divergence throughout the region with evidence for widespread sympatric lineages and a general north versus south clustering. Relaxed molecular clock analysis and tests for demographic expansion identified an initial cladogenesis during the Miocene with subsequent Plio-Pleistocene diversification, with the former corresponding to periods of increased aridity and the onset of monsoonal weather systems. Rates of diversification were relatively constant until the Early Pleistocene when rates increased exponentially. We found equivocal evidence for isolation-by-distance and a potential role of some landscape features as partial barriers to dispersal. Finally, our analyses showed that divergence between insular and mainland populations occurred before Homo sapiens colonized Southeast Asia, suggesting that historical human-mediated dispersal did not drive insular diversification. Our results suggested that demographic expansion in the Late Pleistocene resulted in widespread sympatric lineages in the P. leucomystax complex throughout southern China and Indochina and further clarified the evolutionary history of lineages within P. leucomystax.

Significance of pre-Quaternary climate change for montane species diversity: insights from Asian salamanders (Salamandridae: Pachytriton)

Despite extensive focus on the genetic legacy of Pleistocene glaciation, impacts of earlier climatic change on biodiversity are poorly understood. Because amphibians are highly sensitive to variations in precipitation and temperature, we use a genus of Chinese montane salamanders (Salamandridae: Pachytriton) to study paleoclimatic change in East Asia, which experienced intensification of its monsoon circulation in the late Miocene associated with subsequent Pliocene warming. Using both nuclear and mitochondrial DNA sequences, we reconstruct the species tree under a coalescent model and demonstrate that all major lineages originated before the Quaternary. Initial speciation within the genus occurred after the summer monsoon entered a stage of substantial intensification. Heavy summer precipitation established temporary water connectivity through overflows between adjacent stream systems, which may facilitate geographic range expansion by aquatic species such as Pachytriton. Species were formed in allopatry likely through vicariant isolation during or after range expansion. To evaluate the influence of Pliocene warming on these cold-adapted salamanders, we construct a novel temperature buffer-zone model, which suggests widespread physiological stress or even extinction during the warming period. A significant deceleration of species accumulation rate is consistent with Pliocene range contraction, which affected P. granulosus and P. archospotus the most because they lack large temperature buffer zones. In contrast, demographic growth occurred in species for which refugia persist. The buffer-zone model reveals the Huangshan Mountain as a potential climatic refugium, which is similar to that found for other East Asian organisms. Our approach can incorporate future climatic data to evaluate the potential impact of ongoing global warming on montane species (particularly amphibians) and to predict possible population declines.

Genetic structure and demographic history should inform conservation: Chinese cobras currently treated as homogenous show population divergence

An understanding of population structure and genetic diversity is crucial for wildlife conservation and for determining the integrity of wildlife populations. The vulnerable Chinese cobra (Naja atra) has a distribution from the mouth of the Yangtze River down to northern Vietnam and Laos, within which several large mountain ranges and water bodies may influence population structure. We combined 12 microsatellite loci and 1117 bp of the mitochondrial cytochrome b gene to explore genetic structure and demographic history in this species, using 269 individuals from various localities in Mainland China and Vietnam. High levels of genetic variation were identified for both mtDNA and microsatellites. mtDNA data revealed two main (Vietnam + southern China + southwestern China; eastern + southeastern China) and one minor (comprising only two individuals from the westernmost site) clades. Microsatellite data divided the eastern + southeastern China clade further into two genetic clusters, which include individuals from the eastern and southeastern regions, respectively. The Luoxiao and Nanling Mountains may be important barriers affecting the diversification of lineages. In the haplotype network of cytchrome b, many haplotypes were represented within a “star” cluster and this and other tests suggest recent expansion. However, microsatellite analyses did not yield strong evidence for a recent bottleneck for any population or genetic cluster. The three main clusters identified here should be considered as independent management units for conservation purposes. The release of Chinese cobras into the wild should cease unless their origin can be determined, and this will avoid problems arising from unnatural homogenization.

Population genetic structure and phylogeographical pattern of rice grasshopper, Oxya hyla intricata, across Southeast Asia

The rice grasshopper, Oxya hyla intricata, is a rice pest in Southeast Asia. In this study, population genetic diversity and structure of this Oxya species was examined using both DNA sequences and AFLP technology. The samples of 12 populations were collected from four Southeast Asian countries, among which 175 individuals were analysed using mitochondrial DNA cytochrome c oxidase subunit I (COI) sequences, and 232 individuals were examined using amplified fragment length polymorphisms (AFLP) to test whether the phylogeographical pattern and population genetics of this species are related to past geological events and/or climatic oscillations. No obvious trend of genetic diversity was found along a latitude/longitude gradient among different geographical groups. Phylogenetic analysis indicated three deep monophyletic clades that approximately correspond to three geographical regions separated by high mountains and a deep strait, and TCS analysis also revealed three disconnected networks, suggesting that spatial and temporal separations by vicariance, which were also supported by AMOVA as a source of the molecular variance presented among groups. Gene flow analysis showed that there had been frequent historical gene flow among local populations in different regions, but the networks exhibited no shared haplotype among populations. In conclusion, the past geological events and climatic fluctuations are the most important factor on the phylogeographical structure and genetic patterns of O. hyla intricata in Southeast Asia. Habitat, vegetation, and anthropogenic effect may also contribute to gene flow and introgression of this species. Moreover, temperature, abundant rainfall and a diversity of graminaceous species are beneficial for the migration of O. hyla intricata. High haplotype diversity, deep phylogenetic division, negative Fu's F (s) values and unimodal and multimodal distribution shapes all suggest a complicated demographic expansion pattern of these O. hyla intricata populations, which might have been caused by climatic oscillations during glacial periods in the Quaternary.

Phylogeography of the Alcippe morrisonia (Aves: Timaliidae): long population history beyond late Pleistocene glaciations

Background

The role of Pleistocene glacial oscillations in current biodiversity and distribution patterns varies with latitude, physical topology and population life history and has long been a topic of discussion. However, there had been little phylogeographical research in south China, where the geophysical complexity is associated with great biodiversity. A bird endemic in Southeast Asia, the Grey-cheeked Fulvetta, Alcippe morrisonia, has been reported to show deep genetic divergences among its seven subspecies. In the present study, we investigated the phylogeography of A. morrisonia to explore its population structure and evolutionary history, in order to gain insight into the effect of geological events on the speciation and diversity of birds endemic in south China.

Results

Mitochondrial genes cytochrome b (Cytb) and cytochrome c oxidase I (COI) were represented by 1236 nucleotide sites from 151 individuals from 29 localities. Phylogenetic analysis showed seven monophyletic clades congruent with the geographically separated groups, which were identified as major sources of molecular variance (90.92%) by AMOVA. TCS analysis revealed four disconnected networks, and that no haplotype was shared among the geographical groups. The common ancestor of these populations was dated to 11.6 Mya and several divergence events were estimated along the population evolutionary history. Isolation by distance was inferred by NCPA to be responsible for the current intra-population genetic pattern and gene flow among geographical groups was interrupted. A late Pleistocene demographic expansion was detected in the eastern geographical groups, while the expansion time (0.2–0.4 Mya) was earlier than the Last Glacial Maximum.

Conclusion

It is proposed that the complicated topology preserves high genetic diversity and ancient lineages for geographical groups of A. morrisonia in China mainland and its two major islands, and restricts gene exchange during climate oscillations. Isolation by distance seems to be an important factor of genetic structure formation within geographical populations. Although glacial influence to population fluctuation was observed in late Pleistocene, it seems that populations in eastern China were more susceptible to climate change, and all geographical groups were growing stably through the Last Glacial Maximum. Coalescence analysis suggested that the ancestor of A. morrisonia might be traced back to the late Miocene, and the current phylogeographical structure of A. morrisonia is more likely to be attributable to a series geological events than to Pleistocene glacial cycles.

Pleistocene glacial refugia across the Appalachian Mountains and coastal plain in the millipede genus Narceus: evidence from population genetic, phylogeographic, and paleoclimatic data

Background

Species that are widespread throughout historically glaciated and currently non-glaciated areas provide excellent opportunities to investigate the role of Pleistocene climatic change on the distribution of North American biodiversity. Many studies indicate that northern animal populations exhibit low levels of genetic diversity over geographically widespread areas whereas southern populations exhibit relatively high levels. Recently, paleoclimatic data have been combined with niche-based distribution modeling to locate possible refugia during the Last Glacial Maximum. Using phylogeographic, population, and paleoclimatic data, we show that the distribution and mitochondrial data for the millipede genus Narceus are consistent with classical examples of Pleistocene refugia and subsequent post-glacial population expansion seen in other organismal groups.

Results

The phylogeographic structure of Narceus reveals a complex evolutionary history with signatures of multiple refugia in southeastern North America followed by two major northern expansions. Evidence for refugial populations were found in the southern Appalachian Mountains and in the coastal plain. The northern expansions appear to have radiated from two separate refugia, one from the Gulf Coastal Plain area and the other from the mid-Atlantic coastal region. Distributional models of Narceus during the Last Glacial Maximum show a dramatic reduction from the current distribution, with suitable ecological zones concentrated along the Gulf and Atlantic coastal plain. We found a strong correlation between these zones of ecological suitability inferred from our paleo-model with levels of genetic diversity derived from phylogenetic and population estimates of genetic structuring.

Conclusion

The signature of climatic change, during and after the Pleistocene, on the distribution of the millipede genus Narceus is evident in the genetic data presented. Niche-based historical distribution modeling strengthens the conclusions drawn from the genetic data and proves useful in identifying probable refugia. Such interdisciplinary biogeographic studies provide a comprehensive approach to understanding these processes that generate and maintain biodiversity as well as the framework necessary to explore questions regarding evolutionary diversification of taxa.