How to estimate body condition in large lizards? Argentine black and white tegu (Salvator merianae, Duméril and Bibron, 1839) as a case study

Body condition is a measure of the health and fitness of an organism represented by available energy stores, typically fat. Direct measurements of fat are difficult to obtain non-invasively, thus body condition is usually estimated by calculating body condition indices (BCIs) using mass and length. The utility of BCIs is contingent on the relationship of BCIs and fat, thereby validation studies should be performed to select the best performing BCI before application in ecological investigations. We evaluated 11 BCIs in 883 Argentine black and white tegus (Salvator merianae) removed from their non-native range in South Florida, United States. Because the length-mass relationship in tegus is allometric, a segmented linear regression model was fit to the relationship between mass and length to define size classes. We evaluated percent, residual, and scaled fat and determined percent fat was the best measure of fat, because it was the least-associated with snout-vent length (SVL). We evaluated performance of BCIs with the full dataset and within size classes and identified Fulton's K as the best performing BCI for our sampled population, explaining up to 19% of the variation in fat content. Overall, we found that BCIs: 1) maintained relatively weak relationships with measures of fat and 2) splitting data into size classes reduced the strength of the relationship (i.e., bias) between percent fat and SVL but did not improve the performance of BCIs. We postulate that the weak performance of BCIs in our dataset was likely due to the weak association of fat with SVL, the body plan and life-history traits of tegus, and potentially inadequate accounting of available energy resources. We caution against assuming that BCIs are strong indicators of body condition across species and suggest that validation studies be implemented, or that alternative or complimentary measures of health or fitness should be considered.

Burmese pythons in Florida: A synthesis of biology, impacts, and management tools

Burmese pythons (Python molurus bivittatus) are native to southeastern Asia, however, there is an established invasive population inhabiting much of southern Florida throughout the Greater Everglades Ecosystem. Pythons have severely impacted native species and ecosystems in Florida and represent one of the most intractable invasive-species management issues across the globe. The difficulty stems from a unique combination of inaccessible habitat and the cryptic and resilient nature of pythons that thrive in the subtropical environment of southern Florida, rendering them extremely challenging to detect. Here we provide a comprehensive review and synthesis of the science relevant to managing invasive Burmese pythons. We describe existing control tools and review challenges to productive research, identifying key knowledge gaps that would improve future research and decision making for python control.

Size distribution and reproductive phenology of the invasive Burmese python (Python molurus bivittatus) in the Greater Everglades Ecosystem, Florida, USA

The design of successful invasive species control programs is often hindered by the absence of basic demographic data on the targeted population. Establishment of invasive Burmese pythons (Python molurus bivittatus) in the Greater Everglades Ecosystem, Florida USA has led to local precipitous declines (> 90%) of mesomammal populations and is also a major threat to native populations of reptiles and birds. Efforts to control this species are ongoing but are hampered by the lack of access to and information on the expected biological patterns of pythons in southern Florida. We present data from more than 4,000 wild Burmese pythons that were removed in southern Florida over 26 years (1995–2021), the most robust dataset representing this invasive population to date. We used these data to characterize Burmese python size distribution, size at maturity, clutch size, and seasonal demographic and reproductive trends. We broadened the previously described size ranges by sex and, based on our newly defined size-stage classes, showed that males are smaller than females at sexual maturity, confirmed a positive correlation between maternal body size and potential clutch size, and developed predictive equations to facilitate demographic predictions. We also refined the annual breeding season (approx.100 days December into March), oviposition timing (May), and hatchling emergence and dispersal period (July through October) using correlations of capture morphometrics with observations of seasonal gonadal recrudescence (resurgence) and regression. Determination of reproductive output and timing can inform population models and help managers arrest population growth by targeting key aspects of python life history. These results define characteristics of the species in Florida and provide an enhanced understanding of the ecology and reproductive biology of Burmese pythons in their invasive Everglades range.

Clutch may predict growth of hatchling Burmese pythons better than food availability or sex

