Proliferative strongyloidiasis in a colony of colubrid snakes.

Strongyloides are small rhabditid nematodes primarily associated with enteric disease in a variety of animal species, including reptiles. Strongyloides spp life stages were associated with a disease outbreak in a large breeding colony of snakes. Multiple Pantherophis and Lampropeltis colubrids exhibited respiratory distress, anorexia, stomatitis, facial deformation, and waning body condition that resulted in death or necessitated euthanasia. Postmortem examinations of 13 snakes revealed epithelial hyperplasia and inflammation of the alimentary and respiratory tracts associated with varying numbers of adult and larval nematodes and embryonated or larvated ova. In a subset of snakes, aberrant nematode migration was also observed in the eye, genitourinary system, coelom, and vasculature. Histomorphology and gross examination of parasitic adult female nematodes from host tissues were consistent with a Strongyloides spp. Sedimented fecal material from 101/160 (63%) snakes housed in the affected facility was positive for nematodes and/or larvated ova. Polymerase chain reaction amplification and sequencing of portions of the 18S and 28S ribosomal ribonucleic acid (RNA) genes and the internal transcribed spacer region of adult female parasites and positive fecal samples supported the diagnosis of strongyloidiasis. Strongyloides spp possess a unique life cycle capable of alternating between parasitic (homogonic) and free-living (heterogonic) stages, resulting in the production of directly infective larvae. Commonly utilized husbandry practices in reptile collections can amplify the numbers of infective larvae generated in the captive environment, increasing the risk for rhabditid hyperinfections. This report documents morbidity, mortality, and non-enteric disease manifestations due to Strongyloides hyperinfections in a captive colubrid snake colony.

Reference genome of the California glossy snake, Arizona elegans occidentalis: A declining California Species of Special Concern.

The glossy snake (Arizona elegans) is a polytypic species broadly distributed across southwestern North America. The species occupies habitats ranging from California's coastal chaparral to the shortgrass prairies of Texas and southeastern Nebraska, to the extensive arid scrublands of central México. Three subspecies are currently recognized in California, one of which is afforded state-level protection based on the extensive loss and modification of its preferred alluvial coastal scrub and inland desert habitat. We report the first genome assembly of A. elegans occidentalis as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technologies to produce a de novo assembled genome. The assembly comprises a total of 140 scaffolds spanning 1,842,602,218 base pairs, has a contig NG50 of 61 Mb, a scaffold NG50 of 136 Mb, and a BUSCO complete score of 95.9%, and is one of the most complete snake genome assemblies. The A. e. occidentalis genome will be a key tool for understanding the genomic diversity and the basis of adaptations within this species and close relatives within the hyperdiverse snake family Colubridae.

Acute Coronary Syndrome From Green Snake Envenomation.

BACKGROUND: Snake bite is a grossly underreported public health issue in subtropical, tropical suburban, and rural areas of Africa and South Asia. In literature, ophitoxemia (snake bite envenomation) as a cause of acute coronary syndrome (ACS) is limited to very few case reports. Viper envenomation is the most common cause of ACS among snake bites. We report the first case of unstable angina caused by Colubridae snake bite (Ahaetullanasuta, commonly called green snakes) in a young man without comorbidities. CASE REPORT: A young healthy man had a green snake bite that was camouflaged in the green fodder. He was managed elsewhere with anti-snake serum. He developed acute chest pain and breathlessness on day 3 of his treatment. Electrocardiogram (ECG) showed biphasic T wave inversions suggestive of type A Wellens pattern in the anterior chest leads (V1-V4). He was treated for ACS medically outside and was referred to our institute for further management on the following day. ECG and cardiac enzymes were normal. The echocardiogram showed no regional wall motion abnormality. Computed tomography coronary angiography showed normal epicardial coronaries. He was discharged in stable condition and asymptomatic at 2 months follow-up. WHY SHOULD AN EMERGENCY PHYSICIAN BE AWARE OF THIS?: ACS after a snake bite is not limited to venomous snakes. The diagnosis should be considered promptly even with a nonvenomous snake bite, especially in those with typical symptoms and ECG changes. The time interval between snake bite and development of ACS can be long and warrants prolonged medical supervision.

A new species of Serpentirhabdias Tkach, Kuzmin et Snyder, 2014 (Nematoda: Rhabdiasidae) parasitic in the herald snake, Crotaphopeltis hotamboeia (Laurenti) (Reptilia: Serpentes: Colubridae) in South Africa.

