Insulin Resistance Increases TNBC Aggressiveness and Brain Metastasis via Adipocyte-derived Exosomes.

Patients with triple negative breast cancer (TNBC) and comorbid Type 2 Diabetes (T2D), characterized by insulin resistance of adipose tissue, have higher risk of metastasis and shorter survival. Adipocytes are the main non-malignant cells of the breast tumor microenvironment (TME). However, adipocyte metabolism is usually ignored in oncology and mechanisms that couple T2D to TNBC outcomes are poorly understood. Here we hypothesized that exosomes, small vesicles secreted by TME breast adipocytes, drive epithelial-to-mesenchymal transition (EMT) and metastasis in TNBC via miRNAs. Exosomes were purified from conditioned media of 3T3-L1 mature adipocytes, either insulin-sensitive (IS) or insulin-resistant (IR). Murine 4T1 cells, a TNBC model, were treated with exosomes in vitro (72h). EMT, proliferation and angiogenesis were elevated in IR vs. control and IS. Brain metastases showed more mesenchymal morphology and EMT enrichment in the IR group. MiR-145a-3p is highly differentially expressed between IS and IR, and potentially regulates metastasis. Significance: IR adipocyte exosomes modify TME, increase EMT and promote metastasis to distant organs, likely through miRNA pathways. We suggest metabolic diseases such as T2D reshape the TME, promoting metastasis and decreasing survival. Therefore, TNBC patients with T2D should be closely monitored for metastasis, with metabolic medications considered.

Learning Gaussian Graphical Models from Correlated Data.

Gaussian Graphical Models (GGM) have been widely used in biomedical research to explore complex relationships between many variables. There are well established procedures to build GGMs from a sample of independent and identical distributed observations. However, many studies include clustered and longitudinal data that result in correlated observations and ignoring this correlation among observations can lead to inflated Type I error. In this paper, we propose a Bootstrap algorithm to infer GGM from correlated data. We use extensive simulations of correlated data from family-based studies to show that the Bootstrap method does not inflate the Type I error while retaining statistical power compared to alternative solutions. We apply our method to learn the GGM that represents complex relations between 47 Polygenic Risk Scores generated using genome-wide genotype data from a family-based study known as the Long Life Family Study. By comparing it to the conventional methods that ignore within-cluster correlation, we show that our method controls the Type I error well in this real example.

Cross-platform proteomics signatures of extreme old age.

In previous work we used a Somalogic platform targeting approximately 5000 proteins to generate a serum protein signature of centenarians that we validated in independent studies that used the same technology. We set here to validate and possibly expand the results by profiling the serum proteome of a subset of individuals included in the original study using liquid chromatography tandem mass spectrometry (LC-MS/MS). Following pre-processing, the LC-MS/MS data provided quantification of 398 proteins, with only 266 proteins shared by both platforms. At 1% FDR statistical significance threshold, the analysis of LC-MS/MS data detected 44 proteins associated with extreme old age, including 23 of the original analysis. To identify proteins for which associations between expression and extreme-old age were conserved across platforms, we performed inter-study conservation testing of the 266 proteins quantified by both platforms using a method that accounts for the correlation between the results. From these tests, a total of 80 proteins reached 5% FDR statistical significance, and 26 of these proteins had concordant pattern of gene expression in whole blood. This signature of 80 proteins points to blood coagulation, IGF signaling, extracellular matrix (ECM) organization, and complement cascade as important pathways whose protein level changes provide evidence for age-related adjustments that distinguish centenarians from younger individuals.

The relationship between 11 different polygenic longevity scores, parental lifespan, and disease diagnosis in the UK Biobank.

Large-scale genome-wide association studies (GWAS) strongly suggest that most traits and diseases have a polygenic component. This observation has motivated the development of disease-specific "polygenic scores (PGS)" that are weighted sums of the effects of disease-associated variants identified from GWAS that correlate with an individual's likelihood of expressing a specific phenotype. Although most GWAS have been pursued on disease traits, leading to the creation of refined "Polygenic Risk Scores" (PRS) that quantify risk to diseases, many GWAS have also been pursued on extreme human longevity, general fitness, health span, and other health-positive traits. These GWAS have discovered many genetic variants seemingly protective from disease and are often different from disease-associated variants (i.e., they are not just alternative alleles at disease-associated loci) and suggest that many health-positive traits also have a polygenic basis. This observation has led to an interest in "polygenic longevity scores (PLS)" that quantify the "risk" or genetic predisposition of an individual towards health. We derived 11 different PLS from 4 different available GWAS on lifespan and then investigated the properties of these PLS using data from the UK Biobank (UKB). Tests of association between the PLS and population structure, parental lifespan, and several cancerous and non-cancerous diseases, including death from COVID-19, were performed. Based on the results of our analyses, we argue that PLS are made up of variants not only robustly associated with parental lifespan, but that also contribute to the genetic architecture of disease susceptibility, morbidity, and mortality.