Identifying which environmental and genetic factors affect growth pattern phenotypes can help biologists predict how organisms distribute finite energy resources in response to varying environmental conditions and physiological states. This information may be useful for monitoring and managing populations of cryptic, endangered, and invasive species. Consequently, we assessed the effects of food availability, clutch, and sex on the growth of invasive Burmese pythons (Python bivittatus Kuhl) from the Greater Everglades Ecosystem in Florida, USA. Though little is known from the wild, Burmese pythons have been physiological model organisms for decades, with most experimental research sourcing individuals from the pet trade. Here, we used 60 hatchlings collected as eggs from the nests of two wild pythons, assigned them to High or Low feeding treatments, and monitored growth and meal consumption for 12 weeks, a period when pythons are thought to grow very rapidly. None of the 30 hatchlings that were offered food prior to their fourth week post-hatching consumed it, presumably because they were relying on internal yolk stores. Although only two clutches were used in the experiment, we found that nearly all phenotypic variation was explained by clutch rather than feeding treatment or sex. Hatchlings from clutch 1 (C1) grew faster and were longer, heavier, in better body condition, ate more frequently, and were bolder than hatchlings from clutch 2 (C2), regardless of food availability. On average, C1 and C2 hatchling snout-vent length (SVL) and weight grew 0.15 cm d⁻¹ and 0.10 cm d⁻¹, and 0.20 g d⁻¹ and 0.03 g d⁻¹, respectively. Additional research may be warranted to determine whether these effects remain with larger clutch sample sizes and to identify the underlying mechanisms and fitness implications of this variation to help inform risk assessments and management.

This article has an associated First Person interview with the first author of the paper.

Summary: Hatchling pythons from an invasive population displayed substantial phenotypic variation in morphometrics, growth rates, and behaviors. This information may be useful for managing populations of cryptic, endangered, and invasive species.

Host-specific phenotypic variation of a parasite co-introduced with invasive Burmese pythons

Invasive Burmese pythons (Python bivittatus Kuhl, 1820) have introduced a lung parasite, Raillietiella orientalis, (Hett, 1915) from the python's native range in Southeast Asia to its introduced range in Florida, where parasite spillover from pythons to two families and eight genera of native snakes has occurred. Because these novel host species present a diversity of ecological and morphological traits, and because these parasites attach to their hosts with hooks located on their cephalothorax, we predicted that R. orientalis would exhibit substantial, host-associated phenotypic plasticity in cephalothorax shape. Indeed, geometric morphometric analyses of 39 parasites from five host species revealed significant variation among host taxa in R. orientalis cephalothorax shape. We observed differences associated with host ecology, where parasites from semi-aquatic and aquatic snakes exhibited the greatest morphological similarity. Morphological analyses of R. orientalis recovered from invasive pythons, native pit vipers, and terrestrial snakes each revealed distinct shapes. Our results suggest R. orientalis can exhibit significant differences in morphology based upon host species infected, and this plasticity may facilitate infection with this non-native parasite in a wide array of novel squamate host species.

Parasite spillover: indirect effects of invasive Burmese pythons

Identification of the origin of parasites of nonindigenous species (NIS) can be complex. NIS may introduce parasites from their native range and acquire parasites from within their invaded range. Determination of whether parasites are non-native or native can be complicated when parasite genera occur within both the NIS' native range and its introduced range. We explored potential for spillover and spillback of lung parasites infecting Burmese pythons (Python bivittatus) in their invasive range (Florida). We collected 498 indigenous snakes of 26 species and 805 Burmese pythons during 2004-2016 and examined them for lung parasites. We used morphology to identify three genera of pentastome parasites, Raillietiella, a cosmopolitan form, and Porocephalus and Kiricephalus, both New World forms. We sequenced these parasites at one mitochondrial and one nuclear locus and showed that each genus is represented by a single species, R. orientalis, P. crotali, and K. coarctatus. Pythons are host to R. orientalis and P. crotali, but not K. coarctatus; native snakes are host to all three species. Sequence data show that pythons introduced R. orientalis to North America, where this parasite now infects native snakes. Additionally, our data suggest that pythons are competent hosts to P. crotali, a widespread parasite native to North and South America that was previously hypothesized to infect only viperid snakes. Our results indicate invasive Burmese pythons have affected parasite-host dynamics of native snakes in ways that are consistent with parasite spillover and demonstrate the potential for indirect effects during invasions. Additionally, we show that pythons have acquired a parasite native to their introduced range, which is the initial condition necessary for parasite spillback.

