Influence of increasing acclimation temperature on growth, digestion, antioxidant capacity, liver transcriptome and intestinal microflora of Ussruri whitefish Coregonus ussuriensis Berg.

For effective restoration, conservation of Ussruri whitefish Coregonus ussuriensis Berg and coping with global climate change, effects of environmental temperature on Ussruri whitefish urgently need to be explored. In current study, the effects of different acclimation temperatures on the growth, digestive physiology, antioxidant ability, liver transcriptional responses and intestinal microflora patterns of Ussruri whitefish were investigated. Ussruri whitefish (15.20 g ± 1.23 g) were reared for 42 days under different acclimation temperatures, i.e., 10, 13, 16, 19, 22 and 25 °C, respectively. Result first determined 28 °C as the semi-lethal temperature in order to design the temperature gradient test. Highest main gain rate (MGR) and specific growth rate (SGR) were observed in fish group having acclimation temperature of 19 °C. Significantly decrease (P < 0.05) in triglyceride (TG) content appeared at 19 °C as compared to the 10 °C and 13 °C temperature groups. 19 °C notablely increased protease activities of stomach and intestine and intestinal lipase and amylase activities. 19 °C group obtained the highest activities of chloramphnicol acetyltransferase (CAT) and total antioxidant capacity (T-AOC) and higher activities of superoxide dismutase (SOD). The intestinal microflora composition was most conducive to maintaining overall intestinal health when the temperature was 19 °C, compared to 10 °C and 25 °C. Ussruri whitefish exposed to 10 °C and 25 °C possessed the lower Lactobacillus abundance compared to exposure to 19 °C. Temperature down to 10 °C or up to 25 °C, respectively, triggered cold stress and heat stress, which leading to impairment in intestinal digestion, liver antioxidant capacity and intestinal microflora structure. Liver transcriptome response to 10 °C, 19 °C and 25 °C revealed that Ussruri whitefish might require the initiation of endoplasmic reticulum stress to correct protein damage from cold-temperature and high-temperature stress, and it was speculated that DNAJB11 could be regarded as a biomarker of cold stress response.Based on the quadratic regression analysis of MGR and SGR against temperature, the optimal acclamation temperature were, respectively, 18.0 °C and 18.1 °C. Our findings provide valuable theoretical insights for an in-depth understanding of temperature acclimation mechanisms and laid the foundation for conservation and development of Ussruri whitefish germplasm resources.

Single and joint toxic effects of waterborne exposure to copper and cadmium on Coregonus ussuriensis Berg.

Heavy metal contamination severely affects the aquatic environment and organisms. Copper (Cu) and cadmium (Cd) are two of the most common heavy metal contaminants that impair the survival, development, and reproduction of aquatic organisms. With the growth of agriculture and industry, there is a possibility of heavy metal pollution in Coregonus ussuriensis Berg's water source. However, there are no published studies on the toxicity to C. ussuriensis. Acute toxicity experiments in C. ussuriensis revealed the 96-h median lethal concentrations of copper and cadmium to be 0.492 mg·L-1 (95% confidence interval: 0.452-0.529) and 1.548 mg·L-1 (95% confidence interval: 1.434-1.657), respectively, and safe concentrations of 4.92 µg·L-1 and 15.48 µg·L-1, respectively. C. ussuriensis was then treated for 96 h with Cu (20% of 96 h LC50), Cd (20% of 96 h LC50), and a combination of Cu and Cd (20% of Cu 96 h LC50 + 20% of Cd 96 h LC50). The histological damage caused by the three different exposure modes to the liver and gills of C. ussuriensis was verified using hematoxylin and eosin staining. All three exposure modes caused different degrees of vacuolization, nuclear consolidation, and necrosis in the liver tissue of C. ussuriensis and edema, hyperplasia, laminar fusion, and epithelial elevation in the gill tissue compared with the reference group. The severity of the damage increased with increasing exposure time. Anti-oxidant activity in the gill and liver tissues were measured using enzyme activity assay kits to reflect oxidative stress induced by copper and cadmium exposure alone and in combination. The enzyme activities of superoxide dismutase (SOD), catalase (CAT), and glutathione peroxidase (GSH) were substantially higher than those in the reference groups. However, the activities of the enzymes decreased with increasing exposure time. Malondialdehyde (MDA) activity significantly increased during exposure in relation to that in the reference group. Analysis of immune gene expression in C. ussuriensis gill and liver tissues was executed using real-time inverse transcript polymerase chain response (RT-PCR). The expression levels of the pro-inflammatory cytokines interleukin one beta (IL-1ß) and tumor necrosis factor-alpha (TNF-α) were positively correlated with exposure time and were significantly upregulated with increasing exposure time. Metallothionein (MT) gene expression levels were significantly upregulated in the short term after exposure compared to the reference group but decreased with increasing exposure time. Our results indicate that exposure to aqueous copper and cadmium solutions, either alone or in combination, causes histopathological damage, oxidative stress, and immunotoxicity in C. ussuriensis gill and liver tissue. This study investigated the toxic effects of copper and cadmium on C. ussuriensis to facilitate the monitoring of heavy metals in water sources for healthy aquaculture.

