Introduction

The international wildlife trade facilitates introductions of invasive species and pathogens that threaten native communities. Bullfrogs are native to eastern North America1,2 and are farmed on a large scale in Asia, Central America, and South America to supply the international frog leg trade2,3. Bullfrog farming emerged as an alternative to overharvesting native amphibian species4. However, this practice continues to contribute to the current global amphibian crisis by facilitating biological invasions5,6,7,8,9,10,11. Escapes of bullfrogs from farms have led to the establishment of invasive bullfrog populations6,12,13,14 that negatively influence the local anurofauna by interfering with acoustic communication and jeopardizing reproduction of native anurans15,16, preying on native amphibian species17,18 and competing for resources, thus reducing the fitness of native amphibian populations19,20.

Bullfrogs also play an important role in the dynamics of amphibian chytridiomycosis, a disease caused by the fungus Batrachochytrium dendrobatidis (Bd) that has been linked to amphibian declines worldwide21. Bullfrogs are highly tolerant hosts22,23, meaning that they are able to withstand high Bd infection loads without developing chytridiomycosis (but see exceptions in24,25). Thus, bullfrogs serve as competent pathogen reservoirs12,22,26 and international pathogen vectors11,27,28. Traded bullfrogs have been found to carry multiple genetic lineages of Bd, including the invasive hypervirulent global pandemic lineage (BdGPL) responsible for the decline of amphibians on several continents11,29,30 as well as genotypes from enzootic lineages that tend to be less virulent due to long-term coevolution with their amphibian hosts, as reported for BdASIA-2/BdBRAZIL31,32, BdCAPE, BdCH29 and BdASIA-111.

The combination of bullfrog-Bd interactions and bullfrog farming practices may create ideal conditions for outbreaks of chytridiomycosis and declines of native fauna. High Bd tolerance in bullfrogs suggests that farms could function as reservoirs of particularly virulent Bd genotypes. Novel Bd genotypes introduced through trade may also have higher virulence in native hosts compared to endemic genotypes with which they have co-evolved. Furthermore, translocation of divergent Bd lineages may bring previously isolated Bd genotypes into contact, allowing for the emergence of sexually hybridized strains and possibly giving rise to more virulent hybrids than the parental lineages28,32. In addition, the water used in bullfrog farming is often released into the natural environment without treatment, potentially carrying viable Bd zoospores that could infect native amphibian species. Bullfrogs are also farmed in high densities, providing ideal conditions for (i) Bd transmission among hosts, (ii) hybridization to occur, and (iii) selection to operate. Thus, bullfrog farms may promote continuous spillover of pathogenic zoospores into native anuran communities and propogation of virulent Bd genotypes, but previous studies have not quantified pathogen outflow from bullfrog farms or compared the virulence of genotypes carried by farmed and native host species under controlled experimental conditions.

In Brazil, ranaculture began in the 1930s and escalated throughout the country in the 1970s33, coinciding with sharp historical amphibian declines throughout the Brazilian Atlantic Forest, which have been recently attributed to the emergence of Bd34. Recently, a Brazilian ordinance proposed that introduced aquatic species, including bullfrogs, should be considered native to foster aquaculture development35, which is expected to exacerbate introductions of bullfrogs into native amphibian communities. To test bullfrog farms in Brazil as a source of virulent Bd genotypes for native amphibians, we sampled Bd from farmed bullfrogs across life stages (tadpoles, juveniles and adults), as well as from the outflow water from frog farms. We also isolated and genotyped Bd strains found in farmed bullfrogs and performed infection experiments testing for differences in Bd virulence among isolates from native hosts and farmed bullfrogs. We hypothesized that because of the high Bd tolerance in bullfrogs23,26, isolates from bullfrogs would be more pathogenic to a native Brazilian host species than those isolated from native frogs. Combined, our results provide important information about the dynamics of Bd in the amphibian trade that should be used to guide actions directed at conserving native anurans in Brazil and elsewhere.

Results

Bullfrog farm assessment

