Aquatic Soundscape Recordings Reveal Diverse Vocalizations and Nocturnal Activity of an Endangered Frog.

AbstractAutonomous sensors provide opportunities to observe organisms across spatial and temporal scales that humans cannot directly observe. By processing large data streams from autonomous sensors with deep learning methods, researchers can make novel and important natural history discoveries. In this study, we combine automated acoustic monitoring with deep learning models to observe breeding-associated activity in the endangered Sierra Nevada yellow-legged frog (Rana sierrae), a behavior that current surveys do not measure. By deploying inexpensive hydrophones and developing a deep learning model to recognize breeding-associated vocalizations, we discover three undocumented R. sierrae vocalization types and find an unexpected temporal pattern of nocturnal breeding-associated vocal activity. This study exemplifies how the combination of autonomous sensor data and deep learning can shed new light on species' natural history, especially during times or in locations where human observation is limited or impossible.

The microbiome at the interface between environmental stress and animal health: an example from the most threatened vertebrate group.

Nitrate pollution and global warming are ubiquitous stressors likely to interact and affect the health and survival of wildlife, particularly aquatic ectotherms. Animal health is largely influenced by its microbiome (commensal/symbiotic microorganisms), which responds to such stressors. We used a crossed experimental design including three nitrate levels and five temperature regimes to investigate their interactive and individual effects on an aquatic ectotherm, the European common frog. We associated health biomarkers in larvae with changes in gut bacteria diversity and composition. Larvae experienced higher stress levels and lower body condition under high temperatures and nitrate exposure. Developmental rate increased with temperature but decreased with nitrate pollution. Alterations in bacteria composition but not diversity are likely to correlate with the observed outcomes in larvae health. Leucine degradation decreased at higher temperatures corroborating accelerated development, nitrate degradation increased with nitrate level corroborating reduced body condition and an increase in lysine biosynthesis may have helped larvae deal with the combined effects of both stressors. These results reinforce the importance of associating traditional health biomarkers with underlying microbiome changes. Therefore, we urge studies to investigate the effects of environmental stressors on microbiome composition and consequences for host health in a world threatened by biodiversity loss.

Whole-genome sequencing and analysis of Chryseobacterium arthrosphaerae from Rana nigromaculata.

Chryseobacterium arthrosphaerae strain FS91703 was isolated from Rana nigromaculata in our previous study. To investigate the genomic characteristics, pathogenicity-related genes, antimicrobial resistance, and phylogenetic relationship of this strain, PacBio RS II and Illumina HiSeq 2000 platforms were used for the whole genome sequencing. The genome size of strain FS91703 was 5,435,691 bp and GC content was 37.78%. A total of 4,951 coding genes were predicted; 99 potential virulence factors homologs were identified. Analysis of antibiotic resistance genes revealed that strain FS91703 harbored 10 antibiotic resistance genes in 6 categories and 2 multidrug-resistant efflux pump genes, including adeG and farA. Strain FS91703 was sensitive to ß-lactam combination drugs, cephem, monobactam and carbapenems, intermediately resistant to phenicol, and resistant to penicillin, aminoglycosides, tetracycline, fluoroquinolones, and folate pathway inhibitors. Phylogenetic analysis revealed that strain FS91703 and C. arthrosphaerae CC-VM-7T were on the same branch of the phylogenetic tree based on 16 S rRNA; the ANI value between them was 96.99%; and the DDH values were 80.2, 72.2 and 81.6% by three default calculation formulae. These results suggested that strain FS91703 was a species of C. arthrosphaerae. Pan-genome analysis showed FS91703 had 566 unique genes compared with 13 other C. arthrosphaerae strains, and had a distant phylogenetic relationship with the other C. arthrosphaerae strains of the same branch in phylogenetic tree based on orthologous genes. The results of this study suggest that strain FS91703 is a multidrug-resistant and highly virulent bacterium, that differs from other C. arthrosphaerae strains at the genomic level. The knowledge about the genomic characteristics and antimicrobial resistance of strain FS91703 provides valuable insights into this rare species, as well as guidance for the treatment of the disease caused by FS91703 in Rana nigromaculata.