ARHGAP21 Is Involved in the Carcinogenic Mechanism of Cholangiocarcinoma: A Study Based on Bioinformatic Analyses and Experimental Validation.

Background and Objectives: Rho GTPase-activating protein (RhoGAP) is a negative regulatory element of Rho GTPases and participates in tumorigenesis. Rho GTPase-activating protein 21 (ARHGAP21) is one of the RhoGAPs and its role in cholangiocarcinoma (CCA) has never been disclosed in any publications. Materials and Methods: The bioinformatics public datasets were utilized to investigate the expression patterns and mutations of ARHGAP21 as well as its prognostic significance in CCA. The biological functions of ARHGAP21 in CCA cells (RBE and Hccc9810 cell) were evaluated by scratch assay, cell counting kit-8 assay (CCK8) assay, and transwell migration assay. In addition, the underlying mechanism of ARHGAP21 involved in CCA was investigated by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and the most significant signaling pathway was identified through gene set enrichment analysis (GSEA) and the Western blot method. The ssGSEA algorithm was further used to explore the immune-related mechanism of ARHGAP21 in CCA. Results: The ARHGAP21 expression in CCA tissue was higher than it was in normal tissue, and missense mutation was the main alteration of ARHGAP21 in CCA. Moreover, the expression of ARHGAP21 had obvious differences in patients with different clinical characteristics and it had great prognostic significance. Based on cell experiments, we further observed that the proliferation ability and migration ability of the ARHGAP21-knockdown group was reduced in CCA cells. Several pathological signaling pathways correlated with proliferation and migration were determined by GO and KEGG analysis. Furthermore, the PI3K/Akt signaling pathway was the most significant one. GSEA analysis further verified that ARHGAP21 was highly enriched in PI3K/Akt signaling pathway, and the results of Western blot suggested that the phosphorylated PI3K and Akt were decreased in the ARHGAP21-knockdown group. The drug susceptibility of the PI3K/Akt signaling pathway targeted drugs were positively correlated with ARHGAP21 expression. Moreover, we also discovered that ARHGAP21 was correlated with neutrophil, pDC, and mast cell infiltration as well as immune-related genes in CCA. Conclusions: ARHGAP21 could promote the proliferation and migration of CCA cells by activating the PI3K/Akt signaling pathway, and ARHGAP21 may participate in the immune modulating function of the tumor microenvironment.

Clinical application of traditional Chinese medicine powder in the treatment of acute and chronic wounds.

This study aimed to explore the clinical application and efficacy of traditional Chinese medicine (TCM) powder in the treatment of acute and chronic wounds. Eighty patients with a wound infection were randomly and equally divided into a control group and an observation group. Gauze padding containing furacilin was used to dress the infected wounds of the control group. TCM powder was used to treat the wounds of the observation group. The total response rate of the observation group was significantly higher than the control group (P = .017). The colour and exudate volume scores in the observation group were lower than the control group, and the differences between the two groups were statistically significant (P < .05). The time to the appearance of new epithelium and time to the wound healing of the burns in the observation group were shorter than the control group, and the differences were statistically significant (P < .05). The TCM powder absorbed a large amount of necrotic tissue and exudate from the wound surface, cleared heat and toxins, and activated blood circulation. It also resolved blood stasis, eliminated pus, and allowed for new skin growth, as well as regenerating muscle.

miR-301b-5p and its target gene nfatc2ip regulate inflammatory responses in the liver of rainbow trout (Oncorhynchus mykiss) under high temperature stress.

