Pervasive gene flow despite strong and varied reproductive barriers in swordtails
Data files
Feb 14, 2025 version files 842.47 MB
-
ABC_mito_tdiv_simulations_accepted_November2024.txt
11.05 KB
-
ABC_simulations_accepted_F2s_Feb2024_ch13.txt
114.67 KB
-
ABC_simulations_accepted_F2s_Feb2024_ch14.txt
390.31 KB
-
ABC_simulations_accepted_F2s_Feb2024_ch6.txt
1.70 MB
-
ABC_simulations_accepted_F2s_Feb2024_ch7.txt
595.13 KB
-
all_xmal_xbir_simulated_results_100sim_burnin.txt
400 B
-
all_xmal_xcor_simulated_results_100sim_burnin.txt
400 B
-
all_xmal_xmoz_simulated_results_100sim_burnin.txt
400 B
-
CHPL-contemporary-ancestry.txt
27.41 KB
-
CHPL-historic-ancestry.txt
5.95 KB
-
CHPL-pheromone-trial-keyfile.xlsx
11.92 KB
-
live-male-behavior-data.xlsx
15.50 KB
-
mother-embyro-ancestry.txt
15.56 KB
-
pheromone-trial-keyfile.xlsx
13.47 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____14.xlsx
239.07 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____15.xlsx
225.94 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____16.xlsx
236.11 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____17.xlsx
231.60 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____18.xlsx
197.61 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____19.xlsx
216.67 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____20.xlsx
242.41 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____21.xlsx
232.73 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____22.xlsx
226.68 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____23.xlsx
231.80 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____24.xlsx
236.33 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____25.xlsx
217.48 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____26.xlsx
232.52 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____27.xlsx
225.62 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____28.xlsx
224.28 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____29.xlsx
204.04 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____30.xlsx
193.44 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____31.xlsx
201.57 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____32.xlsx
208.09 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____33.xlsx
225.22 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____34.xlsx
185.42 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____35.xlsx
226.05 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____36.xlsx
239.74 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____37.xlsx
230.21 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____38.xlsx
225.14 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____39.xlsx
233.85 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____40.xlsx
227.89 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____41.xlsx
240.82 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____42.xlsx
243.45 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____43.xlsx
231.29 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____44.xlsx
242.48 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____45.xlsx
243.45 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____46.xlsx
220.83 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____47.xlsx
207.60 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____48.xlsx
219.98 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____49.xlsx
232.64 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____50.xlsx
234.03 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____51.xlsx
248.17 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____52.xlsx
216.12 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____53.xlsx
225.25 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____54.xlsx
213.22 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____55.xlsx
225.96 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____56.xlsx
209.75 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____57.xlsx
224.99 KB
-
Raw_data-CHPL-pheromone-tracking-Trial____58.xlsx
233.16 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____1.xlsx
108.76 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____2.xlsx
100.87 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____3.xlsx
100.76 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____4.xlsx
93.16 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____5.xlsx
94.55 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____6.xlsx
93.79 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____7.xlsx
79 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____8.xlsx
77.21 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial_____9.xlsx
87.06 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____10.xlsx
82.90 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____11.xlsx
93.49 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____12.xlsx
88.89 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____13.xlsx
93.63 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____14.xlsx
94.41 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____15.xlsx
94.87 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____16.xlsx
93.85 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____17.xlsx
76.06 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____18.xlsx
84.79 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____19.xlsx
99.50 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____20.xlsx
101.08 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____21.xlsx
88.33 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____22.xlsx
83.94 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____23.xlsx
95.04 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____24.xlsx
89.80 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____25.xlsx
85.25 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____26.xlsx
86.89 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____27.xlsx
92.90 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____28.xlsx
93.24 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____29.xlsx
72.71 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____30.xlsx
87.56 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____31.xlsx
84.58 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____32.xlsx
93.25 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____33.xlsx
93.48 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____34.xlsx
98.71 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____35.xlsx
74.31 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____36.xlsx
74.58 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____37.xlsx
99.09 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____38.xlsx
95.40 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____39.xlsx
72.77 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____40.xlsx
72.78 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____41.xlsx
102.18 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____42.xlsx
100.44 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____43.xlsx
84.87 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____44.xlsx
87.24 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____45.xlsx
87.09 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____46.xlsx
88.49 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____47.xlsx
85.24 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____48.xlsx
87.78 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____49.xlsx
90.61 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____50.xlsx
91.62 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____51.xlsx
97.92 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____52.xlsx
74.33 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____53.xlsx
98.64 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____54.xlsx
98.56 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____55.xlsx
102.39 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____56.xlsx
90.78 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____57.xlsx
88.18 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____58.xlsx
92.09 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____59.xlsx
94.67 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____60.xlsx
92.64 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____61.xlsx
93.02 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____62.xlsx
93.50 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____63.xlsx
72.73 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____64.xlsx
69.24 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____65.xlsx
87.96 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____66.xlsx
92.24 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____67.xlsx
97.84 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____68.xlsx
84.30 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____69.xlsx
73.64 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____70.xlsx
75.24 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____71.xlsx
98.88 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____72.xlsx
97.76 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____73.xlsx
86 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____74.xlsx
96.15 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____75.xlsx
93.91 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____76.xlsx
82.50 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____77.xlsx
89.97 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____78.xlsx
88.88 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____79.xlsx
84.86 KB
-
Raw_data-cor-bir-pheromone-tracking-Trial____80.xlsx
78.98 KB
-
README.md
7.58 KB
-
sperm-morphology.xlsx
16.06 KB
-
sperm-motilities.csv
5.54 KB
-
testing-priors-xbir.txt
10.66 KB
-
testing-priors-xcor.txt
7.41 KB
-
video-observations-summary.csv
440 B
-
XcorXbir-F2-fullgenotypes.txt_transposed_ancestry_by_site
36.98 MB
-
XcorXbir-PF2-genotypes-110623.txt_transposed
482.45 MB
-
XcorXbir-PF2-I-24-S255-genotypes-021924.txt_transposed
302.81 MB
Abstract
The evolution of reproductive barriers leads to the formation of new species. However, recent research has demonstrated that hybridization has been pervasive across the tree of life even in the presence of strong barriers. Using swordtail fishes (genus Xiphophorus), an emerging model system, we document overlapping mechanisms that act as barriers to gene flow between X. birchmanni and X. cortezi by combining genomic sequencing from natural hybrid populations, experimental laboratory crosses, behavioral assays, sperm performance, and developmental studies. We show that assortative mating plays a role in maintaining subpopulations with distinct ancestry within natural hybrid populations. Using F2 hybrids we identify several genomic regions that strongly impact hybrid viability. Strikingly, two of these regions underlie genetic incompatibilities in hybrids between X. birchmanni and its sister species X. malinche. Our results demonstrate that ancient hybridization has played a role in the origin of this shared genetic incompatibility. Moreover, ancestry mismatch at these incompatible regions has remarkably similar consequences for phenotypes and hybrid survival in X. cortezi × X. birchmanni hybrids as in X. malinche × X. birchmanni hybrids. Our findings identify varied reproductive barriers that shape genetic exchange between naturally hybridizing species and highlight the complex evolutionary outcomes of hybridization.
https://doi.org/10.5061/dryad.dbrv15f8v
Description of the data and file structure
Please see the Main Text and Supplementary Information from Aguillon et al. for details on analyses using these data files.
Note on abbreviations across data files: "CHPL" = Chapulhuacanito (the primary sampling location); "cor" or "xcor" = X. cortezi; "bir" or "xbir" = X. birchmanni; "corXbir" or "XcorXbir" = hybrid between X. cortezi and X. birchmanni; "NA" = indicates a missing value
- Files containing an estimate of overall genomic ancestry:
CHPL-contemporary-ancestry.txt
(contemporary samples from Chapulhuacanito in years 2021-2022),CHPL-historic-ancestry.txt
(historical samples from Chapulhuacanito in years 2003, 2006, 2017),mother-embryo-ancestry.txt
(paired mother/embryo samples),testing-priors-xbir.txt
(testing the use of a uniform versus cluster-specific admixture prior with the ancestry HMM in the Chapulhuacanito X. birchmanni cluster),testing-priors-xcor.txt
(testing the use of a uniform versus cluster-specific admixture prior with the ancestry HMM in the Chapulhuacanito admixed cluster). These files contain columns on # of ancestry informative markers with X. cortezi (cor_count) or X. birchmanni (bir_count) ancestry, the calculated hybrid index, heterozygosity, # of sequencing reads (read_count), and mitochondrial haplotype (mtDNA; 0 = X. birchmanni, 2 = X. cortezi) for each sampled individual (sampleID). For the files testing uniform versus cluster-specific admixture priors, columns end in either _U (uniform prior) or _N (not uniform prior, e.g., cluster-specific prior) to indicate the difference between ancestry HMM runs. - Files containing F2 genotypes:
XcorXbir-PF2-genotypes-110623.txt``transposed
* *andXcorXbir-PF2-I-24-S255-genotypes-021924.txt_transposed
contain hard genotype calls for all ancestry informative markers (AIMs) for X. cortezi x X. birchmanni F2 hybrids. These files are organized with individual F2s as columns and AIMs as rows. Included genotypes are 1 (homozygous X. birchmanni), 2 (heterozygous), or 3 (homozygous X. cortezi). Additionally, the fileXcorXbir-F2-fullgenotypes.txt_transposed_ancestry_by_site
has taken these genotype calls and used them to determine average ancestry (across all the F2s) at all sites along the genome. The files include columns for the chromosome (chr or group), position along the chromosome in basepairs (pos or position), hybrid index, and # of individuals with sufficient sequencing reads at that position (indivs_cov). - Files related to sperm morphology and motility:
sperm-morphology.xlsx
(data on sperm head length, sperm head width, and sperm mid-length) andsperm-motilities.csv
(data on sperm VSL, VCL, VAP, STR, and progressive motility). Both files include data taken from four males each of X. birchmanni, X. cortezi, and their lab-produced F1 and F2 hybrids. Abbreviations: head_l = head length, head_w = head width, mid_l = midpiece length, VSL = straight-line velocity between the start and end of the sperm path, VCL = curvilinear velocity along the actual sperm path, VAP = the velocity over a smoothed sperm path, PM = progressive motility (a measure of the proportion of sperm swimming in a mostly straight line), STR = path straightness (VSL/VAP * 100). - Files related to behavioral assays: Files with the name format
Raw_data-cor-bir-pheromone-tracking-Trial____##.xlsx
orRaw_data-CHPL-pheromone-tracking-Trial____##.xlsx
are raw output files from the tracking program EthoVision from processing video recordings of mate choice trials. Starting at row 35, these raw data files include the X,Y position of the fish within the experimental tank at 0.2 or 0.4 s intervals, along with a 0 (absent) or 1 (present) value for if the tested individual was found within either of the two association zones at the time point. This data is preceded by header information (row 1-33) detailing settings used to process the raw video recordings in EthoVision. The useful variables from the header include Trial name (the # at the end of the file name, and matching the "file_name" column in the associated keyfile), Video file (with info matching the "trial_name" column in the associated keyfile), and Arena settings (which indicates if the processed file was for the A or B trial lane, and matching the "trial_lane" column in the associated keyfile). All associated analyses were conducted using the presence/absence data below the header using custom scripts. The two associated "keyfiles" (pheromone-trial-keyfile.xlsx
andCHPL-pheromone-trial-keyfile.xlsx
) contain columns on individual fish used in trials, date/time trials were conducted, and information connecting trials to EthoVision raw output. Note that the correct date/time of the conducted trials is included in the keyfiles, not the raw Ethovision output (as the header information in the latter reflects the processing of the trial, not the recording of the trial). The filelive-male-behavior-data.xlsx
includes results on behavioral trials that were conducted with live male as stimuli. The filevideo-observations-summary.csv
includes results from video observations in the Chapulhuacanito population. - Files containing output from ABC simulations:
ABC_simulations_accepted_F2s_Feb2024_ch6.txt
(ABC simulation results for the genomic region containing the gene ndufa13)ABC_simulations_accepted_F2s_Feb2024_ch13.txt
(ABC simulation results for the genomic region containing the gene ndufs5),ABC_simulations_accepted_F2s_Feb2024_ch7.txt
(ABC simulation results for a genomic region with distorted ancestry proportions on chromosome 7),ABC_simulations_accepted_F2s_ch14.txt
(ABC simulation results for a genomic region with distorted ancestry proportions on chromosome 14). These simulations were for selection and the relevant columns are s (selection) and h (dominance). Additionally,ABC_mito_tdiv_simulations_accepted_November2024.txt
includes ABC simulation results to estimate the divergence time between X. cortezi and X. malinche mitochondrial haplotypes. The relevant column is Generations. - Files containing output from SLiM simulations:
all_xmal_xcor_simulated_results_100sim_burnin.txt
(results from SLiM simulations comparing X. malinche with X. cortezi),all_xmal_xmoz_simulated_results_100sim_burnin.txt
(results from SLiM simulations comparing X. malinche with X. montezumae),all_xmal_xbir_simulated_results_100sim_burnin.txt
(results from SLiM simulations comparing X. malinche with X. birchmanni).
Genomic Data Availability
New genomic data generated for this project is available through the NCBI BioProject PRJNA1106506.
Code Availability
Data processing scripts to work with all included files are available on GitHub (https://github.com/stepfanie-aguillon/swordtail-reproductive-barriers-NatureEE2025, https://github.com/Schumerlab/Lab_shared_scripts) or Zenodo (https://doi.org/10.5281/zenodo.14736259). Note that data file names in the Dryad Repository are consistent with those used in the GitHub/Zenodo Repository code.