Maintaining the native gut microbiota of bharal (Pseudois nayaur) is crucial in ex situ conservation.

As wildlife protection continue to strengthen, research on the gut microbiota of wildlife is increasing. Carrying out conservation and research on endangered species in the Qinghai Tibet Plateau plays an important role in global biodiversity conservation. This study utilized 16S rRNA sequencing of fecal samples to investigate the composition, function, and changes of the gut microbiota of bharal in different environments, seasons, and genders. The results showed that Firmicutes and Bacteroidota were the dominant phyla and UCG-005, Bacteroides, UCG-010 were the dominant genera of bharal. In the wild, the abundance of Firmicutes increased which was conducive to the decomposition and utilization of cellulose, hemicellulose, and carbohydrate. Due to the variety of food types and nutrition in different seasons, the composition and function of gut microbiota were obviously different between genders. Compared with zoo, higher alpha diversity, a more complex gut microbiota network structure, and stronger metabolic function were conducive bharal to adapting to the wild environment. In the zoo, captive bharals were fed foods rich in high fat and protein, which increased the abundance of Bacteroidota and reduced the alpha diversity of gut microbiota. A fixed diet unified the gut microbiota between genders of bharal. It is very important to pay attention to the impact of captive environments and maintain the native gut microbiota of wildlife.

Gut microbiota non-convergence and adaptations in sympatric Tibetan and Przewalski's gazelles.

Unraveling the connection between gut microbiota and adaptability in wild species in natural habitats is imperative yet challenging. We studied the gut microbiota of sympatric and allopatric populations of two closely related species, the Procapra picticaudata and P. przewalskii, with the latter showing lower adaptability and adaptive potential than the former. Despite shared habitat, sympatric populations showed no convergence in gut microbiota, revealing distinct microbiota-environment relationships between the two gazelle species. Furthermore, the gut microbiota assembly process of the P. przewalskii was shifted toward homogeneous selection processes relative to that of the P. picticaudata. Those taxa which contributed to the shift were mainly from the phyla Firmicutes and Verrucomicrobiota, with functions highly related to micronutrient and macronutrient metabolism. Our study provides new insights into the complex dynamics between gut microbiota, host adaptability, and environment in wildlife adaptation and highlights the need to consider host adaptability when examining wildlife host-microbiome interplay.

Effects of ex situ conservation on diversity and function of the gut microbiota of the Tibetan wild ass (Equus kiang).

Ex situ conservation is the main method for the protection of endangered wildlife. To explore the effect of ex situ conservation on the gut microbiota of the kiang (Equus kiang), metagenomic sequencing combined with bioinformatics analysis was used to investigate the composition and function of the gut microbiota of the kiang. The results showed that ex situ conservation not only protected wildlife, but also affected the composition and function of gut microbiota, as well as the health of animals. In the zoo, the ratio of the relative abundance of Firmicutes to that of Bacteroidetes (F/B) is higher, clusters of potentially pathogenic bacteria (such as Catonella, Catonella, and Mycoplasma) are more numerous, the abundance of resistance genes is higher, and the abundance of metabolic functions is increased. The dynamic changes of the gut microbiota also played an important role in the nutritional absorption, energy metabolism, and environmental adaptation of the kiang. Improving the rearing environment and increasing food diversity play important roles for increasing the diversity of gut microbiota, reducing the spread of potentially pathogenic bacteria, and reducing diseases. In the wild, especially in winter and in food-deficient areas, food supplementation can enhance the gut microbial homeostasis of wild animals and reduce the impact of crises. In depth studies of the gut microbial function of wildlife have important implications for improving ex situ conservation.

Identifying priority reserves favors the sustainable development of wild ungulates and the construction of Sanjiangyuan National Park.