Rainbow trout (Oncorhynchus mykiss) is a typical cold-water aquaculture fish and a high-end aquatic product. When water temperature exceeds its optimal range of 12-18 °C, the immune system of rainbow trout becomes weakened and unbalanced. High temperature in summer and global warming severely impact rainbow trout industry. The focus of this study was to explore the mechanisms regulating the immune response of rainbow trout under high temperature stress and identify molecular elements that account for resistance to high temperature. In this study, individual fish were screened in a high temperature stress experiment and divided into resistant (R) and sensitive (S) groups. The hepatic transcriptome sequencing and analysis of mRNAs and microRNAs of the R, S, and control groups showed that the number of the differentially expressed genes (DEGs) in the S group (9259) was higher than that in the R group (5313). Furthermore, the 1233 genes differentially expressed between S and R groups were mainly enriched in immune-related pathways, including cytokine-cytokine receptor interaction, TNF signaling and IL-17 signaling. Among these DEGs were miR-301b-5p and its target gene that encodes nuclear factor of activated T cells two interacting protein (nfatc2ip). The dual-luciferase reporter system and immunofluorescence experiments verified the relationship between miR-301b-5p and nfatc2ip. We also showed that expression levels of miR-301b-5p and nfatc2ip significantly negatively correlated in the liver of rainbow trout under high temperature stress. By performing functional experiments, we showed that activation of miR-301b-5p expression or inhibition of nfatc2ip expression stimulated the phosphorylation of p65, p38, and JNK in the classical nuclear factor kappa-B and mitogen-activated protein kinase pathways under high temperature stress. These manipulations initially promoted the secretion of the pro-inflammatory factor IL-1ß and then increased the levels of IL-6, IL-12, and TNF-α. In addition, activation of miR-301b-5p expression or inhibition of nfatc2ip expression stimulated the repair of the hepatic ultrastructural damage caused by high temperature stress by activating the inflammatory response in rainbow trout liver.

DOK5 as a Prognostic Biomarker of Gastric Cancer Immunoinvasion: A Bioinformatics Analysis.

BACKGROUND: Docking protein 5 (DOK5) is a member of the docking protein group of membrane proteins and is an adapter protein involved in signal transduction. Nevertheless, the role of DOK5 expression in the prognosis of gastric cancer (GC) remains unclear. METHODS: In this study, clinical prognostic parameters and survival data related to DOK5, in patients with GC, were analyzed using bioinformatics analysis comprising Oncomine and TIMER, UALCAN database, Kaplan-Meier plotter, GEPIA, GSEA, DAVID, and cBioPortal websites. RESULTS: In our study, GC contained various DOK5 expressions, which forecasted poor survival outcomes. Moreover, our research showed that high DOK5 could predict high-level infiltration of several GC immune cells, as evidenced by M1, TAM, M2, B cell, and T cell failure. Hence, DOK5 might become a new gastric cancer biomarker and therapeutic target. In the following analysis, in order to explore the prognostic value of DOK5 in GC, more clinical trials are needed to validate our results. CONCLUSIONS: Through multiple database verifications, DOK5 was found to be part of the pathogenic genes for GC. Thus, it can change the formation and progression of tumors by acting on human immunity.

Molecular characterization and antibacterial immunity functional analysis of the antimicrobial peptide hepcidin from Coregonus ussuriensis berg.

Antimicrobial peptides are immune system molecules existing in different organisms including mollusks, crustaceans and vertebrates. Hepcidins are a group of cysteine rich antimicrobial peptides, which plays an important role in fish response to a variety of pathogens. In this study, we cloned and identified Hepcidin from the Coregonus ussuriensis Berg, and its functions in vivo and in vitro was investigated. Our results showed that, CuHepc contains a 267 bp coding sequence (CDS) region that encodes 88 putative amino acids with a molecular weight of 9.77 kD. Hepcidin transcripts were most abundant in the liver of healthy C. ussuriensis Berg. The synthesized Hepcidin peptide exhibited a wide range of antibacterial activity against Gram-positive and Gram-negative bacteria in vitro, and the results of in vivo bacterial attack assays showed that the CuHepc gene was differentially up-regulated in the six tissues investigated after infection with Aeromonas hydrophila. To analyze the changes in protein levels in C. ussuriensis, we generated Hepc polyclonal antibodies in rabbits and verified that the protein expression was increased after bacterial infection with Western blot assay. MIC assay results showed a geometric mean value of 5.513 µM for CuHepc peptide. In the in vivo experiment, immune-related genes IL-10, NF-κB, TLR3 were up-regulated post-infection CuHepc peptide in liver and intestine. Finally, CuHepc peptide reduced the tissues microbial load compared to infection with Aeromonas hydrophila. The above results indicate that Hepc plays a role in the immune response of C. ussuriensis to exogenous disturbances, indicate that CuHepc might act a candidate for modulation of the innate immune system in C. ussuriensis.

