Gene functions of the Ambystoma altamirani skin microbiome vary across space and time but potential antifungal genes are widespread and prevalent.

Amphibian skin microbiomes can play a critical role in host survival against emerging diseases by protecting their host against pathogens. While a plethora of biotic and abiotic factors have been shown to influence the taxonomic diversity of amphibian skin microbiomes it remains unclear whether functional genomic diversity varies in response to temporal and environmental factors. Here we applied a metagenomic approach to evaluate whether seasonality, distinct elevations/sites, and pathogen presence influenced the functional genomic diversity of the A. altamirani skin microbiome. We obtained a gene catalogue of 92 107 nonredundant annotated genes and a set of 50 unique metagenome assembled genomes (MAGs). Our analysis showed that genes linked to general and potential antifungal traits significantly differed across seasons and sampling locations at different elevations. Moreover, we found that the functional genomic diversity of A. altamirani skin microbiome differed between B. dendrobatidis infected and not infected axolotls only during winter, suggesting an interaction between seasonality and pathogen infection. In addition, we identified the presence of genes and biosynthetic gene clusters (BGCs) linked to potential antifungal functions such as biofilm formation, quorum sensing, secretion systems, secondary metabolite biosynthesis, and chitin degradation. Interestingly genes linked to these potential antifungal traits were mainly identified in Burkholderiales and Chitinophagales MAGs. Overall, our results identified functional traits linked to potential antifungal functions in the A. altamirani skin microbiome regardless of variation in the functional diversity across seasons, elevations/sites, and pathogen presence. Our findings suggest that potential antifungal traits found in Burkholderiales and Chitinophagales taxa could be related to the capacity of A. altamirani to survive in the presence of Bd, although further experimental analyses are required to test this hypothesis.

Environmental RNA can distinguish life stages in amphibian populations.

Applications of environmental DNA (eDNA) analysis methods for biomonitoring have grown exponentially over the last decade and provide a wealth of new information on the distribution of species. However, eDNA methods have limited application for estimating population-level metrics. Environmental RNA (eRNA) has the potential to address ecological questions by gathering population demographic information from environmental media but may be challenging to detect and analyze. We developed gene-specific eRNA assays targeting keratin-associated genes in two focal species, American bullfrogs (Lithobates catesbeianus) and tiger salamanders (Ambystoma mavortium) to answer an important question in amphibian management: whether species detections represent breeding populations versus transitory adults. We performed an extensive laboratory validation with amphibians housed across development stages, where we collected 95 and 127 environmental samples for bullfrogs and salamanders, respectively. Both assays were highly specific to the larval stage and amplified with high sensitivity (90% in bullfrog and 88.4% in tiger salamander samples). We then applied our validated assays to multiple natural systems. When larvae were present, we found 74.1% overall detection in bullfrog field samples and 70.8% and 48.5% overall detection in field samples from ponds with A. macrodactylum and A. californiense larvae, correlating with eDNA detection rates. When only adults were present, we did not detect larvae-specific eRNA in A. macrodactylum ponds, despite high eDNA detection rates. Although much work is ahead for optimizing assay design, sampling and filtering methods, we demonstrate that eRNA can successfully be used to discern life stages with direct application for ecology and conservation management.

Managing invasive hybrids with pond hydroperiod manipulation in an endangered salamander system.

