Background: A comprehensive evaluation of the -omic profiles of venom is important for understanding the potential function and evolution of snake venom. Here, we conducted an integrated multi-omics-analysis to unveil the venom-transcriptomic and venomic profiles in a same group of spine-bellied sea snakes (Hydrophis curtus) from the South China Sea, where the snake is a widespread species and might generate regionally-specific venom potentially harmful to human activities. The capacity of two heterologous antivenoms to immunocapture the H. curtus venom was determined for an in-depth evaluation of their rationality in treatment of H. curtus envenomation. In addition, a phylogenetic analysis by maximum likelihood was used to detect the adaptive molecular evolution of full-length toxin-coding unigenes. Results: A total of 90,909,384 pairs of clean reads were generated via Illumina sequencing from a pooled cDNA library of six specimens, and yielding 148,121 unigenes through de novo assembly. Sequence similarity searching harvested 63,845 valid annotations, including 63,789 non-toxin-coding and 56 toxin-coding unigenes belonging to 22 protein families. Three protein families, three-finger toxins (3-FTx), phospholipase A₂ (PLA₂), and cysteine-rich secretory protein, were detected in the venom proteome. 3-FTx (27.15% in the transcriptome/41.94% in the proteome) and PLA₂ (59.71%/49.36%) were identified as the most abundant families in the venom-gland transcriptome and venom proteome. In addition, 24 unigenes from 11 protein families were shown to have experienced positive selection in their evolutionary history, whereas four were relatively conserved throughout evolution. Commercial Naja atra antivenom exhibited a stronger capacity than Bungarus multicinctus antivenom to immunocapture H. curtus venom components, especially short neurotoxins, with the capacity of both antivenoms to immunocapture short neurotoxins being weaker than that for PLA₂s. Conclusions: Our study clarified the venom-gland transcriptomic and venomic profiles along with the within-group divergence of a H. curtus population from the South China Sea. Adaptive evolution of most venom components driven by natural selection appeared to occur rapidly during evolutionary history. Notably, the utility of commercial N. atra and B. multicinctus antivenoms against H. curtus toxins was not comprehensive; thus, the development of species-specific antivenom is urgently needed.

Numerous studies have demonstrated that multiple intrinsic and extrinsic factors shape the structure and composition of gut microbiota in a host. The disorder of the gut microbiota may trigger various host diseases. Here, we collected fecal samples from wild-caught Japanese geckos (Gekko japonicus) and captive conspecifics fed with mealworms (mealworm-fed geckos) and fruit flies (fly-fed geckos), aiming to examine the dietary and sexual correlates of the gut microbiota. We used 16S rRNA gene sequencing technology to determine the composition of the gut microbiota. The dominant phyla with a mean relative abundance higher than 10% were Verrucomicrobiota, Bacteroidota, and Firmicutes. Gut microbial community richness and diversity were higher in mealworm-fed geckos than in wild geckos. Neither community evenness nor beta diversity of gut microbiota differed among wild, mealworm-fed, and fly-fed geckos. The beta rather than alpha diversity of gut microbiota was sex dependent. Based on the relative abundance of gut bacteria and their gene functions, we concluded that gut microbiota contributed more significantly to the host’s metabolic and immune functions. A higher diversity of gut microbiota in mealworm-fed geckos could result from higher chitin content in insects of the order Coleoptera. This study not only provides basic information about the gut microbiota of G. japonicus but also shows that gut microbiota correlates with dietary habits and sex in the species.

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A₂ (PLA₂) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA₂ were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 μg/g) and long-chain neurotoxin (i.p. 0.14 μg/g) presented fairly high toxicity, whereas both basic and acidic PLA₂s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

The oviparity-viviparity transition is a major evolutionary event, likely altering the reproductive process of the organisms involved. Residual yolk, a portion of yolk remaining unutilized at hatching or birth as parental investment in care, has been investigated in many oviparous amniotes but remained largely unknown in viviparous species. Here, we used data from 20 (12 oviparous and 8 viviparous) species of snakes to see if the oviparity-viviparity transition alters the partitioning of yolk in embryonic snakes. We used ANCOVA to test whether offspring size, mass and components at hatching or birth differed between the sexes in each species. We used both ordinary least squares and phylogenetic generalized least squares regressions to test whether relationships between selected pairs of offspring components were significant. We used phylogenetic ANOVA to test whether offspring components differed between oviparous and viviparous species and, more specifically, the hypothesis that viviparous snakes invest more in the yolk as parental investment in embryogenesis to produce more well developed offspring that are larger in linear size. In none of the 20 species was sex a significant source of variation in any offspring component examined. Newborn viviparous snakes on average contained proportionally more water and, after accounting for body dry mass, had larger carcasses but smaller residual yolks than did newly hatched oviparous snakes. The rates at which carcass dry mass (CDM) and fat body dry mass (FDM) increased with residual yolk dry mass (YDM) did not differ between newborn oviparous and viviparous snakes. Neither CDM nor FDM differed between newborn oviparous and viviparous snakes after accounting for YDM. Our results are not consistent with the hypothesis that the partitioning of yolk between embryonic and post-embryonic stages differs between snakes that differ in parity mode, but instead show that the partitioning of yolk in embryonic snakes is species-specific or phylogenetically related. We conclude that the oviparity-viviparity transition does not alter yolk partitioning in embryonic snakes.