Genome resequencing reveals independent domestication and breeding improvement of naked oat.

As an important cereal crop, common oat, has attracted more and more attention due to its healthy nutritional components and bioactive compounds. Here, high-depth resequencing of 115 oat accessions and closely related hexaploid species worldwide was performed. Based on genetic diversity and linkage disequilibrium analysis, it was found that hulled oat (Avena sativa) experienced a more severe bottleneck than naked oat (Avena sativa var. nuda). Combined with the divergence time of ∼51,200 years ago, the previous speculation that naked oat was a variant of hulled oat was rejected. It was found that the common segments that hulled oat introgressed to naked oat cultivars contained 444 genes, mainly enriched in photosynthetic efficiency-related pathways. Selective sweeps during environmental adaptation and breeding improvement were identified in the naked oat genome. Candidate genes associated with smut resistance and the days to maturity phenotype were also identified. Our study provides genomic resources and new insights into naked oat domestication and breeding.

CDCA3 Is a Novel Prognostic Biomarker Associated with Immune Infiltration in Hepatocellular Carcinoma.

Cell division cycle-associated protein-3 (CDCA3) contributes to the regulation of the cell cycle. CDCA3 plays an important role in the carcinogenesis of various cancers; however, the association between CDCA3 expression, prognosis of patients, and immune infiltration in the tumor microenvironment is still unknown. Here, we demonstrated that CDCA3 was differentially expressed between the tumor tissues and corresponding normal tissues using in silico analysis in the ONCOMINE and Tumor Immune Estimation Resource (TIMER) databases. We analyzed the relationship between the expression of CDCA3 and prognosis of patients with hepatocellular carcinoma (HCC) using the Kaplan-Meier plotter database and Gene Expression Profiling Interactive Analysis (GEPIA). Furthermore, we determined the prognostic value of CDCA3 expression using univariate and multivariate analyses. We observed that CDCA3 expression closely correlated with immune infiltration and gene markers of infiltrating immune cells in the TIMER database. CDCA3 was highly expressed in the tumor tissues than in the adjacent normal tissues in various cancers, including HCC. Increased expression of CDCA3 was accompanied by poorer overall survival (OS), relapse-free survival (RFS), progression-free survival (PFS), and disease-specific survival (DSS). The correlation between CDCA3 expression and OS and disease-free survival (DFS) was also studied using GEPIA. CDCA3 expression was associated with the levels of immune cell infiltration and was positively correlated with tumor purity. Moreover, CDCA3 expression was associated with gene markers such as PD-1, CTLA4, LAG3, and TIM-3 from exhausted T cells, CD3D, CD3E, and CD2 from T cells, and TGFB1 and CCR8 located on the surface of Tregs. Thus, we demonstrated that CDCA3 may be a potential target and biomarker for the management and diagnosis of HCC.

A comparative analysis of the essential oils from two species of garlic seedlings cultivated in China: chemical profile and anticoagulant potential.

Garlic seedlings (GS) and blanched garlic seedlings (BGS) are two kinds of common garlic-derived vegetables in China, but little information is available on their bioactive constituents. In this work, chemical profiles and anticoagulant activities of essential oils from GS (EOGS) and BGS (EOBGS) were disclosed and compared for the first time. Sixteen and fourteen volatile compounds were identified in EOGS and EOBGS by GC-MS analysis, and both of them were rich in sulfur-containing compounds, particularly diallyl sulfides accounting for 74.77% and 85.87%, respectively. EOGS and EOBGS exerted anticoagulant activities via intrinsic, extrinsic, and common coagulation pathways as well as by lowering the content of fibrinogen; EOGS exceeded EOBGS in the activation of intrinsic and extrinsic coagulation pathways, while EOBGS outperformed EOGS on the activation of the common coagulation pathway, which was even superior to that of heparin at the same dose. Herein, the results of the present investigation will give a strong clue that EOGS and EOBGS are more likely to lead to a promising way to vegetable-based anticoagulants.

Preparation of self-assembled collagen fibrillar gel from tilapia skin and its formation in presence of acidic polysaccharides.

