Back

Comparative analysis of cloacal microbiota in Henophidia (non-venomous) and Caenophidia (venomous) snakes

Ghasemian, E.; Nassirnia, S.; Pillonel, T.; Ruegg, S.; Aeby, S.; Bertelli, C.; Borel, N.; Greub, G.

2026-05-14 ecology
10.64898/2026.05.13.724777 bioRxiv
Show abstract

The evolutionary divergence between Henophidia (non-venomous) and Caenophidia (venomous) snakes has produced distinct cranial morphologies, digestive strategies, and presence of specialised venom systems in Caenophidia, yet the extent to which these long-standing diverging trajectories have shaped cloacal microbiota assembly remains poorly understood. We characterised cloacal microbiota in 70 captive snakes (52 Caenophidia, 18 Henophidia) by 16S rRNA amplicon sequencing. Beta diversity was tested by PERMANOVA, differential abundance by ANCOM-BC2, community types by Dirichlet Multinomial Mixture modelling (DMM), and microbial interactions by SparCC co-occurrence networks. Predicted functional potential (PICRUSt2) was analysed by ALDEx2 differential abundance testing and elastic net feature selection. Henophidia exhibited significantly higher bacterial richness and greater compositional variability than Caenophidia. Community composition showed clade-associated differences (PERMANOVA) and partitioned into two distinct DMM community types. The Henophidia network was 11.9-fold denser and more modular, with Burkholderiaceae as a keystone hub, whereas the Caenophidia network was sparse. Henophidia showed predicted enrichment in C1 metabolic pathways (ethylmalonyl-CoA, formaldehyde assimilation I, glycine betaine degradation I, methylaspartate cycle), aromatic compound catabolism, and nitrogen recycling, whilst Caenophidia showed enrichment in allantoin and glucuronate degradation. This multi-method analysis suggests Burkholderiaceae as a candidate keystone taxon in Henophidia and indicates that phylogenetic clade is a major contributor to cloacal microbiota structure. The lower richness in Caenophidia raises a testable hypothesis that broad-spectrum antimicrobial activity of their venom components may selectively filter susceptible microbial lineages, motivating future shotgun metagenomic studies in wild populations of snakes.

Matching journals

The top 6 journals account for 50% of the predicted probability mass.

1
Molecular Ecology
304 papers in training set
Top 0.2%
19.6%
2
Scientific Reports
3102 papers in training set
Top 5%
10.5%
3
Animal Microbiome
26 papers in training set
Top 0.1%
8.5%
4
Microbiome
139 papers in training set
Top 0.5%
6.4%
5
Nature Communications
4913 papers in training set
Top 35%
4.3%
6
Communications Biology
886 papers in training set
Top 1%
4.0%
50% of probability mass above
7
Proceedings of the Royal Society B: Biological Sciences
341 papers in training set
Top 2%
3.6%
8
PLOS ONE
4510 papers in training set
Top 38%
3.6%
9
Frontiers in Ecology and Evolution
60 papers in training set
Top 0.8%
3.6%
10
mSystems
361 papers in training set
Top 4%
2.5%
11
eLife
5422 papers in training set
Top 34%
2.4%
12
Environmental Microbiology
119 papers in training set
Top 1%
1.9%
13
Frontiers in Microbiology
375 papers in training set
Top 4%
1.9%
14
PeerJ
261 papers in training set
Top 7%
1.7%
15
Peer Community Journal
254 papers in training set
Top 2%
1.7%
16
Microbial Ecology
28 papers in training set
Top 0.2%
1.5%
17
Science of The Total Environment
179 papers in training set
Top 3%
1.5%
18
Nature Ecology & Evolution
113 papers in training set
Top 3%
1.2%
19
BMC Genomics
328 papers in training set
Top 4%
1.1%
20
Environmental Microbiome
26 papers in training set
Top 0.4%
0.9%
21
iScience
1063 papers in training set
Top 26%
0.9%
22
Microorganisms
101 papers in training set
Top 2%
0.9%
23
The ISME Journal
194 papers in training set
Top 2%
0.8%
24
FEMS Microbiology Ecology
47 papers in training set
Top 0.4%
0.8%
25
Philosophical Transactions of the Royal Society B: Biological Sciences
53 papers in training set
Top 1%
0.7%
26
Ecology and Evolution
232 papers in training set
Top 4%
0.6%
27
Global Ecology and Biogeography
41 papers in training set
Top 0.7%
0.6%
28
BMC Biology
248 papers in training set
Top 6%
0.6%
29
Royal Society Open Science
193 papers in training set
Top 6%
0.5%
30
The Lancet Microbe
43 papers in training set
Top 2%
0.5%