USGS Science⌕ Search

SEARCH · USGS Science

Results for “Genomics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19Linked to original sources
USGS2019

Putative mitochondrial sex determination in the Bivalvia: Insights from a hybrid transcriptome assembly in freshwater mussels

Bivalves exhibit an astonishing diversity of sexual systems, with genetic and environmental determinants of sex, and possibly the only example of mitochondrial genes influencing sex determination pathways in animals. In contrast to all other animal species in which strict maternal inheritance (SMI) of mitochondria is the rule, bivalves possess a system known as doubly uniparental inheritance (DUI) of mitochondria in which maternal and paternal mitochondria (and their corresponding female-transmitted or F mtDNA and male-transmitted or M mtDNA genomes) are transmitted within a species. Species with DUI also possess sex-associated mtDNA-encoded proteins (in addition to the typical set of 13), which have been hypothesized to play a role in sex determination. In this study, we analyzed the sex-biased transcriptome in gonads of two closely-related freshwater mussel species with different reproductive and mitochondrial transmission modes: the gonochoric, DUI species, Utterbackia peninsularis , and the hermaphroditic, SMI species, Utterbackia imbecillis . Through comparative analysis with other DUI and non-DUI bivalve transcriptomes already available, we identify common male and female-specific genes, as well as SMI and DUI-related genes, that are probably involved in sex determination and mitochondrial inheritance in this animal group. Our results contribute to the understanding of what could be the first animal sex determination system involving the mitochondrial genome

Florida↗
USGS2022

As the goose flies: Migration routes and timing influence patterns of genetic diversity in a circumpolar migratory herbivore

Migration schedules and the timing of other annual events (e.g., pair formation and molt) can affect the distribution of genetic diversity as much as where these events occur. The greater white-fronted goose ( Anser albifrons ) is a circumpolar goose species, exhibiting temporal and spatial variation of events among populations during the annual cycle. Previous range-wide genetic assessments of the nuclear genome based on eight microsatellite loci suggest a single, largely panmictic population despite up to five subspecies currently recognized based on phenotypic differences. We used double digest restriction-site associated DNA (ddRAD-seq) and mitochondrial DNA (mtDNA) sequence data to re-evaluate estimates of spatial genomic structure and to characterize how past and present processes have shaped the patterns of genetic diversity and connectivity across the Arctic and subarctic. We uncovered previously undetected inter-population differentiation with genetic clusters corresponding to sampling locales associated with current management groups. We further observed subtle genetic clustering within each management unit that can be at least partially explained by the timing and directionality of migration events along with other behaviors during the annual cycle. The Tule Goose ( A. a. elgasi ) and Greenland subspecies ( A. a. flavirostris ) showed the highest level of divergence among all sampling locales investigated. The recovery of previously undetected broad and fine-scale spatial structure suggests that the strong cultural transmission of migratory behavior restricts gene flow across portions of the species’ range. Our data further highlight the importance of re-evaluating previous assessments conducted based on a small number of highly variable genetic markers in phenotypically diverse species.

Circumpolar Arctic↗
USGS2020

Novel circoviruses detected in feces of Sonoran felids

Sonoran felids are threatened by drought and habitat fragmentation. Vector range expansion and anthropogenic factors such as habitat encroachment and climate change are altering viral evolutionary dynamics and exposure. However, little is known about the diversity of viruses present in these populations. Small felid populations with lower genetic diversity are likely to be most threatened with extinction by emerging diseases, as with other selective pressures, due to having less adaptive potential. We used a metagenomic approach to identify novel circoviruses, which may have a negative impact on the population viability, from confirmed bobcat ( Lynx rufus ) and puma ( Puma concolor ) scats collected in Sonora, Mexico. Given some circoviruses are known to cause disease in their hosts, such as porcine and avian circoviruses, we took a non-invasive approach using scat to identify circoviruses in free-roaming bobcats and puma. Three circovirus genomes were determined, and, based on the current species demarcation, they represent two novel species. Phylogenetic analyses reveal that one circovirus species is more closely related to rodent associated circoviruses and the other to bat associated circoviruses, sharing highest genome-wide pairwise identity of approximately 70% and 63%, respectively. At this time, it is unknown whether these scat-derived circoviruses infect felids, their prey, or another organism that might have had contact with the scat in the environment. Further studies should be conducted to elucidate the host of these viruses and assess health impacts in felids

