CRISPR-Cas9-mediated barcode insertion into Bacillus thuringiensis for surrogate tracking.

The use of surrogate organisms can enable researchers to safely conduct research on pathogens and in a broader set of conditions. Being able to differentiate between the surrogates used in the experiments and background contamination as well as between different experiments will further improve research efforts. One effective approach is to introduce unique genetic barcodes into the surrogate genome and track their presence using the quantitative polymerase chain reaction (qPCR). In this report, we utilized the CRISPR-Cas9 methodology, which employs a single plasmid and a transformation step to insert five distinct barcodes into Bacillus thuringiensis, a well-established surrogate for Bacillus anthracis when Risk Group 1 organisms are needed. We subsequently developed qPCR assays for barcode detection and successfully demonstrated the stability of the barcodes within the genome through five cycles of sporulation and germination. Additionally, we conducted whole-genome sequencing on these modified strains and analyzed 187 potential Cas9 off-target sites. We found no correlation between the mutations observed in the engineered strains and the predicted off-target sites, suggesting this genome engineering strategy did not directly result in off-target mutations in the genome. This simple approach has the potential to streamline the creation of barcoded B. thuringiensis strains for use in future studies on surrogate genomes. IMPORTANCE: The use of Bacillus anthracis as a biothreat agent poses significant challenges for public health and national security. Bacillus anthracis surrogates, like Bacillus thuringiensis, are invaluable tools for safely understanding Bacillus anthracis properties without the safety concerns that would arise from using a virulent strain of Bacillus anthracis. We report a simple method for barcode insertion into Bacillus thuringiensis using the CRISPR-Cas9 methodology and subsequent tracking by quantitative polymerase chain reaction (qPCR). Moreover, whole-genome sequencing data and CRISPR-Cas9 off-target analyses in Bacillus thuringiensis suggest that this gene-editing method did not directly cause unwanted mutations in the genome. This study should assist in the facile development of barcoded Bacillus thuringiensis surrogate strains, among other biotechnological applications in Bacillus species.

Plant-Derived, Nodule-Specific Cysteine-Rich Peptides as a Novel Source of Biopesticides for Controlling Citrus Greening Disease.

Nodule-specific cysteine-rich (NCR) peptides, encoded in the genome of the Mediterranean legume Medicago truncatula (barrelclover), are known to regulate plant-microbe interactions. A subset of computationally derived 20-mer peptide fragments from 182 NCR peptides was synthesized to identify those with activity against the unculturable vascular pathogen associated with citrus greening disease, 'Candidatus Liberibacter asiaticus' (CLas). Grounded in a design of experiments framework, we evaluated the peptides in a screening pipeline involving three distinct assays: a bacterial culture assay with Liberibacter crescens, a CLas-infected excised citrus leaf assay, and an assay to evaluate effects on bacterial acquisition by the nymphal stage of hemipteran vector Diaphorina citri. A subset of the 20-mer NCR peptide fragments inhibits both CLas growth in citrus leaves and CLas acquisition by D. citri. Two peptides induced higher levels of D. citri mortality. These findings reveal 20-mer NCR peptides as a new class of plant-derived biopesticide molecules to control citrus greening disease.

Strain Tracking of 'Candidatus Liberibacter asiaticus', the Citrus Greening Pathogen, by High-Resolution Microbiome Analysis of Asian Citrus Psyllids.

The Asian citrus psyllid, Diaphorina citri, is an invasive insect and a vector of 'Candidatus Liberibacter asiaticus' (CLas), a bacterium whose growth in Citrus species results in huanglongbing (HLB), also known as citrus greening disease. Methods to enrich and sequence CLas from D. citri often rely on biased genome amplification and nevertheless contain significant quantities of host DNA. To overcome these hurdles, we developed a simple pretreatment DNase and filtration (PDF) protocol to remove host DNA and directly sequence CLas and the complete, primarily uncultivable microbiome from D. citri adults. The PDF protocol yielded CLas abundances upward of 60% and facilitated direct measurement of CLas and endosymbiont replication rates in psyllids. The PDF protocol confirmed our lab strains derived from a progenitor Florida CLas strain and accumulated 156 genetic variants, underscoring the utility of this method for bacterial strain tracking. CLas genetic polymorphisms arising in lab-reared psyllid populations included prophage-encoding regions with key functions in CLas pathogenesis, putative antibiotic resistance loci, and a single secreted effector. These variants suggest that laboratory propagation of CLas could result in different phenotypic trajectories among laboratories and could confound CLas physiology or therapeutic design and evaluation if these differences remain undocumented. Finally, we obtained genetic signatures affiliated with Citrus nuclear and organellar genomes, entomopathogenic fungal mitochondria, and commensal bacteria from laboratory-reared and field-collected D. citri adults. Hence, the PDF protocol can directly inform agricultural management strategies related to bacterial strain tracking, insect microbiome surveillance, and antibiotic resistance screening.

