Data from: Revised taxonomy reveals sustained introgression and secondary contact in ancient lake ricefishes
Data files
Jun 11, 2026 version files 72.53 MB
-
File_01_VCF_9603_SNPs.vcf
31.64 MB
-
File_02_VCF_3947_SNPs.vcf
12.99 MB
-
File_03_VCF_4193_SNPs.vcf
10.75 MB
-
File_04_VCF_9259_SNPs.vcf
17.13 MB
-
File_05_est_and_tpl.zip
27.23 KB
-
README.md
2.93 KB
Abstract
Biotic diversification in ancient lakes is shaped by complex geological histories and genetic exchange among populations. The Malili Lake system on Sulawesi Island represents a classic natural laboratory for studying freshwater fish evolution and harbors multiple endemic Oryzias species that diversified under repeated hydrological reorganizations. Previous genomic analyses inferred that two sympatric species in Lake Towuti (O. profundicola and O. loxolepis) experienced a single ancient introgression event from a “ghost lineage” derived from O. marmoratus inhabiting another lake. However, recent taxonomic re-evaluation has revealed the presence of an extant O. marmoratus population within Lake Towuti itself. This finding suggests that the putative ghost lineage may in fact represent a living population co-occurring in the lake, calling for a re-examination of the introgression history and speciation mode in Lake Towuti. By incorporating newly generated ddRAD-seq data from the true O. marmoratus in Lake Towuti, we reanalyzed phylogenetic relationships and population genetic structure among Malili Lake Oryzias. Previously reported major phylogenetic relationships and inter-lake introgression patterns were largely reproduced. In contrast, TreeMix and f4-statistic analyses revealed that introgression signals previously attributed to a “ghost lineage” into O. profundicola and O. loxolepis instead originated from the extant O. marmoratus population coexisting within Lake Towuti. Demographic model comparisons explicitly incorporating within-lake gene flow further supported a scenario in which O. profundicola and O. loxolepis diverged in allopatry, subsequently came into secondary contact within Lake Towuti, and later experienced additional gene flow following secondary contact with O. marmoratus that entered the lake. Our results demonstrate that introgression from theO. marmoratus lineage intoO. profundicola andO. loxolepis was not a single ancient event, but rather a more sustained process. This finding highlights the critical importance of taxonomic resolution for accurately inferring introgression and divergence history. Comparative studies across other ancient lakes on Sulawesi will be valuable for understanding how the timing and nature of gene flow from third lineages influence patterns of population divergence and the strength of reproductive isolation.
README for data set from:
Yashima et al. 2026 Revised taxonomy reveals sustained introgression and secondary contact in ancient lake ricefishes. BMC Ecology and Evolution.
The data set accompanies five Dryad data files.
Uploaded: December 2025
Revised: June 2026
Contact: Kazunori Yamahira (yamahira@cs.u-ryukyu.ac.jp)
File_01_VCF_9603_SNPs.vcf
The VCF file containing 9,603 SNPs among seven Malili Oryzias populations and two outgroups (O. dopingdopingensis and O. nebulosus), which was used for estimating the maximum-likehood phylogeny by RAxML-NG. REF and ALT alleles represent major and minor alleles, respectively, rather than reference and alternative. This is due to the specification of the populations program of the Stacks pipeline that was used to produce the file. The sequence files are available from DDBJ/EMBL-EBI/NCBI Sequence Read Archive under the accession number DRA010303 and DRA025938.
File_02_VCF_3947_SNPs.vcf
The VCF file containing 3,947 SNPs among seven Malili Oryzias populations and two outgroups (O. dopingdopingensis and O. nebulosus), which was used for estimating the SNAPP species tree and for examining possible introgression among Malili populations by TreeMix. REF and ALT alleles represent major and minor alleles, respectively, rather than reference and alternative. This is due to the specification of the populations program of the Stacks pipeline that was used to produce the file. The sequence files are available from DDBJ/EMBL-EBI/NCBI Sequence Read Archive under the accession number DRA010303 and DRA025938.
File_03_VCF_4193_SNPs.vcf
The VCF file containing 4,193 SNPs among seven Malili Oryzias populations, which was used for estimating the population structure by ADMIXTURE and for building a phylogenetic network by SplitsTree. REF and ALT alleles represent major and minor alleles, respectively, rather than reference and alternative. This is due to the specification of the populations program of the Stacks pipeline that was used to produce the file. The sequence files are available from DDBJ/EMBL-EBI/NCBI Sequence Read Archive under the accession number DRA010303 and DRA025938.
File_04_VCF_9259_SNPs.vcf
The VCF file containing 9,259 biallelic SNPs on 11,330 RAD loci among four Malili populations and one outgroup (O. nebulosus), which was used to create two-dimensional joint minor allele site frequency spectrum for fastsimcoal2 runs. REF and ALT alleles represent major and minor alleles, respectively, rather than reference and alternative. This is due to the specification of the populations program of the Stacks pipeline that was used to produce the file. The sequence files are available from DDBJ/EMBL-EBI/NCBI Sequence Read Archive under the accession number DRA010303 and DRA025938.
File_05_est_and_tpl.zip
The est and tpl files of each of the demographic models compared in fastsimcoal2.