Sonoran Desert↗
USGS2017

Bat white-nose syndrome: a real-time TaqMan polymerase chain reaction test targeting the intergenic spacer region of Geomyces destructans tructans.

The fungus Geomyces destructans is the causative agent of white-nose syndrome (WNS), a disease that has killed millions of North American hibernating bats. We describe a real-time TaqMan PCR test that detects DNA from G. destructans by targeting a portion of the multicopy intergenic spacer region of the rRNA gene complex. The test is highly sensitive, consistently detecting as little as 3.3 fg of genomic DNA from G. destructans . The real-time PCR test specifically amplified genomic DNA from G. destructans but did not amplify target sequence from 54 closely related fungal isolates (including 43 Geomyces spp. isolates) associated with bats. The test was further qualified by analyzing DNA extracted from 91 bat wing skin samples, and PCR results matched histopathology findings. These data indicate the real-time TaqMan PCR method described herein is a sensitive, specific, and rapid test to detect DNA from G. destructans and provides a valuable tool for WNS diagnostics and research.

Mycologia↗
USGS2021

Historical effective population size of North American hoary bat (Lasiurus cinereus) and challenges to estimating trends in contemporary effective breeding population size from archived samples

Background Hoary bats ( Lasiurus cinereus ) are among the bat species most commonly killed by wind turbine strikes in the midwestern United States. The impact of this mortality on species census size is not understood, due in part to the difficulty of estimating population size for this highly migratory and elusive species. Genetic effective population size (Ne) could provide an index of changing census population size if other factors affecting Ne are stable. Methods We used the NeEstimator package to derive effective breeding population size (Nb) estimates for two temporally spaced cohorts: 93 hoary bats collected in 2009–2010 and an additional 93 collected in 2017–2018. We sequenced restriction-site associated polymorphisms and generated a de novo genome assembly to guide the removal of sex-linked and multi-copy loci, as well as identify physically linked markers. Results Analysis of the reference genome with psmc suggested at least a doubling of Ne in the last 100,000 years, likely exceeding Ne = 10,000 in the Holocene. Allele and genotype frequency analyses confirmed that the two cohorts were comparable, although some samples had unusually high or low observed heterozygosities. Additionally, the older cohort had lower mean coverage and greater variability in coverage, and batch effects of sampling locality were observed that were consistent with sample degradation. We therefore excluded samples with low coverage or outlier heterozygosity, as well as loci with sequence coverage far from the mode value, from the final data set. Prior to excluding these outliers, contemporary Nb estimates were significantly higher in the more recent cohort, but this finding was driven by high values for the 2018 sample year and low values for all other years. In the reduced data set, Nb did not differ significantly between cohorts. We found base substitutions to be strongly biased toward cytosine to thymine or the complement, and further partitioning loci by substitution type had a strong effect on Nb estimates. Minor allele frequency and base quality bias thresholds also had strong effects on Nb estimates. Instability of Nb with respect to common data filtering parameters and empirically identified factors prevented robust comparison of the two cohorts. Given that confidence intervals frequently included infinity as the stringency of data filtering increased, contemporary trends in Nb of North American hoary bats may not be tractable with the linkage disequilibrium method, at least using the protocol employed here.

PeerJ↗
USGS2022

Positively selected genes in the hoary bat (Lasiurus cinereus) lineage: Prominence of thymus expression, immune and metabolic function, and regions of ancient synteny