Association between the skin microbiome and MHC class II diversity in an amphibian.

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

Phylotranscriptomic relationships of the Oophila clade of green algae associated to amphibian egg masses.

Green algae usually assigned to the genus Oophila are known to colonize egg capsules of amphibian egg masses across the Nearctic and Palearctic regions. We study the phylogenetic relationships of these algae using a phylotranscriptomic data set of 76 protein-coding single-copy nuclear genes. Our data set includes novel RNAseq data for six amphibian-associated and five free-living green algae, and draft genomes of two of the latter. Within the Oophila clade (nested within Moewusinia), we find samples from two European frogs (Rana dalmatina and R. temporaria) closely related to those of the North American frog R. aurora (Oophila subclade III). An isolate from the North American R. sylvatica (subclade IV) appears to be sister to the Japanese isolate from the salamander Hynobius nigrescens (subclade J1), and subclade I algae from Ambystoma maculatum are sister to all other lineages in the Oophila clade. Two free-living algae (Chlamydomonas nasuta and Cd. pseudogloeogama) are nested within the Oophila clade, and a strain of the type species of Chlorococcum (Cc. infusionum) is related to this assemblage. Our phylotranscriptomic tree suggests that recognition of different species within the Oophila clade ("clade B" of earlier studies) is warranted, and calls for a comprehensive taxonomic revision of Moewusinia.

Rapid molecular identification of Rana dybowskii by species-specific primers.

Oviductus Ranae is the dried oviduct from Rana dybowskii, a forest frog species with medicinal, tonic, and cosmetic properties. Due to the high price and resource shortage, counterfeit varieties of Oviductus Ranae often appear in the market. However, traditional identification methods cannot accurately differentiate between Oviductus Ranae and its adulterants. In this study, a rapid molecular identification method has been established. The method involves extracting genomic DNA in just 30 s using filter paper purification, species-specific rapid polymerase chain reaction (PCR) amplification, and finally, fluorescence detection of the products. It can accurately identify Oviductus Ranae and its three common adulterants in about 30 min, making the process simple, fast, and highly specific.

Living in a multi-stressor world: nitrate pollution and thermal stress interact to affect amphibian larvae.

The interaction of widespread stressors such as nitrate pollution and increasing temperatures associated with climate change are likely to affect aquatic ectotherms such as amphibians. The metamorphic and physiological traits of amphibian larvae during the critical onset of metamorphosis are particularly susceptible to these stressors. We used a crossed experimental design subjecting Rana temporaria larvae to four constant rearing temperatures (18, 22, 26, 28 °C) crossed with three environmentally relevant nitrate concentrations (0, 50, 100 mg×L-1) to investigate the interactive and individual effects of these stressors on metamorphic (i.e., growth and development) and physiological traits (i.e., metabolism and heat tolerance) at the onset of metamorphosis. Larvae exposed to elevated nitrate concentrations and thermal stress displayed increased metabolic rates but decreased developmental rate, highlighting interactive effects of these stressors. However, nitrate pollution alone had no effect on either metamorphic or physiological traits, suggesting that detoxification processes were sufficient to maintain homeostasis but not in combination with increased rearing temperatures. Furthermore, larvae exposed to nitrate displayed diminished abilities to exhibit temperature-induced plasticity in metamorphosis timing and heat tolerance, as well as reduced acclimation capacity in heat tolerance and an increased thermal sensitivity of metabolic rate to higher temperatures. These results highlight the importance of considering the exposure to multiple stressors when investigating how natural populations respond to global change.

Analyses of mRNA Expression Levels of Pituitary Hormones, Their Hypothalamic Regulating Factors, and Receptors Involved in Metamorphosis with Special Reference to the Summer and Winter Seasons.