When invasive and endangered native taxa hybridize, the resulting admixture introduces novel conservation challenges. Across a large region of central California, a hybrid swarm consisting of admixed endangered California tiger salamanders (CTS) (Ambystoma californiense) and introduced barred tiger salamanders (BTS) (Ambystoma mavortium) has replaced native populations, threatening the genetic integrity of CTS and the vernal pool systems they inhabit. We employed a large-scale, genomically informed field experiment to test whether shortening breeding pond hydroperiod would favor native CTS genotypes. We constructed 14 large, seminatural ponds to evaluate the effect of hydroperiod duration on larval survival and mass at metamorphosis. We tracked changes in non-native allele frequencies with a 5237-gene exon capture array and employed a combination of custom Bayesian and generalized linear models to quantify the effect of pond duration on salamander fitness. Earlier work on this system showed hybrid superiority under many conditions and suggested that hybrids are favored in human-modified ponds with artificially long hydroperiods. Consistent with these earlier studies, we found overwhelming evidence for hybrid superiority. Very short hydroperiods substantially reduced the mass (1.1-1.5 fold) and survival probability (10-13 fold) of both native and hybrid larvae, confirming that hydroperiod likely exerts a strong selective pressure in the wild. We identified 86 genes, representing 1.8% of 4723 screened loci, that significantly responded to this hydroperiod-driven selection. In contrast to earlier work, under our more natural experimental conditions, native CTS survival and size at metamorphosis were always less than hybrids, suggesting that hydroperiod management alone will not shift selection to favor native larval genotypes. However, shortening pond hydroperiod may limit productivity of hybrid ponds, complementing other strategies to remove hybrids while maintaining vernal pool ecosystems. This study confirms and expands on previous work that highlights the importance of hydroperiod management to control invasive aquatic species.

Manejo de híbridos invasores mediante la manipulación del hidroperiodo de los estanques en el sistema de una salamandra en peligro de extinción Resumen La hibridación entre un taxón nativo en peligro y uno invasor introduce nuevos retos para la conservación. Una plaga híbrida de salamandras tigre de California (STC) (Ambystoma californiense), especie en peligro, y salamandras tigre barradas (STB) (Ambystoma mavortium) introducidas ha reemplazado a las poblaciones nativas en una región amplia del centro de California, lo que amenaza la integridad genética de las STC y el sistema de estanques vernales que habitan. Realizamos un experimento de campo a gran escala y con información genética para probar si la reducción del hidroperiodo reproductivo del estanque favorecería al genotipo de las STC nativas. Construimos 14 estanques seminaturales grandes para analizar el efecto de la duración del hidroperiodo sobre la supervivencia y masa larval durante la metamorfosis. Monitoreamos los cambios en la frecuencia de alelos no nativos con una matriz de captura de exones de 5,237 genes y utilizamos una combinación de modelos lineales generalizados y bayesianos a medida para cuantificar los efectos de la duración del estanque sobre la adaptabilidad de las salamandras. Los primeros trabajos en este sistema mostraron la superioridad híbrida bajo varias condiciones y sugirieron que los híbridos están favorecidos en los estanques con modificaciones antropogénicas e hidroperiodos de larga duración artificial. En coherencia con estos primeros resultados, encontramos evidencia abrumadora de la superioridad híbrida. Los hidroperiodos muy cortos redujeron sustancialmente la masa (1.1­1.5 más veces) y la probabilidad de supervivencia (10­13 más veces) de las larvas nativas e híbridas, lo que confirma que el hidroperiodo probablemente ejerce una fuerte presión selectiva en vida silvestre. Identificamos 86 genes, que representan el 1.8% de los 4,723 loci examinados, que respondieron significativamente a la selección basada en el hidroperiodo. Con las condiciones más naturales de nuestro experimento, y en contraste a nuestros primeros trabajos, la supervivencia y el tamaño de las STC nativas durante la metamorfosis siempre fueron menores a las de los híbridos, lo que sugiere que el manejo del hidroperiodo por sí solo no cambiará la selección a favor de los genotipos larvales nativos. Sin embargo, la reducción del hidroperiodo del estanque puede limitar la productividad de los estanques híbridos y complementar otras estrategias para extirpar a los híbridos mientras que mantiene el ecosistema del estanque vernal. Este estudio confirma y amplía los trabajos anteriores que resaltan la importancia del manejo del hidroperiodo para controlar las especies acuáticas invasoras.

