Data from: Heterogeneous genomic architecture of skeletal armour traits in sticklebacks
Data files
Jul 31, 2026 version files 2.12 GB
-
Left_side_backcross.zip
1.09 GB
-
qtl_inputs.zip
1.07 MB
-
README.md
3.39 KB
-
Right_side_backcross.zip
1.03 GB
-
TableS1_sample_information_n221.csv
21.65 KB
Abstract
Whether populations adapt to similar selection pressures using the same underlying genetic variants depends on population history and the distribution of standing genetic variation at the metapopulation level. Studies of sticklebacks provide a case in point: when colonising and adapting to freshwater habitats, three-spined sticklebacks (Gasterosteus aculeatus) with high gene flow tend to fix the same adaptive alleles in the same major loci, whereas nine-spined sticklebacks (Pungitius pungitius) with limited gene flow tend to utilize a more heterogeneous set of loci. In accordance with this, we report results of quantitative trait locus (QTL) analyses using a backcross design showing that lateral plate number variation in the western European nine-spined sticklebacks mapped to three moderate-effect QTL, contrary to the major-effect QTL in three-spined sticklebacks and different from the four QTL previously identified in the eastern European nine-spined sticklebacks. Furthermore, several QTL were identified associated with variation in lateral plate size, and three moderate-effect QTL with body size. Together, these findings indicate more heterogenous and polygenic genetic underpinnings of skeletal armour variation in nine-spined than three-spined sticklebacks, indicating limited genetic parallelism underlying armour trait evolution in the family Gasterostidae.
https://doi.org/10.5061/dryad.tmpg4f54x
This Dryad archive includes data generated in the QTL study of the lateral plate variation in nine-spined sticklebacks. A F2 backcross was generated using stickleback lineages having different numbers of lateral plates: the fully-plated Pungitius sinensis and the partially-plated P. pungitius. A total of 221 individuals were prepared into RADseq libraries using the PstI restriction enzyme and sequenced using HiSeqTM 2000 (BGI Hong Kong). The F2 progeny were phenotypically sexed and measured in lateral plate number, lateral plate size, and body size. QTL analyses were conducted to find genetic underpinnings of the putatively adaptive traits (lateral plate variation) and body size.
Description of the data and file structure
Left_side_backcross.zip and Right_side_backcross.zip
These files are the raw digital photographs used for phenotypic measurements. Each F2 individual has two photographs taken on the left and right body sides. Photos are grouped by clutches and named by "clutch ID" - "individual ID", consistent with the IDs in the sample information table.
TableS1_sample_information_n221.csv
This is the sample information of the 221 individuals (also in the supplementary materials of this article), including the raw phenotypic data. The 12 individuals with potentially misidentified or unknown phenotypic sexes were marked. The columns are:
- the full sample ID used in analyses (SampleID)
- The pedigree of the sample (pedigree)
- the individual ID (ID)
- the clutch of the sample (Clutch)
- the identified phenotypic sex (Phenotypic_Sex)
- the measured body length in mm (BodyLength_Stained.fish)
- the measured body height in mm (BodyHeight_Stained.fish)
- the number of myomere having the largest plate on the left (Left_Plate1_Myomere)
- the area size of the largest left plate in square mm (Left_Plate1_Area)
- the width of the largest left plate in mm (Left_Plate1_Width)
- the height of the largest left plate in mm (Left_Plate1_Height)
- the number of myomere having the largest plate on the right (Right_Plate1_Myomere)
- the area size of the largest right plate in square mm (Right_Plate1_Area)
- the width of the largest right plate in mm (Right_Plate1_Width)
- the height of the largest right plate in mm (Right_Plate1_Height)
- the total number of plates on the left side (PlateN_Left)
- the total number of plates on the right side (PlateN_Right)
- the percent of plated myomere on the left (PlateMyom.percent_Left)
- the percent of plated myomere on the right (PlateMyom.percent_Right)
- individuals removed from analyses due to ambiguous sex ID (remove_for_ambiguous_phenotypic_sex)
qtl_inputs.zip
Input files for the QTL analyses, including: the pedigree, the transformed phenotypes (qtl_pheno_203_final.csv), the formatted genotypes (plateQTL_geno_203_SP.csv), the genetic map (plateQTL_gmap.csv), the physical map (plateQTL_pmap.csv), and the R/qtl2 control file (PlateQTL_203_final.yaml).
Additional information
R scripts for these analyses are available on Github and Zenodo (https://github.com/xuelingyi/QTL_plate; DOI: 10.5281/zenodo.1155930).
To identify genomic regions underlying lateral plate variation in Pungitius sticklebacks, a backcross experiment was conducted between the partially plated western European (WL) nine-spined stickleback (Pungitius pungitius) and the fully plated Pungitius sinensis. Phenotypes were measured in ImageJ and genotypes were obtained from the reduced-representation sequencing (RADseq). Linkage maps were constructed using LepMAP3, and quantitative trait loci (QTL) analyses were conducted using R/qtl2.
