Phylogeographic structure of Brachymystax lenok tsinlingensis (Salmonidae) populations in the Qinling Mountains, Shaanxi, based on mtDNA control region

Genetic Structure and Phylogeographic Divergence of Thymallus brevicephalus in the Ob-Irtysh River Headwaters

Clarifying the genetic structure and population history of a species can reveal the impacts of historical climate and geological changes, providing critical insights for developing effective conservation strategies for ecologically significant fish. The Markakol grayling (Thymallus brevicephalus), an endangered species found in the Altai-Sayan Mountain region of Central Asia, serves as an ideal model for studying these factors. In this study, populations of a grayling (Thymallus) species discovered in the upper Irtysh River headwaters in Xinjiang, China, were analyzed to assess genetic diversity and population structure. Mitochondrial DNA sequences (cytochrome b and control region), along with 10 microsatellite markers, were used to examine genetic variation. Phylogenetic and genetic distance analyses confirmed the species, long misidentified as Arctic grayling (T. arcticus), as T. brevicephalus. This species can be divided into two distinct geographic groups: eastern and western, with the Crane River acting as the boundary. The divergence between these groups likely corresponds to refugia formed during the Pleistocene glaciation of the Altai Mountains, approximately 0.48 MA (million years ago) (range: 0.30 to 0.71 Ma). High haplotype diversity (Hd > 0.5) and low nucleotide diversity (π < 0.005) suggest that, despite the species' genetic richness, T. brevicephalus remains vulnerable to genetic drift, which could threaten its long-term survival. This vulnerability may stem from inbreeding within small refugial populations during the glacial period, followed by gradual population expansion. Our study offers novel insights into grayling populations, with results that have direct implications for management by serving as a tool for the identification of conservation units.

The phylogeographical pattern of the Amur minnow Rhynchocypris lagowskii (Cypriniformes: Cyprinidae) in the Qinling Mountains

In this study, the phylogeographical pattern of the Amur minnow (Rhynchocypris lagowskii) widely distributed in the cold freshwaters of the Qinling Mountains was examined. A total of 464 specimens from 48 localities were sequenced at a 540-bp region of the mitochondrial cytochrome b (Cytb) gene, and 69 haplotypes were obtained. The mean ratio of the number of synonymous and nonsynonymous substitutions per site (dN/dS) was 0.028 and indicated purifying selection. Haplotype diversity (h) and nucleotide diversity (π) of natural populations of R. lagowskii varied widely between distinct localities. Phylogenetic trees based on Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and network analysis showed five well-differentiated lineages, but these did not completely correspond to localities and geographic distribution. Meanwhile, analysis of molecular variances (AMOVA) indicated the highest proportion of genetic variation was attributed to the differentiation between populations rather than by our defined lineages. In addition, there was no significant correlation between the pairwise Fst values and geographic distance (p > .05). Based on the molecular clock calibration, the time to the most recent common ancestor (TMRCA) was estimated to have emerged from the Late Miocene to the Early Pleistocene. Finally, the results of demographic history based on the neutrality test, mismatch distribution, and Bayesian skyline plot (BSP) analyses showed that collectively, the populations were stable during the Pleistocene while one lineage (lineage E) probably underwent a slight contraction during the Middle Pleistocene and a rapid expansion from the Middle to the Late Pleistocene. Therefore, the study suggests the current phylogeographical pattern of R. lagowskii was likely shaped by geological events that led to vicariance followed by dispersal and secondary contact, river capture, and climatic oscillation during the Late Miocene to the Early Pleistocene in the Qinling Mountains.

Mitochondrial DNA revealed the validation of Quasipaa robertingeri (Amphibia: Anura: Dicroglossidae) and its population genetic diversity

The spiny frog Quasipaa robertingeri is endemic to a narrow region of southwest China and its taxonomic validation is still controversial. Based on COI gene sequences of 110 individuals from seven populations of Q. robertingeri and its related species, we investigated the phylogenetic position and population genetic structure of the species. Phylogenetic analyses indicated that Q. robertingeri was deeply genetically separated from its closely related species Q. boulengeri, indicating the validation of the species. All samples of Q. robertingeri were clustered into two divergent lineages. Haplotype network, AMOVA, and genetic distance estimations also supported the separations of the two groups. Neutrality tests indicated that one lineage has been likely independently experienced a recent population expansion, leading to a secondary contact area between the two lineages.

Transcriptomic responses to heat stress in gill and liver of endangered Brachymystax lenok tsinlingensis

