In response to the Song et al. comments, which raised concerns about the methodologies and conclusions of our work,1 we would like to clarify some points that should not be misinterpreted. Song et al. partially analysed some of the data described in our paper to support their comments. Despite the criticism, our methodology is based on relevant literature, and we provide a detailed method description to enable replication of our study. By contrast, the letter by Song et al. did not include a methodology that would allow the readers to reproduce the data therein mentioned. The authors suggest a different parameter and cut-off thresholds to re-analyse some data sets used in our study, ultimately providing different results. It is known that changing parameters and analysis workflows produce diverse results. The dilemma of selecting bioinformatics strategies or statistical criteria to analyse gene expression remains a matter of no consensus. The approach suggested by Song et al. is based on the selection of different cut-offs for each cancer type to ensure that different types of cancers present a similar and very high number of differentially expressed genes (DEGs). The authors also suggest that our fixed criterion of DEGs cut-off, used as a universal threshold for DEGs determination among all tested cancers, might not be suitable due to intrinsic differences between tumours with different tissue origins. To the best of our knowledge, no previous pan-cancer study has used the criterion proposed by the authors, which includes, for some tumour types, most of the detected genes as differentially expressed. In contrast, a fixed cut-off to determine DEGs has been widely used to compare gene expression among cancer types in pan-cancer studies2-5 and in tools for large public multi-omics data analyses.6-9 Besides, as previously demonstrated, each type of normal or tumour tissue has a specific and highly variable dimension of the secretome components10, 11 that cannot be shown if we assume the same number of DEGs for all tissues. Accordingly, we found that pancreatic cancer presented the highest number of upregulated secretome genes (1267) across all tumour types. This result is expected for the pancreatic tissues, which have more than 70% of the transcripts encoding secreted proteins in the normal tissues.11 Thus, our data show that these tissue-specific characteristics, captured by the analyses that we performed considering the same statistical cut-off for all tumour types, may influence the prevalence and severity of the cachexia in different tumour types. Song et al. did not mention the data sets and samples tested, as well as the statistical analysis and tools used to re-analyse the data. Tables and figures do not have captions, which hinders a careful examination and interpretation of their study. Regarding the criticism of the statistical analyses conducted in our study, we would like to clarify that we used Gene Expression Profiling Interactive Analysis (GEPIA, http://gepia.cancer-pku.cn/)9 to integrate the TCGA (The Cancer Genome Atlas) and GTEx (Genotype-Tissue Expression) data. In this case, if the authors had used the same tool, their option for analysing the data would be the four-way analysis of variance (ANOVA) or LIMMA. We calculated the DEGs between tumour and normal samples by four-way ANOVA using disease state (tumour or normal), ethnicity, age, and gender as variables. ANOVA measures the strength of the relationship between the log-transformed expression of a gene and all variables mentioned above, between tumour and normal tissues, as previously described.4 We used the standard Benjamini and Hochberg false discovery rate (FDR) method to adjust the P value in each factor to obtain the multiple testing adjusted q value, and not fold change = 2 and P < 0.01 as statistical cut-offs, as reported by Song et al. The DEGs between tumour and normal samples were determined by ANOVA, applying the statistical cut-offs of log2 fold change >1 and q value <0.01. FDR is recommended in high-throughput experiments to correct random events that falsely appear significant, exerting robust control over the error rate even when the hypotheses have dependencies.12 Thus, despite the fixed criterion of DEGs cut-off used by the authors, the number of DEGs identified in their analyses includes many false positives. Moreover, considering the extremely high number of DEGs (~9000) found by the authors for each tumour type, nearly all cachexia-inducing factors (CIFs) are altered. Consequently, their data fail to identify CIFs tumour-specific expression profiles correlated with the prevalence of weight loss and cachexia in different tumour types. The authors also raised concerns about our survival analysis using cancer gene expression data. Access to SurvExpress database13 has been lost since October 2019, and the platform is currently out of funds (https://tinyurl.com/y42v5ovf). However, the data sets have recently been made available for analysis. We analysed three additional pancreatic cancer data sets14, 15 to further validate the expression of CIFs genes as predictors of cancer survival outcomes. The CIFs gene expression profile for each of these data sets was significantly associated with poor prognosis (Figure 1A). Indeed, the prognostic value (survival) of the CIFs gene expression was further confirmed in pancreatic cancer data sets by the area under receiver operating characteristic (ROC) curves statistics (area under the curve, AUC > 0.6 for all data sets) (Figure 1B). We did not carry out the survival analysis for individual CIFs because our data indicate that cachexia outcome is related to a specific set of CIFs rather than a single factor. These results explain our option in performing the survival analysis using the SurvExpress tool, which was developed to assess the prognostic performance of gene expression signatures. The RNA-Seq data from TCGA and GTEx databases constitute a rich source of information for cachexia research. Though we lack a validation set for CIFs gene expression profile, the interpretations of our transcriptomics results are clinically and biologically relevant. Previous pan-cancer studies that used samples exclusively from the TCGA data set have provided a uniquely comprehensive, in-depth, and interconnected understanding of several aspects of human tumours.10, 16-19 Thus, the lack of further validation does not negatively affect the value of our findings. Moreover, one must consider the broad spectrum of tumours that were compared by robust statistical analysis applied to integrate such big data, including 12 tumour types (4651 and 2737 tumour and normal samples, respectively). Finally, the molecular analyses of several primary tumours and normal samples have great potential to impact cachexia research positively. The integration of multi-omics data, including single-cell data, brings the significant potential to understand the syndrome's complexity. However, it should be emphasized that the repertoire of different analyses of such big data might lead to different interpretations and results, as argued before. These aspects should be considered by researchers when interpreting results from systems biology studies that interrogate cachexia comprehensively. This study was supported by the São Paulo Research Foundation—FAPESP (grants 12/13961-6, 13/50343-1, 14/13941-0, and 18/19695-2) and by the National Council for Scientific and Technological Development, CNPq (Process 311530/2019-2 to RFC). The results published here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga. The authors certify that the manuscript complies with the ethical guidelines for publishing in the Journal of Cachexia, Sarcopenia and Muscle.22