Track Yourself, Not the Average: Eight-Year Twin Study Points Toward Personalized Aging Medicine

A landmark longitudinal investigation tracking 335 women across an eight-year span has revealed that human biological aging does not follow a uniform, predetermined script. By analyzing whole-blood gene expression alongside serum metabolites over three distinct clinical visits, researchers uncovered that nearly a third of all monitored genes and an array of critical metabolites exhibit marked temporal shifts. Crucially, the trajectories of individual subjects frequently broke away from population averages, moving in diametrically opposed directions. Biological rhythmicity, environmental toxins, and host genetics interact continuously to shape these molecular paths, revealing that systemic decline is highly individualized rather than a synchronized march of chronological time.

Most of what we know about how human blood chemistry changes with age comes from comparing young people with old people at a single moment. That approach cannot separate real aging from the simple fact that people differ from one another. A study in Science takes the harder route: follow the same people and measure them again and again.

The MultiMuTHER project drew blood from 335 female twins aged 32 to 80 at three or more clinic visits between 2009 and 2017, a median of six years from first to last. Each sample was run for RNA sequencing of whole blood and for 915 usable serum metabolites.

The headline numbers are large. Expression of 5,036 genes, 31% of those measured, shifted consistently over time, split almost evenly between rising and falling. Among metabolites, only 45 (5%) showed a shared direction. Androgenic steroids such as DHEA-S fell, while ceramides and sphingomyelins, lipids tied to cardiometabolic risk, rose.

The more interesting result is how poorly the average describes the individual. Even for the genes with the strongest group trends, some women moved the opposite way and others barely moved. A further 136 metabolites, including serotonin and several acylcarnitines, changed substantially within people but in both directions, so the group average came out flat. The inflammatory chemokine gene CXCL9, a key component of a published inflammatory aging clock, behaved the same way.

Immune cells did not age in step either. Using computational estimates of cell-type-specific expression, the team found that CD4 and CD8 T cells mostly lost expression over time, while natural killer cells mostly gained it. The authors read this as remodeling of the immune system, with adaptive immunity waning and innate immunity turning up, although the cell-level data are inferred and not directly measured.

Context mattered a great deal. About 26% of genes and 39% of metabolites varied with the hour of the blood draw, even within office hours. A quarter of both varied with season, with interferon and immune genes peaking in winter. Genes that declined over time were about twice as long as those that rose, supporting earlier reports that aging blunts the transcription of long genes.

The study also caught an environmental policy at work. Serum levels of the “forever chemicals” PFOA and PFOS dropped steadily, in line with regulatory restrictions on their use. At the first visit PFOS tracked with the expression of 234 genes, including a thyroid hormone receptor gene. By the last visit no associations were detectable.

What the paper does not do is tell anyone which trajectory is good or bad. There are no disease or mortality outcomes, no men, and no non-European participants. Its value is as a reference map and as a warning: a single blood omics snapshot, compared against a population average, says less about a person’s aging than commercial test reports imply.

Actionable Insights

This study tested no treatment, so nothing here shows that changing a marker improves health. The practical lessons concern measurement.

First, one-off blood omics tests are noisy. A woman’s gene expression profile correlated only 0.31 with her own profile at another visit, on a scale where 1.0 is identical. That is about 10% shared variation. Metabolite profiles did better at 0.44 for adjacent visits, falling to 0.32 across the full span.

Second, standardize your blood draws. With 26% of genes and 39% of metabolites shifting by time of day, and about a quarter by season, a morning draw in January and an afternoon draw in July are not comparable. Same hour, same fasting state, and ideally the same season.

Third, track yourself, not the average. For 136 of the 181 changing metabolites (75%), people moved in both directions.

Fourth, body weight and blood sugar dominate. The strongest shared molecular pattern shifted about 0.09 standard deviations per BMI point, and women with type 2 diabetes sat roughly one standard deviation from healthy controls, a large gap.

Finally, PFAS blood levels fell after regulation, and their molecular associations disappeared. That is consistent with exposure reduction being worthwhile but does not prove a health benefit.