Bullfrog (Rana catesbeiana) larvae inhabiting the main island of Japan overwinter as preclimax animals, whereas the larvae that reached climax in summer complete metamorphosis. We analyzed the mRNA expression levels of the adenohypophyseal hormones, hypothalamic hormones, and their receptors that are involved in controlling metamorphosis in tadpoles at various developmental stages available in summer and winter in order to understand the hormonal mechanism regulating metamorphosis progression. Corticotropin-releasing factor (CRF) and thyrotropin ß-subunit (TSHß) mRNA expression was enhanced as they reached the climax stage in metamorphosing summer tadpoles, although type 2 CRF receptor (CRFR2) mRNA levels demonstrated a tendency of elevation, indicating the activation of the hypothalamo-hypophyseal axis for stimulating the release of thyroid hormone in summer. Arginine vasotocin (AVT) mRNA levels were elevated as metamorphosis progressed, but mRNA expression levels were not synchronized with those of proopiomelanocortin (POMC) and V1b-type AVT receptor (V1bR). The elevation of mRNA levels of prolactin (PRL) 1A and type 3 thyrotropin-releasing hormone receptor (TRHR3), but not of thyrotropin-releasing hormone (TRH) precursor mRNA levels, was noted in climactic tadpoles, indicating that PRL mRNA levels are not simply dependent on the expression levels of TRH precursor mRNA. In the preclimactic larvae captured in winter, which are in metamorphic stasis, mRNA levels of pituitary hormones, hypothalamic factors, and their receptors remained low or at levels similar to those of the larvae captured in summer. These results indicate the relationship between the mRNA expression of metamorphosis-related factors and the seasonal progression/stasis of metamorphosis.

Widespread occurrence of the amphibian chytrid panzootic lineage in Uruguay is constrained by climate.

The amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd) causes chytridiomycosis, a disease among the main causes of amphibian declines worldwide. However, Bd studies on Neotropical amphibians from temperate areas are scarce. We present a comprehensive survey of Bd in Uruguay, in temperate central eastern South America, carried out between 2006 and 2014. Skin swabs of 535 specimens of 21 native and exotic frogs were tested by PCR. We used individual-level data to examine the relationship between infection, climatic variables, and their effects on body condition and the number of prey items found in stomach contents. Infection was widespread in free-ranging anurans with an overall prevalence of 41.9%, detected in 15 native species, wild American bullfrogs Aquarana catesbeiana, and captive specimens of Ceratophrys ornata and Xenopus laevis. Three haplotypes of the Bd ITS region were identified in native amphibians, all belonging to the global panzootic lineage (BdGPL), of which only one was present in exotic hosts. Despite high infection frequencies in different anurans, we found no evidence of morbidity or mortality attributable to chytridiomycosis, and we observed no discernible impact on body condition or consumed prey. Climatic conditions at the time of our surveys suggested that the chance of infection is associated with monthly mean temperature, mean humidity, and total precipitation. Temperatures below 21°C combined with moderate humidity and pronounced rainfall may increase the likelihood of infection. Multiple haplotypes of BdGPL combined with high frequencies of infection suggest an enzootic pattern in native species, underscoring the need for continued monitoring.

Cogrowth advantage: Intestinal microbiota analysis of Bufo gargarizans and Rana chensinensis.

Intestinal microbiota has profound effects on host health and adaptation to environmental changes. Bufo gargarizans and Rana chensinensis coexist in the same habitat and have been paid much attention to economically because of their medicinal value. To date, no comparison of differences between single and mixed populations has been made. In our study, differences in the structure and function of the intestinal microbial of B. gargarizans and R. chensinensis in environments of single-species and mixed-species growth were investigated by high-throughput sequencing. Our results suggest that the cogrowth of B. gargarizans and R. chensinensis could lead to the decrease of the abundance of pathogenic bacteria (Bosea) and the introduction or increase of beneficial bacteria (Kaistia, Cetobacterium and Erysipelatoclostridium). The Tax4Fun-based functional predictions revealed that the level of pathways involved in the metabolism of R. chensinensis in mixed-species aquaria is greatly up-regulated. This study provides useful information for ecologists, ecosystem policy makers and wildlife conservationists to promote more effective conservation measures.

Effect of Aeromonas hydrophila infection on leptin receptor overlapping transcript expression in Rana amurensis.