Serpentirhabdias mamlambo n. sp. is described from the lung of the herald snake, Crotaphopeltis hotamboeia (Laurenti) in KwaZulu-Natal Province, South Africa. The new species is characterised by the small lips, the absence of the buccal capsule, the large excretory glands (1.5-2.3 times longer than the oesophagus), and the numerous eggs in uteri. Serpentirhabdias mamlambo n. sp. differs from other Serpentirhabdias spp. in the presence of M-shaped cuticular folds on dorsal and ventral sides of apical surface, the eye-shaped internal labial papillae, and the distinct dilatation of the oesophagus at its mid-length. Morphological differentiation of S. mamlambo n. sp. from its congeners is supported by molecular evidence based on cox1 mitochondrial gene, and ITS and partial 28S region of nuclear rDNA. Serpentirhabdias mamlambo n. sp. is the first species of the genus Serpentirhabdias described from the Afrotropical Realm.

Out of the Hengduan Mountains: Molecular phylogeny and historical biogeography of the Asian water snake genus Trimerodytes (Squamata: Colubridae).

The Asian water snake genus Trimerodytes (formerly Sinonatrix) is endemic to East and Southeast Asia. Although several species have been included in various phylogenetic studies previously, the evolution and relationships among members of this genus as a whole remain unexplored. In this study, we report the sequencing two protein-coding mitochondrial gene fragments (MTCYB and ND2) and three nuclear genes (c-mos, NT3, and Rag1), reconstruct interspecific phylogeny, and explore biogeography for the genus Trimerodytes. Both Bayesian inference and maximum likelihood analyses consistently recover the monophyly of Trimerodytes with strong support, with T. yapingi the sister-group to the remaining species. The divergence date and ancestral area estimation suggest that Trimerodytes likely originated in Hengduan Mountains (eastern Tibetan Plateau) in western China at 23.93 Ma (95% HPD: 17.09-31.30), and intraspecific divergence began at about 4.23 Ma (95% HPD: 2.74-6.10). Analyses support the validity of T. yunnanensis.

Snakes as Novel Biomarkers of Mercury Contamination: A Review.

Mercury (Hg) is an environmental contaminant that has been reported in many wildlife species worldwide. The organic form of Hg bioaccumulates in higher trophic levels, and thus, long-lived predators are at risk for higher Hg exposure. Although ecological risk assessments for contaminants such as Hg include pertinent receptor species, snakes are rarely considered, despite their high trophic status and potential to accumulate high levels of Hg. Our current knowledge of these reptiles suggests that snakes may be useful novel biomarkers to monitor contaminated environments. The few available studies show that snakes can bioaccumulate significant amounts of Hg. However, little is known about the role of snakes in Hg transport in the environment or the individual-level effects of Hg exposure in this group of reptiles. This is a major concern, as snakes often serve as important prey for a variety of taxa within ecosystems (including humans). In this review, we compiled and analyzed the results of over 30 studies to discuss the impact of Hg on snakes, specifically sources of exposure, bioaccumulation, health consequences, and specific scientific knowledge gaps regarding these moderate to high trophic predators.

One-Step Affinity Purification of Leucine-Rich α2-Glycoproteins from Snake Sera and Characterization of Their Phospholipase A2-Inhibitory Activities as ß-Type Phospholipase A2 Inhibitors.

Snakes contain three types of phospholipase A2 (PLA2)-inhibitory proteins in their blood, PLIα, ß, and γ, which protect them from their own venom, PLA2. PLIß is the snake ortholog of leucine-rich α2 glycoprotein (LRG). Since autologous cytochrome c (Cyt c) serves as an endogenous ligand for LRG, in this study, we purified snake LRGs from various snake serum samples using Cyt c affinity chromatography. All purified snake LRGs were found to be dimers linked by disulfide bonds. Laticauda semifasciata and Naja kaouthia LRGs showed no inhibitory activity against L. semifasciata PLA2 and weak inhibitory activity against Gloydius brevicauda basic PLA2. Elaphe climacophora PLIß had weaker inhibitory activity against G. brevicauda basic PLA2 than G. brevicauda and Elaphe quadrivirgata PLIs, which are abundant in blood and known to neutralize G. brevicauda basic PLA2. Protobothrops flavoviridis LRG showed no inhibitory activity against basic venom PLA2, PL-X, or G. brevicauda basic PLA2. Binding analysis of P. flavoviridis LRG using surface plasmon resonance showed very strong binding to snake Cyt c, followed by that to horse Cyt c, weak binding to yeast Cyt c, and no binding to P. flavoviridis PL-X or BPI/II. We also deduced the amino acid sequences of L. semifasciata and P. flavoviridis LRG by means of cDNA sequencing and compared them with those of other known sequences of PLIs and LRGs. This study concluded that snake LRG can potentially inhibit basic PLA2, but, whether it actually functions as a PLA2-inhibitory protein, PLIß, depends on the snake.

