Batrachochytrium salamandrivorans (Bsal) not detected in an intensive survey of wild North American amphibians.

The salamander chytrid fungus (Batrachochytrium salamandrivorans [Bsal]) is causing massive mortality of salamanders in Europe. The potential for spread via international trade into North America and the high diversity of salamanders has catalyzed concern about Bsal in the U.S. Surveillance programs for invading pathogens must initially meet challenges that include low rates of occurrence on the landscape, low prevalence at a site, and imperfect detection of the diagnostic tests. We implemented a large-scale survey to determine if Bsal was present in North America designed to target taxa and localities where Bsal was determined highest risk to be present based on species susceptibility and geography. Our analysis included a Bayesian model to estimate the probability of occurrence of Bsal given our prior knowledge of the occurrence and prevalence of the pathogen. We failed to detect Bsal in any of 11,189 samples from 594 sites in 223 counties within 35 U.S. states and one site in Mexico. Our modeling indicates that Bsal is highly unlikely to occur within wild amphibians in the U.S. and suggests that the best proactive response is to continue mitigation efforts against the introduction and establishment of the disease and to develop plans to reduce impacts should Bsal establish.

Amphibian chytrid prevalence on boreal toads in SE Alaska and NW British Columbia: tests of habitat, life stages, and temporal trends.

Tracking and understanding variation in pathogens such as Batrachochytrium dendrobatidis (Bd), the agent of amphibian chytridiomycosis which has caused population declines globally, is a priority for many land managers. However, relatively little sampling of amphibian communities has occurred at high latitudes. We used skin swabs collected during 2005-2017 from boreal toads Anaxyrus boreas (n = 248), in southeast Alaska (USA; primarily in and near Klondike Gold Rush National Historical Park [KLGO]) and northwest British Columbia (Canada) to determine how Bd prevalence varied across life stages, habitat characteristics, local species richness, and time. Across all years, Bd prevalence peaked in June and was >3 times greater for adult toads (37.5%) vs. juveniles and metamorphs (11.2%). Bd prevalence for toads in the KLGO area, where other amphibian species are rare or absent, was highest from river habitats (55.0%), followed by human-modified upland wetlands (32.3%) and natural upland wetlands (12.7%)-the same rank-order these habitats are used for toad breeding. None of the 12 Columbia spotted frogs Rana luteiventris or 2 wood frogs R. sylvatica from the study area tested Bd-positive, although all were from an area of low host density where Bd has not been detected. Prevalence of Bd on toads in the KLGO area decreased during 2005-2015. This trend from a largely single-species system may be encouraging or concerning, depending on how Bd is affecting vital rates, and emphasizes the need to understand effects of pathogens before translating disease prevalence into management actions.

Spatial capture-recapture reveals age- and sex-specific survival and movement in stream amphibians.

Life-history information sets the foundation for our understanding of ecology and conservation requirements. For many species, this information is lacking even for basic demographic rates such as survival and movement. When survival and movement estimates are available, they are often derived from mixed demographic groups and do not consider differences among life stages or sexes, which is critical, because life stages and sexes often contribute differentially to population dynamics. We used hierarchical models informed with spatial capture-mark-recapture data of Ascaphus montanus (Rocky Mountain tailed frog) in five streams and A. truei (coastal tailed frog) in one stream to estimate variation in survival and movement by sex and age, represented by size. By incorporating survival and movement into a single model, we were able to estimate both parameters with limited bias. Annual survival was similar between sexes of A. montanus [females = 0.885 (95% CI 0.614-1), males = 0.901 (0.657-1)], but was slightly higher for female A. truei [0.836 (0.560-0.993)] than for males [0.664 (0.354-0.962)]. Survival of A. montanus peaked at mid-age, suggesting that lower survival of young and actuarial senescence may influence population demographics. Our models suggest that younger A. montanus moved farther than older individuals, and that females moved farther than males in both species. Our results provide uncommon insight into age- and sex-specific rates of survival and movement that are crucial elements of life-history strategies and are important for modeling population growth and prescribing conservation actions.

Associations between environmental pollutants and larval amphibians in wetlands contaminated by energy-related brines are potentially mediated by feeding traits.