The leptin receptor overlapping transcript (LepROT) has been suggested to play several roles in immunomodulatory mechanisms; however, the understanding of its role in Rana amurensis immunity is still very limited. Here, we performed hematoxylin-eosin staining, quantitative reverse-transcription polymerase chain reaction (qRT-PCR), immunofluorescence and western blotting to investigate the roles of LepROT in the immunomodulatory mechanism and the influence of its expression on the nuclear factor-kappa B (NF-κB) signaling pathway, such as the activation of IκB kinase and NF-кB, in amphibian resistance to infection with Aeromonas hydrophila (Ah). After Ah infection, the liver, lung, kidney, skin, muscle, and stomach of R. amurensis showed cell structure disturbance, bleeding, and texture abnormalities. In addition, the relative expression levels of LepROT, NF-кB, IKKα, and IKKß were all upregulated after Ah infection; however, they showed time-dependent differential expression. The NF-кB signaling pathway exhibited robust expression levels, which might be explained by the positive feedback regulation function of LepROT. Overall, this study provides a basis for further assessment of the biological functions of LepROT and highlights its role in the regulation of immune mechanisms.

Gut dysbiosis of Rana zhenhaiensis tadpoles after lead (Pb) exposure based on integrated analysis of microbiota and gut transcriptome.

Lead (Pb) is a ubiquitously detected heavy metal pollutant in aquatic ecosystems. Previous studies focused mainly on the response of gut microbiota to Pb stress, with less emphasis on gene expression in intestine, thereby limiting the information about impacts of Pb on intestinal homeostasis in amphibians. Here, microbial community and transcriptional response of intestines in Rana zhenhaiensis tadpoles to Pb exposure were evaluated. Our results showed that 10 µg/L Pb significantly decreased bacterial diversity compared to the controls by the Simpson index. Additionally, 1000 µg/L Pb exposure resulted in a significant reduction in the abundance of Fusobacteriota phylum and Cetobacterium genus but a significant expansion in Hafnia-Obesumbacterium genus. Moreover, transcriptome analysis revealed that about 90 % of the DEGs (8458 out of 9450 DEGs) were down-regulated in 1000 µg/L Pb group, mainly including genes annotated with biological functions in fatty acid degradation, and oxidative phosphorylation, while up-regulated DEGs involved in metabolism of xenobiotics by cytochrome P450. The expression of Gsto1, Gsta5, Gstt4, and Nadph showed strong correlation with the abundance of genera Serratia, Lactococcus, and Hafnia-Obesumbacterium. The findings of this study provide important insights into understanding the influence of Pb on intestinal homeostasis in amphibians.

Urea and glucose modulation during freezing exposure in three temperate frogs reveals specific targets in relation to climate.

Amphibian diversity is most prominent in the warm and humid tropical and subtropical regions across the globe. Nonetheless, amphibians also inhabit high-altitude tropical mountains and regions at medium and high latitudes, exposing them to subzero temperatures and requiring behavioural or physiological adaptations to endure freezing events. While freeze tolerance has been predominantly reported in high-latitude zones where species endure prolonged freezing (several weeks or months), less is known about mid-latitudes amphibians exposed to occasional subzero temperatures. In this study, we employed a controlled ecological protocol, subjecting three frog species from the Iberian Peninsula (Rana parvipalmata, Epidalea calamita, and Pelobates cultripes) to a 2-h exposure to temperatures of -2 °C to investigate the accumulation of urea and glucose as physiological mechanisms associated with survival at freezing temperatures. Our results revealed a moderate response in the production of cryoprotectant metabolites under experimental freezing conditions, particularly urea, with notable findings in R. parvipalmata and E. calamita and no response in P. cultripes. However, no significant alterations in glucose concentrations were observed in any of the studied frog species. This relatively weak freezing tolerance response differs from the strong response exhibited by amphibians inhabiting high latitudes and enduring prolonged freezing conditions, suggesting potential reliance on behavioural adaptations to cope with occasional freezing episodes.

Diversity, richness, and evenness of birds in Rana Resort Forest, District Kasur, Punjab, Pakistan.