Conservation of Major Satellite DNAs in Snake Heterochromatin.

Repetitive DNA sequences constitute a sizeable portion of animal genomes, and tandemly organized satellite DNAs are a major part of them. They are usually located in constitutive heterochromatin clusters in or near the centromeres or telomeres, and less frequently in the interstitial parts of chromosome arms. They are also frequently accumulated in sex chromosomes. The function of these clusters is to sustain the architecture of the chromosomes and the nucleus, and to regulate chromosome behavior during mitosis and meiosis. The study of satellite DNA diversity is important for understanding sex chromosome evolution, interspecific hybridization, and speciation. In this work, we identified four satellite DNA families in the genomes of two snakes from different families: Daboia russelii (Viperidae) and Pantherophis guttatus (Colubridae) and determine their chromosomal localization. We found that one family is localized in the centromeres of both species, whereas the others form clusters in certain chromosomes or subsets of chromosomes. BLAST with snake genome assemblies showed the conservation of such clusters, as well as a subtle presence of the satellites in the interspersed manner outside the clusters. Overall, our results show high conservation of satellite DNA in snakes and confirm the "library" model of satellite DNA evolution.

Exotic (non-native) snakebite envenomation in Japan: A review of the literature between 2000 and 2022.

A limited number of studies have investigated the clinical characteristics of snakebite envenomation by exotic (non-native) snakes in Japan. This study reviewed the literature to determine the status and risk of bites by exotic pet snakes in Japan. We reviewed reports of snakebite due to exotic snakes in Japan published between 2000 and 2022, excluding reports of bites by snakes native to Japan, such as Gloydius blomhoffii, Rhabdophis tigrinus, and Protobothrops flavoviridis. During the study period, 11 exotic snakebites were recorded, involving 11 species. The majority of those bitten (10/11 cases) were male, all cases were hand injuries, and there were no fatalities. The snakes responsible belonged to the Colubridae (4/11 cases), Viperidae (4/11 cases), and Elapidae (3/11 cases) families. Cases of envenomation by G. brevicaudus, Bungarus candidus, and Dendroaspis angusticeps were of particular interest. Ten of the eleven patients developed local cytotoxic signs, and three developed "compartment syndrome," in which the surgeons performed decompressive incisions. Two bites from elapid snakes and one from a viperid snake resulted in respiratory failure. Antivenom was given in two cases. Complications were observed, such as acute kidney injury, rhabdomyolysis, coagulopathy, and residual dysfunction of the affected finger. Emergency rooms should be prepared to manage patients who have been bitten by exotic snakes, even though the number of reported cases is not high in Japan. Initial stabilization of patients is crucial, before a definitive diagnosis is made, as with native snakebite envenomation. Finger bites are reported in most cases, which may result in functional impairment of the fingers. In order to collect more comprehensive patient data in Japan, a reporting system for all snakebite envenomations should be considered.

The Endangered Sardinian Grass Snake: Distribution Update, Bioclimatic Niche Modelling, Dorsal Pattern Characterisation, and Literature Review.

The Sardinian grass snake, Natrix helvetica cetti, is an endangered endemic snake subspecies with a restricted and highly fragmented geographic distribution. Information on its ecology and detailed geographic distribution are scarce and may negatively impact on its conservation status. Therefore, a literature review on its taxonomy, morphology, ecology, and conservation is presented here. Moreover, field records from the authors, citizen science and the existing literature provide an updated geographic distribution highlighting its presence within 13 new and 7 historic 10 × 10 km cells. Bioclimatic niche modelling was then applied to explore patterns of habitat suitability and phenotypic variation within N. h. cetti. The geographic distribution of the species was found to be positively correlated with altitude and precipitation values, whereas temperature showed a negative correlation. Taken together, these outcomes may explain the snake's presence, particularly in eastern Sardinia. In addition, analysis of distribution overlap with the competing viperine snake (N. maura) and the urodeles as possible overlooked trophic resources (Speleomantes spp. and Euproctus platycephalus) showed overlaps of 66% and 79%, respectively. Finally, geographical or bioclimatic correlations did not explain phenotypic variation patterns observed in this highly polymorphic taxon. Perspectives on future research to investigate N. h. cetti's decline and support effective conservation measures are discussed.

