Mutation accumulation in a hybrid parthenogenetic vertebrate
Data files
Aug 12, 2026 version files 490.09 MB
-
alignments_nostops.nex
372.93 KB
-
alignments_nostops.nex.RELAX.json
190.94 KB
-
all_orthologs_OUT_0.txt
971.70 KB
-
all_orthologs_OUT_2.txt
972.08 KB
-
all_orthologs_OUT_4.txt
972.14 KB
-
all_phased_orthologs_alignments.zip
232.18 KB
-
all_species_OUT_0.txt
1.27 MB
-
all_species_OUT_2.txt
1.28 MB
-
all_species_OUT_3.txt
1.28 MB
-
all_species_OUT_6.txt
1.28 MB
-
fourratio_tess_parents.zip
1.36 MB
-
JAM_014_cdhit_h.fasta
50.55 MB
-
JOR1048_collapsed_h2.fasta
25.07 MB
-
JOR1112_collapsed_h2.fasta
63.12 MB
-
KLC_066_cdhit_h.fasta
29.85 MB
-
KLC_148_cdhit.fasta
46.29 MB
-
KLC_157_cdhit.fasta
40.63 MB
-
KLC_162_cdhit_h.fasta
50.07 MB
-
KLC_163_header.fasta
31.98 MB
-
KLC_167_cdhit.fasta
40.78 MB
-
KLC_168_cdhit_h.fasta
19.39 MB
-
marm_relaxed_pbsv.vcf
6.49 MB
-
nonparental_sexual_ortholog_alignments.zip
77.11 KB
-
nonparental_sexual_tree.txt
100 B
-
oneratio_tess_parents.zip
1.31 MB
-
paml_annotated_david.xlsx
1 MB
-
README.md
6.87 KB
-
relax_out.zip
1.39 MB
-
sept_relaxed_pbsv.vcf
36.70 MB
-
simulations.zip
33.85 MB
-
tree_phased.txt
55 B
-
tworatio_tess_parents.zip
1.36 MB
Aug 12, 2026 version files 490.09 MB
-
alignments_nostops.nex
372.93 KB
-
alignments_nostops.nex.RELAX.json
190.94 KB
-
all_orthologs_OUT_0.txt
971.70 KB
-
all_orthologs_OUT_2.txt
972.08 KB
-
all_orthologs_OUT_4.txt
972.14 KB
-
all_phased_orthologs_alignments.zip
232.18 KB
-
all_species_OUT_0.txt
1.27 MB
-
all_species_OUT_2.txt
1.28 MB
-
all_species_OUT_3.txt
1.28 MB
-
all_species_OUT_6.txt
1.28 MB
-
fourratio_tess_parents.zip
1.36 MB
-
JAM_014_cdhit_h.fasta
50.55 MB
-
JOR1048_collapsed_h2.fasta
25.07 MB
-
JOR1112_collapsed_h2.fasta
63.12 MB
-
KLC_066_cdhit_h.fasta
29.85 MB
-
KLC_148_cdhit.fasta
46.29 MB
-
KLC_157_cdhit.fasta
40.63 MB
-
KLC_162_cdhit_h.fasta
50.07 MB
-
KLC_163_header.fasta
31.98 MB
-
KLC_167_cdhit.fasta
40.78 MB
-
KLC_168_cdhit_h.fasta
19.39 MB
-
marm_relaxed_pbsv.vcf
6.49 MB
-
nonparental_sexual_ortholog_alignments.zip
77.11 KB
-
nonparental_sexual_tree.txt
100 B
-
oneratio_tess_parents.zip
1.31 MB
-
paml_annotated_david.xlsx
1 MB
-
README.md
6.88 KB
-
relax_out.zip
1.39 MB
-
sept_relaxed_pbsv.vcf
36.70 MB
-
simulations.zip
33.85 MB
-
tree_phased.txt
55 B
-
tworatio_tess_parents.zip
1.36 MB
Abstract
Asexual lineages are thought to experience elevated extinction rates compared with sexual species, yet direct evidence for the underlying genetic causes remains scarce. Muller’s ratchet predicts that the absence of recombination in asexual organisms facilitates the accumulation of deleterious mutations, thereby reducing long-term fitness. Here, we test this hypothesis in the hybrid-origin, parthenogenetic whiptail lizard Aspidoscelis tesselatus by integrating short-read RNA-seq and long-read IsoSeq data from both the asexual lineage and its parental sexual species. We reconstructed phased transcripts for A. tesselatus to quantify mutation accumulation relative to the parental sexual species. Comparative analyses revealed elevated ω ratios in both parental genomic complements (subgenomes) of the parthenogenetic lineage, consistent with accelerated accumulation of nonsynonymous mutations. Structural variant analyses identified multiple indels in expressed transcripts predicted to disrupt protein domains. Functional annotation indicated that genes affected by both single-nucleotide variants and indels were enriched for roles in chromatin organization, apoptosis regulation, and transcriptional control. While both parental subgenomes showed similar evolutionary patterns, the maternal complement exhibited more structural and missense mutations than the paternal complement. Together, these results provide evidence that mutations accumulate in asexual A. tesselatus in genes involved in core cellular functions, supporting theoretical predictions that Muller’s ratchet contributes to mutation accumulation in asexual lineages.
https://doi.org/10.5061/dryad.cnp5hqcjx
Description of the data and file structure
This dataset was generated to investigate whether mutations accumulate in the asexual whiptail lineage, Aspidoscelis tesselatus, when compared with their sexual parental species. All Aspidoscelis tesselatus samples, A. arizonae, A. gularis, and A. sexlineatus are short read RNA-seq transcriptomes from liver tissue, whereas the A. marmoratus and A. septemvittatus fasta files are IsoSeq transcriptomes from liver tissue.
Files and variables
- Samples
The sample names used throughout these files are from collector IDs. The museum accession code for each sample is listed below:- KLC157 = UTA-R-68384: A. sexlineatus
- JOR1048 = UTA-R-68269: A. marmoratus
- KLC163 = UTA-R-68389: A. tesselatus
- KLC167 = UTA-R-68066: A. gularis
- JOR1112 = UTA-R-68355: A. septemvittatus
- KLC168 = UTA-R-68101: A. tesselatus
- KLC162 = UTA-R-68388: A. tesselatus
- KLC066 = UTA-R-68396: A. tesselatus
- JAM014 = UTA-R-68414: A. tesselatus
- KLC148 = UTA-R-68184: A. arizonae
Below is a summary of the included files:
- Tree Files
tree_phased.txt: Tree used for A. tesselatus and parentals PAML and RELAX, generated from homologizer and IQTREE2- Variables:
tes_marm {Asexual}: A. tesselatus-marmoratus subgenomeJOR_1048 {Sexual}: A. marmoratustes_sept {Asexual}: A. tesselatus-septemvittatus subgenomeJOR_1112 {Sexual}: A. septemvittatus
- Variables:
nonparental_sexual_tree.txt: Tree used for parental, A. tesselatus, and nonparental sexual comparison in PAML- Variables:
JOR_1048: A. marmoratustes_marm: A. tesselatus-marmoratus subgenomeKLC_148: A. arizonaeKLC_157: A. sexlineatusKLC_167: A. gularisJOR_1112: A. septemvittatustes_sept: A. tesselatus-septemvittatus subgenome
- Variables:
- Raw Data
all_phased_orthologs_alignments.zip: Ortholog alignments with the phased A. tesselatus and parentals comparison in phylip format. These are output files from the OrthoFinder and homologizer phasing pipeline.nonparental_sexual_ortholog_alignments.zip: Ortholog alignments with the nonparental sexual comparison. This folder includes only the ortholog files in phylip format for PAML.KLC_168_cdhit_h.fasta: A. tesselatus RNAseq transcriptomeKLC_163_header.fasta: A. tesselatus RNAseq transcriptomeKLC_066_cdhit_h.fasta: A. tesselatus RNAseq transcriptomeKLC_162_cdhit_h.fasta: A. tesselatus RNAseq transcriptomeJAM_014_cdhit_h.fasta: A. tesselatus RNAseq transcriptomeKLC_167_cdhit.fasta: A. gularis RNAseq transcriptomeKLC_157_cdhit.fasta: A. sexlineatus RNAseq transcriptomeKLC_148_cdhit.fasta: A. arizonae RNAseq transcriptomeJOR1112_collapsed_h2.fasta: A. septemvittatus IsoSeq transcriptomeJOR1048_collapsed_h2.fasta: A. marmoratus IsoSeq transcriptomealignments_nostops.nex: RELAX input for all concatenated orthologs with phased A. tesselatus and parentals
- Results
relax_out.zip: RELAX output files for each of the orthologs in theall_phased_ortholog_alignments.zipfolder.oneratio_tess_parents.zip: PAML result files for M0 model with A. tesselatus and parentals. This folder includes only the results files in text format for each of the orthologs used in the analysis.tworatio_tess_parents.zip: PAML result files for the M2 two-ratio model with A. tesselatus and parentals. This folder includes only the results files in text format for each of the orthologs used in the analysis.fourratio_tess_parents.zip: PAML result files for M2 four-ratio model with A. tesselatus and parentals. This folder includes only the results files in text format for each of the orthologs used in the analysis.paml_annotated_david.xlsx: Excel spreadsheet with ortholog annotations from InterProScan and DAVID functional annotationall_orthologs_OUT_0.txt: PAML output with M0 model with concatenated orthologs with phased A. tesselatus and parentalsalignments_nostops.nex.RELAX.json: RELAX output for all concatenated orthologs with phased A. tesselatus and parentalsall_orthologs_OUT_4.txt: PAML output with M2 four-ratio model with concatenated orthologs with phased A. tesselatus and parentalsall_orthologs_OUT_2.txt: PAML output with M2 two-ratio model with concatenated orthologs with phased A. tesselatus and parentalsall_species_OUT_0.txt: PAML output file from M0 model with nonparental sexual comparisonall_species_OUT_2.txt: PAML output file from M2 two-ratio model with nonparental sexual comparison with concatenated orthologsall_species_OUT_3.txt: PAML output file from M2 three-ratio model with nonparental sexual comparison with concatenated orthologsall_species_OUT_6.txt: PAML output file from M2 seven-ratio model with nonparental sexual comparison with concatenated orthologs.sept_relaxed_pbsv.vcf: A. tesselatus-septemvittatus phased structural variantsmarm_relaxed_pbsv.vcf: A. tesselatus-marmoratus phased structural variantssimulations.zip: Evolver and SLiM simulation resultsevolver_reps: Alignment output files from evolverreps_paml0: PAML output from M0 model with evolver simulated datasetsreps_paml2: PAML output from M2 model with evolver simulated datasetsreps_paml4: PAML output from M4 model with evolver simulated datasets
unequal_ne: Alignment output files from SLiM simulation of unequal Ne between asexual and sexual populationsreps_paml0: PAML output from M0 model with SLiM unequal Ne datasetsreps_paml2: PAML output from M2 model with SLiM unequal Ne datasetsreps_paml4: PAML output from M4 model with SLiM unequal Ne datasets
equal_ne: Alignment output files from SLiM simulation of equal Ne between asexual and sexual populationsreps_paml0: PAML output from M0 model with SLiM equal Ne datasetsreps_paml2: PAML output from M2 model with SLiM equal Ne datasetsreps_paml4: PAML output from M4 model with SLiM equal Ne datasets
forwardNe.slim: SLiM5 model framework for Ne simulations
- Access information
- Other publicly accessible locations of the data:
- NCBI BioProject: PRJNA1381241
- Scripts at GitHub: https://github.com/zoemullersnakes/MutationAccumulationWhiptails
- Other publicly accessible locations of the data:
Changes after Aug 12, 2026: Edited README.md.
