Gene expression varies within and between enzootic and epizootic lineages of Batrachochytrium dendrobatidis (Bd) in the Americas.

While much research focus is paid to hypervirulent fungal lineages during emerging infectious disease outbreaks, examining enzootic pathogen isolates can be equally fruitful in delineating infection dynamics and determining pathogenesis. The fungal pathogen of amphibians, Batrachochytrium dendrobatidis (Bd), exhibits markedly different patterns of disease in natural populations, where it has caused massive amphibian declines in some regions, yet persists enzootically in others. Here we compare in vitro gene expression profiles of a panel of Bd isolates representing both the enzootic Bd-Brazil lineage, and the more recently diverged, panzootic lineage, Bd-GPL. We document significantly different lineage-specific and intralineage gene expression patterns, with Bd-Brazil upregulating genes with aspartic-type peptidase activity, and Bd-GPL upregulating CBM18 chitin-binding genes, among others. We also find pronounced intralineage variation in membrane integrity and transmembrane transport ability within our Bd-GPL isolates. Finally, we highlight unexpectedly divergent expression profiles in sympatric panzootic isolates, underscoring microgeographic functional variation in a largely clonal lineage. This variation in gene expression likely plays an important role in the relative pathogenesis and host range of Bd-Brazil and Bd-GPL isolates. Together, our results demonstrate that functional genomics approaches can provide information relevant to studies of virulence evolution within the Bd clade.

Land cover and forest connectivity alter the interactions among host, pathogen and skin microbiome.

Deforestation has detrimental consequences on biodiversity, affecting species interactions at multiple scales. The associations among vertebrates, pathogens and their commensal/symbiotic microbial communities (i.e. microbiomes) have important downstream effects for biodiversity conservation, yet we know little about how deforestation contributes to changes in host microbial diversity and pathogen abundance. Here, we tested the effects of landcover, forest connectivity and infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd) on amphibian skin bacterial diversity along deforestation gradients in Brazilian landscapes. If disturbance to natural habitat alters skin microbiomes as it does in vertebrate host communities, then we would expect higher host bacterial diversity in natural forest habitats. Bd infection loads are also often higher in these closed-canopy forests, which may in turn impact skin-associated bacterial communities. We found that forest corridors shaped composition of host skin microbiomes; high forest connectivity predicted greater similarity of skin bacterial communities among host populations. In addition, we found that host skin bacterial diversity and Bd loads increased towards natural vegetation. Because symbiotic bacteria can potentially buffer hosts from Bd infection, we also evaluated the bi-directional microbiome-Bd link but failed to find a significant effect of skin bacterial diversity reducing Bd infections. Although weak, we found support for Bd increasing bacterial diversity and/or for core bacteria dominance reducing Bd loads. Our research incorporates a critical element in the study of host microbiomes by linking environmental heterogeneity of landscapes to the host-pathogen-microbiome triangle.

Genetic diversity and gene flow decline with elevation in montane mayflies.

Montane environments around the globe are biodiversity 'hotspots' and important reservoirs of genetic diversity. Montane species are also typically more vulnerable to environmental change than their low-elevation counterparts due to restricted ranges and dispersal limitations. Here we focus on two abundant congeneric mayflies (Baetis bicaudatus and B. tricaudatus) from montane streams over an elevation gradient spanning 1400 m. Using single-nucleotide polymorphism genotypes, we measured population diversity and vulnerability in these two species by: (i) describing genetic diversity and population structure across elevation gradients to identify mechanisms underlying diversification; (ii) performing spatially explicit landscape analyses to identify environmental drivers of differentiation; and (iii) identifying outlier loci hypothesized to underlie adaptive divergence. Differences in the extent of population structure in these species were evident depending upon their position along the elevation gradient. Heterozygosity, effective population sizes and gene flow all declined with increasing elevation, resulting in substantial population structure in the higher elevation species (B. bicaudatus). At lower elevations, populations of both species are more genetically similar, indicating ongoing gene flow. Isolation by distance was detected at lower elevations only, whereas landscape barriers better predicted genetic distance at higher elevations. At higher elevations, dispersal was restricted due to landscape effects, resulting in greater population isolation. Our results demonstrate differentiation over small spatial scales along an elevation gradient, and highlight the importance of preserving genetic diversity in more isolated high-elevation populations.

