Data from: Exploring incongruence in a plant genus with relatively recent polyploidy and cytoplasmic discordance
Data files
Jun 13, 2026 version files 5.22 GB
-
570_accessions_100_Busco_concatenated.zip
35.08 MB
-
Dataset_5_B-genome.zip
17.98 MB
-
Dataset1_61_accessions.zip
4.02 GB
-
Dataset2_27_accessions.zip
28.05 MB
-
Dataset3_7_accessions.zip
1.06 GB
-
Dataset4_A-genome.zip
18.39 MB
-
Dataset6_DEHI_genome.zip
20.17 MB
-
Divergence_dates.zip
61.48 KB
-
orthofinder.zip
9.22 MB
-
README.md
3.43 KB
-
reference_loci.zip
4.59 MB
Abstract
The eukaryotic genome has been described as a collection of different histories; for any set of taxa one of these histories is the record of cladogenic events that together comprise the species tree. Among the other histories expected to occur are those attributable to deep coalescence/lineage sorting; to biological causes such as introgression and horizontal transfer; or to pseudo-orthology, long branch attraction, and other technical issues. Gene tree topology space is the portion of tree space occupied by the gene trees reconstructed for a particular dataset of sampled genetic loci. Because coalescent theory predicts that the species tree topology will generally be the most frequent among gene trees, a reasonable expectation is that there will be a peak in gene tree topology space at the species tree topology, with secondary peaks present due to trees tracking other histories. Gene tree topology space in the small (~30 species, including the cultivated soybean) legume genus, Glycine should not only have signals from the species tree and from lineage sorting, but also from a likely introgression event that created incongruence between the plastid and nuclear genomes. Additionally, Glycine is the product of a relatively recent (<13 million years) whole genome duplication, raising the possibility of pseudo-orthology. We explored this space using a set of 2389 nuclear genes and representative accessions from a 570-taxon concatenation tree, reconstructing gene trees for all nuclear loci and from complete plastid genomes and partial mitochondrial genomes. Species trees (ASTRAL) and maximum likelihood (ML) concatenation trees were congruent for a 61-taxon dataset but were incongruent with organellar genome trees. Gene tree topology space was flat: No topology was represented by more than one gene tree. This was also true for a reduced dataset of 27 taxa; only when the dataset was reduced to six ingroup taxa were multiple gene trees having the species tree topology observed, along with a topology congruent with the chloroplast genome topology, presumably representing nuclear loci introgressed along with the plastome. Clustering failed to identify any regional differentiation of gene tree topology space populated by loci with similar topologies. Pseudo-orthology did not contribute meaningfully to incongruence, in agreement with recent modeling work that minimizes concerns about this phenomenon. Clearly, different genes have different historical signals, but these signals are complex and exist at the level of clades within trees rather than as entire gene trees.
https://doi.org/10.5061/dryad.2v6wwpzxn
Description of the data and file structure
Files and variables
There are multiple files/folders for the different datasets investigating gene tree space in Glycine.
File: orthofinder.zip
Description: This contains the alignment files, rooted gene trees, and collapsed gene trees (nodes with less than 70% support were collapsed) for the accessions of* Glycine* from Zhuang et al with assembled genomes. This also includes the species tree estimated from OrthoFinder.
File: Dataset6_DEHI_genome.zip
Description: This contains the alignment files and tree files for 29 accessions that make up the DEHI-genome group in Glycine.
File: Dataset3_7_accessions.zip
Description: This contains the alignment files, rooted gene trees, and collapsed gene trees (nodes with less than 70% support were collapsed) for 7 accessions of Glycine from Zhuang et al 2022, plus two outgroups,* G. max and Aphicarpeae edgworthii. The input files for the ABBA-BABA test are also included, along with summary files for each locus and accession across the 7 accessions.
File: reference_loci.zip
Description: Reference fasta file for the BUSCO loci used in the study.
File: Dataset4_A-genome.zip
Description: This contains the alignment files and tree files for 29 accessions that make up the A-genome group in Glycine.
File: Dataset2_27_accessions.zip
Description: This contains the alignment files, rooted gene trees, and collapsed gene trees (nodes with less than 70% support were collapsed) for 27 accessions of Glycine, plus two outgroups, G. max and Aphicarpeae edgworthii.
File: 570_accessions_100_Busco_concatenated.zip
Description: This contains the alignment of 100 BUSCO loci (the first 5 from each of the 20 chromosomes) and the generated newick file from IQTree2
File: Dataset_5_B-genome.zip
Description: This contains the alignment files and tree files for 29 accessions that make up the B-genome group in Glycine.
File: Dataset1_61_accessions.zip
Description: This contains the alignment files, rooted gene trees, and collapsed gene trees (nodes with less than 70% support were collapsed) for 61 accessions of Glycine, plus two outgroups, G. max and Aphicarpeae edgworthii. The plastome assemblies (minus one inverted repeat) and the alignment of said assemblies using progressiveMauve are also provided, along with the newick tree file from IQTree2 and the dated newick file from RelTime in Mega. The input files for the ABBA-BABA test are also included, along with summary files for each locus and accession across the 61 accessions.
File: Divergence_dates.zip
Description: Files with the 95% confidence interval for each node in the phylogeny for the 61 taxon data for nuclear, plastome, and mitogenome.
Code/software
Alignment files can be viewed with a text editor or a multiple sequence alignment program such as AliView. Tree files can be viewed with a text editor or FigTree
Access information
Other publicly accessible locations of the data:
- Raw data is being submitted to the Sequence Read Archive, and the project number will be made available upon acceptance of the manuscript.
Resequencing reads were mapped to BUSCO loci from Glyince, and SNPs were called. The resulting VCF files were converted to fasta format. Gene trees for each locus were inferred with IQTree2, followed by concatenation of loci and species tree inference with Astral. Gene tree clustering was done with TreeTools in R, unique topologies were inferred with phytools in R, and PhyParts was also used to look at topology support.