Energy production in the Williston Basin, located in the Prairie Pothole Region of central North America, has increased rapidly over the last several decades. Advances in recycling and disposal practices of saline wastewaters (brines) co-produced during energy production have reduced ecological risks, but spills still occur often and legacy practices of releasing brines into the environment caused persistent salinization in many areas. Aside from sodium and chloride, these brines contain elevated concentrations of metals and metalloids (lead, selenium, strontium, antimony and vanadium), ammonium, volatile organic compounds, hydrocarbons, and radionuclides. Amphibians are especially sensitive to chloride and some metals, increasing potential effects in wetlands contaminated by brines. We collected bed sediment and larval amphibians (Ambystoma mavortium, Lithobates pipiens and Pseudacris maculata) from wetlands in Montana and North Dakota representing a range of brine contamination history and severity to determine if contamination was associated with metal concentrations in sediments and if metal accumulation in tissues varied by species. In wetland sediments, brine contamination was positively associated with the concentrations of sodium and strontium, both known to occur in oil and gas wastewater, but negatively correlated with mercury. In amphibian tissues, selenium and vanadium were associated with brine contamination. Metal tissue concentrations were higher in tadpoles that graze compared to predatory salamanders; this suggests frequent contact with the sediments could lead to greater ingestion of metal-laden materials. Although many of these metals may not be directly linked with energy development, the potential additive or synergistic effects of exposure along with elevated chloride from brines could have important consequences for aquatic organisms. To effectively manage amphibian populations in wetlands contaminated by saline wastewaters we need a better understanding of how life history traits, species-specific susceptibilities and the physical-chemical properties of metals co-occurring in wetland sediments interact with other stressors like chloride and wetland drying.

Widespread legacy brine contamination from oil production reduces survival of chorus frog larvae.

Advances in drilling techniques have facilitated a rapid increase in hydrocarbon extraction from energy shales, including the Williston Basin in central North America. This area overlaps with the Prairie Pothole Region, a region densely populated with wetlands that provide numerous ecosystem services. Historical (legacy) disposal practices often released saline co-produced waters (brines) with high chloride concentrations, affecting wetland water quality directly or persisting in sediments. Despite the potential threat of brine contamination to aquatic habitats, there has been little research into its ecological effects. We capitalized on a gradient of legacy brine-contaminated wetlands in northeast Montana to conduct laboratory experiments to assess variation in survival of larval Boreal Chorus Frogs (Pseudacris maculata) reared on sediments from 3 local wetlands and a control source. To help provide environmental context for the experiment, we also measured chloride concentrations in 6 brine-contaminated wetlands in our study area, including the 2 contaminated sites used for sediment exposures. Survival of frog larvae during 46- and 55-day experiments differed by up to 88% among sediment sources (Site Model) and was negatively correlated with potential chloride exposure (Chloride Model). Five of the 6 contaminated wetlands exceeded the U.S. EPA acute benchmark for chloride in freshwater (860 mg/L) and all exceeded the chronic benchmark (230 mg/L). However, the Wetland Site model explained more variation in survival than the Chloride Model, suggesting that chloride concentration alone does not fully reflect the threat of contamination to aquatic species. Because the profiles of brine-contaminated sediments are complex, further surveys and experiments are needed across a broad range of conditions, especially where restoration or remediation actions have reduced brine-contamination. Information provided by this study can help quantify potential ecological threats and help land managers prioritize conservation strategies as part of responsible and sustainable energy development.

Interactive effects of wildfire, forest management, and isolation on amphibian and parasite abundance.

Projected increases in wildfire and other climate-driven disturbances will affect populations and communities worldwide, including host-parasite relationships. Research in temperate forests has shown that wildfire can negatively affect amphibians, but this research has occurred primarily outside of managed landscapes where interactions with human disturbances could result in additive or synergistic effects. Furthermore, parasites represent a large component of biodiversity and can affect host fitness and population dynamics, yet they are rarely included in studies of how vertebrate hosts respond to disturbance. To determine how wildfire affects amphibians and their parasites, and whether effects differ between protected and managed landscapes, we compared abundance of two amphibians and two nematodes relative to wildfire extent and severity around wetlands in neighboring protected and managed forests (Montana, USA). Population sizes of adult, male long-toed salamanders (Ambystoma macrodactylum) decreased with increased burn severity, with stronger negative effects on isolated populations and in managed forests. In contrast, breeding population sizes of Columbia spotted frogs (Rana luteiventris) increased with burn extent in both protected and managed protected forests. Path analysis showed that the effects of wildfire on the two species of nematodes were consistent with differences in their life history and transmission strategies and the responses of their hosts. Burn severity indirectly reduced abundance of soil-transmitted Cosmocercoides variabilis through reductions in salamander abundance. Burn severity also directly reduced C. variabilis abundance, possibly though changes in soil conditions. For the aquatically transmitted nematode Gyrinicola batrachiensis, the positive effect of burn extent on density of Columbia spotted frog larvae indirectly increased parasite abundance. Our results show that effects of wildfire on amphibians depend upon burn extent and severity, isolation, and prior land use. Through subsequent effects on the parasites, our results also reveal how changes in disturbance regimes can affect communities across trophic levels.

