Habitat fragmentation effects on annual survival of the federally protected eastern indigo snake

Combining species distribution models and big datasets may provide finer assessments of snakebite impacts

BACKGROUND

Snakebite is a major poverty-related neglected tropical disease. An integrated scientific approach is needed to understand the dynamics of this important health issue. Our objective was to estimate snakebite occurrence in a tropical area by using a blend of ecological modelling and robust statistical analysis.

METHODS

The present study used climatic, environmental, and human population density data to determine the area with snakebite occurrence-probability for the first time in Bangladesh. We also analysed a large, 16-year dataset of hospitalized snakebite cases to reveal the epidemiology of snakebite in the south-eastern zone of the country.

FINDINGS

Our results show that cobra bite is the most commonly occurring venomous snakebite in humans (around ~12% of the total yearly snakebite records), and men are more frequently bitten than women (2/3 of human victims are men). Most bites occur during the rainy season for cobra and green pit viper, while krait bites are not restricted to any particular season. As snakebite incidents are closely related to climate conditions, we can model snakebite risk using temperature and precipitation variables. Whereas there is a lack of snakebite reports from several parts of the study area in official records, our models predict that the entire study area is favourable for snakebite incidents. Based on the combined evidence we estimate that about 200,000 snakebite events occur every year in the south-eastern part of Bangladesh alone. Considering future global climate change, our model projections show that snakebite incidence in Bangladesh might not significantly decrease in the future (- 2070-); however, the distribution of probabilities might change, with a predicted increase of snakebite incidence in the hilly areas of the country.

CONCLUSIONS

Using climatic data to predict snakebite incidence in Bangladesh allowed us to provide estimations of the total annual number of snakebite cases in the study area. As in most countries, the scarcity of accurate epidemiological data in official records might have masked the real magnitude of this problem. Our analysis suggests that the problem of snakebite envenoming in Bangladesh might be worse than currently perceived. A long-term sustainable snakebite program plan should be designed and institutionalized, considering climatic, geographical and human demographic variables, to obtain better data and facilitate the implementation of accurate snakebite management programs for this country.

Health Assessment of Free-Ranging Eastern Indigo Snakes (Drymarchon couperi) from Hydrologic Restoration Construction Sites in South Florida, USA

There is a paucity of information regarding the health status of free-ranging eastern indigo snakes (EIS; Drymarchon couperi) in heavily modified and developing landscapes. As a component of regional Florida Everglades restoration efforts, several areas occupied by EIS are being converted from agricultural lands to reservoirs. From 2020 to 2022, 28 EIS were opportunistically captured at two of these sites and brought into captivity to join a captive breeding colony; however, 11 snakes died within 5 mo of capture. Health assessments were performed on 28 individuals and included hematology and plasma biochemistry analysis, as well as screening for pesticide contaminant levels, parasites, and other pathogens. Overall, the presence of pathogens was relatively high, suggesting immunosuppression secondary to stress: 25/28 (89.4%) Kalicephalus sp.; 12/28 (42.9%) Raillietiella orientalis; 11/28 (39.2%) Ochetosoma validum; 7/28 (25.0%) Cryptosporidium serpentis; 3/28 (10.7%) snake adenovirus 1; and 1/28 (3.6%) Ferlavirus genotype C. Stress may have been caused by physical displacement, habitat modification, and noise pollution. These potential stressors (including the presence of remnant harmful chemicals from previous land use and the impacts on this federally threatened species) should be considered further when making restoration or construction decisions.

Favourable climatic niche in low elevations outside the flood zone characterises the distribution pattern of venomous snakes in Bangladesh

Snakes are sensitive to both environmental and climate gradients. To design conservation plans, a scientific understanding of snake habitats in light of environmental and climatic variables is an essential prerequisite. For venomous snakes, denoting favourable habitats should also be relevant for snakebite management. We have considered 18 spatial variables to portray the range of terrestrial venomous snake distribution in Bangladesh. Our results indicate that the distribution of 29 studied venomous snakes in this country is primarily driven by climatic and environmental variables. We found that especially low elevation and flood risk constrain the distribution of those terrestrial snakes, i.e. regular floods in central Bangladesh push venomous snakes towards the edges of the country. Moreover, none of these species occupies the whole of its anticipated climatically favourable area. Projections into the future indicated that 11 studied species, Amphiesma platyceps, Boiga siamensis, Chrysopelea ornata, Pseudoxenodon macrops, Rhabdophis himalayanus, Rhabdophis subminiatus, Bungarus lividus, Ophiophagus hannah, Daboia russelii, Ovophis monticola and Trimeresurus popeiorum will lose their entire climatically suitable area within the country. Therefore, we suggest establishing more protected areas in the hilly ecosystems in the eastern part and in the mangrove forests in the south-western corner of Bangladesh to mitigate future extinction risks, such as climate change, sea-level rise and increase in flood severity. Conserving village forests and croplands, which are subject to rapid change, will also need to be addressed equally, as these are inhabited by almost one-third of the studied species. The occurrence of the cobras and kraits in village forests and cropland dominant habitats demands more attention to minimise snakebite related mortality and morbidity.