Cryptic species diversity reveals biogeographic support for the 'mountain passes are higher in the tropics' hypothesis.

The 'mountain passes are higher in the tropics' (MPHT) hypothesis posits that reduced climate variability at low latitudes should select for narrower thermal tolerances, lower dispersal and smaller elevational ranges compared with higher latitudes. These latitudinal differences could increase species richness at low latitudes, but that increase may be largely cryptic, because physiological and dispersal traits isolating populations might not correspond to morphological differences. Yet previous tests of the MPHT hypothesis have not addressed cryptic diversity. We use integrative taxonomy, combining morphology (6136 specimens) and DNA barcoding (1832 specimens) to compare the species richness, cryptic diversity and elevational ranges of mayflies (Ephemeroptera) in the Rocky Mountains (Colorado; approx. 40°N) and the Andes (Ecuador; approx. 0°). We find higher species richness and smaller elevational ranges in Ecuador than Colorado, but only after quantifying and accounting for cryptic diversity. The opposite pattern is found when comparing diversity based on morphology alone, underscoring the importance of uncovering cryptic species to understand global biodiversity patterns.

Amphibian-killing chytrid in Brazil comprises both locally endemic and globally expanding populations.

Chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), is the emerging infectious disease implicated in recent population declines and extinctions of amphibian species worldwide. Bd strains from regions of disease-associated amphibian decline to date have all belonged to a single, hypervirulent clonal genotype (Bd-GPL). However, earlier studies in the Atlantic Forest of southeastern Brazil detected a novel, putatively enzootic lineage (Bd-Brazil), and indicated hybridization between Bd-GPL and Bd-Brazil. Here, we characterize the spatial distribution and population history of these sympatric lineages in the Brazilian Atlantic Forest. To investigate the genetic structure of Bd in this region, we collected and genotyped Bd strains along a 2400-km transect of the Atlantic Forest. Bd-Brazil genotypes were restricted to a narrow geographic range in the southern Atlantic Forest, while Bd-GPL strains were widespread and largely geographically unstructured. Bd population genetics in this region support the hypothesis that the recently discovered Brazilian lineage is enzootic in the Atlantic Forest of Brazil and that Bd-GPL is a more recently expanded invasive. We collected additional hybrid isolates that demonstrate the recurrence of hybridization between panzootic and enzootic lineages, thereby confirming the existence of a hybrid zone in the Serra da Graciosa mountain range of Paraná State. Our field observations suggest that Bd-GPL may be more infective towards native Brazilian amphibians, and potentially more effective at dispersing across a fragmented landscape. We also provide further evidence of pathogen translocations mediated by the Brazilian ranaculture industry with implications for regulations and policies on global amphibian trade.

Contact zone dynamics during early stages of speciation in a chorus frog (Pseudacris crucifer).

Characterizing the genetic and behavioural consequences of contact between previously geographically isolated lineages provides insights into the mechanisms underlying diversification and ultimately speciation. The spring peeper (Pseudacris crucifer) is a widespread Nearctic chorus frog with six divergent mitochondrial DNA (mtDNA) lineages, many of which came into secondary contact during the Holocene. We examined genetics, morphology, advertisement calls and female preference for two lineages that began diverging in allopatry in the Pliocene and now overlap in southwestern Ontario, Canada. We found non-coincident clines in mtDNA and nuclear DNA, mirroring directionality of premating isolation barriers. We also found divergence in a range of traits between these two lineages, displacement in male call attributes and female preference for calls of their natal lineage in sympatry. Hybrids were morphologically distinct from both parental lineages, but hybrid male calls were acoustically intermediate. Female hybrids showed asymmetrical preference for Eastern male calls. These results considered together provide evidence of either unidirectional hybridization or selection against hybrids, potentially implying reproductive character displacement. Our work demonstrates the utility of integrated, multi-character approaches to understanding the processes of divergence and the nature of speciation.

Vicariance and marine migration in continental island populations of a frog endemic to the Atlantic Coastal forest.