Cytotoxicity and inflammation induced by Philodryas patagoniensis venom.

The Green racer Philodryas patagoniensis is a snake species from South America and accidents with this genus are often neglected. Therefore, this study aimed to evaluate the toxicological, cytotoxic, and inflammatory potential of P. patagoniensis venom (PpV). The experimental model Artemia salina was used to determine toxicity through the median lethal dose (LD50). Cell viability and genotoxicity were evaluated in human mononuclear cells using the Trypan blue test and the Comet assay, respectively. To assess inflammation, mice had the ventral surface of the right hind paw injected with saline, formalin, and three different concentrations of venom (1, 1.5, and 2 µg. 50 µL-1). LD50 in A. salina was 461 µg. mL-1. PpV caused a significant increase in cell death and genotoxicity in human mononuclear cells at two concentrations (575 and 1150 µg. mL-1). PpV shown also to be a strong agent causing nociception in mice. Paw edema totaled four days at 1.5 µg. 50 µL-1. The hyperalgesia caused by the venom had a long duration in mice, lasting eight days at all concentrations evaluated. Thus, we evaluated for the first time the toxicological potential of PpV in A. salina model and in leukocytes. We concluded that systemic oxidative stress, which we infer to be in the genesis of cytotoxicity and genotoxicity observed in vitro, and the inflammatory process are part of the pathways that trigger the venom damage cascades. Relevant data for both scientific research and clinical medicine. Nonetheless, studies are needed to elucidate these mechanisms.

A Multi-Institutional Collaboration to Understand Neoplasia, Treatment and Survival of Snakes.

This multi-institutional collaborative study of neoplasia in snakes reviewed medical records of snakes at each facility to determine species prevalence, survival, and methods of treatment. Complete species numbers of snakes were also collected at each facility. In total, 65 species, 133 snakes, and 149 unique neoplasias were included in this study. Affected species, age, sex, and their tumor prevalence, tumor type and location, metastasis, treatment, and survival data are reported. The highest species-specific tumor prevalence was in Common or Northern Watersnakes (Nerodia sipedon) (30.8%, n = 4 of 13), Eastern Diamond-Backed Rattlesnakes (Crotalus adamanteus) (26.3%, n = 5 of 19), and Timber rattlesnakes (Crotalus horridus) (22.7%, n = 5 of 22). Malignant tumors predominated (86.6%, n = 129 of 149) with soft tissue sarcomas being the most common (30.2%, n = 45 of 149). Snakes with malignant neoplasia, metastases, or indeterminate presence of metastases were statistically more likely to die from their neoplasms than snakes having either benign neoplasia or no diagnosed metastases (p < 0.05). Gender, taxonomic family, and species of those evaluated did not significantly affect the outcome of snakes with neoplasia. Only 27.1% (n = 36 of 133) of snakes received a reported form of treatment and, for those treated, surgical excision was the most common treatment modality. There was not a significant difference in outcome based on treatment; however, surgery and chemotherapy were associated with death from a cause other than their tumor.

Leaving no stone unturned: three additional new species of Atractus ground snakes (Serpentes, Colubridae) from Ecuador discovered using a biogeographical approach.

The genus Atractus includes 146 species of cryptozoic snakes occurring from Panama to northeastern Argentina. Here, a molecular phylogeny of this genus is presented, which encompasses 29% (= 42; six are included here for the first time) of the species currently recognized. Morphological and phylogenetic support is found for three new species of ground snakes, which are described here based on their unique combination of molecular, meristic, and color pattern characteristics. The name A.arangoi Prado, 1939 is revalidated for a Colombian snake species previously subsumed under A.major Boulenger, 1894 based on new material collected in Ecuador. Reidentifications are provided for Atractus voucher specimens and sequences deposited in GenBank. With these changes, the number of Atractus reported in Ecuador increases from 27 to 31 species. Finally, attention is given to the importance of using a biogeographical framework that includes molecular data and a comprehensive geographic sampling when proposing species limits in complex taxonomic groups.