Background Bats of the genus Lasiurus occur throughout the Americas and have diversified into at least 20 species among three subgenera. The hoary bat ( Lasiurus cinereus ) is highly migratory and ranges farther across North America than any other wild mammal. Despite the ecological importance of this species as a major insect predator, and the particular susceptibility of lasiurine bats to wind turbine strikes, our understanding of hoary bat ecology, physiology, and behavior remains poor. Methods To better understand adaptive evolution in this lineage, we used whole-genome sequencing to identify protein-coding sequence and explore signatures of positive selection. Gene models were predicted with Maker and compared to seven well-annotated and phylogenetically representative species. Evolutionary rate analysis was performed with PAML. Results Of 9,447 single-copy orthologous groups that met evaluation criteria, 150 genes had a significant excess of nonsynonymous substitutions along the L. cinereus branch ( P < 0.001 after manual review of alignments). Selected genes as a group had biased expression, most strongly in thymus tissue. We identified 23 selected genes with reported immune functions as well as a divergent paralog of Steep1 within suborder Yangochiroptera. Seventeen genes had roles in lipid and glucose metabolic pathways, partially overlapping with 15 mitochondrion-associated genes; these adaptations may reflect the metabolic challenges of hibernation, long-distance migration, and seasonal variation in prey abundance. The genomic distribution of positively selected genes differed significantly from background expectation by discrete Kolmogorov–Smirnov test ( P < 0.001). Remarkably, the top three physical clusters all coincided with islands of conserved synteny predating Mammalia, the largest of which shares synteny with the human cat-eye critical region (CECR) on 22q11. This observation coupled with the expansion of a novel Tbx1 -like gene family may indicate evolutionary innovation during pharyngeal arch development: both the CECR and Tbx1 cause dosage-dependent congenital abnormalities in thymus, heart, and head, and craniodysmorphy is associated with human orthologs of other positively selected genes as well.

PeerJ↗
USGS2011

The use of epifluorescent microscopy and quantitative polymerase chain reaction to determine the presence/absence and identification of microorganisms associated with domestic and foreign wallboard samples

Epifluorescent microscopy and quantitative polymerase chain reaction (qPCR) were utilized to determine the presence, concentration and identification of bacteria, and more specifically sulfate reducing bacteria (SRB) in subsamples of Chinese and North American wallboard, and wallboard-mine rock. Bacteria were visible in most subsamples, which included wallboard-lining paper from each side of the wallboard, wallboard filler, wallboard tape and fragments of mined wallboard rock via microscopy. Observed bacteria occurred as single or small clusters of cells and no mass aggregates indicating colonization were noted. Universal 16S qPCR was utilized to directly examine samples and detected bacteria at concentrations ranging from 1.4 x 10 3 to 6.4 x 10 4 genomic equivalents per mm 2 of paper or per gram of wallboard filler or mined rock, in 12 of 41 subsamples. Subsamples were incubated in sulfate reducing broth for ~30 to 60 days (enrichment assay) and then analyzed by universal 16S and SRB qPCR. Enrichment universal 16S qPCR detected bacteria in 32 of 41 subsamples at concentrations ranging from 1.5 x 10 4 to 4.2 x 10 7 genomic equivalents per ml of culture broth. Evaluation of enriched subsamples by SRB qPCR demonstrated that SRB were not detectable in most of the samples and if they were detected, detection was not reproducible (an indication of low concentrations, if present). Enrichment universal 16S and SRB qPCR demonstrated that viable bacteria were present in subsamples (as expected given exposure of the samples following manufacture, transport and use) but that SRB were either not present or present at very low numbers. Further, no differences in trends were noted between the various Chinese and North American wallboard samples. In all, the microscopy and qPCR data indicated that the suspected ‘sulfur emissions’ emanating from suspect wallboard samples is not due to microbial activity.

Report↗
USGS2013

Biogeographic perspective of speciation among desert tortoises in the genus Gopherus : a preliminary evaluation

The enduring processes of time, climate, and adaptation have sculpted the distribution of organisms we observe in the Sonoran Desert. One such organism is Morafka’s desert tortoise, Gopherus morafkai. We apply a genomic approach to identify the evolutionary processes driving diversity in this species and present preliminary findings and emerging hypotheses. The Sonoran Desert form of the tortoise exhibits a continuum of genetic similarity spanning 850 km of Sonoran desertscrub extending from Empalme, Sonora, to Kingman, Arizona. However, at the ecotone between desertscrub and foothills thornscrub we identify a distinct, Sinaloan lineage and this occurrence suggests a more complex evolutionary history for G. morafkai. By using multiple loci from throughout the tortoise’s genome, we aim to determine if divergence between these lineages occurred in allopatry, and further to investigate for signatures of past or current genetic introgression. This international, collaborative project will assist state and federal agencies in developing management strategies that best preserve the evolutionary potential of Morafka’s desert tortoise. Ultimately, an understanding of the evolutionary history of desert tortoises will not only clarify the forces that have driven the divergence in this group, but also contribute to our knowledge of the biogeographic history of the Southwestern deserts and how diversity is maintained within them.