The theory of island biogeography is most often studied in the context of oceanic islands where all island inhabitants are descendants from founding events involving migration from mainland source populations. Far fewer studies have considered predictions of island biogeography in the case of continental islands, where island formation typically splits continuous populations and thus vicariance also contributes to the diversity of island populations. We examined one such case on continental islands in southeastern Brazil, to determine how classic island biogeography predictions and past vicariance explain the population genetic diversity of Thoropa taophora, a frog endemic to the Atlantic Coastal Forest. We used nuclear microsatellite markers to examine the genetic diversity of coastal and island populations of this species. We found that island isolation has a role in shaping the genetic diversity of continental island species, with island populations being significantly less diverse than coastal populations. However, area of the island and distance from coast had no significant effect on genetic diversity. We also found no significant differences between migration among coastal populations and migration to and from islands. We discuss how vicariance and the effects of continued migration between coastal and island populations interact to shape evolutionary patterns on continental islands.

Wildlife disease. Recent introduction of a chytrid fungus endangers Western Palearctic salamanders.

Emerging infectious diseases are reducing biodiversity on a global scale. Recently, the emergence of the chytrid fungus Batrachochytrium salamandrivorans resulted in rapid declines in populations of European fire salamanders. Here, we screened more than 5000 amphibians from across four continents and combined experimental assessment of pathogenicity with phylogenetic methods to estimate the threat that this infection poses to amphibian diversity. Results show that B. salamandrivorans is restricted to, but highly pathogenic for, salamanders and newts (Urodela). The pathogen likely originated and remained in coexistence with a clade of salamander hosts for millions of years in Asia. As a result of globalization and lack of biosecurity, it has recently been introduced into naïve European amphibian populations, where it is currently causing biodiversity loss.

Long-term endemism of two highly divergent lineages of the amphibian-killing fungus in the Atlantic Forest of Brazil.

The recent global spread of the amphibian-killing fungus [Batrachochytrium dendrobatidis (Bd)] has been closely tied to anthropogenic activities; however, regional patterns of spread are not completely understood. Using historical samples, we can test whether Bd was a spreading or endemic pathogen in a region within a particular time frame, because those two disease states provide different predictions for the regional demographic dynamics and population genetics of Bd. Testing historical patterns of pathogen prevalence and population genetics under these predictions is key to understanding the evolution and origin of Bd. Focusing on the Atlantic Forest (AF) of Brazil, we used qPCR assays to determine the presence or absence of Bd on 2799 preserved postmetamorphic anurans collected between 1894 and 2010 and used semi-nested PCRs to determine the frequency of rRNA ITS1 haplotypes from 52 samples. Our earliest date of detection was 1894. A mean prevalence of 23.9% over time and spatiotemporal patterns of Bd clusters indicate that Bd has been enzootic in the Brazilian AF with no evidence of regional spread within the last 116 years. ITS1 haplotypes confirm the long-term presence of two divergent strains of Bd (BdGPL and Bd-Brazil) and three spatiotemporally broad genetic demes within BdGPL, indicating that Bd was not introduced into southeast Brazil by the bullfrog trade. Our data show that the evolutionary history and pathogen dynamics of Bd in Brazil is better explained by the endemic pathogen hypothesis.

Kin-bias, breeding site selection and female fitness in a cannibalistic Neotropical frog.

Resource availability influences sexual selection within populations and determines whether behaviours such as territoriality or resource sharing are adaptive. In Thoropa taophora, a frog endemic to the Atlantic Coastal Rainforest of Brazil, males compete for and defend limited breeding sites while females often share breeding sites with other females; however, sharing breeding sites may involve costs due to cannibalism by conspecific tadpoles. We studied a breeding population of T. taophora to determine (i) whether this species exhibits polygynous mating involving female choice for territorial males and limited breeding resources; (ii) whether limited breeding resources create the potential for male-male cooperation in defence of neighbouring territories; and (iii) whether females sharing breeding sites exhibit kin-biased breeding site choice, possibly driven by fitness losses due to cannibalism among offspring of females sharing sites. We used microsatellites to reconstruct parentage and quantify relatedness at eight breeding sites in our focal population, where these sites are scarce, and in a second population, where sites are abundant. We found that at localities where the appropriate sites for reproduction are spatially limited, the mating system for this species is polygynous, with typically two females sharing a breeding site with a male. We also found that females exhibit negative kin-bias in their choice of breeding sites, potentially to maximize their inclusive fitness by avoiding tadpole cannibalism of highly related kin. Our results indicate that male territorial defence and female site sharing are likely important components of this mating system, and we propose that kinship-dependent avoidance in mating strategies may be more general than previously realized.