The reference genome of the Vernal Pool Tadpole Shrimp, Lepidurus packardi.

In this paper, we report on the scaffold-level assembled genome for the federally endangered, California endemic crustacean Lepidurus packardi (the Vernal Pool Tadpole Shrimp). L. packardi is a key food source for other conserved California species including the California Tiger Salamander Ambystoma californiense. It faces significant habitat loss and fragmentation as vernal pools are threatened by urbanization, agricultural conversion, and climate change. This resource represents the first scaffold-level genome of any Lepidurus species. The assembled genome spans 108.6 Mbps, with 6 chromosome-length scaffolds comprising 71% of total genomic length and 444 total contigs. The BUSCO score for this genome is 97.3%, suggesting a high level of completeness. We produced a predicted gene set for this species trained on the Daphnia magna set of genes and predicted 17,650 genes. These tools can aid researchers in understanding the evolution and adaptive potential of alternative reproductive modes within this species.

Patterns of bacterial diversity in embryonic capsules of the spotted salamander Ambystoma maculatum: an expanding view of a symbiosis.

The unicellular green alga, Oophila amblystomatis, populates egg capsules of the spotted salamander Ambystoma maculatum. This nutrient-exchange mutualism is widely perceived as a bipartite interaction, but the presence and contributing effects of bacteria to this symbiosis are unknown. We used standard cultivation techniques and amplicon sequencing of the V4/V5 region of 16S rRNA gene to identify and compare diversity of bacterial taxa in embryonic capsules with that in the aquatic breeding habitat. Our sampling regime allowed us to investigate diversity among individual capsules of an egg mass and between two ponds and sampling years. Capsules contain much lower diversity of bacteria than pond water, and spatial and temporal variation in intracapsular and pond bacterial diversity was observed. Despite this variation, sequences corresponding to species in the orders Burkholderiales and Oligoflexales were either prevalent or abundant, or both. Isolates most commonly recovered from capsules were closely related to species in the genus Herbaspirillum (Burkholderiaceae); other isolates were pseudomonads, but in all cases are closely related to known vascular plant-associated species. We conclude that, despite observed variation, there are bacterial taxa whose presence is held in common spatially and temporally among capsules and that the symbiosis between O. amblystomatis and A. maculatum may involve these taxa.

Geography is more important than life history in the recent diversification of the tiger salamander complex.

The North American tiger salamander species complex, including its best-known species, the Mexican axolotl, has long been a source of biological fascination. The complex exhibits a wide range of variation in developmental life history strategies, including populations and individuals that undergo metamorphosis; those able to forego metamorphosis and retain a larval, aquatic lifestyle (i.e., paedomorphosis); and those that do both. The evolution of a paedomorphic life history state is thought to lead to increased population genetic differentiation and ultimately reproductive isolation and speciation, but the degree to which it has shaped population- and species-level divergence is poorly understood. Using a large multilocus dataset from hundreds of samples across North America, we identified genetic clusters across the geographic range of the tiger salamander complex. These clusters often contain a mixture of paedomorphic and metamorphic taxa, indicating that geographic isolation has played a larger role in lineage divergence than paedomorphosis in this system. This conclusion is bolstered by geography-informed analyses indicating no effect of life history strategy on population genetic differentiation and by model-based population genetic analyses demonstrating gene flow between adjacent metamorphic and paedomorphic populations. This fine-scale genetic perspective on life history variation establishes a framework for understanding how plasticity, local adaptation, and gene flow contribute to lineage divergence. Many members of the tiger salamander complex are endangered, and the Mexican axolotl is an important model system in regenerative and biomedical research. Our results chart a course for more informed use of these taxa in experimental, ecological, and conservation research.

Allele-specific expression and gene regulation help explain transgressive thermal tolerance in non-native hybrids of the endangered California tiger salamander (Ambystoma californiense).