Polymorphism and gene organization of water buffalo MHC-DQB genes show homology to the BoLA DQB region.

In cattle (Bos taurus), there is evidence of more than 50 alleles of BoLA-DQB (bovine lymphocyte antigen DQB) that are distributed across at least five DQB loci, making this region one of the most complex in the BoLA gene family. In this study, DQB alleles were analysed for the water buffalo (Bubalus bubalis), another economically important bovine species. Twelve alleles for Bubu-DQB (Bubalis bubalis DQB) were determined by nucleotide sequence analysis. A phylogenetic analysis revealed numerous trans-species polymorphisms, with alleles from water buffalo assigned to at least three different loci (BoLA-DQB1, BoLA-DQB3 and BoLA-DQB4) that are also found in cattle. These presumptive loci were analysed for patterns of synonymous (d(S)) and non-synonymous (d(N)) substitution. Like BoLA-DQB1, Bubu-DQB1 was observed to be under strong positive selection for polymorphism. We conclude that water buffalo and cattle share the current arrangement of their DQB region because of their common ancestry.

Assessing substitution patterns, rates and homoplasy at HVRI of Steller sea lions, Eumetopias jubatus.

Despite the widely recognized incidence of homoplasy characterizing this region, the hypervariable region I (HVRI) of the mitochondrial control region is one of the most frequently used genetic markers for population genetic and phylogeographic studies. We present an evolutionary analysis of HVRI and cytochrome b sequences from a range-wide survey of 1031 Steller sea lions, Eumetopias jubatus, to quantify homoplasy and substitution rate at HVRI. Variation in HVRI was distributed across 41 variable sites in the 238-bp segment examined. All variants at HVR1 were found to be transitions. However, our analyses suggest that a minimum of 101 changes have actually occurred within HVRI with as many as 18 substitutions occurring at a single site. By including this hidden variation into our analyses, several instances of apparent long-range dispersal were resolved to be homoplasies and 8.5-12% of observed HVRI haplotypes were found to have geographic distributions descriptive of convergent molecular evolution rather than identity by descent. We estimate the rate of substitution at HVRI in Steller sea lions to be approximately 24 times that of cytochrome b with an absolute rate of HVRI substitution estimated at 27.45% per million years. These findings have direct implications regarding the utility of HVRI data to generate a variety of evolutionary genetic hypotheses.

Glacial refugia and the phylogeography of Steller's sea lion (Eumatopias jubatus) in the North Pacific.

Mitochondrial DNA sequence data were used to examine the phylogeographic history of Steller's sea lions (Eumetopias jubatus) in relation to the presence of Plio-Pleistocene insular refugia. Cytochrome b and control region sequences from 336 Steller's sea lions reveal phylogenetic lineages associated with continental refugia south of the ice sheets in North America and Eurasia. Phylogenetic analysis suggests the genetic structure of E. jubatus is the result of Pleistocene glacial geology, which caused the elimination and subsequent reappearance of suitable rookery habitat during glacial and interglacial periods. The cyclic nature of geological change produced a series of independent population expansions, contractions and isolations that had analogous results on Steller's sea lions and other marine and terrestrial species. Our data show evidence of four glacial refugia in which populations of Steller's sea lions diverged. These events occurred from approximately 60,000 to 180,000 years BP and thus preceded the last glacial maximum.

Gradual and quantum genome size shifts in the hystricognath rodents.

