The technique belongs to a field called quantitative ultrasound, or QUS. Conventional B-mode ultrasound, the grayscale image familiar from maternity wards and emergency rooms, displays anatomy but says little about what the tissue is actually made of at the microscale. QUS goes a step further by treating the raw backscattered radiofrequency signal as data. When an ultrasound pulse travels into tissue, it scatters off microscopic structures such as cells, vessels, and extracellular fibers. The amplitude of the returning echo, known as the envelope, carries a statistical fingerprint of those scatterers: their density, their spatial organization, and their clustering. By fitting probability distributions to the envelope, researchers can extract numbers that correlate directly with tissue microstructure, numbers that a standard grayscale image simply cannot provide.

The team evaluated three complementary statistical models. The Homodyned-K distribution, the most general of the three, separates the echo into a coherent component, arising from periodically arranged scatterers, and a diffuse component, arising from randomly distributed ones. Its parameter alpha describes scatterer clustering, while a derived parameter k captures the ratio of coherent to diffuse signal power. The Nakagami distribution, a simpler model, uses a shape parameter m that reflects the effective number density of scatterers and a scale parameter Omega related to backscattered energy. The Burr distribution, the newest addition to the envelope-statistics toolkit, assumes that scatterers exhibit multi-scale, power-law behavior, making it sensitive to fractal and elongated structures such as the branching vascular networks and fibrillar collagen that characterize tumors. Its shape parameter b reflects the density of fractal scattering structures, while its scale parameter lambda relates to echo amplitude.

To test whether these statistics could distinguish tumor from healthy tissue, the researchers worked with forty human tissue samples, twenty glioblastoma specimens resected from patients aged 15 to 68 and twenty normal brain samples obtained from the NIMHANS brain bank. Because surgeons cannot safely remove healthy brain tissue during tumor surgery, the control specimens came from post-mortem donations collected within 24 hours. Each sample, fixed in formalin and cut into blocks a few millimeters thick, was embedded in a solidified agarose phantom to hold it steady during imaging. A fraction of every sample was processed for histopathology, and two experienced neuropathologists with 10 and 25 years of experience independently confirmed the diagnoses, providing a ground truth against which the ultrasound statistics could be validated.

Ultrasound data were acquired with a Verasonics Vantage 128 research system driving a linear array transducer at a center frequency of 7.6 megahertz, focused 14 millimeters into each sample. Five independent imaging planes were collected per specimen, with the transducer translated by roughly one millimeter between acquisitions, a distance greater than the lateral beamwidth, ensuring that consecutive measurements sampled genuinely different regions of tissue. Before any statistics were computed, the team measured the frequency-dependent acoustic attenuation of the samples using a broadband pulse-echo insertion-loss technique with a separate 20-megahertz single-element transducer. Because formalin fixation and tissue composition alter how sound is absorbed as it propagates, the raw radiofrequency data were compensated for attenuation in the frequency domain before envelope analysis, preventing absorption artifacts from masquerading as microstructural differences.

The results were striking. Every parameter from all three distributions differed significantly between glioblastoma and normal brain. The HK clustering parameter alpha was higher in tumor tissue, while the coherent-to-diffuse ratio k was lower, indicating increased scatterer clustering together with a greater contribution of diffuse scattering, findings that align with the necrotic zones and microvascular proliferation visible under the microscope. The Nakagami parameters m and Omega were both elevated in tumor samples, reflecting a higher effective density of scatterers and greater backscattered energy, consistent with the dense, pleomorphic tumor cells and high nuclear-to-cytoplasmic ratios of glioblastoma. The Burr parameters b and lambda were also higher in tumor tissue, suggesting that the tumor’s fractal vessel networks and its remodeling of the extracellular matrix, with increased deposition of fibrillar collagens, fibronectin, and laminin, create elongated, anisotropic scattering sites that produce heavy-tailed echo amplitudes. Effect sizes were medium to large for most parameters, with Cohen’s d values exceeding 1.0 for several, indicating robust group separation.

