Candidate regions under consistent selection in seed beetle (Acanthoscelides obtectus) lines experimentally evolved for early or late reproduction

Since 1986, replicate lines of the seed beetle Acanthoscelides obtectus have been kept under two experimental evolution regimes: Early (E; lines EI–EIV), where beetles reproduce only during the first two days of adult life, and Late (L; lines LV–LVIII), where egg laying is only possible after day 10. After over 200 generations, late lines have evolved a nearly doubled lifespan. This item contains the processed genome-wide data shown in the external genome browser, in five groups: - Selective sweeps: candidate regions under consistent selection in each regime. - Copy-number variants (CNVs): regime-level coverage profiles and CNV calls. - Quantitative trait locus (QTL) mapping: LOD (logarithm of the odds) profiles for lifespan and body weight, a recombination map, and founder coverage from an Early x Late intercross. - Gene annotation and expression: a gene-level annotation, Drosophila orthologs and differentially expressed genes. - Region quality: masks of unreliable regions, GC content and mappability. Groups 1, 2 and 5 are derived from pooled whole-genome sequencing (Pool-seq) of the eight lines. Each line was sequenced as one female and one male pool of 50 individuals, giving 8 pooled samples per regime and 16 in total. Code related to the analyses is available from the GitHub repository: https://github.com/Goran-Arnqvist/UUBeetleLab (https://github.com/Goran-Arnqvist/UUBeetleLab) 1. Selective sweepsCandidate regions with consistent signatures of a selective sweep across populations: 574 regions in the Early regime and 458 in the Late regime. - Sweep detection per sample. Read pileups were generated for each of the 16 pooled samples with SAMtools v0.1.19, and each sample was analysed independently with Pool-HMM v1.4.4. Pool-HMM identifies regions in which the allele frequency spectrum (AFS) is skewed towards extreme allele frequencies, as expected after a selective sweep, and uses a hidden Markov model to call candidate loci. Parameters: initial theta 0.005, 100 chromosomes per pool, mutation-rate parameter k = 1e-7, every tenth polymorphic position analysed. The X chromosome (CAVLJG010000002.1) was not analysed. - Consistency within each regime. The 8 per-sample call sets of a regime (4 lines x 2 sexes) were intersected with bedtools multiinter (BEDTools v2.31). A region was kept if it was called in at least 6 of the 8 samples. Requiring support across independently evolving replicate lines prioritises regions associated with the selection regime over signals specific to one line. 2. Copy Number Variants CNVs were called by depth of coverage. Standard CNV callers assume a fixed ploidy, which pooled samples do not have, so the two regimes were used as mutual controls: a CNV is only called where coverage differs consistently between all Early and all Late lines. - Coverage. Read depth per non-overlapping 100 base pair (bp) bin was computed for each pooled sample with mosdepth v0.3.3, excluding reads with mapping quality below 20 or with any SAM (Sequence Alignment/Map) flag in 3852. - Filtering and bias correction. Bins were removed if they overlapped the coverage masks of group 5, had GC content below the 5th or above the 99th percentile, or had mappability below 0.10. Depth was corrected for GC and then mappability bias with weighted LOESS (locally estimated scatterplot smoothing) regression. - Regime contrast. Female and male depths were summed per line, and the median across the four lines of each regime was taken per bin. The log2(Late / Early) ratio of the medians was denoised with a stationary wavelet transform (Haar wavelet, 5 levels, SureShrink soft thresholds). - Segmentation. The denoised ratio was segmented with circular binary segmentation (CBS; R package PSCBS v0.68.0, undo threshold 4 standard deviations). - CNV calls. Segments shorter than 1 kb or with more than 60% missing bins were removed. A segment was kept if the fold change between regimes was at least 1.5 and a permutation test over the 8 lines gave P = 0.95, coverage >= 0.80) and intersected with protein-coding genes. The score column is the log2 fold change (positive = higher expression in Early lines); pvalue_bh_fdr is the false discovery rate; color_fold_change is a display colour, mapped directly from the original log2 fold change to a colourmap. The protein accession of a gene (e.g. CAK1619843.1) links the three file types: it is the Name of protein-coding genes in the annotation, the ID in the ortholog file and the geneID in the DEG files. 5. Region qualityMasks of unreliable regions built for this repeat-rich genome (over 60% repeats), used to filter the sweep, CNV and Gene Ontology analyses. Unlike the other files, these cover all 3,796 sequences of the assembly. - Coverage masks. Read depth of the 16 pooled samples in 1,000 bp windows (100 bp step); regions above the 99th or below the 7th percentile, consistently across samples of both regimes. - Mapping-quality mask. Regions where fewer than 40% of reads have mapping quality >= 10. - Assembly mask. Regions within 150 bp of an assembly gap, and contigs shorter than 20 kb. - GC content. Fraction (0–1, not a percentage) of G and C bases per 100 bp bin. - Mappability. Paired-end reads simulated from the reference (dicey v0.3.3, insert size 400 bp) and mapped back with BWA-MEM2 v2.2.1. Values are simulated-read depth (0–300); divide by 300 for mappability on a 0–1 scale. FilesEach .gz file has a tabix index with the same name plus .tbi (not listed below). The complete file list and file sizes are provided in MANIFEST.txt. 1. Selective sweeps - selective_sweeps_consistent_Early.bed.gz: Candidate sweep regions, Early; BED3, no header - selective_sweeps_consistent_Late.bed.gz: Candidate sweep regions, Late; BED3, no header 2. Copy-number variants - coverage_median_corrected_100bp_{Early,Late}.bw: Median across {Early,Late} lines of the summed female + male corrected depth - coverage_log2_ratio_Late_vs_Early_100bp.bw: log2(Late / Early); positive = higher coverage in Late - coverage_log2_ratio_Late_vs_Early_100bp_denoised.bw : As above, wavelet-denoised (used for segmentation) - CNV_segmentation_log2_ratio_Late_vs_Early.bw: All CBS segments (mean log2 ratio), before filtering - CNV_calls_{Early,Late}.bed.gz: Final CNV calls, Early (3,001), Late (3,779). Columns: chrom, start, end, call (deletion/gain), cn_deviation (population copy number minus 2) - CNV_calls_{Early, Late}.bw: Early/Late calls as bigWig (value = cn_deviation) 3. QTL mapping, recombination and founder coverage Files ending in _markers.bw hold the original LOD values at marker positions. Each has a display-only version with the same name ending in _100kb_interpolated.bw (not listed). - QTL_LOD_lifespan_IM_markers.bw: Lifespan, interval mapping - QTL_LOD_bodyweight_IM_markers.bw: Body weight, interval mapping - QTL_LOD_lifespan_MQM_markers.bw: Lifespan, MQM; all six QTLs - QTL_LOD_bodyweight_MQM_chr{2-22.9cM,4-25.3cM,4-28.3cM,6-12.6cM,9-90.5cM}_markers.bw: Body-weight MQM LOD profiles for QTLs on chromosomes 2, 4, 6, and 9 at the indicated cM positions. - recombination_rate_cM_per_Mb_100kb.bw: Local recombination rate per 100 kilobase (kb), unsmoothed - recombination_rate_cM_per_Mb_2Mb_smoothed.bw: As above, 2 Megabase (Mb) moving average - F0_coverage_mean_100bp_{E2,L5}.bw: Mean normalized depth of the four {E2,L5} founders - F0_coverage_log2_ratio_L5_vs_E2_100bp.bw: log2(L5 / E2); positive = higher in L5; capped at -5 to 5 4. Gene annotation and expression - Aobtectus_GCA_963669975.1_genes_eggNOG.gff3.gz: One feature per gene, with eggNOG-mapper annotation - Aobtectus_GCA_963669975.1_Drosophila_orthologs.gff3.gz: D. melanogaster orthologs - DEGs_Early_vs_Late_*.gff3.gz: Differentially expressed genes (DEGs) between Early and Late lines, for female abdomen, female head + thorax, male abdomen, and male head + thorax. 5. Region quality - excluded_regions_high_coverage.bed.gz: Abnormally high coverage; BED3, no header - excluded_regions_low_coverage.bed.gz Abnormally low coverage; BED3, no header - excluded_regions_low_mapping_quality.bed.gz: Low mapping quality; BED3, no header - excluded_regions_assembly_gaps_small_contigs.bed.gz: Near assembly gaps, or small contigs; BED3, no header - GC_fraction_100bp.bw: GC content per 100 bp bin, as a fraction (0–1) - mappability_100bp.bw: Simulated-read depth per 100 bp bin (0–300) Item thumbnail adapted from: Udo Schmidt, CC BY-SA 2.0, via Wikimedia Commons.