We assessed genome size variation by flow cytometry within and among 31 species of nine families of African and South American hystricognath rodents. Interspecific variation was extensive and genome size was relatively high among the South American radiation whereas only moderate variation and smaller estimates of genome size were observed in the African counterparts. The largest genome size, indicating tetraploidy was recorded in the South American octodontid, Tympanoctomys barrerae (16.8 pg DNA). This quantum shift in DNA content represents a novel mechanism of genome evolution in mammals. As expected in polyploid organisms, varying nucleotypic effects were observed in the dimensions of the sperm cells and lymphocytes of T. barrerae. The role of control mechanisms that influence cell dimensions in polyploid organisms is discussed.

Polymorphisms in MHC-DRA and -DRB alleles of water buffalo (Bubalus bubalis) reveal different features from cattle DR alleles.

Seventy-five individuals of Bubalus bubalis belonging to four different breeds, three of river buffalo and one of swamp buffalo, were studied for polymorphism in MHC DRB (Bubu-DRB) and DRA (Bubu-DRA) loci. Eight alleles of Bubu-DRB were found, and all alleles in the swamp type were shared with the three river breeds. All alleles sampled from the breed of European origin (Mediterranean) were present in breeds sampled in Brazil, thus variability of this locus may have been preserved to a great extent in the more recently founded Brazilian population. Bubu-DRB alleles contained higher proportions of synonymous vs. non-synonymous substitutions in the non-peptide-binding sites (PBS) region, in contrast to the pattern of variation found in BoLA-DRB3, the orthologous locus in cattle. This indicated that either the first domain exon (exon 2) of Bubu-DRB has not undergone as much recombination and/or gene conversion as in cattle alleles, or Bubu-DRB may be more ancient than BoLA-DRB3 alleles. Phylogenetic analysis of DRB alleles from Bubalus, Syncerus c. caffer, the Cape buffalo, and domestic cattle demonstrated transspecies polymorphism. Water buffalo contained two alleles of DRA that differed from each other in two amino acid positions, including one in the PBS (alpha22) that was also shared with Anoa depressicornis, the anoa. Discovery of variation in DRA was surprising as the first domain of DRA is a highly conserved polypeptide in mammals in general and especially in ruminants, where no other substitution in PBS was seen.

Historical population size change of bowhead whales inferred from DNA sequence polymorphism data.

Nucleotide sequence data from the mitochondrial control region were used from a phylogenetic context to investigate the long-term history of a population of bowhead whales (Balaena mysticetus). In addition, the coalescence time of these sequences was used to estimate the age of the inferred patterns of population size change. The results indicate that mitochondrial genetic polymorphism was not affected by a recent bottleneck that occurred near the turn of the 20th century, thereby preserving the signature of historical population size change in the mitochondrial genome. Further analysis showed that this population underwent an expansion initiated in the Middle to Late Pleistocene. As such, early Holocene changes in Arctic sea ice distribution appear to have had little influence on patterns of genetic variability in this population.

Molecular phylogeny and divergence time estimates for major rodent groups: evidence from multiple genes.

The order Rodentia contains half of all extant mammal species, and from an evolutionary standpoint, there are persistent controversies surrounding the monophyly of the order, divergence dates for major lineages, and relationships among families. Exons of growth hormone receptor (GHR) and breast cancer susceptibility (BRCA1) genes were sequenced for a wide diversity of rodents and other mammals and combined with sequences of the mitochondrial 12S rRNA gene and previously published sequences of von Willebrand factor (vWF). Rodents exhibit rates of amino acid replacement twice those observed for nonrodents, and this rapid rate of evolution influences estimates of divergence dates. Based on GHR sequences, monophyly is supported, with the estimated divergence between hystricognaths and most sciurognaths dating to about 75 MYA. Most estimated dates of divergence are consistent with the fossil record, including a date of 23 MYA for Mus-Rattus divergence. These dates are considerably later than those derived from some other molecular studies. Among combined and separate analyses of the various gene sequences, moderate to strong support was found for several clades. GHR appears to have greater resolving power than do 12S or vWF. Despite its complete unresponsiveness to growth hormone, Cavia (and other hystricognaths) exhibits a conservative rate of change in the intracellular domain of GHR.

Whence the red panda?