Proceedings↗
USGS2006

Conservation of Toll-like receptor signaling pathways in teleost fish

In mammals, toll-like receptors (TLR) recognize ligands, including pathogen-associated molecular patterns (PAMPs), and respond with ligand-specific induction of genes. In this study, we establish evolutionary conservation in teleost fish of key components of the TLR-signaling pathway that act as switches for differential gene induction, including MYD88, TIRAP, TRIF, TRAF6, IRF3, and IRF7. We further explore this conservation with a molecular phylogenetic analysis of MYD88. To the extent that current genomic analysis can establish, each vertebrate has one ortholog to each of these genes. For molecular tree construction and phylogeny inference, we demonstrate a methodology for including genes with only partial primary sequences without disrupting the topology provided by the high-confidence full-length sequences. Conservation of the TLR-signaling molecules suggests that the basic program of gene regulation by the TLR-signaling pathway is conserved across vertebrates. To test this hypothesis, leukocytes from a model fish, rainbow trout ( Oncorhynchus mykiss ), were stimulated with known mammalian TLR agonists including: diacylated and triacylated forms of lipoprotein, flagellin, two forms of LPS, synthetic double-stranded RNA, and two imidazoquinoline compounds (loxoribine and R848). Trout leukocytes responded in vitro to a number of these agonists with distinct patterns of cytokine expression that correspond to mammalian responses. Our results support the key prediction from our phylogenetic analyses that strong selective pressure of pathogenic microbes has preserved both TLR recognition and signaling functions during vertebrate evolution.

Comparative Biochemistry and Physiology, Part D: G↗
USGS2018

A transcriptome resource for the Arctic Cod (Boreogadus saida)

Arctic Cod ( Boreogadus saida ) serve as an important link in Arctic food webs and are thus considered an important species for environmental monitoring. RNA-Seq was conducted on samples from wild-collected individuals representing various age classes and tissue types to obtain as complete a transcriptome as possible on an Illumina MiSeq, which resulted in a total of 64,457 transcripts with an average length of 295 bp. We identified well-known genes that are associated with temperature change or response to pollutants. This RNA-Seq effort provides the first insight into the B. saida transcriptome, which can be a starting point for investigations identifying genes for local adaptation and genomic responses to future environmental change.

Marine Genomics↗
USGS2018

Deciphering the link between doubly uniparental inheritance of mtDNA and sex determination in bivalves: Clues from comparative transcriptomics

Bivalves exhibit an astonishing diversity of sexual systems and sex-determining mechanisms. They can be gonochoric, hermaphroditic or androgenetic, with both genetic and environmental factors known to determine or influence sex. One unique sex-determining system involving the mitochondrial genome has also been hypothesized to exist in bivalves with doubly uniparental inheritance (DUI) of mtDNA. However, the link between DUI and sex determination remains obscure. In this study, we performed a comparative gonad transcriptomics analysis for two DUI-possessing freshwater mussel species to better understand the mechanisms underlying sex determination and DUI in these bivalves. We used a BLAST reciprocal analysis to identify orthologs between Venustaconcha ellipsiformis and Utterbackia peninsularis and compared our results with previously published sex-specific bivalve transcriptomes to identify conserved sex-determining genes. We also compared our data with other DUI species to identify candidate genes possibly involved in the regulation of DUI. A total of ∼12,000 orthologous relationships were found, with 2,583 genes differentially expressed in both species. Among these genes, key sex-determining factors previously reported in vertebrates and in bivalves (e.g., Sry, Dmrt1, Foxl2 ) were identified, suggesting that some steps of the sex-determination pathway may be deeply conserved in metazoans. Our results also support the hypothesis that a modified ubiquitination mechanism could be responsible for the retention of the paternal mtDNA in male bivalves, and revealed that DNA methylation could also be involved in the regulation of DUI. Globally, our results suggest that sets of genes associated with sex determination and DUI are similar in distantly-related DUI species.

Genome Biology and Evolution↗
USGS2024