Hybridization between native and non-native species is an ongoing global conservation threat. Hybrids that exhibit traits and tolerances that surpass parental values are of particular concern, given their potential to outperform native species. Effective management of hybrid populations requires an understanding of both physiological performance and the underlying mechanisms that drive transgressive hybrid traits. Here, we explore several aspects of the hybridization between the endangered California tiger salamander (Ambystoma californiense; CTS) and the introduced barred tiger salamander (Ambystoma mavortium; BTS). We assayed critical thermal maximum (CTMax) to compare the ability of CTS, BTS and F1 hybrids to tolerate acute thermal stress, and found that hybrids exhibit a wide range of CTMax values, with 33% (4/12) able to tolerate temperatures greater than either parent. We then quantified the genomic response, measured at the RNA transcript level, of each salamander, to explore the mechanisms underlying thermal tolerance strategies. We found that CTS and BTS have strikingly different values and tissue-specific patterns of overall gene expression, with hybrids expressing intermediate values. F1 hybrids display abundant and variable degrees of allele-specific expression (ASE), likely arising from extensive compensatory evolution in gene regulatory mechanisms between CTS and BTS. We found evidence that the proportion of genes with allelic imbalance in individual hybrids correlates with their CTMax, suggesting a link between ASE and expanded thermal tolerance that may contribute to the success of hybrid salamanders in California. Future climate change may further complicate management of CTS if hybrid salamanders are better equipped to deal with rising temperatures.

Depauperate major histocompatibility complex variation in the endangered reticulated flatwoods salamander (Ambystoma bishopi).

Reticulated flatwoods salamander (Ambystoma bishopi) populations began decreasing dramatically in the 1900s. Contemporary populations are small, isolated, and may be susceptible to inbreeding and reduced adaptive potential because of low genetic variation. Genetic variation at immune genes is especially important as it influences disease susceptibility and adaptation to emerging infectious pathogens, a central conservation concern for declining amphibians. We collected samples from across the extant range of this salamander to examine genetic variation at major histocompatibility complex (MHC) class Iα and IIß exons as well as the mitochondrial control region. We screened tail or toe tissue for ranavirus, a pathogen associated with amphibian declines worldwide. Overall, we found low MHC variation when compared to other amphibian species and did not detect ranavirus at any site. MHC class Iα sequencing revealed only three alleles with a nucleotide diversity of 0.001, while MHC class IIß had five alleles with a with nucleotide diversity of 0.004. However, unique variation still exists across this species' range with private alleles at three sites. Unlike MHC diversity, mitochondrial variation was comparable to levels estimated for other amphibians with nine haplotypes observed, including one haplotype shared across all sites. We hypothesize that a combination of a historic disease outbreak and a population bottleneck may have contributed to low MHC diversity while maintaining higher levels of mitochondrial DNA variation. Ultimately, MHC data indicated that the reticulated flatwoods salamander may be at an elevated risk from infectious diseases due to low levels of immunogenetic variation necessary to combat novel pathogens.

Multiscale resistant kernel surfaces derived from inferred gene flow: An application with vernal pool breeding salamanders.

The importance of assessing spatial data at multiple scales when modelling species-environment relationships has been highlighted by several empirical studies. However, no landscape genetics studies have optimized landscape resistance surfaces by evaluating relevant spatial predictors at multiple spatial scales. Here, we model multiscale/layer landscape resistance surfaces to estimate resistance to inferred gene flow for two vernal pool breeding salamander species, spotted (Ambystoma maculatum) and marbled (A. opacum) salamanders. Multiscale resistance surface models outperformed spatial layers modelled at their original spatial scale. A resistance surface with forest land cover at a 500-m Gaussian kernel bandwidth and normalized vegetation index at a 100-m Gaussian kernel bandwidth was the top optimized resistance surface for A. maculatum, while a resistance surface with traffic rate and topographic curvature, both at a 500-m Gaussian kernel bandwidth, was the top optimized resistance surface for A. opacum. Species-specific resistant kernels were fit at all vernal pools in our study area with the optimized multiscale/layer resistance surface controlling kernel spread. Vernal pools were then evaluated and scored based on surrounding upland habitat (local score) and connectivity with other vernal pools on the landscape, with resistant kernels driving vernal pool connectivity scores. As expected, vernal pools that scored highest were in areas within forested habitats and with high vernal pool densities and low species-specific landscape resistance. Our findings highlight the success of using a novel analytical approach in a multiscale framework with applications beyond vernal pool amphibian conservation.

