Geology topics

Kevin P. Mulder

Publications and source records attributed to Kevin P. Mulder.

4 recordsLinked to original sources

USGS2022

Speciation with gene flow in a narrow endemic West Virginia cave salamander (Gyrinophilus subterraneus)

Due to their limited geographic distributions and specialized ecologies, cave species are often highly endemic and can be especially vulnerable to habitat degradation within and surrounding the cave systems they inhabit. We investigated the evolutionary history of the West Virginia Spring Salamander ( Gyrinophilus subterraneus) , estimated the population trend from historic and current survey data, and assessed the current potential for water quality threats to the cave habitat. Our genomic data (mtDNA sequence and ddRADseq-derived SNPs) reveal two, distinct evolutionary lineages within General Davis Cave corresponding to G. subterraneus and its widely distributed sister species, Gyrinophilus porphyriticus , that are also differentiable based on morphological traits. Genomic models of evolutionary history strongly support asymmetric and continuous gene flow between the two lineages, and hybrid classification analyses identify only parental and first generation cross (F1) progeny. Collectively, these results point to a rare case of sympatric speciation occurring within the cave, leading to strong support for continuing to recognize G. subterraneus as a distinct and unique species. Due to its specialized habitat requirements, the complete distribution of G. subterraneus is unresolved, but using survey data in its type locality (and currently the only known occupied site), we find that the population within General Davis Cave has possibly declined over the last 45 years. Finally, our measures of cave and surface stream water quality did not reveal evidence of water quality impairment and provide important baselines for future monitoring. In addition, our unexpected finding of a hybrid zone and partial reproductive isolation between G. subterraneus and G. porphyriticus warrants further attention to better understand the evolutionary and conservation implications of occasional hybridization between the species.

West Virginia↗

USGS2019

Functional variation at an expressed MHC class IIß locus associates with Ranavirus infection intensity in larval anuran populations

Infectious diseases are causing catastrophic losses to biodiversity globally. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian susceptibility to pathogens. MHC genes encode diverse cell-surface molecules that can recognize and bind to a wide array of pathogen peptides, and are divided into two classes. While MHC class I genes are the classic mediators of viral acquired immunity, larval amphibians do not express them. However, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Lithobates sylvaticus), which are highly susceptible to Ranavirus, across 17 ponds and two years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIß locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if physiochemical properties of peptide-binding regions influenced host responses to Ranavirus. Among 334 sampled individuals, 26% were infected with Ranavirus, and among infected individuals the average intensity was 7.12 x 107 virus copies. We recovered 20 unique MHC class IIß alleles that fell into two deeply diverged clades and seven functional supertypes. Variation in MHC supertypes were associated with Ranavirus infection intensity, but not prevalence. MHC supertype heterozygotes and individuals with the MHC supertype genotype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and all other genotypes. We conclude that MHC class IIß functional genetic variation is an important component of Ranavirus susceptibility. Identifying immune system gene signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

Immunogenetics↗

USGS2018

Two-species occupancy modeling accounting for species misidentification and nondetection

In occupancy studies, species misidentification can lead to false‐positive detections, which can cause severe estimator biases. Currently, all models that account for false‐positive errors only consider omnibus sources of false detections and are limited to single‐species occupancy. However, false detections for a given species often occur because of the misidentification with another, closely related species. To exploit this explicit source of false‐positive detection error, we develop a two‐species occupancy model that accounts for misidentifications between two species of interest. As with other false‐positive models, identifiability is greatly improved by the availability of unambiguous detections at a subset of site x occasions. Here, we consider the case where some of the field observations can be confirmed using laboratory or other independent identification methods (“confirmatory data”). We performed three simulation studies to (1) assess the model's performance under various realistic scenarios, (2) investigate the influence of the proportion of confirmatory data on estimator accuracy and (3) compare the performance of this two‐species model with that of the single‐species false‐positive model. The model shows good performance under all scenarios, even when only small proportions of detections are confirmed (e.g. 5%). It also clearly outperforms the single‐species model. We illustrate application of this model using a 4‐year dataset on two sympatric species of lungless salamanders: the US federally endangered Shenandoah salamander Plethodon shenandoah , and its presumed competitor, the red‐backed salamander Plethodon cinereus . Occupancy of red‐backed salamanders appeared very stable across the 4 years of study, whereas the Shenandoah salamander displayed substantial turnover in occupancy of forest habitats among years. Given the extent of species misidentification issues in occupancy studies, this modelling approach should help improve the reliability of estimates of species distribution, which is the goal of many studies and monitoring programmes. Further developments, to account for different forms of state uncertainty, can be readily undertaken under our general approach.

Methods in Ecology and Evolution↗

USGS2017

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Developmental and Comparative Immunology↗