An Excised Leaf Assay to Measure Acquisition of 'Candidatus Liberibacter asiaticus' by Psyllids Associated with Citrus Huanglongbing Disease.

Huanglongbing, or citrus greening disease, is the most serious disease of citrus worldwide and is associated with plant infection by 'Candidatus Liberibacter asiaticus' (CLas) and other Liberibacter species. CLas is transmitted by Diaphorina citri, the Asian citrus psyllid, in a circulative propagative manner. Circulative propagative transmission is a complex process comprising at least three steps: movement of the pathogen into vector tissues, translocation and replication of the pathogen within the vector host, and pathogen inoculation of a new host by the vector. In this work, we describe an excised leaf CLas acquisition assay, which enables precise measurements of CLas acquisition by D. citri in a streamlined laboratory assay. Briefly, healthy fourth and fifth instar D. citri nymphs acquire CLas from excised CLas-positive leaves, where the insects also complete their developmental cycle. CLas titer in the resulting adults is measured using quantitative PCR and CLas-specific 16S rRNA gene primers. We observed positive correlations between CLas titer in each leaf replicate and the CLas titer that developed in the insects after acquisition (rs = 0.78; P = 0.0002). This simple assay could be used to detect CLas acquisition phenotypes and their underlying genotypes, facilitate assessment of plant factors that impact acquisition, and screen for compounds that interfere with CLas acquisition by delivering these compounds through the excised leaf.

Proteomic Polyphenism in Color Morphotypes of Diaphorina citri, Insect Vector of Citrus Greening Disease.

Diaphorina citri is a vector of "Candidatus Liberibacter asiaticus" (CLas), associated with citrus greening disease. D. citri exhibit at least two color morphotypes, blue and non-blue, the latter including gray and yellow morphs. Blue morphs have a greater capacity for long-distance flight and transmit CLas less efficiently as compared to non-blue morphs. Differences in physiology and immunity between color morphs of the insect vector may influence disease epidemiology and biological control strategies. We evaluated the effect of CLas infection on color morph and sex-specific proteomic profiles of D. citri. Immunity-associated proteins were more abundant in blue morphs as compared to non-blue morphs but were upregulated at a higher magnitude in response to CLas infection in non-blue insects. To test for differences in color morph immunity, we measured two phenotypes: (1) survival of D. citri when challenged with the entomopathogenic fungus Beauveria bassiana and (2) microbial load of the surface and internal microbial communities. Non-blue color morphs showed higher mortality at four doses of B. bassinana, but no differences in microbial load were observed. Thus, color morph polyphenism is associated with two distinct proteomic immunity phenotypes in D. citri that may impact transmission of CLas and resistance to B. bassiana under some conditions.

Metagenome-Guided Proteomic Quantification of Reductive Dehalogenases in the Dehalococcoides mccartyi-Containing Consortium SDC-9.

At groundwater sites contaminated with chlorinated ethenes, fermentable substrates are often added to promote reductive dehalogenation by indigenous or augmented microorganisms. Contemporary bioremediation performance monitoring relies on nucleic acid biomarkers of key organohalide-respiring bacteria, such as Dehalococcoides mccartyi (Dhc). Metagenome sequencing of the commercial, Dhc-containing consortium, SDC-9, identified 12 reductive dehalogenase (RDase) genes, including pceA (two copies), vcrA, and tceA, and allowed for specific detection and quantification of RDase peptides using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). Shotgun (i.e., untargeted) proteomics applied to the SDC-9 consortium grown with tetrachloroethene (PCE) and lactate identified 143 RDase peptides, and 36 distinct peptides that covered greater than 99% of the protein-coding sequences of the PceA, TceA, and VcrA RDases. Quantification of RDase peptides using multiple reaction monitoring (MRM) assays with 13C-/15N-labeled peptides determined 1.8 × 103 TceA and 1.2 × 102 VcrA RDase molecules per Dhc cell. The MRM mass spectrometry approach allowed for sensitive detection and accurate quantification of relevant Dhc RDases and has potential utility in bioremediation monitoring regimes.