Tracking climate change in a dispersal-limited species: reduced spatial and genetic connectivity in a montane salamander.

Tropical montane taxa are often locally adapted to very specific climatic conditions, contributing to their lower dispersal potential across complex landscapes. Climate and landscape features in montane regions affect population genetic structure in predictable ways, yet few empirical studies quantify the effects of both factors in shaping genetic structure of montane-adapted taxa. Here, we considered temporal and spatial variability in climate to explain contemporary genetic differentiation between populations of the montane salamander, Pseudoeurycea leprosa. Specifically, we used ecological niche modelling (ENM) and measured spatial connectivity and gene flow (using both mtDNA and microsatellite markers) across extant populations of P. leprosa in the Trans-Mexican Volcanic Belt (TVB). Our results indicate significant spatial and genetic isolation among populations, but we cannot distinguish between isolation by distance over time or current landscape barriers as mechanisms shaping population genetic divergences. Combining ecological niche modelling, spatial connectivity analyses, and historical and contemporary genetic signatures from different classes of genetic markers allows for inference of historical evolutionary processes and predictions of the impacts future climate change will have on the genetic diversity of montane taxa with low dispersal rates. Pseudoeurycea leprosa is one montane species among many endemic to this region and thus is a case study for the continued persistence of spatially and genetically isolated populations in the highly biodiverse TVB of central Mexico.

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication.

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.

Genetic drift and rapid evolution of viviparity in insular fire salamanders (Salamandra salamandra).

Continental islands offer an excellent opportunity to investigate adaptive processes and to time microevolutionary changes that precede macroevolutionary events. We performed a population genetic study of the fire salamander (Salamandra salamandra), a species that displays unique intraspecific diversity of reproductive strategies, to address the microevolutionary processes leading to phenotypic and genetic differentiation of island, coastal and interior populations. We used eight microsatellite markers to estimate genetic diversity, population structure and demographic parameters in viviparous insular populations and ovoviviparous coastal and interior populations. Our results show considerable genetic differentiation (F(ST) range: 0.06-0.27), and no clear signs of gene flow among populations, except between the large and admixed interior populations. We find no support for island colonization by rafting or intentional/accidental anthropogenic introductions, indicating that rising sea levels were responsible for isolation of the island populations approximately 9000 years ago. Our study provides evidence of rapid genetic differentiation between island and coastal populations, and rapid evolution of viviparity driven by climatic selective pressures on island populations, geographic isolation with genetic drift, or a combination of these factors. Studies of these viviparous island populations in early stages of divergence help us better understand the microevolutionary processes involved in rapid phenotypic shifts.

Reproductive success by large, closely related males facilitated by sperm storage in an aggregate breeding amphibian.

The outcome of sexual selection on males may depend on female mate choice and male-male competition as well as postcopulatory processes such as cryptic female choice and sperm competition. We studied the outcome of sexual selection in the spotted salamander (Ambystoma maculatum), specifically examining the role of body size and relatedness on male reproductive success. Using controlled mating experiments in the field, we gave females access to three males of different sizes. We used seven microsatellite loci to determine paternity in the resulting larvae, estimate relatedness (r) between females and their mates, and calculate md(2) (a measure of within-individual genomic divergence), heterozygosity, and standardized heterozygosity in the larvae. Both body size and relatedness to the female were significant predictors of male reproductive success. The relatedness of the males available to a female did not influence the amount of stored sperm she used to sire her larvae. Nonetheless, computer simulations showed that the average md(2), heterozygosity, and standardized heterozygosity of the offspring were lower than expected by random mating. These differences are due to the use of stored sperm to fertilize some eggs; md(2), heterozygosity, and standardized heterozygosity of larvae sired by stored sperm were significantly lower than those of larvae sired by the experimental males. These results suggest that relatedness may further influence a male's long-term reproductive success by determining whether his sperm is stored for later breeding seasons. Sexual selection in this salamander likely involves a complex interaction among many factors and may act over many seasons.

Characterization of microsatellite markers for Thoropa taophora (Anura, Cycloramphidae), a frog endemic to the Brazilian Atlantic Rain Forest.

