The Gobiidae family occupy one of the most diverse habitat ranges of all fishes. One key reason for their successful colonization of different habitats is their ability to adapt to different energy demands. This energy requirement is related to the ability of mitochondria in cells to generate energy via oxidative phosphorylation (OXPHOS). Here, we assembled three complete mitochondrial genomes of Rhinogobius shennongensis, Rhinogobius wuyanlingensis, and Chaenogobius annularis. These mitogenomes are circular and include 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs, and one non-coding control region (CR). We used comparative mitochondrial DNA (mtDNA) genome and selection pressure analyses to explore the structure and evolutionary rates of Gobiidae mitogenomics in different environments. The CmC model showed that the ω ratios of all mtDNA PCGs were <1, and that the evolutionary rate of adenosine triphosphate 8 (atp8) was faster in Gobiidae than in other mitochondrial DNA PCGs. We also found evidence of positive selection for several sites of NADH dehydrogenase (nd) 6 and atp8 genes. Thus, divergent mechanisms appear to underlie the evolution of mtDNA PCGs, which might explain the ability of Gobiidae to adapt to diverse environments. Our study provides new insights on the adaptive evolution of Gobiidae mtDNA genome and molecular mechanisms of OXPHOS.

The gastrointestinal tract of animals contains microbiota, forming a complex microecosystem. Gut microbes and their metabolites can regulate the development of host innate and adaptive immune systems. Animal immune systems maintain intestinal symbiotic microbiota homeostasis. However, relatively few studies have been published on reptiles, particularly snakes, and even fewer studies on different parts of the digestive tracts of these animals. Herein, we used 16S rRNA gene sequencing to investigate the microbial community composition and adaptability in the stomach and small and large intestines of Lycodon rufozonatus. Proteobacteria, Bacteroidetes, and Firmicutes were most abundant in the stomach; Fusobacteria in the small intestine; and Proteobacteria, Bacteroidetes, Fusobacteria, and Firmicutes in the large intestine. No dominant genus could be identified in the stomach; however, dominant genera were evident in the small and large intestines. The microbial diversity index was significantly higher in the stomach than in the small and large intestines. Moreover, the influence of the microbial community structure on function was clarified through function prediction. Collectively, the gut microbes in the different segments of the digestive tract revealed the unique features of the L. rufozonatus gut microbiome. Our results provide insights into the co-evolutionary relationship between reptile gut microbiota and their hosts.

Many mammals develop specific immune responses owing to the changes in their ecological niche and diet that are essential for animal survival. However, pattern recognition receptors (PRRs) serve as the first line of defense in innate immunity and generate immune responses in the host. However, the evolutionary impacts on PRR genes in Carnivora are not well studied. Herein, we explored the evolution of 946 PRR gene sequences in 43 Carnivora species to elucidate the molecular mechanisms of carnivore adaptation to complex habitats. We found that the PRRs were relatively conserved, and different gene families showed different evolutionary patterns. PRRs were highly purified based on their overall roles in Carnivora species but interspersed with positive-selection patterns during evolution. Different niche types may have jointly driven the evolution of PRR genes. In particular, the selection pressure of toll-like receptor (TLR) 10 was relaxed in seven species with pseudogenes, which may have emerged during recent evolutionary events. We speculated that a “functional compensation” mechanism may exist for genes with overlapping functions in the TLR gene family. Additionally, TLR2, TLR4, NLRC5, and DECTIN1 were subject to positive selection in semi-aquatic species, and the adaptive evolution of these genes may have been related to the adaptation to semi-aquatic environments. In summary, our findings offer valuable insights into the molecular and functional evolution of PRR genes, which are important for immune adaptations in Carnivora.

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species’ phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters’ adaptation to semi-aquatic habitats.

Host phylogeny can imprint the gut microbiota, but it is often masked by diet and environment. Leveraging the standardized husbandry of a zoological collection, we profiled fecal microbiota from 55 captive artiodactyls representing 12 species in Bovidae, Cervidae, and Camelidae using 16S rRNA amplicon sequencing targeting the V3–V4 region on the Illumina MiSeq platform. Community composition differed significantly among host families (Bray–Curtis PERMANOVA, R² = 0.1075, p = 0.001). A host–microbiota tanglegram, which juxtaposes the host phylogeny with a dendrogram of microbiota similarity, recovered a topology congruent with the host phylogeny, with camelids forming a distinct branch separate from true ruminants in both trees. The linear discriminant analysis effect size (LEfSe; LDA ≥ 3.5) identified family-specific biomarkers, including enrichment of Acinetobacter/Moraxellaceae in Bovidae, Rikenellaceae (the Rikenellaceae_RC9_gut_group) in Cervidae, and Rummeliibacillus together with the Christensenellaceae_R-7_group in Camelidae. Functional inference with PICRUSt2 revealed significant differences in KEGG level-2 pathways (e.g., carbohydrate metabolism and xenobiotics biodegradation), consistent with taxonomic shifts. Altogether, these findings show that—even under uniform captive conditions—host evolutionary history remains a primary determinant of both the structure and the predicted functions of the artiodactyl gut microbiota, refining the scope of phylosymbiosis and providing actionable baselines for veterinary monitoring and conservation management in zoo settings.

Bacillus paranthracis, a species of the genus Bacillus, is a Gram-positive bacterium classified as an opportunistic pathogen that can cause foodborne diarrhea and other intestinal diseases in humans and various animals. To date, there has been limited research on B. paranthracis, and there are few records of this bacterium being isolated from animal intestines. In this study, a strain named Qf-1 was isolated and purified from faecal samples of mink. Through culturomics, 16S rRNA gene sequencing, whole-genome sequencing, and average nucleotide identity (ANI) analysis, the strain was confirmed to be B. paranthracis. Whole-genome sequencing revealed that the strain has a genome size of 5.27 Mb, comprising one chromosome (5,224,739 bp) and one plasmid (51,964 bp). Functional annotation of its genome identified multiple potential pathogenic factors associated with pneumonia, including the key genes AsbD and AsbF, which facilitate bacterial colonisation of the lungs and trigger inflammatory responses, as well as EsxB and EsxL, which exacerbate lung inflammation and promote infection spread. Comparative genomics analysis revealed that this strain shares a close evolutionary relationship with previously reported B. paranthracis strains. The structure and function of the bacterial genes were analyzed in depth using multi-omics methods. Through mouse pathogenicity experiments, it was found that this bacterium may cause pneumonia and enteritis in mice. We predict that it may also pose a threat to the health of the mink. These research findings contribute to the establishment of a stable experimental model between pathogens and mink hosts, laying the foundation for further elucidating their pathogenicity and pathogenic mechanisms. This is of great significance for the diagnosis and prevention of bacterial diseases in mink in the future.