Metagenomic sequencing sheds light on microbes putatively associated with pneumonia-related fatalities of white-tailed deer (Odocoileus virginianus)

With emerging infectious disease outbreaks in human, domestic and wild animal populations on the rise, improvements in pathogen characterization and surveillance are paramount for the protection of human and animal health, as well as the conservation of ecologically and economically important wildlife. Genomics offers a range of suitable tools to meet these goals, with metagenomic sequencing facilitating the characterization of whole microbial communities associated with emerging and endemic disease outbreaks. Here, we use metagenomic sequencing in a case-control study to identify microbes in lung tissue associated with newly observed pneumonia-related fatalities in 34 white-tailed deer ( Odocoileus virginianus ) in Wisconsin, USA. We identified 20 bacterial species that occurred in more than a single individual. Of these, only Clostridium novyi was found to substantially differ (in number of detections) between case and control sample groups; however, this difference was not statistically significant. We also detected several bacterial species associated with pneumonia and/or other diseases in ruminants ( Mycoplasma ovipneumoniae , Trueperella pyogenes , Pasteurella multocida , Anaplasma phagocytophilum , Fusobacterium necrophorum ); however, these species did not substantially differ between case and control sample groups. On average, we detected a larger number of bacterial species in case samples than controls, supporting the potential role of polymicrobial infections in this system. Importantly, we did not detect DNA of viruses or fungi, suggesting that they are not significantly associated with pneumonia in this system. Together, these results highlight the utility of metagenomic sequencing for identifying disease-associated microbes. This preliminary list of microbes will help inform future research on pneumonia-associated fatalities of white-tailed deer.

Microbial Genomics↗
USGS2016

H9N2 influenza A virus isolated from a Greater White-fronted wild goose (Anser albifrons) in Alaska has a mutation in the PB2 gene, which is associated with pathogenicity in human pandemic 2009 H1N1

We report here the genomic sequence of an H9N2 influenza A virus [A/greater white-fronted goose/Alaska/81081/2008 (H9N2)]. This virus shares ≥99.8% identity with a previously reported virus. Both strains contain a G590S mutation in the polymerase basic 2 (PB2) gene, which is a pathogenicity marker in the pandemic 2009 H1N1 virus when combined with R591.

Genome Announcements↗
USGS2015

Metagenomic analysis of planktonic microbial consortia from a non-tidal urban-impacted segment of James River

Knowledge of the diversity and ecological function of the microbial consortia of James River in Virginia, USA, is essential to developing a more complete understanding of the ecology of this model river system. Metagenomic analysis of James River's planktonic microbial community was performed for the first time using an unamplified genomic library and a 16S rDNA amplicon library prepared and sequenced by Ion PGM and MiSeq, respectively. From the 0.46-Gb WGS library (GenBank:SRR1146621; MG-RAST:4532156.3), 4 × 10 6 reads revealed >3 × 10 6 genes, 240 families of prokaryotes, and 155 families of eukaryotes. From the 0.68-Gb 16S library (GenBank:SRR2124995; MG-RAST:4631271.3; EMB:2184), 4 × 10 6 reads revealed 259 families of eubacteria. Results of the WGS and 16S analyses were highly consistent and indicated that more than half of the bacterial sequences were Proteobacteria , predominantly Comamonadaceae . The most numerous genera in this group were Acidovorax (including iron oxidizers, nitrotolulene degraders, and plant pathogens), which accounted for 10 % of assigned bacterial reads. Polaromonas were another 6 % of all bacterial reads, with many assignments to groups capable of degrading polycyclic aromatic hydrocarbons. Albidiferax (iron reducers) and Variovorax (biodegraders of a variety of natural biogenic compounds as well as anthropogenic contaminants such as polycyclic aromatic hydrocarbons and endocrine disruptors) each accounted for an additional 3 % of bacterial reads. Comparison of these data to other publically-available aquatic metagenomes revealed that this stretch of James River is highly similar to the upper Mississippi River, and that these river systems are more similar to aquaculture and sludge ecosystems than they are to lakes or to a pristine section of the upper Amazon River. Taken together, these analyses exposed previously unknown aspects of microbial biodiversity, documented the ecological responses of microbes to urban effects, and revealed the noteworthy presence of 22 human-pathogenic bacterial genera (e.g., Enterobacteriaceae , pathogenic Pseudomonadaceae , and ‘ Vibrionales' ) and 6 pathogenic eukaryotic genera (e.g., Trypanosomatidae and Vahlkampfiidae). This information about pathogen diversity may be used to promote human epidemiological studies, enhance existing water quality monitoring efforts, and increase awareness of the possible health risks associated with recreational use of James River.

