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.

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-1 and 0.10 cm d-1, and 0.20 g d-1 and 0.03 g d-1, 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.

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.

Modeling the distributions of tegu lizards in native and potential invasive ranges.

Invasive reptilian predators can have substantial impacts on native species and ecosystems. Tegu lizards are widely distributed in South America east of the Andes, and are popular in the international live animal trade. Two species are established in Florida (U.S.A.) - Salvator merianae (Argentine black and white tegu) and Tupinambis teguixin sensu lato (gold tegu) - and a third has been recorded there- S. rufescens (red tegu). We built species distribution models (SDMs) using 5 approaches (logistic regression, multivariate adaptive regression splines, boosted regression trees, random forest, and maximum entropy) based on data from the native ranges. We then projected these models to North America to develop hypotheses for potential tegu distributions. Our results suggest that much of the southern United States and northern México probably contains suitable habitat for one or more of these tegu species. Salvator rufescens had higher habitat suitability in semi-arid areas, whereas S. merianae and T. teguixin had higher habitat suitability in more mesic areas. We propose that Florida is not the only state where these taxa could become established, and that early detection and rapid response programs targeting tegu lizards in potentially suitable habitat elsewhere in North America could help prevent establishment and abate negative impacts on native ecosystems.

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.

A validation of 11 body-condition indices in a giant snake species that exhibits positive allometry.

Body condition is a gauge of the energy stores of an animal, and though it has important implications for fitness, survival, competition, and disease, it is difficult to measure directly. Instead, body condition is frequently estimated as a body condition index (BCI) using length and mass measurements. A desirable BCI should accurately reflect true body condition and be unbiased with respect to size (i.e., mean BCI estimates should not change across different length or mass ranges), and choosing the most-appropriate BCI is not straightforward. We evaluated 11 different BCIs in 248 Burmese pythons (Python bivittatus), organisms that, like other snakes, exhibit simple body plans well characterized by length and mass. We found that the length-mass relationship in Burmese pythons is positively allometric, where mass increases rapidly with respect to length, and this allowed us to explore the effects of allometry on BCI verification. We employed three alternative measures of 'true' body condition: percent fat, scaled fat, and residual fat. The latter two measures mostly accommodated allometry in true body condition, but percent fat did not. Our inferences of the best-performing BCIs depended heavily on our measure of true body condition, with most BCIs falling into one of two groups. The first group contained most BCIs based on ratios, and these were associated with percent fat and body length (i.e., were biased). The second group contained the scaled mass index and most of the BCIs based on linear regressions, and these were associated with both scaled and residual fat but not body length (i.e., were unbiased). Our results show that potential differences in measures of true body condition should be explored in BCI verification studies, particularly in organisms undergoing allometric growth. Furthermore, the caveats of each BCI and similarities to other BCIs are important to consider when determining which BCI is appropriate for any particular taxon.

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.

Host specificity shapes population structure of pinworm parasites in Caribbean reptiles.

Host specificity is one of the potential factors affecting parasite diversification because gene flow may be facilitated or constrained by the number of host species that a parasite can exploit. We test this hypothesis using a costructure approach, comparing two sympatric pinworm parasites that differ in host specificity - Parapharyngodon cubensis and Spauligodon anolis - on the Puerto Rican Bank and St. Croix in the Caribbean. Spauligodon anolis specializes on Anolis lizards, whereas P. cubensis parasitizes Anolis lizards as well as many other species of lizards and snakes. We collected lizards from across the Puerto Rican Bank and St. Croix, sampled them for S. anolis and P. cubensis and generated nuclear and mitochondrial sequence data from the parasites. We used these data to show that P. cubensis is comprised of multiple cryptic species that exhibit limited population structure relative to S. anolis, which is consistent with our prediction based on their host specificity. We also provide evidence that the distribution of P. cubensis species is maintained by competitive exclusion, and in contrast to previous theoretical work, the parasites with the greatest number of host species also reach the highest prevalence rates. Overall, our results are consistent with the hypothesis that host specificity shapes parasite diversification, and suggest that even moderate differences in host specificity may contribute to substantial differences in diversification.

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.