This research was carried out to determine avian diversity in the Rana Resort Forest from June to August 2023. A total of 655 birds were observed; these birds belonged to 36 different species from 11 orders; 27 families and 36 genera were found in the forest. The observed avian species that belong to the Muscicapidae family had the highest number of species, while Strigidae had the lowest number of species. The highest numbers of the order Passeriformes were recorded. Among the bird species observed, 33 were residents, 1 was a winter visitor, and 2 were summer breeders. The species diversity, richness, and evenness indices are 3.5238, 5.3973, and 0.10492, respectively. According to the index results, the Rana Resort Forest is a diverse area, but the species are not evenly distributed due to some of the threats observed during the study, such as habitat destruction, poverty, and unawareness.

A frog-derived antimicrobial peptide as a potential anti-biofilm agent in combating Staphylococcus aureus skin infection.

Staphylococcus aureus (S. aureus), one of the most prevalent bacteria found in atopic dermatitis lesions, can induce ongoing infections and inflammation by downregulating the expression of host defence peptides in the skin. In addition, the emergence of the 'superbug' Methicillin-resistant S. aureus (MRSA) has made the treatment of these infections more challenging. Antimicrobial peptides (AMPs), due to their potent antimicrobial activity, limited evidence of resistance development, and potential immunomodulatory effects, have gained increasing attention as potential therapeutic agents for atopic dermatitis. In this study, we report a novel AMP, brevinin-1E-OG9, isolated from the skin secretions of Odorrana grahami, which shows potent antibacterial activity, especially against S. aureus. Based on the characteristics of the 'Rana Box', we designed a set of brevinin-1E-OG9 analogues to explore its structure-activity relationship. Brevinin-1E-OG9c-De-NH2 exhibited the most potent antimicrobial efficacy in both in vitro and ex vivo studies and attenuated inflammatory responses induced by lipoteichoic acid and heat-killed microbes. As a result, brevinin-1E-OG9c-De-NH2 might represent a promising candidate for the treatment of S. aureus skin infections.

The regulation of m6A-related proteins during whole-body freezing of the freeze-tolerant wood frog.

Rana sylvatica (also known as Boreorana sylvatica) is one of the few vertebrates that spend extreme winters showing no physiological signs of life. Up to 70% of the total body water of the wood frog freezes as extracellular ice. Survival in extreme conditions requires regulation at transcriptional and translational levels to activate prosurvival pathways. N6-methyladenosine (m6A) methylation is one of the most common RNA modifications, regulating transcript processing and translation by executing important functions that affect regulatory pathways in stress conditions. In the study, regulation of m6A-related proteins in the liver of R. sylvatica was analyzed during 24 h frozen and 8 h thaw conditions. Decreases in the activity of demethylases of 28.44 ± 0.4% and 24.1 ± 0.9% of control values in frozen and thaw tissues, respectively, were observed. Total protein levels of m6A methyltransferase complex components methyltransferase-like 14 and Wilm's tumor associated protein were increased by 1.28-fold and 1.42-fold, respectively, during freezing. Demethylase fat mass and obesity, however, showed a decreasing trend, with a significant decrease in abundance during recovery from frozen conditions. Levels of mRNA degraders YTHDF2 and YTHDC2 also decreased under stress. Overall, increased levels of m6A methylation complex components, and suppressed levels of readers/erasers, provide evidence for the potential role of RNA methylation in freezing survival and its regulation in a hypometabolic state.

A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK.

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

Chromosome-level genome assembly of a high-altitude-adapted frog (Rana kukunoris) from the Tibetan plateau provides insight into amphibian genome evolution and adaptation.

BACKGROUND: The high-altitude-adapted frog Rana kukunoris, occurring on the Tibetan plateau, is an excellent model to study life history evolution and adaptation to harsh high-altitude environments. However, genomic resources for this species are still underdeveloped constraining attempts to investigate the underpinnings of adaptation. RESULTS: The R. kukunoris genome was assembled to a size of 4.83 Gb and the contig N50 was 1.80 Mb. The 6555 contigs were clustered and ordered into 12 pseudo-chromosomes covering ~ 93.07% of the assembled genome. In total, 32,304 genes were functionally annotated. Synteny analysis between the genomes of R. kukunoris and a low latitude species Rana temporaria showed a high degree of chromosome level synteny with one fusion event between chr11 and chr13 forming pseudo-chromosome 11 in R. kukunoris. Characterization of features of the R. kukunoris genome identified that 61.5% consisted of transposable elements and expansions of gene families related to cell nucleus structure and taste sense were identified. Ninety-five single-copy orthologous genes were identified as being under positive selection and had functions associated with the positive regulation of proteins in the catabolic process and negative regulation of developmental growth. These gene family expansions and positively selected genes indicate regions for further interrogation to understand adaptation to high altitude. CONCLUSIONS: Here, we reported a high-quality chromosome-level genome assembly of a high-altitude amphibian species using a combination of Illumina, PacBio and Hi-C sequencing technologies. This genome assembly provides a valuable resource for subsequent research on R. kukunoris genomics and amphibian genome evolution in general.