Bites by exotic snakes reported to the UK National Poisons Information Service 2009-2020.

Introduction: Snakebite is recognised as a neglected tropical disease and a cause of substantial morbidity and mortality. Whilst the most medically important snakes are typically native of Asia, Africa, Latin America and Oceania, the possibility of encountering these snakes is no longer limited by geography due to an increasing number of exotic (non-native) snakes being held in captivity.Methods: A retrospective review of snakebite enquiries to the UK National Poisons Information Service (NPIS) between 2009 and 2020. Enquiries about the European adder (Vipera berus) or where the identity of the snake was unknown were excluded.Results: There were 321 exotic snakebites in 300 patients involving 68 different species during this period. Ten patients were bitten on more than one occasion. The majority of patients (64.5%) were male. Most bites were inflicted by snakes of the family Colubridae (184/321, 57.3%); seventeen bites resulted in moderate symptoms (predominantly swelling of the bitten limb). There were 30 (9.3%) bites by Viperidae and 14 (4.3%) bites by Elapidae. All severe cases (n = 15) resulted from bites by either Viperidae (n = 10) or Elapidae (n = 5). Antivenom was given in 17 cases. One fatality was recorded.Conclusions: Despite their low incidence, exotic snakebites present a substantial challenge for UK healthcare professionals. Although rare, these bites typically occur in individuals (usually male) who keep snakes as part of their occupation or hobby and are therefore at risk of multiple bites. Bites can result in venom hypersensitisation and the risk of venom-induced anaphylaxis. Rapid access to expert clinical advice is available in the UK on a 24-hour basis through the National Poisons Information Service and is strongly recommended in all cases of exotic snakebite.

Measuring Annual Variation in Reproductive Output Reveals a Key Role of Maternal Body Condition in Determining the Size of Eggs in Snakes.

Long-term studies are especially suited for disentangling the effects of extrinsic and intrinsic factors on both total reproductive investment and reproductive allocation in offspring number versus offspring size. Female reproductive traits of the red-banded wolf snake (Lycodon rufozonatus) from Zhejiang, East China were studied in four years between 1999 and 2014. Egg-laying dates overall extended from late June to late July, and varied among years. Postpartum body mass, clutch size, clutch mass, and egg size were positively related to female size (snout vent length, SVL) in each year. Postpartum body mass, clutch mass, and egg size differed among years after accounting for female SVL, whereas clutch size did not. Setting female SVL at the same level, postpartum body mass was greater in 2010 than in 2014, clutch mass was greater in 2010 than in 2011 and 2014, and egg size was greater in 2010 than in the other three years. Females did not trade off egg size against number. Egg size was positively related to postpartum body condition in each year. Females laid larger eggs in 2010 than in other three years after removing the influence of maternal body condition. Our study provides evidence for the traditional view that reproductive output is highly linked to maternal body size in snakes, but not following Smith and Fretwell's (1974) classic prediction that females with different amounts of resources to invest in reproduction should give priority to adjusting the number rather than size of their offspring. Maternal body size and condition both are important sources of variation in egg size, but factors other than these two variables may also affect the size of eggs produced by female L. rufozonatus.

Complete mitochondrial genome of the Red Keelback (Pseudagkistrodon rudis Boulenger, 1906).

The Red Keelback (Pseudagkistrodon rudis Boulenger, 1906) is widely distributed in the southern of China. The complete mitochondrial genome (mitogenome) of P. rudis was determined for the first time by using next-generation sequencing. The size of assembled mitogenome for P. rudis was 19,150 bp, which included 13 protein coding genes (PCGs), 22 tRNAs, two rRNAs and two control regions (d-loop1 and d-loop2). The Bayesian tree showed that P. rudis and Rhabdophis tigrinus have a closed relationship. These results can provide data for phylogeny and molecular classification of the genus.

First insights into the biochemical and toxicological characterization of venom from the Banded Cat-eyed Snake Leptodeira annulata pulchriceps.