Context and Source

  • Paywalled Paper: Longitudinal dynamics of gene expression and metabolomics in an aging population cohort, September 2026.
  • Lead institution: Department of Twin Research and Genetic Epidemiology, King’s College London, with the University of Geneva, University of Oxford, and MRC Biostatistics Unit, Cambridge
  • Country: United Kingdom (lead), with Switzerland and Italy
  • Journal: Science, volume 393, issue 6815, 3 September 2026
  • Impact evaluation: The impact score of this journal is 47.3, evaluated against a typical high-end range of 0 to 60+ for top general science, therefore this is an Elite impact journal.

Biomarker Data (Effect Sizes)

The paper reports model coefficients in rank-normalized units, where 1 unit is roughly 1 standard deviation (SD), and does not report absolute concentrations or fold changes. Several figures below are therefore my approximations from the plots and are labeled as such.

How stable is a person’s profile?

  • Gene expression self-correlation between any two visits was 0.31 and did not decay with time. It fell 0.006 per year, and time explained 0.3% of the variation. [Confidence: High]
  • Metabolite self-correlation was 0.44 for adjacent visits and 0.32 from first to last visit. That is an absolute drop of 0.12, or 27% in relative terms, with time explaining 12% of the variation. [Confidence: High]
  • Identical twins resembled each other more than fraternal twins. For gene expression the correlations were about 0.25 versus 0.16 (56% higher), and for metabolites 0.32 versus 0.17 (88% higher).
  • A woman’s expression profile was only modestly closer to her own past self (0.31) than to her identical twin (0.25).

How much changed?

  • Genes with a consistent group trend: 5,036 of 16,292 (31%), with 54% decreasing and 46% increasing.
  • Metabolites with a consistent group trend: 45 of 915 (4.9%), with 24 decreasing and 21 increasing.
  • Features changing within individuals but with no group direction: 136 metabolites (75% of all changing metabolites) but only 25 genes (0.5% of all changing genes).
  • Per-year gene shifts ranged from roughly minus 0.25 to plus 0.30 SD at the extremes (read from the volcano plot). Most significant genes sat well below 0.1 SD per year. [Confidence: Medium]
  • DHEA-S fell by roughly 0.06 SD per year (read from the plot). Over six years that is about 0.35 SD, enough to move a median woman to about the 36th percentile of the baseline distribution. [Confidence: Low to Medium]
  • PFOS and PFOA fell by roughly 0.10 and 0.09 SD per year in the same plot. A second figure plotting levels by calendar year suggests a smaller shift of about 0.2 SD across eight years. The paper does not reconcile these, and no ng/mL values are given. [Confidence: Low]

Gene length

  • Median length of declining genes was 6,602 bp, versus 3,270 bp for rising genes and 4,219 bp for unchanged genes.
  • Declining genes were about twice as long as rising ones and 56% longer than unchanged ones. This is a large and robust separation. [Confidence: High]

Context effects

  • Time of day: 4,265 genes (26%) and 356 metabolites (39%).
  • Season: 4,011 genes (25%) and 220 metabolites (24%).
  • PFOS at first visit: 234 genes (1.4%) and 9 metabolites. PFOA at first visit: 2 genes and 19 metabolites. At last visit: zero for both.

Cross-omic structure

  • 101,171 gene-metabolite associations were found, which is about 0.7% of the roughly 14.9 million pairs tested.
  • Eleven shared factors explained 10% of expression variance and 21% of metabolite variance.
  • Factor 1 shifted 0.09 units per BMI point, so a 5-point BMI difference corresponds to about 0.45 units.
  • For glycemic status, the reported linear contrast of minus 0.69 translates to roughly a 1.0 unit gap between normoglycemic controls and type 2 diabetics. If factor scores have unit variance, that is a Cohen’s d near 1, a large effect. It rests on 11 diabetic women. [Confidence: Low to Medium]

Genetics

  • No genotype-by-time effect on expression passed the 5% false discovery threshold.
  • The 128 “suggestive” genes sit at an 8% false discovery rate, so about 10 are expected to be false.
  • Five variant-metabolite pairs were significant, covering two metabolites: oxidized cysteinylglycine (at DPEP1) and creatinine.

Novelty

  • This is the first study to pair transcriptome and metabolome at three or more time points in the same people over this span, at this sample size.
  • It shows that the transcriptome’s self-similarity holds steady while the metabolome’s decays over years.
  • It quantifies how often individual trajectories contradict the population trend, including for CXCL9.
  • It provides within-person confirmation of the gene-length effect in humans, with NK cells as an apparent exception.
  • It shows PFAS-associated molecular signatures fading as population exposure fell.