Widespread amphibian Perkinsea infections associated with Ranidae hosts, cooler months and Ranavirus co-infection.

Amphibians suffer from large-scale population declines globally, and emerging infectious diseases contribute heavily to these declines. Amphibian Perkinsea (Pr) is a worldwide anuran pathogen associated with mass mortality events, yet little is known about its epidemiological patterns, especially in comparison to the body of literature on amphibian chytridiomycosis and ranavirosis. Here, we establish Pr infection patterns in natural anuran populations and identify important covariates including climate, host attributes and co-infection with Ranavirus (Rv). We used quantitative (q)PCR to determine the presence and intensity of Pr and Rv across 1234 individuals sampled throughout central Florida in 2017-2019. We then implemented random forest ensemble learning models to predict infection with both pathogens based on physiological and environmental characteristics. Perkinsea infected 32% of all sampled anurans, and Pr prevalence was significantly elevated in Ranidae frogs, cooler months, metamorphosed individuals and frogs co-infected with Rv, while Pr intensity was significantly higher in ranid frogs and individuals collected dead. Ranavirus prevalence was 17% overall and was significantly higher in Ranidae frogs, metamorphosed individuals, locations with higher average temperatures, and individuals co-infected with Pr. Perkinsea prevalence was significantly higher than Rv prevalence across months, regions, life stages and species. Among locations, Pr prevalence was negatively associated with crayfish prevalence and positively associated with relative abundance of microhylids, but Rv prevalence did not associate with any tested co-variates. Co-infections were significantly more common than single infections for both pathogens, and we propose that Pr infections may propel Rv infections because seasonal Rv infection peaks followed Pr infection peaks and random forest models found Pr intensity was a leading factor explaining Rv infections. Our study elucidates epidemiological patterns of Pr in Florida and suggests that Pr may be under-recognized as a cause of anuran declines, especially in the context of pathogen co-infection.

Chronic Exposure to a PFAS Mixture Resembling AFFF-Impacted Surface Water Decreases Body Size in Northern Leopard Frogs (Rana pipiens).

Per- and polyfluoroalkyl substances (PFAS) occur in the environment as mixtures, yet mixture toxicity remains poorly understood. Aqueous film-forming foams (AFFFs) are a common source of PFAS. Our objective was to examine chronic effects of a complex PFAS mixture on amphibian growth and development. We tested toxicity of a five-chemical PFAS mixture summing to 10 µg/L and that accounts for >90% of the PFAS in AFFF-affected surface waters: perfluorooctane sulfonate (PFOS, 40%), perfluorohexane sulfonic acid (PFHxS, 30%), perflurooctanoic acid (PFOA, 12.5%), perfluorohexanoic acid (PFHxA, 12.5%), and perfluoropentanoic acid (PFPeA, 5%). We also included treatments to determine whether PFOS drove mixture toxicity and whether PFOS and mixture components act additively. We exposed Northern leopard frog (Rana pipiens) larvae through metamorphosis (∼130 d) in outdoor mesocosms. After 21 days of exposure, the larval body condition fell ∼5% relative to controls in the 4 µg/L PFOS treatment and mixtures lacking PFOS. At metamorphosis, the full 5-component 10 µg/L PFAS mixture reduced mass by 16% relative to controls. We did not observe effects on development. Our results indicate that toxicity of PFOS and other PFAS mixtures typical of AFFF sites act additively and that PFOS is not more inherently toxic than other mixture components.