Fibrillar gel of pepsin-solubilized collagen from tilapia skin was prepared by self-assembly in neutral phosphate buffer at 28 °C. Then effects of acidic polysaccharides, such as sodium alginate (SA), chondroitin sulfate (CS), and hyaluronic acid (HA), on the formation and properties of self-assembled fibrillar gel were investigated. SA and CS prolonged gelling time, whereas HA had no obvious effect. SA made fibril network denser, while CS and HA induced the presence of larger ordered structures. All the acidic polysaccharides broadened the D-periodicity of fibrils. SA and HA increased the maximum mechanical strength of gel to 39.64 and 34.49 kN/m2, respectively, significantly higher than that of pure collagen gel (14.53 kN/m2), while that only 17.20 kN/m2 after CS introduced. HA had no evident effect on enzymatic resistance, while SA and CS decreased. Therefore, tilapia skin collagen with HA has a higher potential as a biomaterial than that with CS or SA.

The association between interleukin-1 polymorphisms and their protein expression in Chinese Han patients with breast cancer.

BACKGROUND: Breast cancer (BC) is the most common cancer in women and the second leading cause of cancer-related deaths among women worldwide. Single nucleotide polymorphisms (SNPs) in cytokine genes have been shown to alter their expressions or functions in patients with BC. In recent years, the molecular structure and function of IL-1 have been studied. Its genetic polymorphism could affect the transcription and expression of the IL-1 gene. Moreover, it is closely related to several diseases. This fact and plethora of gene polymorphism data prompted us to investigate the relationship between IL-1 polymorphisms and IL-1 protein expression in Chinese Han BC patients. METHOD: In total, 298 patients with BC and 287 healthy control women were studied. The genetic polymorphisms for IL-1 were analyzed by the MassARRAY sequencing method. Tumor markers and IL-1ß levels were measured by electrochemiluminescence and ELISA, respectively. All gene selection GRCh38 version. RESULTS: The rs1143623 (NC_000002.12:g.112838252C>G) (GC), rs16944 (NC_000002.12:g.112837290A>G)(AG), and rs10490571 (NC_000002.12:g.102100877C>T) (CC) SNPs were found to be significantly lower in the BC group than in the controls. The variant G/C genotype of rs1143623 was associated with a significantly increased risk for BC (OR = 2.34, p < 0.05). The alleles for rs16944 (A/G; OR = 3.15, p < 0.05) and rs10490571 (T/C; OR = 2.48, p < 0.05) were also significantly associated with BC. Moreover, the genotypes of rs1143623, rs16944, and rs10490571 were significantly correlated with serum IL-1ß levels and other tumor markers. CONCLUSION: Our data reveal the association between genetic polymorphisms of IL-1 and BC susceptibility in the Chinese Han population and indicates that IL-1 polymorphisms are closely associated with tumor markers and IL-1ß protein expression in BC patients.

Genome re-sequencing reveals the evolutionary history of peach fruit edibility.

Peach (Prunus persica) is an economically important fruit crop and a well-characterized model for studying Prunus species. Here we explore the evolutionary history of peach using a large-scale SNP data set generated from 58 high-coverage genomes of cultivated peach and closely related relatives, including 44 newly re-sequenced accessions and 14 accessions from a previous study. Our analyses suggest that peach originated about 2.47 Mya in southwest China in glacial refugia generated by the uplift of the Tibetan plateau. Our exploration of genomic selection signatures and demographic history supports the hypothesis that frugivore-mediated selection occurred several million years before the eventual human-mediated domestication of peach. We also identify a large set of SNPs and/or CNVs, and candidate genes associated with fruit texture, taste, size, and skin color, with implications for genomic-selection breeding in peach. Collectively, this study provides valuable information for understanding the evolution and domestication of perennial fruit tree crops.

Phylogenetic and molecular dating analysis of Taiwan Blue Pheasant (Lophura swinhoii).

The Swinhoe's Pheasant (Lophura swinhoii) is an endemic and most endangered species to Taiwan, China. It belongs to the genus Lophura, family Phasianidae. To further investigate the evolutionary history of L. swinhoii, we determined its complete mitochondrial genome and reconstructed a single, robust phylogenetic tree. Our results showed that L. swinhoii is clustered with Lophura nycthemera and forms a sister group of Lophura ignita. The genus Lophura is strongly supported as the sister taxon of the genus Crossoptilon. The molecular clock analysis showed that the genetic divergence of L. swinhoii occurred in 2.71 (1.31-4.22) Mya. The most common ancestor of L. swinhoii might have migrated from mainland of South East Asia to Taiwan Island by the land bridge at 2.71 Mya ago. Taiwan Island is separated from the mainland by the sea (Taiwan Strait) and formed a separate island at around 2.5 Mya because of the transgression and regression. Therefore, geographical isolation and climate change may accelerate the evolution of L. swinhoii. In this study, we propose a biogeographic hypothesis for speciation of the L. swinhoii based on known events of the geographic and geological history of South East Asia and southeast China, which would benefit the understanding of evolutionary history of L. swinhoii as well as other galliform birds.