The researchers then fed these features into machine learning classifiers, including logistic regression, linear support vector machines, random forests, and multilayer perceptron neural networks. Performance was assessed at two levels: per-sample, where the five imaging planes were averaged into a single observation per specimen, and per-plane, where each image counted as an independent observation. To avoid optimistic bias, the team used leave-one-out cross-validation for per-sample analysis and group K-fold cross-validation for per-plane analysis, ensuring that planes from the same specimen never appeared in both training and test sets. Hyperparameters were tuned with the Optuna framework, and a stricter nested cross-validation was also performed. Among individual distributions, HK-derived features achieved the highest single-model performance at 90.0 percent per-sample accuracy. Combining all six features pushed the linear SVM to 97.5 percent per-sample accuracy with an area under the receiver operating characteristic curve of 0.985, and a reduced four-feature subset, dropping the scale parameters, achieved comparable results, indicating that shape-related parameters carry most of the discriminative information.

Per-plane classification, a tougher test because individual imaging planes capture localized heterogeneity within a sample, still reached 80.5 percent accuracy with an AUC of 0.854 using the six-feature set. Under nested cross-validation, the strictest evaluation, per-sample accuracy decreased to 90.0 percent, a reduction consistent with the optimistic bias expected when model selection and evaluation share the same data splits, but the overall performance trends across feature subsets remained unchanged. Feature importance analysis using random forests identified the HK clustering parameter alpha and the Nakagami scale parameter Omega as the most discriminative features, underscoring the central role of scatterer organization and backscattered energy in telling tumor from healthy tissue. Notably, increasing model complexity did not consistently improve performance, suggesting that the envelope statistics features are largely separable with simple, interpretable decision boundaries, an encouraging property for any future clinical deployment.

The authors are careful to frame this as proof of concept rather than a ready-made clinical tool. Formalin fixation alters acoustic properties relative to fresh tissue, with previous studies reporting attenuation changes ranging from minus 29 to plus 55 percent after fixation, although because both tumor and control samples were fixed under identical conditions, the relative trends are likely to remain valid. The tumor and control cohorts also came from different sources, with potential confounds in post-mortem interval, handling, age distribution, and sex balance, and the dataset of forty samples was imbalanced toward males. Attenuation compensation relied on group-specific average coefficients rather than sample-specific measurements, which could not be obtained for all tissues due to limited dimensions and insufficient signal-to-noise ratio. The authors also note that whole-region parameter averages do not capture spatial heterogeneity within individual samples, a limitation for margin assessment where local variation is exactly what matters.

Even so, the study opens a genuinely exciting translational path. Because ultrasound is cheap, portable, and already standard in neurosurgical operating rooms, a multiparametric envelope-statistics framework could one day complement conventional B-mode imaging by overlaying quantitative microstructural maps onto the surgeon’s view, helping delineate tumor margins in real time without the expense of intraoperative MRI or the limited sensitivity of fluorescence-guided resection. The researchers envision next steps using freshly excised tissues and, ultimately, in-situ examination of brain tumors in patients with a specialized intraoperative probe. If those studies succeed, the statistical whisper of scattered ultrasound echoes, invisible in a grayscale image, could become a powerful new voice in the fight against one of medicine’s most formidable cancers.

Subject of Research: Science

Article Title: Ultrasound echo envelope statistics to differentiate glioblastoma multiforme from normal brain tissue ex vivo: a proof-of-concept study

Article References: Patil, J. M., Agarwal, M., Rao, S., Vazhayil, V., Pandya, H. J., Mahadevan, A., Shekhar, H., & Mercado-Shekhar, K. P. (2026). Ultrasound echo envelope statistics to differentiate glioblastoma multiforme from normal brain tissue ex vivo: a proof-of-concept study. Frontiers in Acoustics, 4, Article 1882624. https://doi.org/10.3389/facou.2026.1882624

Image Credits: AI Generated

DOI: 10.3389/facou.2026.1882624

Keywords: glioblastoma, quantitative ultrasound, envelope statistics, machine learning, tumor margins, medical imaging, brain tumor, Ultrasound, echo, envelope, statistics, differentiate