yQTL Pipeline: a structured computational workflow for large scale quantitative trait loci discovery and downstream visualization.

Quantitative trait loci (QTL) denote regions of DNA whose variation is associated with variations in quantitative traits. QTL discovery is a powerful approach to understand how changes in molecular and clinical phenotypes may be related to DNA sequence changes. However, QTL discovery analysis encompasses multiple analytical steps and the processing of multiple input files, which can be laborious, error prone, and hard to reproduce if performed manually. In order to facilitate and automate large-scale QTL analysis, we developed the yQTL Pipeline, where the 'y' indicates the dependent quantitative variable being modeled. Prior to genome-wide association test, the pipeline supports the calculation or the direct input of pre-defined genome-wide principal components and genetic relationship matrix when applicable. User-specified covariates can also be provided. Depending on whether familial relatedness exists among the subjects, genome-wide association tests will be performed using either a linear mixed-effect model or a linear model. Using the workflow management tool Nextflow, the pipeline parallelizes the analysis steps to optimize run-time and ensure results reproducibility. In addition, a user-friendly R Shiny App is developed to facilitate result visualization. Upon uploading the result file, it can generate Manhattan plots of user-selected phenotype traits and trait-QTL connection networks based on user-specified p-value thresholds. We applied the yQTL Pipeline to analyze metabolomics profiles of blood serum from the New England Centenarians Study (NECS) participants. A total of 9.1M SNPs and 1,052 metabolites across 194 participants were analyzed. Using a p-value cutoff 5e-8, we found 14,983 mQTLs cumulatively associated with 312 metabolites. The built-in parallelization of our pipeline reduced the run time from ~90 min to ~26 min. Visualization using the R Shiny App revealed multiple mQTLs shared across multiple metabolites. The yQTL Pipeline is available with documentation on GitHub at https://github.com/montilab/yQTL-Pipeline.

Total RNA sequencing reveals gene expression and microbial alterations shared by oral pre-malignant lesions and cancer.

Head and neck cancers are a complex malignancy comprising multiple anatomical sites, with cancer of the oral cavity ranking among the deadliest and the most disfiguring cancers globally. Oral cancer (OC) constitutes a subset of head and neck cancer cases, presenting primarily as tobacco- and alcohol-associated oral squamous cell carcinoma (OSCC), with a 5-year survival rate of ~ 65%, partly due to the lack of early detection and effective treatments. OSCC arises from premalignant lesions (PMLs) in the oral cavity through a multi-step series of clinical and histopathological stages, including varying degrees of epithelial dysplasia. To gain insights into the molecular mechanisms associated with the progression of PMLs to OSCC, we profiled the whole transcriptome of 66 human PMLs comprising leukoplakia with dysplasia and hyperkeratosis non-reactive (HkNR) pathologies, alongside healthy controls and OSCC. Our data revealed that PMLs were enriched in gene signatures associated with cellular plasticity, such as partial EMT (p-EMT) phenotypes, and with immune response. Integrated analyses of the host transcriptome and microbiome further highlighted a significant association between differential microbial abundance and PML pathway activity, suggesting a contribution of the oral microbiome toward PML evolution to OSCC. Collectively, this study reveals molecular processes associated with PML progression that may help early diagnosis and disease interception at an early stage.

ß-catenin/CBP activation of mTORC1 signaling promotes partial epithelial-mesenchymal states in head and neck cancer.