The observed Bd prevalence in farmed tadpoles ranged from 0 (2 farms) to 48%. We detected Bd+ juveniles at all farms, with prevalence reaching 100% and infection loads up to 71,029 zoospore genome equivalents (g.e.; Fig. 1a, Table 1). Bd prevalence varied among amphibian developmental stages (F = 4.226; P = 0.027), with juveniles exhibiting higher Bd prevalence than tadpoles (Tukey P = 0.02), but not adults (Fig. 2a). We detected the same pattern after accounting for a conservative estimate of false positive error (Supplementary Table S1, Table S2). Prevalence of Bd was similar across life stage when accounting for a conservative estimate of false negative error (Supplementary Table S1, Table S2). Juveniles presented higher infection loads than adults (t = 6.55, df = 1, P < 0.001; Fig. 2b).

Figure 1
figure 1

Bd prevalence in different developmental stages (tadpole, juvenile and adult) in different sampled farms (a); Bd lineages isolated from bullfrogs (b).

Table 1 Developmental stage of bullfrogs, Bd prevalence [presented as percentage (Bd+/tested individuals)], and infection load [values presented as mean; SD (range)] sampled in farms.
Figure 2
figure 2

Bd prevalence (%) among the three developmental stages (tadpole, juvenile and adult) (a); Bd infection load by developmental stage (juvenile and adult) (b).

We consistently recorded high Bd zoospore concentrations in outflow water, with an average concentration of 114 Bd zoospore g.e. per liter. Outflow of water averaged 60,000 L per day, which is estimated to release approximately 50 million zoospore g.e. per day (Fig. 3).

Figure 3
figure 3

Concentration of Bd zoospores per liter of water that is released from the farms into the natural surrounding waterbodies (grey bars) and volume of outflow water released daily by the farms (black bars).

We cultured 36 Bd isolates [7 from farm #4, 7 from farm #8, 8 from farm #9, and 14 from farm #10 (Fig. 1b)] from tadpoles showing clinical signs of chytridiomycosis (dekeratinized jaw sheath and tooth rows). We detected isolates from BdGPL-2 at all four farms and isolates from BdASIA-2/BdBRAZIL at two farms. One tadpole was infected with both lineages.

Laboratory infection experiment

We found a significant effect of Bd isolate on survival of the native toadlet Brachycephalus ephippium2 = 37.269; df = 5; P < 0.0001; Fig. 4a). Isolates associated with the highest and lowest host mortality rates were both cultured from bullfrogs (Fig. 4a). However, Bd isolates cultured from bullfrogs and native frogs led to similar rates of host mortality (χ2 = 0.082; df = 1; P = 0.774).

Figure 4
figure 4

Brachycephalus ephippium survival curves (%) following inoculation with six Bd isolates (a); Mean infection load (log; ±SE) on day 16 (including of animals that died before day 16) among treatments (b).

Patterns in Bd infection loads were consistent with our survival analysis. Average pathogen loads differed among B. ephippium exposed to different Bd isolates (F(6,49) = 68.918; P < 0.0001; Table 2, Fig. 4b, Supplementary Table S3). However, loads were similar between toadlets exposed to isolates cultured from bullfrogs and native frogs (F(1,46) = 0.002; P = 0.965). Individuals that died during the experiment showed similarly high Bd infection loads (>105 g.e.), independent of isolate (F(5,35) = 1.674; P = 0.168; Table 3).

Table 2 Infection loads by Bd isolate during the experiment (day 16 and day of mortality); values presented as mean; SD (range; sample size). C = control; N = isolates from native amphibian hosts; F = isolates from bullfrogs.
Table 3 Proportional Hazard analysis testing the interactive effects of Bd infection load and Bd isolates on the survival of B. ephippium.

Discussion

Past studies have shown that invasive bullfrog populations harbor high prevalences of Bd and may contribute to the global spread of Bd through the international food trade8,11,12. However, we lack information on how frog farming practices may facilitate the spread of Bd to nearby native host communities and whether bullfrog farms function as a reservoir for highly virulent Bd strains. Our study demonstrates that bullfrog farms constantly release substantial quantities of Bd zoospores into the surrounding natural environment. The high concentration of Bd zoospores in frog pens probably results from movement of water among infected frog pens36,37. Releases of Bd zoospores from bullfrog farms may not only maintain the presence of Bd in the natural environment, but also may introduce the pathogen to new sites.