Challenges to a molecular approach to prey identification in the Burmese python, Python molurus bivittatus

Molecular approaches to prey identification are increasingly useful in elucidating predator-prey relationships, and we aimed to investigate the feasibility of these methods to document the species identities of prey consumed by invasive Burmese pythons in Florida. We were particularly interested in the diet of young snakes, because visual identification of prey from this size class has proven difficult. We successfully extracted DNA from the gastrointestinal contents of 43 young pythons, as well as from several control samples, and attempted amplification of DNA mini-barcodes, a 130-bp region of COX1. Using a PNA clamp to exclude python DNA, we found that prey DNA was not present in sufficient quality for amplification of this locus in 86% of our samples. All samples from the GI tracts of young pythons contained only hair, and the six samples we were able to identify to species were hispid cotton rats. This suggests that young Burmese pythons prey predominantly on small mammals and that prey diversity among snakes of this size class is low. We discuss prolonged gastrointestinal transit times and extreme gastric breakdown as possible causes of DNA degradation that limit the success of a molecular approach to prey identification in Burmese pythons.

Clonal reproduction shapes evolution in the lizard malaria parasite Plasmodium floridense

The preponderant clonal evolution hypothesis (PCE) predicts that frequent clonal reproduction (sex between two clones) in many pathogens capable of sexual recombination results in strong linkage disequilibrium and the presence of discrete genetic subdivisions characterized by occasional gene flow. We expand on the PCE and predict that higher rates of clonal reproduction will result in: (1) morphologically cryptic species that exhibit (2) low within-species variation and (3) recent between-species divergence. We tested these predictions in the Caribbean lizard malaria parasite Plasmodium floridense using 63 single-infection samples in lizards collected from across the parasite's range, and sequenced them at two mitochondrial, one apicoplast, and five nuclear genes. We identified 11 provisionally cryptic species within P. floridense, each of which exhibits low intraspecific variation and recent divergence times between species (some diverged approximately 110,000 years ago). Our results are consistent with the hypothesis that clonal reproduction can profoundly affect diversification of species capable of sexual recombination, and suggest that clonal reproduction may have led to a large number of unrecognized pathogen species. The factors that may influence the rates of clonal reproduction among pathogens are unclear, and we discuss how prevalence and virulence may relate to clonal reproduction.

Tree-based delimitation of morphologically ambiguous taxa: a study of the lizard malaria parasites on the Caribbean island of Hispaniola

Malaria parasites in the genus Plasmodium have been classified primarily on the basis of differences in morphology. These single-celled organisms often lack distinguishing morphological features, and this can encumber both species delimitation and identification. Six saurian malaria parasites have been described from the Caribbean island of Hispaniola. All six infect lizards in the genus Anolis, but only two of these parasites can be distinguished using morphology. The remaining four species overlap in morphology and geography, and cannot be consistently identified using traditional methods. We compared a morphological approach with a molecular phylogenetic approach for assessing the taxonomy of these parasites. We surveyed for blood parasites from 677 Anolis lizards, representing 26 Anolis spp. from a total of 52 sites across Hispaniola. Fifty-five of these lizards were infected with Plasmodium spp., representing several new host records, but only 24 of these infections could be matched to previously described species using traditional morphological criteria. We then estimated the phylogeny of these parasites using both mitochondrial (cytb and coxI) and nuclear (EF2) genes, and included carefully selected GenBank sequences to confirm identities for certain species. Our molecular results unambiguously corroborated our morphology-based species identifications for only the two species previously judged to be morphologically distinctive. The remaining infections fell into two well-supported and reciprocally monophyletic clades, which contained the morphological variation previously reported for all four of the morphologically ambiguous species. One of these clades was identified as Plasmodium floridense and the other as Plasmodium fairchildi hispaniolae. We elevate the latter to Plasmodium hispaniolae comb. nov. because it is polyphyletic with the mainland species Plasmodium fairchildifairchildi and we contribute additional morphological and molecular characters for future species delimitation. Our phylogenetic hypotheses indicate that two currently recognised taxa, Plasmodium minasense anolisi and Plasmodium tropiduri caribbense, are not valid on Hispaniola. These results illustrate that molecular data can improve taxonomic hypotheses in Plasmodium when reliable morphological characters are lacking.