Head and neck cancers, which include oral squamous cell carcinoma (OSCC) as a major subsite, exhibit cellular plasticity that includes features of an epithelial-mesenchymal transition (EMT), referred to as partial-EMT (p-EMT). To identify molecular mechanisms contributing to OSCC plasticity, we performed a multiphase analysis of single cell RNA sequencing (scRNAseq) data from human OSCC. This included a multiresolution characterization of cancer cell subgroups to identify pathways and cell states that are heterogeneously represented, followed by casual inference analysis to elucidate activating and inhibitory relationships between these pathways and cell states. This approach revealed signaling networks associated with hierarchical cell state transitions, which notably included an association between ß-catenin-driven CREB-binding protein (CBP) activity and mTORC1 signaling. This network was associated with subpopulations of cancer cells that were enriched for markers of the p-EMT state and poor patient survival. Functional analyses revealed that ß-catenin/CBP induced mTORC1 activity in part through the transcriptional regulation of a raptor-interacting protein, chaperonin containing TCP1 subunit 5 (CCT5). Inhibition of ß-catenin-CBP activity through the use of the orally active small molecule, E7386, reduced the expression of CCT5 and mTORC1 activity in vitro, and inhibited p-EMT-associated markers and tumor development in a murine model of OSCC. Our study highlights the use of multiresolution network analyses of scRNAseq data to identify targetable signals for therapeutic benefit, thus defining an underappreciated association between ß-catenin/CBP and mTORC1 signaling in head and neck cancer plasticity.

Structure learning for gene regulatory networks.

Inference of biological network structures is often performed on high-dimensional data, yet is hindered by the limited sample size of high throughput "omics" data typically available. To overcome this challenge, often referred to as the "small n, large p problem," we exploit known organizing principles of biological networks that are sparse, modular, and likely share a large portion of their underlying architecture. We present SHINE-Structure Learning for Hierarchical Networks-a framework for defining data-driven structural constraints and incorporating a shared learning paradigm for efficiently learning multiple Markov networks from high-dimensional data at large p/n ratios not previously feasible. We evaluated SHINE on Pan-Cancer data comprising 23 tumor types, and found that learned tumor-specific networks exhibit expected graph properties of real biological networks, recapture previously validated interactions, and recapitulate findings in literature. Application of SHINE to the analysis of subtype-specific breast cancer networks identified key genes and biological processes for tumor maintenance and survival as well as potential therapeutic targets for modulating known breast cancer disease genes.

Multi-modal profiling of peripheral blood cells across the human lifespan reveals distinct immune cell signatures of aging and longevity.

BACKGROUND: Age-related changes in immune cell composition and functionality are associated with multimorbidity and mortality. However, many centenarians delay the onset of aging-related disease suggesting the presence of elite immunity that remains highly functional at extreme old age. METHODS: To identify immune-specific patterns of aging and extreme human longevity, we analyzed novel single cell profiles from the peripheral blood mononuclear cells (PBMCs) of a random sample of 7 centenarians (mean age 106) and publicly available single cell RNA-sequencing (scRNA-seq) datasets that included an additional 7 centenarians as well as 52 people at younger ages (20-89 years). FINDINGS: The analysis confirmed known shifts in the ratio of lymphocytes to myeloid cells, and noncytotoxic to cytotoxic cell distributions with aging, but also identified significant shifts from CD4+ T cell to B cell populations in centenarians suggesting a history of exposure to natural and environmental immunogens. We validated several of these findings using flow cytometry analysis of the same samples. Our transcriptional analysis identified cell type signatures specific to exceptional longevity that included genes with age-related changes (e.g., increased expression of STK17A, a gene known to be involved in DNA damage response) as well as genes expressed uniquely in centenarians' PBMCs (e.g., S100A4, part of the S100 protein family studied in age-related disease and connected to longevity and metabolic regulation). INTERPRETATION: Collectively, these data suggest that centenarians harbor unique, highly functional immune systems that have successfully adapted to a history of insults allowing for the achievement of exceptional longevity. FUNDING: TK, SM, PS, GM, SA, TP are supported by NIH-NIAUH2AG064704 and U19AG023122. MM and PS are supported by NIHNIA Pepper center: P30 AG031679-10. This project is supported by the Flow Cytometry Core Facility at BUSM. FCCF is funded by the NIH Instrumentation grant: S10 OD021587.

Total RNA sequencing reveals gene expression and microbial alterations shared by oral pre-malignant lesions and cancer.

