Data from: Nemo knows: Clownfishes differentiate cryptic host species across fine and broad geographic scales and reveal a potential adaptive radiation in the clownfish-hosting sea anemones
Data files
Apr 24, 2026 version files 142.55 MB
-
dadi_JFS.py
2.13 KB
-
dadiModels.py
57.02 KB
-
EQ_dadi_top_model_Results_Summary.txt
10.68 KB
-
EQ_RADseq_Stats.txt
20.73 KB
-
EQ_RADseq.ustr
1.54 MB
-
EQ_UCE_loci_genetrees.trees
5.35 MB
-
eqJapan-UCE-assemblyStats.csv
12.55 KB
-
eqMantu-75p-msa.fas
135.52 MB
-
README.md
4.82 KB
-
tabS1_EQmorphs.csv
4.83 KB
-
timeCal-unpart.timetree.nex
22.58 KB
Abstract
The symbiosis between clownfish (or anemonefishes) and their host sea anemones ranks among the most recognizable animal interactions on the planet. Found on coral reef habitats across the Indian and Pacific Oceans, 28 recognized species of clownfishes adaptively radiated from a common ancestor to live obligately with only 10 nominal species of host sea anemones. Are the host sea anemones truly less diverse than clownfish? Did the symbiosis with clownfish trigger a reciprocal adaptive radiation in sea anemones, or minimally, a co-evolutionary response to the mutualism? To address these questions, we combined fine- and broad-scale biogeographic sampling with multiple independent genomic datasets for the bubble-tip sea anemone, Entacmaea quadricolor—the most common clownfish host anemone throughout the Indo-West Pacific. Fine-scale sampling and restriction site-associated DNA sequencing (RADseq) throughout the Japanese Archipelago revealed three highly divergent cryptic species: two of which co-occur throughout the Ryukyu Islands and can be differentiated by the clownfish species they host. Remarkably, broader biogeographic sampling and bait-capture sequencing reveal that this pattern is not simply the result of local ecological processes unique to Japan, but part of a deeper evolutionary signal where some species of E. quadricolor serve as host to the generalist clownfish species Amphiprion clarkii and others serve as host to the specialist clownfish A. frenatus. In total, we delimit at least five cryptic species in E. quadricolor that have diversified within the last five million years. The rapid diversification of E. quadricolor, combined with functional ecological and phenotypic differentiation, supports the hypothesis that this may represent an adaptive radiation in response to mutualism with clownfishes. Our data indicate that clownfishes are not merely settling in locally available hosts but recruiting to specialized host lineages with which they have co-evolved. These findings have important implications for understanding how the clownfish-sea anemone symbiosis has evolved and will shape future research agendas on this iconic model system.
Dataset DOI: 10.5061/dryad.0000000fb
Description of the data and file structure
Files and variables
File: dadi_JFS.py
Description: Python script used for generating the joint folded allele frequency spectra used in our Entacmaea quadriclor demographic modeling analyisis in dadi
File: dadiModels.py
Description: Python script used for generating and simulating the 17 demographic models we developed in dadi for Entacmaea quadricolor
File: EQ_RADseq_Stats.txt
Description: This file contains the summary statistics from ipyrad for the final double digest restriction site associated DNA sequencing (ddRADseq) dataset for Entacmaea quadricolor from the Japanese Archipelago.
File: EQ_dadi_top_model_Results_Summary.txt
Description: this file contains results from demographic modeling in dadi. For each model the top five simulated models with the best log-likelihood AIC scores are presented.
Variables
- Model: Name of the simulated demographic model in dadi
- Replicate: the round and replicate of the simulated model with the highest AIC score
- log-likelihood: log-likelihood score for the specified model
- AIC: Akaike Information Criterion score
- chi-squared: chi-squared score calculated by dadi for the specified model
- theta: Theta parameter calculated for each model in dadi. theta=4Nrefmu*L where mu is the mutation rate, L is the length of sequence, and Nref is the effective size of the reference population
- optimized_params: optimized paremeter values calculated in dadi for each model. The parameters for each model are designated in the file dadiModels.py
File: eqJapan-UCE-assemblyStats.csv
Description: This file contains the summary statistics from the bait-capture sequencing dataset assembly as produced by phyluce.
Variables
- samples: Sample ID
- contigs: Number of contigs per sample
- total bp: total number of nucleotide base pairs (bp) sequenced for that individual sample
- mean length: mean per locus sequence length
- 95 CI length: 95% confidence interval sequence length
- min length: minimum sequence length (in bp)
- max length: maximum sequence length
- median legnth:
- contigs >1kb:
File: EQ_RADseq.ustr
Description: This file is contains the double digest restriction site associated DNA sequencing dataset containing unlinked single nucleotide polymorphisms for Entacmaea quadricolor from Japan. The file is formatted for the genetic clustering program STRUCTURE.
File: tabS1_EQmorphs.csv
Description: This file contains the raw data for the Entacmaea quadricolor phenotypic analysis from the Japanese Archipelago.
Variables
- ID number: sample ID
- TSHAPE: Tentacle shape (bubble with tip; digitiform; bubble rounded)
- TLENGTH: Tentacle length (long or small))
- GP: Anemone group size (solitary vs clonal)
- depthCat: Depth (deep or shallow)
- Color: anemone color (bleached, brown, pale brown, green/brown, and green)
- tipColor: anemone tentacle tip color (pink tip-PT; not- noPT)
- tipPattern: anemone tentacle tip pattern (striated, striated/speckled; dull/matte)
- fishSP: clownfish species residing in the anemone
File: timeCal-unpart.timetree.nex
Description: Time calibrated unpartitioned .nexus tree file for our time calibrated phylogenetic analysis of the Entacmaea quadricolor species complex
File: EQ_UCE_loci_genetrees.trees
Description: Individual input gene trees generated from bait-capture DNA sequencing targeting Ultra Conserved Element (UCE) and Exon loci for Entacmaea quadricolor.
File: eqMantu-75p-msa.fas
Description: FASTA file formatted DNA sequence alignment generated from bait-capture sequencing for Entacmaea quadricolor.
Code/software
dadiModels.py: Python script use for generating 17 demographic models in dadi. This script requires Python, numpy, dadi, Numerics, PhiManip, and Integration packages.
dadiJFS.py: Python script used to create the joint-site frequency spectrum used in dadi demographic modeling simulations. This script requires Python, numpy dadi, datetime, Optimize Functions, and Models3D packages
Access information
Other publicly accessible locations of the data:
- Raw sequence data are deposited on GenBank.
- RADseq data accessible via BioProject PRJNA1198303
- Bait-capture sequence data are accessible via BioProject PRJNA1198303
Data was derived from the following sources:
- RADseq and bait-capture sequence data were generated from DNA extracted from Entacmaea quadricolor tentacle tissue.