Fitness Cost of Aflatoxin Production in Aspergillus flavus When Competing with Soil Microbes Could Maintain Balancing Selection.

Selective forces that maintain the polymorphism for aflatoxigenic and nonaflatoxigenic individuals of Aspergillus flavus are largely unknown. As soils are widely considered the natural habitat of A. flavus, we hypothesized that aflatoxin production would confer a fitness advantage in the soil environment. To test this hypothesis, we used A. flavus DNA quantified by quantitative PCR (qPCR) as a proxy for fitness of aflatoxigenic and nonaflatoxigenic field isolates grown in soil microcosms. Contrary to predictions, aflatoxigenic isolates had significantly lower fitness than did nonaflatoxigenic isolates in natural soils across three temperatures (25, 37, and 42°C). The addition of aflatoxin to soils (500 ng/g) had no effect on the growth of A. flavus Amplicon sequencing showed that neither the aflatoxin-producing ability of the fungus nor the addition of aflatoxin had a significant effect on the composition of fungal or bacterial communities in soil. We argue that the fitness disadvantage of aflatoxigenic isolates is most likely explained by the metabolic cost of producing aflatoxin. Coupled with a previous report of a selective advantage of aflatoxin production in the presence of some insects, our findings give an ecological explanation for balancing selection resulting in persistent polymorphisms in aflatoxin production.IMPORTANCE Aflatoxin, produced by the fungus Aspergillus flavus, is an extremely potent hepatotoxin that causes acute toxicosis and cancer, and it incurs hundreds of millions of dollars annually in agricultural losses. Despite the importance of this toxin to humans, it has remained unclear what the fungus gains by producing aflatoxin. In fact, not all strains of A. flavus produce aflatoxin. Previous work has shown an advantage to producing aflatoxin in the presence of some insects. Our current work demonstrates the first evidence of a disadvantage to A. flavus in producing aflatoxin when competing with soil microbes. Together, these opposing evolutionary forces could explain the persistence of both aflatoxigenic and nonaflatoxigenic strains through evolutionary time.

Proteogenomics Reveals Novel Reductive Dehalogenases and Methyltransferases Expressed during Anaerobic Dichloromethane Metabolism.

Dichloromethane (DCM) is susceptible to microbial degradation under anoxic conditions and is metabolized via the Wood-Ljungdahl pathway; however, mechanistic understanding of carbon-chlorine bond cleavage is lacking. The microbial consortium RM contains the DCM degrader "Candidatus Dichloromethanomonas elyunquensis" strain RM, which strictly requires DCM as a growth substrate. Proteomic workflows applied to DCM-grown consortium RM biomass revealed a total of 1,705 nonredundant proteins, 521 of which could be assigned to strain RM. In the presence of DCM, strain RM expressed a complete set of Wood-Ljungdahl pathway enzymes, as well as proteins implicated in chemotaxis, motility, sporulation, and vitamin/cofactor synthesis. Four corrinoid-dependent methyltransferases were among the most abundant proteins. Notably, two of three putative reductive dehalogenases (RDases) encoded within strain RM's genome were also detected in high abundance. Expressed RDase 1 and RDase 2 shared 30% amino acid identity, and RDase 1 was most similar to an RDase of Dehalococcoides mccartyi strain WBC-2 (AOV99960, 52% amino acid identity), while RDase 2 was most similar to an RDase of Dehalobacter sp. strain UNSWDHB (EQB22800, 72% amino acid identity). Although the involvement of RDases in anaerobic DCM metabolism has yet to be experimentally verified, the proteome characterization results implicated the possible participation of one or more reductive dechlorination steps and methyl group transfer reactions, leading to a revised proposal for an anaerobic DCM degradation pathway.IMPORTANCE Naturally produced and anthropogenically released DCM can reside in anoxic environments, yet little is known about the diversity of organisms, enzymes, and mechanisms involved in carbon-chlorine bond cleavage in the absence of oxygen. A proteogenomic approach identified two RDases and four corrinoid-dependent methyltransferases expressed by the DCM degrader "Candidatus Dichloromethanomonas elyunquensis" strain RM, suggesting that reductive dechlorination and methyl group transfer play roles in anaerobic DCM degradation. These findings suggest that the characterized DCM-degrading bacterium Dehalobacterium formicoaceticum and "Candidatus Dichloromethanomonas elyunquensis" strain RM utilize distinct strategies for carbon-chlorine bond cleavage, indicating that multiple pathways evolved for anaerobic DCM metabolism. The specific proteins (e.g., RDases and methyltransferases) identified in strain RM may have value as biomarkers for monitoring anaerobic DCM degradation in natural and contaminated environments.