Genomic data recover previously undetectable fragmentation effects in an endangered amphibian.

A critical consideration when using molecular ecological methods to detect trends and parameterize models at very fine spatial and temporal scales has always been the technical limits of resolution. Key landscape features, including most anthropogenic modifications, can cause biologically important, but very recent changes in gene flow that require substantial statistical power to detect. The problem is one of temporal scale: Human change is rapid and recent, while genetic changes accumulate slowly. We generated SNPs from thousands of nuclear loci to characterize the population structure of New York-endangered eastern tiger salamanders (Ambystoma tigrinum) on Long Island and quantify the impacts of roads on population fragmentation. In stark contrast to a recent microsatellite study, we uncovered highly structured populations over an extremely small spatial scale (approximately 40 km2 ) in an increasingly human-modified landscape. Geographic distance and the presence of roads between ponds were both strong predictors of genetic divergence, suggesting that both natural and anthropogenic factors contribute to the observed patterns of genetic variation. All ponds supported small to modest effective breeding populations, and pond surface area showed a strong positive correlation with population size. None of these patterns emerged in an earlier study of the same system using microsatellite loci, and we determined that at least 300-400 SNPs were needed to recover the fine-scale population structure present in this system. Conservation assessments using earlier genetic techniques in other species may similarly lack the statistical power for small-scale inferences and benefit from reassessments using genomic tools.

Influence of genome and bio-ecology on the prevalence of genome exchange in unisexuals of the Ambystoma complex.

BACKGROUND: Unisexuals of the blue-spotted salamander complex are thought to reproduce by kleptogenesis. Genome exchanges associated with this sperm-dependent mode of reproduction are expected to result in a higher genetic variation and multiple ploidy levels compared to clonality. However, the existence of some populations exclusively formed of genetically identical individuals suggests that factors could prevent genome exchanges. This study aimed at assessing the prevalence of genome exchange among unisexuals of the Ambystoma laterale-jeffersonianum complex from 10 sites in the northern part of their distribution. RESULTS: A total of 235 individuals, including 207 unisexuals, were genotyped using microsatellite loci and AFLP. Unisexual individuals could be sorted in five genetically distinct groups, likely derived from the same paternal A. jeffersonianum haplome. One of these groups exclusively reproduced clonally, even when found in sympatry with lineages presenting signature of genome exchange. Genome exchange was site-dependent for another group. Genome exchange was detected at all sites for the three remaining groups. CONCLUSION: Prevalence of genome exchange appears to be associated with ecological conditions such as availability of effective sperm donors. Intrinsic genomic factors may also affect this process, since different lineages in sympatry present highly variable rate of genome exchange. The coexistence of clonal and genetically diversified lineages opens the door to further research on alternatives to genetic variation.

Natural tissue concentrations in adult Ambystoma maculatum and larval DNA damage from exposure to arsenic and chromium.