The complete mitochondrial genome sequence analysis of Tibetan argali (Ovis ammon hodgsoni): implications of Tibetan argali and Gansu argali as the same subspecies.

The genus Ovis (Bovidae, Artiodactyla) includes six species, i.e. Ovis ammon, Ovis aries, Ovis canadensis, Ovis dalli, Ovis nivicola and Ovis vignei. Based on morphology, geographical location, habitat, etc., the species O. ammon is divided into nine subspecies. The near threatened Tibetan argali is distributed across the Tibetan Plateau and its peripheral mountains, and believed to be one of the O. ammon subspecies (O. a. hodgsoni). However, considering its morphological features and distributions, a question has been proposed by some researchers about the subspecies status of Tibetan argali. In this study, we employed complete mitochondrial DNA (mtDNA) to explore the phylogenetic relationship and population genetic structure of Tibetan argali. The results revealed that the nucleotide composition, gene arrangement and codon usage pattern of the mitochondrial genome of Tibetan argali are similar to those of other caprines. Phylogenetic analyses showed that Tibetan argali was clustered with O. ammon. Interestingly, five Tibetan argali individuals and one of the three Gansu argali (O. a. dalailamae) individuals were clustered in the same branch, which is a sister group to other two Gansu argali individuals. Together with morphological characteristics, our results suggested that Tibetan argali and Gansu argali may belong to the same subspecies (O. a. hodgsoni) of O. ammon, rather than two different subspecies.

Genetic variability of the tokay gecko based on mitochondrial and nuclear DNA.

With largely allopatric distribution, the black tokay and the red tokay are two distinct morphs of the subspecies Gekko gecko gecko. In consideration of their different morphological characteristics, the taxonomic status of G. g. gecko is disputed. Through detailed morphological comparison, it has been proposed that the black tokay should be elevated to species ranking, but without strong genetic evidence. In order to further investigate the taxonomic status of the tokay gecko (G. gecko), we used one mitochondrial marker (ND2) and three nuclear markers (RAG1, c-mos, and ITS2) to explore the phylogenetic and taxonomic relationship of the tokay gecko. Our results revealed a deep phylogeographical divergence in tokay gecko and at the same time provided us with the evidence of possible introgressive hybridization or/and incomplete lineage sorting between the black tokay and the red tokay. The elevation of the black tokay to species level is also supported by our results. However, due to limited sampling and genetic data, this elevation should be further corroborated by more genetic evidence.

Molecular evidence for the subspecific differentiation of blue sheep (Pseudois nayaur) and polyphyletic origin of dwarf blue sheep (Pseudois schaeferi).

Blue sheep (Pseudois nayaur), a Central Asian ungulate with restricted geographic distribution, exhibits unclear variation in morphology and phylogeographic structure. The composition of species and subspecies in the genus Pseudois is controversial, particularly with respect to the taxonomic designation of geographically restricted populations. Here, 26 specimens including 5 dwarf blue sheep (Pseudois schaeferi), which were collected from a broad geographic region in China, were analyzed for 2 mitochondrial DNA fragments (cytochrome b and control region sequences). In a pattern consistent with geographically defined subspecies, we found three deeply divergent mitochondrial lineages restricted to different geographic regions. The currently designated two subspecies of blue sheep, Pseudois nayaur nayaur and Pseudois nayaur szechuanensis, were recognized in the phylogenetic trees. In addition, the Helan Mountain population showed distinct genetic characteristics from other geographic populations, and thus should be classified as a new subspecies. In contrast, dwarf blue sheep clustered closely with some blue sheep from Sichuan Province in the phylogenetic trees. Therefore, dwarf blue sheep appear to be a subset of Pseudois nayaur szechuanensis. After considering both population genetic information and molecular clock analysis, we obtained some relevant molecular phylogeographic information concerning the historical biogeography of blue sheep. These results also indicate that western Sichuan was a potential refugium for blue sheep during the Quaternary period.