Urbanization may limit impacts of an invasive predator on native mammal diversity

Biotic homogenisation and differentiation as directional change in beta diversity: synthesising driver-response relationships to develop conceptual models across ecosystems

Biotic homogenisation is defined as decreasing dissimilarity among ecological assemblages sampled within a given spatial area over time. Biotic differentiation, in turn, is defined as increasing dissimilarity over time. Overall, changes in the spatial dissimilarities among assemblages (termed 'beta diversity') is an increasingly recognised feature of broader biodiversity change in the Anthropocene. Empirical evidence of biotic homogenisation and biotic differentiation remains scattered across different ecosystems. Most meta-analyses quantify the prevalence and direction of change in beta diversity, rather than attempting to identify underlying ecological drivers of such changes. By conceptualising the mechanisms that contribute to decreasing or increasing dissimilarity in the composition of ecological assemblages across space, environmental managers and conservation practitioners can make informed decisions about what interventions may be required to sustain biodiversity and can predict potential biodiversity outcomes of future disturbances. We systematically reviewed and synthesised published empirical evidence for ecological drivers of biotic homogenisation and differentiation across terrestrial, marine, and freshwater realms to derive conceptual models that explain changes in spatial beta diversity. We pursued five key themes in our review: (i) temporal environmental change; (ii) disturbance regime; (iii) connectivity alteration and species redistribution; (iv) habitat change; and (v) biotic and trophic interactions. Our first conceptual model highlights how biotic homogenisation and differentiation can occur as a function of changes in local (alpha) diversity or regional (gamma) diversity, independently of species invasions and losses due to changes in species occurrence among assemblages. Second, the direction and magnitude of change in beta diversity depends on the interaction between spatial variation (patchiness) and temporal variation (synchronicity) of disturbance events. Third, in the context of connectivity and species redistribution, divergent beta diversity outcomes occur as different species have different dispersal characteristics, and the magnitude of beta diversity change associated with species invasions also depends strongly on alpha and gamma diversity prior to species invasion. Fourth, beta diversity is positively linked with spatial environmental variability, such that biotic homogenisation and differentiation occur when environmental heterogeneity decreases or increases, respectively. Fifth, species interactions can influence beta diversity via habitat modification, disease, consumption (trophic dynamics), competition, and by altering ecosystem productivity. Our synthesis highlights the multitude of mechanisms that cause assemblages to be more or less spatially similar in composition (taxonomically, functionally, phylogenetically) through time. We consider that future studies should aim to enhance our collective understanding of ecological systems by clarifying the underlying mechanisms driving homogenisation or differentiation, rather than focusing only on reporting the prevalence and direction of change in beta diversity, per se.

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.

Effect of urbanization on the opossum Didelphis virginiana health and implications for zoonotic diseases

Urban animals can be an important threat to human health as possible hosts of zoonotic diseases and their susceptibility to these diseases can depend on their overall health conditions. Thus, it is important to understand the factors that determine their health conditions. For this, we studied Virginia opossums (Didelphis virginiana) in six locations with different urbanization levels and types in Mexico City, Mexico. We trapped opossums and measured eight health-related characteristics (number of ectoparasites and tartar severity, among others) and estimated the percentage of area covered by the four main types of terrain (natural vegetation, managed vegetation, impermeable terrain and constructions). Data were analyzed by a canonical correspondence analysis. We found that impermeable terrain was related to negative health characteristics, while the constructions were opposite to impermeable terrain and mostly related to good health characteristics. At the same, time constructed areas held a smaller population than the natural areas. This suggests that constructed areas provide few shelters, but opossums are healthier there, while impermeable areas provide more shelter but cause more health problems to the animals, thus increasing the chances of zoonotic diseases. We recommend reducing the impermeable areas in city planning to contribute to a better health of the urban animals and therefore reduce risks of zoonotic diseases with potentially disastrous results.

Invasive snake causes massive reduction of all endemic herpetofauna on Gran Canaria

Invasive snakes represent a serious threat to island biodiversity, being responsible for far-reaching impacts that are noticeably understudied, particularly regarding native reptiles. We analysed the impact of the invasive California kingsnake, Lampropeltis californiae-recently introduced in the Canary Islands-on the abundance of all endemic herpetofauna of the island of Gran Canaria. We quantified the density in invaded and uninvaded sites for the Gran Canaria giant lizard, Gallotia stehlini, the Gran Canaria skink, Chalcides sexlineatus, and Boettger's wall gecko, Tarentola boettgeri. We used spatially explicit capture-recapture and distance-sampling methods for G. stehlini and active searches under rocks for the abundance of the other two reptiles. The abundance of all species was lower in invaded sites, with a reduction in the number of individuals greater than 90% for G. stehlini, greater than 80% for C. sexlineatus and greater than 50% for T. boettgeri in invaded sites. Our results illustrate the severe impact of L. californiae on the endemic herpetofauna of Gran Canaria and highlight the need for strengthened measures to manage this invasion. We also provide further evidence of the negative consequences of invasive snakes on island reptiles and emphasize the need for further research on this matter on islands worldwide.

Niche lability mitigates the impact of invasion but not urbanization