Draft Genome Sequences of the 1,2-Dichloropropane-Respiring Dehalococcoides mccartyi Strains RC and KS.

Dehalococcoides mccartyi strains RC and KS respire toxic 1,2-dichloropropane to environmentally benign propene. Their genomes were sequenced with Ion Torrent technology, assembled, and annotated. The draft genomes of strains RC and KS were 1.50 and 1.49 Mb in size and carried 1,653 and 1,671 genes, respectively.

Phylogenomics Reveal the Dynamic Evolution of Fungal Nitric Oxide Reductases and Their Relationship to Secondary Metabolism.

Fungi expressing P450nor, an unconventional nitric oxide (NO) reducing cytochrome P450, are considered significant contributors to environmental nitrous oxide (N2O) emissions. Despite extensive efforts, fungal contributions to N2O emissions remain uncertain. For example, the majority of N2O emitted from antibiotic-amended soil microcosms is attributed to fungal activity, yet axenic fungal cultures do not couple N-oxyanion respiration to growth and these fungi produce only minor quantities of N2O. To assist in reconciling these conflicting observations and produce a benchmark genomic analysis of fungal denitrifiers, genes underlying denitrification were examined in >700 fungal genomes. Of 167 p450nor-containing genomes identified, 0, 30, and 48 also harbored the denitrification genes narG, napA, or nirK, respectively. Compared with napA and nirK, p450nor was twice as abundant and exhibited 2-5-fold more gene duplications, losses, and transfers, indicating a disconnect between p450nor presence and denitrification potential. Furthermore, cooccurrence of p450nor with genes encoding NO-detoxifying flavohemoglobins (Spearman r = 0.87, p = 1.6e-10) confounds hypotheses regarding P450nor's primary role in NO detoxification. Instead, ancestral state reconstruction united P450nor with actinobacterial cytochrome P450s (CYP105) involved in secondary metabolism (SM) and 19 (11%) p450nor-containing genomic regions were predicted to be SM clusters. Another 40 (24%) genomes harbored genes nearby p450nor predicted to encode hallmark SM functions, providing additional contextual evidence linking p450nor to SM. These findings underscore the potential physiological implications of widespread p450nor gene transfer, support the undiscovered affiliation of p450nor with fungal SM, and challenge the hypothesis of p450nor's primary role in denitrification.

Human impact on the diversity and virulence of the ubiquitous zoonotic parasite Toxoplasma gondii.

A majority of emerging infectious diseases in humans are zoonoses. Understanding factors that influence the emergence and transmission of zoonoses is pivotal for their prevention and control. Toxoplasma gondii is one of the most widespread zoonotic pathogens known today. Whereas only a few genotypes of T. gondii dominate in the Northern Hemisphere, many genotypes coexist in South America. Furthermore, T. gondii strains from South America are more likely to be virulent than those from the Northern Hemisphere. However, it is not clear what factor(s) shaped modern-day genetic diversity and virulence of T. gondii Here, our analysis suggests that the rise and expansion of farming in the past 11,000 years established the domestic cat/mouse transmission cycle for T. gondii, which has undoubtedly played a significant role in the selection of certain linages of T. gondii Our mathematical simulations showed that within the domestic transmission cycle, intermediately mouse-virulent T. gondii genotypes have an adaptive advantage and eventually become dominant due to a balance between lower host mortality and the ability to superinfect mice previously infected with a less virulent T. gondii strain. Our analysis of the global type II lineage of T. gondii suggests its Old World origin but recent expansion in North America, which is likely the consequence of global human migration and trading. These results have significant implications concerning transmission and evolution of zoonotic pathogens in the rapidly expanding anthropized environment demanded by rapid growth of the human population and intensive international trading at present and in the future.

Grape pomace compost harbors organohalide-respiring Dehalogenimonas species with novel reductive dehalogenase genes.

