Contributions
Harvard Forest Long-Term Warming — DNA vs RNA Functional Response analyzed a 25-year +5°C soil-warming experiment at Harvard Forest Barre Woods (NMDC nmdc:sty-11-8ws97026; 42 biosamples) 1. It found a real but modest community-level response that reproduced the published Actinobacteria-up/Acidobacteria-down signal 1. In organic soil, Actinobacteria rose from 0.249 to 0.315 (log2 FC +0.341, q=0.049) and Acidobacteria fell from 0.035 to 0.024 (log2 FC −0.549, q=0.049) 1. Mineral soil showed no phylum-level changes after false-discovery-rate (FDR) correction 1. PERMANOVA (permutational multivariate analysis of variance) on Bray–Curtis genus distances attributed 7.6% of variance to treatment (p=0.069), 30.6% to horizon (p=0.0002), and 41% to the treatment × horizon factor (p=0.0002) 1. The paired n=25 subset confounded horizon with incubation, because every organic sample was incubated and every mineral sample was direct 1. Once that confound was removed, direct-sample analyses gave comparable DNA and RNA treatment R² (10-13%), so the hypothesis that RNA composition shifts more than DNA under warming was not supported 1. The report notes that this is in tension with the larger transcript-level effect reported by Roy Chowdhury et al. (2021) at the same site 1. It also notes that differences in sampling timepoint, replication and analytical pipeline complicate the comparison 1.
The same project tested a curated 62-KO (KEGG Orthology gene family) carbon-cycling list for enrichment among heated-up KOs 1. The list was enriched among heated-up DNA KOs in organic soil: 5 of 57 hits versus 428 of 12,806, Fisher odds ratio 2.78, one-sided p=0.042 1. There was no enrichment in the other horizon × pool combinations, and no individual RNA KO survived FDR (q<0.10) 1. In the RNA pool, pmoA (K10944) and pmoB (K10945) were directionally positive in both horizons; for example, pmoA had log2 FC +0.730 in organic soil (p=0.009) and +0.743 in mineral soil (p=0.029) 1. Glyoxylate-cycle genes aceA/icl and aceB/glcB increased in heated mineral RNA (p=0.037 each) 1. All of these signals were nominal and none survived FDR across 14K KOs 1. Heated mineral samples had fewer detectable metabolites (155 ± 4 versus 167 ± 9 ChEBI identifications; Mann–Whitney p=0.012), while organic samples showed the same direction without significance (p=0.209) 1. No individual metabolite passed Benjamini–Hochberg false-discovery-rate (BH-FDR) correction 1. The report proposes compositional dilution, substrate depletion, and compositional turnover over 25 years as explanations for the RNA-pool results 1. It presents all three as author-proposed explanations rather than tested results 1.
SNIPE Defense System in the KBase Data Lakehouse Pangenome surveyed the SNIPE phage-defense system across the KBase Data Lakehouse pangenome 2. It found DUF4041/PF13250 in 4,572 gene clusters across 1,696 species and 33 bacterial and archaeal phyla, which exceeds the >500 homologues reported in the original paper 2. Of these clusters, 13.3% (604) were core, 30.7% were accessory, and 56.1% (2,566) were singleton 2. The report reads this as consistent with mobile defense-island carriage, but mobility is inferred from pangenome category composition rather than measured directly 2. The project identified PF13455 (Mug113) as the SNIPE nuclease family, distinct from canonical GIY-YIG PF01541 2. It found zero gene clusters containing both DUF4041/PF13250 and PF01541 2. The report records two conflicting species counts for IPR025280: 2,466 in its domain table and 1,373 in its Methods 2. The source gives no way to reconcile them 2. Among species with AlphaEarth embeddings (n=1,069 SNIPE-bearing versus n=13,977 non-SNIPE), 22 of 64 dimensions differed significantly under Welch’s t-tests with Bonferroni correction 2. The largest effect was dimension A19, with Cohen’s d = 0.26 and p = 5.0e-14, and effect sizes of |d| = 0.11–0.26 indicated a modest environmental shift rather than strict habitat segregation 2. AlphaEarth coverage was 28.4% of genomes and was biased toward environmental isolates with latitude/longitude metadata 2.
Using RB-TnSeq (random barcode transposon sequencing) data from the Fitness Browser, the SNIPE project found the strongest E. coli K-12 ManXYZ fitness defects on D-glucosamine and D-mannose 2. On D-mannose, manX, manY and manZ averaged -2.75, -3.00 and -2.74 2. The report reads these results as contradicting the UniProt “fructose-specific” annotation for ManX 2. It frames SNIPE as giving phage resistance without the metabolic cost of losing man transporter genes 2. It offers no direct comparative fitness measurement, so this remains a mechanistic proposal 2. The project also identified one full two-domain SNIPE protein in Methanococcus maripaludis JJ (MMJJ_RS01635) 2. Its minimum fitness was -1.16 on formate/acetate, and it was dispensable under most conditions 2.