At the heart of the controversy is a technique the authors call SLRI, or spectral label-free Raman imaging, which combines spontaneous confocal hyperspectral Raman microscopy with a deterministic mathematical fitting procedure known as non-negative least squares, or NNLS. The method scans tissue point by point with a single 532-nanometer excitation beam, acquiring a full Raman spectrum at every pixel across thousands of spatial positions and multiple depth planes. Each measured spectrum is then fitted against a library of 42 purified reference compounds, ranging from saturated and unsaturated fatty acids to cytochrome c and other biomolecules, producing spatial maps that show where each reference component contributes most strongly to the measured signal.

The critical letter, according to the authors, misread the scope of these claims. The reply is unambiguous on this point: the published article never asserted that the NNLS-derived maps represent absolute molecular concentrations, nor that all 42 analytes had been definitively and independently identified. Instead, the paper explicitly describes the maps as reference-constrained relative distributions, and the small-molecule images as Raman-inferred, model-constrained representations. The authors lay out a clear analytical hierarchy to make the distinction concrete: a measured Raman spectrum feeds into a deterministic, reference-constrained NNLS coefficient, which in turn generates a spatial coefficient map. No step in the Drosophila study converts those coefficients into milligrams per milliliter or moles per liter. A separate concentration-calibrated implementation exists in the team’s earlier qRamanomics work, published in Cell Reports Methods in 2023, but that calibration was claimed only for selected major biomolecular classes, not for all 42 references.

Where the published text used shorthand such as a component being enriched in a region, the authors explain, the intended meaning is a relatively higher distribution of the corresponding reference-constrained spectral contribution, not a calibrated concentration increase. They also address a related criticism about statistics: violin plots in the supplementary figures quantified only three representative components, unsaturated fatty acids, saturated fatty acids, and oxidized cytochrome c, at the whole-testis level. These were chosen as illustrative examples of signal trends, the authors say, not because the imaging was limited to three molecules. Crucially, they acknowledge that pixels within a single specimen are spatially correlated measurements and were never treated as independent biological replicates, so pixel-level statistical tests were not used as evidence of biological replication.

On the question of whether the computational workflow was reconstructible from the public record, the authors concede the Advanced Science article could have linked the code and study-specific data more explicitly, but they reject the charge that the pipeline was opaque. The same MATLAB qRamanomics processing and NNLS codebase used in the study was previously published and archived on Zenodo, and the principal new analytical input was the study-specific 42-reference library. To close the loop, the team has now released a set of files that make the data lineage explicit: the normalized 42 reference spectra, a representative raw hyperspectral Z-plane containing 30,771 measured spectra on a 263 by 117 grid with 1,114 Raman-shift channels per spectrum, a numerical NNLS coefficient vector for a cytochrome c map, and the corresponding rendered TIFF image. Reshaping the coefficient vector to a 600 by 686 array reproduces the published image with a Pearson correlation of approximately 0.9999, a striking demonstration that the colorful figures trace directly back to raw pixel-resolved spectral measurements rather than to stylized graphics.

Perhaps the most technically substantive portion of the reply concerns identifiability, the question of whether 42 overlapping spectra can be reliably separated at all. The authors report new diagnostics on the exact reference library after interpolation onto the measured Raman axis and application of the fitting window from 400 to 3100 inverse centimeters, excluding the spectroscopically silent region between 1800 and 2700. The 42-by-42 design matrix is full rank, meaning no reference spectrum is an exact linear combination of the others, with 500 fitting variables available for 42 components. Yet the team does not oversell this result. The 2-norm condition number of the unit-area-normalized matrix is 391.5, dropping to 242.6 after additional L2 column normalization, and pairwise correlation and hierarchical clustering analyses reveal substantial collinearity among several chemically related references. The authors explicitly state that these diagnostics are not proof that all 42 components are mutually orthogonal or unconditionally identifiable; rather, they quantify the model dependence that motivated the cautious Raman-inferred terminology in the original paper. They also dismiss a tempting shortcut: refitting a noise-free synthetic mixture with the same design matrix would be mathematically circular and therefore meaningless as validation.