Accounting for geographic variation in species-habitat associations during habitat suitability modeling

Range-wide species conservation efforts are facilitated by spatially explicit estimates of habitat suitability. However, species-environment relationships often vary geographically and models assuming geographically constant relationships may result in misleading inferences. We present the first range-wide habitat suitability model (HSM) for the federally threatened eastern indigo snake (Drymarchon couperi) as a case study illustrating an approach to account for known latitudinal variation in habitat associations. Specifically, we modeled habitat suitability using interactive relationships between minimum winter temperature and several a priori environmental covariates and compared our results to those from models assuming geographically constant relationships. We found that multi-scale models including interactive effects with winter temperature outperformed single-scale models and models not including interactive effects with winter temperature. Our top-ranked model had suitable range-wide predictive performance and identified numerous large (i.e., ≥1000 ha) potential habitat patches throughout the indigo snake range. Predictive performance was greatest in southern Georgia and northern Florida likely reflecting more restrictive indigo snake habitat associations in these regions. This study illustrates how modeling interactive effects between temperature and environmental covariates can improve the performance of HSMs across geographically varying environmental gradients.

Multiscale assessment of functional connectivity: Landscape genetics of eastern indigo snakes in an anthropogenically fragmented landscape in central Florida

Landscape features can strongly influence gene flow and the strength and direction of these effects may vary across spatial scales. However, few studies have evaluated methodological approaches for selecting spatial scales in landscape genetics analyses, in part because of computational challenges associated with optimizing landscape resistance surfaces (LRS). We used the federally threatened eastern indigo snake (Drymarchon couperi) in central Florida as a case study with which to compare the importance of landscape features and their scales of effect in influencing gene flow. We used genetic algorithms (ResistanceGA) to empirically optimize LRS using categorical land cover surfaces, multiscale resource selection surfaces (RSS), and four combinations of landscape covariates measured at multiple spatial scales (multisurface multiscale LRS). We compared LRS where scale was selected using pseudo- and full optimization. Multisurface multiscale LRS received more empirical support than LRS optimized from categorical land cover surfaces or RSS. Multiscale LRS with scale selected using full optimization generally outperformed those with scale selected using pseudo-optimization. Multiscale LRS with large spatial scales (1200-1800 m) received the most empirical support. Our results highlight the importance of considering landscape features across multiple spatial scales in landscape genetic analyses, particularly broad scales relative to species movement potential. Different effects of scale on home range-level movements and dispersal could explain weak associations between habitat suitability and gene flow in other studies. Our results also demonstrate the importance of large tracts of undeveloped upland habitat with heterogenous vegetation communities and low urbanization for promoting indigo snake connectivity.

Prioritization of Vulnerable Species Under Scenarios of Anthropogenic-Driven Change in Georgia's Coastal Plain

Effective management of wildlife populations benefits from an understanding of the long-term vulnerability of species to anthropogenic stressors. Exposure to potential habitat change is one measure of vulnerability that wildlife managers often use to assess and prioritize individual species or groups of species for resource allocation or direct management actions. We used species distribution models for 15 species occurring in the coastal plain ecoregion of Georgia to estimate the current amount and distribution of potential habitat and then predict exposure to changes in habitat due to inundation from sea level rise (using the Sea Level Affecting Marshes model) and urban growth (using the Slope Land-use Excluded Urban Topology Hillshade Growth model) for four future time points. Our results predict that all focal species were likely to experience some exposure to habitat change from either sea level rise or urbanization, but few species will experience high exposure to change from both stressors. Species that use salt marsh or beach habitats had the highest predicted exposure from sea level rise (25–69%), while species that use more inland habitats had the highest predicted exposure to urban growth (10–20%). Our models are a resource for managers considering tradeoffs between prioritization schemes under two future stressors. Results suggest that managers may need to prioritize species (or their habitats) based on the predicted magnitude of habitat loss, while also contextualizing prioritization with respect to the current amount of available protected habitat and species global vulnerability.

Sunning themselves in heaps, knots, and snarls: The extraordinary abundance and demography of island watersnakes