Arsenic (As) and chromium (Cr) are two contaminants that are detected in aquatic and terrestrial habitats. Using the spotted salamander, Ambystoma maculatum, to assess impacts from these contaminants may be advantageous as adults live and breed in such environments. Adult amphibians typically exhibit elevated tissue concentrations of contaminants present in their environment, while larval stages were found to exhibit increased sensitivity to pollutants. From January through March of 2015, during the spring breeding season, 5 adults and approximately 32 egg masses were collected from a local breeding site. Field levels of As and Cr ranged from 5.99 to 8.88 µg/L and 1.45 to 2 µg/L, respectively, while mean adult As tissue concentrations were 56.74 µg/g dry weight for heart, 0.92 µg/g for liver, and 1.21 µg/g for tail tissue. Mean tissue concentrations for Cr were 87.64 µg/g for heart, 1.47 µg/g for liver, and 6.92 µg/g for tail. Developing larvae that were collected from the field and exposed in a lab setting for 12 d to 0.2 or 20 mg/L of either As or Cr displayed little DNA damage attributed to As, but marked damage due to exposure to 20 mg/L Cr when assessed using the comet assay. Exposure to a mixture of either 0.25:0.1 or 25:10 mg/L As and Cr resulted in significant DNA damage at the lower concentration of 0.25:0.1 mg/L. As adult spotted salamanders were found to possess high concentrations of these contaminants in cardiac tissue, and larvae were shown to be susceptible to DNA damage from increased exposures, assessing impacts and potential declines of amphibian populations exposed to As and Cr is needed.

Transcriptome analysis illuminates the nature of the intracellular interaction in a vertebrate-algal symbiosis.

During embryonic development, cells of the green alga Oophila amblystomatis enter cells of the salamander Ambystoma maculatum forming an endosymbiosis. Here, using de novo dual-RNA seq, we compared the host salamander cells that harbored intracellular algae to those without algae and the algae inside the animal cells to those in the egg capsule. This two-by-two-way analysis revealed that intracellular algae exhibit hallmarks of cellular stress and undergo a striking metabolic shift from oxidative metabolism to fermentation. Culturing experiments with the alga showed that host glutamine may be utilized by the algal endosymbiont as a primary nitrogen source. Transcriptional changes in salamander cells suggest an innate immune response to the alga, with potential attenuation of NF-κB, and metabolic alterations indicative of modulation of insulin sensitivity. In stark contrast to its algal endosymbiont, the salamander cells did not exhibit major stress responses, suggesting that the host cell experience is neutral or beneficial.

A linkage map for the Newt Notophthalmus viridescens: Insights in vertebrate genome and chromosome evolution.

Genetic linkage maps are fundamental resources that enable diverse genetic and genomic approaches, including quantitative trait locus (QTL) analyses and comparative studies of genome evolution. It is straightforward to build linkage maps for species that are amenable to laboratory culture and genetic crossing designs, and that have relatively small genomes and few chromosomes. It is more difficult to generate linkage maps for species that do not meet these criteria. Here, we introduce a method to rapidly build linkage maps for salamanders, which are known for their enormous genome sizes. As proof of principle, we developed a linkage map with thousands of molecular markers (N=2349) for the Eastern newt (Notophthalmus viridescens). The map contains 12 linkage groups (152.3-934.7cM), only one more than the number of chromosome pairs. Importantly, this map was generated using RNA isolated from a single wild caught female and her 28 offspring. We used the map to reveal chromosome-scale conservation of synteny among N. viridescens, A. mexicanum (Urodela), and chicken (Amniota), and to identify large conserved segments between N. viridescens and Xenopus tropicalis (Anura). We also show that met1, a major effect QTL that regulates the expression of alternate metamorphic and paedomorphic modes of development in Ambystoma, associates with a chromosomal fusion that is not found in the N. viridescens map. Our results shed new light on the ancestral amphibian karyotype and reveal specific fusion and translocation events that shaped the genomes of three amphibian model taxa. The ability to rapidly build linkage maps for large salamander genomes will enable genetic and genomic analyses within this important vertebrate group, and more generally, empower comparative studies of vertebrate biology and evolution.

An approach for identifying cryptic barriers to gene flow that limit species' geographic ranges.