With the aim to widen the current knowledge of toxinological implications of bites from rear-fanged snakes and biological roles of their venoms, this study focuses on the biochemical composition and toxic effects of the venom of Leptodeira annulata pulchriceps from Argentina. We analyzed the protein composition by electrophoresis and mass spectrometry, and enzymatic properties by quantitative assays on different substrates. Additionally, we evaluated local and systemic toxicity in mice, and tested its cross-reactivity with elapid and viperid antivenoms used in Argentina. This venom showed features reminiscent of venoms from snakes of Bothrops genus, containing components ranging from ~17 to 75 kDa, which are mainly tissue-damaging toxins such as proteinases. Although showing low lethality to mice (LD50 = 20 µg/g body weight), prominent hemorrhage developed locally in mice intramuscularly and intradermally injected with the venom, and the minimum hemorrhagic dose was found to be 12.7 µg/mouse. This study is the first comprehensive investigation of the venom of L. a. pulchriceps, and sheds new light on differences between this and those of the other two subspecies of L. annulata. Additionally, the study provides new insights into the venom components of "colubrid" snakes, advocating for considering bites from this rich diversity of snakes as a public health problem that needs to be addressed worldwide.

The complete mitochondrial genome of Plagiopholis styani (Reptilia: Squamata: Colubridae).

In this study, the complete mitochondrial genome of Plagiopholis styani was acquired and described. The mitogenome is a circular molecule of 19,669 base pairs (bp) in length including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and two control regions (D-loop), which is similar to other Colubridae snakes. The total base composition of mitochondrial DNA is A 33.0%, C 24.9%, G 12.8%, and T 29.3%. The mitochondrial genome of Plagiopholis styani contributes to revealing the phylogenetic relationships among species of the Colubridae family.

The first complete mitochondrial genome of the Grey-bellied Green Ratsnake Gonyosoma frenatum (Gray, 1853) and its phylogenetic placement.

We sequenced the complete 17,209 bp mitochondrial genome (mitogenome) of Gonyosoma frenatum (Squamata: Colubridae) using next-generation sequencing. It consists of 13 PCGs, two ribosomal RNA genes, 22 transfer RNA genes, one non-coding region of an L-strand replication origin, and two control regions. The overall nucleotide composition was 34.7% of A, 24.8% of T, 12.3% of G, and 28.1% of C. The result of the phylogenetic analysis showed that G. frenatum, a member of Colubridae, is sister to other New World and Old World ratsnakes. The new data could help better understand the phylogenetic status of the genus Gonyosoma and the evolutionary history of Colubridae species.

Using Google Trends to Determine Current, Past, and Future Trends in the Reptile Pet Trade.

Reptiles are one of the most popular exotic pets in the world, with over a third of all described species currently being traded. However, the most commonly available reptiles are typically non-threatened, captive-bred, and/or domestically obtained, which means they are also largely unregulated and unmonitored, resulting in a large portion of the reptile pet trade remaining unknown. In this study, the past, current, and future trends of the most popular reptiles in the pet trade were examined. Google Trends was used to determine the global popularity of the most popular pets from 2004 to 2020 and compared to the results from an online survey sent to individuals involved in the reptile trade. The most popular pets from the previous five years were also compared globally across regions and countries. The results determined that the most popular reptile species during the last decade is by far bearded dragons, followed by ball pythons and leopard geckos. Although the survey results were similar when asked what the top reptiles were, most respondents named ball pythons as the most popular reptile. However, when asked what reptiles had lost the most popularity during the previous decade, the survey respondents named green iguanas, Burmese pythons, chameleons, red-eared sliders, and green anoles, concurring with what was found with Google Trends. The reptiles thought to be more popular in the upcoming decade by the survey participants were blue-tongued skinks, tegus, uromastyx, crested geckos, and ball pythons-most of which did indeed show an increase in popularity during the last decade, as indicated with Google Trends. The results from Google Trends demonstrated that ball pythons and crested geckos have increased their popularity more than any other reptile in the last two decades. Reptile popularity also differed between countries, with bearded dragons the most popular reptile in Australia, Western Europe, the U.S., and Canada. Leopard geckos were the most popular reptile in Italy and Turkey, and ball pythons were the reptile of choice in Mexico, Indonesia, and India. The general finding of this study is that the reptiles declining in popularity were mostly wild-caught or restricted due to regulations, while current and future species were captive-bred and available in many varieties or morphs. The most popular species were also docile, medium-sized, and easy to handle, with relatively simple care requirements. This study demonstrates that Google Trends can be a useful tool for determining relative popularity among reptiles, or any other pet group, with results closely mirroring those obtained through direct surveying of people involved in the pet trade. However, unlike surveys, this analysis is quick, quantifiable, and can show what is popular and in-demand not only at the global level but at much finer scales. Thus, Google Trends can be a valuable tool in many research applications, especially in topics that may otherwise be difficult to monitor and quantify.