Head and neck cancers are a complex malignancy comprising multiple anatomical sites, with cancer of the oral cavity ranking among the deadliest and most disfiguring cancers globally. Oral cancer (OC) constitutes a subset of head and neck cancer cases, presenting primarily as tobacco-and alcohol-associated oral squamous cell carcinoma (OSCC), with a 5-year survival rate of ∻65%, partly due to the lack of early detection and effective treatments. OSCC arises from premalignant lesions (PMLs) in the oral cavity through a multi-step series of clinical and histopathological stages, including varying degrees of epithelial dysplasia. To gain insights into the molecular mechanisms associated with the progression of PMLs to OSCC, we profiled the whole transcriptome of 66 human PMLs comprising leukoplakia with dysplasia and hyperkeratosis non-reactive (HkNR) pathologies, alongside healthy controls and OSCC. Our data revealed that PMLs were enriched in gene signatures associated with cellular plasticity, such as partial EMT (p-EMT) phenotypes, and with immune response. Integrated analyses of the host transcriptome and microbiome further highlighted a significant association between differential microbial abundance and PML pathway activity, suggesting a contribution of the oral microbiome towards PML evolution to OSCC. Collectively, this study reveals molecular processes associated with PML progression that may help early diagnosis and disease interception at an early stage. AUTHOR SUMMARY: Patients harboring oral premalignant lesions (PMLs) have an increased risk of developing oral squamous cell carcinoma (OSCC), but the underlying mechanisms driving transformation of PMLs to OSCC remain poorly understood. In this study, Khan et al., analyzed a newly generated dataset of gene expression and microbial profiles of oral tissues from patients diagnosed with PMLs from differing histopathological groups, including hyperkeratosis not reactive ( HkNR ) and dysplasia, comparing these profiles with OSCC and normal oral mucosa. Significant similarities between PMLs and OSCC were observed, with PMLs manifesting several cancer hallmarks, including oncogenic and immune pathways. The study also demonstrates associations between the abundance of multiple microbial species and PML groups, suggesting a potential contribution of the oral microbiome to the early stages of OSCC development. The study offers insights into the nature of the molecular, cellular and microbial heterogeneity of oral PMLs and suggests that molecular and clinical refinement of PMLs may provide opportunities for early disease detection and interception.

Mosaic Chromosomal Alterations and Human Longevity.

Mosaic chromosomal alterations (mCAs) are structural alterations associated with aging, cancer, cardiovascular disease, infectious diseases, and mortality. The distribution of mCAs in centenarians and individuals with familial longevity is poorly understood. We used MOsaic CHromosomal Alteration (MoChA) to discover mCAs in 2050 centenarians, offspring, and 248 controls from the New England Centenarian Study (NECS) and in 3 642 subjects with familial longevity and 920 spousal controls from the Long-Life Family Study (LLFS). We analyzed study-specific associations of somatic mCAs with age, familial longevity, the incidence of age-related diseases, and mortality and aggregated the results by meta-analysis. We show that the accumulation of mCAs > 100 KB increased to 102 years and plateaued at older ages. Centenarians and offspring accumulated fewer autosomal mCAs compared with controls (relative risk 0.637, p = .0147). Subjects with the APOE E4 allele had a 35.3% higher risk of accumulating autosomal mCAs (p = .002). Males were at higher risk for mCAs compared to females (male relative risk 1.36, p = 5.15e-05). mCAs were associated with increased hazard for cancer (hazard ratio 1.2) and dementia (hazard ratio 1.259) at a 10% false discovery rate. We observed a borderline significant association between mCAs and risk for mortality (hazard ratio 1.07, p = .0605). Our results show that the prevalence of individuals with mCAs does not continue to increase at ages >102 years and factors promoting familial longevity appear to confer protections from mCAs. These results suggest that limited mCA accumulation could be an important mechanism for extreme human longevity that needs to be investigated.

Bayesian Differential Analysis of Cell Type Proportions.

The analysis of cell type proportions in a biological sample should account for the compositional nature of the data but most analyses ignore this characteristic with the risk of producing misleading conclusions. The recent method scCODA appropriately incorporates these constraints by using a Bayesian Multinomial-Dirichlet model that requires a reference cell type to normalize the distribution of all cell types. However, a reference cell type that is stable across biological conditions may not always be available. Here, we present an approach that uses a Bayesian multinomial regression for the analysis of single cell distribution data without the need for a reference cell type. We show an implementation example using the rjags package within the R software.

A metabolomic signature of the APOE2 allele.