Organohalide-respiring bacteria have key roles in the natural chlorine cycle; however, most of the current knowledge is based on cultures from contaminated environments. We demonstrate that grape pomace compost without prior exposure to chlorinated solvents harbors a Dehalogenimonas (Dhgm) species capable of using chlorinated ethenes, including the human carcinogen and common groundwater pollutant vinyl chloride (VC) as electron acceptors. Grape pomace microcosms and derived solid-free enrichment cultures were able to dechlorinate trichloroethene (TCE) to less chlorinated daughter products including ethene. 16S rRNA gene amplicon and qPCR analyses revealed a predominance of Dhgm sequences, but Dehalococcoides mccartyi (Dhc) biomarker genes were not detected. The enumeration of Dhgm 16S rRNA genes demonstrated VC-dependent growth, and 6.55±0.64 × 108 cells were measured per µmole of chloride released. Metagenome sequencing enabled the assembly of a Dhgm draft genome, and 52 putative reductive dehalogenase (RDase) genes were identified. Proteomic workflows identified a putative VC RDase with 49 and 56.1% amino acid similarity to the known VC RDases VcrA and BvcA, respectively. A survey of 1,173 groundwater samples collected from 111 chlorinated solvent-contaminated sites in the United States and Australia revealed that Dhgm 16S rRNA genes were frequently detected and outnumbered Dhc in 65% of the samples. Dhgm are likely greater contributors to reductive dechlorination of chlorinated solvents in contaminated aquifers than is currently recognized, and non-polluted environments represent sources of organohalide-respiring bacteria with novel RDase genes.

'Candidatus Dichloromethanomonas elyunquensis' gen. nov., sp. nov., a dichloromethane-degrading anaerobe of the Peptococcaceae family.

Taxonomic assignments of anaerobic dichloromethane (DCM)-degrading bacteria remain poorly constrained but are important for understanding the microbial diversity of organisms contributing to DCM turnover in environmental systems. We describe the taxonomic classification of a novel DCM degrader in consortium RM obtained from pristine Rio Mameyes sediment. Phylogenetic analysis of full-length 16S rRNA gene sequences demonstrated that the DCM degrader was most closely related to members of the genera Dehalobacter and Syntrophobotulus, but sequence similarities did not exceed 94% and 93%, respectively. Genome-aggregate average amino acid identities against Peptococcaceae members did not exceed 66%, suggesting that the DCM degrader does not affiliate with any described genus. Phylogenetic analysis of conserved single-copy functional genes supported that the DCM degrader represents a novel clade. Growth strictly depended on the presence of DCM, which was consumed at a rate of 160±3µmolL-1 d-1. The DCM degrader attained 5.25×107±1.0×107 cells per µmol DCM consumed. Fluorescence in situ hybridization revealed rod-shaped cells 4±0.8µm long and 0.4±0.1µm wide. Based on the unique phylogenetic, genomic, and physiological characteristics, we propose that the DCM degrader represents a new genus and species, 'Candidatus Dichloromethanomonas elyunquensis'.

Detection and Diversity of Fungal Nitric Oxide Reductase Genes (p450nor) in Agricultural Soils.

UNLABELLED: Members of the Fungi convert nitrate (NO3 (-)) and nitrite (NO2 (-)) to gaseous nitrous oxide (N2O) (denitrification), but the fungal contributions to N loss from soil remain uncertain. Cultivation-based methodologies that include antibiotics to selectively assess fungal activities have limitations, and complementary molecular approaches to assign denitrification potential to fungi are desirable. Microcosms established with soils from two representative U.S. Midwest agricultural regions produced N2O from added NO3 (-) or NO2 (-) in the presence of antibiotics to inhibit bacteria. Cultivation efforts yielded 214 fungal isolates belonging to at least 15 distinct morphological groups, 151 of which produced N2O from NO2 (-) Novel PCR primers targeting the p450nor gene, which encodes the nitric oxide (NO) reductase responsible for N2O production in fungi, yielded 26 novel p450nor amplicons from DNA of 37 isolates and 23 amplicons from environmental DNA obtained from two agricultural soils. The sequences shared 54 to 98% amino acid identity with reference P450nor sequences within the phylum Ascomycota and expand the known fungal P450nor sequence diversity. p450nor was detected in all fungal isolates that produced N2O from NO2 (-), whereas nirK (encoding the NO-forming NO2 (-) reductase) was amplified in only 13 to 74% of the N2O-forming isolates using two separate nirK primer sets. Collectively, our findings demonstrate the value of p450nor-targeted PCR to complement existing approaches to assess the fungal contributions to denitrification and N2O formation. IMPORTANCE: A comprehensive understanding of the microbiota controlling soil N loss and greenhouse gas (N2O) emissions is crucial for sustainable agricultural practices and addressing climate change concerns. We report the design and application of a novel PCR primer set targeting fungal p450nor, a biomarker for fungal N2O production, and demonstrate the utility of the new approach to assess fungal denitrification potential in fungal isolates and agricultural soils. These new PCR primers may find application in a variety of biomes to assess the fungal contributions to N loss and N2O emissions.