The high prevalence and infection loads observed in farmed bullfrogs can be explained by the high densities of frogs in these farms, which promotes pathogen transmission by direct contact among individuals38 or by circulating through frog pens36,37. Ideally, the density of animals should not surpass 20 tadpoles per liter of water, or 100 juveniles/50 adults per m233,39,40. We observed high bullfrog densities of all frog life stages at sampled farms (Supplementary Fig. S1). Thus, stress caused by high host densities in captivity may lower host immune capacity and increase pathogen loads39.

We found lower Bd prevalence in bullfrog tadpoles than in juveniles and adults. Prevalence of Bd in farmed bullfrog tadpoles approximated or exceeded prevalences found in wild tadpoles in the Brazilian Atlantic Forest and savannah (Cerrado)34, the biomes where our sampled farms were located. Although visual inspection can be used as a proxy for Bd diagnosis34,36,38,41,42, it is possible that mouthpart dekeratinization may not be observed at early stages of infection36; thus, the prevalence we detected in bullfrog farms may be an underestimate. When we applied the error estimation to our raw data, Bd prevalence increased significantly, making our findings even worse. The above method uses different species and environmental conditions, and considers individuals Bd+ with zoospore genomic equivalents (g.e.) ≥0.1, while the present study ≥1 g.e. This may be generating an overestimated false negative value, however, considering the legitimacy of these data, the prevalence in bullfrog tadpoles may be higher than observed, and therefore, therefore we are most concerned about the high prevalence found in the bullfrog farms. Hence, we argue that tadpoles in bullfrog farms may play a key role as reservoirs of the pathogen43, maintaining high concentrations of zoospores in the water and promoting infection or reinfection of individuals.

We showed that juveniles had higher Bd prevalence compared to tadpoles, and higher infection loads compared to adults. Our findings corroborate other studies showing low infection loads in adults43 and juveniles with higher prevalence, infection loads, and mortality rates compared to adults44. During metamorphosis, tadpoles undergo several physiological processes that reshape the immune system, leading to stress that can increase susceptibility to Bd infection45,46. In addition, newly metamorphosed individuals may be more affected by Bd infection, since this is the stage at which keratinization of the skin occurs, providing new substrate and making them an excellent host for keratinophilic pathogens such as Bd47. In the wild, high prevalence and infection loads of juveniles may result in greater spread of Bd, since this stage of development is associated with high dispersion rates48,49.

The presence of Bd at all farms and host life stages suggests that this pathogen has the potential to interfere economically in ranaculture by negatively influencing the commercial production of bullfrogs. Although this species tolerates Bd infection23,26, previous study shows that bullfrog tadpoles displayed multiple cardiac alterations in response to infection50, possibly incurring a high energy cost to developing animals, affecting metamorphosis and potentially reducing growth and post-metamorphic survival. Furthermore, other sublethal effects of Bd have been reported in tadpoles, such as behavioral changes and reduced feeding51. Chytrid infection may also reduce frog immune responses, increasing susceptibility to other infections52. Therefore, controlling chytrid infections in frog farms would potentially increase profits, due to enhanced growth and quality of frogs produced.

In addition to serving as a reservoir of the pathogen, bullfrogs may host high genetic diversity of Bd53. The presence of two Bd lineages within our focal farms, BdGPL and BdASIA-2/BdBRAZIL, even in the same tadpole, creates a high potential for hybridization. Hybrid lineages, such as the one between BdGPL × BdASIA-2/BdBRAZIL in Brazil’s Atlantic forest28,54, may be more virulent than the parental lineages32, in accordance with the hybrid vigor theory55. Hybridization of Bd lineages is apparently rare in the wild11,54, but frog farms may exacerbate the risk of new hybrids32,56.

The evolution of virulence is expected in bullfrog farms due to short generation times of Bd57, high host tolerance23,26, and high host population densities38. However, we found no support for higher virulence of isolates originating from bullfrog farms compared to those originating from wild frogs. Instead, individuals of B. ephippium died when Bd infection load surpassed 100,000 zoospores g.e., independent of isolate or origin (bullfrog or native host)58. In contrast, an infection load threshold of 10,000 zoospores was considered to be fatal for aquatic-breeding species59,60, indicating that these thresholds vary among species with contrasting life histories. The high variation in Bd load that we observed between host individuals could be attributed to multiple factors that can act alone or synergistically. Variation in innate and adaptive host immune responses61,62, gene expression63,64, previous contact with the pathogen65 and host stress levels66 may modulate the intensity of infection. Also, infections from different strains may result in different Bd loads32.