Native species can coexist with invasive congeners by partitioning niche space; however, impacts from invasive species often occur alongside other disturbances. Native species' responses to the interactions of multiple disturbances remain poorly understood. Here we study the impacts of urbanization and an invasive congener on a native species. Using abundance (catch-per-unit effort) and vertical distribution of native green anoles (Anolis carolinensis) and invasive brown anoles (Anolis sagrei) across a gradient of natural-to-urban forests, we ask if niche shifting (lability) is occurring, and if it can mitigate impacts from one or both disturbances. We use generalized linear models to relate species abundances across the landscape to urbanization, forest structural complexity, and congener abundances (i.e., A. sagrei); and test for an interaction between urbanization and congener abundance. Our data show that A. sagrei presence results in a 17-fold upward shift in vertical niche of A. carolinensis-an 8.3 m shift in median perch height, and models reveal urbanization also drives an increase in A. carolinensis perch height. A. carolinensis and A. sagrei abundances negatively and positively correlate with urbanization, respectively, and neither species' abundance correlate with congener abundance. Despite a positive correlation between A. sagrei abundance and urbanization, our results do not show evidence of this interaction affecting A. carolinensis. Instead, niche lability appears to enable the native species to mitigate the impact of one driver of decline (invasive competition) while our data suggest it declines with the second (urbanization).

Invasive Burmese pythons alter host use and virus infection in the vector of a zoonotic virus

The composition of wildlife communities can have strong effects on transmission of zoonotic vector-borne pathogens, with more diverse communities often supporting lower infection prevalence in vectors (dilution effect). The introduced Burmese python, Python bivittatus, is eliminating large and medium-sized mammals throughout southern Florida, USA, impacting local communities and the ecology of zoonotic pathogens. We investigated invasive predator-mediated impacts on ecology of Everglades virus (EVEV), a zoonotic pathogen endemic to Florida that circulates in mosquito-rodent cycle. Using binomial generalized linear mixed effects models of field data at areas of high and low python densities, we show that increasing diversity of dilution host (non-rodent mammals) is associated with decreasing blood meals on amplifying hosts (cotton rats), and that increasing cotton rat host use is associated with increasing EVEV infection in vector mosquitoes. The Burmese python has caused a dramatic decrease in mammal diversity in southern Florida, which has shifted vector host use towards EVEV amplifying hosts (rodents), resulting in an indirect increase in EVEV infection prevalence in vector mosquitoes, putatively elevating human transmission risk. Our results indicate that an invasive predator can impact wildlife communities in ways that indirectly affect human health, highlighting the need for conserving biological diversity and natural communities.

Environmental DNA sampling reveals high occupancy rates of invasive Burmese pythons at wading bird breeding aggregations in the central Everglades

The Burmese python (Python bivittatus) is now established as a breeding population throughout south Florida, USA. However, the extent of the invasion, and the ecological impacts of this novel apex predator on animal communities are incompletely known, in large part because Burmese pythons (hereafter “pythons”) are extremely cryptic and there has been no efficient way to detect them. Pythons are recently confirmed nest predators of long-legged wading bird breeding colonies (orders Ciconiiformes and Pelecaniformes). Pythons can consume large quantities of prey and may not be recognized as predators by wading birds, therefore they could be a particular threat to colonies. To quantify python occupancy rates at tree islands where wading birds breed, we utilized environmental DNA (eDNA) analysis—a genetic tool which detects shed DNA in water samples and provides high detection probabilities. We fitted multi-scale Bayesian occupancy models to test the prediction that pythons occupy islands with wading bird colonies at higher rates compared to representative control islands containing no breeding birds. Our results suggest that pythons are widely distributed across the central Everglades in proximity to active wading bird colonies. In support of our prediction that pythons are attracted to colonies, site-level python eDNA occupancy rates were higher at wading bird colonies (ψ = 0.88, 95% credible interval [0.59–1.00]) than at the control islands (ψ = 0.42 [0.16–0.80]) in April through June (n = 15 colony-control pairs). We found our water temperature proxy (time of day) to be informative of detection probability, in accordance with other studies demonstrating an effect of temperature on eDNA degradation in occupied samples. Individual sample concentrations ranged from 0.26 to 38.29 copies/μL and we generally detected higher concentrations of python eDNA in colony sites. Continued monitoring of wading bird colonies is warranted to determine the effect pythons are having on populations and investigate putative management activities.

Cytonuclear discordance in the Florida Everglades invasive Burmese python (Python bivittatus) population reveals possible hybridization with the Indian python (P. molurus)

The invasive Burmese python (Python bivittatus) has been reproducing in the Florida Everglades since the 1980s. These giant constrictor snakes have caused a precipitous decline in small mammal populations in southern Florida following escapes or releases from the commercial pet trade. To better understand the invasion pathway and genetic composition of the population, two mitochondrial (mtDNA) loci across 1,398 base pairs were sequenced on 426 snakes and 22 microsatellites were assessed on 389 snakes. Concatenated mtDNA sequences produced six haplotypes with an average nucleotide and haplotype diversity of π = 0.002 and h = 0.097, respectively. Samples collected in Florida from morphologically identified P. bivittatus snakes were similar to published cytochrome oxidase 1 and cytochrome b sequences from both P. bivittatus and Python molurus and were highly divergent (genetic distances of 5.4% and 4.3%, respectively). The average number of microsatellite alleles and expected heterozygosity were N A = 5.50 and H E = 0.60, respectively. Nuclear Bayesian assignment tests supported two genetically distinct groups and an admixed group, not geographically differentiated. The effective population size (N E = 315.1) was lower than expected for a population this large, but reflected the low genetic diversity overall. The patterns of genetic diversity between mtDNA and microsatellites were disparate, indicating nuclear introgression of separate mtDNA lineages corresponding to cytonuclear discordance. The introgression likely occurred prior to the invasion, but genetic information on the native range and commercial trade is needed for verification. Our finding that the Florida python population is comprised of distinct lineages suggests greater standing variation for adaptation and the potential for broader areas of suitable habitat in the invaded range.