Candidate-species delimitation in Desmognathus salamanders reveals gene flow across lineage boundaries, confounding phylogenetic estimation and clarifying hybrid zones
Data files
Feb 22, 2022 version files 55.37 MB
-
accessions.csv
31.07 KB
-
conanti_q.mat.csv
22.21 KB
-
concat.treefile
62.13 KB
-
concord.cf.tree
70.49 KB
-
cumberland_q.mat.csv
386 B
-
desmo900_dsuite_species_tree_Fvals.txt
16.54 KB
-
desmo900_formalin.parts
2.92 KB
-
desmo900_formalin.parts.treefile
54.81 KB
-
desmo900_formalin.phy
24.91 MB
-
desmo900.csv
121.78 KB
-
fuscus_cp_q.mat.csv
12.74 KB
-
fuscus_mtn_q.mat.csv
6.71 KB
-
loci.treefile
10.74 MB
-
loci.zip
19.25 MB
-
mountain_q.mat.csv
821 B
-
nantahala_q.mat.csv
5.62 KB
-
ocoee_q.mat.csv
15.68 KB
-
orestes_q.mat.csv
4.07 KB
-
ouachita_q.mat.csv
782 B
-
pigmy_q.mat.csv
1.15 KB
-
pisgah_q.mat.csv
12.75 KB
-
README.txt
1.58 KB
-
seepage_q.mat.csv
5.13 KB
-
USNM468094_q.mat.csv
13.04 KB
-
USNM468095_q.mat.csv
12.83 KB
Abstract
Dusky Salamanders (genus Desmognathus) currently comprise only 22 described, extant species. However, recent mitochondrial and nuclear estimates indicate the presence of up to 49 candidate species based on ecogeographic sampling. Previous studies also suggest a complex history of hybridization between these lineages. Studies in other groups suggest that disregarding admixture may affect both phylogenetic inference and clustering-based species-delimitation. With a dataset comprising 233 Anchored Hybrid Enrichment (AHE) loci sequenced for 896 Desmognathus specimens from all 49 candidate species, we test three hypotheses regarding i) species-level diversity, ii) hybridization and admixture, and iii) misleading phylogenetic inference. Using phylogenetic and population-clustering analyses considering gene flow, we find support for at least 47 candidate species in the phylogenomic dataset, some of which are newly characterized here while others represent combinations of previously named lineages that are collapsed in the current dataset. Within these, we observe significant phylogeographic structure, with up to 64 total geographic genetic lineages, many of which hybridize either narrowly at contact zones or extensively across ecological gradients. We find strong support for both recent admixture between terminal lineages and ancient hybridization across internal branches. This signal appears to distort concatenated phylogenetic inference, wherein more heavily admixed terminal specimens occupy apparently artifactual early diverging topological positions, occasionally to the extent of forming false clades of intermediate hybrids. Additional geographic and genetic sampling and more robust computational approaches will be needed to clarify taxonomy, and to reconstruct a network topology to display evolutionary relationships in a manner that is consistent with their complex history of reticulation. --
Methods
Anchored Hybrid Enrichment
Usage notes
All raw data and output files from the sequencing runs and analyses described in the article are included here.