Canidae is a family of carnivores comprises about 36 extant species that have been defined as three distinct monophyletic groups based on multi-gene data sets. The Tibetan fox ( Vulpes ferrilata ) is a member of the family Canidae that is endemic to the Tibetan Plateau and has seldom been in the focus of phylogenetic analyses. To clarify the phylogenic relationship of V . ferrilata between other canids, we sequenced the mitochondrial genome and firstly attempted to clarify the relative phylogenetic position of V . ferrilata in canids using the complete mitochondrial genome data. The mitochondrial genome of the Tibetan fox was 16,667 bp, including 37 genes (13 protein-coding genes, 2 rRNA, and 22 tRNA) and a control region. A comparison analysis among the sequenced data of canids indicated that they shared a similar arrangement, codon usage, and other aspects. A phylogenetic analysis on the basis of the nearly complete mtDNA genomes of canids agreed with three monophyletic clades, and the Tibetan fox was highly supported as a sister group of the corsac fox within Vulpes . The estimation of the divergence time suggested a recent split between the Tibetan fox and the corsac fox and rapid evolution in canids. There was no genetic evidence for positive selection related to high-altitude adaption for the Tibetan fox in mtDNA and following studies should pay more attention to the detection of positive signals in nuclear genes involved in energy and oxygen metabolisms.

The wolf (Canis lupus) is one of the most widely distributed terrestrial mammals, because it is well adapted to various ecological niches and their corresponding pathogen environments. Immunological competence is a crucial factor involved in adapting to a changing environment and fighting pathogen infection in animals. In this study, the peripheral blood transcriptome of wolves was generated via RNA-seq to advance understanding of the wolf immunome, with a special focus on the major histocompatibility complex class I (MHC I) and toll-like receptor (TLR) gene families, which are involved in pathogen recognition and defense. The blood transcriptomic libraries of eight wolves originating from Tibet and Inner Mongolia were sequenced, and approximately 383 million reads were generated. Using a genome-guided assembly strategy, we obtained 123,851 unigenes, with a mean length of 845 bp and an N50 length of 1121 bp. On the basis of BLAST searches against the NCBI non-redundant protein database (Nr), a total of 36,192 (29.22%) unigenes were annotated. For functional classification, 24,663 unigenes were assigned to 13,016 Gene Ontology (GO) terms belonging to 51 sub-categories of the three main GO categories. Additionally, 7682 unigenes were classified into 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, in which the most represented functional sub-categories were signal transduction and the immune system, and 16,238 unigenes were functionally classified into 25 Eukaryotic Orthologous Groups (KOG) categories. We observed an overall higher ω (d N/d S) value at antigen-binding sites (ABSs) than at non-ABS regions as well as clear evidence of intergenic/intragenic recombination events at wolf MHC I loci. Additionally, our analysis revealed that carnivorous TLRs were dominated by purifying selection, with mean ω values at each TLR locus ranging from 0.173 to 0.527. However, we also found significant instances of positive selection that acted on several codons in pathogen recognition domains and were linked to species-specific differences in pathogen recognition. This study represents the first attempt to characterize the blood transcriptome of the wolf and to highlight the value of investigating the immune system. Balancing selection and recombination have contributed to the historical evolution of wolf MHC I genes. Moreover, TLRs in carnivores have undergone adaptive evolution against the background of purifying selection, and a high level of adaptive evolution was detected in the wolf TLR system.

This study aimed to identify the effects of host species on the gut microbial flora in three species (Hemitragus jemlahicus, Pseudois nayaur, and Ovis orientalis) from the subfamily Caprinae, by excluding the impact of environment factors. We investigated the differences in intestinal flora of three species belonging to Caprinae, which were raised in identical conditions. Fecal samples were collected from tahr, mouflon, and bharal, and the V3–V4 region of the 16S ribosomal RNA gene was analyzed by high‐throughput sequencing. The analysis of 16S rRNA gene sequences reveals that fecal samples were mainly composed of four phyla: Firmicutes, Bacteroidetes, Spirochaetes, and Proteobacteria. The most abundant phyla included Firmicutes and Bacteroidetes accounting for >90% of the bacteria, and a higher Firmicutes/Bacteroidetes ratio was observed in tahrs. Moreover, significant differences existed at multiple levels of classifications in the relative abundance of intestinal flora, differing greatly between species. Phylogenetic analyses based on 16S rRNA gene indicated that mouflon is closely related to bharal, and it is inconsistent with previous reports in the species evolutionary relationships. In this study, we demonstrated that the gut microbiota in tahr had a stronger ability to absorb and store energy from the diet compared with mouflon and bharal, and the characteristics of host–microbiome interactions were not significant.

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

The silver fox and blue fox are economically important fur species and were domesticated by humans from their wild counterparts, the arctic fox and red fox, respectively. Farmed foxes show obvious differences from their wild counterparts, including differences in physiology, body size, energy metabolism, and immunity. However, the molecular mechanisms underlying these differences are presently unclear. In this study, we built transcriptome libraries from multiple pooled tissues for each species of farmed fox, used RNA‐seq to obtain a comprehensive dataset, and performed selection analysis and sequence‐level analyses of orthologous genes to identify the genes that may be influenced by human domestication. More than 153.3, 248.0, 81.6, and 65.8 million clean reads were obtained and assembled into a total of 118,577, 401,520, 79,900, and 186,988 unigenes with an average length range from 521 to 667 bp for AF, BF, RF, and SF, respectively. Selective pressure analysis showed that 11 and 14 positively selected genes were identified, respectively, in the two groups (AF vs. BF and RF vs. SF). Several of these genes were associated with natural immunity (CFI and LRRFIP1), protein synthesis (GOLGA4, CEP19 and SLC35A2), and DNA damage repair (MDC1). Further functional enrichment analyses demonstrated that two positively selected genes (ACO1 and ACAD10) were involved in metabolic process (GO:0008152, p‐value = .032), representing a significant enrichment. Sequence analysis of 117 orthologous genes shared by the two groups showed that the LEMD2, RRBP1, and IGBP1 genes might be affected by artificial selection in farmed foxes, with mutation sites located within sequences that are otherwise highly conserved across most mammals. Our results provide a valuable transcriptomic resource for future genetic studies and improvement in the assisted breeding of foxes and other farmed animals.

The evolutionary patterns of the mitochondrial genome are influenced by both adaptive and nonadaptive forces, with their contributions varying among taxa. There appears to be a correlation linking mutagenesis and latitude, which could be due to differences in metabolic rates. These discrepancies in metabolic rates exhibit a positive connection with mutation pressure. On this basis, we hypothesise that nonadaptive forces play a role in the differences in mutation rates observed along latitudinal gradients. In this study, we selected widely distributed carnivores as representatives of mammals to test our hypothesis. We examined the correlations between the dN/dS ratio (ω), as well as the substitution rates (dS and dN), of 13 PCGs in the mtDNA of 122 carnivores, and the latitude and climatic factors. We found that taxa distributed in higher latitudes tend to have higher substitution rates, but not ω values indicating selective pressure. Notably, dN shows a strong positive correlation with dS, although dS is primarily influenced by mutation pressure, while dN is also influenced by effective population size (N_e). Phylogenetic generalised least squares (PGLS) regression analyses showed that both substitution rates were correlated with climatic factors representing the temperature, precipitation and variability of climate. Based on our findings, we propose that the mutations are primarily influenced by nonadaptive forces (mutation pressure). This forms the fundamental premise for natural selection and speciation. Moreover, the correlation between substitution rates and latitudinal distribution and climate, which are outcomes of nonadaptive factors, can aid in comprehending the global distribution of species diversity.