Here we characterize nine microsatellite loci for Thoropa taophora collected from a coastal and an island population in southeastern Brazil. The number of alleles per locus ranged from three to 11 for the coastal population and from two to eight for the island population. Observed heterozygosities ranged from 0.308 to 0.778 for the coastal population and from zero to 0.737 for the island population. One locus did not conform to Hardy-Weinberg expectations on the coast, while two were in disequilibrium on the island. These markers are currently being used to examine genetic differentiation and patterns of dispersal in T. taophora.

Multiple paternity in an aggregate breeding amphibian: the effect of reproductive skew on estimates of male reproductive success.

Aggregate, or explosive, breeding is widespread among vertebrates and likely increases the probability of multiple paternity. We assessed paternity in seven field-collected clutches of the explosively breeding spotted salamander (Ambystoma maculatum) using 10 microsatellite loci to determine the frequency of multiple paternity and the number of males contributing to a female's clutch. Using the Minimum Method of allele counts, multiple paternity was evident in 70% of these egg masses. Simple allele counts underestimate the number of contributing males because this method cannot distinguish multiple fathers with common or similar alleles. Therefore, we used computer simulations to estimate from the offspring genotypes the most likely number of contributing fathers given the distributions of allele frequencies in this population. We determined that two to eight males may contribute to A. maculatum clutches; therefore, multiple paternity is a common strategy in this aggregate breeding species. In aggregate mating systems competition for mates can be intense, thus differential reproductive success (reproductive skew) among males contributing to a female's clutch could be a probable outcome. We use our data to evaluate the potential effect of reproductive skew on estimates of the number of contributing males. We simulated varying scenarios of differential male reproductive success, ranging from equal contribution to high reproductive skew among contributing sires in multiply sired clutches. Our data suggest that even intermediate levels of reproductive skew decrease confidence substantially in estimates of the number of contributing sires when parental genotypes are unknown.

Early male reproductive advantage, multiple paternity and sperm storage in an amphibian aggregate breeder.

We tested whether the order in which males encounter females affects reproductive fitness in spotted salamanders (Ambystoma maculatum). Using mating chambers in the field, we allowed one male access to a female before a second male. We then used four microsatellite markers in paternity analyses of the resulting larvae. First males sired a significantly larger number of offspring than second males, suggesting that male reproductive success is greatly enhanced by early arrival at breeding ponds. Multiple paternity was common among clutches, and frequently larvae were assigned to unidentified males that had not been in the chambers. Sperm from these males had either been stored by females for a year or obtained more recently at other breeding sites.

The evolution of alternative reproductive strategies: fitness differential, heritability, and genetic correlation between the sexes.

Paternity analyses using molecular markers have become standard in studies of mating systems, parentage, and kinship. In systems where individuals exhibit alternative mating strategies, molecular analyses have been productively used to estimate the reproductive success of each behavioral type and hence the fitness consequences to each individual. Here we review the fitness results in a system of five alternative mating strategies present in one population of side-blotched lizards (Uta stansburiana). Males in this population adopt one of three behavioral strategies that differ in their degree of territoriality and mate guarding. In contrast, females adopt one of two strategies that differ in offspring quantity and quality. We use paternity analyses to estimate the fitness of each morph, the heritability of reproductive strategy, and the correlation in strategy between the sexes and discuss the implications of our findings for the evolution and maintenance of reproductive polymorphism in this and other systems.

Polygyny, mate-guarding, and posthumous fertilization as alternative male mating strategies.

Alternative male mating strategies within populations are thought to be evolutionarily stable because different behaviors allow each male type to successfully gain access to females. Although alternative male strategies are widespread among animals, quantitative evidence for the success of discrete male strategies is available for only a few systems. We use nuclear microsatellites to estimate the paternity rates of three male lizard strategies previously modeled as a rock-paper-scissors game. Each strategy has strengths that allow it to outcompete one morph, and weaknesses that leave it vulnerable to the strategy of another. Blue-throated males mate-guard their females and avoid cuckoldry by yellow-throated "sneaker" males, but mate-guarding is ineffective against aggressive orange-throated neighbors. The ultradominant orange-throated males are highly polygynous and maintain large territories; they overpower blue-throated neighbors and cosire offspring with their females, but are often cuckolded by yellow-throated males. Finally, yellow-throated sneaker males sire offspring via secretive copulations and often share paternity of offspring within a female's clutch. Sneaker males sire more offspring posthumously, indicating that sperm competition may be an important component of their strategy.