Paper Review — Verify Claims Fast

Study aim and novelty: The authors performed paired transcriptome (RNA‑seq) and whole‑genome bisulfite sequencing (WGBS) of the goat pineal gland across prepuberty and puberty and report 2,703 DEGs and 1,294 DMRs, claiming methylation changes may regulate puberty-related pineal transcriptional programs — the first such WGBS+RNA‑seq dataset for goat pineal across puberty (DOI:10.1038/s41598-024-84559-x) This paper (core data)

Strengths — data & methods:

• Paired WGBS (30× mean coverage reported) and RNA‑seq on the same biological timepoints provides direct methylation–expression overlap capability (paper reports ~270 million clean WGBS reads/group and RNA‑seq depth ~72.96M reads/sample) Sequencing depth and QC metrics
• Data publicly deposited (GEO: GSE234298 for RNA-seq, GSE234928 for RRBS/WGBS), enabling re-analysis and reproducibility checks (paper provides accessions) Data availability statement

Limitations & critical weaknesses:

1. Pooling and replication for WGBS: WGBS input was pooled from pairs to produce 1 prepubertal pool (3 pooled samples?) and 1 pubertal pool per group as described; pooling reduces sample-level variance estimates and effectively lowers biological replication for methylation analysis, inflating risk of false positives in DMR calling and making per-animal inference impossible. The methods state samples were paired into 3 pre and 3 pub samples for RNA and pooled for methylome analysis (see Methods). This is a major reproducibility blindspot for epigenetic studies where between-animal variance is large Pooling reduces WGBS replication
2. Statistical thresholds and multiple testing: The paper uses DESeq2 with adjusted P < 0.05 for DEGs (appropriate) but the DMR pipeline relies on swDMR with sliding windows (1 kb windows, step 100 bp) and additional filters; the exact per-region FDR control and false discovery estimates for 1294 DMRs are not fully transparent in the main text (details in supplemental), which impedes independent assessment of DMR reliability. Clearer reporting of per-DMR FDR and the number of CpGs per DMR would improve confidence (paper references swDMR and hierarchical testing approach) DMR calling pipeline description
3. Interpretation: correlation vs causation — authors state methylation "regulated" expression and point to specific genes (e.g., ATP5F1D, CACNB2, PTEN). However, most overlaps between DMRs and DEGs are correlative; directionality is not established and the manuscript lacks targeted validation (e.g., locus‑specific bisulfite validation with multiple animals, chromatin accessibility, or reporter assays). The RT‑qPCR validation confirms expression trends but not causal methylation control RT-qPCR validated expression not causality
4. Biological model & sample size: Study used Wanlin White goats (n=12 animals total) with narrow age windows (2.5–3.0 mo pre, 4.0–4.5 mo pub) and pooled methylome samples; this provides a useful developmental snapshot but may not generalize across breeds or sexes, and the age difference includes large body-weight and physiological changes (weight difference reported ~8 kg vs ~20 kg) that could confound methylation/expression patterns (metabolic pathway enrichment could reflect growth/weight effects) Sample characteristics and potential confounding

Biological plausibility and candidate genes:

• LIN28B — increased in puberty in the pineal gland; LIN28B is a well established regulator of pubertal timing across mammals (GWAS and experimental studies). The paper's finding that LIN28B mRNA increases in pubertal goat pineal is consistent with prior work linking Lin28 family to puberty timing (but cell-type specific roles in pineal are novel and need mechanistic follow-up) Paper LIN28B result & discussion
• ATP5F1D, CACNB2, PTEN — identified as network hubs in STRING and as overlapping DMR/DEG genes. These genes plausibly connect to metabolism (ATP5F1D), calcium signaling (CACNB2), and PI3K/AKT metabolic signaling (PTEN) — consistent with enrichment in OXPHOS/metabolic pathways; again, functional validation is lacking Hub genes identified

Blind spots & alternate explanations:

• Metabolic and mitochondrial pathway enrichment can reflect systemic growth and changing cell composition in the pineal gland between 2.5–3.0 mo and 4.0–4.5 mo, not necessarily direct drivers of puberty timing.
• Pooling and limited WGBS replication prevent estimation of inter-animal variability; reported DMRs may be sensitive to pooling strategy.
• Cell-type heterogeneity: pineal gland contains pinealocytes, astrocytes, vascular cells — bulk RNA/WGBS mixes signals. Single-cell or sorted-cell assays would help define cell-type specific methylation/expression changes.

Opportunities — next experiments that would strengthen claims:

1. Locus‑specific bisulfite pyrosequencing (or targeted bisulfite sequencing) across multiple independent animals (not pooled) for the top DMRs (e.g., promoters/first exons near LIN28B, ATP5F1D, CACNB2, PTEN) to confirm reproducibility and inter-individual variance.
2. Chromatin accessibility (ATAC‑seq) and H3K27ac ChIP in pineal gland at the two stages to test if DMRs overlap active regulatory elements.
3. Cell‑type deconvolution (single‑nucleus RNA‑seq / snATAC) to localize changes to pinealocyte vs glial/vascular compartments.
4. Functional perturbation: e.g., demethylation (TET overexpression) or targeted epi-editing (dCas9‑TET/ DNMT3A) of candidate promoter DMRs in ex vivo pineal cultures to test expression effects.

Reproducibility & transparency: Data are deposited in GEO (GSE234298, GSE234928) enabling re-analysis; authors report software versions (DESeq2 v1.20.0, Hisat2 v2.0.5, Bismark v0.12.5, etc.) which aids reproducibility. However, the methylome pooling strategy and incomplete reporting of DMR counts per region limit independent statistical verification without reprocessing raw FASTQ (which is available) Reproducibility and methods reporting

Confidence in principal claims: Moderate — the sequencing depth and gene-level DE calls are credible; WGBS summary metrics are plausible. But causal links between DNA methylation and expression for specific genes are not demonstrated; pooling reduces confidence in the methylome results as representative of group-level variability.

Transparency note / potential biases considered: I flagged sample pooling, small N, and potential metabolic confounding; funding sources are academic (National Natural Science Foundation of China, etc.) and authors declare no competing interests — no direct COI detected in the manuscript's statements. All claims above are tied to the published paper (DOI provided) and its supplementary materials.