With the goal of identifying metabolites that significantly correlate with the protective e2 allele of the apolipoprotein E (APOE) gene, we established a consortium of five studies of healthy aging and extreme human longevity with 3545 participants. This consortium includes the New England Centenarian Study, the Baltimore Longitudinal Study of Aging, the Arivale study, the Longevity Genes Project/LonGenity studies, and the Long Life Family Study. We analyzed the association between APOE genotype groups E2 (e2e2 and e2e3 genotypes, N = 544), E3 (e3e3 genotypes, N = 2299), and E4 (e3e4 and e4e4 genotypes, N = 702) with metabolite profiles in the five studies and used fixed effect meta-analysis to aggregate the results. Our meta-analysis identified a signature of 19 metabolites that are significantly associated with the E2 genotype group at FDR < 10%. The group includes 10 glycerolipids and 4 glycerophospholipids that were all higher in E2 carriers compared to E3, with fold change ranging from 1.08 to 1.25. The organic acid 6-hydroxyindole sulfate, previously linked to changes in gut microbiome that were reflective of healthy aging and longevity, was also higher in E2 carriers compared to E3 carriers. Three sterol lipids and one sphingolipid species were significantly lower in carriers of the E2 genotype group. For some of these metabolites, the effect of the E2 genotype opposed the age effect. No metabolites reached a statistically significant association with the E4 group. This work confirms and expands previous results connecting the APOE gene to lipid regulation and suggests new links between the e2 allele, lipid metabolism, aging, and the gut-brain axis.

Inactivation of LATS1/2 drives luminal-basal plasticity to initiate basal-like mammary carcinomas.

Basal-like breast cancers, an aggressive breast cancer subtype that has poor treatment options, are thought to arise from luminal mammary epithelial cells that undergo basal plasticity through poorly understood mechanisms. Using genetic mouse models and ex vivo primary organoid cultures, we show that conditional co-deletion of the LATS1 and LATS2 kinases, key effectors of Hippo pathway signaling, in mature mammary luminal epithelial cells promotes the development of Krt14 and Sox9-expressing basal-like carcinomas that metastasize over time. Genetic co-deletion experiments revealed that phenotypes resulting from the loss of LATS1/2 activity are dependent on the transcriptional regulators YAP/TAZ. Gene expression analyses of LATS1/2-deleted mammary epithelial cells notably revealed a transcriptional program that associates with human basal-like breast cancers. Our study demonstrates in vivo roles for the LATS1/2 kinases in mammary epithelial homeostasis and luminal-basal fate control and implicates signaling networks induced upon the loss of LATS1/2 activity in the development of basal-like breast cancer.

Serum Orotidine: A Novel Biomarker of Increased CVD Risk in Type 2 Diabetes Discovered Through Metabolomics Studies.

OBJECTIVE: To identify novel biomarkers of cardiovascular disease (CVD) risk in type 2 diabetes (T2D) via a hypothesis-free global metabolomics study, while taking into account renal function, an important confounder often overlooked in previous metabolomics studies of CVD. RESEARCH DESIGN AND METHODS: We conducted a global serum metabolomics analysis using the Metabolon platform in a discovery set from the Joslin Kidney Study having a nested case-control design comprising 409 individuals with T2D. Logistic regression was applied to evaluate the association between incident CVD events and each of the 671 metabolites detected by the Metabolon platform, before and after adjustment for renal function and other CVD risk factors. Significant metabolites were followed up with absolute quantification assays in a validation set from the Joslin Heart Study including 599 individuals with T2D with and without clinical evidence of significant coronary heart disease (CHD). RESULTS: In the discovery set, serum orotidine and 2-piperidinone were significantly associated with increased odds of incident CVD after adjustment for glomerular filtration rate (GFR) (odds ratio [OR] per SD increment 1.94 [95% CI 1.39-2.72], P = 0.0001, and 1.62 [1.26-2.08], P = 0.0001, respectively). Orotidine was also associated with increased odds of CHD in the validation set (OR 1.39 [1.11-1.75]), while 2-piperidinone did not replicate. Furthermore, orotidine, being inversely associated with GFR, mediated 60% of the effects of declining renal function on CVD risk. Addition of orotidine to established clinical predictors improved (P < 0.05) C statistics and discrimination indices for CVD risk (ΔAUC 0.053, rIDI 0.48, NRI 0.42) compared with the clinical predictors alone. CONCLUSIONS: Through a robust metabolomics approach, with independent validation, we have discovered serum orotidine as a novel biomarker of increased odds of CVD in T2D, independent of renal function. Additionally, orotidine may be a biological mediator of the increased CVD risk associated with poor kidney function and may help improve CVD risk prediction in T2D.

Cell Type Diversity Statistic: An Entropy-Based Metric to Compare Overall Cell Type Composition Across Samples.