Our study emphasizes the prominent role that bullfrog farms may play in the spread of Bd into native host communities. In addition, it suggests that Bd strains that evolve on bullfrog farms are likely to impact susceptible native species through the release of zoospores into the environment, a form of pathogen spread. Our results indicate that bullfrog farms, which are distributed across the world, provide a potential environment for survival, reproduction and dissemination of Bd inside and outside farms. A possible way to reduce the spread of chytrid fungus through farming would be to treat the water that is released to the environment. Bullfrogs should also be treated for Bd infection using previously reported treatment methods67. Additional studies aimed at treatment and control of Bd in amphibian farming systems would be beneficial to native amphibian conservation efforts.

Materials and Methods

Bullfrog farm assessment

We sampled 10 bullfrog farms in the state of São Paulo, located in Brazilian Atlantic Forest and savannah (Cerrado), southeastern Brazil. Our sampling quantified Bd prevalence and infection loads in farmed bullfrogs. We sampled 35 juveniles from each of eight farms and 35 adults from each of nine farms. For each individual, we swabbed the inguinal region and each limb 5 times, specifically between the digits68. We used disposable gloves to handle each individual frog to avoid cross-contamination. Additionally, we sampled 100 tadpoles from each of nine farms for the presence of Bd (Supplementary Fig. S2a). We selected tadpoles at Gosner’s stage 25 and visually inspected jaws sheaths and tooth rows. We did not use molecular methods of detection for tadpoles because dekeratinization of the oral region is a reliable proxy for Bd infection in tadpoles36,41, including amphibians sampled in Brazil34,69. We considered individuals with jaw sheath dekeratinization to be infected with Bd (Bd+), regardless of the level of dekeratinization of the tooth rows (Supplementary Fig. S2). Although visual inspection for mouthpart dekeratinization is used for Bd diagnosis in tadpoles34,41, some studies show that tadpoles with early stages of infection may not present mouthpart dekeratinization36, and that tadpoles with fully dekeratinized mouthparts no longer have substrate for Bd growth34, therefore leading to false positives and negatives42. Thus, in order to validate our method, we ran qPCR analyses with a subsample of 51 bullfrog tadpoles with dekeratinized mouthparts from one of the focal bullfrog farms. We found a very low proportion of false positives (two individuals = 4%), indicating that visual inspection for mouthpart dekeratinization is a realiable method for Bd detection in bullfrogs from the state of São Paulo (Supplementary Table S4). In addition, and to be even more conservative, we included an estimation of false negatives and false positives in our analyses (see below) based on previous studies. We applied as a possible error the largest estimates observed in a recent study42: 12.3% for false positives and 18.7% for false negatives (Supplementary Table S1).

In addition to detecting and quantifying Bd from individual frogs, we collected water samples and standardized a filtering protocol to detect the pathogen in the aquatic environment. We measured the volume of water released by farms into the surrounding natural environment (L/s). We measured outflow rate as the amount of water released per second. Then we collected 1 liter and filtered 500 ml of this water from each farm (Supplementary Fig. S1d). Using a vacuum pump, we filtered the sampled water with a permeable membrane (0.45 μm pore size) (Supplementary Fig. S3). After filtering, we extracted DNA from membranes and quantified Bd zoospores and zoosporangia using qPCR as described below.

qPCR, Bd isolation and sequencing

We extracted Bd DNA from swab samples of juveniles and adult amphibians using PrepMan ULTRA (Life Technologies) and performed quantitative PCR analyses for Bd detection and quantification68. We considered samples with zoospore genomic equivalents (g.e.) ≥1 to be Bd+70. We estimated Bd infection prevalence (for tadpoles, juveniles and adults) as the number of infected individuals divided by the total number of sampled individuals for each farm.

We cultured Bd from bullfrog tadpoles showing mouthpart dekeratinization, following protocols by Vieira & Toledo71 and Fisher et al.72. After isolation, we transferred Bd cultures into Petri dishes containing 1% Tryptone agar and incubated cultures for one week. We then extracted DNA from the culture, following protocols by James et al.73. We genotyped each Bd isolate using a sequence of 6 SNP markers (Supplementary Table S5) as described by Schloegel et al.28 and sequenced them in the Sequencing Core Lab at the University of Michigan.

Laboratory infection experiment