Species' geographic range limits are most often not demarcated by obvious dispersal barriers. Poor-quality habitat at the edge of a species' range can prevent range expansion by preventing outward migration or through reducing adaptive potential resulting from decreased genetic diversity. We identified habitat variables that constrain gene flow across the entire geographic range of an endemic salamander (Ambystoma barbouri) in the eastern United States, and we tested whether increased resistance resulting from these variables provides cryptic dispersal barriers at the range edges. Using polymorphic microsatellite loci, we first identified three genetic clusters that are separated by the Ohio and Kentucky rivers. Through a combination of landscape genetic analyses and generalized dissimilarity modelling, we then classified variables that (i) restrict gene flow in each of the genetic clusters across the geographic distribution of A. barbouri and (ii) become more common towards the peripheries of the distribution. A decrease in limestone availability and an increase in growing season precipitation were correlated with high resistance to gene flow across the range, and both became more common at the edges of the species' distribution. However, other landscape variables were more important for explaining variation in geneflow rates in different portions of the range, such as increased mean annual temperature and frost-free period in the south vs. growing season precipitation in the north. Taken together, these results suggest that there are both range-wide and regionally specific cryptic habitat barriers preventing geographic range expansion. Species 'geographic range limits are probably governed by a set of ecological and evolutionary factors, and our landscape genetic approach could be applied to gain additional insight into many systems.

Cryptic sex? Estimates of genome exchange in unisexual mole salamanders (Ambystoma sp.).

Cryptic sex has been argued to explain the exceptional longevity of certain parthenogenetic vertebrate lineages, yet direct measurements of genetic exchange between sexual and apparently parthenogenetic forms are rare. Female unisexual mole salamanders (Ambystoma sp.) are the oldest known unisexual vertebrate lineage (~5 million years), and one hypothesis for their persistence is that allopolyploid female unisexuals periodically exchange haploid genomes 'genome exchange' during gynogenetic reproduction with males from sympatric sexual species. We test this hypothesis by using genome-specific microsatellite DNA markers to estimate the rates of genome exchange between sexual males and unisexual females in two ponds in NE Ohio. We also test the prediction that levels of gene flow should be higher for 'sympatric' (sexual males present) genomes in unisexuals compared to 'allopatric' (sexual males absent) unisexual genomes. We used a model testing framework in the coalescent-based program MIGRATE-N to compare models where unidirectional gene flow is present and absent between sexual species and unisexuals. As predicted, our results show higher levels of gene flow between sexuals and sympatric unisexual genomes compared to lower (likely artefactual) levels of gene flow between sexuals and allopatric unisexual genomes. Our results provide direct evidence that genome exchange between sexual and unisexual Ambystoma occurs and demonstrate that the magnitude depends on which sexual species are present. The relatively high levels of gene flow suggest that unisexuals must be at a selective advantage over sexual forms so as to avoid extinction due to genetic swamping through genome exchange.

Asymmetric Introgression in a Spotted Salamander Hybrid Zone.

Before the establishment of reproductive isolation, deeply diverged intraspecific lineages can experience complex genetic and behavioral interactions as they come into secondary contact. Divergent selective and demographic processes mediate gene flow among lineages, resulting in hybrid zones with complex biogeographic structure. Discordance in the biogeographic patterns of autosomal and maternally inherited loci provides a useful window to infer the processes mediating admixture and introgression across hybrid zones. Here, we sampled 489 genotypes across a hybrid zone between 2 phylogeographic lineages of the spotted salamander, Ambystoma maculatum, and characterize discordant patterns of nuclear and mitochondrial introgression across the contact boundary. Our results indicate asymmetric introgression of nuclear DNA beyond the contact boundary from the western to eastern lineage, with introgression of eastern mitochondrial DNA into the western lineage. We discuss alternative mechanisms for this pattern and attribute this result to neutral patterns of population expansion of the western lineage into the east in combination with female mate choice for larger-bodied western males. Our results underscore the complexity of interacting mechanisms that give rise to reproductive asymmetries in the earliest stages of the speciation process.