The evolutionary history of the red panda (Ailurus fulgens) plays a pivotal role in the higher-level phylogeny of the "bear-like" arctoid carnivoran mammals. Characters from morphology and molecules have provided inconsistent evidence for placement of the red panda. Whereas it certainly is an arctoid, there has been major controversy about whether it should be placed with the bears (ursids), ursids plus pinnipeds (seals, sea lions, walrus), raccoons (procyonids), musteloids (raccoons plus weasels, skunks, otters, and badgers [mustelids]), or as a monotypic lineage of uncertain phylogenetic affinities. Nucleotide sequence data from three mitochondrial genes and one nuclear intron were analyzed, with more complete taxonomic sampling of relevant taxa (arctoids) than previously available in analyses of primary molecular data, to clarify the phylogenetic relationships of the red panda to other arctoid carnivorans. This study provides detailed phylogenetic analyses (both parsimony and maximum-likelihood) of primary character data for arctoid carnivorans, including bootstrap and decay indices for all arctoid nodes, and three statistical tests of alternative phylogenetic hypotheses for the placement of the red panda. Combined phylogenetic analyses reject the hypotheses that the red panda is most closely related to the bears (ursids) or to the raccoons (procyonids). Rather, evidence from nucleotide sequences strongly support placement of the red panda within a broad Musteloidea (sensu lato) clade, including three major lineages (the red panda, the skunks [mephitids], and a clearly monophyletic clade of procyonids plus mustelids [sensu stricto, excluding skunks]). Within the Musteloidea, interrelationships of the three major lineages are unclear and probably are best considered an unresolved trichotomy. These data provide compelling evidence for the relationships of the red panda and demonstrate that small taxonomic sample sizes can result in misleading or possibly erroneous (based on prior modeling, as well as conflict between the results of our analyses of less and more complete data sets) conclusions about phylogenetic relationships and taxonomy.

Evidence from intron 1 of the nuclear transthyretin (Prealbumin) gene for the phylogeny of African mole-rats (Bathyergidae).

A 900- to 1100-bp fragment encompassing intron 1 of the nuclear transthyretin (prealbumin) gene was examined in 12 taxa of Old World hystricognath rodents of the families Bathyergidae, Petromuridae, Thryonomyidae, and Hystricidae. Within the Bathyergidae, Heterocephalus glaber (naked mole-rat) was basal, and the other East African species, Heliophobius argenteocinereus (silvery mole-rat), was sister to a southern African clade containing Bathyergus, Cryptomys, and Georychus (dune, common, and cape mole-rats). These results are congruent with studies using mitochondrial 12S rRNA gene sequences. A combined analysis of transthyretin and 12S rRNA data resulted in a well-supported topology with better resolution than either gene analyzed separately. These data support the findings by M. W. Allard and R. L. Honeycutt (1992, Mol. Biol. Evol. 9: 27-40) and R. L. Honeycutt (1992, Am. Sci. 80: 43-53) that complex social systems evolved independently at least twice, in the common and naked mole-rats.

Evaluating a putative bottleneck in a population of bowhead whales from patterns of microsatellite diversity and genetic disequilibria.

A size-selected Balaena mysticetus genomic library was screened for clones containing simple sequence repeat, or microsatellite, loci. A total of 11 novel loci was identified. These loci were combined with a set of 9 published loci, for a total of 20 markers, and were scored across a sample of 108 bowhead whales from the Bering-Chukchi-Beaufort Seas population of bowhead whales. Genetic variability was measured in terms of polymorphism information content values and unbiased heterozygosity. From the latter, estimates of long-term effective population size were obtained. In addition, gametic phase disequilibrium among loci was investigated. Moderate to high levels of polymorphism were found overall, and the long-term effective size estimates were large relative to total population size. Tests of heterozygosity excess (Cornuet and Luikart 1996) and allele frequency distribution (Luikart et al. 1998) indicated that the possibility of a recent genetic bottleneck in the Bering-Chukchi-Beaufort Seas population of bowhead whales is highly unlikely. However, the fact that five loci displayed a statistically significant heterozygote deficiency remains to be explained.

Genetic evidence of Dobrava virus in Apodemus agrarius in Hungary.

Using nested polymerase chain reaction, we sequenced Dobrava virus (DOB) from the rodent Apodemus agrarius in Hungary. The samples we isolated group with DOB samples previously isolated from A. flavicollis. This grouping may indicate host switching.