The reply also clarifies sample preparation, which bears directly on whether long acquisitions could have been corrupted by biological motion. The Raman-imaged testes were fixed in 4 percent paraformaldehyde for 20 minutes, washed three times in phosphate-buffered saline, embedded in agarose, and maintained in PBS during acquisition. The specimens were therefore not living tissues undergoing cell migration, differentiation, or metabolite turnover during scans that can nominally run for many hours. The authors further explain that the punctate and ring-like patterns highlighted in the study were identified directly within individual XY planes, comparing measured spectra from peripheral punctate loci and central diffuse loci within the same structure, so variation in z-step size cannot generate the observed patterns. Three-dimensional renderings simply stack consecutively acquired planes recorded in the same microscope coordinate system; no intensity-based registration between separately acquired channels is involved, because every molecular coefficient at a given coordinate derives from the same measured spectrum.

On the modality question, the authors note there is no experimental ambiguity: the study describes spontaneous confocal Raman microscopy with single-beam 532-nanometer excitation, point scanning, full-spectrum acquisition, and EMCCD detection, features that distinguish it from stimulated Raman scattering. The critical letter itself, they observe, reached the same classification, so the objection reduced to a citation-placement issue rather than a claim that the wrong technique was used. The authors also take aim at the conceptual figure accompanying the letter, pointing out that its microscope illustration corresponds to no confocal Raman instrument, its spectral traces are schematic drawings rather than experimentally acquired spectra, and its laser beam is drawn blue when 532-nanometer light is green. If AI-assisted tools generated the figure, they add, appropriate disclosure under journal policies would be important.

The broader stakes of the exchange concern the role of artificial intelligence in spectroscopy. The letter proposed an AI-ready multimodal spatial-biology framework as an alternative, but the authors argue that a conceptual schematic is not evidence of analytical superiority. AI may prove valuable for Raman denoising, representation learning, and multimodal integration, they concede, but it is neither necessary nor evidentiary for validating a deterministic study whose reproducibility rests on fixed inputs, a specified reference matrix, documented preprocessing, and NNLS optimization, all of which yield identical numerical outputs given identical inputs and code. They challenge the correspondent to present one or two core experimental datasets to substantiate the proposed framework.

The authors close by situating the Raman maps within a wider web of orthogonal evidence: two independent Cyp4ae1 RNAi lines consistently produced the spermatogonial differentiation arrest phenotype, while LD540 and TOM20 staining, ATP and NADP+/NADPH assays, and non-targeted metabolomics independently documented lipid accumulation and mitochondrial abnormalities. Bulk assays do not prove microscopic localization of any individual Raman component, the authors acknowledge, but they support the broader metabolic remodeling within which the spatial Raman-inferred patterns are embedded. With the full datasets now posted publicly, including the raw hyperspectral layers, reference library, and diagnostics, the team argues that the study rests on real pixel-resolved data and a documented, deterministic pipeline, and that the appropriate standard for evaluating any spectral imaging claim is a minimum evidentiary framework in which the strength of the interpretation never exceeds the strength of the underlying spectra.

Subject of Research: Label-free Raman imaging methodology and evidentiary standards for spectral map interpretation in biomedical research

Article Title: From Spectral Maps to Biomedical Evidence: A Minimum Evidentiary Framework for Label‐Free Raman Imaging (Reply to Letter‐to‐the‐Editor)

Article References: Li, J., Zhao, X., & Yu, J. (2026). From Spectral Maps to Biomedical Evidence: A Minimum Evidentiary Framework for Label‐Free Raman Imaging (Reply to Letter‐to‐the‐Editor). Advanced Science, Article e77930. https://doi.org/10.1002/advs.77930

Image Credits: AI Generated

DOI: 10.1002/advs.77930

Keywords: Raman spectroscopy, NNLS, label-free imaging, Drosophila, spectral unmixing, reproducibility, qRamanomics, hyperspectral imaging, condition number, biomedical evidence, Cyp4ae1, spatial biology