Test Your Hypothesis

 Quick Analysis Plan Copied

 Long Analysis Plan

Evidence Base

The ~12 My satellite identity half-life is a premise from centromere repeat evolution; after 2 Mya divergence, background sequence decay alone would leave ~89% identity between G. urbanum and G. rivale satellites (0.8908987181403393). The atlas data show this prediction is vastly exceeded: G. urbanum centromeres are 87.5% Gypsy-Athila (~18 Mb) while G. rivale retains 65.7% satellite drGeuRiva.159 (~57 Mb) at syntenic locations Geum sister-species centromere contrast. Additionally, LTR elements inside satellite arrays have higher 5′/3′ identity than those outside (Welch t-test, P<1.55×10⁻¹²), indicating centrophilic insertion is recent and active Centrophilic LTR age evidence.

The ~3× centromere size contraction (57 Mb → 18 Mb) alongside complete satellite replacement cannot be explained by passive sequence decay at a 12 My half-life in 2 My. The ~17% interruption fraction with 70% transposon content across satellite arrays genome-wide further suggests that even in α-state genomes, transposons are actively colonizing centromeric space Recurrent centromere state transitions.

Bioinformatics Analysis Plan

Step 1 – Quantify Athila invasion rate: Extract G. urbanum Athila elements within centromeric boundaries; estimate insertion dates from LTR 5′/3′ divergence using species-specific mutation rate. Determine whether insertion bursts post-date the 2 Mya split.

Step 2 – Model satellite decay vs displacement: Compare observed drGeuRiva.159 retention in G. urbanum (residual copies) against the passive decay expectation of ~89% identity retention. Significant excess loss (array collapse, not just sequence divergence) indicates active displacement.

Step 3 – Test nesting structure: Identify Athila insertions nested inside residual satellite monomers (target-site duplications, insertion breakpoints) as direct evidence of transposon-mediated array fragmentation.

Step 4 – Syntenic reconstruction: Use SyRI/MUMmer alignments to map where G. rivale's ~57 Mb satellite arrays correspond to G. urbanum's ~18 Mb Athila blocks; quantify deletion breakpoints.

Step 5 – Ancestral state timing: Apply chronos-dated phylogenetic branch lengths to constrain whether the α→Ω transition occurred within the 2 Mya window or reflects an older polymorphism.

Conclusion (inferred): The acceleration is plausible but the atlas study performed no direct causal tests (no transposon silencing experiments, no CENH3 binding assays), so the mechanism remains correlative. BGPT infers Athila invasion as the most parsimonious explanation for the observed rapid replacement, but formal falsification requires testing whether Athila insertion causes satellite array contraction (e.g., via recombination suppression at insertion sites) versus a neutral hitchhiking model Authors acknowledge no causal tests.