The importance (or lack thereof) of niche divergence to the maintenance of a northern species complex: the case of the long-toed salamander (Ambystoma macrodactylum Baird).

The relative importance of ecological vs. non-ecological factors for the origin and maintenance of species is an open question in evolutionary biology. Young lineages--such as the distinct genetic groups that make up the ranges of many northern species--represent an opportunity to study the importance of ecological divergence during the early stages of diversification. Yet, few studies have examined the extent of niche divergence between lineages in previously glaciated regions and the role of ecology in maintaining the contact zones between them. In this study, we used tests of niche overlap in combination with ecological niche models to explore the extent of niche divergence between lineages of the long-toed salamander (Ambystoma macrodactylum Baird) species complex and to determine whether contact zones correspond to (divergent) niche limits. We found limited evidence for niche divergence between the different long-toed salamander lineages, substantial overlap in the predicted distribution of suitable climatic space for all lineages and range limits that are independent of niche limits. These results raise questions as to the importance of ecological divergence to the development of this widespread species complex and highlight the potential for non-ecological factors to play a more important role in the maintenance of northern taxa.

Limited influence of local and landscape factors on finescale gene flow in two pond-breeding amphibians.

Dispersal and gene flow within animal populations are influenced by the composition and configuration of the landscape. In this study, we evaluated hypotheses about the impact of natural and anthropogenic factors on genetic differentiation in two amphibian species, the spotted salamander (Ambystoma maculatum) and the wood frog (Lithobates sylvaticus) in a commercial forest in central Maine. We conducted this analysis at two scales: a local level, focused on factors measured at each breeding pond, and a landscape level, focused on factors measured between ponds. We investigated the effects of a number of environmental factors in six categories including Productivity, Physical, Land Composition, Land Configuration, Isolation and Location. Embryos were sampled from 56 spotted salamander breeding ponds and 39 wood frog breeding ponds. We used a hierarchical Bayesian approach in the program GESTE at each breeding pond and a random forest algorithm in conjunction with a network analysis between the ponds. We found overall high genetic connectivity across distances up to 17 km for both species and a limited effect of natural and anthropogenic factors on gene flow. We found the null models best explained patterns of genetic differentiation at a local level and found several factors at the landscape level that weakly influenced gene flow. This research indicates multiscale investigations that incorporate local and landscape factors are valuable for understanding patterns of gene flow. Our findings suggest that dispersal rates in this system are high enough to minimize genetic structuring and that current forestry practices do not significantly impede dispersal.

A test of the central-marginal hypothesis using population genetics and ecological niche modelling in an endemic salamander (Ambystoma barbouri).

The central-marginal hypothesis (CMH) predicts that population size, genetic diversity and genetic connectivity are highest at the core and decrease near the edges of species' geographic distributions. We provide a test of the CMH using three replicated core-to-edge transects that encompass nearly the entire geographic range of the endemic streamside salamander (Ambystoma barbouri). We confirmed that the mapped core of the distribution was the most suitable habitat using ecological niche modelling (ENM) and via genetic estimates of effective population sizes. As predicted by the CMH, we found statistical support for decreased genetic diversity, effective population size and genetic connectivity from core to edge in western and northern transects, yet not along a southern transect. Based on our niche model, habitat suitability is lower towards the southern range edge, presumably leading to conflicting core-to-edge genetic patterns. These results suggest that multiple processes may influence a species' distribution based on the heterogeneity of habitat across a species' range and that replicated sampling may be needed to accurately test the CMH. Our work also emphasizes the importance of identifying the geographic range core with methods other than using the Euclidean centre on a map, which may help to explain discrepancies among other empirical tests of the CMH. Assessing core-to-edge population genetic patterns across an entire species' range accompanied with ENM can inform our general understanding of the mechanisms leading to species' geographic range limits.