Abstract
Understanding the origin and spread of invasive species is critical for predicting when and where new introductions will establish, and impact native species. However, due to the complexity of contributing factors such as multiple introductions, dispersal method, genetic admixture in founding populations, and variable propagule pressure, genetic patterns observed in invasive species may not always conform to a single theoretical expectation. Cuban treefrogs (Osteopilus septentrionalis) are invasive in peninsular Florida and sporadically in the Florida panhandle. Though O. septentrionalis has been occasionally reported in Louisiana since the 1990s, established populations were not present until the discovery of a breeding population in New Orleans in 2017. In this study we investigated the source of this novel population using existing and newly generated cytochrome B (cyt-b) mitochondrial gene sequences from the native and invasive range of O. septentrionalis. We recovered a total of 14 cyt-b haplotypes, nine novel and five previously published. Within the 95 Louisiana invasion samples, we recovered seven haplotypes including five novel haplotypes. The haplotypes most common in Louisiana were shared exclusively with west and east Florida localities in central Florida, indicating a possible source population. The presence of haplotypes private to the Louisiana locality suggests other unsampled localities may also be contributing to the Louisiana settlement. Metrics of genetic diversity across native and invasive localities did not significantly differ. Furthermore, the Louisiana samples had higher genetic diversity than any single location sampled within Florida. Thus, genetic diversity and our haplotype connectivity suggest the Louisiana population is derived from multiple introductions from Florida. Our study highlights how demographic and genetic analyses can be utilized to understand the source and future expansion potential of invasive populations.
1 Introduction
Invasive herpetofauna pose a threat to biodiversity through disruption of native species interactions (Rice et al., 2011; Richter-Boix et al., 2013; ; ), displacement of native species through predation (; Piquet and López-Darias, 2021), and the spread of infectious diseases (Rivera et al., 2019; ). While over 780 introductions of non-native herpetofauna were documented globally as of 2015, only 250 species have established as independent populations (). The disparity between introduced and established invasive species underscores the importance of life history, environmental variation, propagule pressure and genetic diversity on the ability of an invasive population to establish and thrive. Typically, generalist species are more likely to successfully invade and establish compared to specialists that have more specific habitat requirements (Pitt et al., 2005), and environments more like the invasive species native range tend to be conducive to establishment (). Introductions of non-native herpetofauna largely result from anthropomorphic activity, such as the pet trade (Willson et al., 2011), farming (), bio-control efforts (Shanmuganathan et al., 2010), and stowaways in cargo shipments (Richmond et al., 2015). Thus, the numerous factors influencing which species become introduced combined with the complexity of novel ecosystems present a challenge for predicting whether and which non-native amphibian and reptile species will become established.
Invasive populations can arise from novel founder events or expansions of pre-existing invasive ranges, and each of these mechanisms should produce distinct patterns of genetic diversity in the invasive population. Founder events are generally expected to decrease genetic diversity and increase inbreeding depression (), which has been observed in the case of the coqui (Eleutherodactylus coqui) invasion in Hawaii (Peacock et al., 2009) and the American bullfrog (Rana catesbeiana) invasion in China (). Contrary to this expectation, invasive populations founded by multiple introduction events or with high propagule pressure (i.e., many individuals introduced in one or more release events) (Rius and Darling, 2014) can maintain genetic diversity at adaptive loci (Selechnik et al., 2019; ) or even show higher overall genetic diversity relative to native populations (; ). A genetically admixed invasive population has increased adaptive potential and often is capable of faster range expansions (Wagner et al., 2017) and adaptation to novel environmental conditions (; Rius and Darling, 2014).
Mitochondrial DNA (mtDNA) is commonly used to establish phylogenetic relationships between populations and trace the origins of invasive populations (; ; ) because it has a higher mutation rate than nuclear loci, allowing for increased resolution of genetic variation (). Additionally, mtDNA is easy to isolate, has a relatively simple genetic structure (i.e., lacking transposable elements, pseudogenes and introns) and exhibits a straightforward mode of inheritance (). The mitochondrial gene cytochrome b (cyt-b) codes for a central redox catalytic subunit involved in the electron transport chain that is present in nearly all eukaryotic organisms. The cyt-b gene is used to evaluate genetic diversity, trace geographic origins, and distinguish cryptic invasive species across vertebrate groups (; ; ; ). Cyt-b data have been used to investigate the origins of several invasive herpetofauna, including the novel invasion of two subspecies of Boa constrictor in Puerto Rico (Reynolds et al., 2013), the origin and identity of cryptic Xenopus species in Chile () and Italy (), the invasion of the greenhouse frog (Eleutherodactylus planirostris) in Florida () and China (), and the invasive American bullfrog (R. catesbeiana) across the western United States (; ; ).
The Cuban treefrog (Osteopilus septentrionalis) is one of the most successful invasive amphibians in the state of Florida (). Invasive populations negatively influence native Florida treefrogs through predation (Wyatt and Forys, 2004), and by influencing pathogen dynamics (; ). The native O. septentrionalis range spans Cuba, the Cayman Islands and the Bahamas, with invasive populations first establishing in the Florida Keys in the late 1800s to early 1900s (; ). Subsequently, O. septentrionalis spread to peninsular Florida by the mid-1900s via multiple anthropogenic introductions () and has now spread throughout much of central and northern Florida and in some areas of the Florida panhandle (Figure 1; ). Sporadic sightings of O. septentrionalis have been recorded along the southern coastline of Alabama and Mississippi (), but no established populations have been reported. Cyt-b from O. septentrionalis has been used to evaluate the ecological drivers of haplotype and species diversity (Rodríguez et al., 2015). Akin to the previous herpetological invasions mentioned above, the origin of invasive O. septentrionalis in Florida was investigated using cyt-b genetic analysis (). The resulting mitochondrial phylogeny revealed a deep phylogenetic split among haplotypes of O. septentrionalis, with both clades represented throughout the native range in Cuba and the invasive range in Florida. This finding indicated that current populations are descended from at least two source populations ().
Figure 1
In 2017, the first established breeding population of O. septentrionalis in Louisiana, USA, was observed in New Orleans (, ). Although O. septentrionalis had been sporadically documented in Louisiana since the 1990s, no breeding colonies had been observed prior to this event (; , ). In contrast, the discovery of individuals in different age classes, ranging from tadpoles to metamorphs, were found in Audubon Zoo and the adjacent Riverview Park beginning in 2016, leading to further investigation of this invasion event (). Subsequently, a second breeding population was located in St. Charles Parish (). While the novel New Orleans invasion coincided with the arrival of palms that were brought from Lake Placid, Florida for an Audubon Zoo exhibit (), given the occasional observations of O. septentrionalis prior to this event, and that O. septentrionalis are regularly dispersed via anthropogenic activities (Palis et al., 2021), a definitive source for this invasion remains unclear.
Here, we used cyt-b mitochondrial genetics to evaluate genetic variation in the recently established O. septentrionalis population at the Riverview Park and St. Charles Parish sites in Louisiana and assess likely routes of introduction. We sampled and generated cyt-b haplotypes from 31 east Florida individuals and 40 west Florida individuals within central Florida, and 95 individuals in New Orleans, Louisiana, then combined these new data with previously generated sequences from throughout the native and invasive range (; Rodríguez et al., 2015). We conducted phylogenetic and haplotype network analyses and compared metrics of genetic diversity in Louisiana to other native and invasive populations. We used these data to assess whether the Louisiana invasion represents a single introduction event or a more complex history of population establishment.
2 Materials and methods
We collected and processed a total of 166 O. septentrionalis between 2017 and 2022. Of those, 31 were from east Florida, 40 from west Florida, and 95 individuals were from Louisiana populations. East Florida samples were collected from one site in Orange County and west Florida samples were collected from one site in Hillsborough County. Louisianan samples originated from two sampling sites; 85 individuals were from Audubon Riverview Park and the adjacent Audubon Zoo in Orleans Parish, while 10 individuals were from 15 kilometers west in St. Charles Parish (Figure 1).
2.1 Louisiana field sampling
As a means to collect O. septentrionalis, artificial PVC refugia were set at the Riverview Park site in 2019 to capture juvenile metamorphs and adults (). Additional samples were collected by hand via standardized visual surveys. Toe clips from metamorphosed anurans were collected in 99% molecular grade ethanol. Scissors were sterilized after each toe clip by alcohol immersion. O. septentrionalis were euthanized via 20% benzocaine. Samples were placed in 99% molecular grade ethanol and stored at -80 °C prior to DNA extraction.
2.2 Florida field sampling
Juvenile metamorphs and adult frog samples were collected by hand via standardized visual surveys. Tadpoles were collected via standardized dip netting and minnow trapping. Tail clips from tadpoles and toe clips from metamorphosed anurans were collected with flame sterilized scissors and placed in 99% molecular grade ethanol. Per Florida regulations on invasive species, O. septentrionalis were euthanized via MS-222 (tricane mesylate). Tissues were preserved in 99% molecular grade ethanol and stored at -20 °C prior to DNA extraction.
2.3 Laboratory methods
We performed DNA extractions on toe or tail tissue samples using a Qiagen DNeasy Blood and Tissue Kit (Qiagen Corporation, Germantown, Maryland, USA) following manufacturer instructions. Extractions were stored at -20 °C prior to further molecular analysis. PCR was performed on all samples to amplify the cyt-b mitochondrial gene using the forward and reverse primers CBL17 (5’-TAGCCTTYTCATCCGTYGCCCATAT-3’) and CBH18 (5’-GTTGATAATGCAACCCTRACCCGATT-3’) (). Each 25µL reaction consisted of 16.4µL of molecular grade water, 5µL of NEB OneTaq Buffer, 0.5µL of 10mM dNTPs, 0.5µL each of 10uM forward and reverse primers, 0.13µL of NEB One Taq DNA polymerase and 2µL of template DNA or controls. Thermocycling conditions for PCR consisted of an initial denaturation at 94 °C for 5 minutes, followed by 40 cycles of 94 °C for 30 seconds, 46 °C for 30 seconds, and 72 °C for 60 seconds, with a final elongation step of 72 °C for one minute (). PCR products were sent to Eurofins Genomics LLC for Sanger sequencing in the forward and reverse direction. Our newly generated sequences were trimmed of primers and manually inspected and edited to resolve errors and ambiguities by analyzing sanger sequencing graphic outputs. Any sequences that were shorter than 670 bp or had ambiguities, such as stop codons, were excluded from the analyses.
2.4 Phylogenetic analyses
We reconstructed a Bayesian gene tree of O. septentrionalis cyt-b using our sequences (N = 166) and a mix of sequences from native and invasive populations available in GenBank (N = 210; ; Rodríguez et al., 2015). We used Rana septentrionalis (AY083273), Hyla molleri (MK172098), Osteopilus marianae (HQ831741), Osteopilus vastus (HQ831742), Osteopilus dominicensis (HQ831743), and Osteopilus ocellatus (HQ831744) as outgroups. The untruncated sequences of 903 bp were aligned using Clustal Omega (Sievers et al., 2011) with default parameters in Geneious v. 9.0.5 software (). To find a model of evolution that best fit the data and to use in the phylogenetic reconstruction, we used jmodeltest2 v. 2.1.1 () on the CIPRES Science Gateway server (). The best model according to the Akaike information criterion for small sample sizes (AICc) was the Hasegawa-Kishino-Yano (HKY) model with invariant sites and a gamma distribution. We reconstructed a Bayesian gene tree using MrBayes v. 3.2.7a (Ronquist et al., 2012) on the CIPRES Science Gateway server (). The phylogenetic reconstruction consisted of two independent runs of 1.0 × 107 generations and four chains each. We sampled every 500th generation and the first 100,000 generations were discarded as burnin. We used Tracer v. 1.7 (Rambaut et al., 2018) to confirm Markov chain Monte Carlo convergence and adequate sampling of the posterior distribution. We also reconstructed a maximum likelihood gene tree with IQ-TREE with 1000 ultrafast bootstraps and SH-aLRT branch tests on the IQ-TREE webserver (Trifinopoulos et al., 2016) to compare its topology to that of MrBayes. We visualized the Bayesian gene tree and geographic locations where each haplotype has been found using the R package ggtree v3.3.0.900 in RStudio using R v. 4.1.2 (Yu et al., 2017; RStudio Team, 2021).
2.5 Haplotype network
We reconstructed a haplotype network of O. septentrionalis cyt-b utilizing all sequences we generated from our west and east Florida sampling, as well as all additional sequences sampled from Florida O. septentrionalis that were available in GenBank, resulting in the addition of a south Florida locality. Sequences in our alignment were truncated to produce uniform 712 bp long sequences. Cyt-b sequences generated in previous studies that were shorter than 330 bp were excluded from the haplotype network to preserve the quality of our analysis. We used the program PopArt v. 1.7 () and the TCS algorithm () to visualize haplotype connectivity and distribution with a haplotype network analysis.
2.6 Genetic diversity
Untruncated sequence lengths of 903 bp were used to compare genetic diversity among native and invasive locations utilizing our east and west Florida samples, and all O. septentrionalis that were available in GenBank. We separated cyt-b samples into those occurring in Louisiana, central Florida, west Florida, south Florida, west Cuba, central Cuba, east Cuba, and the Caribbean. For each location, we calculated the number of segregating sites (S) and nucleotide diversity (pi) for the preceding locations using pegas v. 1.2 (Paradis, 2010) and ape v. 5.7.1 (Paradis and Schliep, 2019) in RStudio. We generated the distribution of pi with 95% confidence intervals by employing a bootstrap resampling technique with 10,000 iterations using the boot function from boot package v. 1.3.28.1 in R (; ). For reproducibility, we set the random number generation seed to 500 using the set.seed function in R.
3 Results
We generated robust cyt-b sequences from all 166 frogs we sampled in west Florida, east Florida, and Louisiana, recovering 14 total haplotypes. Of these, nine are novel (GenBank Accession numbers: PP320231.1 – PP320239.1) and five were previously published (; Table 1). Novel haplotype se86 was the most common sequence found in our sampling, occurring in 49% of all frogs (Table 1).
Table 1
| Haplotype | Frequency (n) | Localities | ||
|---|---|---|---|---|
| Louisiana | west Florida | east Florida | ||
| se3 | 19 | 18 | 1 | |
| se7 | 1 | 1 | ||
| se13 | 7 | 2 | 5 | |
| se17 | 1 | 1 | ||
| se29 | 8 | 1 | 2 | 5 |
| se86* | 81 | 36 | 29 | 16 |
| se87* | 38 | 37 | 1 | |
| se88* | 4 | 4 | ||
| se89* | 1 | 1 | ||
| se90* | 1 | 1 | ||
| se91* | 2 | 2 | ||
| se92* | 1 | 1 | ||
| se93* | 1 | 1 | ||
| se94* | 1 | 1 | ||
Cytochrome B (cyt-b) mitochondrial haplotypes recovered from our sampling efforts.
Haplotypes with a star (*) are novel to this study.
In Louisiana we recovered seven haplotypes from 95 individuals, including two previously published haplotypes (Table 1) and five novel haplotypes. Among the five novel haplotypes, se89, se90, and se94 were private alleles recovered exclusively from Louisiana, each from a single individual (Table 1). Novel haplotype se87 was the most common haplotype in Louisiana (N = 37 individuals) and it was also found in one east Florida frog. Novel haplotype se86 was the second most common haplotype in Louisiana (N = 36 individuals), and we also found it at high frequency in west and east Florida (Figure 2). All ten individuals collected from the St. Charles Parish population had this haplotype. Haplotypes se3 (N = 17) and se29 (N = 1) were recovered exclusively in Florida in previous studies. None of the haplotypes present in Louisiana occurred in frogs from Cuba or the Caribbean (electronic supplemental material, Supplementary Table S1).
Figure 2
The 31 frogs we sampled from east Florida had seven haplotypes, including four novel haplotypes (se92, se92, se86, se87) and three previously published haplotypes (se13, se17, se29). Only two individuals had the novel private haplotype se91, and one individual had the novel private haplotype se92. Most east Florida individuals had the novel haplotype se86 (N = 16), which was also found in west Florida and Louisiana. One east Florida individual had the novel haplotype se87, which was the most common haplotype in Louisiana. The three haplotypes se13, se17, and se29 were previously recovered from west Florida and west Cuba, south Florida and the Caribbean, and west Florida and south Florida, respectively (electronic supplemental material, Supplementary Table S1). Here, we recovered all three of these haplotypes the first time in east Florida.
In west Florida, we recovered seven haplotypes from 40 individuals: three novel haplotypes (se86, se88, se93) and four previously published haplotypes (se3, se7, se13, se29). A single individual had the private novel haplotype se93 and four individuals had the private novel haplotype se88. The novel haplotype se86 was the most common haplotype in west Florida, similar to its high frequency in Louisiana and east Florida. Se3 was previously recovered only in west Florida, but here we found it in one west Florida frog and in 17 Louisiana frogs. We found se7 in west Florida for the first time, but only in one frog. This haplotype was previously found in a single south Florida individual (electronic supplemental material, Supplementary Table S1). We found two other low-frequency haplotypes (se13 and se29) that had previously been detected in Florida and Cuba (electronic supplemental material, Supplementary Table S1).
Our haplotype network illustrates that almost all frogs sampled in Louisiana had haplotypes se3, se86, or se87 (Figure 2). These sequences are only found in Louisiana, east Florida (se86 and se87), and west Florida (se3 and se86). Only one Louisiana frog had a haplotype that is recovered more broadly throughout Florida (se29), and this sequence was many mutational steps away from all other haplotypes found in Louisiana. Private Louisiana haplotypes se89 and se90 were respectively three and one mutational steps away from se86. Likewise, private Louisiana haplotype se94 is six mutational steps from se87, a haplotype only rarely found outside of Louisiana. Collectively these patterns suggest the evolution of new haplotypes within the invasive Florida and/or Louisiana populations.
Phylogenetic analysis of O. septentrionalis mitochondrial sequences revealed minimal geographic structure (Figure 3). We recovered a well-supported basal split between two sequences only found in the Caribbean and all other sequences. The next split was also well-supported and separated one sequence private to Cuba with everything else. All other sequences were separated into two large clades with moderate support (posterior probability = 0.8917). Both of these clades contained numerous haplotypes found in multiple localities through the native and invasive range. Haplotypes sampled from Louisiana frogs were dispersed throughout the phylogeny.
Figure 3
Among all known cyt-b sequences, Cuba contained the highest number of haplotypes. Frogs sampled in west Cuba had more haplotypes than any other location, followed by east Cuba and central Cuba (Table 2). The Caribbean, Louisiana, and Florida localities each contained similar numbers of haplotypes (Table 2). Native Cuban localities contained similar percentages of private haplotypes, which were high in comparison to the invasive range. Private haplotypes made up 20–30% of the haplotypes recovered in the invasive range localities of east Florida, west Florida, and south Florida whereas 42.86% of Louisiana haplotypes were private (Table 2). Nucleotide diversity in Louisiana was highest among invasive locations in the United States (Figure 4). Nucleotide diversity in Louisiana (0.039; 95% CI = 0.033–0.056) did not overlap with nucleotide diversity of central Cuba (0.017; 95% CI = 0.0092–0.021) or east Cuba (0.018; 95% CI = 0.010–0.023; Table 2; Figure 4). Segregating sites were higher in west Cuba and the Caribbean than within invasive localities in the United States (Table 2).
Table 2
| Location (No. frogs sampled) | No. haplotypes | No. private haplotypes | Percentage of private haplotypes | No. segregating sites | Nucleotide diversity (standard error) |
|---|---|---|---|---|---|
| Louisiana (95) | 7 | 3 | 42.86% | 60 | 0.039 (± 0.02214) |
| east Florida (31) | 7 | 2 | 28.57% | 54 | 0.024 (± 0.01378) |
| west Florida (56) | 8 | 2 | 25.00% | 67 | 0.032 (± 0.01761) |
| south Florida (24) | 5 | 1 | 20.00% | 62 | 0.032 (± 0.01975) |
| west Cuba (70) | 36 | 29 | 80.56% | 86 | 0.040 (± 0.01949) |
| central Cuba (29) | 20 | 15 | 75.00% | 54 | 0.017 (± 0.00900) |
| east Cuba (55) | 28 | 20 | 71.43% | 67 | 0.018 (± 0.00927) |
| Caribbean (16) | 7 | 5 | 71.43% | 99 | 0.048 (± 0.08660) |
Summary of allelic diversity, private alleles, segregating sites, and nucleotide diversity across eight locations.
Figure 4
4 Discussion
Our study is the first to investigate the origins of the establishment of O. septentrionalis in Louisiana, the first documented establishment of O. septentrionalis in the contiguous U.S. outside of Florida (). We aimed to determine whether this invasion was the result of one or more founding events from remote locations. Our results indicate both Louisiana O. septentrionalis populations are the result of multiple introduction events, likely from central Florida localities. Four of seven haplotypes found in Louisiana individuals were present only in Florida and Louisiana, while the remaining three were rare and only found in Louisiana. This pattern excludes Caribbean locations as direct sources of the novel population and strongly points toward Florida origins. One haplotype lending the most support to our introduction hypothesis is se86, which is novel to our study (Table 1). It was the most abundant haplotype we recovered overall and was at high frequency in our Florida and Louisiana sample sites, consistent with the hypothesis that central Florida is a source of the Louisiana O. septentrionalis.
Though we did not sample any localities in the Florida panhandle, the likelihood that frogs spread naturally from Florida to Louisiana in a range expansion is minimal. The high nucleotide diversity of haplotypes in Louisiana is inconsistent with the expectation for a range expansion. Range expansions of invasive species typically result in decreased genetic diversity as the leading edge expands (Wilson et al., 2005; ). Instead, we see higher nucleotide diversity in Louisiana compared to all locations except for west Cuba and the Caribbean (Figure 4), which is where the species originated (; ). Notably, the nucleotide diversity in Louisiana is nearly equivalent to values found in west Cuba. Furthermore, there is no evidence for a leading edge, as established populations of O. septentrionalis between the Florida panhandle and the Louisiana population have not yet been documented (). Future work that samples O. septentrionalis from the Florida panhandle and the Gulf Coast of Alabama and Mississippi, and analyzes their genetic composition could enable assessment of potential new population establishments and identification of their source populations.
High genetic diversity metrics in the Louisiana population suggests a high propagule pressure invasion and/or multiple introduction events (; Simberloff, 2009; Wang et al., 2019). Given that human transport is the most plausible route of invasion, multiple founding events is the most likely scenario explaining the patterns we recovered in Louisiana. For example, repeated transport of one or a few individuals traveling on ornamental plants purchased from locations within Florida and planted in New Orleans is a plausible source of the novel invasion. Anecdotal accounts of O. septentrionalis being found within many purchased plants from Florida (), the scale and economic importance of the nursery industry in Florida (), and low efficacy of state-level biosecurity measures in preventing the spread of established invasive species (Paini et al., 2010) lend support to this hypothesis. Therefore, it is plausible that the known planting of ornamental palms from Lake Placid, Florida in the Audubon Zoo contributed to this invasion (). Sampling in or around the nursery in Lake Placid—which lies farther south than our sampling sites—to verify the presence of haplotypes shared with Louisiana could strengthen support for this hypothesis. Additionally, increased sampling at or near ornamental horticultural nurseries and retailers in top economic contributing regions, such as central and south Florida (), could give greater insight into the impacts of this industry on the spread of O. septentrionalis.
Human traffic near the invasion site has the potential to facilitate introductions from distant source populations. New Orleans is a large city that is visited by nearly 20 million tourists annually, with many visitors arriving via vehicular transportation or large cruise ships that depart from docks approximately 5 miles from the invasion site (). Florida is a frequent point of origin, ranking as the 3rd most common origin of interstate travel entering Louisiana in 2023 (). As vehicular transport is the primary mode of invasion of O. septentrionalis (), it is plausible that multiple instances of individual frogs hitch-hiking within vehicles departing from central Florida have contributed to the New Orleans invasion. Additionally, Interstate 10 (I-10), a major thoroughfare linking Florida to Louisiana, provides direct access to New Orleans and is heavily utilized by Floridians during hurricane evacuations (). While hurricane evacuees likely are minor contributors in comparison to tourists, increasing hurricane severity and resulting evacuation efforts may increase the frequency of these individual translocation events. Additionally, cruise ships departing and arriving at the Port of New Orleans frequent many Caribbean islands and various Floridian ports. Caribbean sampling of O. septentrionalis is minimal, leaving unclear resolution on what haplotypes are present in this region. Increased Caribbean sampling may reveal haplotypes from this study found exclusively in Louisiana may also occur in the Caribbean and may be derived from there. Overall increased sampling efforts throughout Florida, as well as the Caribbean could provide increased clarity on the contribution of translocation events to the founding of the Louisianan population.
previously reported two divergent clades of O. septentrionalis to be phylogeographically grouped based on relation to the Guanahacabibes peninsula in remote west Cuba. They additionally reported Bahamian haplotypes form a basal monophyletic group. With the additional samples of Rodríguez et al. (2015) as well as our novel sequences, our phylogenetic analysis groups together west Cuban haplotypes within and outside of the Guanahacabibes peninsula together, indicating this peninsular boundary does not impede gene flow to other localities in Cuba. Our study does not expand on the finding of Bahamian samples forming a monophyly, as we did not include additional Bahamian sampling outside of providing another reason for additional sampling in the Caribbean and specifically the Bahamas.
We expect that O. septentrionalis will continue to negatively affect native species through resource use and influences on local pathogen dynamics, but changes in the genetic composition of O. septentrionalis populations are a potential cause for concern. Anecdotal sightings of O. septentrionalis outside their established invasive territory are plentiful and range from Georgia to as far north as Vermont along the east coast of the United States (). These sightings are likely a result of anthropogenic transport either via horticulture or other items transported from their invasive range to wider regions of the U.S. and are not yet representative of established breeding colonies. Cold temperatures play a large role in limiting the northern distribution of O. septentrionalis, which have been observed dead at temperatures below 0 °C (), but shifting global temperatures threaten this climatic barrier to expansion. With fewer freeze events occurring in regions north of the current range of O. septentrionalis, these populations may persist and establish farther north than previously expected (Osland et al., 2021).
Admixed populations at range limits may also benefit from an increased evolutionary potential resulting from heterosis, thereby increasing the extent of dispersal and probability of permanent establishment (; Wagner et al., 2017). Novel genotypes resulting from heterosis (Verhoeven et al., 2011) may allow for adaptation to colder temperatures. Evidence shows O. septentrionalis from the northern latitudinal limits of their Florida range had critical thermal minima trending lower than their southern counterparts and that in general O. septentrionalis has a higher tolerance for cold weather than previously believed (Simpson, 2013). Though cold tolerance tends to result from phenotypic plasticity (; ), it is possible that exposure to freeze events and the existence of novel genotypes within northern populations result in adapted cold tolerance. Though we observed patterns suggesting local adaptation of the cyt-b mitochondrial gene in Florida and Louisiana, more work could help to assess adaptation in functionally relevant genes. Future experimental studies on mitochondrial gene transcription associated with thermal regulation across invasive populations could provide deeper insight into thermal adaptation ().
Disparities in size between northern and southern populations of O. septentrionalis have also been recorded (). Studies do not yet indicate if this shift to smaller sizes in colder latitudes is genetically based or if these changes result from environmental effects such as shifts in food availability or metabolic constraints resulting from cold temperatures. It is possible that smaller body size and increasing urban development may aid in their spread, as O. septentrionalis can exploit urban refugia to escape or minimize exposure to thermal temperature limits (; ).
Invasive species present a unique opportunity to study evolution in action by assessing novel selection pressures in situ. These systems can serve as natural experiments, with introduced populations functioning as treatment groups and native populations acting as controls. Invasive O. septentrionalis have the potential to serve as a powerful model for studying rapid adaptation to novel environmental conditions, especially in the recently established Louisiana population. While mitochondrial genes markers have historically served as a valuable tool in determining the origin of invasive species, next-generation sequencing (NGS) techniques offer better resolution to fine scale population structure and more robust findings on this recent population divergence (North et al., 2021). Whole mitochondrial genome analysis can expand on findings derived from single genes for more comprehensive insights as demonstrated with Rhinella marina, the invasive cane toad (). Like many amphibian species, O. septrentrionalis has not yet had its nuclear or mitochondrial genome sequenced, due in part to the difficulty of assemblage (; O’Connell et al., 2024) inhibiting the depth of evolutionary questioning utilizing this species.
Our study concludes that the newly established Louisiana population of O. septentrionalis is derived from multiple introductions, primarily from central Florida. These introductions likely occurred via anthropogenic means (), though we cannot conclude with certainty whether the transport of flora from Florida to Louisiana contributed to this establishment. To fully resolve the origin of Louisiana’s population, future work could prioritize increased sampling across the Florida panhandle, Caribbean, and key Floridian localities. By identifying the source populations and the degree of admixture within the Louisiana population, we gain insight into where management efforts could be prioritized and understanding of the genetic variation available for selection to act upon. In the case of O. septentrionalis, this knowledge can be used not only for reconstructing invasion history, but could also anticipate the species’ capacity to adapt to colder climates, urban environments, and other novel pressures.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The animal study was approved by Institutional Animal Care and Use Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
EB: Investigation, Writing – review & editing, Visualization, Formal analysis, Project administration, Writing – original draft, Data curation. KP: Data curation, Investigation, Writing – review & editing. KM: Investigation, Writing – review & editing, Formal analysis, Visualization, Data curation, Project administration. MA: Data curation, Project administration, Investigation, Funding acquisition, Writing – review & editing, Writing – original draft, Conceptualization. BG: Methodology, Investigation, Resources, Data curation, Project administration, Funding acquisition, Conceptualization, Writing – review & editing. JW: Conceptualization, Methodology, Resources, Data curation, Investigation, Writing – review & editing, Project administration, Funding acquisition. RM: Investigation, Data curation, Writing – review & editing, Resources. AS: Writing – review & editing, Funding acquisition, Investigation, Supervision, Resources, Project administration, Methodology.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. USGS, NSF and FWC contributed funds supporting research. Article processing charges were provided in part by the UCF College of Graduate Studies Open Access Publishing Fund.
Acknowledgments
Thanks to the Conservation, Restoration and Communication REU, and students from the Savage lab including Maddie Thurber, Siena Krullis, Ivan Santana, Patricia Hernandez, and Savannah Freeman who assisted in catching or processing samples. Thanks to Maryam Ghoojaei for her guidance on implementing bootstrapping analyses in RStudio. Animals in Florida were captured under Florida Fish and Wildlife Conservation Commission Scientific Collecting Permits and handled in accordance with the approved IACUC protocols. We thank Lindy J. Muse, Nicole D. Jennings, Brittany R. Maldonado, Raymond P. Kidder, Melanie Litton, CJ Hillard, Katie Everett, and many volunteers and a few private landowners for Louisiana tissue collection. We captured animals under annual Louisiana Department of Wildlife and Fisheries Scientific Collecting Permits and handled all animals in accordance with approved IACUC protocols (USGS WARC FY2008–1). Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. This is contribution number 959 of the U.S. Geological Survey Amphibian Research and Monitoring Initiative (ARMI).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/famrs.2025.1646714/full#supplementary-material
References
1
AdamsA. (2024). “Louisiana visitor profile calendar year 2023,” in MMGY travel intelligence. Mclean, VA: MMGY Travel Intelligence. Available online at: https://www.explorelouisiana.com/sites/default/files/2024-07/2023%20Louisiana%20Visitor%20Profile%20Report.pdf (Accessed March 30, 2025).
2
AllendorfF. W.LundquistL. L. (2003). Introduction: population biology, evolution, and control of invasive species. Conserv. Biol.17, 24–30. doi: 10.1046/j.1523-1739.2003.02365.x
3
AtkinsonM. S.SavageA. E. (2023). Invasive amphibians alter host-pathogen interactions with primarily negative outcomes for native species. Biol. Conserv.286, 110310. doi: 10.1016/j.biocon.2023.110310
4
AviseJ. C.ArnoldJ.BallR. M.BerminghamE.LambT.NeigelJ. E.et al. (1987). Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics. Annu. Rev. Ecol. Syst.18, 489–522. doi: 10.1146/annurev.es.18.110187.002421
5
BaiC.KeZ.ConsuegraS.LiuX.LiY. (2012). The role of founder effects on the genetic structure of the invasive bullfrog (Lithobates catesbeianaus) in China. Biol. Invasions.14, 1785–1796. doi: 10.1007/s10530-012-0189-x
6
BomfordM.KrausF.BarryS. C.LawrenceE. (2009). Predicting establishment success for alien reptiles and amphibians: a role for climate matching. Biol. Invasions.11, 713–724. doi: 10.1007/s10530-008-9285-3
7
CantyA. (2021). boot: Bootstrap R (S-Plus) functions. R package version 1. 3–28.
8
ChatfieldM. W. H.VanceM. (2014). Geographic distribution. Osteopilus septentrionalis (Cuban Treefrog). Herpetol. Rev.45 (2), 278. Available online at: https://ssarherps.org/herpetological-review-pdfs/.
9
CheungK.AmosT. G.ShineR.DeVoreJ. L.DucatezS.EdwardsR. J.et al. (2024). Whole-mitogenome analysis unveils previously undescribed genetic diversity in cane toads across their invasion trajectory. Ecol. Evol.14, e11115. doi: 10.1002/ece3.11115
10
ClementM.PosadaD.CrandallK. A. (2000). TCS: a computer program to estimate gene genealogies. Mol. Ecol.9, 1657–1659. doi: 10.1046/j.1365-294x.2000.01020.x
11
ColangeloP.AbiadhA.AloiseG.AmoriG.CapizziD.VasaE.et al. (2015). Mitochondrial phylogeography of the black rat supports a single invasion of the western Mediterranean basin. Biol. Invasions.17, 1859–1868. doi: 10.1007/s10530-015-0842-2
12
CulbertsonK. A.HerrmannN. C. (2019). Asymmetric interference competition and niche partitioning between native and invasive Anolis lizards. Oecologia.190, 811–820. doi: 10.1007/s00442-019-04466-1
13
DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: more models, new heuristics and high-performance computing. Nat. Methods9, 772. doi: 10.1038/nmeth.2109
14
DavisM. C.NovakS. J.HampikianG. (2011). Mitochondrial DNA analysis of an immigrant Basque population: loss of diversity due to founder effects. Am. J. Phys. Anthropol.144, 516–525. doi: 10.1002/ajpa.21432
15
DavisonA. C.HinkleyD. V. (1997). Bootstrap methods and their Application (Cambridge, UK: Cambridge University Press).
16
DoveC. J.SnowR. W.RochfordM. R.MazzottiF. J. (2011). Birds consumed by the invasive Burmese python (Python molurus bivittatus) in Everglades National Park, Florida, USA. Wilson. J. Ornithol.123, 126–131. doi: 10.1676/10-092.1
17
DuellmanW.SchwartzA. (1958). Amphibians and reptiles of southern Florida. Bull. Fla. Mus. Nat. Hist.3, 181–324. doi: 10.58782/flmnh.gtne5332
18
ForsmanA. (2013). Effects of genotypic and phenotypic variation on establishment are important for conservation, invasion, and infection biology. Proc. Natl. Acad. Sci.111, 302–307. doi: 10.1073/pnas.1317745111
19
FrankhamR. (1998). Inbreeding and extinction: island populations. Conserv. Biol.12, 665–675. doi: 10.1046/j.1523-1739.1998.96456.x
20
GaltN.AtkinsonM.GloriosoB.WaddleH.LittonM.SavageA. (2021). Widespread ranavirus and perkinsea infections in Cuban treefrogs (Osteopilus septentrionalis) invading New Orleans, USA. Herpetol. Conserv. Biol.16, 17–29. doi: 10.5061/DRYAD.7WM37PVS8
21
GhorbanzadehM.BurnsS.RugminiammaL. V. N.Erman OzguvenE.HuangW. (2021). Spatiotemporal analysis of highway traffic patterns in Hurricane Irma evacuation. Transp. Res. Rec.2675, 321–334. doi: 10.1177/03611981211001870
22
GloriosoB. M.MacedoK.MaldonadoB. R.HillardC. J.MorencI. N.GrimesE. S. (2018a). Geographic distribution: Osteopilus septentrionalis (Cuban treefrog). Herpetol. Rev.49 (4), 709. Available online at: https://ssarherps.org/herpetological-review-pdfs/.
23
GloriosoB. M.SteeceA.LemannZ. K.LazareR.BeckJ. W. (2016). Geographic distribution: Osteopilus septentrionalis (Cuban Treefrog). Herpetol. Rev.47 (2), 249.
24
GloriosoB. M.VanbergenP.RoyJ.WalterM.LeonpacherL.FreistakM. (2018b). Geographic distribution: Osteopilus septentrionalis (Cuban Treefrog). Herpetol. Rev.49 (1), 70–71. Available online at: https://ssarherps.org/herpetological-review-pdfs/.
25
GloriosoB. M.WaddleJ. H. (2014). A review of pipe and bamboo artificial refugia as sampling tools in anuran studies. Herpetol. Conserv. Biol.9, 609–625. Available online at: https://www.herpconbio.org/Volume_9/Issue_3/Glorioso_Waddle_2014.pdf.
26
GloriosoB. M.WaddleJ. H.MuseL. J.JenningsN. D.LittonM.HamiltonJ.et al. (2018c). Establishment of the exotic invasive Cuban treefrog (Osteopilus septentrionalis) in Louisiana. Biol. Invasions.20, 2707–2713. doi: 10.1007/s10530-018-1732-1
27
GrayM. W.CedergrenR.AbelY.SankoffD. (1989). On the evolutionary origin of the plant mitochondrion and its genome. Proc. Natl. Acad. Sci.86, 2267–2271. doi: 10.1073/pnas.86.7.2267
28
HaggertyC. J. E.CrismanT. L. (2015). Pulse disturbance impacts from a rare freeze event in Tampa, Florida on the exotic invasive Cuban treefrog, Osteopilus septentrionalis, and native treefrogs. Biol. Invasions.17, 2103–2111. doi: 10.1007/s10530-015-0863-x
29
HamnerR.FreshwaterD. W.WhitfieldP. (2007). Mitochondrial cytochrome b analysis reveals two invasive lionfish species with strong founder effects in the western Atlantic. J. Fish Biol.71, 214–222. doi: 10.1111/j.1095-8649.2007.01575.x
30
HeinickeM. P.DiazL. M.HedgesS. B. (2011). Origin of invasive Florida frogs traced to Cuba. Biol. Lett.7, 407–410. doi: 10.1098/rsbl.2010.1131
31
HillS. A.BeardK. H.SiersS. R.ShielsA. B. (2019). Invasive coqui frogs are associated with differences in mongoose and rat abundances and diets in Hawaii. Biol. Invasions.21, 2177–2190. doi: 10.1007/s10530-019-01965-3
32
HodgesA. W.KhachatryanH.RahmaniM. (2016). “Economic contributions of the environmental horticulture industry in Florida in 2015,” in Electronic data information source, university of Florida institute of food and agricultural sciences. Gainesville, FL: University of Florida-Institute of Food and Agricultural Sciences Food and Resource Economics Department. Available online at: https://web.archive.org/web/20180426111529id_/http://fred.ifas.ufl.edu/pdf/EconContEnvirHortIndFL2015-11-15-16.pdf.
33
HongY.HeY.LinZ.DuY.ChenS.HanL.et al. (2022). Complex origins indicate a potential bridgehead introduction of an emerging amphibian invader (Eleutherodactylus planirostris) in China. NeoBiota.77, 23–37. doi: 10.3897/neobiota.77.83205
34
HongY. H.YuanY.-N.LiK.StoreyK. B.ZhangJ.-Y.ZhangS.-S.et al. (2024). Differential mitochondrial genome expression of four Hylid frog species under low-temperature stress and its relationship with Amphibian temperature adaptation. Int. J. Mol. Sci.25, 5967. doi: 10.3390/ijms25115967
35
JacksonH.StrubbeD.TollingtonS.Prys-JonesR.MatthysenE.GroombridgeJ. J. (2015). Ancestral origins and invasion pathways in a globally invasive bird correlate with climate and influences from bird trade. Mol. Ecol.24, 4269–4285. doi: 10.1111/mec.13307
36
JohnsonS. A. (2023). The Cuban treefrog (Osteopilus septentrionalis) in Florida: WEC218/UW259, 2/2023. EDIS.2023, 1. doi: 10.32473/edis-UW259-2023
37
KamathP. L.SepulvedaA. J.LayheeM. (2016). Genetic reconstruction of a bullfrog invasion to elucidate vectors of introduction and secondary spread. Ecol. Evol.6, 5221–5233. doi: 10.1002/ece3.2278
38
KearseM.MoirR.WilsonA.Stones-HavasS.CheungM.SturrockS.et al (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics28, 1647–1649. doi: 10.1093/bioinformatics/bts199
39
KhachatryanH.KnuthM.HodgesA.HallC. (2022). Florida nursery and landscape industry economic contributions report: FE1114/FE1114, 02/2022. EDIS.2022, 1. doi: 10.32473/edis-fe1114-2022
40
KochziusM.SöllerR.KhalafM. A.BlohmD. (2003). Molecular phylogeny of the lionfish genera Dendrochirus and Pterois (Scorpaenidae, Pteroinae) based on mitochondrial DNA sequences. Mol. Phylogenet. Evol.28, 396–403. doi: 10.1016/S1055-7903(02)00444-X
41
KolbeJ. J.GlorR. E.Rodríguez SchettinoL.LaraA. C.LarsonA.LososJ. B. (2004). Genetic variation increases during biological invasion by a Cuban lizard. Nat.431, 177–181. doi: 10.1038/nature02807
42
KolbeJ. J.LarsonA.LososJ. B.De QueirozK. (2008). Admixture determines genetic diversity and population differentiation in the biological invasion of a lizard species. Biol. Lett.4, 434–437. doi: 10.1098/rsbl.2008.0205
43
KoschT. A.Torres-SánchezM.LiedtkeH. C.SummersK.YunM. H.CrawfordA. J.et al. (2024). The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation. BMC Genomics25, 1025. doi: 10.1186/s12864-024-10899-7
44
KrausF. (2015). Impacts from invasive reptiles and amphibians. Annu. Rev. Ecol. Evol. Syst.46, 75–97. doi: 10.1146/annurev-ecolsys-112414-054450
45
KrehenwinkelH.TautzD. (2013). Northern range expansion of European populations of the wasp spider Argiope bruennichi is associated with global warming–correlated genetic admixture and population-specific temperature adaptations. Mol. Ecol.22, 2232–2248. doi: 10.1111/mec.12223
46
LaFondJ.MartinK. R.DahnH.RichmondJ. Q.MurphyR. W.RollinsonN.et al. (2022). Invasive bullfrogs maintain MHC polymorphism including alleles associated with chytrid fungal infection. Integr. Comp. Biol.62, 262–274. doi: 10.1093/icb/icac044
47
LauferG.CanaveroA.NúñezD.ManeyroR. (2008). Bullfrog (Lithobates catesbeianus) invasion in Uruguay. Biol. Invasions.10, 1183–1189. doi: 10.1007/s10530-007-9178-x
48
LeighJ. W.BryantD. (2015). POPART: full-feature software for haplotype network construction. Methods Ecol. Evol.6, 1110–1116. doi: 10.1111/2041-210X.12410
49
LilloF.DufresnesC.FaraoneF.Lo ValvoM.StöckM. (2013). Identification and potential origin of invasive clawed frogs Xenopus (Anura: Pipidae) in Sicily based on mitochondrial and nuclear DNA. Ital. J. Zool.80, 566–573. doi: 10.1080/11250003.2013.847502
50
LobosG.MendezM. A.CattanP.JaksicF. (2014). Low genetic diversity of the successful invasive African clawed frog Xenopus laevis (Pipidae) in Chile. Stud. Neotrop. Fauna Environ.49, 50–60. doi: 10.1080/01650521.2014.912865
51
McCannS.GreenleesM. J.NewellD.ShineR. (2014). Rapid acclimation to cold allows the cane toad to invade montane areas within its Australian range. Funct. Ecol.28, 1166–1174. doi: 10.1111/1365-2435.12255
52
McGarrityM. E.JohnsonS. A. (2009). Geographic trend in sexual size dimorphism and body size of Osteopilus septentrionalis (Cuban treefrog): implications for invasion of the southeastern United States. Biol. Invasions.11, 1411–1420. doi: 10.1007/s10530-008-9349-4
53
MeshakaW. E. (2001). The Cuban treefrog in Florida: life history of a successful colonizing species (Gainesville, FL: University Press of Florida).
54
MillerM. A.PfeifferW.SchwartzT. (2010). Creating the CIPRES Science Gateway for inference of large phylogenetic trees., in 2010 gateway computing environments workshop (GCE). IEEE., 1–8. doi: 10.1109/GCE.2010.5676129
55
MittanC. S.ZamudioK. R. (2019). Rapid adaptation to cold in the invasive cane toad Rhinella marina. Conserv. Physiol.7, coy075. doi: 10.1093/conphys/coy075
56
MorningstarC. R.DanielW. M.SommaL. A. (2024). Cuban Treefrog (Osteopilus septentrionalis) - Species Profile (Gainesville, FL: U.S. Geological Survey). Available online at: https://nas.er.usgs.gov/queries/factsheet.aspx?SpeciesID=57 ((Accessed April 23, 2025)).
57
NorthH. L.McGaughranA.JigginsC. D. (2021). Insights into invasive species from whole-genome resequencing. Mol. Ecol.30, 6289–6308. doi: 10.1111/mec.15999
58
O’ConnellL. A.RodríguezA.KoschT. A.KwonT.BolsoniL.LourençoH.et al. (2024). ““Genomics: using genomics approaches in amphibian conservation”,” in Amphibian conservation action plan. Eds. WrenS.BorzéeA.R. Marcec-GreavesAnguloA. (Gland, Switzerland: IUCN SSC Amphibian Specialist Group), 310–334. Available online at: https://www.researchgate.net/profile/Simon-Clulow/publication/382486709_Amphibian_conservation_action_plan_a_status_review_and_roadmap_for_global_amphibian_conservation/links/66d90fefb1606e24c2e1a991/Amphibian-conservation-action-plan-a-status-review-and-roadmap-for-global-amphibian-conservation.pdf#page=314.
59
OslandM. J.StevensP. W.LamontM. M.BruscaR. C.HartK. M.WaddleJ. H.et al. (2021). Tropicalization of temperate ecosystems in North America: The northward range expansion of tropical organisms in response to warming winter temperatures. Glob. Change Biol.27, 3009–3034. doi: 10.1111/gcb.15563
60
PainiD. R.WornerS. P.CookD. C.De BarroP. J.ThomasM. B. (2010). Threat of invasive pests from within national borders. Nat. Commun.1, 115. doi: 10.1038/ncomms1118
61
PalisJ. G.JohnsonS. A.MarksJ. R.BoehlerM. B. (2021). Cuban treefrog (Osteopilus septentrionalis) in illinois. Bull. Chic. Herpetol. Soc56 (9), 153–155.
62
ParadisE. (2010). pegas: an R package for population genetics with an integrated–modular approach. Bioinformatics.26, 419–420. doi: 10.1093/bioinformatics/btp696
63
ParadisE.SchliepK. (2019). ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics.35, 526–528. doi: 10.1093/bioinformatics/bty633
64
PeacockM. M.BeardK. H.O’NEILLE. M.KirchoffV. S.PetersM. B. (2009). Strong founder effects and low genetic diversity in introduced populations of Coqui frogs. Mol. Ecol.18, 3603–3615. doi: 10.1111/j.1365-294X.2009.04308.x
65
PiquetJ. C.López-DariasM. (2021). Invasive snake causes massive reduction of all endemic herpetofauna on Gran Canaria. Proc. R. Soc B: Biol. Sci.288, 20211939. doi: 10.1098/rspb.2021.1939
66
PittW.ViceD.PitzlerM. (2005). Challenges of invasive reptiles and amphibians. Wildlife Damage Manage. Conferences - Proc.84, 112– 119. Available online at: https://digitalcommons.unl.edu/icwdm_wdmconfproc/84 (Accessed July 17, 2024).
67
RambautA.DrummondA. J.XieD.BaeleG.SuchardM. A. (2018). Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol.67, 901–904. doi: 10.1093/sysbio/syy032
68
ReynoldsR. G.Puente-RolónA. R.ReedR. N.RevellL. J. (2013). Genetic analysis of a novel invasion of Puerto Rico by an exotic constricting snake. Biol. Invasions.15, 953–959. doi: 10.1007/s10530-012-0354-2
69
RiceK. G.WaddleJ. H.MillerM. W.CrockettM. E.MazzottiF. J.PercivalH. F. (2011). Recovery of native treefrogs after removal of nonindigenous Cuban Treefrogs, Osteopilus septentrionalis. Herpetologica.67, 105–117. doi: 10.1655/HERPETOLOGICA-D-10-00020.1
70
RichmondJ. Q.WoodD. A.StanfordJ. W.FisherR. N. (2015). Testing for multiple invasion routes and source populations for the invasive brown treesnake (Boiga irregularis) on Guam: implications for pest management. Biol. Invasions.17, 337–349. doi: 10.1007/s10530-014-0733-y
71
Richter-BoixA.GarrigaN.MontoriA.FranchM.San SebastiánO.VilleroD.et al. (2013). Effects of the non-native amphibian species Discoglossus pictus on the recipient amphibian community: niche overlap, competition and community organization. Biol. Invasions.15, 799–815. doi: 10.1007/s10530-012-0328-4
72
RiusM.DarlingJ. A. (2014). How important is intraspecific genetic admixture to the success of colonising populations? Trends Ecol. Evol.29, 233–242. doi: 10.1016/j.tree.2014.02.003
73
RiveraB.CookK.AndrewsK.AtkinsonM. S.SavageA. E. (2019). Pathogen dynamics in an invasive frog compared to native species. EcoHealth.16, 222–234. doi: 10.1007/s10393-019-01432-4
74
RodríguezA.BörnerM.PabijanM.GeharaM.HaddadC. F. B.VencesM. (2015). Genetic divergence in tropical anurans: deeper phylogeographic structure in forest specialists and in topographically complex regions. Evol. Ecol.29, 765–785. doi: 10.1007/s10682-015-9774-7
75
RonquistF.TeslenkoM.van der MarkP.AyresD. L.DarlingA.HöhnaS.et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol.61, 539–542. doi: 10.1093/sysbio/sys029
76
RStudio Team (2021). RStudio: integrated development environment for R (Boston, MA: RStudio, PBC). Available online at: http://www.rstudio.com/.
77
SelechnikD.RichardsonM. F.ShineR.DeVoreJ. L.DucatezS.RollinsL. A. (2019). Increased adaptive variation despite reduced overall genetic diversity in a rapidly adapting invader. Front. Genet.10. doi: 10.3389/fgene.2019.01221
78
ShanmuganathanT.PallisterJ.DoodyS.McCallumH.RobinsonT.SheppardA.et al. (2010). Biological control of the cane toad in Australia: a review. Anim. Conserv.13, 16–23. doi: 10.1111/j.1469-1795.2009.00319.x
79
SieversF.WilmA.DineenD.GibsonT. J.KarplusK.LiW.et al. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol.7, 539. doi: 10.1038/msb.2011.75
80
SimberloffD. (2009). The role of propagule pressure in Biol. Invasions. Annu. Rev. Ecol. Evol. Syst.40, 81–102. doi: 10.1146/annurev.ecolsys.110308.120304
81
SimpsonS. E. (2013). Assessing critical thermal minima to determine the thermal limits of the invasive Cuban treefrog (Osteopilus septentrionalis). [master’s thesis]. University of Florida, Gainesville (FL.
82
TrifinopoulosJ.NguyenL.-T.von HaeselerA.MinhB. Q. (2016). W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res.44, W232–W235. doi: 10.1093/nar/gkw256
83
VerhoevenK. J.MacelM.WolfeL. M.BiereA. (2011). Population admixture, Biol. Invasions and the balance between local adaptation and inbreeding depression. Proc. R. Soc B: Biol. Sci.278, 2–8. doi: 10.1098/rspb.2010.1272
84
WagnerN. K.OchockiB. M.CrawfordK. M.CompagnoniA.MillerT. E. (2017). Genetic mixture of multiple source populations accelerates invasive range expansion. J. Anim. Ecol.86, 21–34. doi: 10.1111/1365-2656.12567
85
WangS.LiuC.WuJ.XuC.ZhangJ.BaiC.et al. (2019). Propagule pressure and hunting pressure jointly determine genetic evolution in insular populations of a global frog invader. Sci. Rep.9, 448. doi: 10.1038/s41598-018-37007-6
86
WillsonJ. D.DorcasM. E.SnowR. W. (2011). Identifying plausible scenarios for the establishment of invasive Burmese pythons (Python molurus) in Southern Florida. Biol. Invasions.13, 1493–1504. doi: 10.1007/s10530-010-9908-3
87
WilsonG. A.NishiJ. S.ElkinB. T.StrobeckC. (2005). Effects of a recent founding event and intrinsic population dynamics on genetic diversity in an ungulate population. Conserv. Genet.6, 905–916. doi: 10.1007/s10592-005-9077-6
88
WyattJ. L.ForysE. A. (2004). Conservation implications of predation by Cuban treefrogs (Osteopilus septentrionalis) on native hylids in Florida. Southeast. Nat.3, 695–700. doi: 10.1656/1528-7092(2004)003[0695:CIOPBC]2.0.CO;2
89
YuG.SmithD. K.ZhuH.GuanY.LamT. T.-Y. (2017). ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol.8, 28–36. doi: 10.1111/2041-210X.12628
Summary
Keywords
invasive species, genetics, biological invasion, cytochrome b, Osteopilus septentrionalis
Citation
Brosnan EB, Paniagua Torres KA, Martin KR, Atkinson MS, Glorioso BM, Waddle JH, Mendyk RW and Savage AE (2025) Tracing invasion routes of Cuban treefrogs into Louisiana using mitochondrial DNA. Front. Amphib. Reptile Sci. 3:1646714. doi: 10.3389/famrs.2025.1646714
Received
13 June 2025
Accepted
24 September 2025
Published
17 October 2025
Volume
3 - 2025
Edited by
Amanda Kissel, United States Geological Survey, United States
Reviewed by
Steve A. Johnson, University of Florida, United States
Andrew Parks, US Fish and Wildlife Service, United States
Updates
Copyright
© 2025 Brosnan, Paniagua Torres, Martin, Atkinson, Glorioso, Waddle, Mendyk and Savage.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Erin B. Brosnan, erin.brosnan@ucf.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.