Changes of cell type composition across samples can carry biological significance and provide insight into disease and other conditions. Single cell transcriptomics has made it possible to study cell type composition at a fine resolution. Most single cell studies investigate compositional changes between samples for each cell type independently, not accounting for the fixed number of cells per sample in sequencing data. Here, we provide a metric of the distribution of cell type proportions in a sample that can be used to compare the overall distribution of cell types across multiple samples and biological conditions. This is the first method to measure overall cell type composition at the single cell level. We use the method to assess compositional changes in peripheral blood mononuclear cells (PBMCs) related to aging and extreme old age using multiple single cell datasets from individuals of four age groups across the human lifespan.

Inhibition of LSD1 Attenuates Oral Cancer Development and Promotes Therapeutic Efficacy of Immune Checkpoint Blockade and YAP/TAZ Inhibition.

Lysine-specific demethylase 1 (LSD1) is a histone demethylase that contributes to the etiology of oral squamous cell carcinoma (OSCC) in part by promoting cancer stem cell phenotypes. The molecular signals regulated by LSD1, or acting with LSD1, are poorly understood, particularly in the development of OSSC. In this study, we show that conditional deletion of the Lsd1 gene or pharmacologic inhibition of LSD1 in the tongue epithelium leads to reduced development of OSCC following exposure to the tobacco carcinogen 4NQO. LSD1 inhibition attenuated proliferation and clonogenic survival and showed an additive effect when combined with the YAP inhibitor Verteporfin. Interestingly, LSD1 inhibition upregulated the expression of PD-L1, leading to immune checkpoint inhibitor therapy responses. IMPLICATIONS: Collectively, our studies reveal a critical role for LSD1 in OSCC development and identification of tumor growth targeting strategies that can be combined with LSD1 inhibition for improved therapeutic application.

Novel plasma exosome biomarkers for prostate cancer progression in co-morbid metabolic disease.

Comorbid Type 2 diabetes (T2D), a metabolic complication of obesity, associates with worse cancer outcomes for prostate, breast, head and neck, colorectal and several other solid tumors. However, the molecular mechanisms remain poorly understood. Emerging evidence shows that exosomes carry miRNAs in blood that encode the metabolic status of originating tissues and deliver their cargo to target tissues to modulate expression of critical genes. Exosomal communication potentially connects abnormal metabolism to cancer progression. Here, we hypothesized that T2D plasma exosomes induce epithelial-mesenchymal transition (EMT) and immune checkpoints in prostate cancer cells. We demonstrate that plasma exosomes from subjects with T2D induce EMT features in prostate cancer cells and upregulate the checkpoint genes CD274 and CD155. We demonstrate that specific exosomal miRNAs that are differentially abundant in plasma of T2D adults compared to nondiabetic controls (miR374a-5p, miR-93-5p and let-7b-3p) are delivered to cancer cells, thereby regulating critical target genes. We build on our previous reports showing BRD4 controls migration and dissemination of castration-resistant prostate cancer, and transcription of key EMT genes, to show that T2D exosomes require BRD4 to drive EMT and immune ligand expression. We validate our findings with gene set enrichment analysis of human prostate tumor tissue in TGCA genomic data. These results suggest novel, non-invasive approaches to evaluate and potentially block progression of prostate and other cancers in patients with comorbid T2D.

Evaluating the Homogeneity of Surface Features Induced by Impact-Based Surface Treatments.

Impact surface treatments are well-known for their efficiency in enhancing the mechanical properties of metallic materials, especially under cyclic loadings. These processes, which encompass a wide range of surface treatments based on repetitive impacts of tools of various types, induce surface plastic deformation, compressive residual stresses, and grain refinement alter the surface roughness as a side effect. Thus, it is essential to have suitable indexes to quantify the surface features caused by the typically random nature of these treatments. Herein, we evaluated the rationality of using standard roughness parameters for describing the morphological characteristics of surfaces treated by shot peening as a representative and widely used treatment of the category. A detailed numerical model of the peening process was developed. The output data were elaborated to extract the surface roughness parameters following the standard procedures. The results revealed the validity of the surface roughness parameters to describe the topography of material treated with adequate surface coverage, also highlighting the necessity to use a set of parameters rather than the common practice of relying on single parameters. Not considering a comprehensive set of amplitude and spacing parameters can result in significant, inconsistent, and misleading results while comparing the performance of surfaces.