We conducted experimental inoculations in the laboratory to test for effects of different Bd isolates on amphibian hosts. We used Brachycephalus ephippium (Anura: Brachycephalidae) as an experimental host. Brachycephalus ephippium is a direct-developing species endemic to Brazil’s Atlantic Forest74. Direct developers carry low Bd prevalence in the wild and often show low resistance to chytridiomycosis58. Thus, they are an ideal model organism for infection trials. We collected B. ephippium (about 2 cm in snout-vent length) in the municipality of Mogi das Cruzes, state of São Paulo, Brazil. We individually housed each wild-caught individual in plastic bags with leaf-litter to avoid potential cross-contamination among frogs while in the field. We swabbed all individual frogs and only used those that tested negative for Bd in the field. During the experiment, we individually housed each frog in plastic boxes (22 × 15 × 8 cm) containing autoclaved moist Sphagnum moss. We monitored frogs daily and fed them calcium-fortified pinhead crickets. We carried out the experiment in a temperature-controlled room, with temperatures set at 20 °C and a 12 h day-night cycle.

We exposed frogs to three Bd isolates from bullfrogs sampled at farms for this study and three Bd isolates previously isolated and genotyped from native amphibian hosts. Our experimental design consisted of seven treatments (six Bd isolates and a negative control) with eight frogs per treatment (Supplementary Table S6). All isolates were within the BdGPL-2 clade28, which is the dominant form in the Brazilian Atlantic Forest54. We cultured Bd isolates in Petri dishes with tryptone agar at 17 °C for five days. We then harvested Bd zoospores by flooding Petri dishes with distilled water and waiting for approximately one hour for zoospore release from zoosporangia32,75. We then quantified zoospores in a Neubauer hemocytometer and standardized the inoculum concentration (4.6 × 106 zoospores/ml) among isolates.

For inoculations, we placed frogs in individual Petri dishes containing 1 ml of Bd inoculum (treatment) or 1 ml of autoclaved distilled water (control) for 45 minutes. This procedure occurred only once for each animal. We swabbed each individual 16 and 31 days following inoculation, which is sufficient time for multiple Bd generations76. We monitored amphibians daily and swabbed dead or dying individuals. Experimental protocols were approved by the local animal care committee (CEUA UNICAMP #4688-1/2017). We used the same DNA extraction and qPCR protocols described above to detect and quantify Bd68.

Statistical analyses

Bullfrog farm assessment

We used a General Linear Model (GLM) with a binomial distribution (logit link) to test for differences in Bd prevalence among host developmental stages (tadpole, juvenile and adult); we performed a Tukey HSD a posteriori test for multiple comparisons. In addition, we repeated this analysis taking into account the estimated rates of false negatives and false positives associated with the visual inspection method. We also ran a GLM with normal distribution (identity link) to test for differences in Bd infection loads among host developmental stage (juvenile and adult). We log-transformed (log10) infection load data to obtain GLMs with normally distributed residuals.

Laboratory infection experiment

We built survival curves (Parametric Survival analyses) to test for effects of isolation source (native frogs vs. bullfrogs) on the survival of individual hosts. We performed a similar analysis to test for effects of isolate on host survival. In addition, we used a Generalized Linear Mixed Model (GLMM) to test whether Bd infection loads in B. epphipium exposed to bullfrog Bd isolates were higher than those exposed to isolates from native frogs, including the six Bd isolates as a random effect. For these analyses, we included samples collected halfway through the experiment (day 16) and included swabs of individuals that died before day 16; we excluded the control group from this analysis. We also performed a Tukey HSD a posteriori test for pairwise multiple comparisons among means. We also used a simple Analysis of Variance (ANOVA) to test for differences in average infection loads at the point of mortality among individuals exposed to different isolates. Finally, we used Proportional Hazards analysis, including the interaction between infection load and isolates, to test whether survival depended on these factors. We also excluded the control group from this analysis.

Ethics statement

All experiments were performed in accordance with university guidelines and regulations for animal care and husbandry. Our collecting permit was provided by ICMBio (SISBio #54656-3; 27745-13; 17242-3). Experimental protocols were approved by Universidade Estadual de Campinas (UNICAMP) and the local animal care committee [Comissão de Ética no Uso de Animal – CEUA (#4688-1/2017)]. This research was accessioned at SISGen platform (SISGEN #A1E9E10).