Draft Genome Sequence of a Strictly Anaerobic Dichloromethane-Degrading Bacterium.

An anaerobic, dichloromethane-degrading bacterium affiliated with novel Peptococcaceae was maintained in a microbial consortium. The organism originated from pristine freshwater sediment collected from Rio Mameyes in Luquillo, Puerto Rico, in October 2009 (latitude 18°21'43.9″, longitude -65°46'8.4″). The draft genome sequence is 2.1 Mb and has a G+C content of 43.5%.

Sister Dehalobacter Genomes Reveal Specialization in Organohalide Respiration and Recent Strain Differentiation Likely Driven by Chlorinated Substrates.

The genomes of two closely related Dehalobacter strains (strain CF and strain DCA) were assembled from the metagenome of an anaerobic enrichment culture that reductively dechlorinates chloroform (CF), 1,1,1-trichloroethane (1,1,1-TCA) and 1,1-dichloroethane (1,1-DCA). The 3.1 Mbp genomes of strain CF (that dechlorinates CF and 1,1,1-TCA) and strain DCA (that dechlorinates 1,1-DCA) each contain 17 putative reductive dehalogenase homologous (rdh) genes. These two genomes were systematically compared to three other available organohalide-respiring Dehalobacter genomes (Dehalobacter restrictus strain PER-K23, Dehalobacter sp. strain E1 and Dehalobacter sp. strain UNSWDHB), and to the genomes of Dehalococcoides mccartyi strain 195 and Desulfitobacterium hafniense strain Y51. This analysis compared 42 different metabolic and physiological categories. The genomes of strains CF and DCA share 90% overall average nucleotide identity and >99.8% identity over a 2.9 Mbp alignment that excludes large insertions, indicating that these genomes differentiated from a close common ancestor. This differentiation was likely driven by selection pressures around two orthologous reductive dehalogenase genes, cfrA and dcrA, that code for the enzymes that reduce CF or 1,1,1-TCA and 1,1-DCA. The many reductive dehalogenase genes found in the five Dehalobacter genomes cluster into two small conserved regions and were often associated with Crp/Fnr transcriptional regulators. Specialization is on-going on a strain-specific basis, as some strains but not others have lost essential genes in the Wood-Ljungdahl (strain E1) and corrinoid biosynthesis pathways (strains E1 and PER-K23). The gene encoding phosphoserine phosphatase, which catalyzes the last step of serine biosynthesis, is missing from all five Dehalobacter genomes, yet D. restrictus can grow without serine, suggesting an alternative or unrecognized biosynthesis route exists. In contrast to D. mccartyi, a complete heme biosynthesis pathway is present in the five Dehalobacter genomes. This pathway corresponds to a newly described alternative heme biosynthesis route first identified in Archaea. This analysis of organohalide-respiring Firmicutes and Chloroflexi reveals profound evolutionary differences despite very similar niche-specific metabolism and function.

Urea is not a universal cryoprotectant among hibernating anurans: evidence from the freeze-tolerant boreal chorus frog (Pseudacris maculata).

Freeze-tolerant organisms accumulate a diversity of low molecular weight compounds to combat negative effects of ice formation. Previous studies of anuran freeze tolerance have implicated urea as a cryoprotectant in the wood frog (Lithobates sylvatica). However, a cryoprotective role for urea has been identified only for wood frogs, though urea accumulation is an evolutionarily conserved mechanism for coping with osmotic stress in amphibians. To identify whether multiple solutes are involved in freezing tolerance in the boreal chorus frog (Pseudacris maculata), we examined seasonal and freezing-induced variation in several potential cryoprotectants. We further tested for a cryoprotective role for urea by comparing survival and recovery from freezing in control and urea-loaded chorus frogs. Tissue levels of glucose, urea, and glycerol did not vary significantly among seasons for heart, liver, or leg muscle. Furthermore, no changes in urea or glycerol levels were detected with exposure to freezing temperatures in these tissues. Urea-loading increased tissue urea concentrations, but failed to enhance freezing survival or facilitate recovery from freezing in chorus frogs in this study, suggesting little role for urea as a natural cryoprotectant in this species. These data suggest that urea may not universally serve as a primary cryoprotectant among freeze-tolerant, terrestrially hibernating anurans.