Virginia↗
USGS2020

Subspecies differentiation in an enigmatic chaparral shrub species

Premise Delimiting biodiversity units is difficult in organisms in which differentiation is obscured by hybridization, plasticity, and other factors that blur phenotypic boundaries. Such work is more complicated when the focal units are subspecies, the definition of which has not been broadly explored in the era of modern genetic methods. Eastwood manzanita ( Arctostaphylos glandulosa Eastw.) is a widely distributed and morphologically complex chaparral shrub species with much subspecific variation, which has proven challenging to categorize. Currently 10 subspecies are recognized, however, many of them are not geographically segregated, and morphological intermediates are common. Subspecies delimitation is of particular importance in this species because two of the subspecies are rare. The goal of this study was to apply an evolutionary definition of “subspecies” to characterize structure within Eastwood manzanita. Methods We used publicly available geospatial environmental data and reduced‐representation genome sequencing to characterize environmental and genetic differentiation among subspecies. In addition, we tested whether subspecies could be differentiated by environmentally associated genetic variation. Results Our analyses do not show genetic differentiation among subspecies of Eastwood manzanita, with the exception of one of the two rare subspecies. In addition, our environmental analyses did not show ecological differentiation, though limitations of the analysis prevent strong conclusions. Conclusions Genetic structure within Eastwood manzanita does not correspond to current subspecies circumscriptions, but rather reflects geographic distribution. Our study suggests that subspecies concepts need to be reconsidered in long‐lived plant species, especially in the age of next‐generation sequencing.

American Journal of Botany↗
USGS2016

Assessing models of speciation under different biogeographic scenarios; An empirical study using multi-locus and RNA-seq analyses

Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genus Gopherus . These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNA and four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree. RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

Ecology and Evolution↗
USGS2018

A new tool for studying waterfowl immune and metabolic responses: Molecular level analysis using kinome profiling

Here, we describe the design of an Anas‐ specific kinome peptide array that can be used to study the immunometabolic responses of mallard and American black duck to pathogens, contaminants, and environmental stress. The peptide arrays contain 2,642 unique phosphorylate‐able peptide sequences representing 1,900 proteins. These proteins cover a wide array of metabolic and immunological processes, and 758 Gene Ontology Biological processes are statistically significantly represented on the duck peptide array of those 164 contain the term “metabolic” and 25 “immune.” In addition, we conducted a comparison of mallard to American black duck at a genetic and proteomic level. Our results show a significant genomic and proteomic overlap between these two duck species, so that we have designed a cross‐reactive peptide array capable of studying both species. This is the first reported development of a wildlife species‐specific kinome peptide array.

Ecology and Evolution↗
USGS2019

North-facing slopes and elevation shape asymmetric genetic structure in the range-restricted salamander Plethodon shenandoah

Species with narrow environmental preferences are often distributed across fragmented patches of suitable habitat, and dispersal among subpopulations can be difficult to directly observe. Genetic data collected at population centers can help quantify gene flow, which is especially important for vulnerable species with a disjunct range. Plethodon shenandoah is a Federally Endangered salamander known only from three mountaintops in Virginia, USA. To reconstruct the evolutionary history and population connectivity of this species, we generated both mitochondrial and nuclear data using sequence capture for all three populations and found strong population structure that was independent of geographic distance. Both the nuclear markers and mitochondrial genome indicated a deep split between the most southern population and the combined central and northern population. Although there was some mitochondrial haplotype-splitting between the central and northern populations, there was complete admixture in nuclear markers. This is indicative of either a recent split or current male-biased dispersal among mountain isolates. Models of landscape resistance found that dispersal across north-facing slopes at mid-elevation levels best explain the observed genetic structure among populations. These unexpected results highlight the importance of landscape features in understanding and predicting movement and fragmentation of salamanders across space.

Ecology and Evolution↗