Snakes represent a sizable fraction of vertebrate biodiversity, but until recently, data on their demography have been sparse. Consequently, generalizations regarding patterns of variation are weak and the potential for population projections is limited. We address this information gap through an analysis of spatial and temporal variation in demography (population size, annual survival, and realized population growth) of the Lake Erie Watersnake, Nerodia sipedon insularum, and a review of snake survival more generally. Our study spans a period during which the Lake Erie Watersnake was listed as threatened under the U.S. Endangered Species Act, recovered, and was delisted. We collected capture-mark-recapture data at 14 study sites over 20 years, accruing 20,000 captures of 13,800 individually marked adults. Lake Erie Watersnakes achieve extraordinary abundance, averaging 520 adults per km of shoreline (ca. 260 adult per ha) at our study sites (range = 160-1,600 adults per km; ca. 80-800 adults per ha) and surpassing population recovery and postdelisting monitoring criteria. Annual survival averages 0.68 among adult females and 0.76 among adult males, varies among sites, and is positively correlated with body size among study sites. Temporal process variance in annual survival is low, averaging 0.0011 or less than 4% of total variance; thus, stochasticity in annual survival may be of minor significance to snake extinction risk. Estimates of realized population growth indicate that population size has been stable or increasing over the course of our study. More generally, snake annual survival overlaps broadly across continents, climate zones, families, subfamilies, reproductive modes, body size categories, maturation categories, and parity categories. Differences in survival in relation to size, parity, and maturation are in the directions predicted by life history theory but are of small magnitude with much variation around median values. Overall, annual survival appears to be quite plastic, varying with food availability, habitat quality, and other ecological variables.

HEALTH AND NUTRITIONAL ASSESSMENT OF FREE-RANGING EASTERN INDIGO SNAKES (DRYMARCHON COUPERI) IN GEORGIA, UNITED STATES

Clinical pathology and nutritional parameters are useful in evaluating and monitoring threatened and endangered wildlife populations, but reference ranges for most snake species are lacking. From 2001 to 2005, health assessments were performed on 58 eastern indigo snakes (EIS) (Drymarchon couperi) captured in the wild in southeastern Georgia, United States. Health and nutritional assessments performed included hematology, serum biochemistry, fat-soluble vitamins, heavy metals, pesticide contaminants, parasitology, and surveys of other pathogens. Significant differences in total solids, packed cell volume, glucose, blood urea nitrogen, albumin : globulin ratio, amylase, triglycerides, and bile acids between males and females were observed. Additionally, there was a significant difference between liver and kidney concentrations for vitamins A and E. As previously noted in captive EIS, total Ca was elevated in comparison to concentrations reported in other snake species. Parasitism was a common finding in sampled EIS, but the overall health status of this free-ranging population appeared good. A winter-time dermatitis was found in most snakes, which resolved in the summer months. This study represents the first health and nutritional assessment of free-ranging EIS, and provides needed data to guide monitoring and conservation efforts.

Population dynamics of the critically endangered golden lancehead pitviper, Bothrops insularis: stability or decline?

Little is known about vital rates of snakes generally because of the difficulty in collecting data. Here we used a robust design mark-recapture model to estimate survival, behavioral effects on capture probability, temporary emigration, abundance and test the hypothesis of population decline in the golden lancehead pitviper, Bothrops insularis, an endemic and critically endangered species from southeastern Brazil. We collected data at irregular intervals over ten occasions from 2002 to 2010. Survival was slightly higher in the wet season than in the dry season. Temporal emigration was high, indicating the importance of accounting for this parameter both in the sampling design and modeling. No behavioral effects were detected on capture probability. We detected an average annual population decrease ( = 0.93, CI = 0.47–1.38) during the study period, but estimates included high uncertainty, and caution in interpretation is needed. We discuss the potential effects of the illegal removal of individuals and the implications of the vital rates obtained for the future persistence and conservation of this endemic, endangered species.

Coastal nurseries and their importance for conservation of sea kraits

Destruction and pollution of coral reefs threaten these marine biodiversity hot stops which shelter more than two thirds of sea snake species. Notably, in many coral reef ecosystems of the Western Pacific Ocean, large populations of sea kraits (amphibious sea snakes) have drastically declined during the past three decades. Protecting remaining healthy populations is thus essential. In New Caledonia, coral reefs shelter numerous sea krait colonies spread throughout an immense lagoon (24,000 km²). Sea kraits feed on coral fish but lay their eggs on land. However, ecological information on reproduction and juveniles is extremely fragmentary, precluding protection of key habitats for reproduction. Our 10 years mark recapture study on Yellow sea kraits (L. saintgironsi >8,700 individuals marked) revealed that most neonates aggregate in highly localized coastal sites, where they feed and grow during several months before dispersal. Hundreds of females emigrate seasonally from remote populations (>50 km away) to lay their eggs in these coastal nurseries, and then return home. Protecting these nurseries is a priority to maintain recruitment rate, and to retain sea krait populations in the future.