Sanjiangyuan National Park (SNP), the first national park in China, is one of the most important biodiversity conservation areas in the Sanjiangyuan National Nature Reserve (SNNR) and even the world. The threatened ungulates play an irreplaceable role in maintaining the ecosystem diversity and stability in SNNR. Here, based on 1434 occurrence records of six ungulates, the maximum entropy model, with two different strategies, was utilized to determine the priority reserves. The results indicated that the priority reserves in SNNR was mainly located in and around SNP, which were mainly distributed in the middle east, middle west, and southwest of SNNR. Six ungulates shared preference for altitude ranging 4000-5000 m, the average annual temperature below -3.0°C, and average annual precipitation ranging 200-400 mm on meadow, steppe, and unused land. The proportion of high and medium suitable areas for ungulates in SNP was higher than that in SNNR. As the SNP is not contiguously spaced in space, and some core wildlife habitats are not included, it is suggested to optimize the functional areas and adjust the boundary range on the basis of the pilot scope of SNP, so as to enhance the integrity and connectivity of each functional area.

Comparative Analysis of Gut Microbial Composition and Functions in Przewalski's Gazelle (Procapra przewalskii) From Various Habitats.

Gut microbiota of mammals participates in host nutrient metabolism and plays an important role in host adaptation to the environment. Herein, to understand the relationship between environment differences and the composition and abundance of the gut microbiota of Przewalski's gazelle (Procapra przewalskii) in almost all its habitats, high throughput sequencing of the 16S rRNA gene was used to compared the characteristics of the gut microbiota based on total 120 fecal samples. The results showed that Przewalski's gazelle exhibited different characteristics of microbiota diversity in different habitats. The Jiangxigou Rescue Station (JX), Nongchang (NC), and Ganzihe and Haergai townships (GH) groups had a relatively high microbiota diversity, while the Niaodao scenic area (ND) group had the lowest diversity. This finding seemed to follow a similar pattern of change in the population of Przewalski's gazelle. Bacteroidetes and Actinobacteria were the phyla with significant differences, especially between the Wayu township (WY) and the other groups. The difference in the microbiota mainly included the Ruminococcaceae UCG-005, Christensenellaceae R-7 group, and Bacteroidaceae and was enriched in the ND, WY, and other regions. We speculated that the difference in the gut microbiota was due to a difference in environmental characteristics, particularly the food resources that the host can obtain. We speculated that a similar microbiome has important functions for species survival and represents the evolutionary commonality of Przewalski's gazelle, while a different microbiome plays an important role in the adaptation of Przewalski's gazelle to a different environment. The results of this study illustrate how the same species adapts to different environments from the perspective of gut microbiota plasticity and therefore are of great significance for the protection and restoration of the population of this species.

Genomic insights into zokors' phylogeny and speciation in China.

