What the paper shows. Using 130 curated mammalian sequences (β₯39 species per paralog) from GenBank and Ensembl Compara, the authors analyzed TET1, TET2 and TET3 separately after GARD found no recombination breakpoints (KH test, p = 0.01). Codon-based estimates under model M0 indicate pervasive purifying selection, with a strikingly lower Ο for TET3 (0.115) than TET1 (0.385) or TET2 (0.346). Under the PAML M8 model, TET1 (3.8% of sites, 2ΞlnL = 21.22, p < 0.001) and TET2 (2.2%, 2ΞlnL = 41.27, p < 0.001) showed significant positive selection, while all site-model comparisons for TET3 were non-significant or negative. To reduce false positives, sites were called only when at least three of six methods (PAML M8 plus Datamonkey FEL, FUBAR, MEME, REL, SLAC) agreed, yielding 34 codons in TET1, 29 in TET2, and 13 in TET3. No positively selected sites occur in the CXXC DNA-binding domain, consistent with strong constraint there .
Branch-site and structural evidence. Using Guindon's stochastic branch-site framework (Fitmodel, M3+S2 best-fitting for all paralogs), 25% of TET1 and 17% of TET2 sites were assigned to an Ο3 class >1 (1.34 and 1.54 respectively), whereas TET3's Ο3 was only 1.07 with 9% of sites β interpreted as relaxed rather than positive selection . Mapping five positively selected TET2 catalytic-core residues (1203S, 1204S, 1205S, 1892N, 1312L) onto the human TET2βDNA crystal structure shows all lie on the protein's outer surface, none overlapping catalytic residues β supporting the authors' speculation that selection acted on interaction interfaces rather than enzymatic chemistry . The functional importance of TET constraint is corroborated by knockout work showing TET-mediated demethylation is essential for reprogramming .
Critical caveats. (1) No functional assays test any of the inferred sites β the surface-exposure argument is plausible but unverified; the authors themselves concede the cancer-variant overlap (TET2 positions 123, 817, 1204, 1762, 1811 in myeloproliferative neoplasms) is not established as statistically significant . (2) GARD's null result does not guarantee absence of recombination, which can confound site-model inference. (3) Sampling mirrors genome availability (mammalian-biased), and only coding regions were analyzed, excluding regulatory evolution. (4) Ο values barely above 1 and modest site proportions mean the 'adaptive' signal is moderate, not dramatic. What would falsify the central claim: stricter models (e.g., BUSTED/MEME re-analysis with updated alignments), or mutagenesis showing the five surface residues do not affect TET2 interactions.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.