Global warming significantly affects fish, particularly cold-water fish, because increased temperature adversely impacts their abilities to grow or reproduce, and eventually influences their fitness or even causes death. To survive, fish may alter their distribution or behavior to avoid the stress, and perhaps acclimate or evolve resistance to the elevated temperature. Brachymystax lenok tsinlingensis is an endangered cold-water species in China, and it has been found to alter the altitudinal distribution, decrease swimming efficiency and develop resistance under heat exposure, which badly impact the continuing conservation work. To better protect them, it is essential to understand how they respond to thermal stress behaviorally and physiologically. Therefore, the fish were exposed to 24.5 °C and based on the time taken for them to lose equilibrium, they were separately sampled as sensitive and tolerant groups. Both gill and liver tissues were collected from both groups for transcriptome sequencing. Sequencing results demonstrated that control and tolerant groups were similar in transcriptomic patterns and sensitive groups differentially expressed more genes than tolerant ones, suggesting the gene expression of tolerant groups may return to base levels as exposure time increased. Tissue differences were the major factor affecting gene expression, and they also displayed different physiological responses to heat stress. Consistent with other studies, heat shock response, immune response, metabolic adjustment and ion transport were found to be triggered after exposed to elevated temperature. The findings would contribute to a better understanding of responding mechanisms of fish to thermal stress and provide guidance for future conservation programs.

De novo assembly and microsatellite marker development of the transcriptome of the endangered Brachymystax lenok tsinlingensis

BACKGROUND

Brachymystax lenok tsinlingensis is an endemic freshwater fish in Northeast Asia, but experienced a dramatic population decline due to over-exploitation, deteriorated habitats and global climate change. It has been listed as a threatened or endangered species in South Korea and China, respectively. However, the conservation and restoration work in wild B. lenok tsinlingensis populations require large amount of genetic and molecular data to support effective management of genetic resources, while the corresponding information is very limited.

OBJECTIVE

This study was conducted to generate transcriptome assembly and annotation, as well as to develop novel microsatellite markers for B. lenok tsinlingensis.

METHODS

We collected gill and liver tissues and performed transcriptome sequencing. Then the first transcriptome for B. lenok tsinlingensis was de novo assembled and annotated. Microsatellite markers were searched in the assembled transcripts and characterized within ninety individuals collected from three natural sites.

RESULTS

A total of 110,712 protein-coding transcripts were assembled, of which 82,861 transcripts were successfully annotated. This assembly displayed a high level of completeness with retrieving 94% of the single-copy orthologs conserved across vertebrate species. Furthermore, 75,891 microsatellite loci were identified from this transcriptome assembly and 20 polymorphic markers were randomly selected for characterization.

CONCLUSIONS

The microsatellite markers and the first transcriptome assembly would provide valuable resources for investigating genetic diversity and phylogeographic structure of wild populations and molecular mechanisms responding to stressful environments (e.g. increased water temperature) to guide future conservation studies and breeding programs.

Phylogeography and demographic history of Gyrodactylus konovalovi (Monogenoidea: Gyrodactylidae), an ectoparasite on the East Asia Amur minnow (Cyprinidae) in Central China

Gyrodactylus konovalovi is an ectoparasite on the Amur minnow (Rhynchocypris lagowskii) that is widely distributed in the cold fresh waters of East Asia. In the present study, the phylogeography and demographic history of G. konovalovi and the distribution of its host in the Qinling Mountains are examined. A total of 79 individual parasites was sequenced for a 528 bp region of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene, and 25 haplotypes were obtained. The substitution rate (dN/dS) was 0.068 and indicated purifying selection. Haplotype diversity (h) and nucleotide diversity (π) varied widely in the Qinling Mountains. Phylogenetic trees based on Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods and network analysis revealed that all haplotypes were consistently well-supported in three different lineages, indicating a significant geographic distribution pattern. There was a significant positive correlation between genetic differentiation (F st) and geographic distance. The results of mismatch distribution, neutrality test and Bayesian skyline plot analyses showed that whole populations underwent population contraction during the Pleistocene. Based on the molecular clock calibration, the most common ancestor was estimated to have emerged in the middle Pleistocene. Our study suggests for the first time that a clearly phylogeography of G. konovalovi was shaped by geological events and climate fluctuations, such as orogenesis, drainage capture changes, and vicariance, during the Pleistocene in the Qinling Mountains.

Population genetic diversity of an odorous frog Odorrana grahami (Amphibia: Anura: Ranidae) in relation to conservation based on mitochondrial DNA