The phylogeny and speciation of subterranean zokors in China are unclear, as previous studies on morphology and limited molecular markers have generated conflicting results. This study unraveled the complex evolutionary history of eight zokor species in China based on de novo assembly at chromosome level and whole-genome sequencing of 23 populations. We found extensive phylogenetic discordances between nuclear and mitochondrial phylogenies, and different coalescent phylogenies, which could be explained by introgression and incomplete lineage sorting (ILS). The recent Qinghai-Tibet Plateau uplift (∼3.60 million y ago; Mya) drove Eospalax to speciate into clade A and clade B (∼3.22 Mya), and discordant phylogenies in this node were mainly attributed to introgression rather than ILS. Clade A rapidly diverged into three lineages due to geographical isolation and glaciation, while glaciation and C4 plant expansion contributed to the speciation of clade B. ILS contributed to the discordances of two rapidly radiated nodes rather than introgression. The effective population sizes (Ne's) of all the species of Eospalax were affected by three glaciations. Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of all the species pairs. Positively selected genes putatively related to specific inhabitation adaptations were identified, such as heart development, neurogenesis, DNA repair, and immune response. Climate, geological tectonism, and C4 vegetation shaped the adaptation and speciation of zokors in China.

Comparative analysis of gut microbial composition and potential functions in captive forest and alpine musk deer.

Gut microbiota forms a unique microecosystem and performs various irreplaceable metabolic functions for ruminants. The gut microbiota is important for host health and provides new insight into endangered species conservation. Forest musk deer (FMD) and alpine musk deer (AMD) are typical small ruminants, globally endangered due to excessive hunting and habitat loss. Although nearly 60 years of captive musk deer breeding has reduced the hunting pressure in the wild, fatal gastrointestinal diseases restrict the growth of captive populations. In this study, 16S rRNA high-throughput sequencing revealed the differences in gut microbiota between FMD and AMD based on 166 fecal samples. The alpha diversity was higher in FMD than in AMD, probably helping FMD adapt to different and wider habitats. The ß-diversity was higher between adult FMD and AMD than juveniles and in winter than late spring. The phylum Firmicutes and the genera Christensenellaceae R7 group, Ruminococcus, Prevotellaceae UCG-004, and Monoglobus were significantly higher in abundance in FMD than in AMD. However, the phylum Bacteroidetes and genera Bacteroides, UCG-005, Rikenellaceae RC9 gut group, and Alistipes were significantly higher in AMD than FMD. The expression of metabolic functions was higher in AMD than in FMD, a beneficial pattern for AMD to maintain higher energy and substance metabolism. Captive AMD may be at higher risk of intestinal diseases than FMD, with higher relative abundances of most opportunistic pathogens and the expression of disease-related functions. These results provide valuable data for breeding healthy captive musk deer and assessing their adaptability in the wild. KEY POINTS: • Alpha diversity of gut microbiota was higher in FMD than that in AMD • Expression of metabolic and disease-related functions was higher in AMD than in FMD.

Microbial Biogeography along the Gastrointestinal Tract Segments of Sympatric Subterranean Rodents (Eospalax baileyi and Eospalax cansus).

In this study, based on high-throughput sequencing technology, the biodiversity and the community structure of microbiota in different GIT segments (the stomach, small intestine, cecum and rectum) of plateau zokors and Gansu zokors were studied and compared. A source tracking analysis for the microbial communities of different GIT segments was carried out using the fast expectation-maximization microbial source tracking (FEAST) method. We found that, for both species, the microbial community richness and diversity of the small intestine were almost the lowest while those of the cecum were the highest among the four segments of the GIT. Beta diversity analyses revealed that the bacterial community structures of different GIT segments were significantly different. As for the comparison between species, the bacterial community compositions of the whole GIT, as well as for each segment, were all significantly different. Source tracking conducted on both zokors indicated that the soil has little effect on the bacterial community of the GIT. A fairly high percentage of rectum source for the bacterial community of the stomach indicated that both zokors may engage in coprophagy.

Gut Microbiome Changes in Captive Plateau Zokors (Eospalax baileyi).

Wild-caught animals must cope with drastic lifestyle and dietary changes after being induced to captivity. How the gut microbiome structure of these animals will change in response receives increasing attention. The plateau zokor (Eospalax baileyi), a typic subterranean rodent endemic to the Qinghai-Tibet plateau, spends almost the whole life underground and is well adapted to the environmental pressures of both plateau and underground. However, how the gut microbiome of the plateau zokor will change in response to captivity has not been reported to date. This study compared the microbial community structure and functions of 22 plateau zokors before (the WS group) and after being kept in captivity for 15 days (the LS group, fed on carrots) using the 16S rRNA gene via high-throughput sequencing technology. The results showed that the LS group retained 973 of the 977 operational taxonomic units (OTUs) in the WS group, and no new OTUs were found in the LS group. The dominant bacterial phyla were Bacteroides and Firmicutes in both groups. In alpha diversity analysis, the Shannon, Sobs, and ACE indexes of the LS group were significantly lower than those of the WS group. A remarkable difference (P < 0.01) between groups was also detected in beta diversity analysis. The UPGMA clustering, NMDS, PCoA, and Anosim results all showed that the intergroup difference was significantly greater than the intragroup difference. And compared with the WS group, the intragroup difference of the gut microbiota in the LS group was much larger, which failed to support the assumption that similar diets should drive convergence of gut microbial communities. PICRUSt revealed that although some functional categories displayed significant differences between groups, the relative abundances of these categories were very close in both groups. Based on all the results, we conclude that as plateau zokors enter captivity for a short time, although the relative abundances of different gut microbiota categories shifted significantly, they can maintain almost all the OTUs and the functions of the gut microbiota in the wild. So, the use of wild-caught plateau zokors in gut microbial studies is acceptable if the time in captivity is short.

Characterization of a complete mitochondrial genome of Turdus ruficollis (Passeriformes: Turdidae).

Turdus ruficollis is mainly found in China and Northeast Asia. So far, the mitochondrial genome of more than 20 species from the genus Turdus has been studied. However, the relevant information of T. ruficollis has not been reported. To grasp a better comprehension on the molecular basis of T. ruficollis, we obtained the complete mitochondrial DNA genome sequences of this species. The mitogenome was 16,737 bp in length, which consists of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region. A phylogenetic tree based on complete mitogenome sequences revealed that, within the genus Turdidae, T. ruficollis is closely related to T. naumanni and T. eumomus. The complete mitochondrial genome of T. ruficollis would be of great utility for population genetics and phylogeography of the Turdidae family and would also provide meritorious insights for future studies on conservation, genetics, and phylogeny of the Passeriformes family.

Next generation sequencing yields the complete mitogenome of Smith's zokor (Eospalax smithii).

Sequencing analysis of mitochondrial genomes is useful for understanding the genome structures. In this study, complete mitochondrial genomes of the Eospalax smithii was obtained by using next generation sequencing method. The complete mitogenome of E. smithii was 16,350 bp long, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and 1 non-coding control region (D-loop). The overall base composition of the heavy strand is A (33.65%), C (23.80%), T (30.31%), and G (12.24%). The base compositions present highly biased toward A + T nucleotides. The result of phylogenetic analysis showed the five Eospalax species formed a monophyly with the high bootstrap value and as a sister group of the genus Myospalax. This is the first report of the complete mitochondrial genomes of E. smithii and the mitogenome is potentially important for evolutionary biology, population genetics, and species diagnosis studies of the Mysopalacinae.

Characteristics of the mitochondrial genome of Qinling zokor (Eospalax rufescens).

We report the complete mitochondrial genome for the Eospalax rufescens, a typical subterranean rodent species and endemic in China. The resulting E. rufescens mitogenome is 16,355 bases in size, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and 1 noncoding control region (D-loop). The base compositions present highly biased toward A + T nucleotides. Moreover, twelve of all 13 PCGs initiated with ATN start codon, whereas ND1 began with GTG start codon. Stop codons in 13 PCGs were all typical types except incomplete stop codon T for ND4. We further provide a Maximum Likelihood phylogenetic tree showing relationships among E. rufescens and other common subterranean rodents in family Spalacidae.

Pharmacological inhibition of CXCR2 alleviates neuropathic pain by inactivating microglia in a rat L5 spinal nerve ligation model.

Peripheral nerve injury (PNI)-induced neuropathic pain is a prevalent and severe clinical problem. It has been shown that microglia-mediated neuroinflammation plays a crucial role in neuropathic pain. The present study investigated the abnormal expression of C-X-C motif chemokine receptor type 2 (CXCR2) in a rat L5 spinal nerve ligation (SNL) model and evaluated the role of SB225002, a specific antagonist of CXCR2, in repressing neuroinflammation and neuropathic pain. It was found that CXCR2 expression was significantly upregulated in the dorsal horn of L5-SNL rats compared with sham control. Moreover, CXCR2 expression was increased in spinal microglia of rats after L5-SNL. Based on these results, the present study further examined whether pharmacological inhibition of CXCR2 suppressed microglial activation and neuropathic pain. It was demonstrated that SB225002 treatment inhibited L5-SNL-induced microglia proliferation and activation. Furthermore, SB225002 also significantly suppressed the L5-SNL-induced pro-inflammatory response, as indicated by decreased production of tumor necrosis factor-α, interleukin (IL)-1ß and IL-6 in spinal cord tissues. The results indicated that SB225002 also significantly inhibited microglial cell viability and lipopolysaccharide-induced production of pro-inflammatory cytokines in cultured microglia. Functionally, SB225002 treatment effectively repressed mechanical and cold hypersensitivity after peripheral nerve injury. Collectively, the present results suggested that pharmacological inhibition of CXCR2 by SB225002 suppressed L5-SNL-induced neuroinflammation and neuropathic pain, thus offering a potential therapeutic strategy for neuropathic pain treatment.

Musk deer (Moschus spp.) face redistribution to higher elevations and latitudes under climate change in China.

The population of wild musk deer (Moschus spp.) has declined in recent decades and reached an endangered status in China. Global climate change may drive the extinction rate of these species. To understand the implications of global warming on the future potential space utilization and migration direction of musk deer, both the maximum entropy model and barycenter migration analysis were utilized. Five global climate models and four representative concentration pathway scenarios were considered to simulate the distribution of six species for the years 2050 and 2070. The results indicated that the suitable habitat area would decrease over the next 30 to 50 years. These decreases of suitable habitat were more significant for the Siberian musk deer (reduced by 4.98% of the land area of China), the forest musk deer (1.04%), the black musk deer (0.86%), and the Himalayan musk deer (1.82%) compared with the other two musk deer species. The area with suitable climate for the Siberian musk deer will migrate to the southwest (to higher elevations) while areas suitable for the Alpine musk deer, the Himalayan musk deer, and the Anhui musk deer would all migrate to the northeast (to higher latitudes). However, the forest musk deer and the black musk deer will not migrate in the same direction, but will mainly migrate to the west and the north, respectively. These results provide data in support for in-situ conservation, ex-situ conservation, natural reserve community, and bio-corridor construction of China's musk deer species in response to global warming.

Comparison of gut microbiota diversity between wild and captive bharals (Pseudois nayaur).

BACKGROUND: Gastrointestinal microbiota play an important role in animal host immunity, nutrient metabolism, and energy acquisition, and have therefore drawn increasing attentions. This study compared the diversity of the gut microbiota of both wild and captive bharals, which is an ungulate herbivore of caprid from the Qinghai-Tibet plateau. RESULTS: The sequencing of the V4-V5 region of the 16S rRNA gene via high-throughput sequencing technology showed that the dominant bacterial phyla are Firmicutes and Bacteroides both in wild and captive bharals. However, their abundance differed significantly between groups. Firmicutes were significantly higher in wild bharals, while Bacteroides were significantly higher in captive bharals. Different diets are likely a key influencing factor in the diversity and abundance of gut microbiota in bharals. CONCLUSIONS: Changes in diets affect the diversity of gut microbiota and the relative abundance of pathogenic bacteria, increasing the risk of diseases outbreak in captive bharals. The results of this study suggest that the structure and function of the gut microbiota should be regulated via dietary intervention, accurate provision of an individualized diet, and optimization of the functional network of gut microbiota and its interaction with the host. This will improve the ex situ protection of wild animals.

Maxent modeling for predicting the spatial distribution of three raptors in the Sanjiangyuan National Park, China.

Upland buzzard (Buteo hemilasius), Saker falcon (Falco cherrug), and Himalayan vulture (Gyps himalayensis) are three common large raptors in the Sanjiangyuan National Park (SNP), China's first national park. Among them, Upland buzzard and Saker falcon play a significant role in controlling plateau rodent populations and reducing the transmission of pathogens carried by rodents. The Himalayan vulture can provide services for the redistribution and recycling of nutrients in the ecosystem, and play an irreplaceable role in the celestial burial culture of Tibetans in China. Exploring their habitat suitability is important for the protection of the three raptors. Our research was based on the current distribution of Upland buzzard, Saker falcon, and Himalayan vulture that we had extensively surveyed in the Sanjiangyuan National Park from 2016 to 2017. Combined with the correlation analysis of environmental variables, we utilized maximum entropy model (MaxEnt) to evaluate and compare the habitat suitability of the three species in the Sanjiangyuan National Park. Elevation, climate, and human disturbance factors, which had direct or indirect effects on species survival and reproduction, were all included in the model. Among them, elevation was the most important environmental variables affecting the suitability of habitats of three species. Temperature-related factor was another important predictor. The high (>60%) suitable habitat areas for Upland buzzard, Saker falcon, and Himalayan vulture were 73,017.63, 40,732.78, and 61,654.33 km2, respectively, accounted for 59.32%, 33.09%, and 50.08% of the Sanjiangyuan National Park and their total suitable area (i.e., the sum area of high and moderate habitats) reached 96.07%, 60.59%, and 93.70%, respectively. Besides, the three species have overlapping areas for the suitable habitats, which means that overlapping areas should be highly valued and protected. Therefore, understanding the distribution of suitable habitats of the three raptors can provide useful information and reasonable reference for us to put forward suggestions for their protection and regional management.

Adaptive Transcriptome Profiling of Subterranean Zokor, Myospalax baileyi, to High- Altitude Stresses in Tibet.

Animals living at high altitudes have evolved distinct phenotypic and genotypic adaptations against stressful environments. We studied the adaptive patterns of altitudinal stresses on transcriptome turnover in subterranean plateau zokors (Myospalax baileyi) in the high-altitude Qinghai-Tibetan Plateau. Transcriptomes of zokors from three populations with distinct altitudes and ecologies (Low: 2846 m, Middle: 3282 m, High: 3,714 m) were sequenced and compared. Phylogenetic and principal component analyses classified them into three divergent altitudinal population clusters. Genetic polymorphisms showed that the population at H, approaching the uppermost species boundary, harbors the highest genetic polymorphism. Moreover, 1056 highly up-regulated UniGenes were identified from M to H. Gene ontologies reveal genes like EPAS1 and COX1 were overexpressed under hypoxia conditions. EPAS1, EGLN1, and COX1 were convergent in high-altitude adaptation against stresses in other species. The fixation indices (F ST and G ST )-based outlier analysis identified 191 and 211 genes, highly differentiated among L, M, and H. We observed adaptive transcriptome changes in Myospalax baileyi, across a few hundred meters, near the uppermost species boundary, regardless of their relatively stable underground burrows' microclimate. The highly variant genes identified in Myospalax were involved in hypoxia tolerance, hypercapnia tolerance, ATP-pathway energetics, and temperature changes.

Long-Term Follow-Up of the Fellow Eye in Patients Undergoing Surgery on One Eye for Treating Myopic Traction Maculopathy.

Objective. To observe the fellow eye in patients undergoing surgery on one eye for treating myopic traction maculopathy. Methods. 99 fellow eyes of consecutive patients who underwent unilateral surgery to treat MTM were retrospectively evaluated. All patients underwent thorough ophthalmologic examinations, including age, gender, duration of follow-up, refraction, axial length, intraocular pressure, lens status, presence/absence of a staphyloma, and best-corrected visual acuity (BCVA). Fundus photographs and SD-OCT images were obtained. When feasible, MP-1 microperimetry was performed to evaluate macular sensitivity and fixation stability. Results. At an average follow-up time of 24.7 months, 7% fellow eyes exhibited partial or complete MTM resolution, 68% stabilized, and 25% exhibited progression of MTM. Of the 38 eyes with "normal" macular structure on initial examination, 11% exhibited disease progression. The difference in progression rates in Groups 2, 3, and 4 was statistically significant. Refraction, axial length, the frequency of a posterior staphyloma, chorioretinal atrophy, initial BCVA, final BCVA, and retinal sensitivity all differed significantly among Groups 1-4. Conclusions. Long axial length, chorioretinal atrophy, a posterior staphyloma, and anterior traction contribute to MTM development. Patients with high myopia and unilateral MTM require regular OCT monitoring of the fellow eye to assess progression to myopic pre-MTM. For cases exhibiting one or more potential risk factors, early surgical intervention may maximize the visual outcomes.

Sympatric speciation of spiny mice, Acomys, unfolded transcriptomically at Evolution Canyon, Israel.

Spiny mice, Acomys cahirinus, colonized Israel 30,000 y ago from dry tropical Africa and inhabited rocky habitats across Israel. Earlier, we had shown by mtDNA that A. cahirinus incipiently sympatrically speciates at Evolution Canyon I (EC I) in Mount Carmel, Israel because of microclimatic interslope divergence. The EC I microsite consists of a dry and hot savannoid "African" slope (AS) and an abutting humid and cool-forested "European" slope (ES). Here, we substantiate incipient SS in A. cahirinus at EC I based on the entire transcriptome, showing that multiple slope-specific adaptive complexes across the transcriptome result in two divergent clusters. Tajima's D distribution of the abutting Acomys interslope populations shows that the ES population is under stronger positive selection, whereas the AS population is under balancing selection, harboring higher genetic polymorphisms. Considerable sites of the two populations were differentiated with a coefficient of FST = 0.25-0.75. Remarkably, 24 and 37 putatively adaptively selected genes were detected in the AS and ES populations, respectively. The AS genes involved DNA repair, growth arrest, neural cell differentiation, and heat-shock proteins adapting to the local AS stresses of high solar radiation, drought, and high temperature. In contrast, the ES genes involved high ATP associated with energetics stress. The sharp ecological interslope divergence led to strong slope-specific selection overruling the interslope gene flow. Earlier tests suggested slope-specific mate choice. Habitat interslope-adaptive selection across the transcriptome and mate choice substantiate sympatric speciation (SS), suggesting its prevalence at EC I and commonality in nature.

Molecular identification of Echinococcus species from eastern and southern Qinghai, China, based on the mitochondrial cox1 gene.

The Qinghai-Tibetan Plateau (QTP, in western China), which is the largest and highest plateau on Earth, is a highly epidemic region for Echinococcus spp. We collected 70 Echinococcus samples from humans, dogs, sheep, yaks, plateau pikas, and voles in eastern and southern Qinghai and genotyped them using the mitochondrial DNA marker cytochrome oxidase subunit I gene and maximum parsimony and Bayesian reconstruction methods. Based on the 792-bp sequence matrix, we recorded 124 variable sites, of which, 115 were parsimony-informative. Thirty-four haplotypes (H1-H34) were detected, of which H1-H15, H16-H17, and H18-H34 belonged to Echinococcus shiquicus, Echinococcus multilocularis, and Echinococcus granulosus, respectively. Within 26 human isolates, three were identified as E. multilocularis and 23 were E. granulosus. We also detected a dual infection case in a dog with E. multilocularis and E. granulosus. The intraspecific haplotype (Hd ± SD) and nucleotide (Nd ± SD) diversity of E. shiquicus (0.947 ± 0.021; 0.00441 ± 0.00062) was higher than that for E. granulosus (0.896 ± 0.038; 0.00221 ± 0.00031) and E. multilocularis (0.286 ± 0.196; 0.00036 ± 0.00025). Moreover, the haplotype network of E. shiquicus showed a radial feature rather than a divergent feature in a previous study, indicating this species in the QTP has also evolved with bottleneck effects.