Dataset Information

Base-substitution mutation rate across the nuclear genome of Alpheus snapping shrimp and the timing of isolation by the Isthmus of Panama.

ABSTRACT:

Background

The formation of the Isthmus of Panama and final closure of the Central American Seaway (CAS) provides an independent calibration point for examining the rate of DNA substitutions. This vicariant event has been widely used to estimate the substitution rate across mitochondrial genomes and to date evolutionary events in other taxonomic groups. Nuclear sequence data is increasingly being used to complement mitochondrial datasets for phylogenetic and evolutionary investigations; these studies would benefit from information regarding the rate and pattern of DNA substitutions derived from the nuclear genome.

Results

To estimate the genome-wide neutral mutation rate (µ), genotype-by-sequencing (GBS) datasets were generated for three transisthmian species pairs in Alpheus snapping shrimp. A range of bioinformatic filtering parameters were evaluated in order to minimize potential bias in mutation rate estimates that may result from SNP filtering. Using a Bayesian coalescent approach (G-PhoCS) applied to 44,960 GBS loci, we estimated µ to be 2.64E-9 substitutions/site/year, when calibrated with the closure of the CAS at 3 Ma. Post-divergence gene flow was detected in one species pair. Failure to account for this post-split migration inflates our substitution rate estimates, emphasizing the importance of demographic methods that can accommodate gene flow.

Conclusions

Results from our study, both parameter estimates and bioinformatic explorations, have broad-ranging implications for phylogeographic studies in other non-model taxa using reduced representation datasets. Our best estimate of µ that accounts for coalescent and demographic processes is remarkably similar to experimentally derived mutation rates in model arthropod systems. These results contradicted recent suggestions that the closure of the Isthmus was completed much earlier (around 10 Ma), as mutation rates based on an early calibration resulted in uncharacteristically low genomic mutation rates. Also, stricter filtering parameters resulted in biased datasets that generated lower mutation rate estimates and influenced demographic parameters, serving as a cautionary tale for the adherence to conservative bioinformatic strategies when generating reduced-representation datasets at the species level. To our knowledge this is the first use of transisthmian species pairs to calibrate the rate of molecular evolution from GBS data.

SUBMITTER: Silliman K

PROVIDER: S-EPMC8164322 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BACKGROUND:Lutzomyia gomezi (Nitzulescu, 1931) is one of the main Leishmania (Vianna) panamensis vectors in Panama, and despite its medical significance, there are no population genetic studies regarding this species. In this study, we used the sequences of the mitochondrial gene cytochrome b/start of NADH1 and the nuclear elongation gene ?-1 in order to analyze genetic variation and phylogeographic structure of the Lu. gomezi populations. METHODS:A total of 86 Lu. gomezi individuals were captured in 38 locations where cutaneous leishmaniasis occurred. DNA was extracted with phenol/chloroform methods and amplification of genes was performed using PCR primers for mitochondrial and nuclear markers. RESULTS:We found a total of 37 and 26 haplotypes of mitochondrial and nuclear genes, high haplotype diversity (h) for all three populations were detected with both molecular markers. Nucleotide diversity (?) was estimated to be high for all three populations with the mitochondrial marker, which was opposite to the estimate with the nuclear marker. In the AMOVA ?st recorded moderate (mitochondrial) and small (nuclear) population structure with statistical significance among populations. The analysis of the fixation index (Fst) used to measure the differentiation of populations showed that with the exception of the population located in the region of Bocas del Toro, the other populations presented with minor genetic differentiation. The median-Joining network of the mitochondrial marker reveled three clusters and recorded four haplotypes exclusively of localities sampled from Western Panama, demonstrating strong divergence. We found demographic population expansion with Fu´s Fs neutrality test. In the analysis mismatch distribution was observed as a bimodal curve. CONCLUSION:Lu. gomezi is a species with higher genetic pool or variability and mild population structure, due to possible capacity migration and local adaptation to environmental changes or colonization potential. Thus, knowledge of the genetic population and evolutionary history is useful to understand the implications of different population genetic structures for cutaneous leishmaniasis epidemiology.