The diskless-fingered odorous frog Odorrana grahami is widely distributed in the high-altitude mountains in the southwestern China and northern Indochina regions. In this study, a comparative analysis of the mitochondrial COI gene sequences was performed to examine the population genetic diversity of 76 individuals from 10 localities across the distributional range of the species. Haplotype diversity and nucleotide diversity were 0.605 and 0.00199, respectively, in the total population. An AMOVA indicated that 47.3% of the total variation originated from variation within individual populations and 52.7% came from variation between the 10 geographic populations. Tests of neutral evolution indicated that a recent expansion occurred in total population. The findings provide useful information for the conservation of this species.

Complete Mitochondrial Genomes of the Cherskii's Sculpin Cottus czerskii and Siberian Taimen Hucho taimen Reveal GenBank Entry Errors: Incorrect Species Identification and Recombinant Mitochondrial Genome

The complete mitochondrial (mt) genome is sequenced in 2 individuals of the Cherskii's sculpin Cottus czerskii. A surprisingly high level of sequence divergence (10.3%) has been detected between the 2 genomes of C czerskii studied here and the GenBank mt genome of C czerskii (KJ956027). At the same time, a surprisingly low level of divergence (1.4%) has been detected between the GenBank C czerskii (KJ956027) and the Amur sculpin Cottus szanaga (KX762049, KX762050). We argue that the observed discrepancies are due to incorrect taxonomic identification so that the GenBank accession number KJ956027 represents actually the mt genome of C szanaga erroneously identified as C czerskii. Our results are of consequence concerning the GenBank database quality, highlighting the potential negative consequences of entry errors, which once they are introduced tend to be propagated among databases and subsequent publications. We illustrate the premise with the data on recombinant mt genome of the Siberian taimen Hucho taimen (NCBI Reference Sequence Database NC_016426.1; GenBank accession number HQ897271.1), bearing 2 introgressed fragments (≈0.9 kb [kilobase]) from 2 lenok subspecies, Brachymystax lenok and Brachymystax lenok tsinlingensis, submitted to GenBank on June 12, 2011. Since the time of submission, the H taimen recombinant mt genome leading to incorrect phylogenetic inferences was propagated in multiple subsequent publications despite the fact that nonrecombinant H taimen genomes were also available (submitted to GenBank on August 2, 2014; KJ711549, KJ711550). Other examples of recombinant sequences persisting in GenBank are also considered. A GenBank Entry Error Depositary is urgently needed to monitor and avoid a progressive accumulation of wrong biological information.

Effects of habitat fragmentation on the population genetic diversity of the Amur minnow (Phoxinus lagowskii)

Phoxinus lagowskii is a freshwater fish that is widely distributed in China. In this study, a comparative analysis of the mtDNA control region (D-loop) was performed to analyze the natural population structure and genetic diversity of 54 individuals from eleven locations (T1, T2, T3, T4, T5, T6, T7, T8, T9, T10 and T11) which was divided with reservoirs. The estimated haplotype and nucleotide diversity were 0.734 and 0.03514, respectively. An AMOVA indicated that 79.78% of the total variation originated from individual populations and 20.22% came from variation within the 11 geographic populations, which showed high genetic differentiation among the 11 geographic groups. A test of neutral evolution and mismatch distribution indicated that historical expansion occurred in these populations. However, the findings of low genetic diversity and high genetic differentiation demonstrated that the reproduction isolated by reservoir has showed a certain effect for the development of the populations, and the results should provide new information for the conservation and exploitation of this species.

Phylogeography of chukar partridge (Alectoris chukar) in China based on mtDNA control region

The chukar partridge (Alectoris chukar) is distributed in north and northwest of China, in mountainous areas that were heavily affected by cyclic climate and landscape changes during the last Pleistocene glaciations. Some partridge populations have colonized and expanded their present ranges only after deglaciation and recent deforestation by human. Consequently, partridges from different areas could be genetically differentiated. In this study, a 1152-1154 bp portion of the mitochondrial DNA control region were analyzed for all 279 specimens collected from 28 populations through their distribution in China. A total number of 91 haplotypes were defined by 113 variable sites. The mean haplotype diversity and nucleotide diversity were 0.939 ± 0.008 and 0.0030 ± 0.0017, respectively, for chukar. Haplotype diversity among the 28 populations varied from 0.600 ± 0.215 (HH) to 1.000 ± 0.272 (SBC), and nucleotide diversity ranged from 0.0006 ± 0.0007 (HJ) to (0.0071 ± 0.0041) (HEG). The nucleotide diversity of chukar was descending from west to east and a possible historical scenario might be that A. chukar inhabited central Asia, and then dispersed eastward. Phylogenetic analyses indicated that there was no distinct phylogeographic structure in chukar populations in China. Haplotype network of chukar was star-like with some common haplotypes shared by different samples came from different populations as center. Both unstructured phylogeographic tree and star-like haplotype network are signatures of population expansion.