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.

Microbial symbiotic associations may be beneficial, neutral, or harmful to the host. Symbionts exploit the host space and nutrition or use hosts as carriers to spread to other environments. In order to investigate the fecal bacterial communities of wild sika deer (Cervus nippon) and wapiti (Cervus canadensis), this study aimed to sequence and explore the composition of, and similarity between, the fecal microbiota of sika deer and wapiti using high-throughput sequencing. The composition and relative abundance of fecal microbiota, alpha diversity, and differences in beta diversity between the two species were analyzed. We found that no pathogenic bacteria were present in large quantities in the hosts. The dominant bacterial phyla found in the two deer species were similar and included Firmicutes, Bacteroidetes, Proteobacteria, and Spirochaetes. Moreover, the deer also shared similar dominant genera, including the Rikenellaceae RC9 gut group, Ruminococcaceae_UCG-010, Ruminococcaceae_UCG-005, and Bacteroides. These results demonstrate that the sika deer and wapiti share a similar fecal microbiotal structure, probably due to their common diet and living environment, but there was some evidence of a difference at the species level. These analyses provide new insights into the health status of deer populations outside protected environments and offer a scientific framework for monitoring the health conditions of sika deer and wapiti.

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 risk damage from virus invasion due to frequent environmental changes. The oligoadenylate synthesis (OAS) gene family, which is an important component of the immune system, provides an essential response to the antiviral activities of interferons by regulating immune signal pathways. However, little is known about the evolutionary characteristics of OASs in Laurasiatherian mammals. Here, we examined the evolution of the OAS genes in 64 mammals to explore the accompanying molecular mechanisms of the antiviral ability of Laurasiatherian mammals living in different environments. We found that OAS2 and OAS3 were found to be pseudogenes in Odontoceti species. This may be related to the fact that they live in water. Some Antilopinae, Caprinae, and Cervidae species lacked the OASL gene, which may be related to their habitats being at higher altitudes. The OASs had a high number of positive selection sites in Cetartiodactyla, which drove the expression of strong antiviral ability. The OAS gene family evolved in Laurasiatherian mammals at different rates and was highly correlated with the species’ antiviral ability. The gene evolution rate in Cetartiodactyla was significantly higher than that in the other orders. Compared to other species of the Carnivora family, the higher selection pressure on the OAS gene and the absence of positive selection sites in Canidae may be responsible for its weak resistance to rabies virus. The OAS gene family was relatively conserved during evolution. Conserved genes are able to provide better maintenance of gene function. The rate of gene evolution and the number of positively selected sites combine to influence the resistance of a species to viruses. The positive selection sites demonstrate the adaptive evolution of the OAS gene family to the environment. Adaptive evolution combined with conserved gene function improves resistance to viruses. Our findings offer insights into the molecular and functional evolution of the antiviral ability of Laurasian mammals.

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.