Morphology of the nucho‐dorsal glands and related defensive displays in three species of Asian natricine snakes

Chromosome-level genome assembly and annotation of the Rhabdophis nuchalis (Hubei keelback)

Rhabdophis nuchalis, a snake widely distributed in China, possesses a unique trait: glands beneath the skin on its neck and back, known as nucho-dorsal glands. These features make it a valuable subject for studying genetic diversity and the evolution of complex traits. In this study, we obtained a high-quality chromosome-level reference genome of R. nuchalis using MGI short-read sequencing, PacBio Revio long-read sequencing, and Hi-C sequencing techniques. The final assembly comprised 1.92 Gb of the R. nuchalis genome, anchored to 20 chromosomes (including 9 macrochromosomes and 11 microchromosomes), with a contig N50 of 104.79 Mb, a scaffold N50 of 204.96 Mb, and a BUSCO completeness of 97.50%. Additionally, we annotated a total of 1.09 Gb of repetitive sequences (which constitute 56.51% of the entire genome) and identified 22,057 protein-coding genes. This high-quality reference genome of R. nuchalis furnishes essential genomic data for comprehending the genetic diversity and evolutionary history of the species, as well as for facilitating species conservation efforts and comparative genomics studies.

Dramatic dietary shift maintains sequestered toxins in chemically defended snakes

Unlike other snakes, most species of Rhabdophis possess glands in their dorsal skin, sometimes limited to the neck, known as nucho-dorsal and nuchal glands, respectively. Those glands contain powerful cardiotonic steroids known as bufadienolides, which can be deployed as a defense against predators. Bufadienolides otherwise occur only in toads (Bufonidae) and some fireflies (Lampyrinae), which are known or believed to synthesize the toxins. The ancestral diet of Rhabdophis consists of anuran amphibians, and we have shown previously that the bufadienolide toxins of frog-eating species are sequestered from toads consumed as prey. However, one derived clade, the Rhabdophis nuchalis Group, has shifted its primary diet from frogs to earthworms. Here we confirm that the worm-eating snakes possess bufadienolides in their nucho-dorsal glands, although the worms themselves lack such toxins. In addition, we show that the bufadienolides of R. nuchalis Group species are obtained primarily from fireflies. Although few snakes feed on insects, we document through feeding experiments, chemosensory preference tests, and gut contents that lampyrine firefly larvae are regularly consumed by these snakes. Furthermore, members of the R. nuchalis Group contain compounds that resemble the distinctive bufadienolides of fireflies, but not those of toads, in stereochemistry, glycosylation, acetylation, and molecular weight. Thus, the evolutionary shift in primary prey among members of the R. nuchalis Group has been accompanied by a dramatic shift in the source of the species' sequestered defensive toxins.

Evolution of nuchal glands, unusual defensive organs of Asian natricine snakes (Serpentes: Colubridae), inferred from a molecular phylogeny

A large body of evidence indicates that evolutionary innovations of novel organs have facilitated the subsequent diversification of species. Investigation of the evolutionary history of such organs should provide important clues for understanding the basis for species diversification. An Asian natricine snake, Rhabdophis tigrinus, possesses a series of unusual organs, called nuchal glands, which contain cardiotonic steroid toxins known as bufadienolides. Rhabdophis tigrinus sequesters bufadienolides from its toad prey and stores them in the nuchal glands as a defensive mechanism. Among more than 3,500 species of snakes, only 17 Asian natricine species are known to possess nuchal glands or their homologues. These 17 species belong to three nominal genera, Balanophis, Macropisthodon, and Rhabdophis. In Macropisthodon and Rhabdophis, however, species without nuchal glands also exist. To infer the evolutionary history of the nuchal glands, we investigated the molecular phylogenetic relationships among Asian natricine species with and without nuchal glands, based on variations in partial sequences of Mt-CYB, Cmos, and RAG1 (total 2,767 bp). Results show that all species with nuchal glands belong to a single clade (NGC). Therefore, we infer that the common ancestor of this clade possessed nuchal glands with no independent origins of the glands within the members. Our results also imply that some species have secondarily lost the glands. Given the estimated divergence time of related species, the ancestor of the nuchal gland clade emerged 19.18 mya. Our study shows that nuchal glands are fruitful subjects for exploring the evolution of novel organs. In addition, our analysis indicates that reevaluation of the taxonomic status of the genera Balanophis and Macropisthodon is required. We propose to assign all species belonging to the NGC to the genus Rhabdophis, pending further study.